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.

Heterozygous mutations in p63 are associated with split hand/foot malformations (SHFM), orofacial clefting, and ectodermal abnormalities. Elucidation of the p63 gene network that includes target genes and regulatory elements may reveal new genes for other malformation disorders. We performed genome-wide DNA–binding profiling by chromatin immunoprecipitation (ChIP), followed by deep sequencing (ChIP–seq) in primary human keratinocytes, and identified potential target genes and regulatory elements controlled by p63. We show that p63 binds to an enhancer element in the SHFM1 locus on chromosome 7q and that this element controls expression of DLX6 and possibly DLX5, both of which are important for limb development. A unique micro-deletion including this enhancer element, but not the DLX5/DLX6 genes, was identified in a patient with SHFM. Our study strongly indicates disruption of a non-coding cis-regulatory element located more than 250 kb from the DLX5/DLX6 genes as a novel disease mechanism in SHFM1. These data provide a proof-of-concept that the catalogue of p63 binding sites identified in this study may be of relevance to the studies of SHFM and other congenital malformations that resemble the p63-associated phenotypes.

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).

Dominant congenital spinal muscular atrophy (DCSMA) is a disorder of developing anterior horn cells and shows lower-limb predominance and clinical overlap with hereditary spastic paraplegia (HSP), a lower-limb-predominant disorder of corticospinal motor neurons. We have identified four mutations in bicaudal D homolog 2 (Drosophila) (BICD2) in six kindreds affected by DCSMA, DCSMA with upper motor neuron features, or HSP. BICD2 encodes BICD2, a key adaptor protein that interacts with the dynein-dynactin motor complex, which facilitates trafficking of cellular cargos that are critical to motor neuron development and maintenance. We demonstrate that mutations resulting in amino acid substitutions in two binding regions of BICD2 increase its binding affinity for the cytoplasmic dynein-dynactin complex, which might result in the perturbation of BICD2-dynein-dynactin-mediated trafficking, and impair neurite outgrowth. These findings provide insight into the mechanism underlying both the static and the slowly progressive clinical features and the motor neuron pathology that characterize BICD2-associated diseases, and underscore the importance of the dynein-dynactin transport pathway in the development and survival of both lower and upper motor neurons.

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.

See Cannon (doi: ) for a scientific commentary on this article. Dominant gain-of-function mutations in SCN4A , which encodes the α-subunit of the voltage-gated sodium channel, are a common cause of myotonia and periodic paralysis. Zaharieva et al. now report recessive loss-of-function SCN4A mutations in 11 patents with congenital myopathy. The mutations cause fully non-functional channels or result in reduced channel activity.

Spinal muscular atrophy is a disorder of lower motor neurons, most commonly caused by recessive mutations in SMN1 on chromosome 5q. Cases without SMN1 mutations are subclassified according to phenotype. Spinal muscular atrophy, lower extremity-predominant, is characterized by lower limb muscle weakness and wasting, associated with reduced numbers of lumbar motor neurons and is caused by mutations in DYNC1H1, which encodes a microtubule motor protein in the dynein-dynactin complex and one of its cargo adaptors, BICD2. We have now identified 32 patients with BICD2 mutations from nine different families, providing detailed insights into the clinical phenotype and natural history of BICD2 disease. BICD2 spinal muscular atrophy, lower extremity predominant most commonly presents with delayed motor milestones and ankle contractures. Additional features at presentation include arthrogryposis and congenital dislocation of the hips. In all affected individuals, weakness and wasting is lower-limb predominant, and typically involves both proximal and distal muscle groups. There is no evidence of sensory nerve involvement. Upper motor neuron signs are a prominent feature in a subset of individuals, including one family with exclusively adult-onset upper motor neuron features, consistent with a diagnosis of hereditary spastic paraplegia. In all cohort members, lower motor neuron features were static or only slowly progressive, and the majority remained ambulant throughout life. Muscle MRI in six individuals showed a common pattern of muscle involvement with fat deposition in most thigh muscles, but sparing of the adductors and semitendinosus. Muscle pathology findings were highly variable and included pseudomyopathic features, neuropathic features, and minimal change. The six causative mutations, including one not previously reported, result in amino acid changes within all three coiled-coil domains of the BICD2 protein, and include a possible 'hot spot' mutation, p.Ser107Leu present in four families. We used the recently solved crystal structure of a highly conserved region of the Drosophila orthologue of BICD2 to further-explore how the p.Glu774Gly substitution inhibits the binding of BICD2 to Rab6. Overall, the features of BICD2 spinal muscular atrophy, lower extremity predominant are consistent with a pathological process that preferentially affects lumbar lower motor neurons, with or without additional upper motor neuron involvement. Defining the phenotypic features in this, the largest BICD2 disease cohort reported to date, will facilitate focused genetic testing and filtering of next generation sequencing-derived variants in cases with similar features.

These data support NGS as a first-line tool for genetic evaluation of patients with a clinical phenotype suggestive of CMD, with muscle biopsy reserved as a second-tier investigation. Ann Neurol 2016;80:101-111.