Using the ACMG/AMP framework to capture evidence related to predicted and observed impact on splicing: Recommendations from the ClinGen SVI Splicing Subgroup.

The American College of Medical Genetics and Genomics (ACMG)/Association for Molecular Pathology (AMP) framework for classifying variants uses six evidence categories related to the splicing potential of variants: PVS1, PS3, PP3, BS3, BP4, and BP7. However, the lack of guidance on how to apply such codes has contributed to variation in the specifications developed by different Clinical Genome Resource (ClinGen) Variant Curation Expert Panels. The ClinGen Sequence Variant Interpretation Splicing Subgroup was established to refine recommendations for applying ACMG/AMP codes relating to splicing data and computational predictions. We utilized empirically derived splicing evidence to (1) determine the evidence weighting of splicing-related data and appropriate criteria code selection for general use, (2) outline a process for integrating splicing-related considerations when developing a gene-specific PVS1 decision tree, and (3) exemplify methodology to calibrate splice prediction tools. We propose repurposing the PVS1_Strength code to capture splicing assay data that provide experimental evidence for variants resulting in RNA transcript(s) with loss of function. Conversely, BP7 may be used to capture RNA results demonstrating no splicing impact for intronic and synonymous variants. We propose that the PS3/BS3 codes are applied only for well-established assays that measure functional impact not directly captured by RNA-splicing assays. We recommend the application of PS1 based on similarity of predicted RNA-splicing effects for a variant under assessment in comparison with a known pathogenic variant. The recommendations and approaches for consideration and evaluation of RNA-assay evidence described aim to help standardize variant pathogenicity classification processes when interpreting splicing-based evidence.

APPLICATION OF THE ACMG/AMP FRAMEWORK TO CAPTURE EVIDENCE RELEVANT TO PREDICTED AND OBSERVED IMPACT ON SPLICING: RECOMMENDATIONS FROM THE CLINGEN SVI SPLICING SUBGROUP.

The American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) framework for classifying variants uses six evidence categories related to the splicing potential of variants: PVS1 (null variant in a gene where loss-of-function is the mechanism of disease), PS3 (functional assays show damaging effect on splicing), PP3 (computational evidence supports a splicing effect), BS3 (functional assays show no damaging effect on splicing), BP4 (computational evidence suggests no splicing impact), and BP7 (silent change with no predicted impact on splicing). However, the lack of guidance on how to apply such codes has contributed to variation in the specifications developed by different Clinical Genome Resource (ClinGen) Variant Curation Expert Panels. The ClinGen Sequence Variant Interpretation (SVI) Splicing Subgroup was established to refine recommendations for applying ACMG/AMP codes relating to splicing data and computational predictions. Our study utilised empirically derived splicing evidence to: 1) determine the evidence weighting of splicing-related data and appropriate criteria code selection for general use, 2) outline a process for integrating splicing-related considerations when developing a gene-specific PVS1 decision tree, and 3) exemplify methodology to calibrate bioinformatic splice prediction tools. We propose repurposing of the PVS1_Strength code to capture splicing assay data that provide experimental evidence for variants resulting in RNA transcript(s) with loss of function. Conversely BP7 may be used to capture RNA results demonstrating no impact on splicing for both intronic and synonymous variants, and for missense variants if protein functional impact has been excluded. Furthermore, we propose that the PS3 and BS3 codes are applied only for well-established assays that measure functional impact that is not directly captured by RNA splicing assays. We recommend the application of PS1 based on similarity of predicted RNA splicing effects for a variant under assessment in comparison to a known Pathogenic variant. The recommendations and approaches for consideration and evaluation of RNA assay evidence described aim to help standardise variant pathogenicity classification processes and result in greater consistency when interpreting splicing-based evidence.

Phenotypic continuum of NFU1-related disorders.

Bi-allelic variants in Iron-Sulfur Cluster Scaffold (NFU1) have previously been associated with multiple mitochondrial dysfunctions syndrome 1 (MMDS1) characterized by early-onset rapidly fatal leukoencephalopathy. We report 19 affected individuals from 10 independent families with ultra-rare bi-allelic NFU1 missense variants associated with a spectrum of early-onset pure to complex hereditary spastic paraplegia (HSP) phenotype with a longer survival (16/19) on one end and neurodevelopmental delay with severe hypotonia (3/19) on the other. Reversible or irreversible neurological decompensation after a febrile illness was common in the cohort, and there were invariable white matter abnormalities on neuroimaging. The study suggests that MMDS1 and HSP could be the two ends of the NFU1-related phenotypic continuum.

The phenotypic spectrum and genotype-phenotype correlations in 106 patients with variants in major autism gene CHD8.

CHD8, a major autism gene, functions in chromatin remodelling and has various roles involving several biological pathways. Therefore, unsurprisingly, previous studies have shown that intellectual developmental disorder with autism and macrocephaly (IDDAM), the syndrome caused by pathogenic variants in CHD8, consists of a broad range of phenotypic abnormalities. We collected and reviewed 106 individuals with IDDAM, including 36 individuals not previously published, thus enabling thorough genotype-phenotype analyses, involving the CHD8 mutation spectrum, characterization of the CHD8 DNA methylation episignature, and the systematic analysis of phenotypes collected in Human Phenotype Ontology (HPO). We identified 29 unique nonsense, 25 frameshift, 24 missense, and 12 splice site variants. Furthermore, two unique inframe deletions, one larger deletion (exons 26-28), and one translocation were observed. Methylation analysis was performed for 13 patients, 11 of which showed the previously established episignature for IDDAM (85%) associated with CHD8 haploinsufficiency, one analysis was inconclusive, and one showing a possible gain-of-function signature instead of the expected haploinsufficiency signature was observed. Consistent with previous studies, phenotypical abnormalities affected multiple organ systems. Many neurological abnormalities, like intellectual disability (68%) and hypotonia (29%) were observed, as well as a wide variety of behavioural abnormalities (88%). Most frequently observed behavioural problems included autism spectrum disorder (76%), short attention span (32%), abnormal social behaviour (31%), sleep disturbance (29%) and impaired social interactions (28%). Furthermore, abnormalities in the digestive (53%), musculoskeletal (79%) and genitourinary systems (18%) were noted. Although no significant difference in severity was observed between males and females, individuals with a missense variant were less severely affected. Our study provides an extensive review of all phenotypic abnormalities in patients with IDDAM and provides clinical recommendations, which will be of significant value to individuals with a pathogenic variant in CHD8, their families, and clinicians as it gives a more refined insight into the clinical and molecular spectrum of IDDAM, which is essential for accurate care and counselling.

High-throughput evaluation of epilepsy-associated KCNQ2 variants reveals functional and pharmacological heterogeneity.

Hundreds of genetic variants in KCNQ2 encoding the voltage-gated potassium channel KV7.2 are associated with early onset epilepsy and/or developmental disability, but the functional consequences of most variants are unknown. Absent functional annotation for KCNQ2 variants hinders identification of individuals who may benefit from emerging precision therapies. We employed automated patch clamp recordings to assess at, to our knowledge, an unprecedented scale the functional and pharmacological properties of 79 missense and 2 inframe deletion KCNQ2 variants. Among the variants we studied were 18 known pathogenic variants, 24 mostly rare population variants, and 39 disease-associated variants with unclear functional effects. We analyzed electrophysiological data recorded from 9,480 cells. The functional properties of 18 known pathogenic variants largely matched previously published results and validated automated patch clamp for this purpose. Unlike rare population variants, most disease-associated KCNQ2 variants exhibited prominent loss-of-function with dominant-negative effects, providing strong evidence in support of pathogenicity. All variants responded to retigabine, although there were substantial differences in maximal responses. Our study demonstrated that dominant-negative loss-of-function is a common mechanism associated with missense KCNQ2 variants. Importantly, we observed genotype-dependent differences in the response of KCNQ2 variants to retigabine, a proposed precision therapy for KCNQ2 developmental and epileptic encephalopathy.

A catalogue of new incidence estimates of monogenic neurodevelopmental disorders caused by de novo variants.

A large fraction of rare and severe neurodevelopmental disorders are caused by sporadic de novo variants. Epidemiological disease estimates are not available for the vast majority of these de novo monogenic neurodevelopmental disorders because of phenotypic heterogeneity and the absence of large-scale genomic screens. Yet, knowledge of disease incidence is important for clinicians and researchers to guide health policy planning. Here, we adjusted a statistical method based on genetic data to predict, for the first time, the incidences of 101 known de novo variant-associated neurodevelopmental disorders as well as 3106 putative monogenic disorders. Two corroboration analyses supported the validity of the calculated estimates. First, greater predicted gene-disorder incidences positively correlated with larger numbers of pathogenic variants collected from patient variant databases (Kendall's τ = 0.093, P-value = 6.9 × 10-6). Second, for six of seven (86%) de novo variant associated monogenic disorders for which epidemiological estimates were available (SCN1A, SLC2A1, SALL1, TBX5, KCNQ2, and CDKL5), the predicted incidence estimates matched the reported estimates. We conclude that in the absence of epidemiological data, our catalogue of 3207 incidence estimates for disorders caused by de novo variants can guide patient advocacy groups, clinicians, researchers, and policymakers in strategic decision-making.

Genetic variant pathogenicity prediction trained using disease-specific clinical sequencing data sets.

Recent advances in DNA sequencing have expanded our understanding of the molecular basis of genetic disorders and increased the utilization of clinical genomic tests. Given the paucity of evidence to accurately classify each variant and the difficulty of experimentally evaluating its clinical significance, a large number of variants generated by clinical tests are reported as variants of unknown clinical significance. Population-scale variant databases can improve clinical interpretation. Specifically, pathogenicity prediction for novel missense variants can use features describing regional variant constraint. Constrained genomic regions are those that have an unusually low variant count in the general population. Computational methods have been introduced to capture these regions and incorporate them into pathogenicity classifiers, but these methods have yet to be compared on an independent clinical variant data set. Here, we introduce one variant data set derived from clinical sequencing panels and use it to compare the ability of different genomic constraint metrics to determine missense variant pathogenicity. This data set is compiled from 17,071 patients surveyed with clinical genomic sequencing for cardiomyopathy, epilepsy, or RASopathies. We further use this data set to demonstrate the necessity of disease-specific classifiers and to train PathoPredictor, a disease-specific ensemble classifier of pathogenicity based on regional constraint and variant-level features. PathoPredictor achieves an average precision >90% for variants from all 99 tested disease genes while approaching 100% accuracy for some genes. The accumulation of larger clinical variant training data sets can significantly enhance their performance in a disease- and gene-specific manner.

Diagnostic outcomes for genetic testing of 70 genes in 8565 patients with epilepsy and neurodevelopmental disorders.

OBJECTIVE: We evaluated >8500 consecutive, unselected patients with epilepsy and neurodevelopmental disorders who underwent multigene panel testing to determine the average age at molecular diagnosis and diagnostic yield of 70 genes. METHODS: We reviewed molecular test results for 70 genes known to cause epilepsy and neurodevelopmental disorders using next generation sequencing (NGS) and exon-level array comparative genomic hybridization (aCGH). A positive result was defined as the presence of 1 or 2 pathogenic or likely pathogenic (P/LP) variants in a single gene, depending on the mode of inheritance of the associated disorder. RESULTS: Overall, 22 genes were found to have a high yield of positive findings by genetic testing, with SCN1A and KCNQ2 accounting for the greatest number of positive findings. In contrast, there were no positive findings in 16 genes. Most of the P/LP variants were sequence changes identified by NGS (90.9%), whereas ~9% were gross deletions or duplications detected by exon-level aCGH. The mean age of molecular diagnosis for the cohort was 5 years, 8 months (ranging from 1 week to 47 years). Recurrent P/LP variants were observed in 14 distinct genes, most commonly in MECP2, KCNQ2, SCN1A, SCN2A, STXBP1, and PRRT2. Parental testing was performed in >30% of positive cases. All variants identified in CDKL5, STXBP1, SCN8A, GABRA1, and FOXG1 were de novo, whereas 85.7% of variants in PRRT2 were inherited. SIGNIFICANCE: Using a combined approach of NGS and exon-level aCGH, testing identified a genetic etiology in 15.4% of patients in this cohort and revealed the age at molecular diagnosis for patients. Our study highlights both high- and low-yield genes associated with epilepsy and neurodevelopmental disorders, indicating which genes may be considered for molecular diagnostic testing.

High frequency of mosaic pathogenic variants in genes causing epilepsy-related neurodevelopmental disorders.

PurposeMosaicism probably represents an underreported cause of genetic disorders due to detection challenges during routine molecular diagnostics. The purpose of this study was to evaluate the frequency of mosaicism detected by next-generation sequencing in genes associated with epilepsy-related neurodevelopmental disorders.MethodsWe conducted a retrospective analysis of 893 probands with epilepsy who had a multigene epilepsy panel or whole-exome sequencing performed in a clinical diagnostic laboratory and were positive for a pathogenic or likely pathogenic variant in one of nine genes (CDKL5, GABRA1, GABRG2, GRIN2B, KCNQ2, MECP2, PCDH19, SCN1A, or SCN2A). Parental results were available for 395 of these probands.ResultsMosaicism was most common in the CDKL5, PCDH19, SCN2A, and SCN1A genes. Mosaicism was observed in GABRA1, GABRG2, and GRIN2B, which previously have not been reported to have mosaicism, and also in KCNQ2 and MECP2. Parental mosaicism was observed for pathogenic variants in multiple genes including KCNQ2, MECP2, SCN1A, and SCN2A.ConclusionMosaic pathogenic variants were identified frequently in nine genes associated with various neurological conditions. Given the potential clinical ramifications, our findings suggest that next-generation sequencing diagnostic methods may be utilized when testing these genes in a diagnostic laboratory.

Identification of a novel insertion mutation in FGFR3 that causes thanatophoric dysplasia type 1.

Thanatophoric dysplasia is a type of short-limbed neonatal dwarfism that is usually lethal in the perinatal period. It is characterized by short limbs, a narrow, bell-shaped thorax, macrocephaly with a prominent forehead, and flattened vertebral bodies. These malformations result from autosomal dominant mutations in the fibroblast growth factor receptor 3 (FGFR3) gene. In this report, we describe a novel FGFR3 insertion mutation in a fetus with shortened limbs, curved femurs, and a narrow thorax. The diagnosis of thanatophoric dysplasia type 1 was suspected clinically, and FGFR3 sequencing showed a c.742_743insTGT variant, which predicts p.R248delinsLC. In vivo studies in zebrafish demonstrated that this mutation resulted in the overexpression of zebrafish Fgfr3, leading to the over-activation of downstream signaling and dorsalized embryos. To date, no insertions or deletions in FGFR3 have been reported to cause thanatophoric dysplasia types 1 or 2; therefore, this represents the first report to describe such a mutation. © 2016 Wiley Periodicals, Inc.

Promotion of bone morphogenetic protein signaling by tetraspanins and glycosphingolipids.

Bone morphogenetic proteins (BMPs) belong to the transforming growth factor ß (TGFß) superfamily of secreted molecules. BMPs play essential roles in multiple developmental and homeostatic processes in metazoans. Malfunction of the BMP pathway can cause a variety of diseases in humans, including cancer, skeletal disorders and cardiovascular diseases. Identification of factors that ensure proper spatiotemporal control of BMP signaling is critical for understanding how this pathway is regulated. We have used a unique and sensitive genetic screen to identify the plasma membrane-localized tetraspanin TSP-21 as a key new factor in the C. elegans BMP-like "Sma/Mab" signaling pathway that controls body size and postembryonic M lineage development. We showed that TSP-21 acts in the signal-receiving cells and genetically functions at the ligand-receptor level. We further showed that TSP-21 can associate with itself and with two additional tetraspanins, TSP-12 and TSP-14, which also promote Sma/Mab signaling. TSP-12 and TSP-14 can also associate with SMA-6, the type I receptor of the Sma/Mab pathway. Finally, we found that glycosphingolipids, major components of the tetraspanin-enriched microdomains, are required for Sma/Mab signaling. Our findings suggest that the tetraspanin-enriched membrane microdomains are important for proper BMP signaling. As tetraspanins have emerged as diagnostic and prognostic markers for tumor progression, and TSP-21, TSP-12 and TSP-14 are all conserved in humans, we speculate that abnormal BMP signaling due to altered expression or function of certain tetraspanins may be a contributing factor to cancer development.

TRPV channel-mediated calcium transients in nociceptor neurons are dispensable for avoidance behaviour.

Animals need to sense and react to potentially dangerous environments. TRP ion channels participate in nociception, presumably via Ca(2+) influx, in most animal species. However, the relationship between ion permeation and animals' nocifensive behaviour is unknown. Here we use an invertebrate animal model with relevance for mammalian pain. We analyse the putative selectivity filter of OSM-9, a TRPV channel, in osmotic avoidance behaviour of Caenorhabditis elegans. Using mutagenized OSM-9 expressed in the head nociceptor neuron, ASH, we study nocifensive behaviour and Ca(2+) influx. Within the selectivity filter, M(601)-F(609), Y604G strongly reduces avoidance behaviour and eliminates Ca(2+) transients. Y604F also abolishes Ca(2+) transients in ASH, while sustaining avoidance behaviour, yet it disrupts behavioral plasticity. Homology modelling of the OSM-9 pore suggests that Y(604) may assume a scaffolding role. Thus, aromatic residues in the OSM-9 selectivity filter are critical for pain behaviour and ion permeation. These findings have relevance for understanding evolutionary roots of mammalian nociception.

Truncating mutations in LRP4 lead to a prenatal lethal form of Cenani-Lenz syndrome.

Cenani-Lenz syndrome (CLS) is an autosomal recessive skeletal dysplasia that results in malformations of the distal limb, renal anomalies, and characteristic facies. In 2010, this condition was found to be caused by mutations in LRP4, a member of the low-density lipoprotein family of receptors. LRP4 has been shown to antagonize LRP5/LRP6 activation of WNT and ß-catenin signaling. Loss of LRP4 function leads to excessive Wnt and ß-catenin signaling in the limb bud, which causes abnormal limb development. The large majority of patients with CLS reported in the literature have splicing and missense mutations, which result in syndactyly, oligodactyly, and minor renal malformations. More recently, a patient with CLS has been identified with a homozygous nonsense mutation and a more severe presentation of findings typically associated with this condition. Here we present two sibling fetuses with a prenatal lethal presentation of mesomelic limb reductions, oligosyndactyly, genitourinary malformation and compound heterozygosity for two novel truncating mutations in LRP4. These findings lend further support to the CLS genotype-phenotype correlation presented in recent publications.

An antagonistic role for the C. elegans Schnurri homolog SMA-9 in modulating TGFbeta signaling during mesodermal patterning.

In C. elegans, the Sma/Mab TGFbeta signaling pathway regulates body size and male tail patterning. SMA-9, the C. elegans homolog of Schnurri, has been shown to function as a downstream component to mediate the Sma/Mab TGFbeta signaling pathway in these processes. We have discovered a new role for SMA-9 in dorsoventral patterning of the C. elegans post-embryonic mesoderm, the M lineage. In addition to a small body size, sma-9 mutant animals exhibit a dorsal-to-ventral fate transformation within the M lineage. This M lineage defect of sma-9 mutants is unique in that animals carrying mutations in all other known components of the TGFbeta pathway exhibit no M lineage defects. Surprisingly, mutations in the core components of the Sma/Mab TGFbeta signaling pathway suppressed the M lineage defects of sma-9 mutants without suppressing their body size defects. We show that this suppression specifically happens within the M lineage. Our studies have uncovered an unexpected role of SMA-9 in antagonizing the TGFbeta signaling pathway during mesodermal patterning, suggesting a novel mode of function for the SMA-9/Schnurri family of proteins.