Interpreting de novo Variation in Human Disease Using denovolyzeR

Variants in BSN , encoding the presynaptic protein Bassoon, result in a novel neurodevelopmental disorder with a broad phenotypic range

Disease-causing variants in synaptic function genes are a common cause of neurodevelopmental disorders and epilepsy. Here, we describe 14 individuals with de novo disruptive variants in BSN , which encodes the presynaptic protein Bassoon. To expand the phenotypic spectrum, we identified 15 additional individuals with protein-truncating variants (PTVs) from large biobanks. Clinical features were standardized using the Human Phenotype Ontology (HPO) across all 29 individuals, which revealed common clinical characteristics including epilepsy (13/29 45%), febrile seizures (7/29 25%), generalized tonic-clonic seizures (5/29 17%), and focal onset seizures (3/29 10%). Behavioral phenotypes were present in almost half of all individuals (14/29 48%), which comprised ADHD (7/29 25%) and autistic behavior (5/29 17%). Additional common features included developmental delay (11/29 38%), obesity (10/29 34%), and delayed speech (8/29 28%). In adults with BSN PTVs, milder features were common, suggesting phenotypic variability including a range of individuals without obvious neurodevelopmental features (7/29 24%). To detect gene-specific signatures, we performed association analysis in a cohort of 14,895 individuals with neurodevelopmental disorders (NDDs). A total of 66 clinical features were associated with BSN , including febrile seizures (p=1.26e-06) and behavioral disinhibition (p = 3.39e-17). Furthermore, individuals carrying BSN variants were phenotypically more similar than expected by chance (p=0.00014), exceeding phenotypic relatedness in 179/256 NDD-related conditions. In summary, integrating information derived from community-based gene matching and large data repositories through computational phenotyping approaches, we identify BSN variants as the cause of a new class of synaptic disorder with a broad phenotypic range across the age spectrum.

TRIM71 mutations cause a neurodevelopmental syndrome featuring ventriculomegaly and hydrocephalus

Congenital hydrocephalus, characterized by cerebral ventriculomegaly, is one of the most common reasons for paediatric brain surgery. Recent studies have implicated lin-41 (lineage variant 41)/TRIM71 (tripartite motif 71) as a candidate congenital hydrocephalus risk gene; however, TRIM71 variants have not been systematically examined in a large patient cohort or conclusively linked with an OMIM syndrome. Through cross-sectional analysis of the largest assembled cohort of patients with cerebral ventriculomegaly, including neurosurgically-treated congenital hydrocephalus (totalling 2697 parent-proband trios and 8091 total exomes), we identified 13 protein-altering de novo variants (DNVs) in TRIM71 in unrelated children exhibiting variable ventriculomegaly, congenital hydrocephalus, developmental delay, dysmorphic features and other structural brain defects, including corpus callosum dysgenesis and white matter hypoplasia. Eight unrelated patients were found to harbour arginine variants, including two recurrent missense DNVs, at homologous positions in RPXGV motifs of different NHL domains. Seven patients with rare, damaging, unphased or transmitted variants of uncertain significance were also identified. NHL-domain variants of TRIM71 exhibited impaired binding to the canonical TRIM71 target CDKN1A; other variants failed to direct the subcellular localization of TRIM71 to processing bodies. Single-cell transcriptomic analysis of human embryos revealed expression of TRIM71 in early first-trimester neural stem cells of the brain. These data show TRIM71 is essential for human brain morphogenesis and that TRIM71 mutations cause a novel neurodevelopmental syndrome that we term 'TRIM71-associated developmental disorders (TADD)', featuring variable ventriculomegaly, congenital hydrocephalus and other structural brain defects.

Burden re-analysis of neurodevelopmental disorder cohorts for prioritization of candidate genes

This study aimed to uncover novel genes associated with neurodevelopmental disorders (NDD) by leveraging recent large-scale de novo burden analysis studies to enhance a virtual gene panel used in a diagnostic setting. We re-analyzed historical trio-exome sequencing data from 745 individuals with NDD according to the most recent diagnostic standards, resulting in a cohort of 567 unsolved individuals. Next, we designed a virtual gene panel containing candidate genes from three large de novo burden analysis studies in NDD and prioritized candidate genes by stringent filtering for ultra-rare de novo variants with high pathogenicity scores. Our analysis revealed an increased burden of de novo variants in our selected candidate genes within the unsolved NDD cohort and identified qualifying de novo variants in seven candidate genes: RIF1, CAMK2D, RAB11FIP4, AGO3, PCBP2, LEO1, and VCP. Clinical data were collected from six new individuals with de novo or inherited LEO1 variants and three new individuals with de novo PCBP2 variants. Our findings add additional evidence for LEO1 as a risk gene for autism and intellectual disability. Furthermore, we prioritize PCBP2 as a candidate gene for NDD associated with motor and language delay. In summary, by leveraging de novo burden analysis studies, employing a stringent variant filtering pipeline, and engaging in targeted patient recruitment, our study contributes to the identification of novel genes implicated in NDDs.

Loss-of-Function Variants in CUL3 Cause a Syndromic Neurodevelopmental Disorder

OBJECTIVE

De novo variants in cullin-3 ubiquitin ligase (CUL3) have been strongly associated with neurodevelopmental disorders (NDDs), but no large case series have been reported so far. Here, we aimed to collect sporadic cases carrying rare variants in CUL3, describe the genotype-phenotype correlation, and investigate the underlying pathogenic mechanism.

METHODS

Genetic data and detailed clinical records were collected via multicenter collaboration. Dysmorphic facial features were analyzed using GestaltMatcher. Variant effects on CUL3 protein stability were assessed using patient-derived T-cells.

RESULTS

We assembled a cohort of 37 individuals with heterozygous CUL3 variants presenting a syndromic NDD characterized by intellectual disability with or without autistic features. Of these, 35 have loss-of-function (LoF) and 2 have missense variants. CUL3 LoF variants in patients may affect protein stability leading to perturbations in protein homeostasis, as evidenced by decreased ubiquitin-protein conjugates in vitro. Notably, we show that 4E-BP1 (EIF4EBP1), a prominent substrate of CUL3, fails to be targeted for proteasomal degradation in patient-derived cells.

INTERPRETATION

Our study further refines the clinical and mutational spectrum of CUL3-associated NDDs, expands the spectrum of cullin RING E3 ligase-associated neuropsychiatric disorders, and suggests haploinsufficiency via LoF variants is the predominant pathogenic mechanism. ANN NEUROL 2024.

Genetic constraint at single amino acid resolution in protein domains improves missense variant prioritisation and gene discovery

BACKGROUND

One of the major hurdles in clinical genetics is interpreting the clinical consequences associated with germline missense variants in humans. Recent significant advances have leveraged natural variation observed in large-scale human populations to uncover genes or genomic regions that show a depletion of natural variation, indicative of selection pressure. We refer to this as "genetic constraint". Although existing genetic constraint metrics have been demonstrated to be successful in prioritising genes or genomic regions associated with diseases, their spatial resolution is limited in distinguishing pathogenic variants from benign variants within genes.

METHODS

We aim to identify missense variants that are significantly depleted in the general human population. Given the size of currently available human populations with exome or genome sequencing data, it is not possible to directly detect depletion of individual missense variants, since the average expected number of observations of a variant at most positions is less than one. We instead focus on protein domains, grouping homologous variants with similar functional impacts to examine the depletion of natural variations within these comparable sets. To accomplish this, we develop the Homologous Missense Constraint (HMC) score. We utilise the Genome Aggregation Database (gnomAD) 125 K exome sequencing data and evaluate genetic constraint at quasi amino-acid resolution by combining signals across protein homologues.

RESULTS

We identify one million possible missense variants under strong negative selection within protein domains. Though our approach annotates only protein domains, it nonetheless allows us to assess 22% of the exome confidently. It precisely distinguishes pathogenic variants from benign variants for both early-onset and adult-onset disorders. It outperforms existing constraint metrics and pathogenicity meta-predictors in prioritising de novo mutations from probands with developmental disorders (DD). It is also methodologically independent of these, adding power to predict variant pathogenicity when used in combination. We demonstrate utility for gene discovery by identifying seven genes newly significantly associated with DD that could act through an altered-function mechanism.

CONCLUSIONS

Grouping variants of comparable functional impacts is effective in evaluating their genetic constraint. HMC is a novel and accurate predictor of missense consequence for improved variant interpretation.

Assisting the analysis of insertions and deletions using regional allele frequencies

Accurate estimation of population allele frequency (AF) is crucial for gene discovery and genetic diagnostics. However, determining AF for frameshift-inducing small insertions and deletions (indels) faces challenges due to discrepancies in mapping and variant calling methods. Here, we propose an innovative approach to assess indel AF. We developed CRAFTS-indels (Calculating Regional Allele Frequency Targeting Small indels), an algorithm that combines AF of distinct indels within a given region and provides "regional AF" (rAF). We tested and validated CRAFTS-indels using three independent datasets: gnomAD v2 (n=125,748 samples), an internal dataset (IGM; n=39,367), and the UK BioBank (UKBB; n=469,835). By comparing rAF against standard AF, we identified rare indels with rAF exceeding standard AF (sAF≤10-4 and rAF>10-4) as "rAF-hi" indels. Notably, a high percentage of rare indels were "rAF-hi", with a higher proportion in gnomAD v2 (11-20%) and IGM (11-22%) compared to the UKBB (5-9% depending on the CRAFTS-indels' parameters). Analysis of the overlap of regions based on their rAF with low complexity regions and with ClinVar classification supported the pertinence of rAF. Using the internal dataset, we illustrated the utility of CRAFTS-indel in the analysis of de novo variants and the potential negative impact of rAF-hi indels in gene discovery. In summary, annotation of indels with cohort specific rAF can be used to handle some of the limitations of current annotation pipelines and facilitate detection of novel gene disease associations. CRAFTS-indels offers a user-friendly approach to providing rAF annotation. It can be integrated into public databases such as gnomAD, UKBB and used by ClinVar to revise indel classifications.

Expanding the genetics and phenotypes of ocular congenital cranial dysinnervation disorders

Purpose

To identify genetic etiologies and genotype/phenotype associations for unsolved ocular congenital cranial dysinnervation disorders (oCCDDs).

Methods

We coupled phenotyping with exome or genome sequencing of 467 pedigrees with genetically unsolved oCCDDs, integrating analyses of pedigrees, human and animal model phenotypes, and de novo variants to identify rare candidate single nucleotide variants, insertion/deletions, and structural variants disrupting protein-coding regions. Prioritized variants were classified for pathogenicity and evaluated for genotype/phenotype correlations.

Results

Analyses elucidated phenotypic subgroups, identified pathogenic/likely pathogenic variant(s) in 43/467 probands (9.2%), and prioritized variants of uncertain significance in 70/467 additional probands (15.0%). These included known and novel variants in established oCCDD genes, genes associated with syndromes that sometimes include oCCDDs (e.g., MYH10, KIF21B, TGFBR2, TUBB6), genes that fit the syndromic component of the phenotype but had no prior oCCDD association (e.g., CDK13, TGFB2), genes with no reported association with oCCDDs or the syndromic phenotypes (e.g., TUBA4A, KIF5C, CTNNA1, KLB, FGF21), and genes associated with oCCDD phenocopies that had resulted in misdiagnoses.

Conclusion

This study suggests that unsolved oCCDDs are clinically and genetically heterogeneous disorders often overlapping other Mendelian conditions and nominates many candidates for future replication and functional studies.

Statistical methods for assessing the effects of de novo variants on birth defects

With the development of next-generation sequencing technology, de novo variants (DNVs) with deleterious effects can be identified and investigated for their effects on birth defects such as congenital heart disease (CHD). However, statistical power is still limited for such studies because of the small sample size due to the high cost of recruiting and sequencing samples and the low occurrence of DNVs. DNV analysis is further complicated by genetic heterogeneity across diseased individuals. Therefore, it is critical to jointly analyze DNVs with other types of genomic/biological information to improve statistical power to identify genes associated with birth defects. In this review, we discuss the general workflow, recent developments in statistical methods, and future directions for DNV analysis.

Proteome-Wide Assessment of Clustering of Missense Variants in Neurodevelopmental Disorders Versus Cancer

Missense de novo variants (DNVs) and missense somatic variants contribute to neurodevelopmental disorders (NDDs) and cancer, respectively. Proteins with statistical enrichment based on analyses of these variants exhibit convergence in the differing NDD and cancer phenotypes. Herein, the question of why some of the same proteins are identified in both phenotypes is examined through investigation of clustering of missense variation at the protein level. Our hypothesis is that missense variation is present in different protein locations in the two phenotypes leading to the distinct phenotypic outcomes. We tested this hypothesis in 1D protein space using our software CLUMP. Furthermore, we newly developed 3D-CLUMP that uses 3D protein structures to spatially test clustering of missense variation for proteome-wide significance. We examined missense DNVs in 39,883 parent-child sequenced trios with NDDs and missense somatic variants from 10,543 sequenced tumors covering five TCGA cancer types and two COSMIC pan-cancer aggregates of tissue types. There were 57 proteins with proteome-wide significant missense variation clustering in NDDs when compared to cancers and 79 proteins with proteome-wide significant missense clustering in cancers compared to NDDs. While our main objective was to identify differences in patterns of missense variation, we also identified a novel NDD protein BLTP2. Overall, our study is innovative, provides new insights into differential missense variation in NDDs and cancer at the protein-level, and contributes necessary information toward building a framework for thinking about prognostic and therapeutic aspects of these proteins.

Mutation of key signaling regulators of cerebrovascular development in vein of Galen malformations

To elucidate the pathogenesis of vein of Galen malformations (VOGMs), the most common and most severe of congenital brain arteriovenous malformations, we performed an integrated analysis of 310 VOGM proband-family exomes and 336,326 human cerebrovasculature single-cell transcriptomes. We found the Ras suppressor p120 RasGAP (RASA1) harbored a genome-wide significant burden of loss-of-function de novo variants (2042.5-fold, p = 4.79 x 10-7). Rare, damaging transmitted variants were enriched in Ephrin receptor-B4 (EPHB4) (17.5-fold, p = 1.22 x 10-5), which cooperates with p120 RasGAP to regulate vascular development. Additional probands had damaging variants in ACVRL1, NOTCH1, ITGB1, and PTPN11. ACVRL1 variants were also identified in a multi-generational VOGM pedigree. Integrative genomic analysis defined developing endothelial cells as a likely spatio-temporal locus of VOGM pathophysiology. Mice expressing a VOGM-specific EPHB4 kinase-domain missense variant (Phe867Leu) exhibited disrupted developmental angiogenesis and impaired hierarchical development of arterial-capillary-venous networks, but only in the presence of a "second-hit" allele. These results illuminate human arterio-venous development and VOGM pathobiology and have implications for patients and their families.

Genetic dysregulation of an endothelial Ras signaling network in vein of Galen malformations

To elucidate the pathogenesis of vein of Galen malformations (VOGMs), the most common and severe congenital brain arteriovenous malformation, we performed an integrated analysis of 310 VOGM proband-family exomes and 336,326 human cerebrovasculature single-cell transcriptomes. We found the Ras suppressor p120 RasGAP ( RASA1 ) harbored a genome-wide significant burden of loss-of-function de novo variants (p=4.79×10 ^-7 ). Rare, damaging transmitted variants were enriched in Ephrin receptor-B4 ( EPHB4 ) (p=1.22×10 ^-5 ), which cooperates with p120 RasGAP to limit Ras activation. Other probands had pathogenic variants in ACVRL1 , NOTCH1 , ITGB1 , and PTPN11 . ACVRL1 variants were also identified in a multi-generational VOGM pedigree. Integrative genomics defined developing endothelial cells as a key spatio-temporal locus of VOGM pathophysiology. Mice expressing a VOGM-specific EPHB4 kinase-domain missense variant exhibited constitutive endothelial Ras/ERK/MAPK activation and impaired hierarchical development of angiogenesis-regulated arterial-capillary-venous networks, but only when carrying a "second-hit" allele. These results illuminate human arterio-venous development and VOGM pathobiology and have clinical implications.

Integrated gene analyses of de novo variants from 46,612 trios with autism and developmental disorders

Most genetic studies consider autism spectrum disorder (ASD) and developmental disorder (DD) separately despite overwhelming comorbidity and shared genetic etiology. Here, we analyzed de novo variants (DNVs) from 15,560 ASD (6,557 from SPARK) and 31,052 DD trios independently and also combined as broader neurodevelopmental disorders (NDDs) using three models. We identify 615 NDD candidate genes (false discovery rate [FDR] < 0.05) supported by ≥1 models, including 138 reaching Bonferroni exome-wide significance (P < 3.64e-7) in all models. The genes group into five functional networks associating with different brain developmental lineages based on single-cell nuclei transcriptomic data. We find no evidence for ASD-specific genes in contrast to 18 genes significantly enriched for DD. There are 53 genes that show mutational bias, including enrichments for missense (n = 41) or truncating (n = 12) DNVs. We also find 10 genes with evidence of male- or female-bias enrichment, including 4 X chromosome genes with significant female burden (DDX3X, MECP2, WDR45, and HDAC8). This large-scale integrative analysis identifies candidates and functional subsets of NDD genes.

Germline mosaicism of a missense variant in KCNC2 in a multiplex family with autism and epilepsy characterized by long-read sequencing

Currently, protein-coding de novo variants and large copy number variants have been identified as important for ~30% of individuals with autism. One approach to identify relevant variation in individuals who lack these types of events is by utilizing newer genomic technologies. In this study, highly accurate PacBio HiFi long-read sequencing was applied to a family with autism, epileptic encephalopathy, cognitive impairment, and mild dysmorphic features (two affected female siblings, unaffected parents, and one unaffected male sibling) with no known clinical variant. From our long-read sequencing data, a de novo missense variant in the KCNC2 gene (encodes Kv3.2) was identified in both affected children. This variant was phased to the paternal chromosome of origin and is likely a germline mosaic. In silico assessment revealed the variant was not in controls, highly conserved, and predicted damaging. This specific missense variant (Val473Ala) has been shown in both an ortholog and paralog of Kv3.2 to accelerate current decay, shift the voltage dependence of activation, and prevent the channel from entering a long-lasting open state. Seven additional missense variants have been identified in other individuals with neurodevelopmental disorders (p = 1.03 × 10^-5 ). KCNC2 is most highly expressed in the brain; in particular, in the thalamus and is enriched in GABAergic neurons. Long-read sequencing was useful in discovering the relevant variant in this family with autism that had remained a mystery for several years and will potentially have great benefits in the clinic once it is widely available.

Genetic Influence on Neurodevelopment in Nonsyndromic Craniosynostosis

BACKGROUND

Nonsyndromic craniosynostosis is one of the most common anomalies treated by craniofacial surgeons. Despite optimal surgical management, nearly half of affected children have subtle neurocognitive deficits. Whereas timing and type of surgical intervention have been studied, the possibility of genetic influence on neurodevelopment in nonsyndromic craniosynostosis patients remains unexplored.

METHODS

The authors performed whole-exome sequencing for 404 case-parent trios with sporadic nonsyndromic craniosynostosis. Statistical analyses were performed to assess the burden of de novo mutations in cases compared to both expectation and 1789 healthy control trios. Individuals with and without each mutation class were analyzed, and the presence or absence of various types of neurodevelopmental delay were recorded alongside demographic information.

RESULTS

The authors identified a highly significant burden of damaging de novo mutations in mutation-intolerant [probability of loss of function intolerance (pLI) >0.9] genes in nonsyndromic craniosynostosis probands (p = 5.9 × 10-6). Children with these mutations had a two-fold higher incidence of neurodevelopmental delay (p = 0.001) and a more than 20-fold greater incidence of intellectual disability (p = 7.2 × 10-7), and were 3.6-fold more likely to have delays that persisted past 5 years of age (p = 4.4 × 10-4) in comparison with children with nonsyndromic craniosynostosis without these mutations. Transmitted loss of function mutations in high-pLI genes also conferred a 1.9-fold greater risk of neurodevelopmental delay (p = 4.5 ×10-4).

CONCLUSIONS

These findings implicate genetic lesions concurrently impacting neurodevelopment and cranial morphogenesis in the pathoetiology of nonsyndromic craniosynostosis and identify a strong genetic influence on neurodevelopmental outcomes in affected children. These findings may eventually prove useful in determining which children with nonsyndromic craniosynostosis are most likely to benefit from surgical intervention.

CLINICAL QUESTION/LEVEL OF EVIDENCE

Risk, III.

Computational Genomics in the Era of Precision Medicine: Applications to Variant Analysis and Gene Therapy

Rapid methodological advances in statistical and computational genomics have enabled researchers to better identify and interpret both rare and common variants responsible for complex human diseases. As we continue to see an expansion of these advances in the field, it is now imperative for researchers to understand the resources and methodologies available for various data types and study designs. In this review, we provide an overview of recent methods for identifying rare and common variants and understanding their roles in disease etiology. Additionally, we discuss the strategy, challenge, and promise of gene therapy. As computational and statistical approaches continue to improve, we will have an opportunity to translate human genetic findings into personalized health care.

Schizophrenia, autism spectrum disorders and developmental disorders share specific disruptive coding mutations

People with schizophrenia are enriched for rare coding variants in genes associated with neurodevelopmental disorders, particularly autism spectrum disorders and intellectual disability. However, it is unclear if the same changes to gene function that increase risk to neurodevelopmental disorders also do so for schizophrenia. Using data from 3444 schizophrenia trios and 37,488 neurodevelopmental disorder trios, we show that within shared risk genes, de novo variants in schizophrenia and neurodevelopmental disorders are generally of the same functional category, and that specific de novo variants observed in neurodevelopmental disorders are enriched in schizophrenia (P = 5.0 × 10-6). The latter includes variants known to be pathogenic for syndromic disorders, suggesting that schizophrenia be included as a characteristic of those syndromes. Our findings imply that, in part, neurodevelopmental disorders and schizophrenia have shared molecular aetiology, and therefore likely overlapping pathophysiology, and support the hypothesis that at least some forms of schizophrenia lie on a continuum of neurodevelopmental disorders.

Meta-Analysis of Brain Gene Expression Data from Mouse Model Studies of Maternal Immune Activation Using Poly(I:C)

Maternal immune activation (MIA) is a known risk factor for schizophrenia (SCZ) and autism spectrum disorder (ASD) and is often modelled in animal studies in order to study the effect of prenatal infection on brain function including behaviour and gene expression. Although the effect of MIA on gene expression are highly heterogeneous, combining data from multiple gene expression studies in a robust method may shed light on the true underlying biological effects caused by MIA and this could inform studies of SCZ and ASD. This study combined four RNA-seq and microarray datasets in an overlap analysis and ranked meta-analysis in order to investigate genes, pathways and cell types dysregulated in the MIA mouse models. Genes linked to SCZ and ASD and crucial in neurodevelopmental processes including neural tube folding, regulation of cellular stress and neuronal/glial cell differentiation were among the most consistently dysregulated in these ranked analyses. Gene ontologies including K+ ion channel function, neuron and glial cell differentiation, synaptic structure, axonal outgrowth, cilia function and lipid metabolism were also strongly implicated. Single-cell analysis identified excitatory and inhibitory cell types in the cortex, hippocampus and striatum that may be affected by MIA and are also enriched for genes associated with SCZ, ASD and cognitive phenotypes. This points to the cellular location of molecular mechanisms that may be consistent between the MIA model and neurodevelopmental disease, improving our understanding of its utility to study prenatal infection as an environmental stressor.

A Multi-Omics Approach Using a Mouse Model of Cardiac Malformations for Prioritization of Human Congenital Heart Disease Contributing Genes

Congenital heart disease (CHD) is the most common type of birth defect, affecting ~1% of all live births. Malformations of the cardiac outflow tract (OFT) account for ~30% of all CHD and include a range of CHDs from bicuspid aortic valve (BAV) to tetralogy of Fallot (TOF). We hypothesized that transcriptomic profiling of a mouse model of CHD would highlight disease-contributing genes implicated in congenital cardiac malformations in humans. To test this hypothesis, we utilized global transcriptional profiling differences from a mouse model of OFT malformations to prioritize damaging, de novo variants identified from exome sequencing datasets from published cohorts of CHD patients. Notch1 +/- ; Nos3 -/- mice display a spectrum of cardiac OFT malformations ranging from BAV, semilunar valve (SLV) stenosis to TOF. Global transcriptional profiling of the E13.5 Notch1 +/- ; Nos3 -/- mutant mouse OFTs and wildtype controls was performed by RNA sequencing (RNA-Seq). Analysis of the RNA-Seq dataset demonstrated genes belonging to the Hif1α, Tgf-β, Hippo, and Wnt signaling pathways were differentially expressed in the mutant OFT. Mouse to human comparative analysis was then performed to determine if patients with TOF and SLV stenosis display an increased burden of damaging, genetic variants in gene homologs that were dysregulated in Notch1 +/- ; Nos3 -/- OFT. We found an enrichment of de novo variants in the TOF population among the 1,352 significantly differentially expressed genes in Notch1 +/- ; Nos3 -/- mouse OFT but not the SLV population. This association was not significant when comparing only highly expressed genes in the murine OFT to de novo variants in the TOF population. These results suggest that transcriptomic datasets generated from the appropriate temporal, anatomic and cellular tissues from murine models of CHD may provide a novel approach for the prioritization of disease-contributing genes in patients with CHD.

De novo variation in bipolar disorder

Bipolar disorder (BD) is a common, highly heritable disorder that affects 1-2% of the world's population. To date, most genetic studies of BD have focused on common gene variation, and while robustly associated loci have been identified, a substantial proportion of the heritability remains missing and could be partially attributable to rare variation. In this study, we apply a de novo paradigm in BD to identify newly arisen variants that have yet to undergo natural selection and may represent highly pathogenic variants. We performed whole genome sequencing of 97 trios of Ashkenazi Jewish descent, selecting "simplex" families with no family history of BD and an early age of onset. We found a total of 6882 de novo variants (an average of 70.9 ± 12.9 S.D. variants per trio), including 107 variants within protein-coding genes. We combined our exonic variations with the results of 79 previously published BD trios, identifying 20 loss-of-function (LoF) and 77 missense damaging de novo variants in BD. These variants showed significant enrichment for constrained genes and for genes located to the postsynaptic density (PSD) (all Bonferroni corrected p < 0.05). Pathway analyses showed enrichment in several pathways, including "Phosphoinositides (PI) and their downstream targets" (Bonferroni p = 4.2 × 10-6), a pathway prominently featured in lithium's hypothesized mechanism of action. In addition, while we found overall evidence for transmission of common variant polygenic risk of BD in our full sample (pTDT p = 2.21 × 10-4), specific trios with LoF variants showed no evidence of polygenic transmission. In sum, our findings support the de novo paradigm as a contributor to the genetic architecture of BD and provide evidence that constrained genes, as well as genes within the PSD and PI pathway harbor rare variation associated with BD.

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

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.

Genetic association of FMRP targets with psychiatric disorders

Genes encoding the mRNA targets of fragile X mental retardation protein (FMRP) are enriched for genetic association with psychiatric disorders. However, many FMRP targets possess functions that are themselves genetically associated with psychiatric disorders, including synaptic transmission and plasticity, making it unclear whether the genetic risk is truly related to binding by FMRP or is alternatively mediated by the sampling of genes better characterised by another trait or functional annotation. Using published common variant, rare coding variant and copy number variant data, we examined the relationship between FMRP binding and genetic association with schizophrenia, major depressive disorder and bipolar disorder. High-confidence targets of FMRP, derived from studies of multiple tissue types, were enriched for common schizophrenia risk alleles, as well as rare loss-of-function and de novo nonsynonymous variants in schizophrenia cases. Similarly, through common variation, FMRP targets were associated with major depressive disorder, and we present novel evidence of association with bipolar disorder. These relationships could not be explained by other functional annotations known to be associated with psychiatric disorders, including those related to synaptic structure and function. This study reinforces the evidence that targeting by FMRP captures a subpopulation of genes enriched for genetic association with a range of psychiatric disorders.

Essential genes from genome-wide screenings as a resource for neuropsychiatric disorders gene discovery

Genome-wide screenings of "essential genes", i.e., genes required for an organism or cell survival, have been traditionally conducted in vitro in cancer cell lines, limiting the translation of results to other tissues and non-cancerous cells. Recently, an in vivo screening was conducted in adult mouse striatum tissue, providing the first genome-wide dataset of essential genes in neuronal cells. Here, we aim to investigate the role of essential genes in brain development and disease risk with a comprehensive set of bioinformatics tools, including integration with transcriptomic data from developing human brain, publicly available data from genome-wide association studies, de novo mutation datasets for different neuropsychiatric disorders, and case-control transcriptomic data from postmortem brain tissues. For the first time, we found that the expression of neuronal essential genes (NEGs) increases before birth during the early development of human brain and maintains a relatively high expression after birth. On the contrary, common essential genes from cancer cell line screenings (ACEGs) tend to be expressed at high levels during development but quickly drop after birth. Both gene sets were enriched in neurodevelopmental disorders, but only NEGs were robustly associated with neuropsychiatric disorders risk genes. Finally, NEGs were more likely to show differential expression in the brains of neuropsychiatric disorders patients than ACEGs. Overall, genome-wide central nervous system screening of essential genes can provide new insights into neuropsychiatric diseases.

Analysis workflow to assess de novo genetic variants from human whole-exome sequencing

Here, we present a protocol to analyze de novo genetic variants derived from the whole-exome sequencing (WES) of proband-parent trios. We provide stepwise instructions for using existing pipelines to call de novo mutations (DNMs) and determine whether the observed number of such mutations is enriched relative to the expected number. This protocol may be extended to any human disease trio-based cohort. Cohort size is a limiting determinant to the discovery of high-confidence pathogenic DNMs.

For complete details on the use and execution of this protocol, please refer to Dong et al. (2020).

•

We demonstrate the ability to call de novo mutations from whole-exome sequencing data

•

This protocol is applied to WES from cohorts composed of proband and both parents

•

We demonstrate how to perform enrichment analysis using denovolyzR

•

The size of the trio-based cohort is a limiting determinant of this protocol’s accuracy

SATB2-LEMD2 interaction links nuclear shape plasticity to regulation of cognition-related genes

SATB2 is a schizophrenia risk gene and is genetically associated with human intelligence. How it affects cognition at molecular level is currently unknown. Here, we show that interactions between SATB2, a chromosomal scaffolding protein, and the inner nuclear membrane protein LEMD2 orchestrate the response of pyramidal neurons to neuronal activation. Exposure to novel environment in vivo causes changes in nuclear shape of CA1 hippocampal neurons via a SATB2-dependent mechanism. The activity-driven plasticity of the nuclear envelope requires not only SATB2, but also its protein interactor LEMD2 and the ESCRT-III/VPS4 membrane-remodeling complex. Furthermore, LEMD2 depletion in cortical neurons, similar to SATB2 ablation, affects neuronal activity-dependent regulation of multiple rapid and delayed primary response genes. In human genetic data, LEMD2-regulated genes are enriched for de novo mutations reported in intellectual disability and schizophrenia and are, like SATB2-regulated genes, enriched for common variants associated with schizophrenia and cognitive function. Hence, interactions between SATB2 and the inner nuclear membrane protein LEMD2 influence gene expression programs in pyramidal neurons that are linked to cognitive ability and psychiatric disorder etiology.

Semantic Similarity Analysis Reveals Robust Gene-Disease Relationships in Developmental and Epileptic Encephalopathies

More than 100 genetic etiologies have been identified in developmental and epileptic encephalopathies (DEEs), but correlating genetic findings with clinical features at scale has remained a hurdle because of a lack of frameworks for analyzing heterogenous clinical data. Here, we analyzed 31,742 Human Phenotype Ontology (HPO) terms in 846 individuals with existing whole-exome trio data and assessed associated clinical features and phenotypic relatedness by using HPO-based semantic similarity analysis for individuals with de novo variants in the same gene. Gene-specific phenotypic signatures included associations of SCN1A with “complex febrile seizures” (HP: 0011172; p = 2.1 × 10⁻⁵) and “focal clonic seizures” (HP: 0002266; p = 8.9 × 10⁻⁶), STXBP1 with “absent speech” (HP: 0001344; p = 1.3 × 10⁻¹¹), and SLC6A1 with “EEG with generalized slow activity” (HP: 0010845; p = 0.018). Of 41 genes with de novo variants in two or more individuals, 11 genes showed significant phenotypic similarity, including SCN1A (n = 16, p < 0.0001), STXBP1 (n = 14, p = 0.0021), and KCNB1 (n = 6, p = 0.011). Including genetic and phenotypic data of control subjects increased phenotypic similarity for all genetic etiologies, whereas the probability of observing de novo variants decreased, emphasizing the conceptual differences between semantic similarity analysis and approaches based on the expected number of de novo events. We demonstrate that HPO-based phenotype analysis captures unique profiles for distinct genetic etiologies, reflecting the breadth of the phenotypic spectrum in genetic epilepsies. Semantic similarity can be used to generate statistical evidence for disease causation analogous to the traditional approach of primarily defining disease entities through similar clinical features.

Genes and Pathways Implicated in Tetralogy of Fallot Revealed by Ultra-Rare Variant Burden Analysis in 231 Genome Sequences

Recent genome-wide studies of rare genetic variants have begun to implicate novel mechanisms for tetralogy of Fallot (TOF), a severe congenital heart defect (CHD). To provide statistical support for case-only data without parental genomes, we re-analyzed genome sequences of 231 individuals with TOF (n = 175) or related CHD. We adapted a burden test originally developed for de novo variants to assess ultra-rare variant burden in individual genes, and in gene-sets corresponding to functional pathways and mouse phenotypes, accounting for highly correlated gene-sets and for multiple testing. For truncating variants, the gene burden test confirmed significant burden in FLT4 (Bonferroni corrected p-value < 0.01). For missense variants, burden in NOTCH1 achieved genome-wide significance only when restricted to constrained genes (i.e., under negative selection, Bonferroni corrected p-value = 0.004), and showed enrichment for variants affecting the extracellular domain, especially those disrupting cysteine residues forming disulfide bonds (OR = 39.8 vs. gnomAD). Individuals with NOTCH1 ultra-rare missense variants, all with TOF, were enriched for positive family history of CHD. Other genes not previously implicated in CHD had more modest statistical support in gene burden tests. Gene-set burden tests for truncating variants identified a cluster of pathways corresponding to VEGF signaling (FDR = 0%), and of mouse phenotypes corresponding to abnormal vasculature (FDR = 0.8%); these suggested additional candidate genes not previously identified (e.g., WNT5A and ZFAND5). Results for the most promising genes were driven by the TOF subset of the cohort. The findings support the importance of ultra-rare variants disrupting genes involved in VEGF and NOTCH signaling in the genetic architecture of TOF, accounting for 11–14% of individuals in the TOF cohort. These proof-of-principle data indicate that this statistical methodology could assist in analyzing case-only sequencing data in which ultra-rare variants, whether de novo or inherited, contribute to the genetic etiopathogenesis of a complex disorder.

Genes regulated by BCL11B during T-cell development are enriched for de novo mutations found in schizophrenia patients

While abnormal neurodevelopment contributes to schizophrenia (SCZ) risk, there is also evidence to support a role for immune dysfunction in SCZ. BCL11B, associated with SCZ in genome-wide association study (GWAS), is a transcription factor that regulates the differentiation and development of cells in the central nervous and immune systems. Here, we use functional genomics data from studies of BCL11B to investigate the contribution of neuronal and immune processes to SCZ pathophysiology. We identified the gene targets of BCL11B in brain striatal cells (n = 223 genes), double negative 4 (DN4) developing T cells (n = 114 genes) and double positive (DP) developing T cells (n = 518 genes) using an integrated analysis of RNA-seq and ChIP-seq data. No gene-set was enriched for genes containing common variants associated with SCZ but the DP gene-set was enriched for genes containing missense de novo mutations (DNMs; p = .001) using data from 3,447 SCZ trios. Post hoc analysis revealed the enrichment to be stronger for DP genes negatively regulated by BCL11B. Biological processes enriched for genes negatively regulated by BCL11B in DP gene-set included immune system development and cytokine signaling. These analyses, leveraging a GWAS-identified SCZ risk gene and data on gene expression and transcription factor binding, indicate that DNMs in immune pathways contribute to SCZ risk.

SLC12A2 variants cause a neurodevelopmental disorder or cochleovestibular defect

The SLC12 gene family consists of SLC12A1-SLC12A9, encoding electroneutral cation-coupled chloride co-transporters. SCL12A2 has been shown to play a role in corticogenesis and therefore represents a strong candidate neurodevelopmental disorder gene. Through trio exome sequencing we identified de novo mutations in SLC12A2 in six children with neurodevelopmental disorders. All had developmental delay or intellectual disability ranging from mild to severe. Two had sensorineural deafness. We also identified SLC12A2 variants in three individuals with non-syndromic bilateral sensorineural hearing loss and vestibular areflexia. The SLC12A2 de novo mutation rate was demonstrated to be significantly elevated in the deciphering developmental disorders cohort. All tested variants were shown to reduce co-transporter function in Xenopus laevis oocytes. Analysis of SLC12A2 expression in foetal brain at 16-18 weeks post-conception revealed high expression in radial glial cells, compatible with a role in neurogenesis. Gene co-expression analysis in cells robustly expressing SLC12A2 at 16-18 weeks post-conception identified a transcriptomic programme associated with active neurogenesis. We identify SLC12A2 de novo mutations as the cause of a novel neurodevelopmental disorder and bilateral non-syndromic sensorineural hearing loss and provide further data supporting a role for this gene in human neurodevelopment.

Genome-wide Enrichment of De Novo Coding Mutations in Orofacial Cleft Trios

Although de novo mutations (DNMs) are known to increase an individual's risk of congenital defects, DNMs have not been fully explored regarding orofacial clefts (OFCs), one of the most common human birth defects. Therefore, whole-genome sequencing of 756 child-parent trios of European, Colombian, and Taiwanese ancestry was performed to determine the contributions of coding DNMs to an individual's OFC risk. Overall, we identified a significant excess of loss-of-function DNMs in genes highly expressed in craniofacial tissues, as well as genes associated with known autosomal dominant OFC syndromes. This analysis also revealed roles for zinc-finger homeobox domain and SOX2-interacting genes in OFC etiology.

Rare and de novo duplications containing SHOX in clubfoot

Introduction

Congenital clubfoot is a common birth defect that affects at least 0.1% of all births. Nearly 25% cases are familial and the remaining are sporadic in inheritance. Copy number variants (CNVs) involving transcriptional regulators of limb development, including PITX1 and TBX4, have previously been shown to cause familial clubfoot, but much of the heritability remains unexplained.

Methods

Exome sequence data from 816 unrelated clubfoot cases and 2645 in-house controls were analysed using coverage data to identify rare CNVs. The precise size and location of duplications were then determined using high-density Affymetrix Cytoscan chromosomal microarray (CMA). Segregation in families and de novo status were determined using qantitative PCR.

Results

Chromosome Xp22.33 duplications involving SHOX were identified in 1.1% of cases (9/816) compared with 0.07% of in-house controls (2/2645) (p=7.98×10⁻⁵, OR=14.57) and 0.27% (38/13592) of Atherosclerosis Risk in Communities/the Wellcome Trust Case Control Consortium 2 controls (p=0.001, OR=3.97). CMA validation confirmed an overlapping 180.28 kb duplicated region that included SHOX exons as well as downstream non-coding regions. In four of six sporadic cases where DNA was available for unaffected parents, the duplication was de novo. The probability of four de novo mutations in SHOX by chance in a cohort of 450 sporadic clubfoot cases is 5.4×10^–10.

Conclusions

Microduplications of the pseudoautosomal chromosome Xp22.33 region (PAR1) containing SHOX and downstream enhancer elements occur in ~1% of patients with clubfoot. SHOX and regulatory regions have previously been implicated in skeletal dysplasia as well as idiopathic short stature, but have not yet been reported in clubfoot. SHOX duplications likely contribute to clubfoot pathogenesis by altering early limb development.

BAZ2B haploinsufficiency as a cause of developmental delay, intellectual disability, and autism spectrum disorder

The bromodomain adjacent to zinc finger 2B gene (BAZ2B) encodes a protein involved in chromatin remodeling. Loss of BAZ2B function has been postulated to cause neurodevelopmental disorders. To determine whether BAZ2B deficiency is likely to contribute to the pathogenesis of these disorders, we performed bioinformatics analyses that demonstrated a high level of functional convergence during fetal cortical development between BAZ2B and genes known to cause autism spectrum disorder (ASD) and neurodevelopmental disorder. We also found an excess of de novo BAZ2B loss-of-function variants in exome sequencing data from previously published cohorts of individuals with neurodevelopmental disorders. We subsequently identified seven additional individuals with heterozygous deletions, stop-gain, or de novo missense variants affecting BAZ2B. All of these individuals have developmental delay (DD), intellectual disability (ID), and/or ASD. Taken together, our findings suggest that haploinsufficiency of BAZ2B causes a neurodevelopmental disorder, whose cardinal features include DD, ID, and ASD.

Synaptic and Gene Regulatory Mechanisms in Schizophrenia, Autism, and 22q11.2 Copy Number Variant-Mediated Risk for Neuropsychiatric Disorders

BACKGROUND

22q11.2 copy number variants are among the most highly penetrant genetic risk variants for developmental neuropsychiatric disorders such as schizophrenia (SCZ) and autism spectrum disorder (ASD). However, the specific mechanisms through which they confer risk remain unclear.

METHODS

Using a functional genomics approach, we integrated transcriptomic data from the developing human brain, genome-wide association findings for SCZ and ASD, protein interaction data, and gene expression signatures from SCZ and ASD postmortem cortex to 1) organize genes into the developmental cellular and molecular systems within which they operate, 2) identify neurodevelopmental processes associated with polygenic risk for SCZ and ASD across the allelic frequency spectrum, and 3) elucidate pathways and individual genes through which 22q11.2 copy number variants may confer risk for each disorder.

RESULTS

Polygenic risk for SCZ and ASD converged on partially overlapping neurodevelopmental modules involved in synaptic function and transcriptional regulation, with ASD risk variants additionally enriched for modules involved in neuronal differentiation during fetal development. The 22q11.2 locus formed a large protein network during development that disproportionately affected SCZ-associated and ASD-associated neurodevelopmental modules, including loading highly onto synaptic and gene regulatory pathways. SEPT5, PI4KA, and SNAP29 genes are candidate drivers of 22q11.2 synaptic pathology relevant to SCZ and ASD, and DGCR8 and HIRA are candidate drivers of disease-relevant alterations in gene regulation.

CONCLUSIONS

This approach offers a powerful framework to identify neurodevelopmental processes affected by diverse risk variants for SCZ and ASD and elucidate mechanisms through which highly penetrant, multigene copy number variants contribute to disease risk.

Sex-Based Analysis of De Novo Variants in Neurodevelopmental Disorders

While genes with an excess of de novo mutations (DNMs) have been identified in children with neurodevelopmental disorders (NDDs), few studies focus on DNM patterns where the sex of affected children is examined separately. We considered ∼8,825 sequenced parent-child trios (n ∼26,475 individuals) and identify 54 genes with a DNM enrichment in males (n = 18), females (n = 17), or overlapping in both the male and female subsets (n = 19). A replication cohort of 18,778 sequenced parent-child trios (n = 56,334 individuals) confirms 25 genes (n = 3 in males, n = 7 in females, n = 15 in both male and female subsets). As expected, we observe significant enrichment on the X chromosome for females but also find autosomal genes with potential sex bias (females, CDK13, ITPR1; males, CHD8, MBD5, SYNGAP1); 6.5% of females harbor a DNM in a female-enriched gene, whereas 2.7% of males have a DNM in a male-enriched gene. Sex-biased genes are enriched in transcriptional processes and chromatin binding, primarily reside in the nucleus of cells, and have brain expression. By downsampling, we find that DNM gene discovery is greatest when studying affected females. Finally, directly comparing de novo allele counts in NDD-affected males and females identifies one replicated genome-wide significant gene (DDX3X) with locus-specific enrichment in females. Our sex-based DNM enrichment analysis identifies candidate NDD genes differentially affecting males and females and indicates that the study of females with NDDs leads to greater gene discovery consistent with the female-protective effect.

Convergent Evidence That ZNF804A Is a Regulator of Pre-messenger RNA Processing and Gene Expression

Genome-wide association studies have linked common variation in ZNF804A with an increased risk of schizophrenia. However, little is known about the biology of ZNF804A and its role in schizophrenia. Here, we investigate the function of ZNF804A using a variety of complementary molecular techniques. We show that ZNF804A is a nuclear protein that interacts with neuronal RNA splicing factors and RNA-binding proteins including RBFOX1, which is also associated with schizophrenia, CELF3/4, components of the ubiquitin-proteasome system and the ZNF804A paralog, GPATCH8. GPATCH8 also interacts with splicing factors and is localized to nuclear speckles indicative of a role in pre-messenger RNA (mRNA) processing. Sequence analysis showed that GPATCH8 contains ultraconserved, alternatively spliced poison exons that are also regulated by RBFOX proteins. ZNF804A knockdown in SH-SY5Y cells resulted in robust changes in gene expression and pre-mRNA splicing converging on pathways associated with nervous system development, synaptic contact, and cell adhesion. We observed enrichment (P = 1.66 × 10-9) for differentially spliced genes in ZNF804A-depleted cells among genes that contain RBFOX-dependent alternatively spliced exons. Differentially spliced genes in ZNF804A-depleted cells were also enriched for genes harboring de novo loss of function mutations in autism spectrum disorder (P = 6.25 × 10-7, enrichment 2.16) and common variant alleles associated with schizophrenia (P = .014), bipolar disorder and schizophrenia (P = .003), and autism spectrum disorder (P = .005). These data suggest that ZNF804A and its paralogs may interact with neuronal-splicing factors and RNA-binding proteins to regulate the expression of a subset of synaptic and neurodevelopmental genes.

Disruptive mutations in TANC2 define a neurodevelopmental syndrome associated with psychiatric disorders

Postsynaptic density (PSD) proteins have been implicated in the pathophysiology of neurodevelopmental and psychiatric disorders. Here, we present detailed clinical and genetic data for 20 patients with likely gene-disrupting mutations in TANC2-whose protein product interacts with multiple PSD proteins. Pediatric patients with disruptive mutations present with autism, intellectual disability, and delayed language and motor development. In addition to a variable degree of epilepsy and facial dysmorphism, we observe a pattern of more complex psychiatric dysfunction or behavioral problems in adult probands or carrier parents. Although this observation requires replication to establish statistical significance, it also suggests that mutations in this gene are associated with a variety of neuropsychiatric disorders consistent with its postsynaptic function. We find that TANC2 is expressed broadly in the human developing brain, especially in excitatory neurons and glial cells, but shows a more restricted pattern in Drosophila glial cells where its disruption affects behavioral outcomes.

Mutations in TFAP2B and previously unimplicated genes of the BMP, Wnt, and Hedgehog pathways in syndromic craniosynostosis

Craniosynostosis (CS) is a frequent congenital malformation featuring premature fusion of cranial sutures; 15% of these children have syndromic disease, often due to rare mutations with large effect. While many genes causing Mendelian forms of syndromic CS have been identified, clinical sequencing often fails to identify a likely causative mutation. We performed whole-exome sequencing of 12 case-parent trios with previously negative genetic evaluations. The results identified likely pathogenic mutations in TFAP2B, KAT6A, GLI2, SOX11, CTNNA1, and GPC4 in these families, adding several loci to those known to cause syndromic CS. The findings have implications for determining risk of disease in subsequent offspring and demonstrate that unexplained syndromic CS cases are a particularly rich vein for discovery of CS loci.

Craniosynostosis (CS) is a frequent congenital anomaly featuring the premature fusion of 1 or more sutures of the cranial vault. Syndromic cases, featuring additional congenital anomalies, make up 15% of CS. While many genes underlying syndromic CS have been identified, the cause of many syndromic cases remains unknown. We performed exome sequencing of 12 syndromic CS cases and their parents, in whom previous genetic evaluations were unrevealing. Damaging de novo or transmitted loss of function (LOF) mutations were found in 8 genes that are highly intolerant to LOF mutation (P = 4.0 × 10⁻⁸); additionally, a rare damaging mutation in SOX11, which has a lower level of intolerance, was identified. Four probands had rare damaging mutations (2 de novo) in TFAP2B, a transcription factor that orchestrates neural crest cell migration and differentiation; this mutation burden is highly significant (P = 8.2 × 10⁻¹²). Three probands had rare damaging mutations in GLI2, SOX11, or GPC4, which function in the Hedgehog, BMP, and Wnt signaling pathways; other genes in these pathways have previously been implicated in syndromic CS. Similarly, damaging de novo mutations were identified in genes encoding the chromatin modifier KAT6A, and CTNNA1, encoding catenin α-1. These findings establish TFAP2B as a CS gene, have implications for assessing risk to subsequent children in these families, and provide evidence implicating other genes in syndromic CS. This high yield indicates the value of performing exome sequencing of syndromic CS patients when sequencing of known disease loci is unrevealing.

De novo mutations of the ATP6V1A gene cause developmental encephalopathy with epilepsy

V-type proton (H⁺) ATPase (v-ATPase) is a multi-subunit proton pump that regulates pH homeostasis in all eukaryotic cells; in neurons, v-ATPase plays additional and unique roles in synapse function. Through whole exome sequencing, we identified de novo heterozygous mutations (p.Pro27Arg, p.Asp100Tyr, p.Asp349Asn, p.Asp371Gly) in ATP6V1A, encoding the A subunit of v-ATPase, in four patients with developmental encephalopathy with epilepsy. Early manifestations, observed in all patients, were developmental delay and febrile seizures, evolving to encephalopathy with profound delay, hypotonic/dyskinetic quadriparesis and intractable multiple seizure types in two patients (p.Pro27Arg, p.Asp100Tyr), and to moderate delay with milder epilepsy in the other two (p.Asp349Asn, p.Asp371Gly). Modelling performed on the available prokaryotic and eukaryotic structures of v-ATPase predicted p.Pro27Arg to perturb subunit interaction, p.Asp100Tyr to cause steric hindrance and destabilize protein folding, p.Asp349Asn to affect the catalytic function and p.Asp371Gly to impair the rotation process, necessary for proton transport. We addressed the impact of p.Asp349Asn and p.Asp100Tyr mutations on ATP6V1A expression and function by analysing ATP6V1A-overexpressing HEK293T cells and patients’ lymphoblasts. The p.Asp100Tyr mutant was characterized by reduced expression due to increased degradation. Conversely, no decrease in expression and clearance was observed for p.Asp349Asn. In HEK293T cells overexpressing either pathogenic or control variants, p.Asp349Asn significantly increased LysoTracker® fluorescence with no effects on EEA1 and LAMP1 expression. Conversely, p.Asp100Tyr decreased both LysoTracker® fluorescence and LAMP1 levels, leaving EEA1 expression unaffected. Both mutations decreased v-ATPase recruitment to autophagosomes, with no major impact on autophagy. Experiments performed on patients’ lymphoblasts using the LysoSensor™ probe revealed lower pH of endocytic organelles for p.Asp349Asn and a reduced expression of LAMP1 with no effect on the pH for p.Asp100Tyr. These data demonstrate gain of function for p.Asp349Asn characterized by an increased proton pumping in intracellular organelles, and loss of function for p.Asp100Tyr with decreased expression of ATP6V1A and reduced levels of lysosomal markers. We expressed p.Asp349Asn and p.Asp100Tyr in rat hippocampal neurons and confirmed significant and opposite effects in lysosomal labelling. However, both mutations caused a similar defect in neurite elongation accompanied by loss of excitatory inputs, revealing that altered lysosomal homeostasis markedly affects neurite development and synaptic connectivity. This study provides evidence that de novo heterozygous ATP6V1A mutations cause a developmental encephalopathy with a pathomechanism that involves perturbations of lysosomal homeostasis and neuronal connectivity, uncovering a novel role for v-ATPase in neuronal development.

Missense Variants in the Histone Acetyltransferase Complex Component Gene TRRAP Cause Autism and Syndromic Intellectual Disability

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.

De Novo SOX4 Variants Cause a Neurodevelopmental Disease Associated with Mild Dysmorphism

SOX4, together with SOX11 and SOX12, forms group C of SRY-related (SOX) transcription factors. They play key roles, often in redundancy, in multiple developmental pathways, including neurogenesis and skeletogenesis. De novo SOX11 heterozygous mutations have been shown to cause intellectual disability, growth deficiency, and dysmorphic features compatible with mild Coffin-Siris syndrome. Using trio-based exome sequencing, we here identify de novo SOX4 heterozygous missense variants in four children who share developmental delay, intellectual disability, and mild facial and digital morphological abnormalities. SOX4 is highly expressed in areas of active neurogenesis in human fetuses, and sox4 knockdown in Xenopus embryos diminishes brain and whole-body size. The SOX4 variants cluster in the highly conserved, SOX family-specific HMG domain, but each alters a different residue. In silico tools predict that each variant affects a distinct structural feature of this DNA-binding domain, and functional assays demonstrate that these SOX4 proteins carrying these variants are unable to bind DNA in vitro and transactivate SOX reporter genes in cultured cells. These variants are not found in the gnomAD database of individuals with presumably normal development, but 12 other SOX4 HMG-domain missense variants are recorded and all demonstrate partial to full activity in the reporter assay. Taken together, these findings point to specific SOX4 HMG-domain missense variants as the cause of a characteristic human neurodevelopmental disorder associated with mild facial and digital dysmorphism.

Modeling mutant/wild-type interactions to ascertain pathogenicity of PROKR2 missense variants in patients with isolated GnRH deficiency

A major challenge in human genetics is the validation of pathogenicity of heterozygous missense variants. This problem is well-illustrated by PROKR2 variants associated with Isolated GnRH Deficiency (IGD). Homozygous, loss of function variants in PROKR2 was initially implicated in autosomal recessive IGD; however, most IGD-associated PROKR2 variants are heterozygous. Moreover, while IGD patient cohorts are enriched for PROKR2 missense variants similar rare variants are also found in normal individuals. To elucidate the pathogenic mechanisms distinguishing IGD-associated PROKR2 variants from rare variants in controls, we assessed 59 variants using three approaches: (i) in silico prediction, (ii) traditional in vitro functional assays across three signaling pathways with mutant-alone transfections, and (iii) modified in vitro assays with mutant and wild-type expression constructs co-transfected to model in vivo heterozygosity. We found that neither in silico analyses nor traditional in vitro assessments of mutants transfected alone could distinguish IGD variants from control variants. However, in vitro co-transfections revealed that 15/34 IGD variants caused loss-of-function (LoF), including 3 novel dominant-negatives, while only 4/25 control variants caused LoF. Surprisingly, 19 IGD-associated variants were benign or exhibited LoF that could be rescued by WT co-transfection. Overall, variants that were LoF in ≥ 2 signaling assays under co-transfection conditions were more likely to be disease-associated than benign or 'rescuable' variants. Our findings suggest that in vitro modeling of WT/Mutant interactions increases the resolution for identifying causal variants, uncovers novel dominant negative mutations, and provides new insights into the pathogenic mechanisms underlying heterozygous PROKR2 variants.

De novo variants in congenital diaphragmatic hernia identify MYRF as a new syndrome and reveal genetic overlaps with other developmental disorders

Congenital diaphragmatic hernia (CDH) is a severe birth defect that is often accompanied by other congenital anomalies. Previous exome sequencing studies for CDH have supported a role of de novo damaging variants but did not identify any recurrently mutated genes. To investigate further the genetics of CDH, we analyzed de novo coding variants in 362 proband-parent trios including 271 new trios reported in this study. We identified four unrelated individuals with damaging de novo variants in MYRF (P = 5.3x10(-8)), including one likely gene-disrupting (LGD) and three deleterious missense (D-mis) variants. Eight additional individuals with de novo LGD or missense variants were identified from our other genetic studies or from the literature. Common phenotypes of MYRF de novo variant carriers include CDH, congenital heart disease and genitourinary abnormalities, suggesting that it represents a novel syndrome. MYRF is a membrane associated transcriptional factor highly expressed in developing diaphragm and is depleted of LGD variants in the general population. All de novo missense variants aggregated in two functional protein domains. Analyzing the transcriptome of patient-derived diaphragm fibroblast cells suggest that disease associated variants abolish the transcription factor activity. Furthermore, we showed that the remaining genes with damaging variants in CDH significantly overlap with genes implicated in other developmental disorders. Gene expression patterns and patient phenotypes support pleiotropic effects of damaging variants in these genes on CDH and other developmental disorders. Finally, functional enrichment analysis implicates the disruption of regulation of gene expression, kinase activities, intra-cellular signaling, and cytoskeleton organization as pathogenic mechanisms in CDH.

De Novo Mutation in Genes Regulating Neural Stem Cell Fate in Human Congenital Hydrocephalus

Congenital hydrocephalus (CH), featuring markedly enlarged brain ventricles, is thought to arise from failed cerebrospinal fluid (CSF) homeostasis and is treated with lifelong surgical CSF shunting with substantial morbidity. CH pathogenesis is poorly understood. Exome sequencing of 125 CH trios and 52 additional probands identified three genes with significant burden of rare damaging de novo or transmitted mutations: TRIM71 (p = 2.15 × 10-7), SMARCC1 (p = 8.15 × 10-10), and PTCH1 (p = 1.06 × 10-6). Additionally, two de novo duplications were identified at the SHH locus, encoding the PTCH1 ligand (p = 1.2 × 10-4). Together, these probands account for ∼10% of studied cases. Strikingly, all four genes are required for neural tube development and regulate ventricular zone neural stem cell fate. These results implicate impaired neurogenesis (rather than active CSF accumulation) in the pathogenesis of a subset of CH patients, with potential diagnostic, prognostic, and therapeutic ramifications.

Rare variants in SOX17 are associated with pulmonary arterial hypertension with congenital heart disease

BACKGROUND

Pulmonary arterial hypertension (PAH) is a rare disease characterized by distinctive changes in pulmonary arterioles that lead to progressive pulmonary arterial pressures, right-sided heart failure, and a high mortality rate. Up to 30% of adult and 75% of pediatric PAH cases are associated with congenital heart disease (PAH-CHD), and the underlying etiology is largely unknown. There are no known major risk genes for PAH-CHD.

METHODS

To identify novel genetic causes of PAH-CHD, we performed whole exome sequencing in 256 PAH-CHD patients. We performed a case-control gene-based association test of rare deleterious variants using 7509 gnomAD whole genome sequencing population controls. We then screened a separate cohort of 413 idiopathic and familial PAH patients without CHD for rare deleterious variants in the top association gene.

RESULTS

We identified SOX17 as a novel candidate risk gene (p = 5.5e-7). SOX17 is highly constrained and encodes a transcription factor involved in Wnt/β-catenin and Notch signaling during development. We estimate that rare deleterious variants contribute to approximately 3.2% of PAH-CHD cases. The coding variants identified include likely gene-disrupting (LGD) and deleterious missense, with most of the missense variants occurring in a highly conserved HMG-box protein domain. We further observed an enrichment of rare deleterious variants in putative targets of SOX17, many of which are highly expressed in developing heart and pulmonary vasculature. In the cohort of PAH without CHD, rare deleterious variants of SOX17 were observed in 0.7% of cases.

CONCLUSIONS

These data strongly implicate SOX17 as a new risk gene contributing to PAH-CHD as well as idiopathic/familial PAH. Replication in other PAH cohorts and further characterization of the clinical phenotype will be important to confirm the precise role of SOX17 and better estimate the contribution of genes regulated by SOX17.

Exome Sequencing in Children With Pulmonary Arterial Hypertension Demonstrates Differences Compared With Adults

BACKGROUND

Pulmonary arterial hypertension (PAH) is a rare disease characterized by pulmonary arteriole remodeling, elevated arterial pressure and resistance, and subsequent heart failure. Compared with adult-onset disease, pediatric-onset PAH is more heterogeneous and often associated with worse prognosis. Although BMPR2 mutations underlie ≈70% of adult familial PAH (FPAH) cases, the genetic basis of PAH in children is less understood.

METHODS

We performed genetic analysis of 155 pediatric- and 257 adult-onset PAH patients, including both FPAH and sporadic, idiopathic PAH (IPAH). After screening for 2 common PAH risk genes, mutation-negative FPAH and all IPAH cases were evaluated by exome sequencing.

RESULTS

We observed similar frequencies of rare, deleterious BMPR2 mutations in pediatric- and adult-onset patients: ≈55% in FPAH and 10% in IPAH patients in both age groups. However, there was significant enrichment of TBX4 mutations in pediatric- compared with adult-onset patients (IPAH: 10/130 pediatric versus 0/178 adult-onset), and TBX4 carriers had younger mean age-of-onset compared with BMPR2 carriers. Mutations in other known PAH risk genes were infrequent in both age groups. Notably, among pediatric IPAH patients without mutations in known risk genes, exome sequencing revealed a 2-fold enrichment of de novo likely gene-damaging and predicted deleterious missense variants.

CONCLUSIONS

Mutations in known PAH risk genes accounted for ≈70% to 80% of FPAH in both age groups, 21% of pediatric-onset IPAH, and 11% of adult-onset IPAH. Rare, predicted deleterious variants in TBX4 are enriched in pediatric patients and de novo variants in novel genes may explain ≈19% of pediatric-onset IPAH cases.

SPARK: A US Cohort of 50,000 Families to Accelerate Autism Research

The Simons Foundation Autism Research Initiative (SFARI) has launched SPARKForAutism.org, a dynamic platform that is engaging thousands of individuals with autism spectrum disorder (ASD) and connecting them to researchers. By making all data accessible, SPARK seeks to increase our understanding of ASD and accelerate new supports and treatments for ASD.

De Novo Mutations in Protein Kinase Genes CAMK2A and CAMK2B Cause Intellectual Disability

Calcium/calmodulin-dependent protein kinase II (CAMK2) is one of the first proteins shown to be essential for normal learning and synaptic plasticity in mice, but its requirement for human brain development has not yet been established. Through a multi-center collaborative study based on a whole-exome sequencing approach, we identified 19 exceedingly rare de novo CAMK2A or CAMK2B variants in 24 unrelated individuals with intellectual disability. Variants were assessed for their effect on CAMK2 function and on neuronal migration. For both CAMK2A and CAMK2B, we identified mutations that decreased or increased CAMK2 auto-phosphorylation at Thr286/Thr287. We further found that all mutations affecting auto-phosphorylation also affected neuronal migration, highlighting the importance of tightly regulated CAMK2 auto-phosphorylation in neuronal function and neurodevelopment. Our data establish the importance of CAMK2A and CAMK2B and their auto-phosphorylation in human brain function and expand the phenotypic spectrum of the disorders caused by variants in key players of the glutamatergic signaling pathway.

High Rate of Recurrent De Novo Mutations in Developmental and Epileptic Encephalopathies

Developmental and epileptic encephalopathy (DEE) is a group of conditions characterized by the co-occurrence of epilepsy and intellectual disability (ID), typically with developmental plateauing or regression associated with frequent epileptiform activity. The cause of DEE remains unknown in the majority of cases. We performed whole-genome sequencing (WGS) in 197 individuals with unexplained DEE and pharmaco-resistant seizures and in their unaffected parents. We focused our attention on de novo mutations (DNMs) and identified candidate genes containing such variants. We sought to identify additional subjects with DNMs in these genes by performing targeted sequencing in another series of individuals with DEE and by mining various sequencing datasets. We also performed meta-analyses to document enrichment of DNMs in candidate genes by leveraging our WGS dataset with those of several DEE and ID series. By combining these strategies, we were able to provide a causal link between DEE and the following genes: NTRK2, GABRB2, CLTC, DHDDS, NUS1, RAB11A, GABBR2, and SNAP25. Overall, we established a molecular diagnosis in 63/197 (32%) individuals in our WGS series. The main cause of DEE in these individuals was de novo point mutations (53/63 solved cases), followed by inherited mutations (6/63 solved cases) and de novo CNVs (4/63 solved cases). De novo missense variants explained a larger proportion of individuals in our series than in other series that were primarily ascertained because of ID. Moreover, these DNMs were more frequently recurrent than those identified in ID series. These observations indicate that the genetic landscape of DEE might be different from that of ID without epilepsy.

Contribution of rare inherited and de novo variants in 2,871 congenital heart disease probands

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.

De novo mutations in inhibitors of Wnt, BMP, and Ras/ERK signaling pathways in non-syndromic midline craniosynostosis

Non-syndromic craniosynostosis (NSC) is a frequent congenital malformation in which one or more cranial sutures fuse prematurely. Mutations causing rare syndromic craniosynostoses in humans and engineered mouse models commonly increase signaling of the Wnt, bone morphogenetic protein (BMP), or Ras/ERK pathways, converging on shared nuclear targets that promote bone formation. In contrast, the genetics of NSC is largely unexplored. More than 95% of NSC is sporadic, suggesting a role for de novo mutations. Exome sequencing of 291 parent-offspring trios with midline NSC revealed 15 probands with heterozygous damaging de novo mutations in 12 negative regulators of Wnt, BMP, and Ras/ERK signaling (10.9-fold enrichment, P = 2.4 × 10^-11). SMAD6 had 4 de novo and 14 transmitted mutations; no other gene had more than 1. Four familial NSC kindreds had mutations in genes previously implicated in syndromic disease. Collectively, these mutations contribute to 10% of probands. Mutations are predominantly loss-of-function, implicating haploinsufficiency as a frequent mechanism. A common risk variant near BMP2 increased the penetrance of SMAD6 mutations and was overtransmitted to patients with de novo mutations in other genes in these pathways, supporting a frequent two-locus pathogenesis. These findings implicate new genes in NSC and demonstrate related pathophysiology of common non-syndromic and rare syndromic craniosynostoses. These findings have implications for diagnosis, risk of recurrence, and risk of adverse neurodevelopmental outcomes. Finally, the use of pathways identified in rare syndromic disease to find genes accounting for non-syndromic cases may prove broadly relevant to understanding other congenital disorders featuring high locus heterogeneity.

Genome-wide enrichment of damaging de novo variants in patients with isolated and complex congenital diaphragmatic hernia

Congenital Diaphragmatic Hernia (CDH) is a common and often lethal birth defect characterized by diaphragmatic structural defects and pulmonary hypoplasia. CDH is isolated in 60% of newborns, but may also be part of a complex phenotype with additional anomalies. We performed whole exome sequencing (WES) on 87 individuals with isolated or complex CDH and on their unaffected parents, to assess the contribution of de novo mutations in the etiology of diaphragmatic and pulmonary defects and to identify new candidate genes. A combined analysis with 39 additional trios with complex CDH, previously published, revealed a significant genome-wide burden of de novo variants compared to background mutation rate and 900 control trios. We identified an increased burden of likely gene-disrupting (LGD, i.e. nonsense, frameshift, and canonical splice site) and predicted deleterious missense (D-mis) variants in complex and isolated CDH patients. Overall, an excess of predicted damaging de novo LGD and D-mis variants relative to the expected frequency contributed to 21% of complex cases and 12% of isolated CDH cases. The burden of de novo variants was higher in genes expressed in the developing mouse diaphragm and heart. Some overlap with genes responsible for congenital heart defects and neurodevelopmental disorders was observed in CDH patients within our cohorts. We propose that de novo variants contribute significantly to the development of CDH.