Large-scale discovery of novel neurodevelopmental disorder-related genes through a unified analysis of single-nucleotide and copy number variants

Role of histone modifications in neurogenesis and neurodegenerative disease development

Progressive neuronal dysfunction and death are key features of neurodegenerative diseases; therefore, promoting neurogenesis in neurodegenerative diseases is crucial. With advancements in proteomics and high-throughput sequencing technology, it has been demonstrated that histone post-transcriptional modifications (PTMs) are often altered during neurogenesis when the brain is affected by disease or external stimuli and that the degree of histone modification is closely associated with the development of neurodegenerative diseases. This review aimed to show the regulatory role of histone modifications in neurogenesis and neurodegenerative diseases by discussing the changing patterns and functional significance of histone modifications, including histone methylation, acetylation, ubiquitination, phosphorylation, and lactylation. Finally, we explored the control of neurogenesis and the development of neurodegenerative diseases by artificially modulating histone modifications.

A functional schizophrenia-associated genetic variant near the TSNARE1 and ADGRB1 genes

Recent collaborative genome-wide association studies (GWAS) have identified >200 independent loci contributing to risk for schizophrenia (SCZ). The genes closest to these loci have diverse functions, supporting the potential involvement of multiple relevant biological processes, yet there is no direct evidence that individual variants are functional or directly linked to specific genes. Nevertheless, overlap with certain epigenetic marks suggest that most GWAS-implicated variants are regulatory. Based on the strength of association with SCZ and the presence of regulatory epigenetic marks, we chose one such variant near TSNARE1 and ADGRB1, rs4129585, to test for functional potential and assay differences that may drive the pathogenicity of the risk allele. We observed that the variant-containing sequence drives reporter expression in relevant neuronal populations in zebrafish. Next, we introduced each allele into human induced pluripotent cells and differentiated four isogenic clones homozygous for the risk allele and five clones homozygous for the non-risk allele into neural progenitor cells. Employing RNA sequencing, we found that the two alleles yield significant transcriptional differences in the expression of 109 genes at a false discovery rate (FDR) of <0.05 and 259 genes at a FDR of <0.1. We demonstrate that these genes are highly interconnected in pathways enriched for synaptic proteins, axon guidance, and regulation of synapse assembly. Exploration of genes near rs4129585 suggests that this variant does not regulate TSNARE1 transcripts, as previously thought, but may regulate the neighboring ADGRB1, a regulator of synaptogenesis. Our results suggest that rs4129585 is a functional common variant that functions in specific pathways likely involved in SCZ risk.

Systematic analysis of variants escaping nonsense-mediated decay uncovers candidate Mendelian diseases

Protein-truncating variants (PTVs) near the 3' end of genes may escape nonsense-mediated decay (NMD). PTVs in the NMD-escape region (PTVescs) can cause Mendelian disease but are difficult to interpret given their varying impact on protein function. Previously, PTVesc burden was assessed in an epilepsy cohort, but no large-scale analysis has systematically evaluated these variants in rare disease. We performed a retrospective analysis of 29,031 neurodevelopmental disorder (NDD) parent-offspring trios referred for clinical exome sequencing to identify PTVesc de novo mutations (DNMs). We identified 1,376 PTVesc DNMs and 133 genes that were significantly enriched (binomial p < 0.001). The PTVesc-enriched genes included those with PTVescs previously described to cause dominant Mendelian disease (e.g., SEMA6B, PPM1D, and DAGLA). We annotated ClinVar variants for PTVescs and identified 948 genes with at least one high-confidence pathogenic variant. Twenty-two known Mendelian PTVesc-enriched genes had no prior evidence of PTVesc-associated disease. We found 22 additional PTVesc-enriched genes that are not well established to be associated with Mendelian disease, several of which showed phenotypic similarity between individuals harboring PTVesc variants in the same gene. Four individuals with PTVesc mutations in RAB1A had similar phenotypes including NDD and spasticity. PTVesc mutations in IRF2BP1 were found in two individuals who each had severe immunodeficiency manifesting in NDD. Three individuals with PTVesc mutations in LDB1 all had NDD and multiple congenital anomalies. Using a large-scale, systematic analysis of DNMs, we extend the mutation spectrum for known Mendelian disease-associated genes and identify potentially novel disease-associated genes.

A Functional Schizophrenia-associated genetic variant near the TSNARE1 and ADGRB1 genes

Recent collaborative genome wide association studies (GWAS) have identified >200 independent loci contributing to risk for schizophrenia (SCZ). The genes closest to these loci have diverse functions, supporting the potential involvement of multiple relevant biological processes; yet there is no direct evidence that individual variants are functional or directly linked to specific genes. Nevertheless, overlap with certain epigenetic marks suggest that most GWAS-implicated variants are regulatory. Based on the strength of association with SCZ and the presence of regulatory epigenetic marks, we chose one such variant near TSNARE1 and ADGRB1, rs4129585, to test for functional potential and assay differences that may drive the pathogenicity of the risk allele. We observed that the variant-containing sequence drives reporter expression in relevant neuronal populations in zebrafish. Next, we introduced each allele into human induced pluripotent cells and differentiated 4 isogenic clones homozygous for the risk allele and 5 clones homozygous for the non-risk allele into neural precursor cells. Employing RNA-seq, we found that the two alleles yield significant transcriptional differences in the expression of 109 genes at FDR <0.05 and 259 genes at FDR <0.1. We demonstrate that these genes are highly interconnected in pathways enriched for synaptic proteins, axon guidance, and regulation of synapse assembly. Exploration of genes near rs4129585 suggests that this variant does not regulate TSNARE1 transcripts, as previously thought, but may regulate the neighboring ADGRB1, a regulator of synaptogenesis. Our results suggest that rs4129585 is a functional common variant that functions in specific pathways likely involved in SCZ risk.

Large-scale analysis of de novo mutations identifies risk genes for female infertility characterized by oocyte and early embryo defects

BACKGROUND

Oocyte maturation arrest and early embryonic arrest are important reproductive phenotypes resulting in female infertility and cause the recurrent failure of assisted reproductive technology (ART). However, the genetic etiologies of these female infertility-related phenotypes are poorly understood. Previous studies have mainly focused on inherited mutations based on large pedigrees or consanguineous patients. However, the role of de novo mutations (DNMs) in these phenotypes remains to be elucidated.

RESULTS

To decipher the role of DNMs in ART failure and female infertility with oocyte and embryo defects, we explore the landscape of DNMs in 473 infertile parent-child trios and identify a set of 481 confident DNMs distributed in 474 genes. Gene ontology analysis reveals that the identified genes with DNMs are enriched in signaling pathways associated with female reproductive processes such as meiosis, embryonic development, and reproductive structure development. We perform functional assays on the effects of DNMs in a representative gene Tubulin Alpha 4a (TUBA4A), which shows the most significant enrichment of DNMs in the infertile parent-child trios. DNMs in TUBA4A disrupt the normal assembly of the microtubule network in HeLa cells, and microinjection of DNM TUBA4A cRNAs causes abnormalities in mouse oocyte maturation or embryo development, suggesting the pathogenic role of these DNMs in TUBA4A.

CONCLUSIONS

Our findings suggest novel genetic insights that DNMs contribute to female infertility with oocyte and embryo defects. This study also provides potential genetic markers and facilitates the genetic diagnosis of recurrent ART failure and female infertility.

Stem Cell-Based Organoid Models of Neurodevelopmental Disorders

The past decade has seen an explosion in the identification of genetic causes of neurodevelopmental disorders, including Mendelian, de novo, and somatic factors. These discoveries provide opportunities to understand cellular and molecular mechanisms as well as potential gene-gene and gene-environment interactions to support novel therapies. Stem cell-based models, particularly human brain organoids, can capture disease-associated alleles in the context of the human genome, engineered to mirror disease-relevant aspects of cellular complexity and developmental timing. These models have brought key insights into neurodevelopmental disorders as diverse as microcephaly, autism, and focal epilepsy. However, intrinsic organoid-to-organoid variability, low levels of certain brain-resident cell types, and long culture times required to reach maturity can impede progress. Several recent advances incorporate specific morphogen gradients, mixtures of diverse brain cell types, and organoid engraftment into animal models. Together with nonhuman primate organoid comparisons, mechanisms of human neurodevelopmental disorders are emerging.

The fly homolog of SUPT16H, a gene associated with neurodevelopmental disorders, is required in a cell-autonomous fashion for cell survival

SUPT16H encodes the large subunit of the FAcilitate Chromatin Transcription (FACT) complex, which functions as a nucleosome organizer during transcription. We identified two individuals from unrelated families carrying de novo missense variants in SUPT16H. The probands exhibit global developmental delay, intellectual disability, epilepsy, facial dysmorphism and brain structural abnormalities. We used Drosophila to characterize two variants: p.T171I and p.G808R. Loss of the fly ortholog, dre4, causes lethality at an early developmental stage. RNAi-mediated knockdown of dre4 in either glia or neurons causes severely reduced eclosion and longevity. Tissue-specific knockdown of dre4 in the eye or wing leads to the loss of these tissues, whereas overexpression of SUPT16H has no dominant effect. Moreover, expression of the reference SUPT16H significantly rescues the loss-of-function phenotypes in the nervous system as well as wing and eye. In contrast, expression of SUPT16H p.T171I or p.G808R rescues the phenotypes poorly, indicating that the variants are partial loss-of-function alleles. While previous studies argued that the developmental arrest caused by loss of dre4 is due to impaired ecdysone production in the prothoracic gland, our data show that dre4 is required for proper cell growth and survival in multiple tissues in a cell-autonomous manner. Altogether, our data indicate that the de novo loss-of-function variants in SUPT16H are indeed associated with developmental and neurological defects observed in the probands.

The Weakened Interaction Between HECTD4 and GluN2B in Ischemic Stroke Promotes Calcium Overload and Brain Injury Through a Mechanism Involving the Decrease of GluN2B and MALT1 Ubiquitination

Glutamate receptor ionotropic NMDA 2B (GluN2B) plays an essential role in calcium overload during excitotoxicity. Reverse-phase nano-liquid chromatography-tandem mass spectrometry has revealed an interaction between GluN2B and HECT domain E3 ubiquitin protein ligase 4 (HECTD4), an E3 ubiquitin ligase highly expressed in the brain. As a potential substrate for HECTD4, mucosa-associated lymphoid tissue lymphoma translocation protein 1 (MALT1) acts as a scaffold with hydrolysis activity. This study explores the relationship between HECTD4, GluN2B, and MALT1, focusing on their role in brain injury in ischemic stroke. Rats were subjected to 2 h-ischemia followed by 24-h reperfusion to establish an ischemic stroke model. We observed the downregulation of HECTD4 and the upregulation of MALT1. Additionally, an increased GluN2B phosphorylation was concomitant with weakened interactions between HECTD4 and GluN2B, followed by decreased striatal-enriched protein phosphatase (STEP61). Knockdown of HECTD4 exacerbated hypoxia- or NMDA-induced injury in nerve cells coincident with a decrease in GluN2B and MALT1 ubiquitination, and an increase in GluN2B phosphorylation as well as an increase in intracellular calcium level, which were counteracted by MALT1 siRNA. Blockage of MALT1 with its inhibitor or siRNA reduced STEP61 degradation, accompanied by a decrease in GluN2B phosphorylation, intracellular calcium concentration, and brain cell injury, which were reversed by overexpression of MALT1. Based on these observations, we conclude that the downregulation of HECTD4 in ischemic stroke rat brain accounts for calcium overload and brain injury due to activating GluN2B directly and indirectly through a mechanism involving the reduced ubiquitination of GluN2B and MALT1, respectively.

Nine patients with KCNQ2-related neonatal seizures and functional studies of two missense variants

Mutations in KCNQ2 encoding for voltage-gated K channel subunits underlying the neuronal M-current have been associated with infantile-onset epileptic disorders. The clinical spectrum ranges from self-limited neonatal seizures to epileptic encephalopathy and delayed development. Mutations in KCNQ2 could be either gain- or loss-of-function which require different therapeutic approaches. To better understand genotype-phenotype correlation, more reports of patients and their mutations with elucidated molecular mechanism are needed. We studied 104 patients with infantile-onset pharmacoresistant epilepsy who underwent exome or genome sequencing. Nine patients with neonatal-onset seizures from unrelated families were found to harbor pathogenic or likely pathogenic variants in the KCNQ2 gene. The p.(N258K) was recently reported, and p. (G279D) has never been previously reported. Functional effect of p.(N258K) and p.(G279D) has never been previously studied. The cellular localization study demonstrated that the surface membrane expression of Kv7.2 carrying either variant was decreased. Whole-cell patch-clamp analyses revealed that both variants significantly impaired Kv7.2 M-current amplitude and density, conductance depolarizing shift in voltage dependence of activation, membrane resistance, and membrane time constant (Tau), indicating a loss-of-function in both the homotetrameric and heterotetrameric with Kv7.3 channels. In addition, both variants exerted dominant-negative effects in heterotetrameric with Kv7.3 channels. This study expands the mutational spectrum of KCNQ2- related epilepsy and their functional consequences provide insights into their pathomechanism.

Exposure to Antibiotics and Neurodevelopmental Disorders: Could Probiotics Modulate the Gut-Brain Axis?

In order to develop properly, the brain requires the intricate interconnection of genetic factors and pre-and postnatal environmental events. The gut-brain axis has recently raised considerable interest for its involvement in regulating the development and functioning of the brain. Consequently, alterations in the gut microbiota composition, due to antibiotic administration, could favor the onset of neurodevelopmental disorders. Literature data suggest that the modulation of gut microbiota is often altered in individuals affected by neurodevelopmental disorders. It has been shown in animal studies that metabolites released by an imbalanced gut-brain axis, leads to alterations in brain function and deficits in social behavior. Here, we report the potential effects of antibiotic administration, before and after birth, in relation to the risk of developing neurodevelopmental disorders. We also review the potential role of probiotics in treating gastrointestinal disorders associated with gut dysbiosis after antibiotic administration, and their possible effect in ameliorating neurodevelopmental disorder symptoms.

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.