Genomic rearrangements are a common cause of human congenital abnormalities. However, their origin and consequences are poorly understood. We performed molecular analysis of two patients with congenital disease who carried de novo genomic rearrangements. We found that the rearrangements in both patients hit genes that are recurrently rearranged in cancer (ETV1, FOXP1, and microRNA cluster C19MC) and drive formation of fusion genes similar to those described in cancer. Subsequent analysis of a large set of 552 de novo germline genomic rearrangements underlying congenital disorders revealed enrichment for genes rearranged in cancer and overlap with somatic cancer breakpoints. Breakpoints of common (inherited) germline structural variations also overlap with cancer breakpoints but are depleted for cancer genes. We propose that the same genomic positions are prone to genomic rearrangements in germline and soma but that timing and context of breakage determines whether developmental defects or cancer are promoted.

With recent rapid advances in genomic technologies, precise delineation of structural chromosome rearrangements at the nucleotide level is becoming increasingly feasible. In this era of "next-generation cytogenetics" (i.e., an integration of traditional cytogenetic techniques and next-generation sequencing), a consensus nomenclature is essential for accurate communication and data sharing. Currently, nomenclature for describing the sequencing data of these aberrations is lacking. Herein, we present a system called Next-Gen Cytogenetic Nomenclature, which is concordant with the International System for Human Cytogenetic Nomenclature (2013). This system starts with the alignment of rearrangement sequences by BLAT or BLAST (alignment tools) and arrives at a concise and detailed description of chromosomal changes. To facilitate usage and implementation of this nomenclature, we are developing a program designated BLA(S)T Output Sequence Tool of Nomenclature (BOSToN), a demonstrative version of which is accessible online. A standardized characterization of structural chromosomal rearrangements is essential both for research analyses and for application in the clinical setting.

Persons with neurodevelopmental disorders or autism spectrum disorder (ASD) often harbor chromosomal microdeletions, yet the individual genetic contributors within these regions have not been systematically evaluated. We established a consortium of clinical diagnostic and research laboratories to accumulate a large cohort with genetic alterations of chromosomal region 2q23.1 and acquired 65 subjects with microdeletion or translocation. We sequenced translocation breakpoints; aligned microdeletions to determine the critical region; assessed effects on mRNA expression; and examined medical records, photos, and clinical evaluations. We identified a single gene, methyl-CpG-binding domain 5 (MBD5), as the only locus that defined the critical region. Partial or complete deletion of MBD5 was associated with haploinsufficiency of mRNA expression, intellectual disability, epilepsy, and autistic features. Fourteen alterations, including partial deletions of noncoding regions not typically captured or considered pathogenic by current diagnostic screening, disrupted MBD5 alone. Expression profiles and clinical characteristics were largely indistinguishable between MBD5-specific alteration and deletion of the entire 2q23.1 interval. No copy-number alterations of MBD5 were observed in 7878 controls, suggesting MBD5 alterations are highly penetrant. We surveyed MBD5 coding variations among 747 ASD subjects compared to 2043 non-ASD subjects analyzed by whole-exome sequencing and detected an association with a highly conserved methyl-CpG-binding domain missense variant, p.79Gly>Glu (c.236G>A) (p = 0.012). These results suggest that genetic alterations of MBD5 cause features of 2q23.1 microdeletion syndrome and that this epigenetic regulator significantly contributes to ASD risk, warranting further consideration in research and clinical diagnostic screening and highlighting the importance of chromatin remodeling in the etiology of these complex disorders.

Schwannomatosis has been linked to germline mutations in the SMARCB1 and LZTR1 genes, and is frequently associated with pain.In a cohort study, we assessed the mutation status of 37 patients with clinically diagnosed schwannomatosis and compared to clinical data, whole body MRI (WBMRI), visual analog pain scale, and Short Form 36 (SF-36) bodily pain subscale.We identified a germline mutation in LZTR1 in 5 patients (13.5%) and SMARCB1 in 15 patients (40.5%), but found no germline mutation in 17 patients (45.9%). Peripheral schwannomas were detected in 3 LZTR1-mutant (60%) and 10 SMARCB1-mutant subjects (66.7%). Among those with peripheral tumors, the median tumor number was 4 in the LZTR1 group (median total body tumor volume 30 cc) and 10 in the SMARCB1 group (median volume 85cc), (P=.2915 for tumor number and P = .2289 for volume). mutation was associated with an increased prevalence of spinal schwannomas (100% vs 41%, P = .0197). The median pain score was 3.9/10 in the LZTR1 group and 0.5/10 in the SMARCB1 group (P = .0414), and SF-36 pain-associated quality of life was significantly worse in the LZTR1 group (P = .0106). Pain scores correlated with total body tumor volume (rho = 0.32471, P = .0499), but not with number of tumors (rho = 0.23065, P = .1696).We found no significant difference in quantitative tumor burden between mutational groups, but spinal schwannomas were more common in LZTR1-mutant patients. Pain was significantly higher in LZTR1-mutant than in SMARCB1-mutant patients, though spinal tumor location did not significantly correlate with pain. This suggests a possible genetic association with schwannomatosis-associated pain.

The 6%-9% risk of an untoward outcome previously established by Warburton for prenatally detected de novo balanced chromosomal rearrangements (BCRs) does not account for long-term morbidity. We performed long-term follow-up (mean 17 years) of a registry-based nationwide cohort of 41 individuals carrying a prenatally detected de novo BCR with normal first trimester screening/ultrasound scan. We observed a significantly higher frequency of neurodevelopmental and/or neuropsychiatric disorders than in a matched control group (19.5% versus 8.3%, p = 0.04), which was increased to 26.8% upon clinical follow-up. Chromosomal microarray of 32 carriers revealed no pathogenic imbalances, illustrating a low prognostic value when fetal ultrasound scan is normal. In contrast, mate-pair sequencing revealed disrupted genes (ARID1B, NPAS3, CELF4), regulatory domains of known developmental genes (ZEB2, HOXC), and complex BCRs associated with adverse outcomes. Seven unmappable autosomal-autosomal BCRs with breakpoints involving pericentromeric/heterochromatic regions may represent a low-risk group. We performed independent phenotype-aware and blinded interpretation, which accurately predicted benign outcomes (specificity = 100%) but demonstrated relatively low sensitivity for prediction of the clinical outcome in affected carriers (sensitivity = 45%-55%). This sensitivity emphasizes the challenges associated with prenatal risk prediction for long-term morbidity in the absence of phenotypic data given the still immature annotation of the morbidity genome and poorly understood long-range regulatory mechanisms. In conclusion, we upwardly revise the previous estimates of Warburton to a morbidity risk of 27% and recommend sequencing of the chromosomal breakpoints as the first-tier diagnostic test in pregnancies with a de novo BCR.

Despite improvements in genomics technology, the detection of structural variants (SVs) from short-read sequencing still poses challenges, particularly for complex variation. Here we analyse the genomes of two patients with congenital abnormalities using the MinION nanopore sequencer and a novel computational pipeline-NanoSV. We demonstrate that nanopore long reads are superior to short reads with regard to detection of de novo chromothripsis rearrangements. The long reads also enable efficient phasing of genetic variations, which we leveraged to determine the parental origin of all de novo chromothripsis breakpoints and to resolve the structure of these complex rearrangements. Additionally, genome-wide surveillance of inherited SVs reveals novel variants, missed in short-read data sets, a large proportion of which are retrotransposon insertions. We provide a first exploration of patient genome sequencing with a nanopore sequencer and demonstrate the value of long-read sequencing in mapping and phasing of SVs for both clinical and research applications.

Schizophrenia (SCZ) is a debilitating psychiatric disorder with a large genetic contribution; however, its neurodevelopmental substrates remain largely unknown. Modeling pathogenic processes in SCZ using human induced pluripotent stem cell-derived neurons (iNs) has emerged as a promising strategy. Copy number variants confer high genetic risk for SCZ, with duplication of the 16p11.2 locus increasing the risk 14.5-fold.

Mitochondria dysfunction occurs in the aging brain as well as in several neurodegenerative disorders and predisposes neuronal cells to enhanced sensitivity to neurotoxins. 3-nitropropionic acid (3-NP) is a naturally occurring plant and fungal neurotoxin that causes neurodegeneration predominantly in the striatum by irreversibly inhibiting the tricarboxylic acid respiratory chain enzyme, succinate dehydrogenase (SDH), the main constituent of the mitochondria respiratory chain complex II. Significantly, although 3-NP-induced inhibition of SDH occurs in all brain regions, neurodegeneration occurs primarily and almost exclusively in the striatum for reasons still not understood. In rodents, 3-NP-induced striatal neurodegeneration depends on the strain background suggesting that genetic differences among genotypes modulate toxicant variability and mechanisms that underlie 3-NP-induced neuronal cell death. Using the large BXD family of recombinant inbred (RI) strains we demonstrate that variants in Ccnd1 - the gene encoding cyclin D1 - of the DBA/2 J parent underlie the resistance to 3-NP-induced striatal neurodegeneration. In contrast, the Ccnd1 variant inherited from the widely used C57BL/6 J parental strain confers sensitivity. Given that cellular stress triggers induction of cyclin D1 expression followed by cell-cycle re-entry and consequent neuronal cell death, we sought to determine if the C57BL/6 J and DBA/2 J Ccnd1 variants are differentially modulated in response to 3-NP. We confirm that 3-NP induces cyclin D1 expression in striatal neuronal cells of C57BL/6 J, but this response is blunted in the DBA/2 J. We further show that striatal-specific alternative processing of a highly conserved 3'UTR negative regulatory region of Ccnd1 co-segregates with the C57BL/6 J parental Ccnd1 allele in BXD strains and that its differential processing accounts for sensitivity or resistance to 3-NP. Our results indicate that naturally occurring Ccnd1 variants may play a role in the variability observed in neurodegenerative disorders involving mitochondria complex II dysfunction and point to cyclin D1 as a possible therapeutic target.

At the start of the COVID-19 pandemic, the BNT162b2 (BioNTech-Pfizer) and mRNA-1273 (Moderna) mRNA vaccines were expediently designed and mass produced. Both vaccines produce the full-length SARS-CoV-2 spike protein for gain of immunity and have greatly reduced mortality and morbidity from SARS-CoV-2 infection. The distribution and duration of SARS-CoV-2 mRNA vaccine persistence in human tissues is unclear. Here, we developed specific RT-qPCR-based assays to detect each mRNA vaccine and screened lymph nodes, liver, spleen, and myocardium from recently vaccinated deceased patients. Vaccine was detected in the axillary lymph nodes in the majority of patients dying within 30 days of vaccination, but not in patients dying more than 30 days from vaccination. Vaccine was not detected in the mediastinal lymph nodes, spleen, or liver. Vaccine was detected in the myocardium in a subset of patients vaccinated within 30 days of death. Cardiac ventricles in which vaccine was detected had healing myocardial injury at the time of vaccination and had more myocardial macrophages than the cardiac ventricles in which vaccine was not detected. These results suggest that SARS-CoV-2 mRNA vaccines routinely persist up to 30 days from vaccination and can be detected in the heart.

Familial Dysautonomia (FD) is an autosomal recessive disorder that affects 1/3,600 live births in the Ashkenazi Jewish population, and leads to death before the age of 40. The disease is characterized by abnormal development and progressive degeneration of the sensory and autonomic nervous system. A single base pair substitution in intron 20 of the Ikbkap gene accounts for 98% of FD cases, and results in the expression of low levels of the full-length mRNA with simultaneous expression of an aberrantly spliced mRNA in which exon 20 is missing. To date, there is no animal model for the disease, and the essential cellular functions of IKAP--the protein encoded by Ikbkap--remain unknown. To better understand the normal function of IKAP and in an effort to generate a mouse model for FD, we have targeted the mouse Ikbkap gene by homologous recombination. We created two distinct alleles that result in either loss of Ikbkap expression, or expression of an mRNA lacking only exon 20. Homozygosity for either mutation leads to developmental delay, cardiovascular and brain malformations, accompanied with early embryonic lethality. Our analyses indicate that IKAP is essential for expression of specific genes involved in cardiac morphogenesis, and that cardiac failure is the likely cause of abnormal vascular development and embryonic lethality. Our results also indicate that deletion of exon 20 abolishes gene function. This implies that the truncated IKAP protein expressed in FD patients does not retain any significant biological function.

The splice site mutation in the IKBKAP gene coding for IKAP protein leads to the tissue-specific skipping of exon 20, with concomitant reduction in IKAP protein production. This causes the neurodevelopmental, autosomal-recessive genetic disorder - Familial Dysautonomia (FD). The molecular hallmark of FD is the severe reduction of IKAP protein in the nervous system that is believed to be the main reason for the devastating symptoms of this disease. Our recent studies showed that in the brain of two FD patients, genes linked to oligodendrocyte differentiation and/or myelin formation are significantly downregulated, implicating IKAP in the process of myelination. However, due to the scarcity of FD patient tissues, these results awaited further validation in other models. Recently, two FD mouse models that faithfully recapitulate FD were generated, with two types of mutations resulting in severely low levels of IKAP expression. Here we demonstrate that IKAP deficiency in these FD mouse models affects a similar set of genes as in FD patients' brains. In addition, we identified two new IKAP target genes involved in oligodendrocyte cells differentiation and myelination, further underscoring the essential role of IKAP in this process. We also provide proof that IKAP expression is needed cell-autonomously for the regulation of expression of genes involved in myelin formation since knockdown of IKAP in the Oli-neu oligodendrocyte precursor cell line results in similar deficiencies. Further analyses of these two experimental models will compensate for the lack of human postmortem tissues and will advance our understanding of the role of IKAP in myelination and the disease pathology.

Structural variation (SV) is a significant component of the genetic etiology of both neurodevelopmental and psychiatric disorders; however, routine guidelines for clinical genetic screening have been established only in the former category. Genome-wide chromosomal microarray (CMA) can detect genomic imbalances such as copy-number variants (CNVs), but balanced chromosomal abnormalities (BCAs) still require karyotyping for clinical detection. Moreover, submicroscopic BCAs and subarray threshold CNVs are intractable, or cryptic, to both CMA and karyotyping. Here, we performed whole-genome sequencing using large-insert jumping libraries to delineate both cytogenetically visible and cryptic SVs in a single test among 30 clinically referred youth representing a range of severe neuropsychiatric conditions. We detected 96 SVs per person on average that passed filtering criteria above our highest-confidence resolution (6,305 bp) and an additional 111 SVs per genome below this resolution. These SVs rearranged 3.8 Mb of genomic sequence and resulted in 42 putative loss-of-function (LoF) or gain-of-function mutations per person. We estimate that 80% of the LoF variants were cryptic to clinical CMA. We found myriad complex and cryptic rearrangements, including a "paired" duplication (360 kb, 169 kb) that flanks a 5.25 Mb inversion that appears in 7 additional cases from clinical CNV data among 47,562 individuals. Following convergent genomic profiling of these independent clinical CNV data, we interpreted three SVs to be of potential clinical significance. These data indicate that sequence-based delineation of the full SV mutational spectrum warrants exploration in youth referred for neuropsychiatric evaluation and clinical diagnostic SV screening more broadly.

Rare individuals with inactivating mutations in the Huntington's disease gene (HTT) exhibit variable abnormalities that imply essential HTT roles during organ development. Here we report phenotypes produced when increasingly severe hypomorphic mutations in the murine HTT orthologue Htt, (HdhneoQ20, HdhneoQ50, HdhneoQ111), were placed over a null allele (Hdhex4/5). The most severe hypomorphic allele failed to rescue null lethality at gastrulation, while the intermediate, though still severe, alleles yielded recessive perinatal lethality and a variety of fetal abnormalities affecting body size, skin, skeletal and ear formation, and transient defects in hematopoiesis. Comparative molecular analysis of wild-type and Htt-null retinoic acid-differentiated cells revealed gene network dysregulation associated with organ development that nominate polycomb repressive complexes and miRNAs as molecular mediators. Together these findings demonstrate that Htt is required both pre- and post-gastrulation to support normal development.

Loss-of-function mutations in SGCE, which encodes ε-sarcoglycan (ε-SG), cause myoclonus-dystonia syndrome (OMIM159900, DYT11). A "major" ε-SG protein derived from CCDS5637.1 (NM_003919.2) and a "brain-specific" protein, that includes sequence derived from alternative exon 11b (CCDS47642.1, NM_001099400.1), are reportedly localized in post- and pre-synaptic membrane fractions, respectively. Moreover, deficiency of the "brain-specific" isoform and other isoforms derived from exon 11b may be central to the pathogenesis of DYT11. However, no animal model supports this hypothesis. Gene-trapped ES cells (CMHD-GT_148G1-3, intron 9 of NM_011360) were used to generate a novel Sgce mouse model (C57BL/6J background) with markedly reduced expression of isoforms derived from exons 3' to exon 9 of NM_011360. Among those brain regions analyzed in adult (2month-old) wild-type (WT) mice, cerebellum showed the highest relative expression of isoforms incorporating exon 11b. Homozygotes (SgceGt(148G1)Cmhd/Gt(148G1)Cmhd or SgceGt/Gt) and paternal heterozygotes (Sgcem+/pGt, m-maternal, p-paternal) showed 60 to 70% reductions in expression of total Sgce. Although expression of the major (NM_011360) and brain-specific (NM_001130189) isoforms was markedly reduced, expression of short isoforms was preserved and relatively small amounts of chimeric ε-SG/β-galactosidase fusion protein was produced by the Sgce gene-trap locus. Immunoaffinity purification followed by mass spectrometry assessments of Sgcem+/pGt mouse brain using pan- or brain-specific ε-SG antibodies revealed significant reductions of ε-SG and other interacting sarcoglycans. Genome-wide gene-expression data using RNA derived from adult Sgcem+/pGt mouse cerebellum showed that the top up-regulated genes were involved in cell cycle, cellular development, cell death and survival, while the top down-regulated genes were associated with protein synthesis, cellular development, and cell death and survival. In comparison to WT littermates, Sgcem+/pGt mice exhibited "tiptoe" gait and stimulus-induced appendicular posturing between Postnatal Days 14 to 16. Abnormalities noted in older Sgcem+/pGt mice included reduced body weight, altered gait dynamics, and reduced open-field activity. Overt spontaneous or stimulus-sensitive myoclonus was not apparent on the C57BL/6J background or mixed C57BL/6J-BALB/c and C57BL/6J-129S2 backgrounds. Our data confirm that mouse Sgce is a maternally imprinted gene and suggests that short Sgce isoforms may compensate, in part, for deficiency of major and brain-specific Sgce isoforms.

Huntington's disease (HD) is caused by an autosomal dominant polyglutamine expansion mutation of Huntingtin (HTT). HD patients suffer from progressive motor, cognitive, and psychiatric impairments, along with significant degeneration of the striatal projection neurons (SPNs) of the striatum. HD is widely accepted to be caused by a toxic gain-of-function of mutant HTT. However, whether loss of HTT function, because of dominant-negative effects of the mutant protein, plays a role in HD and whether HTT is required for SPN health and function are not known. Here, we delete Htt from specific subpopulations of SPNs using the Cre-Lox system and find that SPNs require HTT for motor regulation, synaptic development, cell health, and survival during aging. Our results suggest that loss of HTT function in SPNs could play a critical role in HD pathogenesis.

We present the largest exome sequencing study of autism spectrum disorder (ASD) to date (n = 35,584 total samples, 11,986 with ASD). Using an enhanced analytical framework to integrate de novo and case-control rare variation, we identify 102 risk genes at a false discovery rate of 0.1 or less. Of these genes, 49 show higher frequencies of disruptive de novo variants in individuals ascertained to have severe neurodevelopmental delay, whereas 53 show higher frequencies in individuals ascertained to have ASD; comparing ASD cases with mutations in these groups reveals phenotypic differences. Expressed early in brain development, most risk genes have roles in regulation of gene expression or neuronal communication (i.e., mutations effect neurodevelopmental and neurophysiological changes), and 13 fall within loci recurrently hit by copy number variants. In cells from the human cortex, expression of risk genes is enriched in excitatory and inhibitory neuronal lineages, consistent with multiple paths to an excitatory-inhibitory imbalance underlying ASD.

Gene expression levels vary across developmental stage, cell type, and region in the brain. Genomic variants also contribute to the variation in expression, and some neuropsychiatric disorder loci may exert their effects through this mechanism. To investigate these relationships, we present BrainVar, a unique resource of paired whole-genome and bulk tissue RNA sequencing from the dorsolateral prefrontal cortex of 176 individuals across prenatal and postnatal development. Here we identify common variants that alter gene expression (expression quantitative trait loci [eQTLs]) constantly across development or predominantly during prenatal or postnatal stages. Both "constant" and "temporal-predominant" eQTLs are enriched for loci associated with neuropsychiatric traits and disorders and colocalize with specific variants. Expression levels of more than 12,000 genes rise or fall in a concerted late-fetal transition, with the transitional genes enriched for cell-type-specific genes and neuropsychiatric risk loci, underscoring the importance of cataloging developmental trajectories in understanding cortical physiology and pathology.

Gonadotropin-releasing hormone (GnRH) neurons in the hypothalamus play a key role in the regulation of reproductive function. In this study, we sought an efficient method for generating GnRH neurons from human embryonic and induced pluripotent stem cells (hESC and hiPSC, respectively). First, we found that exposure of primitive neuroepithelial cells, rather than neuroprogenitor cells, to fibroblast growth factor 8 (FGF8), was more effective in generating GnRH neurons. Second, addition of kisspeptin to FGF8 further increased the efficiency rates of GnRH neurogeneration. Third, we generated a fluorescent marker mCherry labeled human embryonic GnRH cell line (mCh-hESC) using a CRISPR-Cas9 targeting approach. Fourth, we examined physiological characteristics of GnRH (mCh-hESC) neurons: similar to GnRH neurons in vivo, they released the GnRH peptide in a pulsatile manner at ~60 min intervals; GnRH release increased in response to high potassium, kisspeptin, estradiol, and neurokinin B challenges; and injection of depolarizing current induced action potentials. Finally, we characterized developmental changes in transcriptomes of GnRH neurons using hESC, hiPSC, and mCh-hESC. The developmental pattern of transcriptomes was remarkably similar among the 3 cell lines. Collectively, human stem cell-derived GnRH neurons will be an important tool for establishing disease models to understand diseases, such as idiopathic hypothalamic hypogonadism, and testing contraceptive drugs.

Pre-mRNA splicing is a key controller of human gene expression. Disturbances in splicing due to mutation lead to dysregulated protein expression and contribute to a substantial fraction of human disease. Several classes of splicing modulator compounds (SMCs) have been recently identified and establish that pre-mRNA splicing represents a target for therapy. We describe herein the identification of BPN-15477, a SMC that restores correct splicing of ELP1 exon 20. Using transcriptome sequencing from treated fibroblast cells and a machine learning approach, we identify BPN-15477 responsive sequence signatures. We then leverage this model to discover 155 human disease genes harboring ClinVar mutations predicted to alter pre-mRNA splicing as targets for BPN-15477. Splicing assays confirm successful correction of splicing defects caused by mutations in CFTR, LIPA, MLH1 and MAPT. Subsequent validations in two disease-relevant cellular models demonstrate that BPN-15477 increases functional protein, confirming the clinical potential of our predictions.

Genetic variation within the transcription factor TCF4 locus can cause the intellectual disability and developmental disorder Pitt-Hopkins syndrome (PTHS), whereas single-nucleotide polymorphisms within noncoding regions are associated with schizophrenia. These genetic findings position TCF4 as a link between transcription and cognition; however, the neurobiology of TCF4 remains poorly understood. Here, we quantitated multiple distinct TCF4 transcript levels in human induced pluripotent stem cell-derived neural progenitors and differentiated neurons, and PTHS patient fibroblasts. We identify two classes of pharmacological treatments that regulate TCF4 expression: WNT pathway activation and inhibition of class I histone deacetylases. In PTHS fibroblasts, both of these perturbations upregulate a subset of TCF4 transcripts. Finally, using chromatin immunoprecipitation sequencing in conjunction with genome-wide transcriptome analysis, we identified TCF4 target genes that may mediate the effect of TCF4 loss on neuroplasticity. Our studies identify new pharmacological assays, tools, and targets for the development of therapeutics for cognitive disorders.