Search

As patient derived xenograft (PDX) models are increasingly used for preclinical drug development, strategies to account for the nonhuman component of PDX RNA expression data are critical to its interpretation. A bioinformatics pipeline to separate donor tumor and mouse stroma transcriptome profiles was devised and tested. To examine the molecular fidelity of PDX versus donor tumors, we compared mRNA differences between paired PDX-donor tumors from nine ovarian cancer patients. 1,935 differentially expressed genes were identified between PDX and donor tumors. Over 90% (n = 1767) of these genes were down-regulated in PDX models and enriched in stroma-specific functions. Several protein kinases were also differentially expressed in PDX tumors, e.g. PDGFRA, PDGFRB and CSF1R. Upon in silico removal of these PDX-donor tumor differentially expressed genes, a stronger transcriptional resemblance between PDX-donor tumor pairs was seen (average correlation coefficient increases from 0.91 to 0.95). We devised and validated an effective bioinformatics strategy to separate mouse stroma expression from human tumor expression for PDX RNAseq. In addition, we showed most of the PDX-donor differentially expressed genes were implicated in stromal components. The molecular similarities and differences between PDX and donor tumors have implications in future therapeutic trial designs and treatment response evaluations using PDX models.

PURPOSE: Sudden infant death syndrome (SIDS) is the commonest cause of sudden death of an infant; however, the genetic basis remains poorly understood. We aimed to identify noncardiac genes underpinning SIDS and determine their prevalence compared with ethnically matched controls.METHODS: Using exome sequencing we assessed the yield of ultrarare nonsynonymous variants (minor allele frequency [MAF] ≤0.00005, dominant model; MAF ≤0.01, recessive model) in 278 European SIDS cases (62% male; average age =2.7 ± 2 months) versus 973 European controls across 61 noncardiac SIDS-susceptibility genes. The variants were classified according to American College of Medical Genetics and Genomics criteria. Case-control, gene-collapsing analysis was performed in eight candidate biological pathways previously implicated in SIDS pathogenesis.RESULTS: Overall 43/278 SIDS cases harbored an ultrarare single-nucleotide variant compared with 114/973 controls (15.5 vs. 11.7%, p=0.10). Only 2/61 noncardiac genes were significantly overrepresented in cases compared with controls (ECE1, 3/278 [1%] vs. 1/973 [0.1%] p=0.036; SLC6A4, 2/278 [0.7%] vs. 1/973 [0.1%] p=0.049). There was no difference in yield of pathogenic or likely pathogenic variants between cases and controls (1/278 [0.36%] vs. 4/973 [0.41%]; p=1.0). Gene-collapsing analysis did not identify any specific biological pathways to be significantly associated with SIDS.CONCLUSIONS: A monogenic basis for SIDS amongst the previously implicated noncardiac genes and their encoded biological pathways is negligible.

Angiosarcoma is an aggressive malignancy of endothelial cells that carries a high mortality rate. Cytotoxic chemotherapy can elicit clinical responses, but the duration of response is limited. Sequencing reveals multiple mutations in angiogenesis pathways in angiosarcomas, particularly in vascular endothelial growth factor (VEGFR) and mitogen-activated protein kinase (MAPK) signaling. We aimed to determine the biological relevance of these pathways in angiosarcoma. Tissue microarray consisting of clinical formalin-fixed paraffin embedded tissue archival samples were stained for phospho- extracellular signal-regulated kinase (p-ERK) with immunohistochemistry. Angiosarcoma cell lines were treated with the mitogen-activated protein kinase kinase (MEK) inhibitor trametinib, pan-VEGFR inhibitor cediranib, or combined trametinib and cediranib and viability was assessed. Reverse phase protein array (RPPA) was performed to assess multiple oncogenic protein pathways. SVR angiosarcoma cells were grown in vivo and gene expression effects of treatment were assessed with whole exome RNA sequencing. MAPK signaling was found active in over half of clinical angiosarcoma samples. Inhibition of MAPK signaling with the MEK inhibitor trametinib decreased the viability of angiosarcoma cells. Combined inhibition of the VEGF and MAPK pathways with cediranib and trametinib had an additive effect in in vitro models, and a combinatorial effect in an in vivo model. Combined treatment led to smaller tumors than treatment with either agent alone. RNA-seq demonstrated distinct expression signatures between the trametinib treated tumors and those treated with both trametinib and cediranib. These results indicate a clinical study of combined VEGFR and MEK inhibition in angiosarcoma is warranted.

Background: DNA polymerase epsilon (POLE) is encoded by the POLE gene, and POLE-driven tumors are characterized by high mutational rates. POLE-driven tumors are relatively common in endometrial and colorectal cancer, and their presence is increasingly recognized in ovarian cancer (OC) of endometrioid type. POLE-driven cases possess an abundance of TCT>TAT and TCG>TTG somatic mutations characterized by mutational signature 10 from the Catalog of Somatic Mutations in Cancer (COSMIC). By quantifying the contribution of COSMIC mutational signature 10 in RNA sequencing (RNA-seq) we set out to identify POLE-driven tumors in a set of unselected Mayo Clinic OC. Methods: Mutational profiles were calculated using expressed single-nucleotide variants (eSNV) in the Mayo Clinic OC tumors (n=195), The Cancer Genome Atlas (TCGA) OC tumors (n=419), and the Genotype-Tissue Expression (GTEx) normal ovarian tissues (n=84). Non-negative Matrix Factorization (NMF) of the mutational profiles inferred the contribution per sample of four distinct mutational signatures, one of which corresponds to COSMIC mutational signature 10. Results: In the Mayo Clinic OC cohort we identified six tumors with a predicted contribution from COSMIC mutational signature 10 of over five mutations per megabase. These six cases harbored known POLE hotspot mutations (P286R, S297F, V411L, and A456P) and were of endometrioid histotype (P=5e-04). These six tumors were hypermutated with a higher tumor mutation load (mean, 54.02 mutations per megabase) compared to non-POLE endometrioid OC cases (mean, 7.69 mutations per megabase; P=5e-04), and had an early onset (average age of patients at onset, 48.33 years) when compared to non-POLE endometrioid OC cohort (average age at onset, 60.13 years; P=.008). Samples from TCGA and GTEx had a low COSMIC signature 10 contribution (median 0.16 mutations per megabase; maximum 1.78 mutations per megabase) and carried no POLE hotspot mutations.Conclusions: From the largest cohort of RNA-seq from endometrioid OC to date (n=53), we identified six hypermutated samples likely driven by POLE (frequency, 11%). Our result suggests the clinical need to screen for POLE driver mutations in endometrioid OC, which can guide enrollment in immunotherapy clinical trials.

Background: Sudden infant death syndrome (SIDS) is a leading cause of postneonatal mortality. Genetic heart diseases (GHDs) underlie some cases of SIDS. Objectives: This study aimed to determine the spectrum and prevalence of GHD-associated mutations as a potential monogenic basis for SIDS. Methods: A cohort of 419 unrelated SIDS cases (257 male; average age 2.7 ± 1.9 months) underwent whole exome sequencing and a targeted analysis of 90 GHD-susceptibility genes. The yield of “potentially informative,” ultra-rare variants (minor allele frequency

BACKGROUND: Transplantation of autologous mesenchymal stem cells (MSCs) may be a viable option for treatment of several diseases. MSCs efficacy depends on adequate function of their mitochondria, which might be impaired in a noxious milieu. OBJECTIVES: We hypothesized that obesity compromises MSCs mitochondrial structure and function, possibly via micro-RNA (miRNA)-based mechanisms. METHODS: MSCs were collected from swine abdominal adipose tissue after 16 weeks of Lean or Obese diet (n=7 each). Mitochondrial structure was assessed by electron microscopy and function by membrane potential and cytochrome-c oxidase (COX)-IV activity. Oxidative stress was assessed by Mito-SOX and dihydroethidium staining. Next-generation sequencing (RNA-seq) was performed to identify miRNAs expression in MSCs, and predicted mitochondrial target genes were then identified (MitoCarta). RESULTS: Compared to Lean-MSCs, mitochondria from Obese-MSCs were smaller and showed cristae remodeling and loss. Mitochondrial membrane potential and COX-IV activity decreased in Obese-MSCs, associated with increased mitochondrial oxidative stress. RNA-seq generated reads for 413 miRNAs, of which 5 miRNAs were upregulated in Obese-MSCs (fold change >2, p

Autologous transplantation of mesenchymal stem cells (MSCs) is a viable option for the treatment of several diseases. Evidence indicates that MSCs release extracellular vesicles (EVs) and that EVs shuttle miRNAs to damaged parenchymal cells to activate an endogenous repair program. We hypothesize that comorbidities may interfere with the packaging of cargo in MSC-derived EVs. Therefore, we examined whether metabolic syndrome (MetS) modulates the miRNA content packed within MSC-derived EVs. MSCs were collected from swine abdominal adipose tissue after 16 weeks of lean or obese diet (n = 7 each). Next-generation RNA sequencing of miRNAs (miRNA-seq) was performed to identify miRNAs enriched in MSC-derived EVs and their predicted target genes. Functional pathway analysis of the top 50 target genes of the top 4 miRNAs enriched in each group was performed using gene ontology analysis. Lean- and MetS-EVs were enriched in, respectively, 14 and 8 distinct miRNAs. Target genes of miRNAs enriched in MetS-EVs were implicated in the development of MetS and its complications, including diabetes-related pathways, validated transcriptional targets of AP1 family members Fra1 and Fra2, Class A/1 (Rhodopsin-like receptors), and Peptide ligand-binding receptors. In contrast, miRNAs enriched in Lean EVs target primarily EphrinA-EPHA and the Rho family of GTPases. MetS alters the miRNA content of EVs derived from porcine adipose tissue MSCs. These alterations could impair the efficacy and limit the therapeutic use of autologous MSCs in subjects with MetS. Our findings may assist in developing adequate regenerative strategies to preserve the reparative potency of MSCs in individuals with MetS. © 2017 International Society for Advancement of Cytometry.

Introduction We hypothesized that common Alzheimer's disease (AD)-associated variants within the triggering receptor expressed on myeloid ( TREM ) gene cluster influence disease through gene expression. Methods Expression microarrays on temporal cortex and cerebellum from ∼400 neuropathologically diagnosed subjects and two independent RNAseq replication cohorts were used for expression quantitative trait locus analysis. Results A variant within a DNase hypersensitive site 5′ of TREM2 , rs9357347-C, associates with reduced AD risk and increased TREML1 and TREM2 levels (uncorrected P = 6.3 × 10 −3 and 4.6 × 10 −2 , respectively). Meta-analysis on expression quantitative trait locus results from three independent data sets ( n = 1006) confirmed these associations (uncorrected P = 3.4 × 10 −2 and 3.5 × 10 −3 , Bonferroni-corrected P = 6.7 × 10 −2 and 7.1 × 10 −3 , respectively). Discussion Our findings point to rs9357347 as a functional regulatory variant that contributes to a protective effect observed at the TREM locus in the International Genomics of Alzheimer's Project genome-wide association study meta-analysis and suggest concomitant increase in TREML1 and TREM2 brain levels as a potential mechanism for protection from AD.

Introduction: Acquired resistance to platinum remains a major clinical problem in ovarian cancer since systemic therapies have limited efficacy in this setting. There is growing evidence that the tumor microenvironment (TME) contributes to such resistance, as well as the pathogenesis of the disease. The study of the TME is limited by technical challenges of separating the epithelial from the stromal component. To overcome this limitation, we performed RNA sequencing (RNA-Seq) for three pairs of sensitive and induced resistant patient-derived xenograft models since the epithelial transcripts originate from human tissue and the stroma component is murine. Methods: To determine which TME genes might be contributing to chemotherapy resistance, three platinum-sensitive ovarian cancer PDX lines were treated repeatedly with carboplatin/paclitaxel until resistance occurred. Primary platinum-sensitive and induced platinum resistant tumors were characterized histologically by H&E and Masson Trichrome staining. The three pairs of sensitive and resistant PDX models were then subjected to RNA-Seq. Raw sequencing data were processed through Xenome to classify the reads into epithelial carcinoma (human) or TME (mouse) through Mayo Clinic’s internal RNA-Seq pipeline for transcriptome quantification. Differential expression analysis between induced resistant and sensitive models was carried out for both carcinoma and TME components using edgeR paired analysis. Gene set enrichment analysis was used to identify overrepresented biologic functions for differentially expressed genes between the pairs. Top significantly varied, differentially expressed genes between the pairs were validated using species-specific qPCR. Results: Histologic analysis of the matched sensitive/resistant pairs revealed a marked increase of the stromal component in the resistant tumors evident by H&E and Masson trichrome staining. A total number of 95 and 312 genes were differentially expressed in the epithelial and TME component, respectively (FDR1.5). The stemness marker SOX2 was significantly upregulated in the human component of the resistant tumors as confirmed by qPCR. Genes involved in stroma remodeling were upregulated in the “murine” component of the resistant tumors, and upregulation was confirmed by qPCR (normalized to a murine housekeeper gene). Gene set enrichment analysis showed that the upregulated genes in TME of resistant models were enriched for cytoskeleton related functions, while the downregulated genes were enriched for extracellular matrix-related functions. Conclusion: A bioinformatics strategy able to map human and murine reads separately may identify resistance mechanisms attributed to the two tumor subcompartments and could lead to strategies to potentially target each separately or concurrently. Citation Format: Valentina Zanfagnin, Yuanhang Liu, Pritha Chanana, Xiaonan Pritha, Beverly Long, Joe Pathoulas, Brad Evans, Jaime Davila, Chen Wang, Saravut John Weroha. Transcriptional profiling of tumor stroma using ovarian cancer PDX models with induced platinum-resistance [abstract]. In: Proceedings of the AACR Special Conference on Advances in Ovarian Cancer Research; 2019 Sep 13-16, 2019; Atlanta, GA. Philadelphia (PA): AACR; Clin Cancer Res 2020;26(13_Suppl):Abstract nr B18.

CD4 and CD8 T cells are vital components of the immune system. We found that histone deacetylase 3 (HDAC3) is critical for the development of CD4 T cells, as HDAC3-deficient DP thymocytes generate only CD8SP thymocytes in mice. In the absence of HDAC3, MHC Class II-restricted OT-II thymocytes are redirected to the CD8 cytotoxic lineage, which occurs with accelerated kinetics. Analysis of histone acetylation and RNA-seq reveals that HDAC3-deficient DP thymocytes are biased towards the CD8 lineage prior to positive selection. Commitment to the CD4 or CD8 lineage is determined by whether persistent TCR signaling or cytokine signaling predominates, respectively. Despite elevated IL-21R/γc/STAT5 signaling in HDAC3-deficient DP thymocytes, blocking IL-21R does not restore CD4 lineage commitment. Instead, HDAC3 binds directly to CD8-lineage promoting genes. Thus, HDAC3 is required to restrain CD8-lineage genes in DP thymocytes for the generation of CD4 T cells.

Objective: To determine whether a monogenic basis explains sudden infant death syndrome (SIDS) using an exome-wide focus. Study design: A cohort of 427 unrelated cases of SIDS (257 male; average age = 2.7 ± 1.9 months) underwent whole-exome sequencing. Exome-wide rare variant analyses were carried out with 278 SIDS cases of European ancestry (173 male; average age = 2.7 ± 1.98 months) and 973 ethnic-matched controls based on 6 genetic models. Ingenuity Pathway Analysis also was performed. The cohort was collected in collaboration with coroners, medical examiners, and pathologists by St George's University of London, United Kingdom, and Mayo Clinic, Rochester, Minnesota. Whole-exome sequencing was performed at the Genomic Laboratory, Kings College London, United Kingdom, or Mayo Clinic's Medical Genome Facility, Rochester, Minnesota. Results: Although no exome-wide significant (P < 2.5 × 10−6) difference in burden of ultra-rare variants was detected for any gene, 405 genes had a greater prevalence (P

CD4 and CD8 T cells are vital components of the immune system. We found that histone deacetylase 3 (HDAC3) is critical for the development of CD4 T cells, as HDAC3-deficient DP thymocytes generate only CD8SP thymocytes in mice. In the absence of HDAC3, MHC Class II-restricted OT-II thymocytes are redirected to the CD8 cytotoxic lineage, which occurs with accelerated kinetics. Analysis of histone acetylation and RNA-seq reveals that HDAC3-deficient DP thymocytes are biased towards the CD8 lineage prior to positive selection. Commitment to the CD4 or CD8 lineage is determined by whether persistent TCR signaling or cytokine signaling predominates, respectively. Despite elevated IL-21R/γc/STAT5 signaling in HDAC3-deficient DP thymocytes, blocking IL-21R does not restore CD4 lineage commitment. Instead, HDAC3 binds directly to CD8-lineage promoting genes. Thus, HDAC3 is required to restrain CD8-lineage genes in DP thymocytes for the generation of CD4 T cells.

Sudden infant death syndrome (SIDS) is the commonest cause of sudden death of an infant; however, the genetic basis remains poorly understood. We aimed to identify noncardiac genes underpinning SIDS and determine their prevalence compared with ethnically matched controls.Using exome sequencing we assessed the yield of ultrarare nonsynonymous variants (minor allele frequency [MAF] ≤0.00005, dominant model; MAF ≤0.01, recessive model) in 278 European SIDS cases (62% male; average age =2.7 ± 2 months) versus 973 European controls across 61 noncardiac SIDS-susceptibility genes. The variants were classified according to American College of Medical Genetics and Genomics criteria. Case-control, gene-collapsing analysis was performed in eight candidate biological pathways previously implicated in SIDS pathogenesis.Overall 43/278 SIDS cases harbored an ultrarare single-nucleotide variant compared with 114/973 controls (15.5 vs. 11.7%, p=0.10). Only 2/61 noncardiac genes were significantly overrepresented in cases compared with controls (ECE1, 3/278 [1%] vs. 1/973 [0.1%] p=0.036; SLC6A4, 2/278 [0.7%] vs. 1/973 [0.1%] p=0.049). There was no difference in yield of pathogenic or likely pathogenic variants between cases and controls (1/278 [0.36%] vs. 4/973 [0.41%]; p=1.0). Gene-collapsing analysis did not identify any specific biological pathways to be significantly associated with SIDS.A monogenic basis for SIDS amongst the previously implicated noncardiac genes and their encoded biological pathways is negligible.

Although immune checkpoint inhibitors have resulted in durable clinical benefits in a subset of patients with advanced cancer, some patients who did not respond to initial anti-PD-1 therapy have been found to benefit from the addition of salvage chemotherapy. However, the mechanism responsible for the successful chemoimmunotherapy is not completely understood. Here we show that a subset of circulating CD8+ T cells expressing the chemokine receptor CX3CR1 are able to withstand the toxicity of chemotherapy and are increased in patients with metastatic melanoma who responded to chemoimmunotherapy (paclitaxel and carboplatin plus PD-1 blockade). These CX3CR1+CD8+ T cells have effector memory phenotypes and the ability to efflux chemotherapy drugs via the ABCB1 transporter. In line with clinical observation, our preclinical models identified an optimal sequencing of chemoimmunotherapy that resulted in an increase of CX3CR1+CD8+ T cells. Taken together, we found a subset of PD-1 therapy-responsive CD8+ T cells that were capable of withstanding chemotherapy and executing tumor rejection with their unique abilities of drug efflux (ABCB1), cytolytic activity (granzyme B and perforin), and migration to and retention (CX3CR1 and CD11a) at tumor sites. Future strategies to monitor and increase the frequency of CX3CR1+CD8+ T cells may help to design effective chemoimmunotherapy to overcome cancer resistance to immune checkpoint blockade therapy.

Antibody secreting plasma cells (PCs) play an important role in effective humoral immune responses. The low frequency of bone marrow PCs in humans makes it challenging to obtain sufficient numbers of PCs for biologic studies. Previous studies have employed in vitro model systems to generate cells that morphologically, phenotypically, and functionally resemble normal polyclonal PCs. Gene expression profiles of in vitro generated PCs (IVPCs) mirror their normal counterparts, however to date extensive immunoglobulin (Ig) repertoire analysis of IVPCs is lacking. Here, we used a modified 3-step protocol to generate IVPCs and used RNA-seq to explore the transcriptome with emphasis on the Ig repertoire of plasmablasts and PCs. Total B cells were isolated from 3 normal donors and cultured with various cytokines and the B cell activators CpG ODN and CD40L. RNA was obtained from freshly isolated B cells (Day 0; D0) as well as from Day 4 (D4) plasmablasts, and Day 10 (D10) IVPCs. Morphologically, D10 cells exhibited typical PC morphology, including an eccentric nucleus and perinuclear hof. RNA-seq was performed on total RNA from all 3 donors and time points using the Standard TRuSeq v2 library prep and with paired end sequencing on the Illumina HiSeq 4000 platform. Principle component analysis of gene expression data showed that D0, D4 and D10 cells could be clearly segregated across all 3 normal donors. Of importance, transcripts previously described as distinguishing B cells from PCs were found to be differentially expressed including overexpression of CXCR5, CD19, EBF, CD83, PAX5, IRF8 in D0 B cells and overexpression of IRF4, Blimp-1, XBP1, BCMA, SLAMF7, Syndecan-1, CD38 and CD27 in IVPCs, thus validating our in vitro model for generating PCs. Furthermore, expression of cell cycle related transcripts such as CKS1, CDK1, and CCDN2 followed the pattern of low expression in resting B cells, increased expression in plasmablasts, and decreased expression in IVPCs confirming the cells are actively cycling in a manner comparable to cells in vivo. D10 IVPCs also overexpressed transcripts known to be upregulated during the unfolded protein response. As expected from Ig secreting cells, D10 IVPCs had an over-representation of Ig transcripts. At D0, resting B cells had high levels of IgD and IgM heavy chain (HC) transcripts. At D10, IgM transcripts modestly increased with Log2 fold change (FC) = 3 and as expected, IgD levels decreased significantly (Log2 FC = -2.2). IgA and IgG isotype transcripts significantly increased at D10 (Log2 FC > 6.0) with the IgG4 subtype having the greatest Log2 FC at 8.4. Next we focused on the Ig repertoire of D0, D4, and D10 cells. By aligning to known germline Ig sequences in IMGT/V-Quest (www.imgt.org) and then assembling the paired ends of D0, D4 and D10 Ig transcripts, we were able to analyze the Ig repertoire. Since the Ig HC variable (V) region is encoded by V, diversity (D) and joining (J) segments, only fragments that could be confidently determined were considered. All but 3 IGHV transcripts (IGHV3-35, IGHV3-47 and IGHV7-8) and 2 IGHD transcripts (IGHD4-4 and IGHD5-5) were found and all IGHJ segments were represented across the differentiation spectrum. In D0 cells, the number of unique VDJ combinations ranged from 643 to 863 across all 3 normal samples and increased to a range of 2524 to 2867 in D10 IVPCs. When looking at the differential expression of each VDJ combination from D0 to D10, a pairwise t-test for relative frequency showed that there was no significant change greater than 1%, suggesting the repertoire diversity was not skewed, thus proving the conditions for stimulation were not targeting any one starting B cell. Our data also allowed us to track clonal expansions during differentiation as defined by the increasing frequency of sequences with identical nucleotide sequence in the V region and CDR3 (including D and J regions). Hence, a single sequence could be tracked from D0 to D10. Of interest, in a small sampling of the total available sequences, only those B cells with a mutated IGHV region, characteristic of a memory B cell, went on to expand in this system whereas B cells with an unmutated IGHV did not. Our analysis of the Ig repertoire of IVPCs suggests this system provides a functional model to study Ig repertoire along the B cell differentiation process and further delineate the conditions that may result in a clonal expansion, a hallmark of many hematologic malignancies including multiple myeloma. Disclosures No relevant conflicts of interest to declare.

BackgroundRNA 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.MethodsWe 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.ResultsWe 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.ConclusionsThe 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.