Search

There is evidence that transcription factor (TF) encoding genes, which temporally control development in multiple cell types, can have tens of enhancers that regulate their expression. The NR2F1 TF developmentally promotes caudal and ventral cortical regional fates. Here, we epigenomically compared the activity of Nr2f1's enhancers during mouse cortical development with their activity in a transgenic assay. We identified at least six that are likely to be important in prenatal cortical development, with three harboring de novo mutants identified in ASD individuals. We chose to study the function of two of the most robust enhancers by deleting them singly or together. We found that they have distinct and overlapping functions in driving Nr2f1's regional and laminar expression in the developing cortex. Thus, these two enhancers, probably in combination with the others that we defined epigenetically, precisely tune Nr2f1's regional, cell type, and temporal expression during corticogenesis.

Many autism spectrum disorder (ASD)-associated genes act as transcriptional regulators (TRs). Chromatin immunoprecipitation sequencing (ChIP-seq) was used to identify the regulatory targets of ARID1B, BCL11A, FOXP1, TBR1, and TCF7L2, ASD-associated TRs in the developing human and mouse cortex. These TRs shared substantial overlap in the binding sites, especially within open chromatin. The overlap within a promoter region, 1-2,000 bp upstream of the transcription start site, was highly predictive of brain-expressed genes. This signature was observed in 96 out of 102 ASD-associated genes. In vitro CRISPRi against ARID1B and TBR1 delineated downstream convergent biology in mouse cortical cultures. After 8 days, NeuN+ and CALB+ cells were decreased, GFAP+ cells were increased, and transcriptomic signatures correlated with the postmortem brain samples from individuals with ASD. We suggest that functional convergence across five ASD-associated TRs leads to shared neurodevelopmental outcomes of haploinsufficient disruption.

Identification of Alzheimer's disease (AD) onset risk can facilitate interventions before irreversible disease progression. We demonstrate that electronic health records from the University of California, San Francisco, followed by knowledge networks (for example, SPOKE) allow for (1) prediction of AD onset and (2) prioritization of biological hypotheses, and (3) contextualization of sex dimorphism. We trained random forest models and predicted AD onset on a cohort of 749 individuals with AD and 250,545 controls with a mean area under the receiver operating characteristic of 0.72 (7 years prior) to 0.81 (1 day prior). We further harnessed matched cohort models to identify conditions with predictive power before AD onset. Knowledge networks highlight shared genes between multiple top predictors and AD (for example, APOE, ACTB, IL6 and INS). Genetic colocalization analysis supports AD association with hyperlipidemia at the APOE locus, as well as a stronger female AD association with osteoporosis at a locus near MS4A6A. We therefore show how clinical data can be utilized for early AD prediction and identification of personalized biological hypotheses.

CRISPR-based gene activation (CRISPRa) is a strategy for upregulating gene expression by targeting promoters or enhancers in a tissue/cell-type specific manner. Here, we describe an experimental framework that combines highly multiplexed perturbations with single-cell RNA sequencing (sc-RNA-seq) to identify cell-type-specific, CRISPRa-responsive cis-regulatory elements and the gene(s) they regulate. Random combinations of many gRNAs are introduced to each of many cells, which are then profiled and partitioned into test and control groups to test for effect(s) of CRISPRa perturbations of both enhancers and promoters on the expression of neighboring genes. Applying this method to a library of 493 gRNAs targeting candidate cis-regulatory elements in both K562 cells and iPSC-derived excitatory neurons, we identify gRNAs capable of specifically upregulating intended target genes and no other neighboring genes within 1 Mb, including gRNAs yielding upregulation of six autism spectrum disorder (ASD) and neurodevelopmental disorder (NDD) risk genes in neurons. A consistent pattern is that the responsiveness of individual enhancers to CRISPRa is restricted by cell type, implying a dependency on either chromatin landscape and/or additional trans-acting factors for successful gene activation. The approach outlined here may facilitate large-scale screens for gRNAs that activate genes in a cell type-specific manner.

ObjectiveThe aim of this study was to catalog and evaluate response biomarkers correlated with autism spectrum disorder (ASD) symptoms to improve clinical trials.MethodsA systematic review of MEDLINE, Embase, and Scopus was conducted in April 2020. Seven criteria were applied to focus on original research that includes quantifiable response biomarkers measured alongside ASD symptoms. Interventional studies or human studies that assessed the correlation between biomarkers and ASD-related behavioral measures were included.ResultsA total of 5,799 independent records yielded 280 articles for review that reported on 940 biomarkers, 755 of which were unique to a single publication. Molecular biomarkers were the most frequently assayed, including cytokines, growth factors, measures of oxidative stress, neurotransmitters, and hormones, followed by neurophysiology (e.g., EEG and eye tracking), neuroimaging (e.g., functional MRI), and other physiological measures. Studies were highly heterogeneous, including in phenotypes, demographic characteristics, tissues assayed, and methods for biomarker detection. With a median total sample size of 64, almost all of the reviewed studies were only powered to identify biomarkers with large effect sizes. Reporting of individual-level values and summary statistics was inconsistent, hampering mega- and meta-analysis. Biomarkers assayed in multiple studies yielded mostly inconsistent results, revealing a "replication crisis."ConclusionsThere is currently no response biomarker with sufficient evidence to inform ASD clinical trials. This review highlights methodological imperatives for ASD biomarker research necessary to make definitive progress: consistent experimental design, correction for multiple comparisons, formal replication, sharing of sample-level data, and preregistration of study designs. Systematic "big data" analyses of multiple potential biomarkers could accelerate discovery.

Autism spectrum disorder (ASD) is a common, complex, and highly heritable condition with contributions from both common and rare genetic variations. While disruptive, rare variants in protein-coding regions clearly contribute to symptoms, the role of rare non-coding remains unclear. Variants in these regions, including promoters, can alter downstream RNA and protein quantity; however, the functional impacts of specific variants observed in ASD cohorts remain largely uncharacterized. Here, we analyzed 3600 de novo mutations in promoter regions previously identified by whole-genome sequencing of autistic probands and neurotypical siblings to test the hypothesis that mutations in cases have a greater functional impact than those in controls. We leveraged massively parallel reporter assays (MPRAs) to detect transcriptional consequences of these variants in neural progenitor cells and identified 165 functionally high confidence de novo variants (HcDNVs). While these HcDNVs are enriched for markers of active transcription, disruption to transcription factor binding sites, and open chromatin, we did not identify differences in functional impact based on ASD diagnostic status.

Posttranscriptional RNA modifications by adenosine-to-inosine (A-to-I) editing are abundant in the brain, yet elucidating functional sites remains challenging. To bridge this gap, we investigate spatiotemporal and genetically regulated A-to-I editing sites across prenatal and postnatal stages of human brain development. More than 10,000 spatiotemporally regulated A-to-I sites were identified that occur predominately in 3' UTRs and introns, as well as 37 sites that recode amino acids in protein coding regions with precise changes in editing levels across development. Hyper-edited transcripts are also enriched in the aging brain and stabilize RNA secondary structures. These features are conserved in murine and non-human primate models of neurodevelopment. Finally, thousands of cis-editing quantitative trait loci (edQTLs) were identified with unique regulatory effects during prenatal and postnatal development. Collectively, this work offers a resolved atlas linking spatiotemporal variation in editing levels to genetic regulatory effects throughout distinct stages of brain maturation.

Some individuals with autism spectrum disorder (ASD) carry functional mutations rarely observed in the general population. We explored the genes disrupted by these variants from joint analysis of protein-truncating variants (PTVs), missense variants and copy number variants (CNVs) in a cohort of 63,237 individuals. We discovered 72 genes associated with ASD at false discovery rate (FDR) ≤ 0.001 (185 at FDR ≤ 0.05). De novo PTVs, damaging missense variants and CNVs represented 57.5%, 21.1% and 8.44% of association evidence, while CNVs conferred greatest relative risk. Meta-analysis with cohorts ascertained for developmental delay (DD) (n = 91,605) yielded 373 genes associated with ASD/DD at FDR ≤ 0.001 (664 at FDR ≤ 0.05), some of which differed in relative frequency of mutation between ASD and DD cohorts. The DD-associated genes were enriched in transcriptomes of progenitor and immature neuronal cells, whereas genes showing stronger evidence in ASD were more enriched in maturing neurons and overlapped with schizophrenia-associated genes, emphasizing that these neuropsychiatric disorders may share common pathways to risk.

lternative DNA conformations, termed non-B DNA structures, can affect transcription, but the underlying mechanisms and their functional impact have not been systematically characterized. Here, we used computational genomic analyses coupled with massively parallel reporter assays (MPRAs) to show that certain non-B DNA structures have a substantial effect on gene expression. Genomic analyses found that non-B DNA structures at promoters harbor an excess of germline variants. Analysis of multiple MPRAs, including a promoter library specifically designed to perturb non-B DNA structures, functionally validated that Z-DNA can significantly affect promoter activity. We also observed that biophysical properties of non-B DNA motifs, such as the length of Z-DNA motifs and the orientation of G-quadruplex structures relative to transcriptional direction, have a significant effect on promoter activity. Combined, their higher mutation rate and functional effect on transcription implicate a subset of non-B DNA motifs as major drivers of human gene-expression-associated phenotypes.

Knowledge representation and reasoning (KR&R) has been successfully implemented in many fields to enable computers to solve complex problems with AI methods. However, its application to biomedicine has been lagging in part due to the daunting complexity of molecular and cellular pathways that govern human physiology and pathology. In this article we describe concrete uses of SPOKE, an open knowledge network that connects curated information from 37 specialized and human-curated databases into a single property graph, with 3 million nodes and 15 million edges to date. Applications discussed in this article include drug discovery, COVID-19 research and chronic disease diagnosis and management.

White matter lesions (WML) underlie multiple brain disorders, and automatic WML segmentation is crucial to evaluate the natural disease course and effectiveness of clinical interventions, including drug discovery. Although recent research has achieved tremendous progress in WML segmentation, accurate detection of subtle WML present early in the disease course remains particularly challenging. Here we propose an approach to automatic WML segmentation of mild WML loads using an intensity standardisation technique, gray level co-occurrence matrix (GLCM) embedded clustering technique, and random forest (RF) classifier to extract texture features and identify morphology specific to true WML. We precisely define their boundaries through a local outlier factor (LOF) algorithm that identifies edge pixels by local density deviation relative to its neighbors. The automated approach was validated on 32 human subjects, demonstrating strong agreement and correlation (excluding one outlier) with manual delineation by a neuroradiologist through Intra-Class Correlation (ICC = 0.881, 95% CI 0.769, 0.941) and Pearson correlation (r = 0.895, p-value

Neurodevelopmental neuropsychiatric disorders, such as autism spectrum disorder and schizophrenia, have strong genetic risk components, but the underlying mechanisms have proven difficult to decipher. Rare, high-risk variants may offer an opportunity to delineate the biological mechanisms responsible more clearly for more common idiopathic diseases. Indeed, different rare variants can cause the same behavioral phenotype, demonstrating genetic heterogeneity, while the same rare variant can cause different behavioral phenotypes, demonstrating variable expressivity. These observations suggest convergent underlying biological and neurological mechanisms; identification of these mechanisms may ultimately reveal new therapeutic targets. At the 2021 Keystone eSymposium "Neuropsychiatric and Neurodevelopmental Disorders: Harnessing Rare Variants" a panel of experts in the field described significant progress in genomic discovery and human phenotyping and raised several consistent issues, including the need for detailed natural history studies of rare disorders, the challenges in cohort recruitment, and the importance of viewing phenotypes as quantitative traits that are impacted by rare variants.

BackgroundGenetic variants in the voltage-gated sodium channels SCN1A, SCN2A, SCN3A, and SCN8A are leading causes of epilepsy, developmental delay, and autism spectrum disorder. The mRNA splicing patterns of all four genes vary across development in the rodent brain, including mutually exclusive copies of the fifth protein-coding exon detected in the neonate (5N) and adult (5A). A second pair of mutually exclusive exons is reported in SCN8A only (18N and 18A). We aimed to quantify the expression of individual exons in the developing human brain.MethodsRNA-seq data from 783 human brain samples across development were analyzed to estimate exon-level expression. Developmental changes in exon utilization were validated by assessing intron splicing. Exon expression was also estimated in RNA-seq data from 58 developing mouse neocortical samples.ResultsIn the mature human neocortex, exon 5A is consistently expressed at least 4-fold higher than exon 5N in all four genes. For SCN2A, SCN3A, and SCN8A, a brain-wide synchronized 5N to 5A transition occurs between 24 post-conceptual weeks (2nd trimester) and 6 years of age. In mice, the equivalent 5N to 5A transition begins at or before embryonic day 15.5. In SCN8A, over 90% of transcripts in the mature human cortex include exon 18A. Early in fetal development, most transcripts include 18N or skip both 18N and 18A, with a transition to 18A inclusion occurring from 13 post-conceptual weeks to 6 months of age. No other protein-coding exons showed comparably dynamic developmental trajectories.ConclusionsExon usage in SCN1A, SCN2A, SCN3A, and SCN8A changes dramatically during human brain development. These splice isoforms, which alter the biophysical properties of the encoded channels, may account for some of the observed phenotypic differences across development and between specific variants. Manipulation of the proportion of splicing isoforms at appropriate stages of development may act as a therapeutic strategy for specific mutations or even epilepsy in general.

BackgroundRecent large-scale initiatives have led to systematically collected phenotypic data for several rare genetic conditions implicated in autism spectrum disorder (ASD). The onset of developmentally expected skills (e.g. walking, talking) serve as readily quantifiable aspects of the behavioral phenotype. This study's aims were: (a) describe the distribution of ages of attainment of gross motor and expressive language milestones in several rare genetic conditions, and (b) characterize the likelihood of delays in these conditions compared with idiopathic ASD.MethodsParticipants aged 3 years and older were drawn from two Simons Foundation Autism Research Initiative registries that employed consistent phenotyping protocols. Inclusion criteria were a confirmed genetic diagnosis of one of 16 genetic conditions (Simons Searchlight) or absence of known pathogenic genetic findings in individuals with ASD (SPARK). Parent-reported age of acquisition of three gross motor and two expressive language milestones was described and categorized as on-time or delayed, relative to normative expectations.ResultsDevelopmental milestone profiles of probands with genetic conditions were marked by extensive delays (including nonattainment), with highest severity in single gene conditions and more delays than idiopathic ASD in motor skills. Compared with idiopathic ASD, the median odds of delay among the genetic groups were higher by 8.3 times (IQR 5.8-16.3) for sitting, 12.4 times (IQR 5.3-19.5) for crawling, 26.8 times (IQR 7.7-41.1) for walking, 2.7 times (IQR 1.7-5.5) for single words, and 5.7 times (IQR 2.7-18.3) for combined words.ConclusionsDelays in developmental milestones, particularly in gross motor skills, are frequent and may be among the earliest indicators of differentially affected developmental processes in specific genetically defined conditions associated with ASD, as compared with those with clinical diagnoses of idiopathic ASD. The possibility of different developmental pathways leading to ASD-associated phenotypes should be considered when deciding how to employ specific genetic conditions as models for ASD.

Lay abstractWhat is Already Known about This Subject: Genetics, (including de novo mutations), environmental factors (including toxic exposures), and their interactions impact autism spectrum disorder etiology. Paternal smoking is a candidate risk for autism spectrum disorder due to biological plausibility, high prevalence, and potential intervention.What This Study Adds: This original study and its replication confirms that paternal factors can substantially contribute to autism spectrum disorder risk for their offspring. It specifically indicates that paternal smoking both before and during pregnancy contributes significantly to autism spectrum disorder risk.Implications for practice, research, or policy: Smoking prevention, especially in pregnancy planning, may decrease autism spectrum disorder risk in offspring.

Loss-of-function variants in the gene SCN2A, which encodes the sodium channel NaV1.2, are strongly associated with autism spectrum disorder and intellectual disability. An estimated 20%-30% of children with these variants also suffer from epilepsy, with altered neuronal activity originating in neocortex, a region where NaV1.2 channels are expressed predominantly in excitatory pyramidal cells. This is paradoxical, as sodium channel loss in excitatory cells would be expected to dampen neocortical activity rather than promote seizure. Here, we examined pyramidal neurons lacking NaV1.2 channels and found that they were intrinsically hyperexcitable, firing high-frequency bursts of action potentials (APs) despite decrements in AP size and speed. Compartmental modeling and dynamic-clamp recordings revealed that NaV1.2 loss prevented potassium channels from properly repolarizing neurons between APs, increasing overall excitability by allowing neurons to reach threshold for subsequent APs more rapidly. This cell-intrinsic mechanism may, therefore, account for why SCN2A loss-of-function can paradoxically promote seizure.

Healthy pregnancy requires tolerance to fetal alloantigens as well as syngeneic embryonic and placental antigens. Given the importance of the autoimmune regulator (Aire) gene in self-tolerance, we investigated the role of Aire-expressing cells in maternal-fetal tolerance. We report that maternal ablation of Aire-expressing (Aire +) cells during early mouse pregnancy caused intrauterine growth restriction (IUGR) in both allogeneic and syngeneic pregnancies. This phenotype is immune mediated, as IUGR was rescued in Rag1-deficient mice, and involved a memory response, demonstrated by recurrence of severe IUGR in second pregnancies. Single-cell RNA sequencing demonstrated that Aire + cell depletion in pregnancy results in expansion of activated T cells, particularly T follicular helper cells. Unexpectedly, selective ablation of either Aire-expressing medullary thymic epithelial cells or extrathymic Aire-expressing cells (eTACs) mapped the IUGR phenotype exclusively to eTACs. Thus, we report a previously undescribed mechanism for the maintenance of maternal-fetal immune homeostasis and demonstrate that eTACs protect the conceptus from immune-mediated IUGR.

Interactions between genetic variants-epistasis-is pervasive in model systems and can profoundly impact evolutionary adaption, population disease dynamics, genetic mapping, and precision medicine efforts. In this work, we develop a model for structured polygenic epistasis, called coordinated epistasis (CE), and prove that several recent theories of genetic architecture fall under the formal umbrella of CE. Unlike standard epistasis models that assume epistasis and main effects are independent, CE captures systematic correlations between epistasis and main effects that result from pathway-level epistasis, on balance skewing the penetrance of genetic effects. To test for the existence of CE, we propose the even-odd (EO) test and prove it is calibrated in a range of realistic biological models. Applying the EO test in the UK Biobank, we find evidence of CE in 18 of 26 traits spanning disease, anthropometric, and blood categories. Finally, we extend the EO test to tissue-specific enrichment and identify several plausible tissue-trait pairs. Overall, CE is a dimension of genetic architecture that can capture structured, systemic forms of epistasis in complex human traits.

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.

3D imaging data necessitate 3D reference atlases for accurate quantitative interpretation. Existing computational methods to generate 3D atlases from 2D-derived atlases result in extensive artifacts, while manual curation approaches are labor-intensive. We present a computational approach for 3D atlas construction that substantially reduces artifacts by identifying anatomical boundaries in the underlying imaging data and using these to guide 3D transformation. Anatomical boundaries also allow extension of atlases to complete edge regions. Applying these methods to the eight developmental stages in the Allen Developing Mouse Brain Atlas (ADMBA) led to more comprehensive and accurate atlases. We generated imaging data from 15 whole mouse brains to validate atlas performance and observed qualitative and quantitative improvement (37% greater alignment between atlas and anatomical boundaries). We provide the pipeline as the MagellanMapper software and the eight 3D reconstructed ADMBA atlases. These resources facilitate whole-organ quantitative analysis between samples and across development.

MagellanMapper is a software suite designed for visual inspection and end-to-end automated processing of large-volume, 3D brain imaging datasets in a memory-efficient manner. The rapidly growing number of large-volume, high-resolution datasets necessitates visualization of raw data at both macro- and microscopic levels to assess the quality of data, as well as automated processing to quantify data in an unbiased manner for comparison across a large number of samples. To facilitate these analyses, MagellanMapper provides both a graphical user interface for manual inspection and a command-line interface for automated image processing. At the macroscopic level, the graphical interface allows researchers to view full volumetric images simultaneously in each dimension and to annotate anatomical label placements. At the microscopic level, researchers can inspect regions of interest at high resolution to build ground truth data of cellular locations such as nuclei positions. Using the command-line interface, researchers can automate cell detection across volumetric images, refine anatomical atlas labels to fit underlying histology, register these atlases to sample images, and perform statistical analyses by anatomical region. MagellanMapper leverages established open-source computer vision libraries and is itself open source and freely available for download and extension. © 2020 Wiley Periodicals LLC. Basic Protocol 1: MagellanMapper installation Alternate Protocol: Alternative methods for MagellanMapper installation Basic Protocol 2: Import image files into MagellanMapper Basic Protocol 3: Region of interest visualization and annotation Basic Protocol 4: Explore an atlas along all three dimensions and register to a sample brain Basic Protocol 5: Automated 3D anatomical atlas construction Basic Protocol 6: Whole-tissue cell detection and quantification by anatomical label Support Protocol: Import a tiled microscopy image in proprietary format into MagellanMapper.

Clinical exome sequencing is frequently used to identify gene-disrupting variants in individuals with neurodevelopmental disorders. While splice-disrupting variants are known to contribute to these disorders, clinical interpretation of cryptic splice variants outside of the canonical splice site has been challenging. Here, we discuss papers that improve such detection.

BackgroundThe cause of most fetal anomalies is not determined prenatally. Exome sequencing has transformed genetic diagnosis after birth, but its usefulness for prenatal diagnosis is still emerging. Nonimmune hydrops fetalis (NIHF), a fetal abnormality that is often lethal, has numerous genetic causes; the extent to which exome sequencing can aid in its diagnosis is unclear.MethodsWe evaluated a series of 127 consecutive unexplained cases of NIHF that were defined by the presence of fetal ascites, pleural or pericardial effusions, skin edema, cystic hygroma, increased nuchal translucency, or a combination of these conditions. The primary outcome was the diagnostic yield of exome sequencing for detecting genetic variants that were classified as either pathogenic or likely pathogenic according to the criteria of the American College of Medical Genetics and Genomics. Secondary outcomes were the percentage of cases associated with specific genetic disorders and the proportion of variants that were inherited.ResultsIn 37 of the 127 cases (29%), we identified diagnostic genetic variants, including those for disorders affecting the RAS-MAPK cell-signaling pathway (known as RASopathies) (30% of the genetic diagnoses); inborn errors of metabolism and musculoskeletal disorders (11% each); lymphatic, neurodevelopmental, cardiovascular, and hematologic disorders (8% each); and others. Prognoses ranged from a relatively mild outcome to death during the perinatal period. Overall, 68% of the cases (25 of 37) with diagnostic variants were autosomal dominant (of which 12% were inherited and 88% were de novo), 27% (10 of 37) were autosomal recessive (of which 95% were inherited and 5% were de novo), 1 was inherited X-linked recessive, and 1 was of uncertain inheritance. We identified potentially diagnostic variants in an additional 12 cases.ConclusionsIn this large case series of 127 fetuses with unexplained NIHF, we identified a diagnostic genetic variant in approximately one third of the cases. (Funded by the UCSF Center for Maternal-Fetal Precision Medicine and others; ClinicalTrials.gov number, NCT03412760.).

To discover regulatory elements driving the specificity of gene expression in different cell types and regions of the developing human brain, we generated an atlas of open chromatin from nine dissected regions of the mid-gestation human telencephalon, as well as microdissected upper and deep layers of the prefrontal cortex. We identified a subset of open chromatin regions (OCRs), termed predicted regulatory elements (pREs), that are likely to function as developmental brain enhancers. pREs showed temporal, regional, and laminar differences in chromatin accessibility and were correlated with gene expression differences across regions and gestational ages. We identified two functional de novo variants in a pRE for autism risk gene SLC6A1, and using CRISPRa, demonstrated that this pRE regulates SCL6A1. Additionally, mouse transgenic experiments validated enhancer activity for pREs proximal to FEZF2 and BCL11A. Thus, this atlas serves as a resource for decoding neurodevelopmental gene regulation in health and disease.

We identify a set of common phenotypic modifiers that interact with five independent autism gene orthologs (RIMS1, CHD8, CHD2, WDFY3, ASH1L) causing a common failure of presynaptic homeostatic plasticity (PHP) in Drosophila. Heterozygous null mutations in each autism gene are demonstrated to have normal baseline neurotransmission and PHP. However, PHP is sensitized and rendered prone to failure. A subsequent electrophysiology-based genetic screen identifies the first known heterozygous mutations that commonly genetically interact with multiple ASD gene orthologs, causing PHP to fail. Two phenotypic modifiers identified in the screen, PDPK1 and PPP2R5D, are characterized. Finally, transcriptomic, ultrastructural and electrophysiological analyses define one mechanism by which PHP fails; an unexpected, maladaptive up-regulation of CREG, a conserved, neuronally expressed, stress response gene and a novel repressor of PHP. Thus, we define a novel genetic landscape by which diverse, unrelated autism risk genes may converge to commonly affect the robustness of synaptic transmission.

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.

Abstract Gene expression differs between cell types and regions within complex tissues such as the developing brain. To discover regulatory elements underlying this specificity, we generated genome-wide maps of chromatin accessibility in eleven anatomically-defined regions of the developing human telencephalon, including upper and deep layers of the prefrontal cortex. We predicted a subset of open chromatin regions (18%) that are most likely to be active enhancers, many of which are dynamic with 26% differing between early and late mid-gestation and 28% present in only one brain region. These region-specific predicted regulatory elements (pREs) are enriched proximal to genes with expression differences across regions and developmental stages and harbor distinct sequence motifs that suggest potential upstream regulators of regional and temporal transcription. We leverage this atlas to identify regulators of genes associated with autism spectrum disorder (ASD) including an enhancer of BCL11A , validated in mouse, and two functional de novo mutations in individuals with ASD in an enhancer of SLC6A1 , validated in neuroblastoma cells. These applications demonstrate the utility of this atlas for decoding neurodevelopmental gene regulation in health and disease. Summary To discover regulatory elements driving the specificity of gene expression in different cell types and regions of the developing human brain, we generated an atlas of open chromatin from eleven dissected regions of the mid-gestation human telencephalon, including upper and deep layers of the prefrontal cortex. We identified a subset of open chromatin regions (OCRs), termed predicted regulatory elements (pREs), that are likely to function as developmental brain enhancers. pREs showed regional differences in chromatin accessibility, including many specific to one brain region, and were correlated with gene expression differences across the same regions and gestational ages. pREs allowed us to map neurodevelopmental disorder risk genes to developing telencephalic regions, and we identified three functional de novo noncoding variants in pREs that alter enhancer function. In addition, transgenic experiments in mouse validated enhancer activity for a pRE proximal to BCL11A , showing how this atlas serves as a resource for decoding neurodevelopmental gene regulation in health and disease.

De novo and inherited rare genetic disorders (RGDs) are a major cause of human morbidity, frequently involving neuropsychiatric symptoms. Recent advances in genomic technologies and data sharing have revolutionized the identification and diagnosis of RGDs, presenting an opportunity to elucidate the mechanisms underlying neuropsychiatric disorders by investigating the pathophysiology of high-penetrance genetic risk factors. Here we seek out the best path forward for achieving these goals. We think future research will require consistent approaches across multiple RGDs and developmental stages, involving both the characterization of shared neuropsychiatric dimensions in humans and the identification of neurobiological commonalities in model systems. A coordinated and concerted effort across patients, families, researchers, clinicians and institutions, including rapid and broad sharing of data, is now needed to translate these discoveries into urgently needed therapies.

Autism spectrum disorder (ASD) is strongly associated with de novo gene mutations. One of the most commonly affected genes is SCN2A. ASD-associated SCN2A mutations impair the encoded protein NaV1.2, a sodium channel important for action potential initiation and propagation in developing excitatory cortical neurons. The link between an axonal sodium channel and ASD, a disorder typically attributed to synaptic or transcriptional dysfunction, is unclear. Here we show that NaV1.2 is unexpectedly critical for dendritic excitability and synaptic function in mature pyramidal neurons in addition to regulating early developmental axonal excitability. NaV1.2 loss reduced action potential backpropagation into dendrites, impairing synaptic plasticity and synaptic strength, even when NaV1.2 expression was disrupted in a cell-autonomous fashion late in development. These results reveal a novel dendritic function for NaV1.2, providing insight into cellular mechanisms probably underlying circuit and behavioral dysfunction in ASD.

Autism spectrum disorder (ASD) is a common disorder that causes substantial distress. Heritability studies consistently show a strong genetic contribution, raising the hope that identifying ASD-associated genetic variants will offer insights into neurobiology and ultimately therapeutics. Next-generation sequencing (NGS) enabled the identification of disruptive variants throughout protein-coding regions of the genome. Alongside large cohorts and novel statistical methods, these NGS methods revolutionized ASD gene discovery. NGS methods have also contributed substantially to functional genetic data, such as gene expression, used to understand the neurobiological consequences of disrupting these ASD-associated genes. These functional data are also critical for annotating the noncoding genome as whole-genome sequencing (WGS) begins to provide initial insights outside of protein-coding regions. NGS methods still have a major role to play, as do similarly transformative advances in stem cell and gene-editing methods, in translating genetic discoveries into a first generation of ASD therapeutics.

Acetylation of the lysine residues in histones and other DNA-binding proteins plays a major role in regulation of eukaryotic gene expression. This process is controlled by histone acetyltransferases (HATs/KATs) found in multiprotein complexes that are recruited to chromatin by the scaffolding subunit transformation/transcription domain-associated protein (TRRAP). TRRAP is evolutionarily conserved and is among the top five genes intolerant to missense variation. Through an international collaboration, 17 distinct de novo or apparently de novo variants were identified in TRRAP in 24 individuals. A strong genotype-phenotype correlation was observed with two distinct clinical spectra. The first is a complex, multi-systemic syndrome associated with various malformations of the brain, heart, kidneys, and genitourinary system and characterized by a wide range of intellectual functioning; a number of affected individuals have intellectual disability (ID) and markedly impaired basic life functions. Individuals with this phenotype had missense variants clustering around the c.3127G>A p.(Ala1043Thr) variant identified in five individuals. The second spectrum manifested with autism spectrum disorder (ASD) and/or ID and epilepsy. Facial dysmorphism was seen in both groups and included upslanted palpebral fissures, epicanthus, telecanthus, a wide nasal bridge and ridge, a broad and smooth philtrum, and a thin upper lip. RNA sequencing analysis of skin fibroblasts derived from affected individuals skin fibroblasts showed significant changes in the expression of several genes implicated in neuronal function and ion transport. Thus, we describe here the clinical spectrum associated with TRRAP pathogenic missense variants, and we suggest a genotype-phenotype correlation useful for clinical evaluation of the pathogenicity of the variants.

Whole-genome sequencing (WGS) has facilitated the first genome-wide evaluations of the contribution of de novo noncoding mutations to complex disorders. Using WGS, we identified 255,106 de novo mutations among sample genomes from members of 1902 quartet families in which one child, but not a sibling or their parents, was affected by autism spectrum disorder (ASD). In contrast to coding mutations, no noncoding functional annotation category, analyzed in isolation, was significantly associated with ASD. Casting noncoding variation in the context of a de novo risk score across multiple annotation categories, however, did demonstrate association with mutations localized to promoter regions. We found that the strongest driver of this promoter signal emanates from evolutionarily conserved transcription factor binding sites distal to the transcription start site. These data suggest that de novo mutations in promoter regions, characterized by evolutionary and functional signatures, contribute to ASD.

To broaden our understanding of human neurodevelopment, we profiled transcriptomic and epigenomic landscapes across brain regions and/or cell types for the entire span of prenatal and postnatal development. Integrative analysis revealed temporal, regional, sex, and cell type-specific dynamics. We observed a global transcriptomic cup-shaped pattern, characterized by a late fetal transition associated with sharply decreased regional differences and changes in cellular composition and maturation, followed by a reversal in childhood-adolescence, and accompanied by epigenomic reorganizations. Analysis of gene coexpression modules revealed relationships with epigenomic regulation and neurodevelopmental processes. Genes with genetic associations to brain-based traits and neuropsychiatric disorders (including MEF2C, SATB2, SOX5, TCF4, and TSHZ3) converged in a small number of modules and distinct cell types, revealing insights into neurodevelopment and the genomic basis of neuropsychiatric risks.

Advances in gene discovery for neurodevelopmental disorders have identified SCN2A dysfunction as a leading cause of infantile seizures, autism spectrum disorder, and intellectual disability. SCN2A encodes the neuronal sodium channel NaV1.2. Functional assays demonstrate strong correlation between genotype and phenotype. This insight can help guide therapeutic decisions and raises the possibility that ligands that selectively enhance or diminish channel function may improve symptoms. The well-defined function of sodium channels makes SCN2A an important test case for investigating the neurobiology of neurodevelopmental disorders more generally. Here, we discuss the progress made, through the concerted efforts of a diverse group of academic and industry scientists as well as policy advocates, in understanding and treating SCN2A-related disorders.

In the version of this article initially published, the consortium authorship and corresponding authors were not presented correctly. In the PDF and print versions, the Whole Genome Sequencing for Psychiatric Disorders (WGSPD) consortium was missing from the author list at the beginning of the paper, where it should have appeared as the seventh author; it was present in the author list at the end of the paper, but the footnote directing readers to the Supplementary Note for a list of members was missing. In the HTML version, the consortium was listed as the last author instead of as the seventh, and the line directing readers to the Supplementary Note for a list of members appeared at the end of the paper under Author Information but not in association with the consortium name itself. Also, this line stated that both member names and affiliations could be found in the Supplementary Note; in fact, only names are given. In all versions of the paper, the corresponding author symbols were attached to A. Jeremy Willsey, Steven E. Hyman, Anjene M. Addington and Thomas Lehner; they should have been attached, respectively, to Steven E. Hyman, Anjene M. Addington, Thomas Lehner and Nelson B. Freimer. As a result of this shift, the respective contact links in the HTML version did not lead to the indicated individuals. The errors have been corrected in the HTML and PDF versions of the article.

Genomic association studies of common or rare protein-coding variation have established robust statistical approaches to account for multiple testing. Here we present a comparable framework to evaluate rare and de novo noncoding single-nucleotide variants, insertion/deletions, and all classes of structural variation from whole-genome sequencing (WGS). Integrating genomic annotations at the level of nucleotides, genes, and regulatory regions, we define 51,801 annotation categories. Analyses of 519 autism spectrum disorder families did not identify association with any categories after correction for 4,123 effective tests. Without appropriate correction, biologically plausible associations are observed in both cases and controls. Despite excluding previously identified gene-disrupting mutations, coding regions still exhibited the strongest associations. Thus, in autism, the contribution of de novo noncoding variation is probably modest in comparison to that of de novo coding variants. Robust results from future WGS studies will require large cohorts and comprehensive analytical strategies that consider the substantial multiple-testing burden.

RORα, the RAR-related orphan nuclear receptor alpha, is essential for cerebellar development. The spontaneous mutant mouse staggerer, with an ataxic gait caused by neurodegeneration of cerebellar Purkinje cells, was discovered two decades ago to result from homozygous intragenic Rora deletions. However, RORA mutations were hitherto undocumented in humans. Through a multi-centric collaboration, we identified three copy-number variant deletions (two de novo and one dominantly inherited in three generations), one de novo disrupting duplication, and nine de novo point mutations (three truncating, one canonical splice site, and five missense mutations) involving RORA in 16 individuals from 13 families with variable neurodevelopmental delay and intellectual disability (ID)-associated autistic features, cerebellar ataxia, and epilepsy. Consistent with the human and mouse data, disruption of the D. rerio ortholog, roraa, causes significant reduction in the size of the developing cerebellum. Systematic in vivo complementation studies showed that, whereas wild-type human RORA mRNA could complement the cerebellar pathology, missense variants had two distinct pathogenic mechanisms of either haploinsufficiency or a dominant toxic effect according to their localization in the ligand-binding or DNA-binding domains, respectively. This dichotomous direction of effect is likely relevant to the phenotype in humans: individuals with loss-of-function variants leading to haploinsufficiency show ID with autistic features, while individuals with de novo dominant toxic variants present with ID, ataxia, and cerebellar atrophy. Our combined genetic and functional data highlight the complex mutational landscape at the human RORA locus and suggest that dual mutational effects likely determine phenotypic outcome.

As technology advances, whole genome sequencing (WGS) is likely to supersede other genotyping technologies. The rate of this change depends on its relative cost and utility. Variants identified uniquely through WGS may reveal novel biological pathways underlying complex disorders and provide high-resolution insight into when, where, and in which cell type these pathways are affected. Alternatively, cheaper and less computationally intensive approaches may yield equivalent insights. Understanding the role of rare variants in the noncoding gene-regulating genome, through pilot WGS projects, will be critical to determine which of these two extremes best represents reality. With large cohorts, well-defined risk loci, and a compelling need to understand the underlying biology, psychiatric disorders have a role to play in this preliminary WGS assessment. The WGSPD consortium will integrate data for 18,000 individuals with psychiatric disorders, beginning with autism spectrum disorder, schizophrenia, bipolar disorder, and major depressive disorder, along with over 150,000 controls.

Congenital heart disease (CHD) is the leading cause of mortality from birth defects. Here, exome sequencing of a single cohort of 2,871 CHD probands, including 2,645 parent-offspring trios, implicated rare inherited mutations in 1.8%, including a recessive founder mutation in GDF1 accounting for ∼5% of severe CHD in Ashkenazim, recessive genotypes in MYH6 accounting for ∼11% of Shone complex, and dominant FLT4 mutations accounting for 2.3% of Tetralogy of Fallot. De novo mutations (DNMs) accounted for 8% of cases, including ∼3% of isolated CHD patients and ∼28% with both neurodevelopmental and extra-cardiac congenital anomalies. Seven genes surpassed thresholds for genome-wide significance, and 12 genes not previously implicated in CHD had >70% probability of being disease related. DNMs in ∼440 genes were inferred to contribute to CHD. Striking overlap between genes with damaging DNMs in probands with CHD and autism was also found.

BackgroundVariants in the SCN2A gene that disrupt the encoded neuronal sodium channel NaV1.2 are important risk factors for autism spectrum disorder (ASD), developmental delay, and infantile seizures. Variants observed in infantile seizures are predominantly missense, leading to a gain of function and increased neuronal excitability. How variants associated with ASD affect NaV1.2 function and neuronal excitability are unclear.MethodsWe examined the properties of 11 ASD-associated SCN2A variants in heterologous expression systems using whole-cell voltage-clamp electrophysiology and immunohistochemistry. Resultant data were incorporated into computational models of developing and mature cortical pyramidal cells that express NaV1.2.ResultsIn contrast to gain of function variants that contribute to seizure, we found that all ASD-associated variants dampened or eliminated channel function. Incorporating these electrophysiological results into a compartmental model of developing excitatory neurons demonstrated that all ASD variants, regardless of their mechanism of action, resulted in deficits in neuronal excitability. Corresponding analysis of mature neurons predicted minimal change in neuronal excitability.ConclusionsThis functional characterization thus identifies SCN2A mutation and NaV1.2 dysfunction as the most frequently observed ASD risk factor detectable by exome sequencing and suggests that associated changes in neuronal excitability, particularly in developing neurons, may contribute to ASD etiology.

Autism spectrum disorder (ASD) risk is influenced by common polygenic and de novo variation. We aimed to clarify the influence of polygenic risk for ASD and to identify subgroups of ASD cases, including those with strongly acting de novo variants, in which polygenic risk is relevant. Using a novel approach called the polygenic transmission disequilibrium test and data from 6,454 families with a child with ASD, we show that polygenic risk for ASD, schizophrenia, and greater educational attainment is over-transmitted to children with ASD. These findings hold independent of proband IQ. We find that polygenic variation contributes additively to risk in ASD cases who carry a strongly acting de novo variant. Lastly, we show that elements of polygenic risk are independent and differ in their relationship with phenotype. These results confirm that the genetic influences on ASD are additive and suggest that they create risk through at least partially distinct etiologic pathways.