The Role of Microtubule Associated Serine/Threonine Kinase 3 Variants in Neurodevelopmental Diseases: Genotype-Phenotype Association

Objective: To prove microtubule associated serine/threonine kinase 3 (MAST3) gene is associated with neurodevelopmental diseases (NDD) and the genotype-phenotype correlation.Methods: Trio exome sequencing (trio ES) was performed on four NDD trios. Bioinformatic analysis was conducted based on large-scale genome sequencing data and human brain transcriptomic data. Further in vivo zebrafish studies were performed.Results: In our study, we identified four de novo MAST3 variants (NM_015016.1: c.302C > T:p.Ser101Phe; c.311C > T:p.Ser104Leu; c.1543G > A:p.Gly515Ser; and c.1547T > C:p.Leu516Pro) in four patients with developmental and epileptic encephalopathy (DEE) separately. Clinical heterogeneities were observed in patients carrying variants in domain of unknown function (DUF) and serine-threonine kinase (STK) domain separately. Using the published large-scale exome sequencing data, higher CADD scores of missense variants in DUF domain were found in NDD cohort compared with gnomAD database. In addition, we obtained an excess of missense variants in DUF domain when compared autistic spectrum disorder (ASD) cohort with gnomAD database, similarly an excess of missense variants in STK domain when compared DEE cohort with gnomAD database. Based on Brainspan datasets, we showed that MAST3 expression was significantly upregulated in ASD and DEE-related brain regions and was functionally linked with DEE genes. In zebrafish model, abnormal morphology of central nervous system was observed in mast3a/b crispants.Conclusion: Our results support the possibility that MAST3 is a novel gene associated with NDD which could expand the genetic spectrum for NDD. The genotype-phenotype correlation may contribute to future genetic counseling.

Novel biallelic variants expand the SLC5A6-related phenotypic spectrum

The sodium (Na+):multivitamin transporter (SMVT), encoded by SLC5A6, belongs to the sodium:solute symporter family and is required for the Na+-dependent uptake of biotin (vitamin B7), pantothenic acid (vitamin B5), the vitamin-like substance α-lipoic acid, and iodide. Compound heterozygous SLC5A6 variants have been reported in individuals with variable multisystemic disorder, including failure to thrive, developmental delay, seizures, cerebral palsy, brain atrophy, gastrointestinal problems, immunodeficiency, and/or osteopenia. We expand the phenotypic spectrum associated with biallelic SLC5A6 variants affecting function by reporting five individuals from three families with motor neuropathies. We identified the homozygous variant c.1285 A > G [p.(Ser429Gly)] in three affected siblings and a simplex patient and the maternally inherited c.280 C > T [p.(Arg94*)] variant and the paternally inherited c.485 A > G [p.(Tyr162Cys)] variant in the simplex patient of the third family. Both missense variants were predicted to affect function by in silico tools. 3D homology modeling of the human SMVT revealed 13 transmembrane helices (TMs) and Tyr162 and Ser429 to be located at the cytoplasmic facing region of TM4 and within TM11, respectively. The SLC5A6 missense variants p.(Tyr162Cys) and p.(Ser429Gly) did not affect plasma membrane localization of the ectopically expressed multivitamin transporter suggesting reduced but not abolished function, such as lower catalytic activity. Targeted therapeutic intervention yielded clinical improvement in four of the five patients. Early molecular diagnosis by exome sequencing is essential for timely replacement therapy in affected individuals.

Missense variants in the TNFA epitopes and their effects on interaction with therapeutic antibodies—in silico analysis

Background Tumor necrosis factor alpha (TNFA) is an important cytokine that influences multiple biological processes. It is one of the key mediators of acute and chronic systemic inflammatory reactions and plays a central role in several autoimmune diseases. A number of approved monoclonal antibodies (mAbs) are widely used to subside these autoimmune diseases. However, there is a high rate of non-responsiveness to treatments with these mAbs. Therefore, it is important to be able to predict responses of the patients in an individualistic manner to these therapeutic antibodies before administration. In the present study, we used in silico tools to explore the effects of missense variants in the respective epitopes of four therapeutic anti-TNFA mAbs—adalimumab (ADA), certolizumab pegol (CZP), golimumab (GLM), and infliximab (IFX)—on their interactions with TNFA. Results The binding affinities of CZP and ADA to corresponding epitopes appear to be reduced by four (TNFAR131Q, TNFAE135G, TNFAR138Q, and TNFAR138W) and two (TNFAG66C and TNFAG66S) variants, respectively. The binding of GLM and IFX appears to be affected by TNFAR141S and TNFAR138W, respectively. TNFAG66C and TNFAG66S may be associated with autoimmune diseases, whereas TNFAE135G, TNFAR138W, and TNFAR141S may be pathogenic per se. Conclusion These variants may contribute to the observed inter-individual variability in response to anti-TNFA mAbs treatments and be used as markers to predict responses, and thus optimize therapeutic benefits to the patients.

In-silico phenotype prediction by normal mode variant analysis in TUBB4A-related disease

TUBB4A-associated disorder is a rare condition affecting the central nervous system. It displays a wide phenotypic spectrum, ranging from isolated late-onset torsion dystonia to a severe early-onset disease with developmental delay, neurological deficits, and atrophy of the basal ganglia and cerebellum, therefore complicating variant interpretation and phenotype prediction in patients carrying TUBB4A variants. We applied entropy-based normal mode analysis (NMA) to investigate genotype–phenotype correlations in TUBB4A-releated disease and to develop an in-silico approach to assist in variant interpretation and phenotype prediction in this disorder. Variants included in our analysis were those reported prior to the conclusion of data collection for this study in October 2019. All TUBB4A pathogenic missense variants reported in ClinVar and Pubmed, for which associated clinical information was available, and all benign/likely benign TUBB4A missense variants reported in ClinVar, were included in the analysis. Pathogenic variants were divided into five phenotypic subgroups. In-silico point mutagenesis in the wild-type modeled protein structure was performed for each variant. Wild-type and mutated structures were analyzed by coarse-grained NMA to quantify protein stability as entropy difference value (ΔG) for each variant. Pairwise ΔG differences between all variant pairs in each structural cluster were calculated and clustered into dendrograms. Our search yielded 41 TUBB4A pathogenic variants in 126 patients, divided into 11 partially overlapping structural clusters across the TUBB4A protein. ΔG-based cluster analysis of the NMA results revealed a continuum of genotype–phenotype correlation across each structural cluster, as well as in transition areas of partially overlapping structural clusters. Benign/likely benign variants were integrated into the genotype–phenotype continuum as expected and were clearly separated from pathogenic variants. We conclude that our results support the incorporation of the NMA-based approach used in this study in the interpretation of variant pathogenicity and phenotype prediction in TUBB4A-related disease. Moreover, our results suggest that NMA may be of value in variant interpretation in additional monogenic conditions.

L-Serine Treatment is Associated with Improvements in Behavior, EEG, and Seizure Frequency in Individuals with GRIN-Related Disorders Due to Null Variants

Pathogenic missense variants in GRIN2A and GRIN2B may result in gain or loss of function (GoF/LoF) of the N-methyl-D-aspartate receptor (NMDAR). This observation gave rise to the hypothesis of successfully treating GRIN-related disorders due to LoF variants with co-agonists of the NMDAR. In this respect, we describe a retrospectively collected series of ten individuals with GRIN2A- or GRIN2B-related disorders who were treated with L-serine, each within an independent n-of-1 trial. Our cohort comprises one individual with a LoF missense variant with clinical improvements confirming the above hypothesis and replicating a previous n-of-1 trial. A second individual with a GoF missense variant was erroneously treated with L-serine and experienced immediate temporary behavioral deterioration further supporting the supposed functional pathomechanism. Eight additional individuals with null variants (that had been interpreted as loss-of-function variants despite not being missense) again showed clinical improvements. Among all nine individuals with LoF missense or null variants, L-serine treatment was associated with improvements in behavior in eight (89%), in development in four (44%), and/or in EEG or seizure frequency in four (44%). None of these nine individuals experienced side effects or adverse findings in the context of L-serine treatment. In summary, we describe the first evidence that L-serine treatment may not only be associated with clinical improvements in GRIN-related disorders due to LoF missense but particularly also null variants.

MvPPT: a highly efficient and sensitive pathogenicity prediction tool for missense variants

Next generation sequencing technologies both boost the discovery of variants in the human genome and exacerbate the challenges of pathogenic variant identification. In this study, we developed mvPPT (Pathogenicity Prediction Tool for missense variants), a highly sensitive and accurate missense variant classifier based on gradient boosting. MvPPT adopts high-confidence training sets with a wide spectrum of variant profiles, and extracts three categories of features, including scores from existing prediction tools, allele, amino acid and genotype frequencies, and genomic context. Compared with established predictors, mvPPT achieved superior performance in all test sets, regardless of data source. In addition, our study also provides guidance for training set and feature selection strategies, as well as reveals highly relevant features, which may further provide biological insights of variant pathogenicity.

Missense variants causing Wiedemann-Steiner syndrome preferentially occur in the KMT2A-CXXC domain and are accurately classified using AlphaFold2

Wiedemann-Steiner syndrome (WSS) is a neurodevelopmental disorder caused by de novo variants in KMT2A, which encodes a multi–domain histone methyltransferase. To gain insight into the currently unknown pathogenesis of WSS, we examined the spatial distribution of likely WSS–causing variants across the 15 different domains of KMT2A. Compared to variants in healthy controls, WSS variants exhibit a 64.1–fold overrepresentation within the CXXC domain – which mediates binding to unmethylated CpGs – suggesting a major role for this domain in mediating the phenotype. In contrast, we find no significant overrepresentation within the catalytic SET domain. Corroborating these results, we find that hippocampal neurons from Kmt2a–deficient mice demonstrate disrupted H3K4me1 preferentially at CpG-rich regions, but this has no systematic impact on gene expression. Motivated by these results, we combine accurate prediction of the CXXC domain structure by AlphaFold2 with prior biological knowledge to develop a classification scheme for missense variants in the CXXC domain. Our classifier achieved 96.0% positive and 92.3% negative predictive value on a hold–out test set. This classification performance enabled us to subsequently perform an in silico saturation mutagenesis and classify a total of 445 variants according to their functional effects. Our results yield a novel insight into the mechanistic basis of WSS and provide an example of how AlphaFold2 can contribute to the in silico characterization of variant effects with very high accuracy, establishing a paradigm potentially applicable to many other Mendelian disorders.

Missense variants in ANKRD11 cause KBG syndrome by impairment of stability or transcriptional activity of the encoded protein

Purpose: Although haploinsufficiency of ANKRD11 is among the most common genetic causes of neurodevelopmental disorders, the role of rare ANKRD11 missense variation remains unclear. We characterized the clinical, molecular and functional spectra of ANKRD11 missense variants. Methods: We collected clinical information of individuals with ANKRD11 missense variants and evaluated phenotypic fit to KBG syndrome. We assessed pathogenicity of variants by in silico analyses and cell-based experiments. Results: We identified 29 individuals with (mostly de novo) ANKRD11 missense variants, who presented with syndromic neurodevelopmental disorders and were phenotypically similar to individuals with KBG syndrome caused by ANKRD11 protein truncating variants or 16q24.3 microdeletions. Missense variants significantly clustered in Repression Domain 2. Cellularly, most variants caused reduced ANKRD11 stability. One variant resulted in decreased proteasome degradation and loss of ANKRD11 transcriptional activity. Conclusion: Our study indicates that pathogenic heterozygous missense variants in ANKRD11 cause the clinically recognizable KBG syndrome. Disrupted transrepression capacity and reduced protein stability each independently lead to ANKRD11 loss-of-function, consistent with haploinsufficiency. This highlights the diagnostic relevance of ANKRD11 missense variants, but also poses diagnostic challenges, as the KBG-associated phenotype may be mild and inherited pathogenic ANKRD11 (missense) variants are increasingly observed, warranting stringent variant classification and careful phenotyping.

Rare coding variation illuminates the allelic architecture, risk genes, cellular expression patterns, and phenotypic context of autism

Individuals with autism spectrum disorder (ASD) or related neurodevelopmental disorders (NDDs) often carry disruptive mutations in genes that are depleted of functional variation in the broader population. We build upon this observation and exome sequencing from 154,842 individuals to explore the allelic diversity of rare protein-coding variation contributing risk for ASD and related NDDs. Using an integrative statistical model, we jointly analyzed rare protein-truncating variants (PTVs), damaging missense variants, and copy number variants (CNVs) derived from exome sequencing of 63,237 individuals from ASD cohorts. We discovered 71 genes associated with ASD at a false discovery rate (FDR) ≤ 0.001, a threshold approximately equivalent to exome-wide significance, and 183 genes at FDR ≤ 0.05. Associations were predominantly driven by de novo PTVs, damaging missense variants, and CNVs: 57.4%, 21.2%, and 8.32% of evidence, respectively. Though fewer in number, CNVs conferred greater relative risk than PTVs, and repeat-mediated de novo CNVs exhibited strong maternal bias in parent-of-origin (e.g., 92.3% of 16p11.2 CNVs), whereas all other CNVs showed a paternal bias. To explore how genes associated with ASD and NDD overlap or differ, we analyzed our ASD cohort alongside a developmental delay (DD) cohort from the deciphering developmental disorders study (DDD; n=91,605 samples). We first reanalyzed the DDD dataset using the same models as the ASD cohorts, then performed joint analyses of both cohorts and identified 373 genes contributing to NDD risk at FDR ≤ 0.001 and 662 NDD risk genes at FDR ≤ 0.05. Of these NDD risk genes, 54 genes (125 genes at FDR ≤ 0.05) were unique to the joint analyses and not significant in either cohort alone. Our results confirm overlap of most ASD and DD risk genes, although many differ significantly in frequency of mutation. Analyses of single-cell transcriptome datasets showed that genes associated predominantly with DD were strongly enriched for earlier neurodevelopmental cell types, whereas genes displaying stronger evidence for association in ASD cohorts were more enriched for maturing neurons. The ASD risk genes were also enriched for genes associated with schizophrenia from a separate rare coding variant analysis of 121,570 individuals, emphasizing that these neuropsychiatric disorders share common pathways to risk.

Variant-specific effects define the phenotypic spectrum of HNRNPH2-associated neurodevelopmental disorders in males

Bain type of X-linked syndromic intellectual developmental disorder, caused by pathogenic missense variants in HRNRPH2, was initially described in six female individuals affected by moderate-to-severe neurodevelopmental delay. Although it was initially postulated that the condition would not be compatible with life in males, several affected male individuals harboring pathogenic variants in HNRNPH2 have since been documented. However, functional in-vitro analyses of identified variants have not been performed and, therefore, possible genotype–phenotype correlations remain elusive. Here, we present eight male individuals, including a pair of monozygotic twins, harboring pathogenic or likely pathogenic HNRNPH2 variants. Notably, we present the first individuals harboring nonsense or frameshift variants who, similarly to an individual harboring a de novo p.(Arg29Cys) variant within the first quasi-RNA-recognition motif (qRRM), displayed mild developmental delay, and developed mostly autistic features and/or psychiatric co-morbidities. Additionally, we present two individuals harboring a recurrent de novo p.(Arg114Trp), within the second qRRM, who had a severe neurodevelopmental delay with seizures. Functional characterization of the three most common HNRNPH2 missense variants revealed dysfunctional nucleocytoplasmic shuttling of proteins harboring the p.(Arg206Gln) and p.(Pro209Leu) variants, located within the nuclear localization signal, whereas proteins with p.(Arg114Trp) showed reduced interaction with members of the large assembly of splicing regulators (LASR). Moreover, RNA-sequencing of primary fibroblasts of the individual harboring the p.(Arg114Trp) revealed substantial alterations in the regulation of alternative splicing along with global transcriptome changes. Thus, we further expand the clinical and variant spectrum in HNRNPH2-associated disease in males and provide novel molecular insights suggesting the disorder to be a spliceopathy on the molecular level.