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 current study characterizes a cohort of limb-girdle muscular dystrophy (LGMD) in the United States using whole-exome sequencing. Fifty-five families affected by LGMD were recruited using an institutionally approved protocol. Exome sequencing was performed on probands and selected parental samples. Pathogenic mutations and cosegregation patterns were confirmed by Sanger sequencing. Twenty-two families (40%) had novel and previously reported pathogenic mutations, primarily in LGMD genes, and also in genes for Duchenne muscular dystrophy, facioscapulohumeral muscular dystrophy, congenital myopathy, myofibrillar myopathy, inclusion body myopathy and Pompe disease. One family was diagnosed via clinical testing. Dominant mutations were identified in COL6A1, COL6A3, FLNC, LMNA, RYR1, SMCHD1 and VCP, recessive mutations in ANO5, CAPN3, GAA, LAMA2, SGCA and SGCG, and X-linked mutations in DMD. A previously reported variant in DMD was confirmed to be benign. Exome sequencing is a powerful diagnostic tool for LGMD. Despite careful phenotypic screening, pathogenic mutations were found in other muscle disease genes, largely accounting for the increased sensitivity of exome sequencing. Our experience suggests that broad sequencing panels are useful for these analyses because of the phenotypic overlap of many neuromuscular conditions. The confirmation of a benign DMD variant illustrates the potential of exome sequencing to help determine pathogenicity.

Mendel’s Law of equal segregation states that each of the two alleles in a diploid organism has an equal probability of being transmitted into...Meiotic drive occurs when a selfish element increases its transmission frequency above the Mendelian ratio by hijacking the asymmetric divisions of female meiosis. Meiotic drive causes genomic conflict and potentially has a major impact on genome evolution, but only a few drive loci of large effect have been described. New methods to reliably detect meiotic drive are therefore needed, particularly for discovering moderate-strength drivers that are likely to be more prevalent in natural populations than strong drivers. Here, we report an efficient method that uses sequencing of large pools of backcross (BC1) progeny to test for deviations from Mendelian segregation genome-wide with single-nucleotide polymorphisms (SNPs) that distinguish the parental strains. We show that meiotic drive can be detected by a characteristic pattern of decay in distortion of SNP frequencies, caused by recombination unlinking the driver from distal loci. We further show that control crosses allow allele-frequency distortion caused by meiotic drive to be distinguished from distortion resulting from developmental effects. We used this approach to test whether chromosomes with extreme telomere-length differences segregate at Mendelian ratios, as telomeric regions are a potential hotspot for meiotic drive due to their roles in meiotic segregation and multiple observations of high rates of telomere sequence evolution. Using four different pairings of long and short telomere strains, we find no evidence that extreme telomere-length variation causes meiotic drive in Drosophila. However, we identify one candidate meiotic driver in a centromere-linked region that shows an ∼8% increase in transmission frequency, corresponding to a ∼54:46 segregation ratio. Our results show that candidate meiotic drivers of moderate strength can be readily detected and localized in pools of BC1 progeny.

Mutations in MEGF10 cause early onset myopathy, areflexia, respiratory distress, and dysphagia (EMARDD), a rare congenital muscle disease, but the pathogenic mechanisms remain largely unknown. We demonstrate that short hairpin RNA (shRNA)-mediated knockdown of Megf10, as well as overexpression of the pathogenic human p.C774R mutation, leads to impaired proliferation and migration of C2C12 cells. Myoblasts from Megf10^−/− mice and Megf10^−/−/mdx double knockout (dko) mice also show impaired proliferation and migration compared to myoblasts from wild type and mdx mice, whereas the dko mice show histological abnormalities that are not observed in either single mutant mouse. Cell proliferation and migration are known to be regulated by the Notch receptor, which plays an essential role in myogenesis. Reciprocal co-immunoprecipitation studies show that Megf10 and Notch1 interact via their respective intracellular domains. These interactions are impaired by the pathogenic p.C774R mutation. Megf10 regulation of myoblast function appears to be mediated at least in part via interactions with key components of the Notch signaling pathway, and defects in these interactions may contribute to the pathogenesis of EMARDD.

Next-generation sequencing is commonly used to screen for pathogenic mutations in families with Mendelian disorders, but due to the pace of discoveries, gaps have widened for some diseases between genetic and pathophysiological knowledge. We recruited and analyzed 16 families with limb-girdle muscular dystrophy (LGMD) of Arab descent from Saudi Arabia and Sudan who did not have confirmed genetic diagnoses. The analysis included both traditional and next-generation sequencing approaches. Cellular and metabolic studies were performed on Pyroxd1 siRNA C2C12 myoblasts and controls. Pathogenic mutations were identified in eight of the 16 families. One Sudanese family of Arab descent residing in Saudi Arabia harbored a homozygous c.464A>G, p.Asn155Ser mutation in PYROXD1, a gene recently reported in association with myofibrillar myopathy and whose protein product reduces thiol residues. Pyroxd1 deficiency in murine C2C12 myoblasts yielded evidence for impairments of cellular proliferation, migration, and differentiation, while CG10721 (Pyroxd1 fly homolog) knockdown in Drosophila yielded a lethal phenotype. Further investigations indicated that Pyroxd1 does not localize to mitochondria, yet Pyroxd1 deficiency is associated with decreased cellular respiration. This study identified pathogenic mutations in half of the LGMD families from the cohort, including one in PYROXD1. Developmental impairments were demonstrated in vitro for Pyroxd1 deficiency and in vivo for CG10721 deficiency, with reduced metabolic activity in vitro for Pyroxd1 deficiency.

Background: Pathogenic mutations causing aberrant splicing are often difficult to detect. Standard variant analysis of next‐generation sequence (NGS) data focuses on canonical splice sites. Noncanonical splice sites are more difficult to ascertain.Methods: We developed a bioinformatics pipeline that screens existing NGS data for potentially aberrant novel essential splice sites (PANESS) and performed a pilot study on a family with a myotonic disorder. Further analyses were performed via qRT‐PCR, immunoblotting, and immunohistochemistry. RNAi knockdown studies were performed in Drosophila to model the gene deficiency.Results: The PANESS pipeline identified a homozygous ATP2A1 variant (NC_000016.9:g.28905928G>A; NM_004320.4:c.1287G>A:p.(Glu429=)) that was predicted to cause the omission of exon 11. Aberrant splicing of ATP2A1 was confirmed via qRT‐PCR, and abnormal expression of the protein product sarcoplasmic/endoplasmic reticulum Ca⁺⁺ ATPase 1 (SERCA1) was demonstrated in quadriceps femoris tissue from the proband. Ubiquitous knockdown of SERCA led to lethality in Drosophila, as did knockdown targeting differentiating or fusing myoblasts.Conclusions: This study confirms the potential of novel in silico algorithms to detect cryptic mutations in existing NGS data; expands the phenotypic spectrum of ATP2A1 mutations beyond classic Brody myopathy; and suggests that genetic testing of ATP2A1 should be considered in patients with clinical myotonia.