Skip to main content

Coding and noncoding variants in EBF3 are involved in HADDS and simplex autism

Abstract

Background

Previous research in autism and other neurodevelopmental disorders (NDDs) has indicated an important contribution of protein-coding (coding) de novo variants (DNVs) within specific genes. The role of de novo noncoding variation has been observable as a general increase in genetic burden but has yet to be resolved to individual functional elements. In this study, we assessed whole-genome sequencing data in 2671 families with autism (discovery cohort of 516 families, replication cohort of 2155 families). We focused on DNVs in enhancers with characterized in vivo activity in the brain and identified an excess of DNVs in an enhancer named hs737.

Results

We adapted the fitDNM statistical model to work in noncoding regions and tested enhancers for excess of DNVs in families with autism. We found only one enhancer (hs737) with nominal significance in the discovery (p = 0.0172), replication (p = 2.5 × 10−3), and combined dataset (p = 1.1 × 10−4). Each individual with a DNV in hs737 had shared phenotypes including being male, intact cognitive function, and hypotonia or motor delay. Our in vitro assessment of the DNVs showed they all reduce enhancer activity in a neuronal cell line. By epigenomic analyses, we found that hs737 is brain-specific and targets the transcription factor gene EBF3 in human fetal brain. EBF3 is genome-wide significant for coding DNVs in NDDs (missense p = 8.12 × 10−35, loss-of-function p = 2.26 × 10−13) and is widely expressed in the body. Through characterization of promoters bound by EBF3 in neuronal cells, we saw enrichment for binding to NDD genes (p = 7.43 × 10−6, OR = 1.87) involved in gene regulation. Individuals with coding DNVs have greater phenotypic severity (hypotonia, ataxia, and delayed development syndrome [HADDS]) in comparison to individuals with noncoding DNVs that have autism and hypotonia.

Conclusions

In this study, we identify DNVs in the hs737 enhancer in individuals with autism. Through multiple approaches, we find hs737 targets the gene EBF3 that is genome-wide significant in NDDs. By assessment of noncoding variation and the genes they affect, we are beginning to understand their impact on gene regulatory networks in NDDs.

Background

Large-scale whole-genome sequencing (WGS) is becoming instrumental in assessing the contribution of protein-coding (coding), but more importantly noncoding variants in complex diseases [1]. Unlike coding exons, the boundaries of noncoding regions are not well defined and hence different types of annotations including but not limited to evolutionary conservation [2], sequence constraint [3], and epigenetic marks [4] are useful guides. Though genome-wide association studies (GWAS) have identified multiple common noncoding variants associated with human disorders [5,6,7], WGS has now provided access to rare de novo variants (DNVs) which are difficult to associate with phenotype without using aggregation methods [8,9,10,11]. These aggregation methods have been successfully used for rare coding variants [12,13,14,15,16], but have been challenging for noncoding regions because of the lack of clearly defined, discrete genomic boundaries and sequence-based models of variant effect.

Autism is a complex neurodevelopmental disorder with a heritability of ~ 80% [17]. Large copy number variants [18,19,20,21] and coding DNVs contribute to ~ 30% of cases with higher enrichment in females with autism and those with intellectual disability [12,13,14,15,16, 22, 23]. Recently, we and others have identified an overall enrichment of de novo [8,9,10,11, 24] or paternally inherited variants [25] within the regulatory sequence of individuals with autism. However, these studies have mostly assessed this aggregation of genetic burden across a large panel of pooled regulatory elements. To begin to parse out the underlying biology of autism DNVs in individual regulatory regions, we turned to VISTA, which is a database of functionally characterized developmental enhancers [26,27,28]. These enhancers were identified based on multiple strategies including sequence conservation and epigenetic signatures. Each enhancer has been tested in transgenic mouse assays providing information on the spatial-temporal dynamics of their activity during mammalian development. We adapted the fitDNM model [29] that was previously used to test for excess DNV load in coding regions so that it would work in noncoding regions. We then applied our updated version of the fitDNM model to VISTA enhancers with known ability to drive expression in the embryonic brain. Application of this test in 2671 families with autism (n = 9831 individuals) revealed one VISTA enhancer (named hs737) with nominal significance for excess of DNVs in autism in our discovery cohort (516 families), replication cohort (2155 families), and the combined dataset. We extensively tested enhancer hs737 in follow-up genomic, epigenomic, phenotypic, in silico, and in vitro analyses. Our analysis revealed this enhancer targets the transcription factor gene EBF3 which is enriched for coding DNVs in the hypotonia, ataxia, and delayed development syndrome (HADDS). We also identified new patients with coding variants in EBF3 and performed extensive phenotypic analysis. Combining this with phenotype data we collected for individuals with the hs737 enhancer DNVs, we found marked increases in phenotypic severity in individuals with coding than noncoding variants. This work provides critical insights into coding and noncoding DNVs at EBF3 and more generally in neurodevelopmental disorders.

Results

Statistical assessment of DNVs in VISTA elements

To assess for DNVs in individuals with autism, we aggregated DNV data from a WGS study of 2671 families with autism [30]. To test the enrichment of DNVs in noncoding regions, we modified the existing fitDNM [29]. In particular, we focused on 544 VISTA human noncoding enhancers (Supplemental Table S1, Supplemental Table S2) previously shown to have enhancer activity in the brain using a lacZ transgenic assay at embryonic day 11.5 in mice [27]. We assessed these same enhancers in our previous paper (Turner et al. [8]) studying 516 families (discovery cohort) and we wanted to test whether any of these enhancers replicated in a new set of 2155 families (replication cohort). Of these enhancers, we identified one (hs737, see Table 1) reaching nominal significance in both cohorts (discovery p = 0.0172, replication p = 2.5 × 10−3) and in the combined dataset (p = 1.1 × 10−4). The hs737 enhancer drives expression in the midbrain and hindbrain at embryonic day 11.5 (E11.5) [8].

Table 1 VISTA enhancers with an excess of de novo mutation based on fitDNM analysis

In each individual with a DNV in hs737, we also assessed the rest of the genome for other potentially relevant genomic variants that could be contributing to autism in each of the individuals (Supplemental Table S3, S4). Individual 13396.p1 had no other large de novo CNVs or coding de novo SNVs/indels. Individual 12975.p1 had three total de novo missense variants with one in each of the following genes CHD6, FAM129B, and KCNC1 and also had a 1.6-Mbp de novo deletion at 11q24.1 containing the following genes BLID, BSX, C11orf63, CRTAM, SORL1, and UBASH3B. Individual 11257.p1 had two de novo missense variants in each of the following genes RUVBL1 and VKORC1L1. We scored each of the variants using a clinical variant scoring program (https://franklin.genoox.com/clinical-db/home) and all of the variants were classified as variants of uncertain significance. We checked each individual’s polygenic risk score [31] for autism spectrum disorder, schizophrenia, and educational attainment and find no significant contribution for any of the three individuals (Supplementary Figure 1).

In vitro assessment of hs737 DNVs

In order to quantify the in vitro transcriptional effects of hs737 DNVs, we transfected the neuronal cell line Neuro2a with a reporter construct that had both the non-risk and risk allele of the enhancer individually cloned upstream of a luciferase gene and a minimal promoter [32]. As a control, we examined the expression of the known RET+3 enhancer, shown to be functionally active in this cell line [32]. Our data shows that all DNVs in hs737 led to a significant reduction in reporter gene expression when compared to their respective non-risk allele and the promoter-only construct and the control enhancer (RET+3) had high transcriptional activity (Fig. 1C). Thus, these DNVs in hs737 can individually affect the transcription of their cognate gene.

Fig. 1
figure 1

Characterization of DNVs in hs737. A Pedigrees of families with de novo variants in hs737. Lightning symbols indicate de novo variants with red = regulatory, purple = missense, and blue = deletion. Family identifiers are shown above the pedigree and the full-scale IQ is shown below each proband. B Sequence analysis of each of the three hs737 de novo mutations, identified in individuals with autism, including transcription factor binding site analysis results. C Results of luciferase assays in neuroblastoma (Neuro2a) cell lines with rs2435357 (RET+3) as a positive control for enhancer activity, promoter only (Basal), the wild type sequence of hs737 (hs737wt), and each of the three DNVs identified in individuals with autism. Error bars represent standard error (SE). D log2 normalized expression of genes from the transcription factor binding site analysis in the brain throughout development and adulthood. E Correlogram of candidate genes and EBF3 after performing regression, with positive control MECP2 and negative control CFTR

In silico transcription factor binding assessment

To better understand functional consequences of the hs737 variants we observed, we analyzed the reference and variant enhancer sequences using QBiC [33, 34] (http://qbic.gcb.duke.edu), a program designed to predict the impact of noncoding mutations on transcription factor binding sites. We found that the DNVs identified in individuals with autism were each predicted to impact transcription factor binding (Fig. 1B). Here we report the most significant hits, realizing that it could be any of the transcription factors in the family causing the functional impact on expression dynamics. All three variants are predicted to impact transcription factor binding via a transition mutation at a highly preferential base that is also highly conserved to frog (Fig. 1B). Each mutation occurs at a location within the position weight matrix that is predicted to completely change the binding status of the transcription factors (Fig. 1B). The first two variants are predicted to each respectively cause SOX30 and ARX to go from the bound state to unbound (SOX30 p = 2.48 × 10−189, z-score = −29.35; ARX p < 2.2 × 10−16, z-score = −40.27) while the third variant is predicted to create a new binding site for ARID5 (p < 2.2 × 10−16, z-score = 31.87). Mouse RNAseq at day E11.5 (see the “Methods” section) provides further support for these transcription factors being impactful as ARX and 12 members of the SOX family are in the 90th percentile and the ARID family members are in the 80th percentile of all genes that are expressed in the brain (Supplemental Table S6). We specifically assessed RNAseq data at E11.5 since that is the timepoint at which lacZ reporter expression was tested and observed in the VISTA enhancer database for the hs737 enhancer (http://enhancer.lbl.gov/cgi-bin/imagedb3.pl?form=presentation&show=1&experiment_id=737&organism_id=1).

Next, we turned to human RNAseq data from Brainspan (http://brainspan.org/rnaseq/search/index.html), which contains expression profiles from various timepoints ranging from early in development (8 post-conception weeks [pcw]) to adulthood (40 years) and examines many regions of the brain. We wanted to determine what transcription factor may be binding to the enhancer. As we describe below, the target gene of hs737 is EBF3. We hypothesized that the transcription factor that is binding may be correlated in expression with EBF3. We find that EBF3 is expressed throughout the brain prenatally and is most highly expressed through the hindbrain in the cerebellum and cerebellar cortex (Fig. 1D). We find further evidence for the transcription factors implicated by the in silico analysis. EBF3 expression is significantly correlated with SOX10 and the ARID1 family after performing a linear regression for age and brain region (Fig. 1E).

Dosage sensitivity of hs737 in the human genome

To further assess the impact of enhancer hs737 on NDDs, we measured the effect of its dosage on the NDD phenotypes. We hypothesize that if heterozygous point mutations in hs737 alter phenotype strongly, the enhancer would be dosage sensitive and show copy number variation (CNV) only in individuals with an NDD. To test this, we first applied two approaches to assess CNVs in hs737. We assessed the morbidity map [35, 36] database containing 29,085 individuals with neurodevelopmental disorders (NDDs) and 19,584 controls. In particular, we looked at the window analysis in Coe et al. [35] across the genome to identify the window containing hs737. This enhancer resided in a genomic window with an excess of deletions (case counts = 27, control counts = 0, p = 9.14 × 10−7) and duplications (case counts = 6, control counts = 0, p = 4.55 × 10−2) in individuals with NDDs. None of the 19,584 control individuals contained a CNV in this enhancer (Fig. 2). On average, these CNVs containing hs737 were 10,570,885 ± 5,789,023 bp long and overlapped 70 genes (Supplemental Table S7). The smallest CNV was 161 kbp and was in completely noncoding space.

Fig. 2
figure 2

Copy number variation over hs737. Displays the counts for both deletions and duplications over hs737 in individuals with neurodevelopmental disorders and controls

We also applied a tool to determine copy number (with paralog-specific sensitivity) in 1-kbp windows across the whole genome [37] to our WGS data and identified only one CNV in this enhancer. It was a deletion and occurred in proband 14091.p1 and upon further inspection was found to be part of a larger known deletion (hg38: chr10:126450330–133655780) [38]. To determine the frequency of deletion/duplication in this enhancer in individuals without autism, we also ran this copy number approach on the newly generated 3202 individuals (high-coverage WGS) from the 1000 genomes project [39]. Combining our WGS parental data, 1000 genomes project data, and morbidity map, there are no deletions or duplications in this enhancer in 28,128 non-NDD individuals (56,256 alleles). We also surveyed the gnomAD [40] database (v2.1) and observed no CNVs in the 10,847 individuals contained there (n = 21,694 alleles). To avoid possible double counting between gnomAD and other datasets, we do not present the aggregate data. Taken together, these results suggest that in addition to point mutations, CNVs involving hs737 may play a significant role in NDDs.

Epigenetic characteristics of hs737

Enhancers have well-characterized epigenetic signatures that are predictive of their activity in specific biological contexts. Thus, to examine the activity of hs737 in its native genomic context, we took advantage of available epigenomic datasets from relevant human samples. We found that hs737 has several hallmarks of neuronal enhancer activity in humans including H3K27ac enrichment in fetal brain tissue [41], DNaseI hypersensitive sites (DHS) indicative of chromatin accessibility in CNS tissues [8, 24], and conserved transcription factor binding sites (TFBS) [8] (HMR conserved transcription factor binding sites track in the UCSC Genome Browser [42]) (Supplemental Figure S2).

Hs737 was assessed for reporter activity in mouse embryos at E11.5 in the VISTA database. To further examine the activity of this element during development, we analyzed the orthologous region in the mouse genome (mouse ortholog of hs737, or mo-hs737) using epigenomic data from a mouse embryonic developmental time series recently published [43,44,45]. Consistent with the reporter expression pattern of hs737, we found that mo-hs737 has the chromatin signature of an active enhancer in the midbrain and hindbrain at E11.5, but not in the non-neuronal tissues assayed (Fig. 3A).

Fig. 3
figure 3

hs737 is a prenatal, brain-specific enhancer. A Genome browser view (chr7:136,079,964–136,087,591; mm10) of chromatin states in mouse from [43] called by chromHMM [44] based on eight histone modifications: H3K4me1, H3K4me2, H3K4me3, H3K27ac, H3K27me3, H3K36me3, H3K9me3, and H3K9ac. B Genome browser view (chr7:136,079,964–136,087,591; mm10) of ATAC-seq and H3K27ac ChIP-seq signal in midbrain, hindbrain, and forebrain at multiple developmental mouse stages from E11.5 to the day of birth (P0)

Strikingly, we found that characteristics of enhancer activity at mo-hs737 such as H3K27ac and chromatin accessibility (as measured by ATAC-seq) reach their height in brain tissue in mid to late gestation and decline at birth (Fig. 3B). This suggests that hs737/mo-hs737 may exert its regulatory influence specifically during embryonic development, which could explain its involvement in NDDs, and may point to developmental stages and model systems that are most appropriate for future studies of this element.

Gene target of hs737

To determine the potential gene target of the hs737 enhancer, we first looked at all of the genes residing within the same topologically associating domain (TAD) (hg38, chr10:128151746–130191746) [46]. We focused on these genes since they are the most likely to be the targets of this enhancer. The genes included C10orf143, CTAGE7P, EBF3, GLRX3, LINC01163, LOC728327, MGMT, and MIR4297; the only gene in this region that is constrained in the human population with a gnomAD [3] o/e value for loss-of-function SNVs of 0.03 is EBF3. To examine physical interactions between hs737/mo-hs737 and genes in this region, we analyzed two separate high-resolution Hi-C datasets examining the 3D architecture of the genome throughout mouse neuronal differentiation and human fetal corticogenesis. The mouse dataset was generated in a differentiation course from mouse embryonic stem cells (ESCs) to neural progenitor cells (NPCs) and then to cortical-like neurons (CNs) [47] with Hi-C at each of these major cell transition stages. We found that the Ebf3 promoter makes strong contacts across a large region encompassing mouse hs737 (mo-hs737) and that these interactions become stronger during differentiation to NPCs and CNs (Fig. 4B). We used HiCCUPS [48] to call loops at each stage and found that in CNs there is a ~1.3-Mbp loop that brings mo-hs737 into close proximity with the Ebf3 promoter (loop anchors chr7:136,050,000–136,075,000 and chr7:137,300,000–137,325,000). We did not observe loops between mo-hs737 and any other genes on the chromosome.

Fig. 4
figure 4

EBF3 is the gene target of hs737. A Schematic of hs737 and target genes. Gray boxes represent promoters and colored boxes represent gene bodies and red box represents hs737. Hi-C contact map generated using data from Won et al. [49] visualized with Juicebox [48] at 25 kbp. Heatmaps are symmetrical across the diagonal, except that HiCCUPS loop calls are shown as black boxes in the upper right half of each heatmap. B Hi-C contact maps from Bonev et al. [47] visualized with Juicebox [48] at 5-kbp resolution. C EBF3 protein diagram (plotted using the DOG protein plotter [67]) with DNVs identified in NDDs. Shown in blue are the missense variants and in red are the loss-of-function variants. D 3D model (plotted using the MuPIT program [68, 69]) of the EBF3 protein with DNVs identified in individuals with NDDs shown in green. E Genes with promoters bound by EBF3 based on ChIP sequencing in SK-N-SH cells. Enrichment is seen for the promoters of known NDD genesets

Next, analyzing the human corticogenesis data which was generated by dissecting three fetal brain samples into cortical plate (CP) and germinal zone (GZ) layers [49], we created chromatin contact maps using Juicer [48] for each dissection layer and called loops again using HiCCUPS [48]. We find that the EBF3 promoter interacts with the region of hs737 in both the GZ (loop anchors chr10:128,535,000–128,540,000) and CP (loop anchors chr10:128,525,001–128,550,00 and chr10:129,995,001–129,975,000) samples. We note that there is a loop in CP between hs737 and a second gene GRLX3 (loop anchors ch10:128,525,001–128,550,000 and chr10:130,125,001–130,150,000), a gene known to be involved in multiple mitochondrial dysfunctions syndrome (Fig. 4A). However, the loop between hs737 and GLRX3 is both weaker and less significant than the one formed between hs737 and EBF3.

Thus, narrowing in on EBF3, we searched recent literature on coding DNVs and found EBF3 is the only gene in this TAD that has known statistical enrichment for coding DNVs in NDDs [50, 51]. Combining data from these two previous studies (Fig. 4C, D; Supplemental Table S8), we also saw genome-wide significance (chimpanzee-human [50] missense p = 8.12 × 10−35, chimpanzee-human loss-of-function p = 2.26 × 10−13, denovolyzeR [52] missense p = 4.79 × 10−13, denovolyzeR loss-of-function p = 7.97 × 10−22) for coding DNVs using two different statistical tests. This is also the gene for the Mendelian phenotype hypotonia, ataxia, and delayed development syndrome (HADDS [53]) and has been characterized in detail in a set of ten patients [54]. Searching the GTEx database, we find that EBF3 is widely expressed throughout the body and in the brain (Fig. 6A) and similarly find in the Human Protein Atlas that EBF3 is detected in many tissues including the brain (Fig. 6B), whereas hs737 displays much more restricted activity only being active throughout the midbrain and hindbrain (Fig. 6C).

EBF3 regulates many NDD-significant genes

To assess the global transcriptional control of EBF3, we analyzed its genome-wide bonding profile from chromatin immunoprecipitation sequencing (ChIP-seq) in the human neuroblastoma cell line SK-N-SH [55]. We mapped all EBF3 peaks to promoters to detect EBF3 binding and identified 3100 genes (16% of all genes in the genome) bound by EBF3. We then focused in on genes with statistical enrichment for coding DNVs in NDDs and that were bound by EBF3 (Supplemental Table S9). Of the 253 significant NDD genes in Coe et al. [50], 26.9% of them were bound by EBF3 at their promoter (p = 8.95 × 10−6, OR = 1.93). Of the 285 significant genes in another study [51], 26.3% of them were bound by EBF3 (p = 7.43 × 10−6, OR = 1.87) (Fig. 4E). Many of these bound NDD genes (Supplemental Table S8,9) are involved in chromatin regulation (Chromatin Binding Gene Ontology p = 3.2 × 10−7) and/or transcription factor activity (DNA Binding Gene Ontology p = 6.4 × 10−11) (e.g., CHD8, CHD2, ARID1B) indicating that EBF3 may be a master-regulator of many NDD genes. Chromatin binding genes account for a sizeable fraction of DNM attributable cases of autism [12, 14, 56], suggesting that EBF3 disruption could result in a milder phenotype in the spectrum, as we observe in these cases. We find further support for this by observing that EBF3 expression is highly correlated (r > |0.6|) with 116 high confidence SFARI genes (SFARI score < 3) (Fig. 5A). Within the cluster of genes positively correlated with EBF3, there is a significant enrichment of genes involved in nucleosome organization (FDR p = 3.95 × 10−2), regulation of histone modification (FDR p = 2.12 × 10−2), chromatin remodeling (FDR p = 2.23 × 10−2), chromosome organization (FDR p = 3.00 × 10−7), chromatin organization (FDR p = 6.58 × 10−7), and positive regulation of the cell cycle process (FDR p = 3.68 × 10−2). Performing a network analysis using String-db, we find that the genes positively correlated with EBF3 have significantly more interactions than expected (PPI enrichment p-value < 1.0 × 10−16, expected edges 95, observed edges 202) (Supplemental Figure 3). However, EBF3 is only connected in this network using a low confidence threshold for connectivity, so we also performed a second network analysis using GeneMania (https://genemania.org) and find EBF3 is part of the network with evidence from genetic interactions (Fig. 5B).

Fig. 5
figure 5

EBF3 gene network analysis. A Correlation matrix of EBF3 with other high scoring SFARI genes (score < 3) after performing regression and hierarchical clustering and having an absolute correlation greater than 0.6. There is a significant enrichment of genes involved in chromatin binding. B Network analysis of genes from cluster 2 from Genemania where each gene is a node and different forms of supporting evidence for an interaction are edges

Phenotypes of individuals with DNV in hs737

In order to understand more about the phenotypic consequences of variation in the enhancer, we reviewed de-identified phenotype information for each of the individuals with autism that had a DNV in hs737. As gene discovery in simplex autism, based on large CNVs and coding DNVs, has yielded the most findings in females and in individuals who have intellectual disability (full-scale IQ < 70) [12], the first two phenotypes that we assessed were sex and the full-scale IQ scores. All three individuals were male and their full-scale IQs were 103, 89, and 132, respectively (Fig. 1A), suggesting that none of the individuals had an intellectual disability. We also found that all three individuals had evidence of motor problems and/or hypotonia.

Comparing the individuals with noncoding variants to 14 previously published individuals (13 probands) [54, 55] and 7 new individuals with coding DNVs in EBF3, we find that those with coding mutations in EBF3 typically have more severe phenotypic consequences (Fig. 6D). All individuals with a coding DNV in EBF3 regardless of position within the protein presented with an intellectual disability or global developmental delay, while no individuals with a DNV in hs737 had either of these phenotypes. Individuals with DNVs in EBF3 had higher rates of ataxia compared to all individuals with DNVs in EBF3, while this phenotype was also absent in individuals with hs737 DNVs. Diagnoses of cerebellar vermis hypoplasia were also found in 4/10 individuals with EBF3 DNVs who had a brain MRI. Within the EBF3 mutation group, there are significantly higher rates of intellectual disability or global developmental delay compared to the hs737 group (p = 0.00088) and we find significantly higher rates of autism within the hs737 DNV group when compared to the EBF3 group (p = 0.0088) (Fig. 6D) (Supplemental Table S10).

Fig. 6
figure 6

Consequence of coding and noncoding variation in EBF3. A GTEx expression data of EBF3 with color corresponding to the organ system. B Human Protein Atlas highlighting where EBF3 expression is detected in the human body. C LacZ staining assay for reporter activity driven by the hs737 enhancer at mouse E11.5. D Phenotypic analysis comparing the frequency of ataxia, hypotonia, ID and GDD, autism, and having 7 or more symptoms between all patients, individuals with EBF3 mutations, individuals with mutations specifically in the EBF3 DNA binding domain, and in the hs737 enhancer. E Gene regulatory network encompassing EBF3 built using current molecular biological knowledge

Discussion

Eleven years ago, the Simons Simplex Collection began [57] and started its efforts to understand the role of genomic variation in simplex autism. It was hypothesized that there is a contribution from DNVs in these simplex autism families [57]. Over the 11 years, microarray [18, 20, 58], WES [12,13,14, 16], and now WGS data [8,9,10, 24, 25] have been generated to fuel this discovery. The first fruits of these efforts were large copy number variants and coding DNVs. They have turned out to be critical for explaining ~ 30% of individuals with autism [12]. Intriguingly, these variants have been found to be enriched more in females and/or individuals who also had intellectual disability [12]. In particular from this work, coupled with the study of DNVs in neurodevelopmental disorders more broadly, there are now > 100 genes with genome-wide significance for excess of coding variation [50, 51, 59] including the gene EBF3. Recent ongoing efforts looking at common variation are providing insights into other aspects of the genetics of autism and explain ~ 50% of autism risk [60]. In our study, we focused on the elusive noncoding DNVs for which we and others have seen aggregate evidence for enrichment in promoters and enhancers [8,9,10]. We assessed VISTA enhancers in a discovery cohort of 516 families previously published [8] and a replication cohort of 2155 new families. Recent work [11] has indicated the importance of these enhancer regions and we find one reaching nominal significance in the discovery, replication, and combined cohorts (hs737).

We found a hs737 DNV in three different individuals with autism (one in the discovery cohort, two in the replication cohort). The individuals with DNVs had shared phenotypes including being male, intact cognitive function, and all had hypotonia or motor delay (unlike coding DNVs which are enriched for females with intellectual disability [12]). Previous work examining individuals with mutations in the same gene show shared phenotypes [61, 62] much like our individuals with enhancer DNVs. Each DNV had a quantitative effect on reporter gene expression in our in vitro assessment. These DNVs are at highly conserved nucleotides and are predicted to affect binding of transcription factors at the enhancer. Beyond de novo single-nucleotide variants, CNVs encompassing hs737 are also enriched in individuals with NDDs. In our assessment of 28,128 non-NDD individuals, there are none with a deletion or duplication of this enhancer suggesting that it is dosage sensitive in the human population. These lines of genetic evidence support the finding that this enhancer has an important role in the human genome.

Analysis of epigenetic data shows that this enhancer is active in the embryonic brain. A major hurdle in the study of enhancers is determining which gene they regulate. This is especially relevant when enhancers are very distal from the promoter they target. Since hs737 resides in a large noncoding region, we utilized an innovative approach of combining constraint information for genes within the TAD, enrichment of coding DNVs in genes in the TAD, and chromatin contact data. This analysis led to our identification of EBF3 as the gene targeted by hs737, in both embryonic mouse and human fetal brain, with an interaction across a distance of ~ 1.4 Mbp in the human fetal brain. EBF3 is a well-established NDD gene with genome-wide significant enrichment of coding DNVs and an established syndrome called HADDS that shares the phenotype of hypotonia with the individuals with autism in our study that had hs737 noncoding DNVs. HADDS is a severe phenotype and affects many parts of the body. This is likely because EBF3 is a transcription factor that is ubiquitously expressed in humans [63] (Fig. 6A, B). It affects many genes in the genome and in particular, in neuronal cells, is enriched for regulating other NDD genes involved in further regulation. We speculate that since hs737 is a brain-specific enhancer that is why we see a less severe phenotype in the individuals with autism in our study than in HADDS (Fig. 6C–E). Aggregating information from the literature and from this current study, we can begin to build the gene regulatory network containing EBF3 (Fig. 6E) and in this network we can see the importance of genotypic and phenotypic assessment of individuals. In particular, the higher the location of the variant in the network, the more severe the phenotypic consequence (e.g., mutation in ARX), and the lower in the network, the less severe the phenotypic consequence (e.g., mutation in the hs737 enhancer). Along with this observation, we point out the critical work of completing the allelic series for each of the genome-wide significant NDD genes. For EBF3, we show that coding and noncoding DNVs result in more and less severe phenotypic outcomes. This type of gene regulatory network building that incorporates coding and noncoding variation is essential for understanding the etiology of autism. For continued delineation of this gene regulatory network, it will be critical to move beyond the in vitro work to in vivo models as in vitro models can be a limitation when studying noncoding elements with specific spatiotemporal activity.

Conclusions

We identify hs737 as an enhancer with excess DNVs in autism and find through several approaches that it is brain-specific and targets the gene EBF3. This study provides critical insights into noncoding DNVs in autism and how they can have similar and differential effects on phenotypic outcomes. This work provides a framework for both future studies of noncoding DNVs and considerations of effects at the level of gene regulatory networks.

Methods

DNVs in 2671 autism families

We accessed DNV data from Wilfert et al. [30] through SFARI Base (accession: SFARI_SSC_WGS_2a, https://base.sfari.org/).

Statistical assessment of DNVs

A list of VISTA enhancers driving the expression of their target genes in the brain was downloaded from the VISTA enhancer browser [27]. DNVs were annotated to each enhancer using bedtools [64]. We modified the fitDNM statistical approach (https://github.com/TNTurnerLab/fitDNM) [29], a method to assess the excess mutational load of DNVs using variant-specific mutation rates calculated based on local sequence context, now applied to noncoding variants in the VISTA brain enhancers.

Copy number assessment of hs737

To test copy number variant enrichment in morbidity map [35], we downloaded Supplementary Dataset 1 from Coe et al. [35] and identified the window in the genome containing hs737. We report in this paper the case counts, control counts, and p-values for deletions and duplications in this window.

The QuicK-mer2 [37] (https://github.com/KiddLab/QuicK-mer2) workflow was run on WGS data to generate copy number estimates, in 1-kbp windows across the genome, in each individual. Briefly, this method utilizes a kmer-based approach to perform copy number estimation. After running QuicK-mer2, we utilized the bedtools [64] map function to calculate the average copy number across the copy number windows covering hs737 (b38, chr10:128,568,604–128,569,741). If the copy number was less than 1.3, we called it as a deletion and if it was greater than 2.7, we called it as a duplication. QuicK-mer2 was run on the autism families in this study and also the high-coverage 1000 genomes project data available as described at http://ftp.1000genomes.ebi.ac.uk/vol1/ftp/data_collections/1000G_2504_high_coverage/.

To assess structural variation in gnomAD v2.1 [40], we queried for our enhancer region on hg19 (10-130366868-130368005) and also available at this link https://gnomad.broadinstitute.org/region/10-130366868-130368005?dataset=gnomad_sv_r2_1. There were “No variants found” in this region.

Statistical testing of coding DNVs in EBF3

DNV data was collected from two recent papers on NDDs [50, 51]. After overlapping samples between the two studies were removed, there were a total of 37,692 sequenced parent-child trios. To test for enrichment of coding DNVs that were loss-of-function or missense, we applied the chimpanzee-human [50] and denovolyzeR [52] models as previously described [50, 65].

Mouse ENCODE chromatin state and interaction tracks

Chromatin state and interaction data from mouse developmental timepoints were assessed in the ENCODE Regulation, ENC+EPD Enhc-Gene, ENCODE cCREs, and EPDnew Promoters tracks in the mm10 genome browser at UCSC [42,43,44,45].

ChIP sequencing assessment of data from Harms et al. [55]

We downloaded ChIP sequencing data from Harms et al. [55] at the following GEO link https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/geo/query/acc.cgi?acc=GSE90682. To identify the promoter locations in the human genome, we looked at sequence 5 kbp upstream of the transcription start site using the Table Browser feature of the UCSC Genome Browser [42]. We then used bedtools [64] intersect to identify which ChIP peaks overlapped with promoters in the human genome. To determine which NDD genes were bound at their promoter, we pulled the genome-wide significant gene lists from Coe et al. [50] and Kaplanis et al. [51] and compared to our EBF3 bound promoter list. Gene Ontology enrichment was performed using the Database for Annotation, Visualization and Integrated Discovery tool version 6.8 (https://david.ncifcrf.gov/) [66].

RNA sequencing at mouse embryonic day 11.5 from ENCODE

Mouse E11.5 forebrain RNAseq data was downloaded from https://www.encodeproject.org/files/ENCFF465SNB/@@download/ENCFF465SNB.tsv, mouse E11.5 midbrain RNAseq data was downloaded from https://www.encodeproject.org/files/ENCFF359ZOA/@@download/ENCFF359ZOA.tsv, and mouse E11.5 hindbrain data was downloaded from https://www.encodeproject.org/files/ENCFF750FTK/@@download/ENCFF750FTK.tsv. For each file, we retained all Ensembl gene identifiers and annotated them to HGNC-approved identifiers using biomart (https://m.ensembl.org/info/data/biomart/index.html). We used an expression cutoff of > 2 to call a gene as expressed and < 2 as not expressed in each region of the brain.

Cell lines

Neuro2a (ATCC CCL-131) were purchased from ATCC and grown under standard conditions (DMEM + 10% FBS and 1% penicillin-streptomycin).

Luciferase assays

Five hundred nanograms of firefly luciferase vector (pGL 4.23, Promega Corporation) containing the enhancer sequence cloned upstream of luc2 and 9 ng of Renilla luciferase vector (transfection control) were transiently transfected into the Neuro2A cell line (1 × 105 cells/well) using 3μl FuGene HD transfection reagent in 100μl of OPTI-MEM medium. Neuro2A cells were incubated for 48 h and luminescence measured using a Dual-Luciferase Reporter Assay System on a Promega GloMax luminometer. All assays were performed in triplicate for a total of six independent readings of each construct. Significance was calculated by using a two-sided t-test assuming unequal variance and used in two situations, first to compare the wild type construct to basal construct and second to compare the variant constructs to the wild type (Supplemental Table S5).

Transcription factor binding predictions

Transcription factor binding analysis was performed using QBiC-Pred [34] and selecting all transcription factor families and using a p-value threshold of 0.0001 and output to a VCF format. Once predictions were obtained, transcription factors were then cross referenced with RNA sequencing from mouse embryonic brains at day 11.5 to identify which transcription factors are highly expressed.

Brainspan RNAseq analysis

RNAseq data (RNA-Seq Gencode v10 summarized to gene (n = 52,377 genes)) was downloaded from http://brainspan.org/static/download.html on 4 February 2021 for the developmental transcriptome which contains samples from 8 weeks after conception to 40 years. Linear regression was performed on the data set for age and brain region using R. Spearman correlations were calculated for each gene (X) and EBF3 (Y) and SFARI genes with an absolute correlation greater than 0.6 were retained for further analysis. Pathway analysis was performed using GeneMania (https://genemania.org) and String-db (https://string-db.org). GOTERM analysis was performed using http://geneontology.org where the input genes were those from cluster positively correlated with EBF3 and SFARI gene list (1/13/21 release) was used as the background.

Statistical assessment of phenotypes

Phenotypes were assessed by counting the number of individuals with a given phenotype. Specifically, we assessed all probands from previous studies for a total of 13 probands in addition to 7 new probands. The 7 or more symptoms category was calculated as the total number of neurological abnormalities, other diagnoses, and the presence of craniofacial abnormalities, where 1 was added to the total number of symptoms if the individual had at least one craniofacial abnormality. A minus (−) value was taken to be the absence of the phenotype whereas NA was taken to be that the assessment was missing. Two probands were missing assessments for ataxia and were excluded from that analysis for a total of 18 individuals rather than 20. A one-sided Fisher’s exact test was used to calculate significance and odds ratio for the phenotypes and only probands who had an assessment for that phenotype were included in the calculation.

Availability of data and materials

The DNV callset is available at SFARI Base accession: SFARI_SSC_WGS_2a, https://base.sfari.org/. We are grateful to all of the families at the participating SSC sites, as well as the principal investigators (A. Beaudet, R. Bernier, J. Constantino, E. Cook, E. Fombonne, D. Geschwind, R. Goin-Kochel, E. Hanson, D. Grice, A. Klin, D. Ledbetter, C. Lord, C. Martin, D. Martin, R. Maxim, J. Miles, O. Ousley, K. Pelphrey, B. Peterson, J. Piggot, C. Saulnier, M. State, W. Stone, J. Sutcliffe, C. Walsh, Z. Warren, and E. Wijsman). We appreciate obtaining access to phenotypic and genetic data on SFARI Base. Approved researchers can obtain the SSC population dataset described in this study (https://www.sfari.org/resource/simons-simplex-collection/) by applying at https://base.sfari.org.

References.

  1. Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, Hunter DJ, et al. Finding the missing heritability of complex diseases. Nature. 2009;461(7265):747–53. https://0-doi-org.brum.beds.ac.uk/10.1038/nature08494.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Pollard KS, Hubisz MJ, Rosenbloom KR, Siepel A. Detection of nonneutral substitution rates on mammalian phylogenies. Genome Res. 2010;20(1):110–21. https://0-doi-org.brum.beds.ac.uk/10.1101/gr.097857.109.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature. 2020;581(7809):434–43. https://0-doi-org.brum.beds.ac.uk/10.1038/s41586-020-2308-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Moore JE, Purcaro MJ, Pratt HE, Epstein CB, Shoresh N, Adrian J, et al. Expanded encyclopaedias of DNA elements in the human and mouse genomes. Nature. 2020;583(7818):699–710. https://0-doi-org.brum.beds.ac.uk/10.1038/s41586-020-2493-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Grove J, Ripke S, Als TD, Mattheisen M, Walters RK, Won H, et al. Identification of common genetic risk variants for autism spectrum disorder. Nat Genet. 2019;51(3):431–44. https://0-doi-org.brum.beds.ac.uk/10.1038/s41588-019-0344-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Yang J, Manolio TA, Pasquale LR, Boerwinkle E, Caporaso N, Cunningham JM, et al. Genome partitioning of genetic variation for complex traits using common SNPs. Nat Genet. 2011;43(6):519–25. https://0-doi-org.brum.beds.ac.uk/10.1038/ng.823.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Maurano MT, Humbert R, Rynes E, Thurman RE, Haugen E, Wang H, et al. Systematic localization of common disease-associated variation in regulatory DNA. Science (New York, NY). 2012;337(6099):1190–5.

    Article  CAS  Google Scholar 

  8. Turner TN, Coe BP, Dickel DE, Hoekzema K, Nelson BJ, Zody MC, et al. Genomic patterns of de novo mutation in simplex autism. Cell. 2017;171(3):710–722.e712.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. An JY, Lin K, Zhu L, Werling DM, Dong S, Brand H, et al. Genome-wide de novo risk score implicates promoter variation in autism spectrum disorder. Science (New York, NY). 2018;362(6420).

  10. Zhou J, Park C, Theesfeld C, Yuan Y, Sawicka K, Darnell J, et al. Whole-genome deep learning analysis reveals causal role of noncoding mutations in autism. bioRxiv. 2018.

  11. Markenscoff-Papadimitriou E, Whalen S, Przytycki P, Thomas R, Binyameen F, Nowakowski TJ, et al. A chromatin accessibility atlas of the developing human telencephalon. Cell. 2020;182(3):754–769.e718.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Iossifov I, O'Roak BJ, Sanders SJ, Ronemus M, Krumm N, Levy D, et al. The contribution of de novo coding mutations to autism spectrum disorder. Nature. 2014;515(7526):216–21. https://0-doi-org.brum.beds.ac.uk/10.1038/nature13908.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. O'Roak BJ, Deriziotis P, Lee C, Vives L, Schwartz JJ, Girirajan S, et al. Exome sequencing in sporadic autism spectrum disorders identifies severe de novo mutations. Nat Genet. 2011;43(6):585–9. https://0-doi-org.brum.beds.ac.uk/10.1038/ng.835.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. O'Roak BJ, Vives L, Girirajan S, Karakoc E, Krumm N, Coe BP, et al. Sporadic autism exomes reveal a highly interconnected protein network of de novo mutations. Nature. 2012;485(7397):246–50. https://0-doi-org.brum.beds.ac.uk/10.1038/nature10989.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. De Rubeis S, He X, Goldberg AP, Poultney CS, Samocha K, Cicek AE, et al. Synaptic, transcriptional and chromatin genes disrupted in autism. Nature. 2014;515(7526):209–15. https://0-doi-org.brum.beds.ac.uk/10.1038/nature13772.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Sanders SJ, Murtha MT, Gupta AR, Murdoch JD, Raubeson MJ, Willsey AJ, et al. De novo mutations revealed by whole-exome sequencing are strongly associated with autism. Nature. 2012;485(7397):237–41. https://0-doi-org.brum.beds.ac.uk/10.1038/nature10945.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Sandin S, Lichtenstein P, Kuja-Halkola R, Hultman C, Larsson H, Reichenberg A. The heritability of autism spectrum disorder. Jama. 2017;318(12):1182–4. https://0-doi-org.brum.beds.ac.uk/10.1001/jama.2017.12141.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Levy D, Ronemus M, Yamrom B, Lee YH, Leotta A, Kendall J, et al. Rare de novo and transmitted copy-number variation in autistic spectrum disorders. Neuron. 2011;70(5):886–97. https://0-doi-org.brum.beds.ac.uk/10.1016/j.neuron.2011.05.015.

    Article  CAS  PubMed  Google Scholar 

  19. Sebat J, Lakshmi B, Malhotra D, Troge J, Lese-Martin C, Walsh T, et al. Strong association of de novo copy number mutations with autism. Science (New York, NY). 2007;316(5823):445–9.

    Article  CAS  Google Scholar 

  20. Sanders SJ, Ercan-Sencicek AG, Hus V, Luo R, Murtha MT, Moreno-De-Luca D, et al. Multiple recurrent de novo CNVs, including duplications of the 7q11.23 Williams syndrome region, are strongly associated with autism. Neuron. 2011;70(5):863–85. https://0-doi-org.brum.beds.ac.uk/10.1016/j.neuron.2011.05.002.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Girirajan S, Dennis MY, Baker C, Malig M, Coe BP, Campbell CD, et al. Refinement and discovery of new hotspots of copy-number variation associated with autism spectrum disorder. Am J Hum Genet. 2013;92(2):221–37. https://0-doi-org.brum.beds.ac.uk/10.1016/j.ajhg.2012.12.016.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Iossifov I, Ronemus M, Levy D, Wang Z, Hakker I, Rosenbaum J, et al. De novo gene disruptions in children on the autistic spectrum. Neuron. 2012;74(2):285–99. https://0-doi-org.brum.beds.ac.uk/10.1016/j.neuron.2012.04.009.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Neale BM, Kou Y, Liu L, Ma'ayan A, Samocha KE, Sabo A, et al. Patterns and rates of exonic de novo mutations in autism spectrum disorders. Nature. 2012;485(7397):242–5. https://0-doi-org.brum.beds.ac.uk/10.1038/nature11011.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Turner TN, Hormozdiari F, Duyzend MH, McClymont SA, Hook PW, Iossifov I, et al. Genome sequencing of autism-affected families reveals disruption of putative noncoding regulatory DNA. Am J Hum Genet. 2016;98(1):58–74. https://0-doi-org.brum.beds.ac.uk/10.1016/j.ajhg.2015.11.023.

    Article  CAS  PubMed  Google Scholar 

  25. Brandler WM, Antaki D, Gujral M, Kleiber ML, Whitney J, Maile MS, et al. Paternally inherited cis-regulatory structural variants are associated with autism. Science (New York, NY). 2018;360(6386):327–31.

    Article  CAS  Google Scholar 

  26. Frazer KA, Pachter L, Poliakov A, Rubin EM, Dubchak I. VISTA: computational tools for comparative genomics. Nucleic Acids Res. 2004;32(Web Server issue):W273–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Visel A, Minovitsky S, Dubchak I, Pennacchio LA. VISTA Enhancer Browser--a database of tissue-specific human enhancers. Nucleic Acids Res. 2007;35(Database issue):D88–92. https://0-doi-org.brum.beds.ac.uk/10.1093/nar/gkl822.

    Article  CAS  PubMed  Google Scholar 

  28. Kvon EZ, Zhu Y, Kelman G, Novak CS, Plajzer-Frick I, Kato M, et al. Comprehensive in vivo interrogation reveals phenotypic impact of human enhancer variants. Cell. 2020;180(6):1262–1271.e1215.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Jiang Y, Han Y, Petrovski S, Owzar K, Goldstein DB, Allen AS. Incorporating functional information in tests of excess de novo mutational load. Am J Hum Genet. 2015;97(2):272–83. https://0-doi-org.brum.beds.ac.uk/10.1016/j.ajhg.2015.06.013.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Wilfert AB, Turner TN, Murali SC, Hsieh P, Sulovari A, Wang T, et al. Recent ultra-rare inherited mutations identify novel autism candidate risk genes. Nat Genet. in press.

  31. Weiner DJ, Wigdor EM, Ripke S, Walters RK, Kosmicki JA, Grove J, et al. Polygenic transmission disequilibrium confirms that common and rare variation act additively to create risk for autism spectrum disorders. Nat Genet. 2017;49(7):978–85. https://0-doi-org.brum.beds.ac.uk/10.1038/ng.3863.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Chatterjee S, Kapoor A, Akiyama JA, Auer DR, Lee D, Gabriel S, et al. Enhancer variants synergistically drive dysfunction of a gene regulatory network in Hirschsprung disease. Cell. 2016;167(2):355–368.e310.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Zhao J, Li D, Seo J, Allen AS, Gordân R. Quantifying the impact of non-coding variants on transcription factor-DNA binding. Res Comput Mol Biol. 2017;10229:336–52. https://0-doi-org.brum.beds.ac.uk/10.1007/978-3-319-56970-3_21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Martin V, Zhao J, Afek A, Mielko Z, Gordân R. QBiC-Pred: quantitative predictions of transcription factor binding changes due to sequence variants. Nucleic Acids Res. 2019;47(W1):W127–w135. https://0-doi-org.brum.beds.ac.uk/10.1093/nar/gkz363.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Coe BP, Witherspoon K, Rosenfeld JA, van Bon BW, Vulto-van Silfhout AT, Bosco P, et al. Refining analyses of copy number variation identifies specific genes associated with developmental delay. Nat Genet. 2014;46(10):1063–71. https://0-doi-org.brum.beds.ac.uk/10.1038/ng.3092.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Cooper GM, Coe BP, Girirajan S, Rosenfeld JA, Vu TH, Baker C, et al. A copy number variation morbidity map of developmental delay. Nat Genet. 2011;43(9):838–46. https://0-doi-org.brum.beds.ac.uk/10.1038/ng.909.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Shen F, Kidd JM. Rapid, paralog-sensitive CNV analysis of 2457 human genomes using QuicK-mer2. Genes (Basel). 2020;11(2).

  38. Sanders SJ, He X, Willsey AJ, Ercan-Sencicek AG, Samocha KE, Cicek AE, et al. Insights into autism spectrum disorder genomic architecture and biology from 71 risk loci. Neuron. 2015;87(6):1215–33. https://0-doi-org.brum.beds.ac.uk/10.1016/j.neuron.2015.09.016.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Byrska-Bishop M, Evani US, Zhao X, Basile AO, Abel HJ, Regier AA, et al. High coverage whole genome sequencing of the expanded 1000 Genomes Project cohort including 602 trios. bioRxiv. 2021; 2021.2002.2006.430068.

  40. Collins RL, Brand H, Karczewski KJ, Zhao X, Alföldi J, Francioli LC, et al. A structural variation reference for medical and population genetics. Nature. 2020;581(7809):444–51. https://0-doi-org.brum.beds.ac.uk/10.1038/s41586-020-2287-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Reilly SK, Yin J, Ayoub AE, Emera D, Leng J, Cotney J, et al. Evolutionary genomics. Evolutionary changes in promoter and enhancer activity during human corticogenesis. Science (New York, NY). 2015;347(6226):1155–9.

    Article  CAS  Google Scholar 

  42. Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, et al. The human genome browser at UCSC. Genome Res. 2002;12(6):996–1006. https://0-doi-org.brum.beds.ac.uk/10.1101/gr.229102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Gorkin DU, Barozzi I, Zhao Y, Zhang Y, Huang H, Lee AY, et al. An atlas of dynamic chromatin landscapes in mouse fetal development. Nature. 2020;583(7818):744–51. https://0-doi-org.brum.beds.ac.uk/10.1038/s41586-020-2093-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Ernst J, Kellis M. ChromHMM: automating chromatin-state discovery and characterization. Nat Methods. 2012;9(3):215–6. https://0-doi-org.brum.beds.ac.uk/10.1038/nmeth.1906.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Sloan CA, Chan ET, Davidson JM, Malladi VS, Strattan JS, Hitz BC, et al. ENCODE data at the ENCODE portal. Nucleic Acids Res. 2016;44(D1):D726–32. https://0-doi-org.brum.beds.ac.uk/10.1093/nar/gkv1160.

    Article  CAS  PubMed  Google Scholar 

  46. Dixon JR, Selvaraj S, Yue F, Kim A, Li Y, Shen Y, et al. Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature. 2012;485(7398):376–80. https://0-doi-org.brum.beds.ac.uk/10.1038/nature11082.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Bonev B, Mendelson Cohen N, Szabo Q, Fritsch L, Papadopoulos GL, Lubling Y, et al. Multiscale 3D genome rewiring during mouse neural development. Cell. 2017;171(3):557–572.e524.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Durand NC, Shamim MS, Machol I, Rao SS, Huntley MH, Lander ES, et al. Juicer provides a one-click system for analyzing loop-resolution Hi-C experiments. Cell Syst. 2016;3(1):95–8. https://0-doi-org.brum.beds.ac.uk/10.1016/j.cels.2016.07.002.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Won H, de la Torre-Ubieta L, Stein JL, Parikshak NN, Huang J, Opland CK, et al. Chromosome conformation elucidates regulatory relationships in developing human brain. Nature. 2016;538(7626):523–7. https://0-doi-org.brum.beds.ac.uk/10.1038/nature19847.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Coe BP, Stessman HAF, Sulovari A, Geisheker MR, Bakken TE, Lake AM, et al. Neurodevelopmental disease genes implicated by de novo mutation and copy number variation morbidity. Nat Genet. 2019;51(1):106–16. https://0-doi-org.brum.beds.ac.uk/10.1038/s41588-018-0288-4.

    Article  CAS  PubMed  Google Scholar 

  51. Kaplanis J, Samocha KE, Wiel L, Zhang Z, Arvai KJ, Eberhardt RY, et al. Evidence for 28 genetic disorders discovered by combining healthcare and research data. Nature. 2020;586(7831):757–62. https://0-doi-org.brum.beds.ac.uk/10.1038/s41586-020-2832-5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Ware JS, Samocha KE, Homsy J, Daly MJ. Interpreting de novo variation in human disease using denovolyzeR. Current protocols in human genetics. 2015;87:7.25.21–15.

    Article  Google Scholar 

  53. Sleven H, Welsh SJ, Yu J, Churchill MEA, Wright CF, Henderson A, et al. De novo mutations in EBF3 cause a neurodevelopmental syndrome. Am J Hum Genet. 2017;100(1):138–50. https://0-doi-org.brum.beds.ac.uk/10.1016/j.ajhg.2016.11.020.

    Article  CAS  PubMed  Google Scholar 

  54. Chao HT, Davids M, Burke E, Pappas JG, Rosenfeld JA, McCarty AJ, et al. A syndromic neurodevelopmental disorder caused by de novo variants in EBF3. Am J Hum Genet. 2017;100(1):128–37. https://0-doi-org.brum.beds.ac.uk/10.1016/j.ajhg.2016.11.018.

    Article  CAS  PubMed  Google Scholar 

  55. Harms FL, Girisha KM, Hardigan AA, Kortüm F, Shukla A, Alawi M, et al. Mutations in EBF3 disturb transcriptional profiles and cause intellectual disability, ataxia, and facial dysmorphism. Am J Hum Genet. 2017;100(1):117–27. https://0-doi-org.brum.beds.ac.uk/10.1016/j.ajhg.2016.11.012.

    Article  CAS  PubMed  Google Scholar 

  56. Hormozdiari F, Penn O, Borenstein E, Eichler EE. The discovery of integrated gene networks for autism and related disorders. Genome Res. 2015;25(1):142–54. https://0-doi-org.brum.beds.ac.uk/10.1101/gr.178855.114.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Fischbach GD, Lord C. The Simons Simplex Collection: a resource for identification of autism genetic risk factors. Neuron. 2010;68(2):192–5. https://0-doi-org.brum.beds.ac.uk/10.1016/j.neuron.2010.10.006.

    Article  CAS  PubMed  Google Scholar 

  58. Celestino-Soper PB, Shaw CA, Sanders SJ, Li J, Murtha MT, Ercan-Sencicek AG, et al. Use of array CGH to detect exonic copy number variants throughout the genome in autism families detects a novel deletion in TMLHE. Hum Mol Genet. 2011;20(22):4360–70. https://0-doi-org.brum.beds.ac.uk/10.1093/hmg/ddr363.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Satterstrom FK, Kosmicki JA, Wang J, Breen MS, De Rubeis S, An JY, et al. Large-scale exome sequencing study implicates both developmental and functional changes in the neurobiology of autism. Cell. 2020;180(3):568–584.e523.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Gaugler T, Klei L, Sanders SJ, Bodea CA, Goldberg AP, Lee AB, et al. Most genetic risk for autism resides with common variation. Nat Genet. 2014;46(8):881–5. https://0-doi-org.brum.beds.ac.uk/10.1038/ng.3039.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Bernier R, Golzio C, Xiong B, Stessman HA, Coe BP, Penn O, et al. Disruptive CHD8 mutations define a subtype of autism early in development. Cell. 2014;158(2):263–76. https://0-doi-org.brum.beds.ac.uk/10.1016/j.cell.2014.06.017.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Earl RK, Turner TN, Mefford HC, Hudac CM, Gerdts J, Eichler EE, et al. Clinical phenotype of ASD-associated DYRK1A haploinsufficiency. Molecular autism. 2017;8(1):54. https://0-doi-org.brum.beds.ac.uk/10.1186/s13229-017-0173-5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Lonsdale J, Thomas J, Salvatore M, Phillips R, Lo E, Shad S, et al. The Genotype-Tissue Expression (GTEx) project. Nat Genet. 2013;45(6):580–5.

    Article  CAS  Google Scholar 

  64. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics (Oxford, England). 2010;26(6):841–2.

    Article  CAS  Google Scholar 

  65. Turner TN, Wilfert AB, Bakken TE, Bernier RA, Pepper MR, Zhang Z, et al. Sex-based analysis of de novo variants in neurodevelopmental disorders. Am J Hum Genet. 2019;105(6):1274–85. https://0-doi-org.brum.beds.ac.uk/10.1016/j.ajhg.2019.11.003.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.

    Article  PubMed  Google Scholar 

  67. Ren J, Wen L, Gao X, Jin C, Xue Y, Yao X. DOG 1.0: illustrator of protein domain structures. Cell Res. 2009;19(2):271–3. https://0-doi-org.brum.beds.ac.uk/10.1038/cr.2009.6.

    Article  CAS  PubMed  Google Scholar 

  68. Tokheim C, Bhattacharya R, Niknafs N, Gygax DM, Kim R, Ryan M, et al. Exome-scale discovery of hotspot mutation regions in human cancer using 3D protein structure. Cancer Res. 2016;76(13):3719–31. https://0-doi-org.brum.beds.ac.uk/10.1158/0008-5472.CAN-15-3190.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Niknafs N, Kim D, Kim R, Diekhans M, Ryan M, Stenson PD, et al. MuPIT interactive: webserver for mapping variant positions to annotated, interactive 3D structures. Hum Genet. 2013;132(11):1235–43. https://0-doi-org.brum.beds.ac.uk/10.1007/s00439-013-1325-0.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Thank you to Dr. Evan Eichler, Dr. Jay Shendure, Dr. Barak Cohen, and Dr. Jeffrey Kidd for helpful discussions on this work.

Code availability

Codes are available from previous publications as described in the methods. Our modified fitDNM code that works on noncoding regions is available at https://github.com/TNTurnerLab/fitDNM.

Funding

This work was supported by grants from the National Institutes of Health (R00MH117165 to T.N.T., UM1HG008901 to M.C.Z.). The Centers for Common Disease Genomics are funded by the National Human Genome Research Institute and the National Heart, Lung, and Blood Institute and the GSP Coordinating Center (U24HG008956) contributed to cross-program scientific initiatives and provided logistical and general study coordination. Research conducted at the E.O. Lawrence Berkeley National Laboratory was additionally supported by NIH grants (R01HG003988) (to L.A.P. and D.E.D.) and performed under a Department of Energy Contract (DE-AC02-05CH11231), University of California.

Author information

Authors and Affiliations

Authors

Contributions

E.M.P. and T.N.T. designed the study; E.M.P. and T.N.T. wrote the paper; E.M.P., T.J.H., M.B., N.S., A.S.A., M.C.Z., and T.N.T. assessed the DNV data; T.J.H., N.S., and A.S.A. optimized the fitDNM model for noncoding regions; E.M.P. and T.N.T. assessed the morbidity map data; E.M.P., R.M., G.N., A.A., J.K., J.N., M.C.Z., and T.N.T. generated the QuicK-mer2 files and performed genotyping; E.M.P. and T.N.T. performed the statistical assessment of protein-coding DNVs in EBF3; E.M.P., Z.C., D.U.G., and T.N.T. assessed the mouse ENCODE data; E.M.P. and T.N.T. assessed the ChIP sequencing data; E.M.P. and T.N.T. assessed the RNAseq data; E.M.P., S.C., L.F., and T.N.T. designed and implemented the luciferase assays; E.M.P. and T.N.T. performed transcription factor binding analysis; E.M.P., Z.C., and D.U.G. performed the Hi-C analysis; E.M.P, R.A.B., and T.N.T. assessed phenotypes of individuals with hs737 variants; E.M.P., M.W., L.P., S.R., B.C., U.K., A.S., L.M.C., F.C., S.E., F.R, S.M., R.M., H.H., R.A.B., and T.N.T. assessed phenotypes of individuals with protein-coding DNVs in EBF3. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Tychele N. Turner.

Ethics declarations

Ethics approval and consent to participate

This study was reviewed by the Washington University in St. Louis Human Research Protection Office and, since all of the human genomic data is completely de-identified, determined it was not subject to Institutional Review Board (IRB) oversight.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1

: Coordinates in the human genome (build 38) of VISTA enhancers driving expression in the brain.

Additional file 2: Fig S1

. Polygenic risk scores (PRS) for the three individuals with hs737 mutations. Fig. S2. Zoom in on the hs737 enhancer with annotations from other datasets. The enhancer is in a PsychEncode fetal enhancer, contains central nervous system DNaseI hypersensitive sites, contains conserved transcription factor binding sites, and is highly conserved across the vertebrate lineage. Also shown are the locations of the hs737 de novo mutations identified in individuals with autism. Fig S3. String-db network analysis, predicts interactions between EBF3 and SPAST, CSNK2A1, HNRNPU, CHD7, KDM6A, KDM6B. Table S2. fitDNM results for de novo mutations in VISTA enhancers driving brain expression. Table S3. Other de novo SNVs/indels seen in individuals with hs737 enhancer mutations. Table S4. Other copy number variation in individuals with hs737 de novo mutations. Table S5. Statistical significance calculations for the luciferase assays. Table S6. Expression of transcription factors potentially binding at variant locations in the hs737 enhancers.

Additional file 3: Table S7:  

should have the label "CNVs over hs737 in the morbidity map dataset.

Additional file 4: Table S8

: Protein-coding de novo mutations in EBF3 from Coe et al. 2019, Nature Genetics and Kaplanis et al. 2020, bioRxiv.

Additional file 5: Table S9: 

should have the label "NDD genes bound by EBF3.

Additional file 6: Table S10: 

should have the label "Phenotypes of individuals with coding and noncoding variants affecting EBF3.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Padhi, E.M., Hayeck, T.J., Cheng, Z. et al. Coding and noncoding variants in EBF3 are involved in HADDS and simplex autism. Hum Genomics 15, 44 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s40246-021-00342-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s40246-021-00342-3

Keywords