Identification of skewed X chromosome inactivation using exome and transcriptome sequencing in patients with suspected rare genetic disease.

BACKGROUND: X-chromosome inactivation (XCI) is an epigenetic process that occurs during early development in mammalian females by randomly silencing one of two copies of the X chromosome in each cell. The preferential inactivation of either the maternal or paternal copy of the X chromosome in a majority of cells results in a skewed or non-random pattern of X inactivation and is observed in over 25% of adult females. Identifying skewed X inactivation is of clinical significance in patients with suspected rare genetic diseases due to the possibility of biased expression of disease-causing genes present on the active X chromosome. The current clinical test for the detection of skewed XCI relies on the methylation status of the methylation-sensitive restriction enzyme (Hpall) binding site present in proximity of short tandem polymorphic repeats on the androgen receptor (AR) gene. This approach using one locus results in uninformative or inconclusive data for 10-20% of tests. Further, recent studies have shown inconsistency between methylation of the AR locus and the state of inactivation of the X chromosome. Herein, we develop a method for estimating X inactivation status, using exome and transcriptome sequencing data derived from blood in 227 female samples. We built a reference model for evaluation of XCI in 135 females from the GTEx consortium. We tested and validated the model on 11 female individuals with different types of undiagnosed rare genetic disorders who were clinically tested for X-skew using the AR gene assay and compared results to our outlier-based analysis technique. RESULTS: In comparison to the AR clinical test for identification of X inactivation, our method was concordant with the AR method in 9 samples, discordant in 1, and provided a measure of X inactivation in 1 sample with uninformative clinical results. We applied this method on an additional 81 females presenting to the clinic with phenotypes consistent with different hereditary disorders without a known genetic diagnosis. CONCLUSIONS: This study presents the use of transcriptome and exome sequencing data to provide an accurate and complete estimation of X-inactivation and skew status in a cohort of female patients with different types of suspected rare genetic disease.

Semiautomated approach focused on new genomic information results in time and effort-efficient reannotation of negative exome data.

Most rare disease patients (75-50%) undergoing genomic sequencing remain unsolved, often due to lack of information about variants identified. Data review over time can leverage novel information regarding disease-causing variants and genes, increasing this diagnostic yield. However, time and resource constraints have limited reanalysis of genetic data in clinical laboratories setting. We developed RENEW, (REannotation of NEgative WES/WGS) an automated reannotation procedure that uses relevant new information in on-line genomic databases to enable rapid review of genomic findings. We tested RENEW in an unselected cohort of 1066 undiagnosed cases with a broad spectrum of phenotypes from the Mayo Clinic Center for Individualized Medicine using new information in ClinVar, HGMD and OMIM between the date of previous analysis/testing and April of 2022. 5741 variants prioritized by RENEW were rapidly reviewed by variant interpretation specialists. Mean analysis time was approximately 20 s per variant (32 h total time). Reviewed cases were classified as: 879 (93.0%) undiagnosed, 63 (6.6%) putatively diagnosed, and 4 (0.4%) definitively diagnosed. New strategies are needed to enable efficient review of genomic findings in unsolved cases. We report on a fast and practical approach to address this need and improve overall diagnostic success in patient testing through a recurrent reannotation process.

Bi-allelic variants in INTS11 are associated with a complex neurological disorder.

The Integrator complex is a multi-subunit protein complex that regulates the processing of nascent RNAs transcribed by RNA polymerase II (RNAPII), including small nuclear RNAs, enhancer RNAs, telomeric RNAs, viral RNAs, and protein-coding mRNAs. Integrator subunit 11 (INTS11) is the catalytic subunit that cleaves nascent RNAs, but, to date, mutations in this subunit have not been linked to human disease. Here, we describe 15 individuals from 10 unrelated families with bi-allelic variants in INTS11 who present with global developmental and language delay, intellectual disability, impaired motor development, and brain atrophy. Consistent with human observations, we find that the fly ortholog of INTS11, dIntS11, is essential and expressed in the central nervous systems in a subset of neurons and most glia in larval and adult stages. Using Drosophila as a model, we investigated the effect of seven variants. We found that two (p.Arg17Leu and p.His414Tyr) fail to rescue the lethality of null mutants, indicating that they are strong loss-of-function variants. Furthermore, we found that five variants (p.Gly55Ser, p.Leu138Phe, p.Lys396Glu, p.Val517Met, and p.Ile553Glu) rescue lethality but cause a shortened lifespan and bang sensitivity and affect locomotor activity, indicating that they are partial loss-of-function variants. Altogether, our results provide compelling evidence that integrity of the Integrator RNA endonuclease is critical for brain development.

Cre/lox regulated conditional rescue and inactivation with zebrafish UFlip alleles generated by CRISPR-Cas9 targeted integration.

The ability to regulate gene activity spatially and temporally is essential to investigate cell-type-specific gene function during development and in postembryonic processes and disease models. The Cre/lox system has been widely used for performing cell and tissue-specific conditional analysis of gene function in zebrafish. However, simple and efficient methods for isolation of stable, Cre/lox regulated zebrafish alleles are lacking. Here, we applied our GeneWeld CRISPR-Cas9 targeted integration strategy to generate floxed alleles that provide robust conditional inactivation and rescue. A universal targeting vector, UFlip, with sites for cloning short homology arms flanking a floxed 2A-mRFP gene trap, was integrated into an intron in rbbp4 and rb1. rbbp4off and rb1off integration alleles resulted in strong mRFP expression,>99% reduction of endogenous gene expression, and recapitulated known indel loss-of-function phenotypes. Introduction of Cre led to stable inversion of the floxed cassette, loss of mRFP expression, and phenotypic rescue. rbbp4on and rb1on integration alleles did not cause phenotypes in combination with a loss-of-function mutation. Addition of Cre led to conditional inactivation by stable inversion of the cassette, gene trapping and mRFP expression, and the expected mutant phenotype. Neural progenitor Cre drivers were used for conditional inactivation and phenotypic rescue to showcase how this approach can be used in specific cell populations. Together these results validate a simplified approach for efficient isolation of Cre/lox-responsive conditional alleles in zebrafish. Our strategy provides a new toolkit for generating genetic mosaics and represents a significant advance in zebrafish genetics.

Heterozygous variants in PRPF8 are associated with neurodevelopmental disorders.

The pre-mRNA-processing factor 8, encoded by PRPF8, is a scaffolding component of a spliceosome complex involved in the removal of introns from mRNA precursors. Previously, heterozygous pathogenic variants in PRPF8 have been associated with autosomal dominant retinitis pigmentosa. More recently, PRPF8 was suggested as a candidate gene for autism spectrum disorder due to the enrichment of sequence variants in this gene in individuals with neurodevelopmental disorders. We report 14 individuals with various forms of neurodevelopmental conditions, found to have heterozygous, predominantly de novo, missense, and loss-of-function variants in PRPF8. These individuals have clinical features that may represent a new neurodevelopmental syndrome.

Rbbp4 loss disrupts neural progenitor cell cycle regulation independent of Rb and leads to Tp53 acetylation and apoptosis.

BACKGROUND: Retinoblastoma binding protein 4 (Rbbp4) is a component of transcription regulatory complexes that control cell cycle gene expression. Previous work indicated that Rbbp4 cooperates with the Rb tumor suppressor to block cell cycle entry. Here, we use genetic analysis to examine the interactions of Rbbp4, Rb, and Tp53 in zebrafish neural progenitor cell cycle regulation and survival. RESULTS: Rbbp4 is upregulated across the spectrum of human embryonal and glial brain cancers. Transgenic rescue of rbbp4 mutant embryos shows Rbbp4 is essential for zebrafish neurogenesis. Rbbp4 loss leads to apoptosis and γ-H2AX in the developing brain that is suppressed by tp53 knockdown or maternal zygotic deletion. Mutant retinal neural precursors accumulate in M phase and fail to initiate G0 gene expression. rbbp4; rb1 mutants show an additive effect on the number of M phase cells. In rbbp4 mutants, Tp53 acetylation is detected; however, Rbbp4 overexpression did not rescue DNA damage-induced apoptosis. CONCLUSION: Rbbp4 is necessary for neural progenitor cell cycle progression and initiation of G0 independent of Rb. Tp53-dependent apoptosis in the absence of Rbpb4 correlates with Tp53 acetylation. Together these results suggest that Rbbp4 is required for cell cycle exit and contributes to neural progenitor survival through the regulation of Tp53 acetylation.

Childhood-onset dystonia-causing KMT2B variants result in a distinctive genomic hypermethylation profile.

BACKGROUND: Dystonia is a clinically and genetically heterogeneous movement disorder characterized by sustained or intermittent muscle contractions causing abnormal, often repetitive, movements and/or postures. Heterozygous variants in lysine methyltransferase 2B (KMT2B), encoding a histone H3 methyltransferase, have been associated with a childhood-onset, progressive and complex form of dystonia (dystonia 28, DYT28). Since 2016, more than one hundred rare KMT2B variants have been reported, including frameshift, nonsense, splice site, missense and other in-frame changes, many having an uncertain clinical impact. RESULTS: We characterize the genome-wide peripheral blood DNA methylation profiles of a cohort of 18 patients with pathogenic and unclassified KMT2B variants. We resolve the "episignature" associated with KMT2B haploinsufficiency, proving that this approach is robust in diagnosing clinically unsolved cases, properly classifying them with respect to other partially overlapping dystonic phenotypes, other rare neurodevelopmental disorders and healthy controls. Notably, defective KMT2B function in DYT28 causes a non-random DNA hypermethylation across the genome, selectively involving promoters and other regulatory regions positively controlling gene expression. CONCLUSIONS: We demonstrate a distinctive DNA hypermethylation pattern associated with DYT28, provide an epigenetic signature for this disorder enabling accurate diagnosis and reclassification of ambiguous genetic findings and suggest potential therapeutic approaches.

CSNK2B: A broad spectrum of neurodevelopmental disability and epilepsy severity.

CSNK2B has recently been implicated as a disease gene for neurodevelopmental disability (NDD) and epilepsy. Information about developmental outcomes has been limited by the young age and short follow-up for many of the previously reported cases, and further delineation of the spectrum of associated phenotypes is needed. We present 25 new patients with variants in CSNK2B and refine the associated NDD and epilepsy phenotypes. CSNK2B variants were identified by research or clinical exome sequencing, and investigators from different centers were connected via GeneMatcher. Most individuals had developmental delay and generalized epilepsy with onset in the first 2 years. However, we found a broad spectrum of phenotypic severity, ranging from early normal development with pharmacoresponsive seizures to profound intellectual disability with intractable epilepsy and recurrent refractory status epilepticus. These findings suggest that CSNK2B should be considered in the diagnostic evaluation of patients with a broad range of NDD with treatable or intractable seizures.

Truncating SRCAP variants outside the Floating-Harbor syndrome locus cause a distinct neurodevelopmental disorder with a specific DNA methylation signature.

Truncating variants in exons 33 and 34 of the SNF2-related CREBBP activator protein (SRCAP) gene cause the neurodevelopmental disorder (NDD) Floating-Harbor syndrome (FLHS), characterized by short stature, speech delay, and facial dysmorphism. Here, we present a cohort of 33 individuals with clinical features distinct from FLHS and truncating (mostly de novo) SRCAP variants either proximal (n = 28) or distal (n = 5) to the FLHS locus. Detailed clinical characterization of the proximal SRCAP individuals identified shared characteristics: developmental delay with or without intellectual disability, behavioral and psychiatric problems, non-specific facial features, musculoskeletal issues, and hypotonia. Because FLHS is known to be associated with a unique set of DNA methylation (DNAm) changes in blood, a DNAm signature, we investigated whether there was a distinct signature associated with our affected individuals. A machine-learning model, based on the FLHS DNAm signature, negatively classified all our tested subjects. Comparing proximal variants with typically developing controls, we identified a DNAm signature distinct from the FLHS signature. Based on the DNAm and clinical data, we refer to the condition as "non-FLHS SRCAP-related NDD." All five distal variants classified negatively using the FLHS DNAm model while two classified positively using the proximal model. This suggests divergent pathogenicity of these variants, though clinically the distal group presented with NDD, similar to the proximal SRCAP group. In summary, for SRCAP, there is a clear relationship between variant location, DNAm profile, and clinical phenotype. These results highlight the power of combined epigenetic, molecular, and clinical studies to identify and characterize genotype-epigenotype-phenotype correlations.

SPEN haploinsufficiency causes a neurodevelopmental disorder overlapping proximal 1p36 deletion syndrome with an episignature of X chromosomes in females.

Deletion 1p36 (del1p36) syndrome is the most common human disorder resulting from a terminal autosomal deletion. This condition is molecularly and clinically heterogeneous. Deletions involving two non-overlapping regions, known as the distal (telomeric) and proximal (centromeric) critical regions, are sufficient to cause the majority of the recurrent clinical features, although with different facial features and dysmorphisms. SPEN encodes a transcriptional repressor commonly deleted in proximal del1p36 syndrome and is located centromeric to the proximal 1p36 critical region. Here, we used clinical data from 34 individuals with truncating variants in SPEN to define a neurodevelopmental disorder presenting with features that overlap considerably with those of proximal del1p36 syndrome. The clinical profile of this disease includes developmental delay/intellectual disability, autism spectrum disorder, anxiety, aggressive behavior, attention deficit disorder, hypotonia, brain and spine anomalies, congenital heart defects, high/narrow palate, facial dysmorphisms, and obesity/increased BMI, especially in females. SPEN also emerges as a relevant gene for del1p36 syndrome by co-expression analyses. Finally, we show that haploinsufficiency of SPEN is associated with a distinctive DNA methylation episignature of the X chromosome in affected females, providing further evidence of a specific contribution of the protein to the epigenetic control of this chromosome, and a paradigm of an X chromosome-specific episignature that classifies syndromic traits. We conclude that SPEN is required for multiple developmental processes and SPEN haploinsufficiency is a major contributor to a disorder associated with deletions centromeric to the previously established 1p36 critical regions.

Impact of integrated translational research on clinical exome sequencing.

PURPOSE: Exome sequencing often identifies pathogenic genetic variants in patients with undiagnosed diseases. Nevertheless, frequent findings of variants of uncertain significance necessitate additional efforts to establish causality before reaching a conclusive diagnosis. To provide comprehensive genomic testing to patients with undiagnosed disease, we established an Individualized Medicine Clinic, which offered clinical exome testing and included a Translational Omics Program (TOP) that provided variant curation, research activities, or research exome sequencing. METHODS: From 2012 to 2018, 1101 unselected patients with undiagnosed diseases received exome testing. Outcomes were reviewed to assess impact of the TOP and patient characteristics on diagnostic rates through descriptive and multivariate analyses. RESULTS: The overall diagnostic yield was 24.9% (274 of 1101 patients), with 174 (15.8% of 1101) diagnosed on the basis of clinical exome sequencing alone. Four hundred twenty-three patients with nondiagnostic or without access to clinical exome sequencing were evaluated by the TOP, with 100 (9% of 1101) patients receiving a diagnosis, accounting for 36.5% of the diagnostic yield. The identification of a genetic diagnosis was influenced by the age at time of testing and the disease phenotype of the patient. CONCLUSION: Integration of translational research activities into clinical practice of a tertiary medical center can significantly increase the diagnostic yield of patients with undiagnosed disease.

Novel loss-of-function variants in TRIO are associated with neurodevelopmental disorder: case report.

BACKGROUND: Damaging variants in TRIO have been associated with moderate to severe neurodevelopmental disorders in humans. While recent work has delineated the positional effect of missense variation on the resulting phenotype, the clinical spectrum associated with loss-of-function variation has yet to be fully defined. CASE PRESENTATION: We report on two probands with novel loss-of-function variants in TRIO. Patient 1 presents with a severe neurodevelopmental disorder and macrocephaly. The TRIO variant is inherited from his affected mother. Patient 2 presents with moderate developmental delays, microcephaly, and cutis aplasia with a frameshift variant of unknown inheritance. CONCLUSIONS: We describe two patients with neurodevelopmental disorder, macro/microcephaly, and cutis aplasia in one patient. Both patients have loss-of-function variants, helping to further characterize how these types of variants affect the phenotypic spectrum associated with TRIO. We also present the third reported case of autosomal dominant inheritance of a damaging variant in TRIO.

Haploinsufficiency as a disease mechanism in GNB1-associated neurodevelopmental disorder.

BACKGROUND: GNB1 encodes a subunit of a heterotrimeric G-protein complex that transduces intracellular signaling cascades. Disruptions to the gene have previously been shown to be embryonic lethal in knockout mice and to cause complex neurodevelopmental disorders in humans. To date, the majority of variants associated with disease in humans have been missense variants in exons 5-7. METHODS: Genetic sequencing was performed on two patients presenting with complex neurological phenotypes including intellectual disability, hypotonia, and in one patient seizures. Reported variants were assessed using RNA sequencing and functional BRET/BiFC assays. RESULTS: A splice variant reported in patient 1 was confirmed to cause usage of a cryptic splice site leading to a truncated protein product. Patient 2 was reported to have a truncating variant. BRET and BiFC assays of both patient variants confirmed both were deficient in inducing GPCR-induced G protein activation due to lack of dimer formation with the Gγ subunit. CONCLUSION: Here, we report two patients with functionally confirmed loss of function variants in GNB1 and neurodevelopmental phenotypes including intellectual disability, hypotonia, and seizures in one patient. These results suggest haploinsufficiency of GNB1 is a mechanism for neurodevelopmental disorders in humans.

SPECC1L regulates palate development downstream of IRF6.

SPECC1L mutations have been identified in patients with rare atypical orofacial clefts and with syndromic cleft lip and/or palate (CL/P). These mutations cluster in the second coiled-coil and calponin homology domains of SPECC1L and severely affect the ability of SPECC1L to associate with microtubules. We previously showed that gene-trap knockout of Specc1l in mouse results in early embryonic lethality. We now present a truncation mutant mouse allele, Specc1lΔC510, that results in perinatal lethality. Specc1lΔC510/ΔC510 homozygotes showed abnormal palate rugae but did not show cleft palate. However, when crossed with a gene-trap allele, Specc1lcGT/ΔC510 compound heterozygotes showed a palate elevation delay with incompletely penetrant cleft palate. Specc1lcGT/ΔC510 embryos exhibit transient oral epithelial adhesions at E13.5, which may delay shelf elevation. Consistent with oral adhesions, we show periderm layer abnormalities, including ectopic apical expression of adherens junction markers, similar to Irf6 hypomorphic mutants and Arhgap29 heterozygotes. Indeed, SPECC1L expression is drastically reduced in Irf6 mutant palatal shelves. Finally, we wanted to determine if SPECC1L deficiency also contributed to non-syndromic (ns) CL/P. We sequenced 62 Caucasian, 89 Filipino, 90 Ethiopian, 90 Nigerian and 95 Japanese patients with nsCL/P and identified three rare coding variants (p.Ala86Thr, p.Met91Iso and p.Arg546Gln) in six individuals. These variants reside outside of SPECC1L coiled-coil domains and result in milder functional defects than variants associated with syndromic clefting. Together, our data indicate that palate elevation is sensitive to deficiency of SPECC1L dosage and function and that SPECC1L cytoskeletal protein functions downstream of IRF6 in palatogenesis.

Protein-elongating mutations in MYH11 are implicated in a dominantly inherited smooth muscle dysmotility syndrome with severe esophageal, gastric, and intestinal disease.

Gastrointestinal motility disorders include a spectrum of mild to severe clinical phenotypes that are caused by smooth muscle dysfunction. We investigated the genetic etiology of severe esophageal, gastric, and colonic dysmotility in two unrelated families with autosomal dominant disease presentation. Using exome sequencing, we identified a 2 base pair insertion at the end of the myosin heavy chain 11 (MYH11) gene in all affected members of Family 1 [NM_001040113:c.5819_5820insCA(p.Gln1941Asnfs*91)] and a 1 base pair deletion at the same genetic locus in Proband 2 [NM_001040113:c.5819del(p.Pro1940Hisfs*91)]. Both variants are predicted to result in a similarly elongated protein product. Heterozygous dominant negative MYH11 pathogenic variants have been associated with thoracic aortic aneurysm and dissection while biallelic null alleles have been associated with megacystis microcolon intestinal hypoperistalsis syndrome. This report highlights heterozygous protein-elongating MYH11 variants affecting the SM2 isoforms of MYH11 as a cause for severe gastrointestinal dysmotility, and we hypothesize that the mechanistic pathogenesis of this disease, dominant hypercontractile loss-of-function, is distinct from those implicated in other diseases involving MYH11 dysfunction.

Clinical and Biochemical Phenotypes in a Family With ENPP1 Mutations.

Inactivating mutations of the ENPP1 gene are associated with generalized arterial calcification of infancy (GACI) and less often autosomal-recessive hypophosphatemic rickets type 2 (ARHR2). We aimed to investigate the spectrum of phenotypes in a family with monoallelic and biallelic mutations of ENPP1 after identification through whole exome sequencing of a 54-year-old female with biallelic mutation of ENPP1, c.323G > T; p.Cys108Phe and c.1441C > T; p.Arg481Trp. Including the proband, 2 subjects had biallelic mutations, 5 had monoallelic mutations, and 2 had no mutation of ENPP1. The maternal mutation, a known pathogenic variant associated with GACI, was found in 3 subjects with monoallelic mutations, while the paternal mutation, which was not previously reported, was present in 2 subjects with monoallelic mutations. Both subjects with biallelic mutations had bowing of bilateral femurs, periarticular mineral deposition, normocalcemic primary hyperparathyroidism with multigland parathyroidectomy, increased carotid intima-media thickness, and enthesopathy was also noted in one subject. Intact FGF23 was elevated in both subjects with biallelic mutations, while C-terminal FGF23 was only elevated in one and PPi was reduced in one. Subjects with monoallelic mutations did not have periarticular calcifications or bone deformities. To conclude, patients with biallelic GACI causing mutations can survive well into adulthood, and despite the same biallelic ENPP1 pathogenic variants, clinical and biochemical manifestations can significantly differ, and include enthesopathy and primary hyperparathyroidism, which have not been previously described. Although carriers of monoallelic ENPP1 variants appear unaffected by classic disease manifestations, there may be subtle biochemical and clinical findings that warrant further investigation. © 2019 American Society for Bone and Mineral Research.

Novel biallelic variants in MSTO1 associated with mitochondrial myopathy.

Mitochondrial disorders are caused by nuclear and mitochondrial pathogenic variants leading to defects in mitochondrial function and cellular respiration. Recently, the nuclear-encoded mitochondrial fusion gene MSTO1 (Misato 1) has been implicated in mitochondrial myopathy and ataxia. Here we report on a 30-yr-old man presenting with a maternally inherited NM_018116.3:c.651C>G, p.F217L missense variant as well as a paternally inherited arr[GRCh37] 1q22(155581773_155706887) × 1 deletion encompassing exons 7-14 of MSTO1 His phenotype included muscle weakness, hypotonia, early motor developmental delay, pectus excavatum, and scoliosis. Testing revealed elevated plasma creatine kinase, and electromyogram results were consistent with longstanding generalized myopathy. These phenotypic features overlap well with previously reported patients harboring biallelic MSTO1 variants. Additionally, our patient presents with dysphagia and restrictive lung disease, not previously reported for MSTO1-associated disorders. The majority of patients with disease-associated variants in MSTO1 present with biallelic variants suggesting autosomal recessive inheritance; however, one family has been reported with a single variant and presumed autosomal dominant inheritance. The pattern of inheritance we observed is consistent with the majority of previous reports suggesting an autosomal recessive disorder. We add to our knowledge of the syndrome caused by variants in MSTO1 and provide additional evidence supporting autosomal recessive inheritance. We also describe phenotypic features not reported in previous cases, although further research is needed to confirm they are associated with defects in MSTO1.

A tailored approach to fusion transcript identification increases diagnosis of rare inherited disease.

BACKGROUND: RNA sequencing has been proposed as a means of increasing diagnostic rates in studies of undiagnosed rare inherited disease. Recent studies have reported diagnostic improvements in the range of 7.5-35% by profiling splicing, gene expression quantification and allele specific expression. To-date however, no study has systematically assessed the presence of gene-fusion transcripts in cases of germline disease. Fusion transcripts are routinely identified in cancer studies and are increasingly recognized as having diagnostic, prognostic or therapeutic relevance. Isolated reports exist of fusion transcripts being detected in cases of developmental and neurological phenotypes, and thus, systematic application of fusion detection to germline conditions may further increase diagnostic rates. However, current fusion detection methods are unsuited to the investigation of germline disease due to performance biases arising from their development using tumor, cell-line or in-silico data. METHODS: We describe a tailored approach to fusion candidate identification and prioritization in a cohort of 47 undiagnosed, suspected inherited disease patients. We modify an existing fusion transcript detection algorithm by eliminating its cell line-derived filtering steps, and instead, prioritize candidates using a custom workflow that integrates genomic and transcriptomic sequence alignment, biological and technical annotations, customized categorization logic, and phenotypic prioritization. RESULTS: We demonstrate that our approach to fusion transcript identification and prioritization detects genuine fusion events excluded by standard analyses and efficiently removes phenotypically unimportant candidates and false positive events, resulting in a reduced candidate list enriched for events with potential phenotypic relevance. We describe the successful genetic resolution of two previously undiagnosed disease cases through the detection of pathogenic fusion transcripts. Furthermore, we report the experimental validation of five additional cases of fusion transcripts with potential phenotypic relevance. CONCLUSIONS: The approach we describe can be implemented to enable the detection of phenotypically relevant fusion transcripts in studies of rare inherited disease. Fusion transcript detection has the potential to increase diagnostic rates in rare inherited disease and should be included in RNA-based analytical pipelines aimed at genetic diagnosis.