HD CAGnome: A Search Tool for Huntingtin CAG Repeat Length-Correlated Genes

Background: The length of the huntingtin (HTT) CAG repeat is strongly correlated with both age at onset of Huntington’s disease (HD) symptoms and age at death of HD patients. Dichotomous analysis comparing HD to controls is widely used to study the effects of HTT CAG repeat expansion. However, a potentially more powerful approach is a continuous analysis strategy that takes advantage of all of the different CAG lengths, to capture effects that are expected to be critical to HD pathogenesis. Methodology/Principal Findings We used continuous and dichotomous approaches to analyze microarray gene expression data from 107 human control and HD lymphoblastoid cell lines. Of all probes found to be significant in a continuous analysis by CAG length, only 21.4% were so identified by a dichotomous comparison of HD versus controls. Moreover, of probes significant by dichotomous analysis, only 33.2% were also significant in the continuous analysis. Simulations revealed that the dichotomous approach would require substantially more than 107 samples to either detect 80% of the CAG-length correlated changes revealed by continuous analysis or to reduce the rate of significant differences that are not CAG length-correlated to 20% (n = 133 or n = 206, respectively). Given the superior power of the continuous approach, we calculated the correlation structure between HTT CAG repeat lengths and gene expression levels and created a freely available searchable website, “HD CAGnome,” that allows users to examine continuous relationships between HTT CAG and expression levels of ∼20,000 human genes. Conclusions/Significance: Our results reveal limitations of dichotomous approaches compared to the power of continuous analysis to study a disease where human genotype-phenotype relationships strongly support a role for a continuum of CAG length-dependent changes. The compendium of HTT CAG length-gene expression level relationships found at the HD CAGnome now provides convenient routes for discovery of candidates influenced by the HD mutation

A modifier of Huntington's disease onset at the MLH1 locus

Huntington’s disease (HD) is a dominantly inherited neurodegenerative disease caused by an expanded CAG repeat in HTT. Many clinical characteristics of HD such as age at motor onset are determined largely by the size of HTT CAG repeat. However, emerging evidence strongly supports a role for other genetic factors in modifying the disease pathogenesis driven by mutant huntingtin. A recent genome-wide association analysis to discover genetic modifiers of HD onset age provided initial evidence for modifier loci on chromosomes 8 and 15 and suggestive evidence for a locus on chromosome 3. Here, genotyping of candidate single nucleotide polymorphisms in a cohort of 3,314 additional HD subjects yields independent confirmation of the former two loci and moves the third to genome-wide significance at MLH1, a locus whose mouse orthologue modifies CAG length-dependent phenotypes in a Htt-knock-in mouse model of HD. Both quantitative and dichotomous association analyses implicate a functional variant on 32% of chromosomes with the beneficial modifier effect that delays HD motor onset by 0.7 years/allele. Genomic DNA capture and sequencing of a modifier haplotype localize the functional variation to a 78 kb region spanning the 3’end of MLH1 and the 5’end of the neighboring LRRFIP2, and marked by an isoleucinevaline missense variant in MLH1. Analysis of expression Quantitative Trait Loci (eQTLs) provides modest support for altered regulation of MLH1 and LRRFIP2, raising the possibility that the modifier affects regulation of both genes. Finally, polygenic modification score and heritability analyses suggest the existence of additional genetic modifiers, supporting expanded, comprehensive genetic analysis of larger HD datasets

Huntington's disease pathogenesis: two sequential components

Historically, Huntington’s disease (HD; OMIM #143100) has played an important role in the enormous advances in human genetics seen over the past four decades. This familial neurodegenerative disorder involves variable onset followed by consistent worsening of characteristic abnormal movements along with cognitive decline and psychiatric disturbances. HD was the first autosomal disease for which the genetic defect was assigned to a position on the human chromosomes using only genetic linkage analysis with common DNA polymorphisms. This discovery set off a multitude of similar studies in other diseases, while the HD gene, later renamed HTT, and its vicinity in chromosome 4p16.3 then acted as a proving ground for development of technologies to clone and sequence genes based upon their genomic location, with the growing momentum of such advances fueling the Human Genome Project. The identification of the HD gene has not yet led to an effective treatment, but continued human genetic analysis of genotype-phenotype relationships in large HD subject populations, first at the HTT locus and subsequently genome-wide, has provided insights into pathogenesis that divide the course of the disease into two sequential, mechanistically distinct components

Quantification of Age-Dependent Somatic CAG Repeat Instability in Hdh CAG Knock-In Mice Reveals Different Expansion Dynamics in Striatum and Liver

Age at onset of Huntington's disease (HD) is largely determined by the CAG trinucleotide repeat length in the HTT gene. Importantly, the CAG repeat undergoes tissue-specific somatic instability, prevalent in brain regions that are disease targets, suggesting a potential role for somatic CAG repeat instability in modifying HD pathogenesis. Thus, understanding underlying mechanisms of somatic CAG repeat instability may lead to discoveries of novel therapeutics for HD. Investigation of the dynamics of the CAG repeat size changes over time may provide insights into the mechanisms underlying CAG repeat instability.To understand how the HTT CAG repeat length changes over time, we quantified somatic instability of the CAG repeat in Huntington's disease CAG knock-in mice from 2-16 months of age in liver, striatum, spleen and tail. The HTT CAG repeat in spleen and tail was very stable, but that in liver and striatum expanded over time at an average rate of one CAG per month. Interestingly, the patterns of repeat instability were different between liver and striatum. Unstable CAG repeats in liver repeatedly gained similar sizes of additional CAG repeats (approximately two CAGs per month), maintaining a distinct population of unstable repeats. In contrast, unstable CAG repeats in striatum gained additional repeats with different sizes resulting in broadly distributed unstable CAG repeats. Expanded CAG repeats in the liver were highly enriched in polyploid hepatocytes, suggesting that the pattern of liver instability may reflect the restriction of the unstable repeats to a unique cell type.Our results are consistent with repeat expansion occurring as a consequence of recurrent small repeat insertions that differ in different tissues. Investigation of the specific mechanisms that underlie liver and striatal instability will contribute to our understanding of the relationship between instability and disease and the means to intervene in this process

Haplotype-based stratification of Huntington's disease

Huntington’s disease (HD) is an autosomal dominant neurodegenerative disease caused by expansion of a CAG trinucleotide repeat in HTT, resulting in an extended polyglutamine tract in huntingtin. We and others have previously determined that the HD-causing expansion occurs on multiple different haplotype backbones, reflecting more than one ancestral origin of the same type of mutation. In view of the therapeutic potential of mutant allele-specific gene silencing, we have compared and integrated two major systems of HTT haplotype definition, combining data from 74 sequence variants to identify the most frequent disease-associated and control chromosome backbones and revealing that there is potential for additional resolution of HD haplotypes. We have used the large collection of 4078 heterozygous HD subjects analyzed in our recent genome-wide association study of HD age at onset to estimate the frequency of these haplotypes in European subjects, finding that common genetic variation at HTT can distinguish the normal and CAG-expanded chromosomes for more than 95% of European HD individuals. As a resource for the HD research community, we have also determined the haplotypes present in a series of publicly available HD subject-derived fibroblasts, induced pluripotent cells, and embryonic stem cells in order to facilitate efforts to develop inclusive methods of allele-specific HTT silencing applicable to most HD patients. Our data providing genetic guidance for therapeutic gene-based targeting will significantly contribute to the developments of rational treatments and implementation of precision medicine in HD

Population-specific genetic modification of Huntington\u27s disease in Venezuela.

Modifiers of Mendelian disorders can provide insights into disease mechanisms and guide therapeutic strategies. A recent genome-wide association (GWA) study discovered genetic modifiers of Huntington\u27s disease (HD) onset in Europeans. Here, we performed whole genome sequencing and GWA analysis of a Venezuelan HD cluster whose families were crucial for the original mapping of the HD gene defect. The Venezuelan HD subjects develop motor symptoms earlier than their European counterparts, implying the potential for population-specific modifiers. The main Venezuelan HD family inherits HTT haplotype hap.03, which differs subtly at the sequence level from European HD hap.03, suggesting a different ancestral origin but not explaining the earlier age at onset in these Venezuelans. GWA analysis of the Venezuelan HD cluster suggests both population-specific and population-shared genetic modifiers. Genome-wide significant signals at 7p21.2-21.1 and suggestive association signals at 4p14 and 17q21.2 are evident only in Venezuelan HD, but genome-wide significant association signals at the established European chromosome 15 modifier locus are improved when Venezuelan HD data are included in the meta-analysis. Venezuelan-specific association signals on chromosome 7 center on SOSTDC1, which encodes a bone morphogenetic protein antagonist. The corresponding SNPs are associated with reduced expression of SOSTDC1 in non-Venezuelan tissue samples, suggesting that interaction of reduced SOSTDC1 expression with a population-specific genetic or environmental factor may be responsible for modification of HD onset in Venezuela. Detection of population-specific modification in Venezuelan HD supports the value of distinct disease populations in revealing novel aspects of a disease and population-relevant therapeutic strategies

Haplotype-based stratification of Huntington's disease

Huntington’s disease (HD) is an autosomal dominant neurodegenerative disease caused by expansion of a CAG trinucleotide repeat in HTT, resulting in an extended polyglutamine tract in huntingtin. We and others have previously determined that the HD-causing expansion occurs on multiple different haplotype backbones, reflecting more than one ancestral origin of the same type of mutation. In view of the therapeutic potential of mutant allele-specific gene silencing, we have compared and integrated two major systems of HTT haplotype definition, combining data from 74 sequence variants to identify the most frequent disease-associated and control chromosome backbones and revealing that there is potential for additional resolution of HD haplotypes. We have used the large collection of 4078 heterozygous HD subjects analyzed in our recent genome-wide association study of HD age at onset to estimate the frequency of these haplotypes in European subjects, finding that common genetic variation at HTT can distinguish the normal and CAG-expanded chromosomes for more than 95% of European HD individuals. As a resource for the HD research community, we have also determined the haplotypes present in a series of publicly available HD subject-derived fibroblasts, induced pluripotent cells, and embryonic stem cells in order to facilitate efforts to develop inclusive methods of allele-specific HTT silencing applicable to most HD patients. Our data providing genetic guidance for therapeutic gene-based targeting will significantly contribute to the developments of rational treatments and implementation of precision medicine in HD

Mismatch Repair Genes Mlh1 and Mlh3 Modify CAG Instability in Huntington's Disease Mice: Genome-Wide and Candidate Approaches

The Huntington's disease gene (HTT) CAG repeat mutation undergoes somatic expansion that correlates with pathogenesis. Modifiers of somatic expansion may therefore provide routes for therapies targeting the underlying mutation, an approach that is likely applicable to other trinucleotide repeat diseases. Huntington's disease HdhQ111 mice exhibit higher levels of somatic HTT CAG expansion on a C57BL/6 genetic background (B6.HdhQ111) than on a 129 background (129.HdhQ111). Linkage mapping in (B6x129).HdhQ111 F2 intercross animals identified a single quantitative trait locus underlying the strain-specific difference in expansion in the striatum, implicating mismatch repair (MMR) gene Mlh1 as the most likely candidate modifier. Crossing B6.HdhQ111 mice onto an Mlh1 null background demonstrated that Mlh1 is essential for somatic CAG expansions and that it is an enhancer of nuclear huntingtin accumulation in striatal neurons. HdhQ111 somatic expansion was also abolished in mice deficient in the Mlh3 gene, implicating MutLγ (MLH1–MLH3) complex as a key driver of somatic expansion. Strikingly, Mlh1 and Mlh3 genes encoding MMR effector proteins were as critical to somatic expansion as Msh2 and Msh3 genes encoding DNA mismatch recognition complex MutSβ (MSH2–MSH3). The Mlh1 locus is highly polymorphic between B6 and 129 strains. While we were unable to detect any difference in base-base mismatch or short slipped-repeat repair activity between B6 and 129 MLH1 variants, repair efficiency was MLH1 dose-dependent. MLH1 mRNA and protein levels were significantly decreased in 129 mice compared to B6 mice, consistent with a dose-sensitive MLH1-dependent DNA repair mechanism underlying the somatic expansion difference between these strains. Together, these data identify Mlh1 and Mlh3 as novel critical genetic modifiers of HTT CAG instability, point to Mlh1 genetic variation as the likely source of the instability difference in B6 and 129 strains and suggest that MLH1 protein levels play an important role in driving of the efficiency of somatic expansions

Mismatch Repair Genes Mlh1 and Mlh3 Modify CAG Instability in Huntington\u27s Disease Mice: Genome-Wide and Candidate Approaches

The Huntington\u27s disease gene (HTT) CAG repeat mutation undergoes somatic expansion that correlates with pathogenesis. Modifiers of somatic expansion may therefore provide routes for therapies targeting the underlying mutation, an approach that is likely applicable to other trinucleotide repeat diseases. Huntington\u27s disease Hdh(Q111) mice exhibit higher levels of somatic HTT CAG expansion on a C57BL/6 genetic background (B6.Hdh(Q111) ) than on a 129 background (129.Hdh(Q111) ). Linkage mapping in (B6x129).Hdh(Q111) F2 intercross animals identified a single quantitative trait locus underlying the strain-specific difference in expansion in the striatum, implicating mismatch repair (MMR) gene Mlh1 as the most likely candidate modifier. Crossing B6.Hdh(Q111) mice onto an Mlh1 null background demonstrated that Mlh1 is essential for somatic CAG expansions and that it is an enhancer of nuclear huntingtin accumulation in striatal neurons. Hdh(Q111) somatic expansion was also abolished in mice deficient in the Mlh3 gene, implicating MutLγ (MLH1-MLH3) complex as a key driver of somatic expansion. Strikingly, Mlh1 and Mlh3 genes encoding MMR effector proteins were as critical to somatic expansion as Msh2 and Msh3 genes encoding DNA mismatch recognition complex MutSβ (MSH2-MSH3). The Mlh1 locus is highly polymorphic between B6 and 129 strains. While we were unable to detect any difference in base-base mismatch or short slipped-repeat repair activity between B6 and 129 MLH1 variants, repair efficiency was MLH1 dose-dependent. MLH1 mRNA and protein levels were significantly decreased in 129 mice compared to B6 mice, consistent with a dose-sensitive MLH1-dependent DNA repair mechanism underlying the somatic expansion difference between these strains. Together, these data identify Mlh1 and Mlh3 as novel critical genetic modifiers of HTT CAG instability, point to Mlh1 genetic variation as the likely source of the instability difference in B6 and 129 strains and suggest that MLH1 protein levels play an important role in driving of the efficiency of somatic expansions