Search

Targeted DNA sequencing approaches will improve how the size of short tandem repeats is measured for diagnostic tests and preclinical studies. The expansion of these sequences causes dozens of disorders, with longer tracts generally leading to a more severe disease. Interrupted alleles are sometimes present within repeats and can alter disease manifestation. Determining repeat size mosaicism and identifying interruptions in targeted sequencing datasets remains a major challenge. This is in part because standard alignment tools are ill-suited for repetitive and unstable sequences. To address this, we have developed Repeat Detector (RD), a deterministic profile weighting algorithm for counting repeats in targeted sequencing data. We tested RD using blood-derived DNA samples from Huntington’s disease and Fuchs endothelial corneal dystrophy patients sequenced using either Illumina MiSeq or Pacific Biosciences single-molecule, real-time sequencing platforms. RD was highly accurate in determining repeat sizes of 609 blood-derived samples from Huntington’s disease individuals and did not require prior knowledge of the flanking sequences. Furthermore, RD can be used to identify alleles with interruptions and provide a measure of repeat instability within an individual. RD is therefore highly versatile and may find applications in the diagnosis of expanded repeat disorders and in the development of novel therapies.

We aimed to determine the genetic diversity and molecular characteristics of the Huntington disease (HD) gene (HTT) in Spain. We performed an extended haplotype and exon one deep sequencing analysis of the HTT gene in a nationwide cohort of population-based controls (n = 520) and families with symptomatic individuals referred for HD genetic testing. This group included 331 HD cases and 140 carriers of intermediate alleles. Clinical and family history data were obtained when available. Spanish normal alleles are enriched in C haplotypes (40.1%), whereas A1 (39.8%) and A2 (31.6%) prevail among intermediate and expanded alleles, respectively. Alleles ≥ 50 CAG repeats are primarily associated with haplotypes A2 (38.9%) and C (32%), which are also present in 50% and 21.4%, respectively, of HD families with large intergenerational expansions. Non-canonical variants of exon one sequence are less frequent, but much more diverse, in alleles of ≥27 CAG repeats. The deletion of CAACAG, one of the six rare variants not observed among smaller normal alleles, is associated with haplotype C and appears to correlate with larger intergenerational expansions and early onset of symptoms. Spanish HD haplotypes are characterized by a high genetic diversity, potentially admixed with other non-Caucasian populations, with a higher representation of A2 and C haplotypes than most European populations. Differences in haplotype distributions across the CAG length range support differential germline expansion dynamics, with A2 and C showing the largest intergenerational expansions. This haplotype-dependent germline instability may be driven by specific cis-elements, such as the CAACAG deletion. This study was funded by the Instituto de Salud Carlos III (grant no. FIS PI15/02227), and the CHDI Foundation (Enroll -HD study). A. Ruiz de Sabando was supported by a Fellowship of Departamento de Salud, Gobierno de Navarra (ref. 0011-4809-2018-000001).

Genome-wide association studies (GWAS) of Huntington’s disease (HD) have identified six DNA maintenance gene loci (among others) as modifiers and implicated a two step-mechanism of pathogenesis: somatic instability of the causative HTT CAG repeat with subsequent triggering of neuronal damage. The largest studies have been limited to HD individuals with a rater-estimated age at motor onset. To capitalize on the wealth of phenotypic data in several large HD natural history studies, we have performed algorithmic prediction using common motor and cognitive measures to predict age at other disease landmarks as additional phenotypes for GWAS. Combined with imputation using the Trans-Omics for Precision Medicine reference panel, predictions using integrated measures provided objective landmark phenotypes with greater power to detect most modifier loci. Importantly, substantial differences in the relative modifier signal across loci, highlighted by comparing common modifiers at MSH3 and FAN1, revealed that individual modifier effects can act preferentially in the motor or cognitive domains. Individual components of the DNA maintenance modifier mechanisms may therefore act differentially on the neuronal circuits underlying the corresponding clinical measures. In addition, we identified new modifier effects at the PMS1 and PMS2 loci and implicated a potential new locus on chromosome 7. These findings indicate that broadened discovery and characterization of HD genetic modifiers based on additional quantitative or qualitative phenotypes offers not only the promise of in-human validated therapeutic targets, but also a route to dissecting the mechanisms and cell types involved in both the somatic instability and toxicity components of HD pathogenesis.