RNA variant assessment using transactivation and transdifferentiation.

Understanding the impact of splicing and nonsense variants on RNA is crucial for the resolution of variant classification as well as their suitability for precision medicine interventions. This is primarily enabled through RNA studies involving transcriptomics followed by targeted assays using RNA isolated from clinically accessible tissues (CATs) such as blood or skin of affected individuals. Insufficient disease gene expression in CATs does however pose a major barrier to RNA based investigations, which we show is relevant to 1,436 Mendelian disease genes. We term these "silent" Mendelian genes (SMGs), the largest portion (36%) of which are associated with neurological disorders. We developed two approaches to induce SMG expression in human dermal fibroblasts (HDFs) to overcome this limitation, including CRISPR-activation-based gene transactivation and fibroblast-to-neuron transdifferentiation. Initial transactivation screens involving 40 SMGs stimulated our development of a highly multiplexed transactivation system culminating in the 6- to 90,000-fold induction of expression of 20/20 (100%) SMGs tested in HDFs. Transdifferentiation of HDFs directly to neurons led to expression of 193/516 (37.4%) of SMGs implicated in neurological disease. The magnitude and isoform diversity of SMG expression following either transactivation or transdifferentiation was comparable to clinically relevant tissues. We apply transdifferentiation and/or gene transactivation combined with short- and long-read RNA sequencing to investigate the impact that variants in USH2A, SCN1A, DMD, and PAK3 have on RNA using HDFs derived from affected individuals. Transactivation and transdifferentiation represent rapid, scalable functional genomic solutions to investigate variants impacting SMGs in the patient cell and genomic context.

Lmod2 is necessary for effective skeletal muscle contraction.

Muscle contraction is a regulated process driven by the sliding of actin-thin filaments over myosin-thick filaments. Lmod2 is an actin filament length regulator and essential for life since human mutations and complete loss of Lmod2 in mice lead to dilated cardiomyopathy and death. To study the little-known role of Lmod2 in skeletal muscle, we created a mouse model with Lmod2 expressed exclusively in the heart but absent in skeletal muscle. Loss of Lmod2 in skeletal muscle results in decreased force production in fast- and slow-twitch muscles. Soleus muscle from rescued Lmod2 knockout mice have shorter thin filaments, increased Lmod3 levels, and present with a myosin fiber type switch from fast myosin heavy chain (MHC) IIA to the slower MHC I isoform. Since Lmod2 regulates thin-filament length in slow-twitch but not fast-twitch skeletal muscle and force deficits were observed in both muscle types, this work demonstrates that Lmod2 regulates skeletal muscle contraction, independent of its role in thin-filament length regulation.

Refining clinically relevant parameters for mis-splicing risk in shortened introns with donor-to-branchpoint space constraint.

Intronic deletions that critically shorten donor-to-branchpoint (D-BP) distance of a precursor mRNA impose biophysical space constraint on assembly of the U1/U2 spliceosomal complex, leading to canonical splicing failure. Here we use a series of ß-globin (HBB) gene constructs with intron 1 deletions to define D-BP lengths that present low/no risk of mis-splicing and lengths which are critically short and likely elicit clinically relevant mis-splicing. We extend our previous observation in EMD intron 5 of 46 nt as the minimal productive D-BP length, demonstrating spliceosome assembly constraint persists at D-BP lengths of 47-56 nt. We exploit the common HBB exon 1 ß-thalassemia variant that strengthens a cryptic donor (NM_000518.5(HBB):c.79G > A) to provide a simple barometer for the earliest signs of space constraint, via cryptic donor activation. For clinical evaluation of intronic deletions, we assert D-BP lengths > 60 nt present low mis-splicing risk while space constraint increases exponentially with D-BP lengths < 55 nt, with critical risk and profound splicing abnormalities with D-BP lengths < 50 nt.

Genome sequencing enables diagnosis and treatment of SLC5A6 neuropathy.

The sodium-dependent multivitamin transporter encoded by SLC5A6 is responsible for uptake of biotin, pantothenic acid, and α-lipoic acid. Thirteen individuals from eight families are reported with pathogenic biallelic SLC5A6 variants. Phenotype ranges from multisystem metabolic disorder to childhood-onset peripheral motor neuropathy. We report three additional affected individuals with biallelic SLC5A6 variants. In Family A, a male proband (AII:1) presenting in early childhood with gross motor regression, motor axonal neuropathy, recurrent cytopenia and infections, and failure to thrive was diagnosed at 12 years of age via genome sequencing (GS) with a paternal NM_021095.4:c.393+2T>C variant and a maternal c.1285A>G p.(Ser429Gly) variant. An uncle with recurrent cytopenia and peripheral neuropathy was subsequently found to have the same genotype. We also report an unrelated female with peripheral neuropathy homozygous for the c.1285A>G p.(Ser429Gly) recurrent variant identified in seven reported cases, including this study. RT-PCR studies on blood mRNA from AII:1 showed c.393+2T>C caused mis-splicing with all canonically spliced transcripts in AII:1 containing the c.1285A>G variant. SLC5A6 mRNA expression in AII:1 fibroblasts was ~50% of control levels, indicative of nonsense-mediated decay of mis-spliced transcripts. Biotin uptake studies on AII:1 fibroblasts, expressing the p.(Ser429Gly) variant, showed an ~90% reduction in uptake compared to controls. Targeted treatment of AII:1 with biotin, pantothenic acid, and lipoic acid resulted in clinical improvement. Health Economic analyses showed implementation of GS as an early investigation could have saved $ AUD 105,988 and shortened diagnostic odyssey and initiation of treatment by up to 7 years.

Expert Panel Curation of 31 Genes in Relation to Limb Girdle Muscular Dystrophy.

Introduction: Limb girdle muscular dystrophies (LGMDs) are a group of genetically heterogeneous autosomal conditions with some degree of phenotypic homogeneity. LGMD is defined as having onset >2 years of age with progressive proximal weakness, elevated serum creatine kinase levels and dystrophic features on muscle biopsy. Advances in massively parallel sequencing have led to a surge in genes linked to LGMD. Methods: The ClinGen Muscular Dystrophies and Myopathies gene curation expert panel (MDM GCEP, formerly Limb Girdle Muscular Dystrophy GCEP) convened to evaluate the strength of evidence supporting gene-disease relationships (GDR) using the ClinGen gene-disease clinical validity framework to evaluate 31 genes implicated in LGMD. Results: The GDR was exclusively LGMD for 17 genes, whereas an additional 14 genes were related to a broader phenotype encompassing congenital weakness. Four genes (CAPN3, COL6A1, COL6A2, COL6A3) were split into two separate disease entities, based on each displaying both dominant and recessive inheritance patterns, resulting in curation of 35 GDRs. Of these, 30 (86%) were classified as Definitive, 4 (11%) as Moderate and 1 (3%) as Limited. Two genes, POMGNT1 and DAG1, though definitively related to myopathy, currently have insufficient evidence to support a relationship specifically with LGMD. Conclusions: The expert-reviewed assertions on the clinical validity of genes implicated in LGMDs form an invaluable resource for clinicians and molecular geneticists. We encourage the global neuromuscular community to publish case-level data that help clarify disputed or novel LGMD associations.

Genome and RNA sequencing boost neuromuscular diagnoses to 62% from 34% with exome sequencing alone.

OBJECTIVE: Most families with heritable neuromuscular disorders do not receive a molecular diagnosis. Here we evaluate diagnostic utility of exome, genome, RNA sequencing, and protein studies and provide evidence-based recommendations for their integration into practice. METHODS: In total, 247 families with suspected monogenic neuromuscular disorders who remained without a genetic diagnosis after standard diagnostic investigations underwent research-led massively parallel sequencing: neuromuscular disorder gene panel, exome, genome, and/or RNA sequencing to identify causal variants. Protein and RNA studies were also deployed when required. RESULTS: Integration of exome sequencing and auxiliary genome, RNA and/or protein studies identified causal or likely causal variants in 62% (152 out of 247) of families. Exome sequencing alone informed 55% (83 out of 152) of diagnoses, with remaining diagnoses (45%; 69 out of 152) requiring genome sequencing, RNA and/or protein studies to identify variants and/or support pathogenicity. Arrestingly, novel disease genes accounted for <4% (6 out of 152) of diagnoses while 36.2% of solved families (55 out of 152) harbored at least one splice-altering or structural variant in a known neuromuscular disorder gene. We posit that contemporary neuromuscular disorder gene-panel sequencing could likely provide 66% (100 out of 152) of our diagnoses today. INTERPRETATION: Our results emphasize thorough clinical phenotyping to enable deep scrutiny of all rare genetic variation in phenotypically consistent genes. Post-exome auxiliary investigations extended our diagnostic yield by 81% overall (34-62%). We present a diagnostic algorithm that details deployment of genomic and auxiliary investigations to obtain these diagnoses today most effectively. We hope this provides a practical guide for clinicians as they gain greater access to clinical genome and transcriptome sequencing.

Digenic inheritance involving a muscle-specific protein kinase and the giant titin protein causes a skeletal muscle myopathy.

In digenic inheritance, pathogenic variants in two genes must be inherited together to cause disease. Only very few examples of digenic inheritance have been described in the neuromuscular disease field. Here we show that predicted deleterious variants in SRPK3, encoding the X-linked serine/argenine protein kinase 3, lead to a progressive early onset skeletal muscle myopathy only when in combination with heterozygous variants in the TTN gene. The co-occurrence of predicted deleterious SRPK3/TTN variants was not seen among 76,702 healthy male individuals, and statistical modeling strongly supported digenic inheritance as the best-fitting model. Furthermore, double-mutant zebrafish (srpk3-/-; ttn.1+/-) replicated the myopathic phenotype and showed myofibrillar disorganization. Transcriptome data suggest that the interaction of srpk3 and ttn.1 in zebrafish occurs at a post-transcriptional level. We propose that digenic inheritance of deleterious changes impacting both the protein kinase SRPK3 and the giant muscle protein titin causes a skeletal myopathy and might serve as a model for other genetic diseases.