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.

Rare diseases collectively affect 1 in 10 individuals, yet current genetic testing fails to identify a causal variant for most cases. At present, cytogenetic methods and/or sequencing approaches such as exome (ES) or short-read genome sequencing (srGS) represent the state-of-the-art for comprehensive clinical discovery of sequence and structural variants (SVs), including copy number variants, balanced SVs, complex SVs, and tandem repeats (TRs). Recently, long-read genome sequencing (lrGS), coupled with multiomics data, has presented great promise to resolve variation in genomic regions recalcitrant to characterization by srGS such as highly repetitive simple repeat sequences and segmental duplications. However, there are few guidelines to enable clinical interpretation of genetic variation in these highly repetitive genomic regions, and the enthusiasm of the field in adopting lrGS has made it difficult to assess the true added diagnostic yield of this technology due to widely variable and inconsistently applied analytic pipelines and variable degrees of pre-screening by ES or srGS. Here, we investigated the contribution of SVs to rare diseases using srGS as a front-line strategy when paired with highly sensitive SV discovery and evaluate the added diagnostic yield of incorporating lrGS for a subset of cases. Our srGS analysis encompassed 1,462 families (3,450 individuals) recruited through the Broad Institute Center for Mendelian Genetics and the Genomics Research to Elucidate the Genetics of Rare Diseases (GREGoR) programs. Diagnostic SVs were identified in 5.4% of cases (79/1,462), of which 80% were uniquely detectable by srGS compared to standard cytogenetic techniques. For 96 families (including 10 families with a heterozygous variant observed in a known recessive gene of clinical relevance), we performed lrGS with methylation profiling, as well as long-read transcriptomic analyses in a subset of 20 trios. Analyses with lrGS yielded over 25,000 SVs per genome, 63% of which were not captured by srGS, along with an additional ~200 rare SNV/indels per genome not previously captured and 12 differentially methylated regions per genome. Among these, we identified only one diagnostic variant not interpreted by srGS, an apparently mosaic de novo SNV in CASK that was absent in the srGS callset due to allelic imbalance. No new diagnoses were supported by long-read transcriptomics or episignatures. In this well characterized rare disease cohort, the added diagnostic yield was thus 1.04% (1/96 families). Following a systematic literature review of prior lrGS studies, we find that most reported diagnoses were detectable by srGS and that our added diagnostic yield is consistent with those prior studies. These studies emphasize the significant impact of comprehensive SV discovery in rare disease cases and further demonstrate the power for increased discovery of novel genomic variation and episignatures from lrGS. Nonetheless, they also serve to temper expectations of dramatic diagnostic advances in rare disease patients until there is more extensive annotation of the functional and clinical impact of all coding and noncoding variation uniquely accessible to lrGS with extensive reference databases spanning highly repetitive genomic sequencing that could be enabled by this transformative technology.

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.

Tauopathies encompass diverse neurodegenerative diseases unified by aberrant patterns of tau deposition in brain. Although most appear sporadic, some are linked to genetic etiologies that offer unique mechanistic insights. Here we report that X-linked Dystonia-Parkinsonism (XDP), caused by a non-coding retrotransposon-associated repeat insertion in TAF1 , involves a significant imbalance of tau isoforms and the accumulation of hyperphosphorylated, four-repeat tau in the brain. In striatal tissue, both misfolded tau accumulation, predominantly in astrocytes, and MAPT exon 10 inclusion correlated with repeat length within the causal insertion. Transcriptomic profiling across brain regions revealed dysregulation of known tau-related pathways. Levels of phosphorylated tau181, glial fibrillary acidic protein, and neurofilament light chain were elevated in patient plasma and discriminated XDP from controls. These findings implicate defective tau proteostasis as a key pathogenic mechanism and position XDP as a genetic model for uncovering cellular drivers that may disrupt tau in other more common neurodegenerative diseases.

X-linked dystonia parkinsonism (XDP) is a progressive adult-onset neurogenerative disorder caused by the insertion of a SINE-VNTR-Alu (SVA) retrotransposon in TAF1. One element of the SVA is a tandem polymorphic CCCTCT repeat tract whose length inversely correlates with the age of disease onset. Previous observations that the repeat exhibits length-dependent somatic expansion and that XDP onset is modified by variation in DNA repair gene MSH3 indicated that somatic repeat expansion is an important disease driver. Here, we sought to uncover genetic modifiers of CCCTCT instability in XDP individuals and to provide a mechanistic link between somatic instability and disease. We determined quantitative metrics of both repeat expansion and repeat contraction in blood. Using genetic association analyses of exome sequencing data and directed sequencing of a variant MSH3 repeat, we found that MSH3 modifies repeat expansion and contraction in blood as well as age at onset. MSH3 alleles associated with earlier disease onset were associated with more expansion and less contraction. Conversely, alleles associated with later disease onset were associated with less expansion and more contraction. Notably, MSH3 repeat alleles were also similarly associated with expansion and contraction in brain tissues. Our findings provide key evidence that the role of MSH3 in CCCTCT repeat dynamics underlies its impact on clinical disease and indicate that therapeutic strategies to lower or inhibit MSH3 are predicted to both slow CCCTCT expansion and promote CCCTCT contraction, impacting the disease course prior to clinical onset.

Understanding how common and rare genetic variation raises risk for neuropsychiatric disease remains a major challenge. We identify the long arm of chromosome 22 (chr22q) as a region where common polygenic risk for schizophrenia, autism, ADHD, and lower IQ is associated with coordinated downregulation of gene expression in postmortem human brain tissue. The effects are strikingly consistent between neuropsychiatric diagnoses and across brain cell types, and appear to be specific to brain-related traits. We observe that common variant risk for neuropsychiatric diseases has remarkably diffuse expression associations across chr22q, including long-range aggregate associations between genetic variants and genes over 10 Mb away. Polygenic risk for psychiatric disease at chr22q is more strongly associated with lower cognitive ability than elsewhere in the genome, suggesting phenotypic convergence with the 22q11.2 deletion, a rare genetic disorder that causes intellectual disability, schizophrenia, autism, and ADHD. Using human iPSC data, we show that the 22q11.2del induces similarly broad expression downregulation across chr22q in multiple neural cell types and experimental settings. Altogether, our results nominate chr22q as a regulatory hub in neuropsychiatric disease, where common and rare genetic risk factors converge both functionally and phenotypically. ### Competing Interest Statement The authors have declared no competing interest. ### Funding Statement This work was supported by a grant from SFARI (704413, E.B.R). ### Author Declarations I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained. Yes The details of the IRB/oversight body that provided approval or exemption for the research described are given below: All individual level data have been de-identified prior to use in our study. I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals. Yes I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance). Yes I have followed all appropriate research reporting guidelines, such as any relevant EQUATOR Network research reporting checklist(s) and other pertinent material, if applicable. Yes Raw sequencing data for novel datasets will be uploaded to GEO prior to publication; data can be shared with reviewers upon reviewer request.

BACKGROUND:Laryngeal dystonia (LD) is isolated task-specific focal dystonia, predominantly impairing speech production. Clinical observations and population survey studies have reported that up to 58% of patients with LD may have symptom improvement following alcohol intake. OBJECTIVES:To determine the objective characteristics of alcohol responsiveness in LD using a standardized alcohol challenge test and genetic testing. METHODS:A total of 109 patients with isolated focal LD participated in the study. Patients were administered two non-diluted drinks of 40-proof vodka 30 min apart, followed by assessments of voice symptoms, breath alcohol content, and side effects. Patients were considered alcohol-responsive (EtOH+) if their symptoms changed by ≥10% from baseline. Whole-exome sequencing was performed to identify genetic variants associated with the alcohol responsiveness of LD. RESULTS:All patients tolerated the standardized alcohol challenge test without major adverse events. Fifty-two patients (47.7%) had an average of 44.4 ± 25.0% improvement of LD symptoms about 45 min after alcohol intake. Five genetic variants in GABAergic pathway-related genes were enriched in EtOH+ patients, and two of these, rs11644926 (ADCY7) and rs2230741 (ADCY9), were associated with the improvement of dystonic voice symptoms. CONCLUSIONS:Alcohol responsiveness of LD symptoms is a robust feature related to the genes regulating the GABAergic synapses. This finding provides support for the evaluation of novel oral medications with mechanisms of action similar to alcohol for the treatment of patients with alcohol-responsive dystonia.

Recursive splice sites are rare motifs postulated to facilitate splicing across massive introns and shape isoform diversity, especially for long, brain-expressed genes. The necessity of this unique mechanism remains unsubstantiated, as does the role of recursive splicing (RS) in human disease. From analyses of rare copy number variants (CNVs) from almost one million individuals, we previously identified large, heterozygous deletions eliminating an RS site (RS1) in the first intron of CADM2 that conferred substantial risk for attention deficit hyperactivity disorder (ADHD) and other neurobehavioral traits. CADM2 encodes a neuronally expressed cell adhesion molecule that has repeatedly been associated with ADHD and numerous similar traits. To explore the molecular impact of RS ablation in CADM2 , we used CRISPR to model patient deletions and to target a smaller region (~500 base pairs) containing RS1 in both human induced neurons (iNs) and rats. Transcriptome analyses in unedited iNs provided a catalog of CADM2 transcripts, including novel transcripts that retained RS exons. Intriguingly, ablating RS1 altered the gradient of RNA abundance across the first intron of CADM2 , decreased the level of CADM2 expression, and impacted transcript usage. Decreased CADM2 expression was reflected in reduced exon usage downstream of the RS1 site and global alteration to genes involved in neuronal processes including synapse and axon development. Given the scale of our analyses and the widespread association of CADM2 with neurobehavioral traits, we sought to validate these findings using in vivo models and found that rodent models harboring Cadm2 RS1 deletions exhibited significant changes in relevant behaviors and functional brain connectivity. In summary, our analyses demonstrate a functional role for RS as a noncoding regulatory mechanism in a gene associated with a spectrum of neuropsychiatric and behavioral traits.

X-linked Dystonia-Parkinsonism (XDP) is a lethal adult-onset neurodegenerative disorder that exhibits features of dystonia and parkinsonism and is exclusively associated with a causal founder haplotype that is indigenous to the Philippines and affects Filipino males. Using patient-specific fibroblasts, neural stem cells (NSC), and other neuronal models, we discovered that cryptic alternative splicing caused by a novel SINE-VNTR-Alu (SVA) mobile element insertion into intron 32 of TAF1 is a mechanistic hallmark of XDP. We leveraged postmortem brain samples from an XDP-specific brain bank to demonstrate that the molecular hallmarks of XDP observed in neural stem cells (NSCs) mirror abnormalities observed in brain tissues from affected patients. Based on these findings that patient-specific NSCs reproduce mechanistic signatures found in the brain, we sought to develop a bespoke precision therapeutic for XDP and evaluate its relative efficacy in ameliorating transcriptomic signatures in neuronal models. We first used CRISPR-based excision of the SVA and demonstrated ablation of all aberrant splicing and dysregulation of TAF1 expression in NSCs across 30 independent clones. CRISPR-based correction of the XDP haplotype also restored the expression of 424 of 1,490 (30%) differentially expressed genes (DEGs) that were altered in XDP patient lines and greatly exceeded what would be expected by chance (p-value = 9.89e-87). While in vivo delivery of a gold standard CRISPR therapy is currently not feasible for XDP, we evaluated a tractable approach for Filipino patients by exploring the potential to modulate alternative splicing in XDP patients using antisense oligonucleotides (ASOs). To accomplish this, we developed a large-scale and well controlled functional genomics platform that screened eighty ASOs targeting intron 32 of XDP patients, followed by prioritization of lead ASOs based on attenuation of the alternative splicing signature. In transcriptome analyses across 1,550 libraries, we found that 8 of the 12 lead ASOs ameliorated the targeted XDP aberrant splicing. Moreover, we found that the two lead ASOs exhibited 38% and 43% rescue of XDP-specific DEGs that were also rescued by CRISPR excision of the SVA (enrichment p-values = 2.06e-13 and 2.27e-05, respectively). These rescues represented restoration of key molecular functions previously implicated in XDP, such as synaptic function, DNA-binding transcription factor activity, and gliogenesis. This study highlights a path to a potential targeted therapeutic for XDP and the capacity to exploit functional genomic signatures in patient-derived neural models to develop a scalable precision therapeutic platform for rare genetic disorders.

Familial dysautonomia (FD) is a rare recessive neurodevelopmental disease caused by a splice mutation in the Elongator acetyltransferase complex subunit 1 ( ELP1 ) gene. This mutation results in a tissue-specific reduction of ELP1 protein, with the lowest levels in the central and peripheral nervous systems (CNS and PNS, respectively). FD patients exhibit complex neurological phenotypes due to the loss of sensory and autonomic neurons. Disease symptoms include decreased pain and temperature perception, impaired or absent myotatic reflexes, proprioceptive ataxia, and progressive retinal degeneration. While the involvement of the PNS in FD pathogenesis has been clearly recognized, the underlying mechanisms responsible for the preferential neuronal loss remain unknown. In this study, we aimed to elucidate the molecular mechanisms underlying FD by conducting a comprehensive transcriptome analysis of neuronal tissues from the phenotypic mouse model TgFD9 ; Elp1 Δ 20/flox . This mouse recapitulates the same tissue-specific ELP1 mis-splicing observed in patients while modeling many of the disease manifestations. Comparison of FD and control transcriptomes from dorsal root ganglion (DRG), trigeminal ganglion (TG), medulla (MED), cortex, and spinal cord (SC) showed significantly more differentially expressed genes (DEGs) in the PNS than the CNS. We then identified genes that were tightly co-expressed and functionally dependent on the level of full-length ELP1 transcript. These genes, defined as ELP1 dose-responsive genes, were combined with the DEGs to generate tissue-specific dysregulated FD signature genes and networks. Within the PNS networks, we observed direct connections between Elp1 and genes involved in tRNA synthesis and genes related to amine metabolism and synaptic signaling. Importantly, transcriptomic dysregulation in PNS tissues exhibited enrichment for neuronal subtype markers associated with peptidergic nociceptors and myelinated sensory neurons, which are known to be affected in FD. In summary, this study has identified critical tissue-specific gene networks underlying the etiology of FD and provides new insights into the molecular basis of the disease.

Glycine is an obligatory co-agonist at excitatory NMDA receptors in the brain, especially in the dentate gyrus, which has been postulated to be crucial for the development of psychotic associations and memories with psychotic content. Drugs modulating glycine levels are in clinical development for improving cognition in schizophrenia. However, the functional relevance of the regulation of glycine metabolism by endogenous enzymes is unclear. Using a chromosome-engineered allelic series in mice, we report that a triplication of the gene encoding the glycine-catabolizing enzyme glycine decarboxylase (GLDC) - as found on a small supernumerary marker chromosome in patients with psychosis - reduces extracellular glycine levels as determined by optical fluorescence resonance energy transfer (FRET) in dentate gyrus (DG) and suppresses long-term potentiation (LTP) in mPP-DG synapses but not in CA3-CA1 synapses, reduces the activity of biochemical pathways implicated in schizophrenia and mitochondrial bioenergetics, and displays deficits in schizophrenia-like behaviors which are in part known to be dependent on the activity of the dentate gyrus, e.g., prepulse inhibition, startle habituation, latent inhibition, working memory, sociability and social preference. Our results demonstrate that Gldc negatively regulates long-term synaptic plasticity in the dentate gyrus in mice, suggesting that an increase in GLDC copy number possibly contributes to the development of psychosis in humans.

There is a paucity of genetic characterization in people with Parkinson’s disease (PD) of Latino and Afro-Caribbean descent. Screening LRRK2 and GBA variants in 32 New Yorkers of Puerto Rican ethnicity with PD and in 119 non-Hispanic-non-Jewish European PD cases revealed that Puerto Rican participants were more likely to harbor the LRRK2-p.G2019S variant (15.6% vs. 4.2%, respectively). Additionally, whole exome sequencing of twelve Puerto Rican and Dominican PD participants was performed as an exploratory study.

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.

The biological significance of a small supernumerary marker chromosome that results in dosage alterations to chromosome 9p24.1, including triplication of the GLDC gene encoding glycine decarboxylase, in two patients with psychosis is unclear. In an allelic series of copy number variant mouse models, we identify that triplication of Gldc reduces extracellular glycine levels as determined by optical fluorescence resonance energy transfer (FRET) in dentate gyrus (DG) but not in CA1, suppresses long-term potentiation (LTP) in mPP-DG synapses but not in CA3-CA1 synapses, reduces the activity of biochemical pathways implicated in schizophrenia and mitochondrial bioenergetics, and displays deficits in prepulse inhibition, startle habituation, latent inhibition, working memory, sociability and social preference. Our results thus provide a link between a genomic copy number variation, biochemical, cellular and behavioral phenotypes, and further demonstrate that GLDC negatively regulates long-term synaptic plasticity at specific hippocampal synapses, possibly contributing to the development of neuropsychiatric disorders.

The canonical paradigm for converting genetic association to mechanism involves iteratively mapping individual associations to the proximal genes through which they act. In contrast, in the present study we demonstrate the feasibility of extracting biological insights from a very large region of the genome and leverage this strategy to study the genetic influences on autism. Using a new statistical approach, we identified the 33-Mb p-arm of chromosome 16 (16p) as harboring the greatest excess of autism’s common polygenic influences. The region also includes the mechanistically cryptic and autism-associated 16p11.2 copy number variant. Analysis of RNA-sequencing data revealed that both the common polygenic influences within 16p and the 16p11.2 deletion were associated with decreased average gene expression across 16p. The transcriptional effects of the rare deletion and diffuse common variation were correlated at the level of individual genes and analysis of Hi-C data revealed patterns of chromatin contact that may explain this transcriptional convergence. These results reflect a new approach for extracting biological insight from genetic association data and suggest convergence of common and rare genetic influences on autism at 16p.