A novel approach to investigate tissue-specific trinucleotide repeat instability
© Lee et al; licensee BioMed Central Ltd. 2010
Received: 25 September 2009
Accepted: 19 March 2010
Published: 19 March 2010
In Huntington's disease (HD), an expanded CAG repeat produces characteristic striatal neurodegeneration. Interestingly, the HD CAG repeat, whose length determines age at onset, undergoes tissue-specific somatic instability, predominant in the striatum, suggesting that tissue-specific CAG length changes could modify the disease process. Therefore, understanding the mechanisms underlying the tissue specificity of somatic instability may provide novel routes to therapies. However progress in this area has been hampered by the lack of sensitive high-throughput instability quantification methods and global approaches to identify the underlying factors.
Here we describe a novel approach to gain insight into the factors responsible for the tissue specificity of somatic instability. Using accurate genetic knock-in mouse models of HD, we developed a reliable, high-throughput method to quantify tissue HD CAG repeat instability and integrated this with genome-wide bioinformatic approaches. Using tissue instability quantified in 16 tissues as a phenotype and tissue microarray gene expression as a predictor, we built a mathematical model and identified a gene expression signature that accurately predicted tissue instability. Using the predictive ability of this signature we found that somatic instability was not a consequence of pathogenesis. In support of this, genetic crosses with models of accelerated neuropathology failed to induce somatic instability. In addition, we searched for genes and pathways that correlated with tissue instability. We found that expression levels of DNA repair genes did not explain the tissue specificity of somatic instability. Instead, our data implicate other pathways, particularly cell cycle, metabolism and neurotransmitter pathways, acting in combination to generate tissue-specific patterns of instability.
Our study clearly demonstrates that multiple tissue factors reflect the level of somatic instability in different tissues. In addition, our quantitative, genome-wide approach is readily applicable to high-throughput assays and opens the door to widespread applications with the potential to accelerate the discovery of drugs that alter tissue instability.
Expansions of trinucleotide repeat sequences over certain thresholds cause more than 30 human diseases including Huntington's disease (HD), a number of spinocerebellar ataxias (SCAs), myotonic dystrophy 1 (DM1), and fragile X syndrome. Interestingly, expanded trinucleotide repeat sequences undergo progressive, expansion-biased tissue-specific somatic instability [1–6]. As the severity of these disorders is highly dependent on repeat length, somatic instability in tissues that are the pathogenic targets is predicted to contribute to disease. Notably, in HD, striking somatic expansion of the HD CAG repeat occurs in the striatum and cortex, brain regions that are major targets of the pathogenic process. Furthermore, studies both in HD patients and in a knock-in mouse model of HD provide compelling evidence indicating that somatic expansion in these brain regions accelerates the ongoing pathogenic process [7–9]. Therefore, understanding the mechanisms underlying tissue-specific somatic instability in HD may provide novel routes to therapies.
Somatic instability is critically dependent on DNA repair genes and is also influenced by cis-factors [7, 8, 10–16]. However, it is unknown what determines its tissue specificity. It has been proposed that the expression levels of DNA repair genes and/or the pathogenic process itself may underlie tissue patterns of instability . Given that somatic HD CAG instability occurs in many tissues to varying extents [3, 6, 17], we reasoned firstly, that tissue specificity may governed by many factors, and secondly, that studying a large cross-section of tissues with different instabilities would provide the most information concerning the major factors underlying tissue instability patterns. Therefore, in order to gain insight into the factors that govern the tissue specificity CAG instability in HD, we have taken quantitative, global and unbiased approaches.
Using accurate genetic knock-in mouse models of HD [6, 18] that exhibit similar tissue-specific patterns of somatic instability to those seen in HD patients [3, 6], we developed a novel instability quantification method that is sensitive and applicable to high-throughput assays. We then integrated this methodology with unbiased and global bioinformatic approaches to identify a gene expression "signature" and biological pathways that correlate with tissue instability. Using these methods we have, a) tested the role played by factors previously proposed to contribute to the tissue specificity of somatic instability, and b) uncovered novel pathways that may be important in determining the tissue specificity of instability in HD.
Validation of the relative peak height threshold quantification method
We first determined the reproducibility of our method by quantifying instability index in 17 tissues from 2-6 different HdhQ 111/+mice at 5 months of age. As shown in Figure 2, the shift in the highest peak compared to tail (panel A) and the instability index (panel B) were highly reproducible between mice for all tissues tested. Note that the instability indices of stable tissues (i.e. lung, heart, spleen) were negative because stable tissue GeneMapper traces were biased toward contraction likely due to the increased amplification efficiency of shorter CAG alleles. Instability indices of 17 tissues ranged from -1.03 (testis) to 6.37 (striatum).
The relative peak height threshold method can also be applied to different types of instability quantification depending on the focus of the biological question. Thus, after applying the relative peak height threshold, we can determine contraction and expansion indices (Figure 2C), the number of contracted and expanded peaks (Figure 2D) or the relative composition (%) of contracted, expanded and unchanged peaks (Figure 2E). Importantly, these measurements of different aspects of instability may be useful to capture the complexity of tissue instability. In all cases, measurements were reproducible for all tissues across multiple mice. To represent the levels of instability of tissues for further analysis, we used the instability index (Figure 2B).
We then compared instability indices using our relative peak height threshold method to somatic instability quantified using SP-PCR on genomic DNA of tissues from the same mouse (9 tissues, 5 month, HdhQ 111/+). Figure 3B shows examples of tissues exhibiting high, medium and low instability indices, and the corresponding CAG repeat length frequency distributions obtained by SP-PCR. These data indicated that the instability index broadly captured the bulk of the somatic variation detected by SP-PCR, but not the rare large expansions. However, there was a highly significant correlation between the instability index obtained using the bulk DNA method and an instability index quantified from the small pool data (Figure 3C, p value, 0.00015), suggesting that although instability index using bulk DNA may not be sensitive enough to detect rare molecules, it can give a good estimate of overall instability.
Together, our analyses suggest that the instability index, determined from GeneMapper traces of bulk genomic DNA, is a reproducible measurement, relatively insensitive to input DNA amount and well suited for high-throughput analyses where SP-PCR may be impractical.
Genome-wide identification of an instability-correlated gene expression signature
With the aim of investigating the tissue specificity of somatic instability in a global and unbiased manner we then took a bioinformatics approach. Using 16 different tissues from 5-month HdhQ 111/+mice as our training set (Figure 2B, excluding tail), with instability index as a quantitative phenotype, we analyzed mouse tissue gene expression data (Mouse Gene Expression Atlas GSE11339, C57BL/6J, 10 weeks) to identify a gene expression signature that correlated with tissue repeat instability. HdhQ 111somatic instability (and therefore instability index) increases over time . We chose 5 months as this represents a time-point at which tissue differences in instability can be readily resolved. Notably, the Gene Expression data is derived from mice that differ in age and genetic background (B6 versus CD1, absence versus presence of HD CAG knock-in allele) to the HdhQ 111mice in this study. While age and genetic background-related gene expression changes will increase the noise in our system, this broad, tissue-based analysis allows us to pull out major tissue-specific gene expression differences that occur over and above age- and genetic background-related effects.
We then confirmed the predictive power of this instability-correlated gene expression signature by comparing measured instability indices with predicted instability indices from our regression model in new independent samples. For this, 1) we measured instability indices of four new independent HdhQ 111/+tissues (muscle, olfactory bulb, white adipose tissue and adrenal gland) and compared these with instability indices predicted from the regression model in the same tissues (Figure 4, blue), and 2) we predicted instability indices using independent microarray data from HdhQ 111striatum and cerebellum and compared these with measured instability indices (Figure 4, red). As shown in Figure 4, the predicted and measured instability indices matched closely in all cases (test set RMSEP, 0.5444) with a significant correlation (Pearson correlation coefficient, 0.9783; p value, 9.6 × 10-7), indicating that instability index can be relatively precisely predicted from the gene expression signature. Furthermore, these data demonstrate that although the model was based on gene expression data and instability index data from mice that differed in age and genetic background, it nevertheless has significant predictive power. This indicates the presence of tissue-specific factors related to instability independent of age and genetic background.
Tissue instability prediction
Tissue instability index predicted by a PLSR model.
common myeloid progenitor
dendritic plasmacytoid B220+
granulo mono progenitor
macrophage peri LPS thio 0 hrs
mega erythrocyte progenitor
dendritic cells lymphoid CD8a+
macrophage bone marrow 6 hr LPS
B-cells marginal zone
dorsal root ganglia
macrophage bone marrow 24 h LPS
macrophage bone marrow 2 hr LPS
mammary gland lact
macrophage peri LPS thio 1 hrs
osteoblast day 14
cerebral cortex prefrontal
macrophage peri LPS thio 7 hrs
mammary gland non-lactating
mast cells IgE
thymocyte DP CD4+CD8+
macrophage bone marrow 0 hr
thymocyte SP CD8+
osteoblast day 21
mast cells IgE+antigen 1 hr
dendritic cells myeloid CD8a-
mast cells IgE+antigen 6 hr
retinal pigment epithelium
thymocyte SP CD4+
Pathogenesis and instability
To test the prediction that somatic instability does not occur as a consequence of ongoing pathogenesis, we performed two genetic experiments. Since the expanded Hdh CAG repeat is both a source of a pathogenic process and a target of instability, it is very difficult to delineate the relationship between the HD pathogenic process and somatic instability. Therefore, we used genetic mouse models in which neurodegenerative processes are modulated or caused by factors independent of the HD CAG repeat. We first investigated HdhQ 92mice lacking the dopamine transporter (DAT), which show accelerated HD pathogenesis in the striatum . As shown in Figure 5B, striatal instability indices of HdhQ 92/+DAT-/- and HdhQ 92/+DAT+/+ mice were not different, indicating that HD CAG instability is not contributed by the disease process. We also tested whether inducing neurodegeneration in the cerebellum, a normally stable tissue, would cause instability in the cerebellum by crossing HdhQ 111mice to Harlequin (Hq) mice, a model of cerebellar granule cell degeneration . As shown in Figure 5C, HdhQ 111/+Hq/Y mice and HdhQ 111/++/Y control mice exhibited similar low cerebellar instability indices, indicating that neurodegeneration per se is insufficient to induce instability.
Taken together, these results support the prediction from our mathematical model, that the HD CAG disease process is not responsible for the striatal specificity of HD CAG repeat instability, arguing against the sequestration of DNA repair proteins or other factors, as a contributor to somatic instability as previously suggested . Our results are also in agreement with similar levels of instability seen in knock-in and fragment transgenic models of HD that exhibit different rates of inclusion formation , and with the observation that striatal instability occurs in SCA1 and DM1, although the striatum is not the target of pathogenesis in these disorders [2, 5].
DNA repair and repeat instability
Selective neuronal expression of Msh3 was recently proposed to contribute to the greater levels of instability in neurons compared to glia , and therefore we explored this further. Analyses of gene expression data revealed nearly identical Msh3 expression levels in purified neurons and glia (Figure 6C). Together with the lack of correlation between instability index and Msh3 expression levels across 16 tissues (Figure 6A, Additional file 2), the data argue against a major role for Msh3 expression levels in determining tissue- or cell type-specific instability.
Therefore, our results suggest that although certain DNA repair genes are absolutely critical for somatic instability [7, 8, 11–16], their expression levels are unlikely to be the primary determinants of tissue specificity. Clearly, posttranscriptional and/or posttranslational regulation of DNA repair genes could still play a tissue-specific role. It would therefore be of further interest to determine whether there is a correlation between DNA repair enzyme activities and tissue instability.
Genome-wide survey for pathways that correlate with tissue instability
Pathways significantly correlated with the instability index.
G1 to S cell cycle reactome
Negative regulation of progression through cell cycle
M phase of mitotic cell cycle
Protein kinase inhibitor activity
Cell cycle pathway
Mitotic cell cycle
Kinase inhibitor activity
Protein amino acid-ribosylation
Integrin mediated cell adhesion
RNA helicase activity
UDP-galactose beta-N-acetylglucosamine beta-1,3-galactosyltransferase activity
Amine receptor activity
Mono amine GPCRS
Neuromuscular junction development
Serotonin receptor activity
Oxidoreductase activity, acting on the CH-CH groups of donors, oxygen as acceptor
It is possible that as striatum is particularly unstable, the highly correlated pathways are simply those that are predominantly present or absent in this tissue, and that the correlation with instability is coincidental. However, pathways significantly up-regulated or down-regulated in striatum compared to cerebellum (data not shown) showed little overlap with those that correlated with instability; for example, the dopamine pathway is strongly up-regulated in striatum, but does not correlate with instability. This suggests that the instability-correlated pathways are directly related to instability rather than simply being striatal-specific.
Test of prediction from GSEA
Contribution of multiple processes to somatic instability
We have developed a novel approach for use in investigations of tissue-specific somatic HD CAG repeat instability that combines a reliable, high-throughput method for quantifying somatic instability with mathematical modeling based on gene expression data. Predictions based on our modeling were confirmed using genetic, biochemical and cell culture-based experiments, indicating the validity of our bioinformatics approach.
It has been proposed that somatic instability may be a consequence of disease pathogenesis , potentially explaining the striatal specificity of somatic expansion in HD. Our results directly demonstrate that HD pathogenesis does not explain the tissue specificity of HD CAG instability. In addition, DNA repair proteins have been found to be essential factors for somatic instability of trinucleotide repeats [7, 8, 11–16]. However, here we demonstrate that differences in expression levels of DNA repair genes do not underlie the tissue-specific differences in HD CAG instability. In addition, Hdh expression levels did not correlate with instability index in tissues (data not shown), confirming observations that although transcription through expanded repeats may be important in somatic instability , tissue-specific patterns are not reflected in the steady state levels of Hdh mRNA. Alternatively, our study suggests new pathways, notably metabolism, neurotransmitter, and cell cycle that may contribute, in combination, to the level of somatic instability in different tissues, providing a starting point to identify additional factors that contribute to somatic instability. Notably, there was no predominant factor that could explain the tissue-specificity of HD CAG instability, suggesting that patterns of instability are determined by the combined effects/interactions of many genes.
Somatic instability of trinucleotide repeats not only requires trans-acting factors, but has also been shown to depend on cis-acting sequences . Thus, while certain tissues are more predisposed to somatic expansion, the expandability of a particular repeat in a particular tissue is further modified by its context. This could at least in part explain differences in the precise tissue patterns of somatic expansion in different diseases [2, 29]. It would therefore be of interest to determine instability-correlated gene expression signatures and instability-correlated biological pathways for other trinucleotide repeat diseases. Instability-correlated genes/pathways that are shared between diseases would provide further insight into fundamental aspects of tissue-specific instability.
Our bioinformatics method based on gene expression data can only address aspects of tissue instability that are related to steady-state mRNA levels. In principle, however, a similar bioinformatic approach could be also applied to proteomics data. Irrespective of the particular method however, the strength of our approach is in its high-throughput, global and predictive nature, facilitating a number of important applications. Our GeneMapper quantification method is readily applicable to high-throughput assays such as screening small molecules that modulate instability in cells, or screening for genetic modifiers in mice. A powerful application of our bioinformatics approach is that the instability-correlated gene expression signature can be used as a surrogate marker for instability in situations where repeat instability cannot be directly measured. For example, gene expression databases can be screened to identify cell or tissue states that have the propensity for somatic instability, even in the absence of an expanded CAG repeat target as a read-out. Similarly, databases can be screened for compounds that reduce the instability propensity. Together, these approaches promise to accelerate the discovery of drugs that modulate instability and that are therefore candidate modifiers of disease.
Our study demonstrates that multiple tissue factors including metabolism, neurotransmitter, and cell cycle combine to reflect the level of somatic instability in different tissues. Our findings also indicate that DNA repair proteins act largely in a non tissue-specific manner. In addition, the combination of our instability quantification method and mathematical modeling is a powerful strategy that has allowed us, in an unbiased manner, to gain critical new insights into the tissue specificity of trinucleotide repeat instability in HD. It opens the door to widespread downstream applications with the potential to make significant advances in novel avenues for therapeutic intervention in both Huntington's disease and trinucleotide expansion disorders in general.
HdhQ 111knock-in mice with 109 CAGs  were used for quantification of tissue instability and for microarray gene expression analyses (Affymetrix MG 430 2.0). Mice were genotyped as previously described . For accelerated pathology models in cerebellum or striatum, HdhQ 111/+(CD1) and HdhQ 92/+mice (CD1)  were crossed with Harlequin (Hq) mutant (B6CBACa-AW-J/A)  and dopamine transporter (DAT) knockout mice (C57Bl/6J) , respectively. HdhQ 92mice were crossed with DAT knock-out mice and progeny intercrossed to generate HdhQ 92/Q92 DAT-/- mice and HdhQ 92DAT+/+ control littermates for comparisons of instability. HdhQ 111males were crossed with Hq/+ females, and HdhQ 111/+Hq/Y males and control HdhQ 111/++/Y littermate males used for comparisons of instability. All animal experiments were performed to minimize pain and discomfort, under an approved Institutional Animal Care and Use Committee protocol.
CAG length determination and instability quantification
Genomic DNA, isolated from mouse tissues and cell lines (DNeasy, Qiagen), was used for PCR amplification using HD CAG repeat-specific primers as previously described . The forward primer was fluorescently labeled with 6-FAM (Perkin Elmer) and PCR products were resolved using the ABI 3730 DNA analyzer (Applied Biosystems) using GeneMapper v.3.7 and GeneScan 500-LIZ as internal size standard to assign repeat size. GeneMapper traces were used to determine an instability index as described (Figure 1).
Genomic DNA was digested with Eco RV and diluted in 10 mM Tris-HCl, pH 8.0, 1 mM EDTA containing 0.1 μM carrier primer (MD16) to a final concentration of approximately 10 ng/μl. The amount of input DNA equivalent to a single amplifiable mutant Hdh allele was determined empirically using Poisson analysis, and for each tissue between 32 and 117 single mutant amplifiable molecules were analyzed. A nested PCR protocol was used, in which only the mutant (knock-in) Hdh allele is amplified. Mutant Hdh alleles were amplified using 0.5 μM MD16 primer 5'-CCCATTCATTGCCTTGCTGCTAAG (forward)  and 0.5 μM LKH5 primer 5'-TGGGTTGCTGGGTCACTCTGTC (reverse)  in 1× Thermo Scientific Custom PCR mix (containing 45 mM Tris-HCl pH 8.8, 11 mM ammonium sulfate, 4.5 mM MgCl2, 6.7 mM 2-mercaptoethanol, 4.4 μM EDTA, 1 mM dNTPs and 113 μg/ml BSA), 10% DMSO and 0.5 U units Taq polymerase (Fisher). Cycling conditions were 94°C 5 min, 35 cycles of 94°C 30 sec, 58°C 30 sec, 72°C 3 min, followed by 10 minutes at 72°C. PCR products were diluted 100-fold in TE and amplified in a second round using 0.8 μM Hu4 primer 5'-CCTGGAAAAGCTGATGAAGG (forward) and 0.8 μM Hu3 primer 5'-GGCGGCTGAGGAAGCTGAGGA (reverse) in a PCR buffer containing 67 mM Tris-HCl pH 8.8, 16.7 mM (NH4)2SO4, 2 mM MgCl2, 0.17 mg/mg BSA, 10 mM 2-mercaptoethanol, 10% DMSO, 200 μM dNTPs, with 0.5 U Taq polymerase (Fisher). Cycling conditions were 94°C 90 sec, 25 cycles of 94°C 30 sec, 65°C 30 sec, 72°C 90 sec, followed by 10 minutes at 72°. Hu4 was fluorescently labelled with 6-FAM (Applied Biosystems). PCR products were resolved using the ABI 3730 automated DNA analyzer (Applied Biosystems) using GeneMapper v.3.7 and GeneScan 500-LIZ as internal size standard to assign repeat size. HD CAG size was assigned as the highest peak. All PCR reactions were set up in a laminar flow hood and 20% of zero DNA control PCR reactions were included per run. To determine a small pool instability index we determined the frequency of each CAG repeat length, and multiplied each frequency by the number of repeats (+ or -) from the modal CAG length. These values were then summed.
Analysis of GNF mouse Gene Expression Atlas and regression modeling
We used the mouse tissue gene expression database of Genomics Institute of the Novartis Research Foundation (mouse Gene Expression Atlas, GSE11339). All microarrays were background corrected and normalized by gcRMA. To identify an instability-correlated gene expression signature, Pearson correlation coefficients and corresponding p values between gene expression levels and instability indices of training samples (16 tissues, 2 gene expression replicates) were calculated for each probe, and the gene expression data was sorted by p values. We used Pearson correlation coefficients only as a ranking metric and this linear correlation information has not been used in actual modeling. Therefore, our models capture not just linear relationship but covariance between instability and expression. To identify an instability-correlated gene expression signature, we sequentially introduced the top n most highly correlated probes into the regression algorithms in a forward selection procedure, and calculated root mean squared error of prediction (RMSEP) by leave one out cross validation (LOO CV) of training samples (R, 2.4.1 and 'pls' package, 2.5.0). In addition to LOO CV of training samples, we further tested our model using 2 different test set samples. Firstly, we measured instability indices in additional tissues (muscle, olfactory bulb, white adipose tissue and adrenal gland (HdhQ 111/+, 5 months, n = 4-6 mice for each tissue) and compared them with instability indices predicted by our model. Secondly, we additionally analyzed gene expression profiles of striatum and cerebellum (HdhQ 111/+, 5 months, n = 1) and used these to predict instability indices for comparison to previously measured instability indices in these tissues. Test set RMSEP was calculated based on the difference between measured and predicted instability indices. Prediction of instability index for each of the tissues analyzed in mouse Gene Expression Atlas was based on our regression model and the instability-correlated signature.
Gene set enrichment analysis
Using all probes, gene set enrichment analysis  was performed to sensitively identify significantly correlated pathways with instability index. Measured instability indices of training samples (16 tissues, Figure 2A) were used as continuous phenotype labels, and Pearson correlation was selected for a ranking metric. Our gene set database included pathways annotated by Gene Ontology, KEGG, GenMAPP, and the Molecular Signature Database from the Broad Institute. Significant gene sets were identified by permutation-based nominal p value (p < 0.01).
- HD :
myotonic dystrophy type 1
partial least square regression
- LOO CV:
leave-one-out cross validation.
We thank James W. MacDonald, Danh V. Nguyen, Jason M. Laramie, and Partners Research Computing (Jerry Xu and Dennis Gurgle) for technical assistance. Supported by NINDS grants NS049206 (to VCW), NS16367 (to JFG and MEM, HD Center Without Walls), and NS32765 (to MEM), NCBC grant LM008748 (to IK), the Hereditary Disease Foundation (to MC), Canada Research Chair in Molecular Neuropharmacology (to MC), and the Huntington's Disease Society of America (Coalition for the Cure, to JFG and MEM). We also thank Delta Squared for supporting this project.
- Ashizawa T, Dubel JR, Harati Y: Somatic instability of CTG repeat in myotonic dystrophy. Neurology. 1993, 43: 2674-2678.View ArticlePubMedGoogle Scholar
- Fortune MT, Vassilopoulos C, Coolbaugh MI, Siciliano MJ, Monckton DG: Dramatic, expansion-biased, age-dependent, tissue-specific somatic mosaicism in a transgenic mouse model of triplet repeat instability. Hum Mol Genet. 2000, 9: 439-445. 10.1093/hmg/9.3.439View ArticlePubMedGoogle Scholar
- Kennedy L, Evans E, Chen CM, Craven L, Detloff PJ, Ennis M, Shelbourne PF: Dramatic tissue-specific mutation length increases are an early molecular event in Huntington disease pathogenesis. Hum Mol Genet. 2003, 12: 3359-3367. 10.1093/hmg/ddg352View ArticlePubMedGoogle Scholar
- Kennedy L, Shelbourne PF: Dramatic mutation instability in HD mouse striatum: does polyglutamine load contribute to cell-specific vulnerability in Huntington's disease?. Hum Mol Genet. 2000, 9: 2539-2544. 10.1093/hmg/9.17.2539View ArticlePubMedGoogle Scholar
- Watase K, Venken KJ, Sun Y, Orr HT, Zoghbi HY: Regional differences of somatic CAG repeat instability do not account for selective neuronal vulnerability in a knock-in mouse model of SCA1. Hum Mol Genet. 2003, 12: 2789-2795. 10.1093/hmg/ddg300View ArticlePubMedGoogle Scholar
- Wheeler VC, Auerbach W, White JK, Srinidhi J, Auerbach A, Ryan A, Duyao MP, Vrbanac V, Weaver M, Gusella JF, et al.: Length-dependent gametic CAG repeat instability in the Huntington's disease knock-in mouse. Hum Mol Genet. 1999, 8: 115-122. 10.1093/hmg/8.1.115View ArticlePubMedGoogle Scholar
- Dragileva E, Hendricks A, Teed A, Gillis T, Lopez ET, Friedberg EC, Kucherlapati R, Edelmann W, Lunetta KL, MacDonald ME, Wheeler VC: Intergenerational and striatal CAG repeat instability in Huntington's disease knock-in mice involve different DNA repair genes. Neurobiol Dis. 2009, 33: 37-47. 10.1016/j.nbd.2008.09.014PubMed CentralView ArticlePubMedGoogle Scholar
- Wheeler VC, Lebel LA, Vrbanac V, Teed A, te Riele H, MacDonald ME: Mismatch repair gene Msh2 modifies the timing of early disease in Hdh(Q111) striatum. Hum Mol Genet. 2003, 12: 273-281. 10.1093/hmg/ddg056View ArticlePubMedGoogle Scholar
- Swami M, Hendricks AE, Gillis T, Massood T, Mysore J, Myers RH, Wheeler VC: Somatic expansion of the Huntington's disease CAG repeat in the brain is associated with an earlier age of disease onset. Hum Mol Genet. 2009, 18: 3039-3047. 10.1093/hmg/ddp242PubMed CentralView ArticlePubMedGoogle Scholar
- Cleary JD, Pearson CE: The contribution of cis-elements to disease-associated repeat instability: clinical and experimental evidence. Cytogenet Genome Res. 2003, 100: 25-55. 10.1159/000072837View ArticlePubMedGoogle Scholar
- Foiry L, Dong L, Savouret C, Hubert L, te Riele H, Junien C, Gourdon G: Msh3 is a limiting factor in the formation of intergenerational CTG expansions in DM1 transgenic mice. Hum Genet. 2006, 119: 520-526. 10.1007/s00439-006-0164-7View ArticlePubMedGoogle Scholar
- Gomes-Pereira M, Fortune MT, Ingram L, McAbney JP, Monckton DG: Pms2 is a genetic enhancer of trinucleotide CAG.CTG repeat somatic mosaicism: implications for the mechanism of triplet repeat expansion. Hum Mol Genet. 2004, 13: 1815-1825. 10.1093/hmg/ddh186View ArticlePubMedGoogle Scholar
- Manley K, Shirley TL, Flaherty L, Messer A: Msh2 deficiency prevents in vivo somatic instability of the CAG repeat in Huntington disease transgenic mice. Nat Genet. 1999, 23: 471-473. 10.1038/70598View ArticlePubMedGoogle Scholar
- Owen BA, Yang Z, Lai M, Gajek M, Badger JD, Hayes JJ, Edelmann W, Kucherlapati R, Wilson TM, McMurray CT: (CAG)(n)-hairpin DNA binds to Msh2-Msh3 and changes properties of mismatch recognition. Nat Struct Mol Biol. 2005, 12: 663-670. 10.1038/nsmb965View ArticlePubMedGoogle Scholar
- Savouret C, Brisson E, Essers J, Kanaar R, Pastink A, te Riele H, Junien C, Gourdon G: CTG repeat instability and size variation timing in DNA repair-deficient mice. Embo J. 2003, 22: 2264-2273. 10.1093/emboj/cdg202PubMed CentralView ArticlePubMedGoogle Scholar
- Broek van den WJ, Nelen MR, Wansink DG, Coerwinkel MM, te Riele H, Groenen PJ, Wieringa B: Somatic expansion behaviour of the (CTG)n repeat in myotonic dystrophy knock-in mice is differentially affected by Msh3 and Msh6 mismatch-repair proteins. Hum Mol Genet. 2002, 11: 191-198. 10.1093/hmg/11.2.191View ArticlePubMedGoogle Scholar
- Telenius H, Kremer B, Goldberg YP, Theilmann J, Andrew SE, Zeisler J, Adam S, Greenberg C, Ives EJ, Clarke LA, et al.: Somatic and gonadal mosaicism of the Huntington disease gene CAG repeat in brain and sperm. Nat Genet. 1994, 6: 409-414. 10.1038/ng0494-409View ArticlePubMedGoogle Scholar
- Wheeler VC, White JK, Gutekunst CA, Vrbanac V, Weaver M, Li XJ, Li SH, Yi H, Vonsattel JP, Gusella JF, et al.: Long glutamine tracts cause nuclear localization of a novel form of huntingtin in medium spiny striatal neurons in HdhQ92 and HdhQ111 knock-in mice. Hum Mol Genet. 2000, 9: 503-513. 10.1093/hmg/9.4.503View ArticlePubMedGoogle Scholar
- Gomes-Pereira M, Bidichandani SI, Monckton DG: Analysis of unstable triplet repeats using small-pool polymerase chain reaction. Methods Mol Biol. 2004, 277: 61-76.PubMedGoogle Scholar
- Mevik B-H, Cederkvist HR: Mean squared error of prediction (MSEP) estimates for principal component regression (PCR) and partial least squares regression (PLSR). Journal of Chemometrics. 2004, 18: 422-429. 10.1002/cem.887.View ArticleGoogle Scholar
- Jung J, Bonini N: CREB-binding protein modulates repeat instability in a Drosophila model for polyQ disease. Science. 2007, 315: 1857-1859. 10.1126/science.1139517View ArticlePubMedGoogle Scholar
- Cyr M, Sotnikova TD, Gainetdinov RR, Caron MG: Dopamine enhances motor and neuropathological consequences of polyglutamine expanded huntingtin. Faseb J. 2006, 20: 2541-2543. 10.1096/fj.06-6533fjeView ArticlePubMedGoogle Scholar
- Klein JA, Longo-Guess CM, Rossmann MP, Seburn KL, Hurd RE, Frankel WN, Bronson RT, Ackerman SL: The harlequin mouse mutation downregulates apoptosis-inducing factor. Nature. 2002, 419: 367-374. 10.1038/nature01034View ArticlePubMedGoogle Scholar
- Gonitel R, Moffitt H, Sathasivam K, Woodman B, Detloff PJ, Faull RL, Bates GP: DNA instability in postmitotic neurons. Proc Natl Acad Sci USA. 2008, 105: 3467-3472. 10.1073/pnas.0800048105PubMed CentralView ArticlePubMedGoogle Scholar
- Kovtun IV, Liu Y, Bjoras M, Klungland A, Wilson SH, McMurray CT: OGG1 initiates age-dependent CAG trinucleotide expansion in somatic cells. Nature. 2007, 447: 447-452. 10.1038/nature05778PubMed CentralView ArticlePubMedGoogle Scholar
- Pearson CE, Ewel A, Acharya S, Fishel RA, Sinden RR: Human MSH2 binds to trinucleotide repeat DNA structures associated with neurodegenerative diseases. Hum Mol Genet. 1997, 6: 1117-1123. 10.1093/hmg/6.7.1117View ArticlePubMedGoogle Scholar
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102: 15545-15550. 10.1073/pnas.0506580102PubMed CentralView ArticlePubMedGoogle Scholar
- Shelbourne PF, Keller-McGandy C, Bi WL, Yoon SR, Dubeau L, Veitch NJ, Vonsattel JP, Wexler NS, Arnheim N, Augood SJ: Triplet repeat mutation length gains correlate with cell-type specific vulnerability in Huntington disease brain. Hum Mol Genet. 2007, 16: 1133-1142. 10.1093/hmg/ddm054View ArticlePubMedGoogle Scholar
- Gomes-Pereira M, Fortune MT, Monckton DG: Mouse tissue culture models of unstable triplet repeats: in vitro selection for larger alleles, mutational expansion bias and tissue specificity, but no association with cell division rates. Hum Mol Genet. 2001, 10: 845-854. 10.1093/hmg/10.8.845View ArticlePubMedGoogle Scholar
- Trettel F, Rigamonti D, Hilditch-Maguire P, Wheeler VC, Sharp AH, Persichetti F, Cattaneo E, MacDonald ME: Dominant phenotypes produced by the HD mutation in STHdh(Q111) striatal cells. Hum Mol Genet. 2000, 9: 2799-2809. 10.1093/hmg/9.19.2799View ArticlePubMedGoogle Scholar
- Lin Y, Hubert L, Wilson JH: Transcription destabilizes triplet repeats. Mol Carcinog. 2009, 48: 350-361. 10.1002/mc.20488PubMed CentralView ArticlePubMedGoogle Scholar