Exome and whole-genome sequencing are becoming increasingly routine approaches in Mendelian disease diagnosis. Despite their success, the current diagnostic rate for genomic analyses across a variety of rare diseases is approximately 25 to 50%. We explore the utility of transcriptome sequencing [RNA sequencing (RNA-seq)] as a complementary diagnostic tool in a cohort of 50 patients with genetically undiagnosed rare muscle disorders. We describe an integrated approach to analyze patient muscle RNA-seq, leveraging an analysis framework focused on the detection of transcript-level changes that are unique to the patient compared to more than 180 control skeletal muscle samples. We demonstrate the power of RNA-seq to validate candidate splice-disrupting mutations and to identify splice-altering variants in both exonic and deep intronic regions, yielding an overall diagnosis rate of 35%. We also report the discovery of a highly recurrent de novo intronic mutation in COL6A1 that results in a dominantly acting splice-gain event, disrupting the critical glycine repeat motif of the triple helical domain. We identify this pathogenic variant in a total of 27 genetically unsolved patients in an external collagen VI–like dystrophy cohort, thus explaining approximately 25% of patients clinically suggestive of having collagen VI dystrophy in whom prior genetic analysis is negative. Overall, this study represents a large systematic application of transcriptome sequencing to rare disease diagnosis and highlights its utility for the detection and interpretation of variants missed by current standard diagnostic approaches.

The T-box family transcription factor gene TBX20 acts in a conserved regulatory network, guiding heart formation and patterning in diverse species. Mouse Tbx20 is expressed in cardiac progenitor cells, differentiating cardiomyocytes, and developing valvular tissue, and its deletion or RNA interference-mediated knockdown is catastrophic for heart development. TBX20 interacts physically, functionally, and genetically with other cardiac transcription factors, including NKX2-5, GATA4, and TBX5, mutations of which cause congenital heart disease (CHD). Here, we report nonsense (Q195X) and missense (I152M) germline mutations within the T-box DNA-binding domain of human TBX20 that were associated with a family history of CHD and a complex spectrum of developmental anomalies, including defects in septation, chamber growth, and valvulogenesis. Biophysical characterization of wild-type and mutant proteins indicated how the missense mutation disrupts the structure and function of the TBX20 T-box. Dilated cardiomyopathy was a feature of the TBX20 mutant phenotype in humans and mice, suggesting that mutations in developmental transcription factors can provide a sensitized template for adult-onset heart disease. Our findings are the first to link TBX20 mutations to human pathology. They provide insights into how mutation of different genes in an interactive regulatory circuit lead to diverse clinical phenotypes, with implications for diagnosis, genetic screening, and patient follow-up.

Nemaline myopathy (NM) is a common form of congenital myopathy, affecting approximately 1 in 50,000 individuals, and is defined by the presence of nonprogressive generalized muscle weakness and numerous electron-dense protein inclusions (nemaline bodies or rods) in skeletal myofibers (1). The most severely affected Nemaline myopathy (NM) is a genetic muscle disorder characterized by muscle dysfunction and electron-dense protein accumulations (nemaline bodies) in myofibers. Pathogenic mutations have been described in 9 genes to date, but the genetic basis remains unknown in many cases. Here, using an approach that combined whole-exome sequencing (WES) and Sanger sequencing, we identified homozygous or compound heterozygous variants in LMOD3 in 21 patients from 14 families with severe, usually lethal, NM. LMOD3 encodes leiomodin-3 (LMOD3), a 65-kDa protein expressed in skeletal and cardiac muscle. LMOD3 was expressed from early stages of muscle differentiation; localized to actin thin filaments, with enrichment near the pointed ends; and had strong actin filament-nucleating activity. Loss of LMOD3 in patient muscle resulted in shortening and disorganization of thin filaments. Knockdown of lmod3 in zebrafish replicated NM-associated functional and pathological phenotypes. Together, these findings indicate that mutations in the gene encoding LMOD3 underlie congenital myopathy and demonstrate that LMOD3 is essential for the organization of sarcomeric thin filaments in skeletal muscle.

Congenital muscular dystrophies with hypoglycosylation of α-dystroglycan (α-DG) are a heterogeneous group of disorders often associated with brain and eye defects in addition to muscular dystrophy. Causative variants in 14 genes thought to be involved in the glycosylation of α-DG have been identified thus far. Allelic mutations in these genes might also cause milder limb-girdle muscular dystrophy phenotypes. Using a combination of exome and Sanger sequencing in eight unrelated individuals, we present evidence that mutations in guanosine diphosphate mannose (GDP-mannose) pyrophosphorylase B (GMPPB) can result in muscular dystrophy variants with hypoglycosylated α-DG. GMPPB catalyzes the formation of GDP-mannose from GTP and mannose-1-phosphate. GDP-mannose is required for O-mannosylation of proteins, including α-DG, and it is the substrate of cytosolic mannosyltransferases. We found reduced α-DG glycosylation in the muscle biopsies of affected individuals and in available fibroblasts. Overexpression of wild-type GMPPB in fibroblasts from an affected individual partially restored glycosylation of α-DG. Whereas wild-type GMPPB localized to the cytoplasm, five of the identified missense mutations caused formation of aggregates in the cytoplasm or near membrane protrusions. Additionally, knockdown of the GMPPB ortholog in zebrafish caused structural muscle defects with decreased motility, eye abnormalities, and reduced glycosylation of α-DG. Together, these data indicate that GMPPB mutations are responsible for congenital and limb-girdle muscular dystrophies with hypoglycosylation of α-DG.

One Sentence Summary: Transcriptome sequencing improves the diagnostic rate for Mendelian disease in patients for whom genetic analysis has not returned a diagnosis. AbstractExome and whole-genome sequencing are becoming increasingly routine approaches in Mendelian disease diagnosis. Despite their success, the current diagnostic rate for genomic analyses across a variety of rare diseases is approximately 25-50%. Here, we explore the utility of transcriptome sequencing (RNA-seq) as a complementary diagnostic tool in a cohort of 50 patients with genetically undiagnosed rare muscle disorders. We describe an integrated approach to analyze patient muscle RNA-seq, leveraging an analysis framework focused on the detection of transcript-level changes that are unique to the patient compared to over 180 control skeletal muscle samples. We demonstrate the power of RNA-seq to validate candidate splice-disrupting mutations and to identify splice-altering variants in both exonic and deep intronic regions, yielding an overall diagnosis rate of 35%. We also report the discovery of a highly recurrent de novo intronic mutation in COL6A1 that results in a dominantly acting splice-gain event, disrupting the critical glycine repeat motif of the triple helical domain. We identify this pathogenic variant in a total of 27 genetically unsolved patients in an external collagen VI-like dystrophy cohort, thus explaining approximately 25% of patients clinically suggestive of collagen VI dystrophy in whom prior genetic analysis is negative. Overall, this study represents a large systematic application of transcriptome sequencing to rare disease diagnosis and highlights its utility for the detection and interpretation of variants missed by current standard diagnostic approaches. IntroductionThe advent of exome (WES) and whole genome (WGS) sequencing has greatly accelerated our capacity to identify variants that explain many Mendelian diseases in both known and new disease genes. While these technologies are mainstays in Mendelian disease diagnosis, their success rate for detecting causal variants is far from complete, ranging from 25-50% (1-4). The primary challenge of these genome-based diagnostics is that the capacity of WES and WGS to discover genetic variants substantially exceeds our ability to interpret their functional and clinical impact (5-7).

IMPORTANCETo our knowledge, the efficacy of transferring next-generation sequencing from a research setting to neuromuscular clinics has never been evaluated.OBJECTIVE To translate whole-exome sequencing (WES) to clinical practice for the genetic diagnosis of a large cohort of patients with limb-girdle muscular dystrophy (LGMD) for whom protein-based analyses and targeted Sanger sequencing failed to identify the genetic cause of their disorder. DESIGN, SETTING, AND PARTICIPANTSWe performed WES on 60 families with LGMDs (100 exomes). Data analysis was performed between January 6 and December 19, 2014, using the xBrowse bioinformatics interface (Broad Institute). Patients with LGMD were ascertained retrospectively through the Institute for Neuroscience and Muscle Research Biospecimen Bank between 2006 and 2014. Enrolled patients had been extensively investigated via protein studies and candidate gene sequencing and remained undiagnosed. Patients presented with more than 2 years of muscle weakness and with dystrophic or myopathic changes present in muscle biopsy specimens. MAIN OUTCOMES AND MEASURESThe diagnostic rate of LGMD in Australia and the relative frequencies of the different LGMD subtypes. Our central goals were to improve the genetic diagnosis of LGMD, investigate whether the WES platform provides adequate coverage of known LGMD-related genes, and identify new LGMD-related genes.RESULTS With WES, we identified likely pathogenic mutations in known myopathy genes for 27 of 60 families. Twelve families had mutations in known LGMD-related genes. However, 15 families had variants in disease-related genes not typically associated with LGMD, highlighting the clinical overlap between LGMD and other myopathies. Common causes of phenotypic overlap were due to mutations in congenital muscular dystrophy-related genes (4 families) and collagen myopathy-related genes (4 families). Less common myopathies included metabolic myopathy (2 families), congenital myasthenic syndrome (DOK7), congenital myopathy (ACTA1), tubular aggregate myopathy (STIM1), myofibrillar myopathy (FLNC), and mutation of CHD7, usually associated with the CHARGE syndrome. Inclusion of family members increased the diagnostic efficacy of WES, with a diagnostic rate of 60% for "trios" (an affected proband with both parents) vs 40% for single probands. A follow-up screening of patients whose conditions were undiagnosed on a targeted neuromuscular disease-related gene panel did not improve our diagnostic yield. CONCLUSIONS AND RELEVANCEWith WES, we achieved a diagnostic success rate of 45.0% in our difficult-to-diagnose cohort of patients with LGMD. We expand the clinical phenotypes associated with known myopathy genes, and we stress the importance of accurate clinical examination and histopathological results for interpretation of WES, with many diagnoses requiring follow-up review and ancillary investigations of biopsy specimens or serum samples.

This detailed clinical reference dataset will greatly facilitate diagnostic confirmation and management of patients, and has provided important insights into disease pathogenesis. Ann Neurol 2018;83:1105-1124.

Mutations in ANO5 are a frequent cause of undetermined muscular dystrophy, with both distal and proximal presentation. Other types include high hyperCKemia, myalgia, or calf hypertrophy over decades without significant weakness, especially in female patients. Mutations are distributed all over the gene, indicating that muscular dystrophy caused by ANO5 can be expected to occur in all populations.