Targeted sequencing and integrative analysis to prioritize candidate genes in neurodevelopmental disorders

MSL2 variants lead to a neurodevelopmental syndrome with lack of coordination, epilepsy, specific dysmorphisms, and a distinct episignature

Epigenetic dysregulation has emerged as an important etiological mechanism of neurodevelopmental disorders (NDDs). Pathogenic variation in epigenetic regulators can impair deposition of histone post-translational modifications leading to aberrant spatiotemporal gene expression during neurodevelopment. The male-specific lethal (MSL) complex is a prominent multi-subunit epigenetic regulator of gene expression and is responsible for histone 4 lysine 16 acetylation (H4K16ac). Using exome sequencing, here we identify a cohort of 25 individuals with heterozygous de novo variants in MSL complex member MSL2. MSL2 variants were associated with NDD phenotypes including global developmental delay, intellectual disability, hypotonia, and motor issues such as coordination problems, feeding difficulties, and gait disturbance. Dysmorphisms and behavioral and/or psychiatric conditions, including autism spectrum disorder, and to a lesser extent, seizures, connective tissue disease signs, sleep disturbance, vision problems, and other organ anomalies, were observed in affected individuals. As a molecular biomarker, a sensitive and specific DNA methylation episignature has been established. Induced pluripotent stem cells (iPSCs) derived from three members of our cohort exhibited reduced MSL2 levels. Remarkably, while NDD-associated variants in two other members of the MSL complex (MOF and MSL3) result in reduced H4K16ac, global H4K16ac levels are unchanged in iPSCs with MSL2 variants. Regardless, MSL2 variants altered the expression of MSL2 targets in iPSCs and upon their differentiation to early germ layers. Our study defines an MSL2-related disorder as an NDD with distinguishable clinical features, a specific blood DNA episignature, and a distinct, MSL2-specific molecular etiology compared to other MSL complex-related disorders.

Novel protein-truncating variants of a chromatin-modifying gene MSL2 in syndromic neurodevelopmental disorders

Numerous large scale genomic studies have uncovered rare but recurrent pathogenetic variants in a significant number of genes encoding epigenetic machinery in cases with neurodevelopmental disorders (NDD) especially autism spectrum disorder (ASD). These findings provide strong support for the functional importance of epigenetic regulators in neurodevelopment. After the clinical genomics evaluation of the patients using exome sequencing, we have identified, three novel protein-truncating variants (PTVs) in the MSL2 gene (OMIM: 614802) which encodes a chromatin modifying enzyme. MSL2 modifies chromatin through both mono-ubiquitination of histone 2B on lysine 34 (K34) and acetylation of histone H4 on lysine 16 (K16). We reported first time the detailed clinical features associated with 3 MSL2 PTVs. There are 15 PTVs (13 de novo) reported from the large genomics studies (12 cases) or ClinVar (3 cases) of NDD, ASD, and developmental disorders (DD) but the specific clinical features for these cases are not described. Taken together, our descriptions of dysmorphic face and other features support the causal role of MSL2 in a likely syndromic neurodevelopmental disorder and add MSL2 to a growing list of epigenetic genes implicated in ASD.

Shared and divergent contribution of vitamin A and oxytocin to the aetiology of autism spectrum disorder

Rare genetic variations contribute to the heterogeneity of autism spectrum disorder (ASD) and the responses to various interventions for ASD probands. However, the associated molecular underpinnings remain unclear. Herein, we estimated the association between rare genetic variations in 410 vitamin A (VA)-related genes (VARGs) and ASD aetiology using publicly available de novo mutations (DNMs), rare inherited variants, and copy number variations (CNVs) from about 50,000 ASD probands and 20,000 normal controls (discovery and validation cohorts). Additionally, given the functional relevance of VA and oxytocin, we systematically compared the similarities and differences between VA and oxytocin with respect to ASD aetiology and evaluated their potential for clinical applications. Functional DNMs and pathogenic CNVs in VARGs contributed to ASD pathogenesis in the discovery and validation cohorts. Additionally, 324 potential VA-related biomarkers were identified, 243 of which were shared with previously identified oxytocin-related biomarkers, while 81 were unique VA biomarkers. Moreover, multivariable logistic regression analysis revealed that both VA- and oxytocin-related biomarkers were able to predict ASD aetiology for individuals carrying functional DNM in corresponding biomarkers with an average precision of 0.94. As well as, convergent and divergent functions were also identified between VA- and oxytocin-related biomarkers. The findings of this study provide a basis for future studies aimed at understanding the pathophysiological mechanisms underlying ASD while also defining a set of potential molecular biomarkers for adjuvant diagnosis and intervention in ASD.

Relating enhancer genetic variation across mammals to complex phenotypes using machine learning

Protein-coding differences between species often fail to explain phenotypic diversity, suggesting the involvement of genomic elements that regulate gene expression such as enhancers. Identifying associations between enhancers and phenotypes is challenging because enhancer activity can be tissue-dependent and functionally conserved despite low sequence conservation. We developed the Tissue-Aware Conservation Inference Toolkit (TACIT) to associate candidate enhancers with species' phenotypes using predictions from machine learning models trained on specific tissues. Applying TACIT to associate motor cortex and parvalbumin-positive interneuron enhancers with neurological phenotypes revealed dozens of enhancer-phenotype associations, including brain size-associated enhancers that interact with genes implicated in microcephaly or macrocephaly. TACIT provides a foundation for identifying enhancers associated with the evolution of any convergently evolved phenotype in any large group of species with aligned genomes.

A rigorous in silico genomic interrogation at 1p13.3 reveals 16 autosomal dominant candidate genes in syndromic neurodevelopmental disorders

Genome-wide chromosomal microarray is extensively used to detect copy number variations (CNVs), which can diagnose microdeletion and microduplication syndromes. These small unbalanced chromosomal structural rearrangements ranging from 1 kb to 10 Mb comprise up to 15% of human mutations leading to monogenic or contiguous genomic disorders. Albeit rare, CNVs at 1p13.3 cause a variety of neurodevelopmental disorders (NDDs) including development delay (DD), intellectual disability (ID), autism, epilepsy, and craniofacial anomalies (CFA). Most of the 1p13.3 CNV cases reported in the pre-microarray era encompassed a large number of genes and lacked the demarcating genomic coordinates, hampering the discovery of positional candidate genes within the boundaries. In this study, we present four subjects with 1p13.3 microdeletions displaying DD, ID, autism, epilepsy, and CFA. In silico comparative genomic mapping with three previously reported subjects with CNVs and 22 unreported DECIPHER CNV cases has resulted in the identification of four different sub-genomic loci harboring five positional candidate genes for DD, ID, and CFA at 1p13.3. Most of these genes have pathogenic variants reported, and their interacting genes are involved in NDDs. RT-qPCR in various human tissues revealed a high expression pattern in the brain and fetal brain, supporting their functional roles in NDDs. Interrogation of variant databases and interacting protein partners led to the identification of another set of 11 potential candidate genes, which might have been dysregulated by the position effect of these CNVs at 1p13.3. Our studies define 1p13.3 as a genomic region harboring 16 NDD candidate genes and underscore the critical roles of small CNVs in in silico comparative genomic mapping for disease gene discovery. Our candidate genes will help accelerate the isolation of pathogenic heterozygous variants from exome/genome sequencing (ES/GS) databases.

Integrative analysis prioritised oxytocin-related biomarkers associated with the aetiology of autism spectrum disorder

Background

Autism spectrum disorder (ASD) is a neurodevelopmental disorder with high phenotypic and genetic heterogeneity. The common variants of specific oxytocin-related genes (OTRGs), particularly OXTR, are associated with the aetiology of ASD. The contribution of rare genetic variations in OTRGs to ASD aetiology remains unclear.

Methods

We catalogued publicly available de novo mutations (DNMs) [from 6,511 patients with ASD and 3,391 controls], rare inherited variants (RIVs) [from 1,786 patients with ASD and 1,786 controls], and both de novo copy number variations (dnCNVs) and inherited CNVs (ihCNVs) [from 15,581 patients with ASD and 6,017 controls] in 963 curated OTRGs to explore their contribution to ASD pathology, respectively. Finally, a combined model was designed to prioritise the contribution of each gene to ASD aetiology by integrating DNMs and CNVs.

Findings

The rare genetic variations of OTRGs were significantly associated with ASD aetiology, in the order of dnCNVs > ihCNVs > DNMs. Furthermore, 172 OTRGs and their connected 286 ASD core genes were prioritised to positively contribute to ASD aetiology, including top-ranked MAPK3. Probands carrying rare disruptive variations in these genes were estimated to account for 10∼11% of all ASD probands.

Interpretation

Our findings suggest that rare disruptive variations in 172 OTRGs and their connected 286 ASD core genes are associated with ASD aetiology and may be potential biomarkers predicting the effects of oxytocin treatment.

Funding

Guangdong Key Project, National Natural Science Foundation of China, Key Laboratory of Clinical Laboratory Diagnosis and Translational Research of Zhejiang Province.