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.

Connective tissue presentation in two families expands the phenotypic spectrum of PYROXD1 disorders.

Recessive variants in the oxidoreductase PYROXD1 are reported to cause a myopathy in 22 affected individuals from 15 families. Here, we describe two female probands from unrelated families presenting with features of a congenital connective tissue disorder including osteopenia, blue sclera, soft skin, joint hypermobility and neuromuscular junction dysfunction in addition to known features of PYROXD1 myopathy including respiratory difficulties, weakness, hypotonia and oromotor dysfunction. Proband AII:1 is compound heterozygous for the recurrent PYROXD1 variant Chr12(GRCh38):g.21452130A>G;NM_024854.5:c.464A>G;p.(N155S) and Chr12(GRCh38):g.21462019_21462022del;NM_024854.5:c.892_895del;p.(V298Mfs*4) and proband BII:1 is compound heterozygous for Chr12(GRCh38):g.21468739-21468741del;NM_024854.5:c.1488_1490del;p.(E496del) and Chr12(GRCh38):g.21467619del;NM_024854.5:c.1254+1del. RNA studies demonstrate c.892_895del;p.(V298Mfs*4) is targeted by nonsense mediated decay and c.1254+1delG elicits in-frame skipping of exon-11. Western blot from cultured fibroblasts shows reduced PYROXD1 protein levels in both probands. Testing urine from BII:1 and six individuals with PYROXD1 myopathy showed elevated levels of deoxypyridinoline, a mature collagen crosslink, correlating with PYROXD1-disorder severity. Urine and serum amino acid testing of the same individuals revealed no reportable changes. In contrast to PYROXD1 knock-out, we find no evidence for disrupted tRNA ligase activity, as measured via XBP1 splicing, in fibroblasts expressing PYROXD1 variants. In summary, we expand the clinical spectrum of PYROXD1-related disorders to include an overlapping connective tissue and myopathy presentation, identify three novel, pathogenic PYROXD1 variants, and provide preliminary evidence that elevated urine DPD crosslinks may provide a clinical biomarker for PYROXD1 disorders. Our results advocate consideration of PYROXD1 variants in the differential diagnosis for undiagnosed individuals presenting with a connective tissue disorder and myopathy.

Case report: Adult-onset limb girdle muscular dystrophy in sibling pair due to novel homozygous LAMA2 missense variant.

Recessive pathogenic variants in the laminin subunit alpha 2 (LAMA2) gene cause a spectrum of disease ranging from severe congenital muscular dystrophy to later-onset limb girdle muscular dystrophy (LGMDR23). The phenotype of LGMDR23 is characterized by slowly progressive proximal limb weakness, contractures, raised creatine kinase, and sometimes distinctive cerebral white matter changes and/or epilepsy. We present two siblings, born to consanguineous parents, who developed adult-onset LGMDR23 associated with typical cerebral white matter changes and who both later developed dementia. The male proband also had epilepsy and upper motor neuron signs when he presented at age 72. Merosin immunohistochemistry and Western blot on muscle biopsies taken from both subjects was normal. Whole exome sequencing revealed a previously unreported homozygous missense variant in LAMA2 [Chr6(GRCh38):g.129297734G>A; NM_000426.3:c.2906G>A; p.(Cys969Tyr)] in the proband. The same homozygous LAMA2 variant was confirmed by Sanger sequencing in the proband's affected sister. These findings expand the genotypic and phenotypic spectrum of LGMDR23.

Compound heterozygous splicing variants expand the genotypic spectrum of EMC1-related disorders.

EMC1 encodes subunit 1 of the endoplasmic reticulum (ER) membrane protein complex (EMC), a transmembrane domain insertase involved in membrane protein biosynthesis. Variants in EMC1 are described as a cause of global developmental delay, hypotonia, cortical visual impairment, and commonly, cerebral atrophy on MRI scan. We report an individual with severe global developmental delay and progressive cerebellar atrophy in whom exome sequencing identified a heterozygous essential splice-site variant in intron-3 of EMC1 (NM_015047.3:c.287-1G>A). Whole genome sequencing (WGS) identified a deep intronic variant in intron-20 of EMC1 (NM_015047.3:c.2588-771C>G) that was poorly predicted by in silico programs to disrupt pre-mRNA splicing. Reverse Transcription-PCR (RT-PCR) revealed stochastic activation of a pseudo-exon associated with the c.2588-771C>G variant and mis-splicing arising from the c.287-1G>A variant. This case highlights the utility of WGS and RNA studies to identify and assess likely pathogenicity of deep intronic variants and expands the genotypic and phenotypic spectrum of EMC1-related disorders.

SpliceVault predicts the precise nature of variant-associated mis-splicing.

Even for essential splice-site variants that are almost guaranteed to alter mRNA splicing, no current method can reliably predict whether exon-skipping, cryptic activation or multiple events will result, greatly complicating clinical interpretation of pathogenicity. Strikingly, ranking the four most common unannotated splicing events across 335,663 reference RNA-sequencing (RNA-seq) samples (300K-RNA Top-4) predicts the nature of variant-associated mis-splicing with 92% sensitivity. The 300K-RNA Top-4 events correctly identify 96% of exon-skipping events and 86% of cryptic splice sites for 140 clinical cases subject to RNA testing, showing higher sensitivity and positive predictive value than SpliceAI. Notably, RNA re-analyses showed we had missed 300K-RNA Top-4 events for several clinical cases tested before the development of this empirical predictive method. Simply, mis-splicing events that happen around a splice site in RNA-seq data are those most likely to be activated by a splice-site variant. The SpliceVault web portal allows users easy access to 300K-RNA for informed splice-site variant interpretation and classification.

A genetic basis is identified in 74% cases of paediatric hyperCKaemia without weakness presenting to a tertiary paediatric neuromuscular centre.

Paediatric hyperCKaemia without weakness presents a clinical conundrum. Invasive investigations with low diagnostic yields, including muscle biopsy, may be considered unjustifiable. Improved access to genome-wide genetic testing has shifted first-line investigations towards genetic studies in neuromuscular disease. This research aims to provide an evidence-based diagnostic approach to paediatric hyperCKaemia without weakness, a current gap in the literature. We identified 47 individuals (10-months to 16-years-old; 34 males, 13 females) from 43 families presenting with hyperCKaemia on two or more occasions, without weakness, from The Children's Hospital at Westmead Neuromuscular Clinic Database. Clinical features, investigations and outcomes were analysed via retrospective chart review. Genetic testing has been performed in 34/43. Genetic variants explaining hyperCKaemia were identified in 25/34 (74%) using multiplex ligation-dependent probe amplification, massive parallel sequencing, single gene testing and exome sequencing. Pathogenic/likely pathogenic variants were identified in 19 neuromuscular disease genes and six metabolic myopathy genes. Individuals with metabolic diagnoses had higher peak creatine kinase levels that sometimes normalized. Conversely, creatine kinase levels remained persistently elevated those with neuromuscular diagnoses. In summary, a genetic cause is found in most paediatric patients with hyperCKaemia without weakness informing clinical management and counselling. Thus, we propose a diagnostic algorithm for this cohort.

Prevalence, parameters, and pathogenic mechanisms for splice-altering acceptor variants that disrupt the AG exclusion zone.

Predicting the pathogenicity of acceptor splice-site variants outside the essential AG is challenging, due to high sequence diversity of the extended splice-site region. Critical analysis of 24,445 intronic extended acceptor splice-site variants reported in ClinVar and the Leiden Open Variation Database (LOVD) demonstrates 41.9% of pathogenic variants create an AG dinucleotide between the predicted branchpoint and acceptor (AG-creating variants in the AG exclusion zone), 28.4% result in loss of a pyrimidine at the -3 position, and 15.1% result in loss of one or more pyrimidines in the polypyrimidine tract. Pathogenicity of AG-creating variants was highly influenced by their position. We define a high-risk zone for pathogenicity: > 6 nucleotides downstream of the predicted branchpoint and >5 nucleotides upstream from the acceptor, where 93.1% of pathogenic AG-creating variants arise and where naturally occurring AG dinucleotides are concordantly depleted (5.8% of natural AGs). SpliceAI effectively predicts pathogenicity of AG-creating variants, achieving 95% sensitivity and 69% specificity. We highlight clinical examples showing contrasting mechanisms for mis-splicing arising from AG variants: (1) cryptic acceptor created; (2) splicing silencer created: an introduced AG silences the acceptor, resulting in exon skipping, intron retention, and/or use of an alternative existing cryptic acceptor; and (3) splicing silencer disrupted: loss of a deep intronic AG activates inclusion of a pseudo-exon. In conclusion, we establish AG-creating variants as a common class of pathogenic extended acceptor variant and outline factors conferring critical risk for mis-splicing for AG-creating variants in the AG exclusion zone, between the branchpoint and acceptor.

Standardized practices for RNA diagnostics using clinically accessible specimens reclassifies 75% of putative splicing variants.

PURPOSE: Genetic variants causing aberrant premessenger RNA splicing are increasingly being recognized as causal variants in genetic disorders. In this study, we devise standardized practices for polymerase chain reaction (PCR)-based RNA diagnostics using clinically accessible specimens (blood, fibroblasts, urothelia, biopsy). METHODS: A total of 74 families with diverse monogenic conditions (31% prenatal-congenital onset, 47% early childhood, and 22% teenage-adult onset) were triaged into PCR-based RNA testing, with comparative RNA sequencing for 19 cases. RESULTS: Informative RNA assay data were obtained for 96% of cases, enabling variant reclassification for 75% variants that can be used for genetic counseling (71%), to inform clinical care (32%) and prenatal counseling (41%). Variant-associated mis-splicing was highly reproducible for 28 cases with samples from ≥2 affected individuals or heterozygotes and 10 cases with ≥2 biospecimens. PCR amplicons encompassing another segregated heterozygous variant was vital for clinical interpretation of 22 of 79 variants to phase RNA splicing events and discern complete from partial mis-splicing. CONCLUSION: RNA diagnostics enabled provision of a genetic diagnosis for 64% of recruited cases. PCR-based RNA diagnostics has capacity to analyze 81.3% of clinically significant genes, with long amplicons providing an advantage over RNA sequencing to phase RNA splicing events. The Australasian Consortium for RNA Diagnostics (SpliceACORD) provide clinically-endorsed, standardized protocols and recommendations for interpreting RNA assay data.

WGS and RNA Studies Diagnose Noncoding DMD Variants in Males With High Creatine Kinase.

OBJECTIVE: To describe the diagnostic utility of whole-genome sequencing and RNA studies in boys with suspected dystrophinopathy, for whom multiplex ligation-dependent probe amplification and exomic parallel sequencing failed to yield a genetic diagnosis, and to use remnant normal DMD splicing in 3 families to define critical levels of wild-type dystrophin bridging clinical spectrums of Duchenne to myalgia. METHODS: Exome, genome, and/or muscle RNA sequencing was performed for 7 males with elevated creatine kinase. PCR of muscle-derived complementary DNA (cDNA) studied consequences for DMD premessenger RNA (pre-mRNA) splicing. Quantitative Western blot was used to determine levels of dystrophin, relative to control muscle. RESULTS: Splice-altering intronic single nucleotide variants or structural rearrangements in DMD were identified in all 7 families. Four individuals, with abnormal splicing causing a premature stop codon and nonsense-mediated decay, expressed remnant levels of normally spliced DMD mRNA. Quantitative Western blot enabled correlation of wild-type dystrophin and clinical severity, with 0%-5% dystrophin conferring a Duchenne phenotype, 10% ± 2% a Becker phenotype, and 15% ± 2% dystrophin associated with myalgia without manifesting weakness. CONCLUSIONS: Whole-genome sequencing relied heavily on RNA studies to identify DMD splice-altering variants. Short-read RNA sequencing was regularly confounded by the effectiveness of nonsense-mediated mRNA decay and low read depth of the giant DMD mRNA. PCR of muscle cDNA provided a simple, yet informative approach. Highly relevant to genetic therapies for dystrophinopathies, our data align strongly with previous studies of mutant dystrophin in Becker muscular dystrophy, with the collective conclusion that a fractional increase in levels of normal dystrophin between 5% and 20% is clinically significant.

Pathogenic deep intronic MTM1 variant activates a pseudo-exon encoding a nonsense codon resulting in severe X-linked myotubular myopathy.

X-linked myotubular myopathy (XLMTM) is a severe congenital myopathy characterised by generalised weakness and respiratory insufficiency. XLMTM is associated with pathogenic variants in MTM1; a gene encoding the lipid phosphatase myotubularin. Whole genome sequencing (WGS) of an exome-negative male proband with severe hypotonia, respiratory insufficiency and centralised nuclei on muscle biopsy identified a deep intronic MTM1 variant NG_008199.1(NM_000252.2):c.1468-577A>G, which strengthened a cryptic 5' splice site (A>G substitution at the +5 position). Muscle RNA sequencing was non-diagnostic due to low read depth. Reverse transcription PCR (RT-PCR) of muscle RNA confirmed the c.1468-577A>G variant activates inclusion of a pseudo-exon encoding a premature stop codon into all detected MTM1 transcripts. Western blot analysis establishes deficiency of myotubularin protein, consistent with the severe XLMTM phenotype. We expand the genotypic spectrum of XLMTM and highlight benefits of screening non-coding regions of MTM1 in male probands with phenotypically concordant XLMTM who remain undiagnosed following exome sequencing.

Recurrent TTN metatranscript-only c.39974-11T>G splice variant associated with autosomal recessive arthrogryposis multiplex congenita and myopathy.

We present eight families with arthrogryposis multiplex congenita and myopathy bearing a TTN intron 213 extended splice-site variant (NM_001267550.1:c.39974-11T>G), inherited in trans with a second pathogenic TTN variant. Muscle-derived RNA studies of three individuals confirmed mis-splicing induced by the c.39974-11T>G variant; in-frame exon 214 skipping or use of a cryptic 3' splice-site effecting a frameshift. Confounding interpretation of pathogenicity is the absence of exons 213-217 within the described skeletal muscle TTN N2A isoform. However, RNA-sequencing from 365 adult human gastrocnemius samples revealed that 56% specimens predominantly include exons 213-217 in TTN transcripts (inclusion rate ≥66%). Further, RNA-sequencing of five fetal muscle samples confirmed that 4/5 specimens predominantly include exons 213-217 (fifth sample inclusion rate 57%). Contractures improved significantly with age for four individuals, which may be linked to decreased expression of pathogenic fetal transcripts. Our study extends emerging evidence supporting a vital developmental role for TTN isoforms containing metatranscript-only exons.

Importance of muscle biopsy to establish pathogenicity of DMD missense and splice variants.

A precise genetic diagnosis of a dystrophinopathy has far-reaching implications for affected boys and their families. We present three boys with DMD single nucleotide variants associated with Becker muscular dystrophy presenting with myalgia, reduced exercise capacity, neurodevelopmental symptoms and elevated creatine kinase. The DMD variants were difficult to classify: AIII:1 a synonymous variant in exon 13 c.1602G>A, p.Lys534Lys; BIII:1 an essential splice-site variant in intron 33 c.4674+1G>A, and CII:1 a missense mutation within the cysteine-rich domain, exon 66 c.9619T>C, p.Cys3207Arg. Complementary DNA (cDNA) analysis using muscle-derived mRNA established splice-altering effects of variants for AIII:1 and BIII:1, and normal splicing in CII:1. Western blot analysis demonstrated mildly to moderately reduced dystrophin levels (17.6 - 36.1% the levels of controls), supporting dystrophinopathy as a probable diagnosis. These three cases highlight the diagnostic utility of muscle biopsy for mRNA studies and western blot to investigate DMD variants of uncertain pathogenicity, by exploring effects on splicing and dystrophin protein levels.

Pathogenic Abnormal Splicing Due to Intronic Deletions that Induce Biophysical Space Constraint for Spliceosome Assembly.

A precise genetic diagnosis is the single most important step for families with genetic disorders to enable personalized and preventative medicine. In addition to genetic variants in coding regions (exons) that can change a protein sequence, abnormal pre-mRNA splicing can be devastating for the encoded protein, inducing a frameshift or in-frame deletion/insertion of multiple residues. Non-coding variants that disrupt splicing are extremely challenging to identify. Stemming from an initial clinical discovery in two index Australian families, we define 25 families with genetic disorders caused by a class of pathogenic non-coding splice variant due to intronic deletions. These pathogenic intronic deletions spare all consensus splice motifs, though they critically shorten the minimal distance between the 5' splice-site (5'SS) and branchpoint. The mechanistic basis for abnormal splicing is due to biophysical constraint precluding U1/U2 spliceosome assembly, which stalls in A-complexes (that bridge the 5'SS and branchpoint). Substitution of deleted nucleotides with non-specific sequences restores spliceosome assembly and normal splicing, arguing against loss of an intronic element as the primary causal basis. Incremental lengthening of 5'SS-branchpoint length in our index EMD case subject defines 45-47 nt as the critical elongation enabling (inefficient) spliceosome assembly for EMD intron 5. The 5'SS-branchpoint space constraint mechanism, not currently factored by genomic informatics pipelines, is relevant to diagnosis and precision medicine across the breadth of Mendelian disorders and cancer genomics.

Recessive DES cardio/myopathy without myofibrillar aggregates: intronic splice variant silences one allele leaving only missense L190P-desmin.

We establish autosomal recessive DES variants p.(Leu190Pro) and a deep intronic splice variant causing inclusion of a frameshift-inducing artificial exon/intronic fragment, as the likely cause of myopathy with cardiac involvement in female siblings. Both sisters presented in their twenties with slowly progressive limb girdle weakness, severe systolic dysfunction, and progressive, severe respiratory weakness. Desmin is an intermediate filament protein typically associated with autosomal dominant myofibrillar myopathy with cardiac involvement. However a few rare cases of autosomal recessive desminopathy are reported. In this family, a paternal missense p.(Leu190Pro) variant was viewed unlikely to be causative of autosomal dominant desminopathy, as the father and brothers carrying this variant were clinically unaffected. Clinical fit with a DES-related myopathy encouraged closer scrutiny of all DES variants, identifying a maternal deep intronic variant within intron-7, predicted to create a cryptic splice site, which segregated with disease. RNA sequencing and studies of muscle cDNA confirmed the deep intronic variant caused aberrant splicing of an artificial exon/intronic fragment into maternal DES mRNA transcripts, encoding a premature termination codon, and potently activating nonsense-mediate decay (92% paternal DES transcripts, 8% maternal). Western blot showed 60-75% reduction in desmin levels, likely comprised only of missense p.(Leu190Pro) desmin. Biopsy showed fibre size variation with increased central nuclei. Electron microscopy showed extensive myofibrillar disarray, duplication of the basal lamina, but no inclusions or aggregates. This study expands the phenotypic spectrum of recessive DES cardio/myopathy, and emphasizes the continuing importance of muscle biopsy for functional genomics pursuit of 'tricky' variants in neuromuscular conditions.