Search

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.

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.

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.

Autism spectrum disorder (ASD) is a heritable neurodevelopmental disorder characterized by deficits in social interactions and communication. Protein-altering variants in many genes have been shown to contribute to ASD; however, understanding the convergence across many genes remains a challenge. We demonstrate that coexpression patterns from 993 human postmortem brains are significantly correlated with the transcriptional consequences of CRISPR perturbations in human neurons. Across 71 ASD risk genes, there was significant tissue-specific convergence implicating synaptic pathways. Tissue-specific convergence was further demonstrated across schizophrenia and atrial fibrillation risk genes. The degree of ASD convergence was significantly correlated with ASD association from rare variation and differential expression in ASD brains. Positively convergent genes showed intolerance to functional mutations and had shorter coding lengths than known risk genes even after removing association with ASD. These results indicate that convergent coexpression can identify potentially novel genes that are unlikely to be discovered by sequencing studies.

SUMMARYCopy number variants (CNVs) are major contributors to genetic diversity and disease. To date, exome sequencing (ES) has been generated for millions of individuals in international biobanks, human disease studies, and clinical diagnostic screening. While standardized methods exist for detecting short variants (single nucleotide and insertion/deletion variants) using tools such as the Genome Analysis ToolKit (GATK), technical challenges have confounded similarly uniform large-scale CNV analyses from ES data. Given the profound impact of rare and de novo coding CNVs on genome organization and human disease, the lack of widely-adopted and robustly benchmarked rare CNV discovery tools has presented a barrier to routine exome-wide assessment of this critical class of variation. Here, we introduce GATK-gCNV, a flexible algorithm to discover rare CNVs from genome sequencing read-depth information, which we distribute as an open-source tool packaged in GATK. GATK-gCNV uses a probabilistic model and inference framework that accounts for technical biases while simultaneously predicting CNVs, which enables self-consistency between technical read-depth normalization and variant calling. We benchmarked GATK-gCNV in 7,962 exomes from individuals in quartet families with matched genome sequencing and microarray data. These analyses demonstrated 97% recall of rare (≤1% site frequency) coding CNVs detected by microarrays and 95% recall of rare coding CNVs discovered by genome sequencing at a resolution of more than two exons. We applied GATK-gCNV to generate a reference catalog of rare coding CNVs in 197,306 individuals with ES from the UK Biobank. We observed strong correlations between CNV rates per gene and measures of mutational constraint, as well as rare CNV associations with multiple traits. In summary, GATK-gCNV is a tunable approach for sensitive and specific CNV discovery in ES, which can easily be applied across trait association and clinical screening.

Autism spectrum disorder (ASD) is diagnosed three to four times more frequently in males than in females. Genetic studies of rare variants support a female protective effect (FPE) against ASD. However, sex differences in common inherited genetic risk for ASD are less studied, particularly within families. Leveraging the Danish iPSYCH resource, we found siblings of female ASD cases (n = 1,707) had higher rates of ASD than siblings of male ASD cases (n = 6,270; p < 1.0 × 10-10). In the Simons Simplex and SPARK collections, mothers of ASD cases (n = 7,436) carried more polygenic risk for ASD than fathers of ASD cases (n = 5,926; 0.08 polygenic risk score [PRS] SD; p = 7.0 × 10-7). Further, male unaffected siblings under-inherited polygenic risk (n = 1,519; p = 0.03). Using both epidemiologic and genetic approaches, our findings strongly support an FPE against ASD's common inherited influences.

Virtually all genome sequencing efforts in national biobanks, complex and Mendelian disease programs, and medical genetic initiatives are reliant upon short-read whole-genome sequencing (srWGS), which presents challenges for the detection of structural variants (SVs) relative to emerging long-read WGS (lrWGS) technologies. Given this ubiquity of srWGS in large-scale genomics initiatives, we sought to establish expectations for routine SV detection from this data type by comparison with lrWGS assembly, as well as to quantify the genomic properties and added value of SVs uniquely accessible to each technology. Analyses from the Human Genome Structural Variation Consortium (HGSVC) of three families captured ~11,000 SVs per genome from srWGS and ~25,000 SVs per genome from lrWGS assembly. Detection power and precision for SV discovery varied dramatically by genomic context and variant class: 9.7% of the current GRCh38 reference is defined by segmental duplication (SD) and simple repeat (SR), yet 91.4% of deletions that were specifically discovered by lrWGS localized to these regions. Across the remaining 90.3% of reference sequence, we observed extremely high (93.8%) concordance between technologies for deletions in these datasets. In contrast, lrWGS was superior for detection of insertions across all genomic contexts. Given that non-SD/SR sequences encompass 95.9% of currently annotated disease-associated exons, improved sensitivity from lrWGS to discover novel pathogenic deletions in these currently interpretable genomic regions is likely to be incremental. However, these analyses highlight the considerable added value of assembly-based lrWGS to create new catalogs of insertions and transposable elements, as well as disease-associated repeat expansions in genomic sequences that were previously recalcitrant to routine assessment.

Summary Each human genome includes de novo mutations that arose during gametogenesis. While these germline mutations represent a fundamental source of new genetic diversity, they can also create deleterious alleles that impact fitness. Whereas the rate and patterns of point mutations in the human germline are now well understood, far less is known about the frequency and features that impact de novo structural variants (dnSVs). We report a family-based study of germline mutations among 9,599 human genomes from 33 multigenerational CEPH-Utah families and 2,384 families from the Simons Foundation Autism Research Initiative. We find that de novo structural mutations detected by alignment-based, short-read WGS occur at an overall rate of at least 0.160 events per genome in unaffected individuals, and we observe a significantly higher rate (0.206 per genome) in ASD-affected individuals. In both probands and unaffected samples, nearly 73% of de novo structural mutations arose in paternal gametes, and we predict most de novo structural mutations to be caused by mutational mechanisms that do not require sequence homology. After multiple testing correction, we did not observe a statistically significant correlation between parental age and the rate of de novo structural variation in offspring. These results highlight that a spectrum of mutational mechanisms contribute to germline structural mutations and that these mechanisms most likely have markedly different rates and selective pressures than those leading to point mutations.

ObjectiveBoth the gut microbiome and host genetics are known to play significant roles in the pathogenesis of IBD. However, the interaction between these two factors and its implications in the aetiology of IBD remain underexplored. Here, we report on the influence of host genetics on the gut microbiome in IBD.DesignTo evaluate the impact of host genetics on the gut microbiota of patients with IBD, we combined whole exome sequencing of the host genome and whole genome shotgun sequencing of 1464 faecal samples from 525 patients with IBD and 939 population-based controls. We followed a four-step analysis: (1) exome-wide microbial quantitative trait loci (mbQTL) analyses, (2) a targeted approach focusing on IBD-associated genomic regions and protein truncating variants (PTVs, minor allele frequency (MAF) >5%), (3) gene-based burden tests on PTVs with MAF ResultsWe identified 12 mbQTLs, including variants in the IBD-associated genes IL17REL, MYRF, SEC16A and WDR78. For example, the decrease of the pathway acetyl-coenzyme A biosynthesis, which is involved in short chain fatty acids production, was associated with variants in the gene MYRF (false discovery rate CYP2D6, GPR151 and CD160 genes. These genes are known for their function in the immune system. Moreover, interaction analyses confirmed previously known IBD disease-specific mbQTLs in TNFSF15.ConclusionThis study highlights that both common and rare genetic variants affecting the immune system are key factors in shaping the gut microbiota in the context of IBD and pinpoints towards potential mechanisms for disease treatment.

Over the last two years, remarkable gene discovery efforts have implicated disruption of pathways involving gene regulatory functions and neuronal processes in autism spectrum disorder (ASD), and more broadly defined neurodevelopmental disorders (NDDs). Functional studies in the developing brain and across cell types demonstrate that the spatiotemporal expression patterns of many of these genes coalesce on subnetworks with distinct developmental trajectories. Here, we review the convergent biological processes derived from gene discovery and functional genomics in ASD and NDD from 2018-2020. We further probe the mechanistic insights that suggest these frequently perturbed pathways are interconnected and, ultimately, converge on specific functional deficits in human neurodevelopment.

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.To dissect the contribution of induced excitatory neurons (iENs) versus GABAergic (gamma-aminobutyric acidergic) neurons (iGNs) to SCZ pathophysiology, we induced iNs from CRISPR (clustered regularly interspaced short palindromic repeats)-Cas9 isogenic and SCZ patient-derived induced pluripotent stem cells and analyzed SCZ-related phenotypes in iEN monocultures and iEN/iGN cocultures.In iEN/iGN cocultures, neuronal firing and synchrony were reduced at later, but not earlier, stages of in vitro development. These were fully recapitulated in iEN monocultures, indicating a primary role for iENs. Moreover, isogenic iENs showed reduced dendrite length and deficits in calcium handling. iENs from 16p11.2 duplication-carrying patients with SCZ displayed overlapping deficits in network synchrony, dendrite outgrowth, and calcium handling. Transcriptomic analysis of both iEN cohorts revealed molecular markers of disease related to the glutamatergic synapse, neuroarchitecture, and calcium regulation.Our results indicate the presence of 16p11.2 duplication-dependent alterations in SCZ patient-derived iENs. Transcriptomics and cellular phenotyping reveal overlap between isogenic and patient-derived iENs, suggesting a central role of glutamatergic, morphological, and calcium dysregulation in 16p11.2 duplication-mediated pathogenesis. Moreover, excitatory dysfunction during early neurodevelopment is implicated as the basis of SCZ pathogenesis in 16p11.2 duplication carriers. Our results support network synchrony and calcium handling as outcomes directly linked to this genetic risk variant.

The influence of genetic background on driver mutations is well established; however, the mechanisms by which the background interacts with Mendelian loci remain unclear. We performed a systematic secondary-variant burden analysis of two independent cohorts of patients with Bardet-Biedl syndrome (BBS) with known recessive biallelic pathogenic mutations in one of 17 BBS genes for each individual. We observed a significant enrichment of trans-acting rare nonsynonymous secondary variants in patients with BBS compared with either population controls or a cohort of individuals with a non-BBS diagnosis and recessive variants in the same gene set. Strikingly, we found a significant over-representation of secondary alleles in chaperonin-encoding genes-a finding corroborated by the observation of epistatic interactions involving this complex in vivo. These data indicate a complex genetic architecture for BBS that informs the biological properties of disease modules and presents a model for secondary-variant burden analysis in recessive disorders.

Structural variants (SVs) contribute to many disorders, yet, functionally annotating them remains a major challenge. Here, we integrate SVs with RNA-sequencing from human post-mortem brains to quantify their dosage and regulatory effects. We show that genic and regulatory SVs exist at significantly lower frequencies than intergenic SVs. Functional impact of copy number variants (CNVs) stems from both the proportion of genic and regulatory content altered and loss-of-function intolerance of the gene. We train a linear model to predict expression effects of rare CNVs and use it to annotate regulatory disruption of CNVs from 14,891 independent genome-sequenced individuals. Pathogenic deletions implicated in neurodevelopmental disorders show significantly more extreme regulatory disruption scores and if rank ordered would be prioritized higher than using frequency or length alone. This work shows the deleteriousness of regulatory SVs, particularly those altering CTCF sites and provides a simple approach for functionally annotating the regulatory consequences of CNVs., Structural variants (SVs) contribute to the genetic architecture of many brain-related disorders. Here, the authors integrate SV calls from genome sequencing (n = 755) with RNA-seq data (n = 629) from post-mortem dorsal lateral prefrontal cortex to annotate the gene regulatory effects of SVs in the human brain and their potential to contribute to disease.

Background Self-reported consanguinity is associated with risk for schizophrenia (SZ) in several inbred populations, but estimates using DNA-based coefficients of inbreeding are unavailable. Further, it is not known whether recessively inherited risk mutations can be identified through homozygosity by descent (HBD) mapping. Methods We studied self-reported and DNA-based estimates of inbreeding among Egyptian patients with SZ (n = 421, DSM IV criteria) and adult controls without psychosis (n = 301), who were evaluated using semi-structured diagnostic interview schedules and genotyped using the Illumina Infinium PsychArray. Following quality control checks, coefficients of inbreeding (F) and regions of homozygosity (ROH) were estimated using PLINK software for HBD analysis. Exome sequencing was conducted in selected cases. Results Inbreeding was associated with schizophrenia based on self-reported consanguinity (χ2 = 4.506, 1 df, p = 0.034) and DNA-based estimates for inbreeding (F); the latter with a significant F × age interaction (β = 32.34, p = 0.0047). The association was most notable among patients older than age 40 years. Eleven ROH were over-represented in cases on chromosomes 1, 3, 6, 11, and 14; all but one region is novel for schizophrenia risk. Exome sequencing identified six recessively-acting genes in ROH with loss-of-function variants; one of which causes primary hereditary microcephaly. Conclusions We propose consanguinity as an age-dependent risk factor for SZ in Egypt. HBD mapping is feasible for SZ in adequately powered samples.

Background Tuberous sclerosis complex (TSC) is a neurodevelopmental disorder with frequent occurrence of epilepsy, autism spectrum disorder (ASD), intellectual disability (ID), and tumors in multiple organs. The aberrant activation of mTORC1 in TSC has led to treatment with mTORC1 inhibitor rapamycin as a lifelong therapy for tumors, but TSC-associated neurocognitive manifestations remain unaffected by rapamycin. Methods Here, we generated patient-specific, induced pluripotent stem cells (iPSCs) from a TSC patient with a heterozygous, germline, nonsense mutation in exon 15 of TSC1 and established an isogenic set of heterozygous (Het), null and corrected wildtype (Corr-WT) iPSCs using CRISPR/Cas9-mediated gene editing. We differentiated these iPSCs into neural progenitor cells (NPCs) and examined neurodevelopmental phenotypes, signaling and changes in gene expression by RNA-seq. Results Differentiated NPCs revealed enlarged cell size in TSC1-Het and Null NPCs, consistent with mTORC1 activation. TSC1-Het and Null NPCs also revealed enhanced proliferation and altered neurite outgrowth in a genotype-dependent manner, which was not reversed by rapamycin. Transcriptome analyses of TSC1-NPCs revealed differentially expressed genes that display a genotype-dependent linear response, i.e., genes upregulated/downregulated in Het were further increased/decreased in Null. In particular, genes linked to ASD, epilepsy, and ID were significantly upregulated or downregulated warranting further investigation. In TSC1-Het and Null NPCs, we also observed basal activation of ERK1/2, which was further activated upon rapamycin treatment. Rapamycin also increased MNK1/2-eIF4E signaling in TSC1-deficient NPCs. Conclusion MEK-ERK and MNK-eIF4E pathways regulate protein translation, and our results suggest that aberrant translation distinct in TSC1/2-deficient NPCs could play a role in neurodevelopmental defects. Our data showing upregulation of these signaling pathways by rapamycin support a strategy to combine a MEK or a MNK inhibitor with rapamycin that may be superior for TSC-associated CNS defects. Importantly, our generation of isogenic sets of NPCs from TSC patients provides a valuable platform for translatome and large-scale drug screening studies. Overall, our studies further support the notion that early developmental events such as NPC proliferation and initial process formation, such as neurite number and length that occur prior to neuronal differentiation, represent primary events in neurogenesis critical to disease pathogenesis of neurodevelopmental disorders such as ASD.

Familial dysautonomia (FD), a hereditary sensory and autonomic neuropathy, is caused by a mutation in the Elongator complex protein 1 (ELP1) gene that leads to a tissue-specific reduction of ELP1 protein. Our work to generate a phenotypic mouse model for FD headed to the discovery that homozygous deletion of the mouse Elp1 gene leads to embryonic lethality prior to mid-gestation. Given that FD is caused by a reduction, not loss, of ELP1, we generated two new mouse models by introducing different copy numbers of the human FD ELP1 transgene into the Elp1 knockout mouse (Elp1−/−) and observed that human ELP1 expression rescues embryonic development in a dose-dependent manner. We then conducted a comprehensive transcriptome analysis in mouse embryos to identify genes and pathways whose expression correlates with the amount of ELP1. We found that ELP1 is essential for the expression of genes responsible for nervous system development. Further, gene length analysis of the differentially expressed genes showed that the loss of Elp1 mainly impacts the expression of long genes and that by gradually restoring Elongator, their expression is progressively rescued. Finally, through evaluation of co-expression modules, we identified gene sets with unique expression patterns that depended on ELP1 expression.

Background Genomic structural variants (SVs) can affect many genes and regulatory elements. Therefore, the molecular mechanisms driving the phenotypes of patients carrying de novo SVs are frequently unknown. Methods We applied a combination of systematic experimental and bioinformatic methods to improve the molecular diagnosis of 39 patients with multiple congenital abnormalities and/or intellectual disability harboring apparent de novo SVs, most with an inconclusive diagnosis after regular genetic testing. Results In 7 of these cases (18%), whole-genome sequencing analysis revealed disease-relevant complexities of the SVs missed in routine microarray-based analyses. We developed a computational tool to predict the effects on genes directly affected by SVs and on genes indirectly affected likely due to the changes in chromatin organization and impact on regulatory mechanisms. By combining these functional predictions with extensive phenotype information, candidate driver genes were identified in 16/39 (41%) patients. In 8 cases, evidence was found for the involvement of multiple candidate drivers contributing to different parts of the phenotypes. Subsequently, we applied this computational method to two cohorts containing a total of 379 patients with previously detected and classified de novo SVs and identified candidate driver genes in 189 cases (50%), including 40 cases whose SVs were previously not classified as pathogenic. Pathogenic position effects were predicted in 28% of all studied cases with balanced SVs and in 11% of the cases with copy number variants. Conclusions These results demonstrate an integrated computational and experimental approach to predict driver genes based on analyses of WGS data with phenotype association and chromatin organization datasets. These analyses nominate new pathogenic loci and have strong potential to improve the molecular diagnosis of patients with de novo SVs.

Reciprocal copy number variations (CNVs) of 16p11.2 are associated with a wide spectrum of neuropsychiatric and neurodevelopmental disorders. Here, we use human induced pluripotent stem cells (iPSCs)-derived dopaminergic (DA) neurons carrying CNVs of 16p11.2 duplication (16pdup) and 16p11.2 deletion (16pdel), engineered using CRISPR-Cas9. We show that 16pdel iPSC-derived DA neurons have increased soma size and synaptic marker expression compared to isogenic control lines, while 16pdup iPSC-derived DA neurons show deficits in neuronal differentiation and reduced synaptic marker expression. The 16pdel iPSC-derived DA neurons have impaired neurophysiological properties. The 16pdel iPSC-derived DA neuronal networks are hyperactive and have increased bursting in culture compared to controls. We also show that the expression of RHOA is increased in the 16pdel iPSC-derived DA neurons and that treatment with a specific RHOA-inhibitor, Rhosin, rescues the network activity of the 16pdel iPSC-derived DA neurons. Our data suggest that 16p11.2 deletion-associated iPSC-derived DA neuron hyperactivation can be rescued by RHOA inhibition., 16p11.2 CNVs are associated with neurodevelopmental disorders. Here, the authors show that 16p11.2 deletion is associated with hyperactivation of human iPSC-derived dopaminergic neuron networks and is rescued by RHOA inhibition in vitro.

The 1000 Genomes Project (1kGP) is the largest fully open resource of whole-genome sequencing (WGS) data consented for public distribution without access or use restrictions. The final, phase 3 release of the 1kGP included 2,504 unrelated samples from 26 populations and was based primarily on low-coverage WGS. Here, we present a high-coverage 3,202-sample WGS 1kGP resource, which now includes 602 complete trios, sequenced to a depth of 30X using Illumina. We performed single-nucleotide variant (SNV) and short insertion and deletion (INDEL) discovery and generated a comprehensive set of structural variants (SVs) by integrating multiple analytic methods through a machine learning model. We show gains in sensitivity and precision of variant calls compared to phase 3, especially among rare SNVs as well as INDELs and SVs spanning frequency spectrum. We also generated an improved reference imputation panel, making variants discovered here accessible for association studies.

Background The Autism Sequencing Consortium identified 102 high-confidence autism spectrum disorder (ASD) genes, showing that individuals with ASD and with potentially damaging single nucleotide variation (pdSNV) in these genes had lower cognitive levels and delayed age at walking, when compared to ASD participants without pdSNV. Here, we made use of a Swedish sample of individuals with ASD (called PAGES, for Population-Based Autism Genetics & Environment Study) to evaluate the frequency of pdSNV and their impact on medical and psychiatric phenotypes, using an epidemiological frame and universal health reporting. We then combine findings with those for potentially damaging copy number variation (pdCNV). Methods SNV and CNV calls were generated from whole-exome sequencing and chromosome microarray data, respectively. Birth and medical register data were used to collect phenotypes. Results Of 808 individuals assessed by sequencing, 69 (9%) had pdSNV in the 102 ASC genes, and 144 (18%) had pdSNV in the 102 ASC genes or in a larger set of curated neurodevelopmental genes (from the Deciphering Developmental Disorders study, the gene2phenotype database, and the Radboud University gene lists). Three or more individuals had pdSNV in GRIN2B, POGZ, SATB1, DYNC1H1, SCN8A, or CREBBP. In comparison, out of the 996 individuals from whom CNV were called, 105 (11%) carried one or more pdCNV, including four or more individuals with CNV in the recurrent 15q11q13, 22q11.2, and 16p11.2 loci. Carriers of pdSNV were more likely to have intellectual disability (ID) and epilepsy, while carriers of pdCNV showed increased rates of congenital anomalies and scholastic skill disorders. Carriers of either pdSNV or pdCNV were more likely to have ID, scholastic skill disorders, and epilepsy. Limitations The cohort only included individuals with autistic disorder, the more severe form of ASD, and phenotypes are defined from medical registers. Not all genes studied are definitively ASD genes, and we did not have de novo information to aid in classification. Conclusions In this epidemiological sample, rare pdSNV were more common than pdCNV and the combined yield of potentially damaging variation was substantial at 27%. The results provide compelling rationale for the use of high-throughout sequencing as part of routine clinical workup for ASD and support the development of precision medicine in ASD.

Megabase-scale intervals of active, gene-rich and inactive, gene-poor chromatin are known to segregate, forming the A and B compartments. Fine mapping of the contents of these A and B compartments has been hitherto impossible, owing to the extraordinary sequencing depths required to distinguish between the long-range contact patterns of individual loci, and to the computational complexity of the associated calculations. Here, we generate the largest published in situ Hi-C map to date, spanning 33 billion contacts. We also develop a computational method, dubbed PCA of Sparse, SUper Massive Matrices (POSSUMM), that is capable of efficiently calculating eigenvectors for sparse matrices with millions of rows and columns. Applying POSSUMM to our Hi-C dataset makes it possible to assign loci to the A and B compartment at 500 bp resolution. We find that loci frequently alternate between compartments as one moves along the contour of the genome, such that the median compartment interval is only 12.5 kb long. Contrary to the findings in coarse-resolution compartment profiles, we find that individual genes are not uniformly positioned in either the A compartment or the B compartment. Instead, essentially all (95%) active gene promoters localize in the A compartment, but the likelihood of localizing in the A compartment declines along the body of active genes, such that the transcriptional termini of long genes (>60 kb) tend to localize in the B compartment. Similarly, nearly all active enhancers elements (95%) localize in the A compartment, even when the flanking sequences are comprised entirely of inactive chromatin and localize in the B compartment. These results are consistent with a model in which DNA-bound regulatory complexes give rise to phase separation at the scale of individual DNA elements.

Genetic variants that inactivate protein-coding genes are a powerful source of information about the phenotypic consequences of gene disruption: genes that are crucial for the function of an organism will be depleted of such variants in natural populations, whereas non-essential genes will tolerate their accumulation. However, predicted loss-of-function variants are enriched for annotation errors, and tend to be found at extremely low frequencies, so their analysis requires careful variant annotation and very large sample sizes

Triggers of innate immune signaling in the CNS of patients with amyotrophic lateral sclerosis and frontotemporal degeneration (ALS/FTD) remain elusive. We report the presence of cytoplasmic double-stranded RNA (cdsRNA), an established trigger of innate immunity, in ALS-FTD brains carrying C9ORF72 intronic hexanucleotide expansions that included genomically encoded expansions of the G4C2 repeat sequences. The presence of cdsRNA in human brains was coincident with cytoplasmic TAR DNA binding protein 43 (TDP-43) inclusions, a pathologic hallmark of ALS/FTD. Introducing cdsRNA into cultured human neural cells induced type I interferon (IFN-I) signaling and death that was rescued by FDA-approved JAK inhibitors. In mice, genomically encoded dsRNAs expressed exclusively in a neuronal class induced IFN-I and death in connected neurons non-cell-autonomously. Our findings establish that genomically encoded cdsRNAs trigger sterile, viral-mimetic IFN-I induction and propagated death within neural circuits and may drive neuroinflammation and neurodegeneration in patients with ALS/FTD.

Dystonia is a neurologic disorder associated with an increasingly large number of variants in many genes, resulting in characteristic disturbances in volitional movement. Dissecting the relationships between these mutations and their functional outcomes is a critical step in understanding the key pathways that drive dystonia pathogenesis. Here we established a pipeline for characterizing an allelic series of dystonia-specific mutations in isogenic induced pluripotent stem cells (iPSCs). We used this strategy to investigate the molecular consequences of variation in THAP1, which encodes a transcription factor that has been linked to neural differentiation. Multiple pathogenic mutations that have been associated with dystonia cluster within distinct THAP1 functional domains and are predicted to alter its DNA binding properties and/or protein interactions differently, yet the relative impact of these varied changes on molecular signatures and neural deficits is unclear. To determine the effects of these mutations on THAP1 transcriptional activity, we engineered an allelic series of eight mutations in a common iPSC background and differentiated these lines into a panel of near-isogenic neural stem cells (n = 94 lines). Transcriptome profiling of these neural derivatives followed by joint analysis of the most robust individual signatures across mutations identified a convergent pattern of dysregulated genes functionally related to neurodevelopment, lysosomal lipid metabolism, and myelin. Based on these observations, we examined mice bearing Thap1-disruptive alleles and detected significant changes in myelin gene expression and reduction of myelin structural integrity relative to tissue from control mice. These results suggest that deficits in neurodevelopment and myelination are common consequences of dystonia-associated THAP1 mutations and highlight the potential role of neuron-glial interactions in the pathogenesis of dystonia.

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., Drugs that modify RNA splicing are promising treatments for many genetic diseases. Here the authors show that deep learning strategies can predict drug targets, strongly supporting the use of in silico approaches to expand the therapeutic potential of drugs that modulate RNA splicing.

Context Kallmann syndrome (KS) is a rare, genetically heterogeneous Mendelian disorder. Structural defects in KS patients have helped define the genetic architecture of gonadotropin-releasing hormone (GnRH) neuronal development in this condition. Objective Examine the functional role a novel structural defect affecting a long noncoding RNA (lncRNA), RMST, found in a KS patient. Design Whole genome sequencing, induced pluripotent stem cells and derived neural crest cells (NCC) from the KS patient were contrasted with controls. Setting The Harvard Reproductive Sciences Center, Massachusetts General Hospital Center for Genomic Medicine, and Singapore Genome Institute. Patient A KS patient with a unique translocation, t(7;12)(q22;q24). Interventions/Main Outcome Measure/Results A novel translocation was detected affecting the lncRNA, RMST, on chromosome 12 in the absence of any other KS mutations. Compared with controls, the patient’s induced pluripotent stem cells and NCC provided functional information regarding RMST. Whereas RMST expression increased during NCC differentiation in controls, it was substantially reduced in the KS patient’s NCC coincident with abrogated NCC morphological development and abnormal expression of several “downstream” genes essential for GnRH ontogeny (SOX2, PAX3, CHD7, TUBB3, and MKRN3). Additionally, an intronic single nucleotide polymorphism in RMST was significantly implicated in a genome-wide association study associated with age of menarche. Conclusions A novel deletion in RMST implicates the loss of function of a lncRNA as a unique cause of KS and suggests it plays a critical role in the ontogeny of GnRH neurons and puberty.

Familial dysautonomia (FD) is a recessive neurodegenerative disease caused by a splice mutation in Elongator complex protein 1 (ELP1, also known as IKBKAP); this mutation leads to variable skipping of exon 20 and to a drastic reduction of ELP1 in the nervous system. Clinically, many of the debilitating aspects of the disease are related to a progressive loss of proprioception; this loss leads to severe gait ataxia, spinal deformities, and respiratory insufficiency due to neuromuscular incoordination. There is currently no effective treatment for FD, and the disease is ultimately fatal. The development of a drug that targets the underlying molecular defect provides hope that the drastic peripheral neurodegeneration characteristic of FD can be halted. We demonstrate herein that the FD mouse TgFD9;Ikbkap(Δ20/flox) recapitulates the proprioceptive impairment observed in individuals with FD, and we provide the in vivo evidence that postnatal correction, promoted by the small molecule kinetin, of the mutant ELP1 splicing can rescue neurological phenotypes in FD. Daily administration of kinetin starting at birth improves sensory-motor coordination and prevents the onset of spinal abnormalities by stopping the loss of proprioceptive neurons. These phenotypic improvements correlate with increased amounts of full-length ELP1 mRNA and protein in multiple tissues, including in the peripheral nervous system (PNS). Our results show that postnatal correction of the underlying ELP1 splicing defect can rescue devastating disease phenotypes and is therefore a viable therapeutic approach for persons with FD.

Background: The relative prevalence and clinical importance of monogenic mutations related to familial hypercholesterolemia and of high polygenic score (cumulative impact of many common variants) pathways for early-onset myocardial infarction remain uncertain. Whole-genome sequencing enables simultaneous ascertainment of both monogenic mutations and polygenic score for each individual. Methods: We performed deep-coverage whole-genome sequencing of 2081 patients from 4 racial subgroups hospitalized in the United States with early-onset myocardial infarction (age ≤55 years) recruited with a 2:1 female-to-male enrollment design. We compared these genomes with those of 3761 population-based control subjects. We first identified individuals with a rare, monogenic mutation related to familial hypercholesterolemia. Second, we calculated a recently developed polygenic score of 6.6 million common DNA variants to quantify the cumulative susceptibility conferred by common variants. We defined high polygenic score as the top 5% of the control distribution because this cutoff has previously been shown to confer similar risk to that of familial hypercholesterolemia mutations. Results: The mean age of the 2081 patients presenting with early-onset myocardial infarction was 48 years, and 66% were female. A familial hypercholesterolemia mutation was present in 36 of these patients (1.7%) and was associated with a 3.8-fold (95% CI, 2.1-6.8; P3-fold increased odds of early-onset myocardial infarction. However, high polygenic score has a 10-fold higher prevalence among patients presents with early-onset myocardial infarction. Clinical Trial Registration: URL: https://www.clinicaltrials.gov. Unique identifier: NCT00597922. ©2019, National Institutes of Health (no. TR001100), American Heart Association (no. 17SDG33680041), National Human Genome Research Institute (no. 5UM1HG008895)