Neurodevelopmental disease genes implicated by de novo mutation and copy number variation morbidity

NASP gene contributes to autism by epigenetic dysregulation of neural and immune pathways

BACKGROUND

Epigenetics makes substantial contribution to the aetiology of autism spectrum disorder (ASD) and may harbour a unique opportunity to prevent the development of ASD. We aimed to identify novel epigenetic genes involved in ASD aetiology.

METHODS

Trio-based whole exome sequencing was conducted on ASD families. Genome editing technique was used to knock out the candidate causal gene in a relevant cell line. ATAC-seq, ChIP-seq and RNA-seq were performed to investigate the functional impact of knockout (KO) or mutation in the candidate gene.

RESULTS

We identified a novel candidate gene NASP (nuclear autoantigenic sperm protein) for epigenetic dysregulation in ASD in a Chinese nuclear family including one proband with autism and comorbid atopic disease. The de novo likely gene disruptive variant tNASP(Q289X) subjects the expression of tNASP to nonsense-mediated decay. tNASP KO increases chromatin accessibility, promotes the active promoter state of genes enriched in synaptic signalling and leads to upregulated expression of genes in the neural signalling and immune signalling pathways. Compared with wild-type tNASP, tNASP(Q289X) enhances chromatin accessibility of the genes with enriched expression in the brain. RNA-seq revealed that genes involved in neural and immune signalling are affected by the tNASP mutation, consistent with the phenotypic impact and molecular effects of nasp-1 mutations in Caenorhabditis elegans. Two additional patients with ASD were found carrying deletion or deleterious mutation in the NASP gene.

CONCLUSION

We identified novel epigenetic mechanisms mediated by tNASP which may contribute to the pathogenesis of ASD and its immune comorbidity.

Human-derived fecal microbiota transplantation alleviates social deficits of the BTBR mouse model of autism through a potential mechanism involving vitamin B6 metabolism

Autism spectrum disorder (ASD) is a heterogeneous neurodevelopmental condition characterized by social communication deficiencies and stereotypic behaviors influenced by hereditary and/or environmental risk factors. There are currently no approved medications for treating the core symptoms of ASD. Human fecal microbiota transplantation (FMT) has emerged as a potential intervention to improve autistic symptoms, but the underlying mechanisms are not fully understood. In this study, we evaluated the effects of human-derived FMT on behavioral and multi-omics profiles of the BTBR mice, an established model for ASD. FMT effectively alleviated the social deficits in the BTBR mice and normalized their distinct plasma metabolic profile, notably reducing the elevated long-chain acylcarnitines. Integrative analysis linked these phenotypic changes to specific Bacteroides species and vitamin B6 metabolism. Indeed, vitamin B6 supplementation improved the social behaviors in BTBR mice. Collectively, these findings shed new light on the interplay between FMT and vitamin B6 metabolism and revealed a potential mechanism underlying the therapeutic role of FMT in ASD.IMPORTANCEAccumulating evidence supports the beneficial effects of human fecal microbiota transplantation (FMT) on symptoms associated with autism spectrum disorder (ASD). However, the precise mechanism by which FMT induces a shift in the microbiota and leads to symptom improvement remains incompletely understood. This study integrated data from colon-content metagenomics, colon-content metabolomics, and plasma metabolomics to investigate the effects of FMT treatment on the BTBR mouse model for ASD. The analysis linked the amelioration of social deficits following FMT treatment to the restoration of mitochondrial function and the modulation of vitamin B6 metabolism. Bacterial species and compounds with beneficial roles in vitamin B6 metabolism and mitochondrial function may further contribute to improving FMT products and designing novel therapies for ASD treatment.

X-chromosome inactivation in human iPSCs provides insight into X-regulated gene expression in autosomes

BACKGROUND

Variation in X chromosome inactivation (XCI) in human-induced pluripotent stem cells (hiPSCs) can impact their ability to model biological sex biases. The gene-wise landscape of X chromosome gene dosage remains unresolved in female hiPSCs. To characterize patterns of de-repression and escape from inactivation, we performed a systematic survey of allele specific expression in 165 female hiPSC lines.

RESULTS

XCI erosion is non-random and primarily affects genes that escape XCI in human tissues. Individual genes and cell lines vary in the frequency and degree of de-repression. Bi-allelic expression increases gradually after modest decrease of XIST in cultures, whose loss is commonly used to mark lines with eroded XCI. We identify three clusters of female lines at different stages of XCI. Increased XCI erosion amplifies female-biased expression at hypomethylated sites and regions normally occupied by repressive histone marks, lowering male-biased differences in the X chromosome. In autosomes, erosion modifies sex differences in a dose-dependent way. Male-biased genes are enriched for hypermethylated regions, and de-repression of XIST-bound autosomal genes in female lines attenuates normal male-biased gene expression in eroded lines. XCI erosion can compensate for a dominant loss of function effect in several disease genes.

CONCLUSIONS

We present a comprehensive view of X chromosome gene dosage in hiPSCs and implicate a direct mechanism for XCI erosion in regulating autosomal gene expression in trans. The uncommon and variable reactivation of X chromosome genes in female hiPSCs can provide insight into X chromosome's role in regulating gene expression and sex differences in humans.

Prenatal diagnosis and pregnancy outcomes in fetuses with ventriculomegaly

Objective

Genetic etiology plays a critical role in fetal ventriculomegaly (VM). However, the studies on chromosomal copy number variants (CNVs) in fetal VM are limited. This study aimed to investigate the chromosomal CNVs in fetuses with mild to moderate VM, and explore its genotype-phenotype correlation.

Methods

A total of 242 fetuses with mild to moderate VM detected by prenatal ultrasound were enrolled in our study from October 2018 to October 2022. All cases underwent chromosomal microarray analysis (CMA) and G-banding simultaneously. All VM cases were classified different subgroups according to the maternal age, severity, VM distribution and presence/absence of other ultrasound abnormalities. The pregnancy outcomes and health conditions after birth were followed up. We also performed a pooled analysis regarding likely pathogenic and pathogenic CNVs (LP/P CNVs) for VM.

Results

The detection rate of chromosomal abnormalities by karyotyping was 9.1% (22/242), whereas it was 16.5% (40/242) when CMA was conducted (P < 0.05). The total detection rate of chromosomal abnormalities by karyotyping and CMA was 21.1% (51/242). A 12.0% incremental yield of CMA over karyotyping was observed. The detection rate of total genetic variants in fetuses with bilateral VM was significantly higher than in fetuses with unilateral VM (30.0% vs. 16.7%, P = 0.017). No significant differences were discovered between isolated VM and non-isolated VM, or between mild and moderate VM, or between advanced maternal age (AMA) and non-AMA (all P > 0.05). 28 fetuses with VM were terminated and 214 fetuses were delivered: one presented developmental delay and one presented congenital heart disease. The VM cases with both positive CMA and karyotypic results had a higher rate of termination of pregnancy than those with either a positive CMA or karyotypic result, or both negative testing results (P < 0.001).

Conclusion

The combination of CMA and karyotyping should be adopted to improve the positive detection rate of chromosomal abnormalities for VM. The total genetic abnormalities detected using both techniques would affect the final pregnancy outcomes. LP/P CNVs at 16p11.2, 17p13, and 22q11.21 were identified as the top three chromosomal hotspots associated with VM, which would enable genetic counselors to provide more precise genetic counseling for VM pregnancies.

Whole exome sequencing and proteomics-based investigation of the pathogenesis of coronary artery disease with diffuse long lesion

OBJECTIVES

The long-term prognosis of patients with coronary artery disease (CAD) with diffuse long lesion underwent coronary artery bypass graft (CABG) or percutaneous coronary intervention (PCI) remains worse. Here, we aimed to identify distinctive genes involved and offer novel insights into the pathogenesis of diffuse long lesion.

MATERIALS AND METHODS

Whole exome sequencing was performed on peripheral blood samples from 20 CAD patients with diffuse long lesion (CAD-DLL) and from 10 controls with focal lesion (CAD-FL) through a uniform pipeline. Proteomics analysis was conducted on the serum samples from 10 CAD-DLL patients and from 10 controls with CAD-FL by mass spectrometry. Bioinformatics analysis was performed to elucidate the involved genes, including functional annotation and protein-protein interaction analysis.

RESULTS

A total of 742 shared variant genes were found in CAD-DLL patients but not in controls. Of these, 46 genes were identified as high-frequency variant genes (≥ 4/20) distinctive genes. According to the consensus variant site, 148 shared variant sites were found in the CAD-DLL group. The lysosome and cellular senescence-related pathway may be the most significant pathway in diffuse long lesion. Following the DNA-protein combined analysis, eight genes were screened whose expression levels were altered at both DNA and protein levels. Among these genes, the MAN2A2 gene, the only one that was highly expressed at the protein level, was associated with metabolic and immune-inflammatory dysregulation.

CONCLUSIONS

Compared to individuals with CAD-FL, patients with CAD-DLL show additional variants. These findings contribute to the understanding of the mechanism of CAD-DLL and provide potential targets for the diagnosis and treatment of CAD-DLL.

Comparative analysis of the nucleus accumbens transcriptional features in multiple depressive animal models

Chronic stress is deemed a significant clinical contributor to depression. The use of animal models of chronic stress can fully reveal the complex pathological mechanisms and their changing trends in the pathogenesis of depression, which is crucial for both disease prevention and therapy. It is also unknown how various forms of stress differ in their impact on animal physiology and behavior. The nucleus accumbens (NAc), an essential brain area for the pathophysiology of depression, and its underlying neural mechanisms remain unclear. Here, we systematically compared transcriptional signatures in the NAc of four chronic stress models in rats: chronic unpredictable mild stress (CUMS), chronic social defeat stress (CSDS), learned helplessness (LH), chronic restraint stress (CRS). The majority of differentially expressed genes (DEGs) were unique to a single depression model, while the rank-rank hypergeometric overlap analysis showed that the CSDS and CRS models had the greatest overlap, and the CRS and CUMS models had the least. Then, we performed pathway analysis of the differential genes and found that the neuroactive ligand-receptor interaction pathway was significantly enriched not only in the LH, CRS and CSDS stress models, but also significantly enriched in stress genes that were also altered in at least two stress models. Finally, we found three hub genes (Dcx, Tnc and Wdfy4) by constructing co-expression networks for stress genes. In summary, our research has the potential to offer fresh insights into the molecular mechanisms underlying depression induced by different types of stress, highlighting both their similarities and differences. It may provide valuable clues for understanding the pathogenesis of depression.

KCTD10 regulates brain development by destabilizing brain disorder-associated protein KCTD13

KCTD10 belongs to the KCTD (potassiumchannel tetramerization domain) family, many members of which are associated with neuropsychiatric disorders. However, the biological function underlying the association with brain disorders remains to be explored. Here, we reveal that Kctd10 is highly expressed in neuronal progenitors and layer V neurons throughout brain development. Kctd10 deficiency triggers abnormal proliferation and differentiation of neuronal progenitors, reduced deep-layer (especially layer V) neurons, increased upper-layer neurons, and lowered brain size. Mechanistically, we screened and identified a unique KCTD10-interacting protein, KCTD13, associated with neurodevelopmental disorders. KCTD10 mediated the ubiquitination-dependent degradation of KCTD13 and KCTD10 ablation resulted in a considerable increase of KCTD13 expression in the developing cortex. KCTD13 overexpression in neuronal progenitors led to reduced proliferation and abnormal cell distribution, mirroring KCTD10 deficiency. Notably, mice with brain-specific Kctd10 knockout exhibited obvious motor deficits. This study uncovers the physiological function of KCTD10 and provides unique insights into the pathogenesis of neurodevelopmental disorders.

Heterozygous truncating variant of TAOK1 in a boy with periventricular nodular heterotopia: a case report and literature review of TAOK1-related neurodevelopmental disorders

BACKGROUND

Thousand and one amino-acid kinase 1 (TAOK1) encodes the MAP3K protein kinase TAO1, which has recently been displayed to be essential for neuronal maturation and cortical differentiation during early brain development. Heterozygous variants in TAOK1 have been reported in children with neurodevelopmental disorders, with or without macrocephaly, hypotonia and mild dysmorphic traits. Literature reports lack evidence of neuronal migration disorders in TAOK1 patients, although studies in animal models suggest this possibility.

CASE PRESENTATION

We provide a clinical description of a child with a neurodevelopmental disorder due to a novel TAOK1 truncating variant, whose brain magnetic resonance imaging displays periventricular nodular heterotopia.

CONCLUSIONS

To our knowledge, this is the first report of a neuronal migration disorder in a patient with a TAOK1-related neurodevelopmental disorder, thus supporting the hypothesized pathogenic mechanisms of TAOK1 defects.

Joint, multifaceted genomic analysis enables diagnosis of diverse, ultra-rare monogenic presentations

Genomics for rare disease diagnosis has advanced at a rapid pace due to our ability to perform "N-of-1" analyses on individual patients. The increasing sizes of ultra-rare, "N-of-1" disease cohorts internationally newly enables cohort-wide analyses for new discoveries, but well-calibrated statistical genetics approaches for jointly analyzing these patients are still under development.1,2 The Undiagnosed Diseases Network (UDN) brings multiple clinical, research and experimental centers under the same umbrella across the United States to facilitate and scale N-of-1 analyses. Here, we present the first joint analysis of whole genome sequencing data of UDN patients across the network. We apply existing and introduce new, well-calibrated statistical methods for prioritizing disease genes with de novo recurrence and compound heterozygosity. We also detect pathways enriched with candidate and known diagnostic genes. Our computational analysis, coupled with a systematic clinical review, recapitulated known diagnoses and revealed new disease associations. We make our gene-level findings and variant-level information across the cohort available in a public-facing browser (https://dbmi-bgm.github.io/udn-browser/). These results show that N-of-1 efforts should be supplemented by a joint genomic analysis across cohorts.

Topologically associating domains define the impact of de novo promoter variants on autism spectrum disorder risk

Whole-genome sequencing (WGS) studies of autism spectrum disorder (ASD) have demonstrated the roles of rare promoter de novo variants (DNVs). However, most promoter DNVs in ASD are not located immediately upstream of known ASD genes. In this study analyzing WGS data of 5,044 ASD probands, 4,095 unaffected siblings, and their parents, we show that promoter DNVs within topologically associating domains (TADs) containing ASD genes are significantly and specifically associated with ASD. An analysis considering TADs as functional units identified specific TADs enriched for promoter DNVs in ASD and indicated that common variants in these regions also confer ASD heritability. Experimental validation using human induced pluripotent stem cells (iPSCs) showed that likely deleterious promoter DNVs in ASD can influence multiple genes within the same TAD, resulting in overall dysregulation of ASD-associated genes. These results highlight the importance of TADs and gene-regulatory mechanisms in better understanding the genetic architecture of ASD.

Loss of NSD2 causes dysregulation of synaptic genes and altered H3K36 dimethylation in mice

Background: Epigenetic disruptions have been implicated in neurodevelopmental disorders. NSD2 is associated with developmental delay/intellectual disability; however, its role in brain development and function remains unclear. Methods: We performed transcriptomic and epigenetic analyses using Nsd2 knockout mice to better understand the role of NSD2 in the brain. Results and discussion: Transcriptomic analysis revealed that the loss of NSD2 caused dysregulation of genes related to synaptic transmission and formation. By analyzing changes in H3 lysine 36 dimethylation (H3K36me2), NSD2-mediated H3K36me2 mainly marked quiescent state regions and the redistribution of H3K36me2 occurred at transcribed genes and enhancers. By integrating transcriptomic and epigenetic data, we observed that H3K36me2 changes in a subset of dysregulated genes related to synaptic transmission and formation. These results suggest that NSD2 is involved in the regulation of genes important for neural function through H3K36me2. Our findings provide insights into the role of NSD2 and improve our understanding of epigenetic regulation in the brain.

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.

Utilizing hiPSC-derived oligodendrocytes to study myelin pathophysiology in neuropsychiatric and neurodegenerative disorders

Oligodendrocytes play a crucial role in our central nervous system (CNS) by myelinating axons for faster action potential conduction, protecting axons from degeneration, structuring the position of ion channels, and providing nutrients to neurons. Oligodendrocyte dysfunction and/or dysmyelination can contribute to a range of neurodegenerative diseases and neuropsychiatric disorders such as Multiple Sclerosis (MS), Leukodystrophy (LD), Schizophrenia (SCZ), and Autism Spectrum Disorder (ASD). Common characteristics identified across these disorders were either an inability of oligodendrocytes to remyelinate after degeneration or defects in oligodendrocyte development and maturation. Unfortunately, the causal mechanisms of oligodendrocyte dysfunction are still uncertain, and therapeutic targets remain elusive. Many studies rely on the use of animal models to identify the molecular and cellular mechanisms behind these disorders, however, such studies face species-specific challenges and therefore lack translatability. The use of human induced pluripotent stem cells (hiPSCs) to model neurological diseases is becoming a powerful new tool, improving our understanding of pathophysiology and capacity to explore therapeutic targets. Here, we focus on the application of hiPSC-derived oligodendrocyte model systems to model disorders caused by oligodendrocyte dysregulation.

An evolutionary perspective on complex neuropsychiatric disease

The forces of evolution-mutation, selection, migration, and genetic drift-shape the genetic architecture of human traits, including the genetic architecture of complex neuropsychiatric illnesses. Studying these illnesses in populations that are diverse in genetic ancestry, historical demography, and cultural history can reveal how evolutionary forces have guided adaptation over time and place. A fundamental truth of shared human biology is that an allele responsible for a disease in anyone, anywhere, reveals a gene critical to the normal biology underlying that condition in everyone, everywhere. Understanding the genetic causes of neuropsychiatric disease in the widest possible range of human populations thus yields the greatest possible range of insight into genes critical to human brain development. In this perspective, we explore some of the relationships between genes, adaptation, and history that can be illuminated by an evolutionary perspective on studies of complex neuropsychiatric disease in diverse populations.

Prenatal prevalence and postnatal manifestations of 16p11.2 deletions: A new insights into neurodevelopmental disorders based on clinical investigations combined with multi-omics analysis

BACKGROUND

The 16p11.2 deletion is one of the most common genetic aetiologies of neurodevelopmental disorders (NDDs). The prenatal phenotype of 16p11.2 deletion and the potential mechanism associated with postnatal clinical manifestations were largely unknow. We revealed the developmental trajectories of 16p11.2 deletion from the prenatal to postnatal periods and to identify key signaling pathways and candidate genes contributing to neurodevelopmental abnormalities.

METHODS

In this 5-y retrospective cohort study, women with singleton pregnancies who underwent amniocentesis for chromosomal abnormalities were included. Test of copy-number variations (CNVs) involved single nucleotide polymorphism-array and CNV-seq was performed to detected 16p11.2 deletion. For infants born carrying the 16p11.2 deletion, neurological and intellectual evaluations using the Chinese version of the Gesell Development Scale. For patients observed to have vertebral malformations, Sanger sequencing for T-C-A haplotype of TBX6 was performed. For those infants with clinical manifestations, whole-exome sequencing was consecutively performed in trios to rule out single-gene diseases, and transcriptomics combined with untargeted metabolomics were performed.

RESULTS

The prevalence of 16p11.2 deletion was 0.063% (55/86,035) in the prenatal period. Up to 80% (20/25) of the 16p11.2 deletions were proven de novo by parental confirmation. Approximately half of 16p11.2 deletions (28/55) were detected with prenatal abnormal ultrasound findings. Vertebral malformations were identified as the most distinctive structural malformations and were enriched in fetuses with 16p11.2 deletions compared with controls (90.9‰ [5/55] vs. 8.4‰ [72/85,980]; P < 0.001). All 5 fetuses with vertebral malformations were confirmed to have the TBX6 haplotype of T-C-A. Overall, 47.6% (10/21) infants birthed were diagnosed with NDDs of different degrees. Language impairment was the predominant manifestation (7/10; 70.0%), followed by motor delay (5/10; 50%). Multi-omics analysis indicated that MAPK3 was the central hub of the differentially expressed gene (DEG) network. We firstly reported that histidine-associated metabolism may be the core metabolic pathway related to the 16p11.2 deletion.

CONCLUSION

We demonstrated the prenatal presentation, incomplete penetrance and variable expressivity of the 16p11.2 deletion. We identified vertebral malformations were the most distinctive prenatal phenotypes, and language impairment was the predominant postnatal manifestation. Most of the 16p11.2 deletion was de novo. Meanwhile, we suggested that MAPK3 and histidine-associated metabolism may contribute to neurodevelopmental abnormalities of 16p11.2 deletion.

Astrocytic β-catenin signaling via TCF7L2 regulates synapse development and social behavior

The Wnt/β-catenin pathway contains multiple high-confidence risk genes that are linked to neurodevelopmental disorders, including autism spectrum disorder. However, its ubiquitous roles across brain cell types and developmental stages have made it challenging to define its impact on neural circuit development and behavior. Here, we show that TCF7L2, which is a key transcriptional effector of the Wnt/β-catenin pathway, plays a cell-autonomous role in postnatal astrocyte maturation and impacts adult social behavior. TCF7L2 was the dominant Wnt effector that was expressed in both mouse and human astrocytes, with a peak during astrocyte maturation. The conditional knockout of Tcf7l2 in postnatal astrocytes led to an enlargement of astrocytes with defective tiling and gap junction coupling. These mice also exhibited an increase in the number of cortical excitatory and inhibitory synapses and a marked increase in social interaction by adulthood. These data reveal an astrocytic role for developmental Wnt/β-catenin signaling in restricting excitatory synapse numbers and regulating adult social behavior.

Prioritizing de novo potential non-canonical splicing variants in neurodevelopmental disorders

BACKGROUND

Genomic variants outside of the canonical splicing site (±2) may generate abnormal mRNA splicing, which are defined as non-canonical splicing variants (NCSVs). However, the clinical interpretation of NCSVs in neurodevelopmental disorders (NDDs) is largely unknown.

METHODS

We investigated the contribution of NCSVs to NDDs from 345,787 de novo variants (DNVs) in 47,574 patients with NDDs. We performed functional enrichment and protein-protein interaction analysis to assess the association between genes carrying prioritised NCSVs and NDDs. Minigene was used to validate the impact of NCSVs on mRNA splicing.

FINDINGS

We observed significantly more NCSVs (p = 0.02, odds ratio [OR] = 2.05) among patients with NDD than in controls. Both canonical splicing variants (CSVs) and NCSVs contributed to an equal proportion of patients with NDD (0.76% vs. 0.82%). The candidate genes carrying NCSVs were associated with glutamatergic synapse and chromatin remodelling. Minigene successfully validated 59 of 79 (74.68%) NCSVs that led to abnormal splicing in 40 candidate genes, and 9 of the genes (ARID1B, KAT6B, TCF4, SMARCA2, SHANK3, PDHA1, WDR45, SCN2A, SYNGAP1) harboured recurrent NCSVs with the same variant present in more than two unrelated patients with NDD. Moreover, 36 of 59 (61.02%) NCSVs are novel clinically relevant variants, including 34 unreported and 2 clinically conflicting interpretations or of uncertain significance NCSVs in the ClinVar database.

INTERPRETATION

This study highlights the common pathology and clinical importance of NCSVs in unsolved patients with NDD.

FUNDING

The present study was funded by grants from the National Natural Science Foundation of China, China Postdoctoral Science Foundation, the Hunan Youth Science and Technology Innovation Talent Project, the Provincial Natural Science Foundation of Hunan, The Scientific Research Program of FuRong laboratory, and the Natural Science Project of the University of Anhui Province.

Whole Genome Analysis of SNV and Indel Polymorphism in Common Marmosets (Callithrix jacchus)

The common marmoset (Callithrix jacchus) is one of the most widely used nonhuman primate models of human disease. Owing to limitations in sequencing technology, early genome assemblies of this species using short-read sequencing suffered from gaps. In addition, the genetic diversity of the species has not yet been adequately explored. Using long-read genome sequencing and expert annotation, we generated a high-quality genome resource creating a 2.898 Gb marmoset genome in which most of the euchromatin portion is assembled contiguously (contig N50 = 25.23 Mbp, scaffold N50 = 98.2 Mbp). We then performed whole genome sequencing on 84 marmosets sampling the genetic diversity from several marmoset research centers. We identified a total of 19.1 million single nucleotide variants (SNVs), of which 11.9 million can be reliably mapped to orthologous locations in the human genome. We also observed 2.8 million small insertion/deletion variants. This dataset includes an average of 5.4 million SNVs per marmoset individual and a total of 74,088 missense variants in protein-coding genes. Of the 4956 variants orthologous to human ClinVar SNVs (present in the same annotated gene and with the same functional consequence in marmoset and human), 27 have a clinical significance of pathogenic and/or likely pathogenic. This important marmoset genomic resource will help guide genetic analyses of natural variation, the discovery of spontaneous functional variation relevant to human disease models, and the development of genetically engineered marmoset disease models.

Toward a CRISPR understanding of gene function in human brain development

Genome-wide association studies pinpoint genetic risk factors for neurodevelopmental disorders (NDDs), but the next challenge is to understand the mechanisms through which these genes affect brain development. Two recent CRISPR screens in human brain organoids1,2 interrogate the function of risk genes for autism spectrum disorder and other NDDs.

The functions of FOXP transcription factors and their regulation by post-translational modifications

The forkhead box subfamily P (FOXP) of transcription factors, consisting of FOXP1, FOXP2, FOXP3, and FOXP4, is involved in the regulation of multisystemic functioning. Disruption of the transcriptional activity of FOXP proteins leads to neurodevelopmental disorders and immunological diseases, as well as the suppression or promotion of carcinogenesis. The transcriptional activities of FOXP proteins are directly or indirectly regulated by diverse post-translational modifications, including phosphorylation, ubiquitination, SUMOylation, acetylation, O-GlcNAcylation, and methylation. Here, we discuss how post-translational modifications modulate the multiple functions of FOXP proteins and examine the implications for tumorigenesis and cancer therapy.

Distinct disease mutations in DNMT3A result in a spectrum of behavioral, epigenetic, and transcriptional deficits

Phenotypic heterogeneity in monogenic neurodevelopmental disorders can arise from differential severity of variants underlying disease, but how distinct alleles drive variable disease presentation is not well understood. Here, we investigate missense mutations in DNA methyltransferase 3A (DNMT3A), a DNA methyltransferase associated with overgrowth, intellectual disability, and autism, to uncover molecular correlates of phenotypic heterogeneity. We generate a Dnmt3aP900L/+ mouse mimicking a mutation with mild to moderate severity and compare phenotypic and epigenomic effects with a severe R878H mutation. P900L mutants exhibit core growth and behavioral phenotypes shared across models but show subtle epigenomic changes, while R878H mutants display extensive disruptions. We identify mutation-specific dysregulated genes that may contribute to variable disease severity. Shared transcriptomic disruption identified across mutations overlaps dysregulation observed in other developmental disorder models and likely drives common phenotypes. Together, our findings define central drivers of DNMT3A disorders and illustrate how variable epigenomic disruption contributes to phenotypic heterogeneity in neurodevelopmental disease.

Quantitative Relaxometry Assessment of Brain Microstructural Abnormality of Preschool Children With Autism Spectrum Disorder With Synthetic Magnetic Resonance Imaging

OBJECTIVE

This study aimed to perform an assessment of brain microstructure in children with autism aged 2 to 5 years using relaxation times acquired by synthetic magnetic resonance imaging.

MATERIALS AND METHODS

Thirty-four children with autism spectrum disorder (ASD) (ASD group) and 17 children with global developmental delay (GDD) (GDD group) were enrolled, and synthetic magnetic resonance imaging was performed to obtain T1 and T2 relaxation times. The differences in brain relaxation times between the 2 groups of children were compared, and the correlation between significantly changed T1/T2 and clinical neuropsychological scores in the ASD group was analyzed.

RESULTS

Compared with the GDD group, shortened T1 relaxation times in the ASD group were distributed in the genu of corpus callosum (GCC) ( P = 0.003), splenium of corpus callosum ( P = 0.002), and right thalamus (TH) ( P = 0.014), whereas shortened T2 relaxation times in the ASD group were distributed in GCC ( P = 0.011), left parietal white matter ( P = 0.035), and bilateral TH (right, P = 0.014; left, P = 0.016). In the ASD group, the T2 of the left parietal white matter is positively correlated with gross motor (developmental quotient [DQ] 2) and personal-social behavior (DQ5), respectively ( r = 0.377, P = 0.028; r = 0.392, P = 0.022); the T2 of the GCC was positively correlated with DQ5 ( r = 0.404, P = 0.018); and the T2 of the left TH is positively correlated with DQ2 and DQ5, respectively ( r = 0.433, P = 0.009; r = 0.377, P = 0.028). All significantly changed relaxation values were not significantly correlated with Childhood Autism Rating Scale scores.

CONCLUSIONS

The shortened relaxometry times in the brain of children with ASD may be associated with the increased myelin content and decreased water content in the brain of children with ASD in comparison with GDD, contributing the understanding of the pathophysiology of ASD. Therefore, the T1 and T2 relaxometry may be used as promising imaging markers for ASD diagnosis.

Assembloid CRISPR screens reveal impact of disease genes in human neurodevelopment

The assembly of cortical circuits involves the generation and migration of interneurons from the ventral to the dorsal forebrain1-3, which has been challenging to study at inaccessible stages of late gestation and early postnatal human development4. Autism spectrum disorder and other neurodevelopmental disorders (NDDs) have been associated with abnormal cortical interneuron development5, but which of these NDD genes affect interneuron generation and migration, and how they mediate these effects remains unknown. We previously developed a platform to study interneuron development and migration in subpallial organoids and forebrain assembloids6. Here we integrate assembloids with CRISPR screening to investigate the involvement of 425 NDD genes in human interneuron development. The first screen aimed at interneuron generation revealed 13 candidate genes, including CSDE1 and SMAD4. We subsequently conducted an interneuron migration screen in more than 1,000 forebrain assembloids that identified 33 candidate genes, including cytoskeleton-related genes and the endoplasmic reticulum-related gene LNPK. We discovered that, during interneuron migration, the endoplasmic reticulum is displaced along the leading neuronal branch before nuclear translocation. LNPK deletion interfered with this endoplasmic reticulum displacement and resulted in abnormal migration. These results highlight the power of this CRISPR-assembloid platform to systematically map NDD genes onto human development and reveal disease mechanisms.

Rare variation in noncoding regions with evolutionary signatures contributes to autism spectrum disorder risk

Little is known about the role of noncoding regions in the etiology of autism spectrum disorder (ASD). We examined three classes of noncoding regions: Human Accelerated Regions (HARs), which show signatures of positive selection in humans; experimentally validated neural Vista Enhancers (VEs); and conserved regions predicted to act as neural enhancers (CNEs). Targeted and whole genome analysis of >16,600 samples and >4900 ASD probands revealed that likely recessive, rare, inherited variants in HARs, VEs, and CNEs substantially contribute to ASD risk in probands whose parents share ancestry, which enriches for recessive contributions, but modestly, if at all, in simplex family structures. We identified multiple patient variants in HARs near IL1RAPL1 and in a VE near SIM1 and showed that they change enhancer activity. Our results implicate both human-evolved and evolutionarily conserved noncoding regions in ASD risk and suggest potential mechanisms of how changes in regulatory regions can modulate social behavior.

Developing a cluster-based approach for deciphering complexity in individuals with neurodevelopmental differences

Objective

Individuals with neurodevelopmental disorders such as global developmental delay (GDD) present both genotypic and phenotypic heterogeneity. This diversity has hampered developing of targeted interventions given the relative rarity of each individual genetic etiology. Novel approaches to clinical trials where distinct, but related diseases can be treated by a common drug, known as basket trials, which have shown benefits in oncology but have yet to be used in GDD. Nonetheless, it remains unclear how individuals with GDD could be clustered. Here, we assess two different approaches: agglomerative and divisive clustering.

Methods

Using the largest cohort of individuals with GDD, which is the Deciphering Developmental Disorders (DDD), characterized using a systematic approach, we extracted genotypic and phenotypic information from 6,588 individuals with GDD. We then used a k-means clustering (divisive) and hierarchical agglomerative clustering (HAC) to identify subgroups of individuals. Next, we extracted gene network and molecular function information with regard to the clusters identified by each approach.

Results

HAC based on phenotypes identified in individuals with GDD revealed 16 clusters, each presenting with one dominant phenotype displayed by most individuals in the cluster, along with other minor phenotypes. Among the most common phenotypes reported were delayed speech, absent speech, and seizure. Interestingly, each phenotypic cluster molecularly included several (3-12) gene sub-networks of more closely related genes with diverse molecular function. k-means clustering also segregated individuals harboring those phenotypes, but the genetic pathways identified were different from the ones identified from HAC.

Conclusion

Our study illustrates how divisive (k-means) and agglomerative clustering can be used in order to group individuals with GDD for future basket trials. Moreover, the result of our analysis suggests that phenotypic clusters should be subdivided into molecular sub-networks for an increased likelihood of successful treatment. Finally, a combination of both agglomerative and divisive clustering may be required for developing of a comprehensive treatment.

An unsupervised deep learning framework for predicting human essential genes from population and functional genomic data

BACKGROUND

The ability to accurately predict essential genes intolerant to loss-of-function (LOF) mutations can dramatically improve the identification of disease-associated genes. Recently, there have been numerous computational methods developed to predict human essential genes from population genomic data. While the existing methods are highly predictive of essential genes of long length, they have limited power in pinpointing short essential genes due to the sparsity of polymorphisms in the human genome.

RESULTS

Motivated by the premise that population and functional genomic data may provide complementary evidence for gene essentiality, here we present an evolution-based deep learning model, DeepLOF, to predict essential genes in an unsupervised manner. Unlike previous population genetic methods, DeepLOF utilizes a novel deep learning framework to integrate both population and functional genomic data, allowing us to pinpoint short essential genes that can hardly be predicted from population genomic data alone. Compared with previous methods, DeepLOF shows unmatched performance in predicting ClinGen haploinsufficient genes, mouse essential genes, and essential genes in human cell lines. Notably, at a false positive rate of 5%, DeepLOF detects 50% more ClinGen haploinsufficient genes than previous methods. Furthermore, DeepLOF discovers 109 novel essential genes that are too short to be identified by previous methods.

CONCLUSION

The predictive power of DeepLOF shows that it is a compelling computational method to aid in the discovery of essential genes.

Acorn: an R package for de novo variant analysis

BACKGROUND

The study of de novo variation is important for assessing biological characteristics of new variation and for studies related to human phenotypes. Software programs exist to call de novo variants and programs also exist to test the burden of these variants in genomic regions; however, I am unaware of a program that fits in between these two aspects of de novo variant assessment. This intermediate space is important for assessing the quality of de novo variants and to understand the characteristics of the callsets. For this reason, I developed an R package called acorn.

RESULTS

Acorn is an R package that examines various features of de novo variants including subsetting the data by individual(s), variant type, or genomic region; calculating features including variant change counts, variant lengths, and presence/absence at CpG sites; and characteristics of parental age in relation to de novo variant counts.

CONCLUSIONS

Acorn is an R package that fills a critical gap in assessing de novo variants and will be of benefit to many investigators studying de novo variation.

Sex difference contributes to phenotypic diversity in individuals with neurodevelopmental disorders

Objective

Gain a better understanding of sex-specific differences in individuals with global developmental delay (GDD), with a focus on phenotypes and genotypes.

Methods

Using the Deciphering Developmental Disorders (DDD) dataset, we extracted phenotypic information from 6,588 individuals with GDD and then identified statistically significant variations in phenotypes and genotypes based on sex. We compared genes with pathogenic variants between sex and then performed gene network and molecular function enrichment analysis and gene expression profiling between sex. Finally, we contrasted individuals with autism as an associated condition.

Results

We identified significantly differentially expressed phenotypes in males vs. females individuals with GDD. Autism and macrocephaly were significantly more common in males whereas microcephaly and stereotypies were more common in females. Importantly, 66% of GDD genes with pathogenic variants overlapped between both sexes. In the cohort, males presented with only slightly increased X-linked genes (9% vs. 8%, respectively). Individuals from both sexes harbored a similar number of pathogenic variants overall (3) but females presented with a significantly higher load for GDD genes with high intolerance to loss of function. Sex difference in gene expression correlated with genes identified in a sex specific manner. While we identified sex-specific GDD gene mutations, their pathways overlapped. Interestingly, individuals with GDD but also co-morbid autism phenotypes, we observed distinct mutation load, pathways and phenotypic presentation.

Conclusion

Our study shows for the first time that males and females with GDD present with significantly different phenotypes. Moreover, while most GDD genes overlapped, some genes were found uniquely in each sex. Surprisingly they shared similar molecular functions. Sorting genes by predicted tolerance to loss of function (pLI) led to identifying an increased mutation load in females with GDD, suggesting potentially a tolerance to GDD genes of higher pLI compared to overall GDD genes. Finally, we show that considering associated conditions (for instance autism) may influence the genomic underpinning found in individuals with GDD and highlight the importance of comprehensive phenotyping.

VBASS enables integration of single cell gene expression data in Bayesian association analysis of rare variants

Rare or de novo variants have substantial contribution to human diseases, but the statistical power to identify risk genes by rare variants is generally low due to rarity of genotype data. Previous studies have shown that risk genes usually have high expression in relevant cell types, although for many conditions the identity of these cell types are largely unknown. Recent efforts in single cell atlas in human and model organisms produced large amount of gene expression data. Here we present VBASS, a Bayesian method that integrates single-cell expression and de novo variant (DNV) data to improve power of disease risk gene discovery. VBASS models disease risk prior as a function of expression profiles, approximated by deep neural networks. It learns the weights of neural networks and parameters of Gamma-Poisson likelihood models of DNV counts jointly from expression and genetics data. On simulated data, VBASS shows proper error rate control and better power than state-of-the-art methods. We applied VBASS to published datasets and identified more candidate risk genes with supports from literature or data from independent cohorts. VBASS can be generalized to integrate other types of functional genomics data in statistical genetics analysis.

Methyl-CpG binding domain 2 (Mbd2) is an epigenetic regulator of autism-risk genes and cognition

The Methyl-CpG-Binding Domain Protein family has been implicated in neurodevelopmental disorders. The Methyl-CpG-binding domain 2 (Mbd2) binds methylated DNA and was shown to play an important role in cancer and immunity. Some evidence linked this protein to neurodevelopment. However, its exact role in neurodevelopment and brain function is mostly unknown. Here we show that Mbd2-deficiency in mice (Mbd2-/-) results in deficits in cognitive, social and emotional functions. Mbd2 binds regulatory DNA regions of neuronal genes in the hippocampus and loss of Mbd2 alters the expression of hundreds of genes with a robust down-regulation of neuronal gene pathways. Further, a genome-wide DNA methylation analysis found an altered DNA methylation pattern in regulatory DNA regions of neuronal genes in Mbd2-/- mice. Differentially expressed genes significantly overlap with gene-expression changes observed in brains of Autism Spectrum Disorder (ASD) individuals. Notably, downregulated genes are significantly enriched for human ortholog ASD risk genes. Observed hippocampal morphological abnormalities were similar to those found in individuals with ASD and ASD rodent models. Hippocampal Mbd2 knockdown partially recapitulates the behavioral phenotypes observed in Mbd2-/- mice. These findings suggest that Mbd2 is a novel epigenetic regulator of genes that are associated with ASD in humans. Mbd2 loss causes behavioral alterations that resemble those found in ASD individuals.

Phenotypic effects of genetic variants associated with autism

While over 100 genes have been associated with autism, little is known about the prevalence of variants affecting them in individuals without a diagnosis of autism. Nor do we fully appreciate the phenotypic diversity beyond the formal autism diagnosis. Based on data from more than 13,000 individuals with autism and 210,000 undiagnosed individuals, we estimated the odds ratios for autism associated to rare loss-of-function (LoF) variants in 185 genes associated with autism, alongside 2,492 genes displaying intolerance to LoF variants. In contrast to autism-centric approaches, we investigated the correlates of these variants in individuals without a diagnosis of autism. We show that these variants are associated with a small but significant decrease in fluid intelligence, qualification level and income and an increase in metrics related to material deprivation. These effects were larger for autism-associated genes than in other LoF-intolerant genes. Using brain imaging data from 21,040 individuals from the UK Biobank, we could not detect significant differences in the overall brain anatomy between LoF carriers and non-carriers. Our results highlight the importance of studying the effect of the genetic variants beyond categorical diagnosis and the need for more research to understand the association between these variants and sociodemographic factors, to best support individuals carrying these variants.

NeuroCNVscore: a tissue-specific framework to prioritise the pathogenicity of CNVs in neurodevelopmental disorders

BACKGROUND

Neurodevelopmental disorders (NDDs) are associated with altered development of the brain especially in childhood. Copy number variants (CNVs) play a crucial role in the genetic aetiology of NDDs by disturbing gene expression directly at linear sequence or remotely at three-dimensional genome level in a tissue-specific manner. Despite the substantial increase in NDD studies employing whole-genome sequencing, there is no specific tool for prioritising the pathogenicity of CNVs in the context of NDDs.

METHODS

Using an XGBoost classifier, we integrated 189 features that represent genomic sequences, gene information and functional/genomic segments for evaluating genome-wide CNVs in a neuro/brain-specific manner, to develop a new tool, neuroCNVscore. We used Human Phenotype Ontology to construct an independent NDD-related set.

RESULTS

Our neuroCNVscore framework (https://github.com/lxsbch/neuroCNVscore) achieved high predictive performance (precision recall=0.82; area under curve=0.85) and outperformed an existing reference method SVScore. Notably, the predicted pathogenic CNVs showed enrichment in known genes associated with autism.

CONCLUSIONS

NeuroCNVscore prioritises functional, deleterious and pathogenic CNVs in NDDs at whole genome-wide level, which is important for genetic studies and clinical genomic screening of NDDs as well as for providing novel biological insights into NDDs.

Loss-of-function variants in CUL3 cause a syndromic neurodevelopmental disorder

Purpose

De novo variants in CUL3 (Cullin-3 ubiquitin ligase) 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 multi-center 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 35 individuals with heterozygous CUL3 variants presenting a syndromic NDD characterized by intellectual disability with or without autistic features. Of these, 33 have loss-of-function (LoF) and two 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 . Specifically, we show that cyclin E1 (CCNE1) and 4E-BP1 (EIF4EBP1), two prominent substrates of CUL3, fail to be targeted for proteasomal degradation in patient-derived cells.

Conclusion

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.

Jansen-de Vries syndrome: Expansion of the PPM1D clinical and phenotypic spectrum in 34 families

Jansen-de Vries syndrome (JdVS) is a neurodevelopmental condition attributed to pathogenic variants in Exons 5 and 6 of PPM1D. As the full phenotypic spectrum and natural history remain to be defined, we describe a large cohort of children and adults with JdVS. This is a retrospective cohort study of 37 individuals from 34 families with disease-causing variants in PPM1D leading to JdVS. Clinical data were provided by treating physicians and/or families. Of the 37 individuals, 27 were male and 10 female, with median age 8.75 years (range 8 months to 62 years). Four families document autosomal dominant transmission, and 32/34 probands were diagnosed via exome sequencing. The facial gestalt, including a broad forehead and broad mouth with a thin and tented upper lip, was most recognizable between 18 and 48 months of age. Common manifestations included global developmental delay (35/36, 97%), hypotonia (25/34, 74%), short stature (14/33, 42%), constipation (22/31, 71%), and cyclic vomiting (6/35, 17%). Distinctive personality traits include a hypersocial affect (21/31, 68%) and moderate-to-severe anxiety (18/28, 64%). In conclusion, JdVS is a clinically recognizable neurodevelopmental syndrome with a characteristic personality and distinctive facial features. The association of pathogenic variants in PPM1D with cyclic vomiting bears not only medical attention but also further pathogenic and mechanistic evaluation.

The Central Noradrenergic System in Neurodevelopmental Disorders: Merging Experimental and Clinical Evidence

The aim of this article is to highlight the potential role of the locus-coeruleus-noradrenergic (LC-NA) system in neurodevelopmental disorders (NdDs). The LC is the main brain noradrenergic nucleus, key in the regulation of arousal, attention, and stress response, and its early maturation and sensitivity to perinatal damage make it an interesting target for translational research. Clinical data shows the involvement of the LC-NA system in several NdDs, suggesting a pathogenetic role in the development of such disorders. In this context, a new neuroimaging tool, LC Magnetic Resonance Imaging (MRI), has been developed to visualize the LC in vivo and assess its integrity, which could be a valuable tool for exploring morphological alterations in NdD in vivo in humans. New animal models may be used to test the contribution of the LC-NA system to the pathogenic pathways of NdD and to evaluate the efficacy of NA-targeting drugs. In this narrative review, we provide an overview of how the LC-NA system may represent a common pathophysiological and pathogenic mechanism in NdD and a reliable target for symptomatic and disease-modifying drugs. Further research is needed to fully understand the interplay between the LC-NA system and NdD.

Prediction of Neurodevelopmental Disorders Based on De Novo Coding Variation

The early detection of neurodevelopmental disorders (NDDs) can significantly improve patient outcomes. The differential burden of non-synonymous de novo mutation among NDD cases and controls indicates that de novo coding variation can be used to identify a subset of samples that will likely display an NDD phenotype. Thus, we have developed an approach for the accurate prediction of NDDs with very low false positive rate (FPR) using de novo coding variation for a small subset of cases. We use a shallow neural network that integrates de novo likely gene-disruptive and missense variants, measures of gene constraint, and conservation information to predict a small subset of NDD cases at very low FPR and prioritizes NDD risk genes for future clinical study.

Translating precision medicine for autism spectrum disorder: A pressing need

Autism spectrum disorder (ASD) is a heterogenous group of neurodevelopmental disorders (NDDs) with a high unmet medical need. Currently, ASD is diagnosed according to behavior-based criteria that overlook clinical and genomic heterogeneity, thus repeatedly resulting in failed clinical trials. Here, we summarize the scientific evidence pointing to the pressing need to create a precision medicine framework for ASD and other NDDs. We discuss the role of omics and systems biology to characterize more homogeneous disease subtypes with different underlying pathophysiological mechanisms and to determine corresponding tailored treatments. Finally, we provide recent initiatives towards tackling the complexity in NDDs for precision medicine and cost-effective drug discovery.

Distinct disease mutations in DNMT3A result in a spectrum of behavioral, epigenetic, and transcriptional deficits

Phenotypic heterogeneity is a common feature of monogenic neurodevelopmental disorders that can arise from differential severity of missense variants underlying disease, but how distinct alleles impact molecular mechanisms to drive variable disease presentation is not well understood. Here, we investigate missense mutations in the DNA methyltransferase DNMT3A associated with variable overgrowth, intellectual disability, and autism, to uncover molecular correlates of phenotypic heterogeneity in neurodevelopmental disease. We generate a DNMT3A P900L/+ mouse model mimicking a disease mutation with mild-to-moderate severity and compare phenotypic and epigenomic effects with a severe R878H mutation. We show that the P900L mutation leads to disease-relevant overgrowth, obesity, and social deficits shared across DNMT3A disorder models, while the R878H mutation causes more extensive epigenomic disruption leading to differential dysregulation of enhancers elements. We identify distinct gene sets disrupted in each mutant which may contribute to mild or severe disease, and detect shared transcriptomic disruption that likely drives common phenotypes across affected individuals. Finally, we demonstrate that core gene dysregulation detected in DNMT3A mutant mice overlaps effects in other developmental disorder models, highlighting the importance of DNMT3A-deposited methylation in neurodevelopment. Together, these findings define central drivers of DNMT3A disorders and illustrate how variable disruption of transcriptional mechanisms can drive the spectrum of phenotypes in neurodevelopmental disease.

UBE3A expression during early postnatal brain development is required for proper dorsomedial striatal maturation

Angelman syndrome (AS) is a severe neurodevelopmental disorder (NDD) caused by loss of functional ubiquitin protein ligase E3A (UBE3A). Previous studies showed that UBE3A plays an important role in the first postnatal weeks of mouse brain development, but its precise role is unknown. Since impaired striatal maturation has been implicated in several mouse models for NDDs, we studied the importance of UBE3A in striatal maturation. We used inducible Ube3a mouse models to investigate the maturation of medium spiny neurons (MSNs) from dorsomedial striatum. MSNs of mutant mice matured properly till postnatal day 15 (P15) but remained hyperexcitable with fewer excitatory synaptic events at later ages, indicative of stalled striatal maturation in Ube3a mice. Reinstatement of UBE3A expression at P21 fully restored MSN excitability but only partially restored synaptic transmission and the operant conditioning behavioral phenotype. Gene reinstatement at P70 failed to rescue both electrophysiological and behavioral phenotypes. In contrast, deletion of Ube3a after normal brain development did not result in these electrophysiological and behavioral phenotypes. This study emphasizes the role of UBE3A in striatal maturation and the importance of early postnatal reinstatement of UBE3A expression to obtain a full rescue of behavioral phenotypes associated with striatal function in AS.

Dysfunction of Prkcaa Links Social Behavior Defects with Disturbed Circadian Rhythm in Zebrafish

Protein kinase Cα (PKCα/PRKCA) is a crucial regulator of circadian rhythm and is associated with human mental illnesses such as autism spectrum disorders and schizophrenia. However, the roles of PRKCA in modulating animal social behavior and the underlying mechanisms remain to be explored. Here we report the generation and characterization of prkcaa-deficient zebrafish (Danio rerio). The results of behavioral tests indicate that a deficiency in Prkcaa led to anxiety-like behavior and impaired social preference in zebrafish. RNA-sequencing analyses revealed the significant effects of the prkcaa mutation on the expression of the morning-preferring circadian genes. The representatives are the immediate early genes, including egr2a, egr4, fosaa, fosab and npas4a. The downregulation of these genes at night was attenuated by Prkcaa dysfunction. Consistently, the mutants demonstrated reversed day-night locomotor rhythm, which are more active at night than in the morning. Our data show the roles of PRKCA in regulating animal social interactions and link the social behavior defects with a disturbed circadian rhythm.

Comprehensive multi-omic profiling of somatic mutations in malformations of cortical development

Malformations of cortical development (MCD) are neurological conditions involving focal disruptions of cortical architecture and cellular organization that arise during embryogenesis, largely from somatic mosaic mutations, and cause intractable epilepsy. Identifying the genetic causes of MCD has been a challenge, as mutations remain at low allelic fractions in brain tissue resected to treat condition-related epilepsy. Here we report a genetic landscape from 283 brain resections, identifying 69 mutated genes through intensive profiling of somatic mutations, combining whole-exome and targeted-amplicon sequencing with functional validation including in utero electroporation of mice and single-nucleus RNA sequencing. Genotype-phenotype correlation analysis elucidated specific MCD gene sets associated with distinct pathophysiological and clinical phenotypes. The unique single-cell level spatiotemporal expression patterns of mutated genes in control and patient brains indicate critical roles in excitatory neurogenic pools during brain development and in promoting neuronal hyperexcitability after birth.

Cholinergic neurons in the basal forebrain are involved in behavioral abnormalities associated with Cul3 deficiency: Role of prefrontal cortex projections in cognitive deficits

Loss-of-function mutations of the gene Cul3 have been identified as a risk factor for autism-spectrum disorder (ASD), but the pathogenic mechanisms are not well understood. Conditional Cul3 ablation in cholinergic neurons of mice (ChatCRECul3F/+) recapitulated ASD-like social and sensory gating phenotypes and caused significant cognitive impairments, with diminished activity of cholinergic neurons in the basal forebrain (BF). Chemogenetic inhibition of BF cholinergic neurons in healthy mice induced similar social and cognitive deficits. Conversely, chemogenetic stimulation of BF cholinergic neurons in ChatCRECul3F/+ mice reversed abnormalities in sensory gating and cognition. Cortical hypofunction was also found after ChAT-specific Cul3 ablation and stimulation of cholinergic projections from the BF to the prefrontal cortex (PFC) mitigated cognitive deficits. Overall, we demonstrate that cholinergic dysfunction due to Cul3 deficiency is involved in ASD-like behavioral abnormalities, and that BF cholinergic neurons are particularly critical for cognitive component through their projections to the PFC.

Is bRaQCing bad? New roles for ribosome associated quality control factors in stress granule regulation

Maintenance of proteostasis is of utmost importance to cellular viability and relies on the coordination of many post-transcriptional processes to respond to stressful stimuli. Stress granules (SGs) are RNA–protein condensates that form after translation initiation is inhibited, such as during the integrated stress response (ISR), and may facilitate cellular adaptation to stress. The ribosome-associated quality control (RQC) pathway is a critical translation monitoring system that recognizes aberrant mRNAs encoding potentially toxic nascent peptides to target them for degradation. Both SG regulation and the RQC pathway are directly associated with translation regulation, thus it is of no surprise recent developments have demonstrated a connection between them. VCP's function in the stress activated RQC pathway, ribosome collisions activating the ISR, and the regulation of the 40S ribosomal subunit by canonical SG proteins during the RQC all connect SGs to the RQC pathway. Because mutations in genes that are involved in both SG and RQC regulation are associated with degenerative and neurological diseases, understanding the coordination and interregulation of SGs and RQC may shed light on disease mechanisms. This minireview will highlight recent advances in understanding how SGs and the RQC pathway interact in health and disease contexts.

de novo variant calling identifies cancer mutation signatures in the 1000 Genomes Project

Detection of de novo variants (DNVs) is critical for studies of disease-related variation and mutation rates. To accelerate DNV calling, we developed a graphics processing units-based workflow. We applied our workflow to whole-genome sequencing data from three parent-child sequenced cohorts including the Simons Simplex Collection (SSC), Simons Foundation Powering Autism Research (SPARK), and the 1000 Genomes Project (1000G) that were sequenced using DNA from blood, saliva, and lymphoblastoid cell lines (LCLs), respectively. The SSC and SPARK DNV callsets were within expectations for number of DNVs, percent at CpG sites, phasing to the paternal chromosome of origin, and average allele balance. However, the 1000G DNV callset was not within expectations and contained excessive DNVs that are likely cell line artifacts. Mutation signature analysis revealed 30% of 1000G DNV signatures matched B-cell lymphoma. Furthermore, we found variants in DNA repair genes and at Clinvar pathogenic or likely-pathogenic sites and significant excess of protein-coding DNVs in IGLL5; a gene known to be involved in B-cell lymphomas. Our study provides a new rapid DNV caller for the field and elucidates important implications of using sequencing data from LCLs for reference building and disease-related projects.

Top2a promotes the development of social behavior via PRC2 and H3K27me3

Little is understood about the embryonic development of sociality. We screened 1120 known drugs and found that embryonic inhibition of topoisomerase IIα (Top2a) resulted in lasting social deficits in zebrafish. In mice, prenatal Top2 inhibition caused defects in social interaction and communication, which are behaviors that relate to core symptoms of autism. Mutation of Top2a in zebrafish caused down-regulation of a set of genes highly enriched for genes associated with autism in humans. Both the Top2a-regulated and autism-associated gene sets have binding sites for polycomb repressive complex 2 (PRC2), a regulatory complex responsible for H3K27 trimethylation (H3K27me3). Moreover, both gene sets are highly enriched for H3K27me3. Inhibition of the PRC2 component Ezh2 rescued social deficits caused by Top2 inhibition. Therefore, Top2a is a key component of an evolutionarily conserved pathway that promotes the development of social behavior through PRC2 and H3K27me3.

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.

Estimating the Prevalence of De Novo Monogenic Neurodevelopmental Disorders from Large Cohort Studies

Rare diseases impact up to 400 million individuals globally. Of the thousands of known rare diseases, many are rare neurodevelopmental disorders (RNDDs) impacting children. RNDDs have proven to be difficult to assess epidemiologically for several reasons. The rarity of them makes it difficult to observe them in the population, there is clinical overlap among many disorders, making it difficult to assess the prevalence without genetic testing, and data have yet to be available to have accurate counts of cases. Here, we utilized large sequencing cohorts of individuals with rare, de novo monogenic disorders to estimate the prevalence of variation in over 11,000 genes among cohorts with developmental delay, autism spectrum disorder, and/or epilepsy. We found that the prevalence of many RNDDs is positively correlated to the previously estimated incidence. We identified the most often mutated genes among neurodevelopmental disorders broadly, as well as developmental delay and autism spectrum disorder independently. Finally, we assessed if social media group member numbers may be a valuable way to estimate prevalence. These data are critical for individuals and families impacted by these RNDDs, clinicians and geneticists in their understanding of how common diseases are, and for researchers to potentially prioritize research into particular genes or gene sets.

From bugs to bedside: functional annotation of human genetic variation for neurological disorders using invertebrate models

The exponential accumulation of DNA sequencing data has opened new avenues for discovering the causative roles of single-nucleotide polymorphisms (SNPs) in neurological diseases. The opportunities emerging from this are staggering, yet only as good as our abilities to glean insights from this surplus of information. Whereas computational biology continues to improve with respect to predictions and molecular modeling, the differences between in silico and in vivo analysis remain substantial. Invertebrate in vivo model systems represent technically advanced, experimentally mature, high-throughput, efficient and cost-effective resources for investigating a disease. With a decades-long track record of enabling investigators to discern function from DNA, fly (Drosophila) and worm (Caenorhabditis elegans) models have never been better poised to serve as living engines of discovery. Both of these animals have already proven useful in the classification of genetic variants as either pathogenic or benign across a range of neurodevelopmental and neurodegenerative disorders-including autism spectrum disorders, ciliopathies, amyotrophic lateral sclerosis, Alzheimer's and Parkinson's disease. Pathogenic SNPs typically display distinctive phenotypes in functional assays when compared with null alleles and frequently lead to protein products with gain-of-function or partial loss-of-function properties that contribute to neurological disease pathogenesis. The utility of invertebrates is logically limited by overt differences in anatomical and physiological characteristics, and also the evolutionary distance in genome structure. Nevertheless, functional annotation of disease-SNPs using invertebrate models can expedite the process of assigning cellular and organismal consequences to mutations, ascertain insights into mechanisms of action, and accelerate therapeutic target discovery and drug development for neurological conditions.

Tissue- and cell-type-specific molecular and functional signatures of 16p11.2 reciprocal genomic disorder across mouse brain and human neuronal models

Chromosome 16p11.2 reciprocal genomic disorder, resulting from recurrent copy-number variants (CNVs), involves intellectual disability, autism spectrum disorder (ASD), and schizophrenia, but the responsible mechanisms are not known. To systemically dissect molecular effects, we performed transcriptome profiling of 350 libraries from six tissues (cortex, cerebellum, striatum, liver, brown fat, and white fat) in mouse models harboring CNVs of the syntenic 7qF3 region, as well as cellular, transcriptional, and single-cell analyses in 54 isogenic neural stem cell, induced neuron, and cerebral organoid models of CRISPR-engineered 16p11.2 CNVs. Transcriptome-wide differentially expressed genes were largely tissue-, cell-type-, and dosage-specific, although more effects were shared between deletion and duplication and across tissue than expected by chance. The broadest effects were observed in the cerebellum (2,163 differentially expressed genes), and the greatest enrichments were associated with synaptic pathways in mouse cerebellum and human induced neurons. Pathway and co-expression analyses identified energy and RNA metabolism as shared processes and enrichment for ASD-associated, loss-of-function constraint, and fragile X messenger ribonucleoprotein target gene sets. Intriguingly, reciprocal 16p11.2 dosage changes resulted in consistent decrements in neurite and electrophysiological features, and single-cell profiling of organoids showed reciprocal alterations to the proportions of excitatory and inhibitory GABAergic neurons. Changes both in neuronal ratios and in gene expression in our organoid analyses point most directly to calretinin GABAergic inhibitory neurons and the excitatory/inhibitory balance as targets of disruption that might contribute to changes in neurodevelopmental and cognitive function in 16p11.2 carriers. Collectively, our data indicate the genomic disorder involves disruption of multiple contributing biological processes and that this disruption has relative impacts that are context specific.

Phenotypic Impact of Rare Potentially Damaging Copy Number Variation in Obsessive-Compulsive Disorder and Chronic Tic Disorders

BACKGROUND

Recent studies report an important-and previously underestimated-role of rare variation in risk of obsessive-compulsive disorder (OCD) and chronic tic disorders (CTD). Using data from a large epidemiological study, we evaluate the distribution of potentially damaging copy number variation (pdCNV) in OCD and CTD, examining associations between pdCNV and the phenotypes of probands, including a consideration of early- vs. late-diagnoses.

METHOD

The Obsessive-Compulsive Inventory-Revised (OCI-R) questionnaire was used to ascertain psychometric profiles of OCD probands. CNV were identified genome-wide using chromosomal microarray data.

RESULTS

For 993 OCD cases, 86 (9%) were identified as pdCNV carriers. The most frequent pdCNV found was at the 16p13.11 region. There was no significant association between pdCNV and the OCI-R total score. However, pdCNV was associated with Obsessing and Checking subscores. There was no significant difference in pdCNV frequency between early- vs. late-diagnosed OCD probands. Of the 217 CTD cases, 18 (8%) were identified as pdCNV carriers. CTD probands with pdCNV were significantly more likely to have co-occurring autism spectrum disorder (ASD).

CONCLUSIONS

pdCNV represents part of the risk architecture for OCD and CTD. If replicated, our findings suggest pdCNV impact some OCD symptoms. Genes within the 16p13.11 region are potential OCD risk genes.