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.

Timothy syndrome (TS) is a multisystem disorder that causes syncope and sudden death from cardiac arrhythmias. Prominent features include congenital heart disease, immune deficiency, intermittent hypoglycemia, cognitive abnormalities, and autism. All TS individuals have syndactyly (webbing of fingers and toes). We discovered that TS resulted from a recurrent, de novo cardiac L-type calcium channel (Ca _V 1.2) mutation, G406R. G406 is located in alternatively spliced exon 8A, encoding transmembrane segment S6 of domain I. Here, we describe two individuals with a severe variant of TS (TS2). Neither child had syndactyly. Both individuals had extreme prolongation of the QT interval on electrocardiogram, with a QT interval corrected for heart rate ranging from 620 to 730 ms, causing multiple arrhythmias and sudden death. One individual had severe mental retardation and nemaline rod skeletal myopathy. We identified de novo missense mutations in exon 8 of Ca _V 1.2 in both individuals. One was an analogous mutation to that found in exon 8A in classic TS, G406R. The other mutation was G402S. Exon 8 encodes the same region as exon 8A, and the two are mutually exclusive. The spliced form of Ca _V 1.2 containing exon 8 is highly expressed in heart and brain, accounting for ≈80% of Ca _V 1.2 mRNAs. G406R and G402S cause reduced channel inactivation, resulting in maintained depolarizing L-type calcium currents. Computer modeling showed prolongation of cardiomyocyte action potentials and delayed afterdepolarizations, factors that increase risk of arrhythmia. These data indicate that gain-of-function mutations of Ca _V 1.2 exons 8 and 8A cause distinct forms of TS.