One of the seminal discoveries from genetic studies of autism spectrum disorder and related neurodevelopmental disorders (NDDs) has been that loss-of-function (LoF) mutations in genes that impact transcriptional regulation confer substantial liability to NDDs. Haploinsufficiency of the epigenetic regulator POGZ represents one of the strongest such associations; however, little is known about the mechanisms by which POGZ LoF alters early neuronal development. Here, we created an allelic series of CRISPR-engineered human induced pluripotent stem cell (hiPSC) clones harboring mono- and bi-allelic POGZ deletions. In hiPSC-derived neural stem cells (NSCs) and Neurogenin-2-induced neurons (iNs), POGZ LoF altered the expression of genes associated with synaptic and intracellular signaling and extracellular matrix organization. Our multiomics profiling also showed altered footprinting of critical transcription factors (e.g., activator protein 1 complexes) that were enriched at promoters of differentially expressed genes associated with synaptic function. To further interrogate the shared molecular changes associated with NDDs, we compared our results to deletions of the transcription factor MEF2C and the sodium channel gene SCN2A that we generated in these same isogenic iNs. These analyses revealed strong enrichment of extracellular matrix and intracellular signaling disruption associated with POGZ and MEF2C deletion, whereas POGZ and SCN2A haploinsufficiency exhibited shared transcriptional effects on gene modules enriched for NDD-associated genes with opposing regulatory effects. Notably, we also observed alterations to synaptic firing rate and neurite extension with bi-allelic deletions. These shared molecular consequences suggest key points of convergence that connect gene regulation to neuronal function in the etiology of neurodevelopmental pathologies.

Cohesin orchestrates gene expression via three-dimensional chromosome folding. Genes encoding cohesin and cohesin loaders have been associated with Mendelian disorders, whereas genes encoding cohesin release factors, including WAPL and its binding partners PDS5A and PDS5B, have not. We explored the relevance of cohesin release factors in Mendelian disease by phenotyping individuals with heterozygous predicted damaging variants in WAPL (n = 27), PDS5A (n = 8), and PDS5B (n = 8), by modeling WAPL deficiency in human cells and mice, and by aggregating disease association statistics from consortia studies. We identified a WAPL-related disorder featuring developmental delay, intellectual disability, and risk of other developmental anomalies. Similarities between individuals with damaging WAPL variants and those with large, recurrent 10q22.3q23.2 (10q) deletions encompassing WAPL nominate WAPL as a driver gene within this genomic disorder region. While individuals with PDS5A or PDS5B variants exhibited features of developmental disorders, neither cohort-based statistics nor subject phenotyping associated these genes with specific phenotypes. We used CRISPR to generate truncating variants in WAPL and 10q deletion or duplication in human induced pluripotent stem cells (iPSCs) and induced neurons. Transcriptomics identified significant overlap between WAPL haploinsufficiency and 10q deletion differentially expressed genes. Mice with 50% Wapl expression exhibited mild deficits of growth and learning/memory, whereas those with 25% residual Wapl displayed birth defects and postnatal lethality, revealing a dosage liability threshold below the level of heterozygosity. In summary, we delineated a genetic condition caused by cohesin release factor deficiency, nominated WAPL as a driver gene within a genomic disorder region, and further illuminated dosage sensitivity of human cohesin.

Cohesin is a fundamental genome-organizing complex that orchestrates three-dimensional chromosome folding and gene expression via DNA loop extrusion. Alterations to genes encoding cohesin subunits and cohesin loaders cause Mendelian disorders, including Cornelia de Lange syndrome (CdLS). By contrast, disruption of factors that remove cohesin from DNA, including WAPL and its binding partners PDS5A and PDS5B, have not yet been associated with human disease. Here, we explored the relevance of these cohesin release factors in Mendelian disease by establishing a rare disease cohort of deeply phenotyped individuals with heterozygous, predicted damaging variants in WAPL (n=27), PDS5A (n=8), and PDS5B (n=8), by modeling WAPL deficiency in human cell lines and mice, and by aggregating rare disease association statistics from consortia studies. We identified a WAPL-related disorder characterized by developmental delay, intellectual disability, and risk of other developmental anomalies including clubfoot. Similarities between individuals with damaging WAPL variants and those with large, recurrent 10q22.3q23.2 (10q) deletions (which encompass WAPL) nominate WAPL as a driver gene within this genomic disorder region. While carriers of PDS5A or PDS5B variants exhibited features of developmental disorders, neither cohort-based statistics nor case phenotyping associated these genes with specific phenotypes. We used CRISPR engineering to generate truncating variants in WAPL, as well the 7.8 Mb 10q deletion or duplication in human iPSCs and induced neurons. Transcriptomic analyses identified differentially expressed genes in both models, with highly significant overlap between WAPL haploinsufficiency and 10q deletion signatures. Mice with 50% residual Wapl expression exhibited mild deficits of growth and learning/memory, whereas those with 25% residual Wapl expression displayed birth defects and postnatal lethality, revealing a dosage liability threshold below the level of heterozygosity. In summary, we delineated a novel genetic condition caused by cohesin release factor deficiency, nominated WAPL as a driver gene within a genomic disorder region, and further illuminated dosage sensitivity of human cohesin.

Huntington's disease (HD) is a fatal neurodegenerative disorder caused by inheriting an expanded CAG repeat tract in the huntingtin gene (HTT) that further expands in somatic cells over an individual's lifetime. Genome-wide association studies have provided critical insight into factors that modify the course of disease. These include DNA repair genes that alter the rate of somatic expansion and other genes that do not appear to directly influence this process. One modifier gene is DNA ligase 1 (LIG1), in which a variant specifying a lysine to asparagine substitution (K845N) is associated with a profound (7 to 8 y) delay in the onset of motor signs. Here, we have taken a multifaceted approach to gain insight into the protective nature of this variant in HD. We demonstrate using in vitro ligase assays and enzyme kinetics that K845N enhances discrimination toward mismatched substrates and increases repair fidelity. Consistent with increased ligation fidelity, K845N confers protection against oxidative stress in cell-based assays. Finally, we demonstrate that the mouse LIG1 K843N orthologue suppresses somatic CAG expansion in HD knock-in mice. Overall, our data provide evidence that altered LIG1 function due to the K845N substitution may contribute to HD clinical delay by slowing somatic expansion in the brain and protecting the genome globally against damage. Significantly, our results provide a mechanistic foundation for considering DNA ligase fidelity as a therapeutic target in HD and potentially in other trinucleotide repeat disorders.

Huntington's disease (HD) is driven by somatic expansion of the HTT CAG repeat, with onset modified by genetic factors. One such modifier, 8AM1, maps to chromosome 8 near RRM2B, a gene not directly involved in the machinery that lengthens the repeat. To investigate this locus, we performed capture sequencing and identified variants at both the 5' and 3' ends of RRM2B with expected minor allele frequencies. A polymorphic frameshift variant (rs1037699) in an alternate exon 1 disrupts expression of a previously uncharacterized RRM2B isoform 2, but not isoform 1. Functional analyses in RRM2B knock-out cells and 8AM1 heterozygous LCLs suggest that isoform 2 may function at mitochondria. Several 3' variants, including a 21 bp 3'UTR deletion (rs200678743) and peak tag-SNV (rs79136984), act as cis expression quantitative trait loci. Analysis of HD onset data (n = 12,982) revealed that 5' and 3' variants contribute independently to the 8AM1 modifier effect, with full impact observed only in the absence of the frameshift variant. Knockdown of both isoforms increased neurodegeneration in HD neurons derived from pre-symptomatic patient fibroblasts, supporting an intersection of RRM2B biology and HD pathogenesis. We conclude that the 8AM1 haplotype, present in ~ 14% of Europeans, modifies RRM2B expression in a cell- and context-dependent manner, thereby accelerating HD onset in mutation carriers.

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 onsetdelaying 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 onsetdelaying 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.

Delineation of structural variants (SVs) at sequence resolution in highly repetitive genomic regions has long been intractable. The sequence properties, origins, and functional effects of classes of genomic rearrangements such as ring chromosomes and Robertsonian translocations thus remain unknown. To resolve these complex structures, we leveraged several recent milestones in the field, including (1) the emergence of long-read sequencing, (2) the gapless telomere-to-telomere (T2T) assembly, and (3) a tool (BigClipper) to discover chromosomal rearrangements from long reads. We applied these technologies across 13 cases with ring chromosomes, Robertsonian translocations, and complex SVs that were unresolved by short reads, followed by validation using optical genome mapping (OGM). Our analyses resolved 10 of 13 cases, including a Robertsonian translocation and all ring chromosomes. Multiple breakpoints were localized to genomic regions previously recalcitrant to sequencing such as acrocentric p-arms, ribosomal DNA arrays, and telomeric repeats, and involved complex structures such as a deletion-inversion and interchromosomal dispersed duplications. We further performed methylation profiling from long-read data to discover phased differential methylation in a gene promoter proximal to a ring fusion, suggesting a long-range position effect (LRPE) with heterochromatin spreading. Breakpoint sequences suggested mechanisms of SV formation such as microhomology-mediated and non-homologous end-joining, as well as non-allelic homologous recombination. These methods provide some of the first glimpses into the sequence resolution of Robertsonian translocations and illuminate the structural diversity of ring chromosomes and complex chromosomal rearrangements with implications for genome biology, prediction of LRPEs from integrated multi-omics technologies, and molecular diagnostics in rare disease cases.

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.

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.

Annexin A11 mutations are a rare cause of amyotrophic lateral sclerosis (ALS), wherein replicated protein variants P36R, G38R, D40G and D40Y are located in a small-alpha helix within the long, disordered N-terminus. To elucidate disease mechanisms, we characterised the phenotypes induced by a genetic loss of function (LoF) and by misexpression of G38R and D40G in vivo. Loss of Annexin A11 results in a low-penetrant behavioural phenotype and aberrant axonal morphology in zebrafish homozygous knockout larvae, which is rescued by human WT Annexin A11. Both Annexin A11 knockout/down and ALS variants trigger nuclear dysfunction characterised by Lamin B2 mis-localisation. The Lamin B2 signature also presented in anterior horn, spinal cord neurons from post-mortem ALS+/-FTD patient tissue possessing G38R and D40G protein variants. These findings suggest mutant Annexin A11 acts as a dominant negative, revealing a potential early nucleopathy highlighting nuclear envelope abnormalities preceding behavioural abnormality in animal models.

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.

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.

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.

Short-read genome sequencing (GS) holds the promise of becoming the primary diagnostic approach for the assessment of autism spectrum disorder (ASD) and fetal structural anomalies (FSAs). However, few studies have comprehensively evaluated its performance against current standard-of-care diagnostic tests: karyotype, chromosomal microarray (CMA), and exome sequencing (ES). To assess the clinical utility of GS, we compared its diagnostic yield against these three tests in 1,612 quartet families including an individual with ASD and in 295 prenatal families. Our GS analytic framework identified a diagnostic variant in 7.8% of ASD probands, almost 2-fold more than CMA (4.3%) and 3-fold more than ES (2.7%). However, when we systematically captured copy-number variants (CNVs) from the exome data, the diagnostic yield of ES (7.4%) was brought much closer to, but did not surpass, GS. Similarly, we estimated that GS could achieve an overall diagnostic yield of 46.1% in unselected FSAs, representing a 17.2% increased yield over karyotype, 14.1% over CMA, and 4.1% over ES with CNV calling or 36.1% increase without CNV discovery. Overall, GS provided an added diagnostic yield of 0.4% and 0.8% beyond the combination of all three standard-of-care tests in ASD and FSAs, respectively. This corresponded to nine GS unique diagnostic variants, including sequence variants in exons not captured by ES, structural variants (SVs) inaccessible to existing standard-of-care tests, and SVs where the resolution of GS changed variant classification. Overall, this large-scale evaluation demonstrated that GS significantly outperforms each individual standard-of-care test while also outperforming the combination of all three tests, thus warranting consideration as the first-tier diagnostic approach for the assessment of ASD and FSAs.

Blood neurofilament light chain (NfL) is a robust predictor of phenoconversion in familial frontotemporal lobar degeneration (fFTLD). The comparative value of CSF NfL and other CSF markers of degeneration in fFTLD remains unexplored.