Genetic modifiers of somatic expansion and clinical phenotypes in Huntington's disease reveal shared and tissue-specific effects.

Huntington's disease (HD), due to expansion of a CAG repeat in HTT , is representative of a growing number of disorders involving somatically unstable short tandem repeats. We find that overlapping and distinct genetic modifiers of clinical landmarks and somatic expansion in blood DNA reveal an underlying complexity and cell-type specificity to the mismatch repair-related processes that influence disease timing. Differential capture of non-DNA-repair gene modifiers by multiple measures of cognitive and motor dysfunction argues additionally for cell-type specificity of pathogenic processes. Beyond trans modifiers, differential effects are also illustrated at HTT by a 5'-UTR variant that promotes somatic expansion in blood without influencing clinical HD, while, even after correcting for uninterrupted CAG length, a synonymous sequence change at the end of the CAG repeat dramatically hastens onset of motor signs without increasing somatic expansion. Our findings are directly relevant to therapeutic suppression of somatic expansion in HD and related disorders and provide a route to define the individual neuronal cell types that contribute to different HD clinical phenotypes.

Identification of genetic modifiers of Huntington's disease somatic CAG repeat instability by in vivo CRISPR-Cas9 genome editing.

Huntington's disease (HD), one of >50 inherited repeat expansion disorders (Depienne and Mandel, 2021), is a dominantly-inherited neurodegenerative disease caused by a CAG expansion in HTT (The Huntington's Disease Collaborative Research Group, 1993). Inherited CAG repeat length is the primary determinant of age of onset, with human genetic studies underscoring that the property driving disease is the CAG length-dependent propensity of the repeat to further expand in brain (Swami et al ., 2009; GeM-HD, 2015; Hensman Moss et al ., 2017; Ciosi et al ., 2019; GeM-HD, 2019; Hong et al ., 2021). Routes to slowing somatic CAG expansion therefore hold great promise for disease-modifying therapies. Several DNA repair genes, notably in the mismatch repair (MMR) pathway, modify somatic expansion in HD mouse models (Wheeler and Dion, 2021). To identify novel modifiers of somatic expansion, we have used CRISPR-Cas9 editing in HD knock-in mice to enable in vivo screening of expansion-modifier candidates at scale. This has included testing of HD onset modifier genes emerging from human genome-wide association studies (GWAS), as well as interactions between modifier genes, thereby providing new insight into pathways underlying CAG expansion and potential therapeutic targets.

Single nuclei RNA-seq reveals a medium spiny neuron glutamate excitotoxicity signature prior to the onset of neuronal death in an ovine Huntington's disease model.

Huntington's disease (HD) is a neurodegenerative genetic disorder caused by an expansion in the CAG repeat tract of the huntingtin (HTT) gene resulting in behavioural, cognitive, and motor defects. Current knowledge of disease pathogenesis remains incomplete, and no disease course-modifying interventions are in clinical use. We have previously reported the development and characterisation of the OVT73 transgenic sheep model of HD. The 73 polyglutamine repeat is somatically stable and therefore likely captures a prodromal phase of the disease with an absence of motor symptomatology even at 5-years of age and no detectable striatal cell loss. To better understand the disease-initiating events we have undertaken a single nuclei transcriptome study of the striatum of an extensively studied cohort of 5-year-old OVT73 HD sheep and age matched wild-type controls. We have identified transcriptional upregulation of genes encoding N-methyl-D-aspartate (NMDA), α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) and kainate receptors in medium spiny neurons, the cell type preferentially lost early in HD. Further, we observed an upregulation of astrocytic glutamate uptake transporters and medium spiny neuron GABAA receptors, which may maintain glutamate homeostasis. Taken together, these observations support the glutamate excitotoxicity hypothesis as an early neurodegeneration cascade-initiating process but the threshold of toxicity may be regulated by several protective mechanisms. Addressing this biochemical defect early may prevent neuronal loss and avoid the more complex secondary consequences precipitated by cell death.

Posttranscriptional regulation of FAN1 by miR-124-3p at rs3512 underlies onset-delaying genetic modification in Huntington's disease.

Many Mendelian disorders, such as Huntington's disease (HD) and spinocerebellar ataxias, arise from expansions of CAG trinucleotide repeats. Despite the clear genetic causes, additional genetic factors may influence the rate of those monogenic disorders. Notably, genome-wide association studies discovered somewhat expected modifiers, particularly mismatch repair genes involved in the CAG repeat instability, impacting age at onset of HD. Strikingly, FAN1, previously unrelated to repeat instability, produced the strongest HD modification signals. Diverse FAN1 haplotypes independently modify HD, with rare genetic variants diminishing DNA binding or nuclease activity of the FAN1 protein, hastening HD onset. However, the mechanism behind the frequent and the most significant onset-delaying FAN1 haplotype lacking missense variations has remained elusive. Here, we illustrated that a microRNA acting on 3'-UTR (untranslated region) SNP rs3512, rather than transcriptional regulation, is responsible for the significant FAN1 expression quantitative trait loci signal and allelic imbalance in FAN1 messenger ribonucleic acid (mRNA), accounting for the most significant and frequent onset-delaying modifier haplotype in HD. Specifically, miR-124-3p selectively targets the reference allele at rs3512, diminishing the stability of FAN1 mRNA harboring that allele and consequently reducing its levels. Subsequent validation analyses, including the use of antagomir and 3'-UTR reporter vectors with swapped alleles, confirmed the specificity of miR-124-3p at rs3512. Together, these findings indicate that the alternative allele at rs3512 renders the FAN1 mRNA less susceptible to miR-124-3p-mediated posttranscriptional regulation, resulting in increased FAN1 levels and a subsequent delay in HD onset by mitigating CAG repeat instability.

Splice modulators target PMS1 to reduce somatic expansion of the Huntington's disease-associated CAG repeat.

Huntington's disease (HD) is a dominant neurological disorder caused by an expanded HTT exon 1 CAG repeat that lengthens huntingtin's polyglutamine tract. Lowering mutant huntingtin has been proposed for treating HD, but genetic modifiers implicate somatic CAG repeat expansion as the driver of onset. We find that branaplam and risdiplam, small molecule splice modulators that lower huntingtin by promoting HTT pseudoexon inclusion, also decrease expansion of an unstable HTT exon 1 CAG repeat in an engineered cell model. Targeted CRISPR-Cas9 editing shows this effect is not due to huntingtin lowering, pointing instead to pseudoexon inclusion in PMS1. Homozygous but not heterozygous inactivation of PMS1 also reduces CAG repeat expansion, supporting PMS1 as a genetic modifier of HD and a potential target for therapeutic intervention. Although splice modulation provides one strategy, genome-wide transcriptomics also emphasize consideration of cell-type specific effects and polymorphic variation at both target and off-target sites.

Modification of Huntington's disease by short tandem repeats.

Expansions of glutamine-coding CAG trinucleotide repeats cause a number of neurodegenerative diseases, including Huntington's disease and several of spinocerebellar ataxias. In general, age-at-onset of the polyglutamine diseases is inversely correlated with the size of the respective inherited expanded CAG repeat. Expanded CAG repeats are also somatically unstable in certain tissues, and age-at-onset of Huntington's disease corrected for individual HTT CAG repeat length (i.e. residual age-at-onset), is modified by repeat instability-related DNA maintenance/repair genes as demonstrated by recent genome-wide association studies. Modification of one polyglutamine disease (e.g. Huntington's disease) by the repeat length of another (e.g. ATXN3, CAG expansions in which cause spinocerebellar ataxia 3) has also been hypothesized. Consequently, we determined whether age-at-onset in Huntington's disease is modified by the CAG repeats of other polyglutamine disease genes. We found that the CAG measured repeat sizes of other polyglutamine disease genes that were polymorphic in Huntington's disease participants but did not influence Huntington's disease age-at-onset. Additional analysis focusing specifically on ATXN3 in a larger sample set (n = 1388) confirmed the lack of association between Huntington's disease residual age-at-onset and ATXN3 CAG repeat length. Additionally, neither our Huntington's disease onset modifier genome-wide association studies single nucleotide polymorphism data nor imputed short tandem repeat data supported the involvement of other polyglutamine disease genes in modifying Huntington's disease. By contrast, our genome-wide association studies based on imputed short tandem repeats revealed significant modification signals for other genomic regions. Together, our short tandem repeat genome-wide association studies show that modification of Huntington's disease is associated with short tandem repeats that do not involve other polyglutamine disease-causing genes, refining the landscape of Huntington's disease modification and highlighting the importance of rigorous data analysis, especially in genetic studies testing candidate modifiers.

Parallelized engineering of mutational models using piggyBac transposon delivery of CRISPR libraries.

New technologies and large-cohort studies have enabled novel variant discovery and association at unprecedented scale, yet functional characterization of these variants remains paramount to deciphering disease mechanisms. Approaches that facilitate parallelized genome editing of cells of interest or induced pluripotent stem cells (iPSCs) have become critical tools toward this goal. Here, we developed an approach that incorporates libraries of CRISPR-Cas9 guide RNAs (gRNAs) together with inducible Cas9 into a piggyBac (PB) transposon system to engineer dozens to hundreds of genomic variants in parallel against isogenic cellular backgrounds. This method empowers loss-of-function (LoF) studies through the introduction of insertions or deletions (indels) and copy-number variants (CNVs), though generating specific nucleotide changes is possible with prime editing. The ability to rapidly establish high-quality mutational models at scale will facilitate the development of isogenic cellular collections and catalyze comparative functional genomic studies investigating the roles of hundreds of genes and mutations in development and disease.

PMS1 as a target for splice modulation to prevent somatic CAG repeat expansion in Huntington's disease.

Huntington's disease (HD) is a dominantly inherited neurodegenerative disorder whose motor, cognitive, and behavioral manifestations are caused by an expanded, somatically unstable CAG repeat in the first exon of HTT that lengthens a polyglutamine tract in huntingtin. Genome-wide association studies (GWAS) have revealed DNA repair genes that influence the age-at-onset of HD and implicate somatic CAG repeat expansion as the primary driver of disease timing. To prevent the consequent neuronal damage, small molecule splice modulators (e.g., branaplam) that target HTT to reduce the levels of huntingtin are being investigated as potential HD therapeutics. We found that the effectiveness of the splice modulators can be influenced by genetic variants, both at HTT and other genes where they promote pseudoexon inclusion. Surprisingly, in a novel hTERT-immortalized retinal pigment epithelial cell (RPE1) model for assessing CAG repeat instability, these drugs also reduced the rate of HTT CAG expansion. We determined that the splice modulators also affect the expression of the mismatch repair gene PMS1, a known modifier of HD age-at-onset. Genome editing at specific HTT and PMS1 sequences using CRISPR-Cas9 nuclease confirmed that branaplam suppresses CAG expansion by promoting the inclusion of a pseudoexon in PMS1, making splice modulation of PMS1 a potential strategy for delaying HD onset. Comparison with another splice modulator, risdiplam, suggests that other genes affected by these splice modulators also influence CAG instability and might provide additional therapeutic targets.

A Cdk5-derived peptide inhibits Cdk5/p25 activity and improves neurodegenerative phenotypes.

Aberrant activity of cyclin-dependent kinase (Cdk5) has been implicated in various neurodegenerative diseases. This deleterious effect is mediated by pathological cleavage of the Cdk5 activator p35 into the truncated product p25, leading to prolonged Cdk5 activation and altered substrate specificity. Elevated p25 levels have been reported in humans and rodents with neurodegeneration, and the benefit of genetically blocking p25 production has been demonstrated previously in rodent and human neurodegenerative models. Here, we report a 12-amino-acid-long peptide fragment derived from Cdk5 (Cdk5i) that is considerably smaller than existing peptide inhibitors of Cdk5 (P5 and CIP) but shows high binding affinity toward the Cdk5/p25 complex, disrupts the interaction of Cdk5 with p25, and lowers Cdk5/p25 kinase activity. When tagged with a fluorophore (FITC) and the cell-penetrating transactivator of transcription (TAT) sequence, the Cdk5i-FT peptide exhibits cell- and brain-penetrant properties and confers protection against neurodegenerative phenotypes associated with Cdk5 hyperactivity in cell and mouse models of neurodegeneration, highlighting Cdk5i's therapeutic potential.

Haplotype-specific MAPK3 expression in 16p11.2 deletion contributes to variable neurodevelopment.

Recurrent proximal 16p11.2 deletion (16p11.2del) is a risk factor for diverse neurodevelopmental disorders with incomplete penetrance and variable expressivity. Although investigation with human induced pluripotent stem cell models has confirmed disruption of neuronal development in 16p11.2del neuronal cells, which genes are responsible for abnormal cellular phenotypes and what determines the penetrance of neurodevelopmental abnormalities are unknown. We performed haplotype phasing of the 16p11.2 region in a 16p11.2del neurodevelopmental disorders cohort and generated human induced pluripotent stem cells for two 16p11.2del families with distinct residual haplotypes and variable neurodevelopmental disorder phenotypes. Using transcriptomic profiles and cellular phenotypes of the human induced pluripotent stem cell-differentiated cortex neuronal cells, we revealed MAPK3 to be a contributor to dysfunction in multiple pathways related to early neuronal development, with altered soma and electrophysiological properties in mature neuronal cells. Notably, MAPK3 expression in 16p11.2del neuronal cells varied on the basis of a 132 kb 58 single nucleotide polymorphism (SNP) residual haplotype, with the version composed entirely of minor alleles associated with reduced MAPK3 expression. Ten SNPs on the residual haplotype were mapped to enhancers of MAPK3. We functionally validated six of these SNPs by luciferase assay, implicating them in the residual haplotype-specific differences in MAPK3 expression via cis-regulation. Finally, the analysis of three different cohorts of 16p11.2del subjects showed that this minor residual haplotype is associated with neurodevelopmental disorder phenotypes in 16p11.2del carriers.

CHD8 suppression impacts on histone H3 lysine 36 trimethylation and alters RNA alternative splicing.

Disruptive mutations in the chromodomain helicase DNA-binding protein 8 gene (CHD8) have been recurrently associated with autism spectrum disorders (ASDs). Here we investigated how chromatin reacts to CHD8 suppression by analyzing a panel of histone modifications in induced pluripotent stem cell-derived neural progenitors. CHD8 suppression led to significant reduction (47.82%) in histone H3K36me3 peaks at gene bodies, particularly impacting on transcriptional elongation chromatin states. H3K36me3 reduction specifically affects highly expressed, CHD8-bound genes and correlates with altered alternative splicing patterns of 462 genes implicated in 'regulation of RNA splicing' and 'mRNA catabolic process'. Mass spectrometry analysis uncovered a novel interaction between CHD8 and the splicing regulator heterogeneous nuclear ribonucleoprotein L (hnRNPL), providing the first mechanistic insights to explain the CHD8 suppression-derived splicing phenotype, partly implicating SETD2, a H3K36me3 methyltransferase. In summary, our results point toward broad molecular consequences of CHD8 suppression, entailing altered histone deposition/maintenance and RNA processing regulation as important regulatory processes in ASD.

Transcriptional and functional consequences of alterations to MEF2C and its topological organization in neuronal models.

Point mutations and structural variants that directly disrupt the coding sequence of MEF2C have been associated with a spectrum of neurodevelopmental disorders (NDDs). However, the impact of MEF2C haploinsufficiency on neurodevelopmental pathways and synaptic processes is not well understood, nor are the complex mechanisms that govern its regulation. To explore the functional changes associated with structural variants that alter MEF2C expression and/or regulation, we generated an allelic series of 204 isogenic human induced pluripotent stem cell (hiPSC)-derived neural stem cells and glutamatergic induced neurons. These neuronal models harbored CRISPR-engineered mutations that involved direct deletion of MEF2C or deletion of the boundary points for topologically associating domains (TADs) and chromatin loops encompassing MEF2C. Systematic profiling of mutation-specific alterations, contrasted to unedited controls that were exposed to the same guide RNAs for each edit, revealed that deletion of MEF2C caused differential expression of genes associated with neurodevelopmental pathways and synaptic function. We also discovered significant reduction in synaptic activity measured by multielectrode arrays (MEAs) in neuronal cells. By contrast, we observed robust buffering against MEF2C regulatory disruption following deletion of a distal 5q14.3 TAD and loop boundary, whereas homozygous loss of a proximal loop boundary resulted in down-regulation of MEF2C expression and reduced electrophysiological activity on MEA that was comparable to direct gene disruption. Collectively, these studies highlight the considerable functional impact of MEF2C deletion in neuronal cells and systematically characterize the complex interactions that challenge a priori predictions of regulatory consequences from structural variants that disrupt three-dimensional genome organization.

Huntington's disease age at motor onset is modified by the tandem hexamer repeat in TCERG1.

Huntington's disease is caused by an expanded CAG tract in HTT. The length of the CAG tract accounts for over half the variance in age at onset of disease, and is influenced by other genetic factors, mostly implicating the DNA maintenance machinery. We examined a single nucleotide variant, rs79727797, on chromosome 5 in the TCERG1 gene, previously reported to be associated with Huntington's disease and a quasi-tandem repeat (QTR) hexamer in exon 4 of TCERG1 with a central pure repeat. We developed a method for calling perfect and imperfect repeats from exome-sequencing data, and tested association between the QTR in TCERG1 and residual age at motor onset (after correcting for the effects of CAG length in the HTT gene) in 610 individuals with Huntington's disease via regression analysis. We found a significant association between age at onset and the sum of the repeat lengths from both alleles of the QTR (p = 2.1 × 10-9), with each added repeat hexamer reducing age at onset by one year (95% confidence interval [0.7, 1.4]). This association explained that previously observed with rs79727797. The association with age at onset in the genome-wide association study is due to a QTR hexamer in TCERG1, translated to a glutamine/alanine tract in the protein. We could not distinguish whether this was due to cis-effects of the hexamer repeat on gene expression or of the encoded glutamine/alanine tract in the protein. These results motivate further study of the mechanisms by which TCERG1 modifies onset of HD.

Tissue- and cell-type-specific molecular and functional signatures of 16p11.2 reciprocal genomic disorder across mouse brain and human neuronal models.

Chromosome 16p11.2 reciprocal genomic disorder, resulting from recurrent copy-number variants (CNVs), involves intellectual disability, autism spectrum disorder (ASD), and schizophrenia, but the responsible mechanisms are not known. To systemically dissect molecular effects, we performed transcriptome profiling of 350 libraries from six tissues (cortex, cerebellum, striatum, liver, brown fat, and white fat) in mouse models harboring CNVs of the syntenic 7qF3 region, as well as cellular, transcriptional, and single-cell analyses in 54 isogenic neural stem cell, induced neuron, and cerebral organoid models of CRISPR-engineered 16p11.2 CNVs. Transcriptome-wide differentially expressed genes were largely tissue-, cell-type-, and dosage-specific, although more effects were shared between deletion and duplication and across tissue than expected by chance. The broadest effects were observed in the cerebellum (2,163 differentially expressed genes), and the greatest enrichments were associated with synaptic pathways in mouse cerebellum and human induced neurons. Pathway and co-expression analyses identified energy and RNA metabolism as shared processes and enrichment for ASD-associated, loss-of-function constraint, and fragile X messenger ribonucleoprotein target gene sets. Intriguingly, reciprocal 16p11.2 dosage changes resulted in consistent decrements in neurite and electrophysiological features, and single-cell profiling of organoids showed reciprocal alterations to the proportions of excitatory and inhibitory GABAergic neurons. Changes both in neuronal ratios and in gene expression in our organoid analyses point most directly to calretinin GABAergic inhibitory neurons and the excitatory/inhibitory balance as targets of disruption that might contribute to changes in neurodevelopmental and cognitive function in 16p11.2 carriers. Collectively, our data indicate the genomic disorder involves disruption of multiple contributing biological processes and that this disruption has relative impacts that are context specific.

A cross-disorder dosage sensitivity map of the human genome.

Rare copy-number variants (rCNVs) include deletions and duplications that occur infrequently in the global human population and can confer substantial risk for disease. In this study, we aimed to quantify the properties of haploinsufficiency (i.e., deletion intolerance) and triplosensitivity (i.e., duplication intolerance) throughout the human genome. We harmonized and meta-analyzed rCNVs from nearly one million individuals to construct a genome-wide catalog of dosage sensitivity across 54 disorders, which defined 163 dosage sensitive segments associated with at least one disorder. These segments were typically gene dense and often harbored dominant dosage sensitive driver genes, which we were able to prioritize using statistical fine-mapping. Finally, we designed an ensemble machine-learning model to predict probabilities of dosage sensitivity (pHaplo & pTriplo) for all autosomal genes, which identified 2,987 haploinsufficient and 1,559 triplosensitive genes, including 648 that were uniquely triplosensitive. This dosage sensitivity resource will provide broad utility for human disease research and clinical genetics.

Orgo-Seq integrates single-cell and bulk transcriptomic data to identify cell type specific-driver genes associated with autism spectrum disorder.

Cerebral organoids can be used to gain insights into cell type specific processes perturbed by genetic variants associated with neuropsychiatric disorders. However, robust and scalable phenotyping of organoids remains challenging. Here, we perform RNA sequencing on 71 samples comprising 1,420 cerebral organoids from 25 donors, and describe a framework (Orgo-Seq) to integrate bulk RNA and single-cell RNA sequence data. We apply Orgo-Seq to 16p11.2 deletions and 15q11-13 duplications, two loci associated with autism spectrum disorder, to identify immature neurons and intermediate progenitor cells as critical cell types for 16p11.2 deletions. We further applied Orgo-Seq to identify cell type-specific driver genes. Our work presents a quantitative phenotyping framework to integrate multi-transcriptomic datasets for the identification of cell types and cell type-specific co-expressed driver genes associated with neuropsychiatric disorders.

Exome sequencing of individuals with Huntington's disease implicates FAN1 nuclease activity in slowing CAG expansion and disease onset.

The age at onset of motor symptoms in Huntington's disease (HD) is driven by HTT CAG repeat length but modified by other genes. In this study, we used exome sequencing of 683 patients with HD with extremes of onset or phenotype relative to CAG length to identify rare variants associated with clinical effect. We discovered damaging coding variants in candidate modifier genes identified in previous genome-wide association studies associated with altered HD onset or severity. Variants in FAN1 clustered in its DNA-binding and nuclease domains and were associated predominantly with earlier-onset HD. Nuclease activities of purified variants in vitro correlated with residual age at motor onset of HD. Mutating endogenous FAN1 to a nuclease-inactive form in an induced pluripotent stem cell model of HD led to rates of CAG expansion similar to those observed with complete FAN1 knockout. Together, these data implicate FAN1 nuclease activity in slowing somatic repeat expansion and hence onset of HD.

Genetic modifiers of Huntington disease differentially influence motor and cognitive domains.

Genome-wide association studies (GWASs) of Huntington 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 by using common motor and cognitive measures to predict age at other disease landmarks as additional phenotypes for GWASs. Combined with imputation with 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 additional modifier effects at the PMS1 and PMS2 loci and implicated a potential second 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.