Lipoblastoma-Like Tumor and Fibrosarcoma-Like Lipomatous Neoplasm Represent the Same Entity: A Clinicopathologic and Molecular Genetic Study of 23 Cases Occurring in Both Men and Women at Diverse Locations.

Lipoblastoma-like tumor (LLT) is a benign soft tissue tumor demonstrating mixed morphologic features of lipoblastoma, myxoid liposarcoma, and spindle cell lipoma but lacking genetic alterations associated with those tumors. LLT was originally thought to be specific to the vulva but has since been reported in the paratesticular region. The morphologic features of LLT overlap with those of "fibrosarcoma-like lipomatous neoplasm" (FLLN), a rare, indolent adipocytic neoplasm considered by some to form part of the spectrum of atypical spindle cell and pleomorphic lipomatous tumor. We compared the morphologic, immunohistochemical, and genetic features of 23 tumors previously classified as LLT (n = 17) and FLLN (n = 6). The 23 tumors occurred in 13 women and 10 men (mean age, 42 years; range, 17 to 80 years). Eighteen (78%) cases arose in the inguinogenital region, whereas 5 tumors (22%) involved noninguinogenital soft tissue, including the flank (n = 1), shoulder (n = 1), foot (n = 1), forearm (n = 1), and chest wall (n = 1). Microscopically, the tumors were lobulated and septated, with variably collagenized fibromyxoid stroma, prominent thin-walled vessels, scattered univacuolated or bivacuolated lipoblasts, and a minor component of mature adipose tissue. Using immunohistochemistry, 5 tumors (42%) showed complete RB1 loss, with partial loss in 7 cases (58%). RNA sequencing, chromosomal microarray, and DNA next-generation sequencing study results were negative for significant alterations. There were no clinical, morphologic, immunohistochemical, or molecular genetic differences between cases previously classified as LLT or FLLN. Clinical follow-up (11 patients [48%]; range, 2-276 months; mean, 48.2 months) showed all patients were alive without disease, and only one patient had experienced a single local recurrence. We conclude that LLT and FLLN represent the same entity, for which "LLT" seems most appropriate. LLT may occur in either sex and any superficial soft tissue location. Careful morphologic study and appropriate ancillary testing should allow for the distinction of LLT from its potential mimics.

Multi-Omic Architecture of Obstructive Hypertrophic Cardiomyopathy.

BACKGROUND: Hypertrophic cardiomyopathy (HCM) is characterized by asymmetric left ventricular hypertrophy. Currently, hypertrophy pathways responsible for HCM have not been fully elucidated. Their identification could serve as a nidus for the generation of novel therapeutics aimed at halting disease development or progression. Herein, we performed a comprehensive multi-omic characterization of hypertrophy pathways in HCM. METHODS: Flash-frozen cardiac tissues were collected from genotyped HCM patients (n=97) undergoing surgical myectomy and tissue from 23 controls. RNA sequencing and mass spectrometry-enabled deep proteome and phosphoproteomic assessment were performed. Rigorous differential expression, gene set enrichment, and pathway analyses were performed to characterize HCM-mediated alterations with emphasis on hypertrophy pathways. RESULTS: We identified transcriptional dysregulation with 1246 (8%) differentially expressed genes and elucidated downregulation of 10 hypertrophy pathways. Deep proteomic analysis identified 411 proteins (9%) that differed between HCM and controls with strong dysregulation of metabolic pathways. Seven hypertrophy pathways were upregulated with antagonistic upregulation of 5 of 10 hypertrophy pathways shown to be downregulated in the transcriptome. Most upregulated hypertrophy pathways encompassed the rat sarcoma-mitogen-activated protein kinase signaling cascade. Phosphoproteomic analysis demonstrated hyperphosphorylation of the rat sarcoma-mitogen-activated protein kinase system suggesting activation of this signaling cascade. There was a common transcriptomic and proteomic profile regardless of genotype. CONCLUSIONS: At time of surgical myectomy, the ventricular proteome, independent of genotype, reveals widespread upregulation and activation of hypertrophy pathways, mainly involving the rat sarcoma-mitogen-activated protein kinase signaling cascade. In addition, there is a counterregulatory transcriptional downregulation of the same pathways. Rat sarcoma-mitogen-activated protein kinase activation may serve a crucial role in hypertrophy observed in HCM.

Polymorphous Low-Grade Neuroepithelial Tumor of the Young (PLNTY): Molecular Profiling Confirms Frequent MAPK Pathway Activation.

Polymorphous low-grade neuroepithelial tumor of the young (PLNTY) is a recently described epileptogenic tumor characterized by oligodendroglioma-like components, aberrant CD34 expression, and frequent mitogen-activated protein kinase (MAPK) pathway activation. We molecularly profiled 13 cases with diagnostic histopathological features of PLNTY (10 female; median age, 16 years; range, 5-52). Patients frequently presented with seizures (9 of 12 with available history) and temporal lobe tumors (9 of 13). MAPK pathway activating alterations were identified in all 13 cases. Fusions were present in the 7 youngest patients: FGFR2-CTNNA3 (n = 2), FGFR2-KIAA1598 (FGFR2-SHTN1) (n = 1), FGFR2-INA (n = 1), FGFR2-MPRIP (n = 1), QKI-NTRK2 (n = 1), and KIAA1549-BRAF (n = 1). BRAF V600E mutation was present in 6 patients (17 years or older). Two fusion-positive cases additionally harbored TP53/RB1 abnormalities suggesting biallelic inactivation. Copy number changes predominantly involving whole chromosomes were observed in all 10 evaluated cases, with losses of chromosome 10q occurring with FGFR2-KIAA1598 (SHTN1)/CTNNA3 fusions. The KIAA1549-BRAF and QKI-NTRK2 fusions were associated respectively with a 7q34 deletion and 9q21 duplication. This study shows that despite its name, PLNTY also occurs in older adults, who frequently show BRAF V600E mutation. It also expands the spectrum of the MAPK pathway activating alterations associated with PLNTY and demonstrates recurrent chromosomal copy number changes consistent with chromosomal instability.

EWSR1-WT1 gene fusions in neoplasms other than desmoplastic small round cell tumor: a report of three unusual tumors involving the female genital tract and review of the literature.

Desmoplastic small round cell tumor (DSRCT) is a high-grade round cell sarcoma that typically arises in the abdominopelvic cavity of young males, co-expresses keratins and desmin, and carries a pathognomonic EWSR1-WT1 gene fusion. The EWSR1-WT1 gene fusion is generally considered specific for DSRCT, although there are two reports of this fusion in tumors otherwise lacking features of DSRCT. We report three female genital tract tumors with EWSR1-WT1 fusions but showing morphologic and immunohistochemical features incompatible with DSRCT. The tumors occurred in the uterine cervix, uterine corpus/ovaries, and vagina, respectively, of 46, 30, and 20-year-old women. Two tumors consisted of a sheet-like to fascicular proliferation of relatively uniform spindled to occasionally more epithelioid cells arrayed about thick-walled, hyalinized, and capillary-sized vessels, with distinctive areas of pseudovascular change, and absence of desmoplastic stroma. The third tumor resembled a monomorphic spindle cell sarcoma with necrosis. All had diffuse desmin and variable but more limited keratin expression, two of three expressed smooth muscle actin, and all were negative for h-caldesmon, CD10, estrogen receptor, myogenin, N-terminus WT-1, and S100 protein. One patient received neoadjuvant chemotherapy and radiation therapy followed by resection and is disease-free 42 months after diagnosis. Another patient was managed by resection only and is disease-free 9 months after initial diagnosis. The remaining patient recently underwent resection of multifocal pelvic disease. Comprehensive differential gene expression analysis on two tumors compared to two classic DSRCTs with known EWSR1-WT1 fusions resulted in 1726 genes that were differentially expressed (log2 fold change >2 or < -2) and statistically significant (FDR < 5%). In combination with previous reports, our findings suggest pleiotropy of the EWSR1-WT1 fusion is possible and not limited to DSRCT. Subsets of non-DSRCT EWSR1-WT1 positive tumors may represent discrete entities, but further study is necessary.

RNA-Seq Reveals Differences in Expressed Tumor Mutation Burden in Colorectal and Endometrial Cancers with and without Defective DNA-Mismatch Repair.

Tumor mutation burden (TMB) is an emerging biomarker of immunotherapy response. RNA sequencing in FFPE tissue samples was used for determining TMB in microsatellite-stable (MSS) and microsatellite instability-high (MSI-H) tumors in patients with colorectal or endometrial cancer. Tissue from tumors and paired normal tissue from 46 MSI-H and 12 MSS cases were included. Of the MSI-H tumors, 29 had defective DNA mismatch-repair mutations, and 17 had MLH1 promoter hypermethylation. TMB was measured using the expressed somatic nucleotide variants (eTMB). A method of accurate measurement of eTMB was developed that removes FFPE-derived artifacts by leveraging mutation signatures. There was a significant difference in the median eTMB values observed between MSI-H and MSS cases: 27.3 versus 6.7 mutations/megabase (mut/Mb) (P = 3.5 × 10-9). Among tumors with defective DNA-mismatch repair, those with mismatch-repair mutations had a significantly higher median eTMB than those with hypermethylation: 28.1 versus 17.5 mut/Mb (P = 0.037). Multivariate analysis showed that MSI status, tumor type (endometrial or colorectal), and age were significantly associated with eTMB. Additionally, using whole-exome sequencing in a subset of these patients, it was determined that DNA TMB correlated well with eTMB (Spearman correlation coefficient, 0.83). These results demonstrate that RNA sequencing can be used for measuring eTMB in FFPE tumor specimens.

Frequency of MicroRNA Response Elements Identifies Pathologically Relevant Signaling Pathways in Triple-Negative Breast Cancer.

Complex interactions between mRNAs and microRNAs influence cellular functions. The mRNA-microRNA interactions also determine the post-transcriptional availability of mRNAs and unbound microRNAs. MicroRNAs binds to one or more microRNA response elements (MREs) located on the 3'UTR of mRNAs. In this study, we leveraged MREs and their frequencies in cancer and matched normal tissues to obtain insights into disease-specific interactions between mRNAs and microRNAs. We developed a bioinformatics method "ReMIx" that utilizes RNA sequencing (RNA-Seq) data to quantify MRE frequencies across the transcriptome. We applied ReMIx to triple-negative (TN) breast cancer tumor-normal adjacent pairs and identified MREs specific to TN tumors. ReMIx identified candidate mRNAs and microRNAs in the MAPK signaling cascade. Further analysis of MAPK gene regulatory networks revealed microRNA partners that influence and modulate MAPK signaling. In conclusion, we demonstrate a novel method of using MREs in the identification of functionally relevant mRNA-microRNA interactions in TN breast cancer.

Desmoplastic Infantile Ganglioglioma: A MAPK Pathway-Driven and Microglia/Macrophage-Rich Neuroepithelial Tumor.

MAPK pathway activation has been recurrently observed in desmoplastic infantile ganglioglioma/astrocytoma (DIG/DIA) with reported disproportionally low mutation allele frequencies relative to the apparent high tumor content, suggesting that MAPK pathway alterations may be subclonal. We sought to expand the number of molecularly profiled cases and investigate if tumor cell composition could account for the observed low mutation allele frequencies. Molecular (targeted neuro-oncology next-generation sequencing/RNA sequencing and OncoScan microarray) and immunohistochemical (CD68-PGM1/CD163/CD14/CD11c/lysozyme/CD3/CD20/CD34/PD-L1) studies were performed in 7 DIG. Activating MAPK pathway alterations were identified in 4 (57%) cases: 3 had a BRAF mutation (V600E/V600D/V600_W604delinsDQTDG, at 8%-27% variant allele frequency) and 1 showed a TPM3-NRTK1 fusion. Copy number changes were infrequent and nonrecurrent. All tumors had at least 30% of cells morphologically and immunophenotypically consistent with microglial/macrophage lineage. Two subtotally resected tumors regrew; 1 was re-excised and received adjuvant treatment (chemotherapy/targeted therapy), with clinical response to targeted therapy only. Even with residual tumor, all patients are alive (median follow-up, 83 months; 19-139). This study further supports DIG as another MAPK pathway-driven neuroepithelial tumor, thus expanding potential treatment options for tumors not amenable to surgical cure, and suggests that DIG is a microglia/macrophage-rich neuroepithelial tumor with frequent low driver mutation allele frequencies.

An expanded variant list and assembly annotation identifies multiple novel coding and noncoding genes for prostate cancer risk using a normal prostate tissue eQTL data set.

Prostate cancer (PrCa) is highly heritable; 284 variants have been identified to date that are associated with increased prostate cancer risk, yet few genes contributing to its development are known. Expression quantitative trait loci (eQTL) studies link variants with affected genes, helping to determine how these variants might regulate gene expression and may influence prostate cancer risk. In the current study, we performed eQTL analysis on 471 normal prostate epithelium samples and 249 PrCa-risk variants in 196 risk loci, utilizing RNA sequencing transcriptome data based on ENSEMBL gene definition and genome-wide variant data. We identified a total of 213 genes associated with known PrCa-risk variants, including 141 protein-coding genes, 16 lncRNAs, and 56 other non-coding RNA species with differential expression. Compared to our previous analysis, where RefSeq was used for gene annotation, we identified an additional 130 expressed genes associated with known PrCa-risk variants. We detected an eQTL signal for more than half (n = 102, 52%) of the 196 loci tested; 52 (51%) of which were a Group 1 signal, indicating high linkage disequilibrium (LD) between the peak eQTL variant and the PrCa-risk variant (r2>0.5) and may help explain how risk variants influence the development of prostate cancer.

RNA sequencing identifies a novel USP9X-USP6 promoter swap gene fusion in a primary aneurysmal bone cyst.

Primary aneurysmal bone cyst (ABC) is a benign multiloculated cystic lesion of bone that is defined cytogenetically by USP6 gene rearrangements. Rearrangements involving USP6 are promoter swaps, usually generated by fusion of the noncoding upstream exons of different partner genes with exon 1 or 2 of USP6, thus leading to transcriptional upregulation of full-length USP6 coding sequence. Testing for USP6 rearrangements is used diagnostically to distinguish it from secondary ABC and other giant cell-rich primary bone tumors. In this report, we present a case of a 16-year-old male with a primary ABC of the left distal femur. USP6 break apart fluorescence in situ hybridization was positive for a rearrangement and conventional chromosome analysis identified a reciprocal X;17 translocation. In order to identify the putative USP6 fusion partner, we performed RNA sequencing and uncovered a novel USP9X-USP6 promoter swap fusion. This result was confirmed by reverse transcription-polymerase chain reaction (RT-PCR) and by mate pair sequencing thus showing the utility of these alternative methodologies in identifying novel fusion candidates. Ubiquitin-specific protease 9X (USP9X), like USP6, encodes a highly conserved substrate-specific deubiquitylating enzyme. USP9X is highly expressed in a number of tissue types and acts as both an oncogene and tumor suppressor in several human cancers. We conclude that oncogenic activation of USP6 via USP9X promoter exchange represents a novel driver of primary ABC formation.

Characterization of three ciliopathy pedigrees expands the phenotype associated with biallelic C2CD3 variants.

Whole exome sequencing (WES) is utilized in diagnostic odyssey cases to identify the underlying genetic cause associated with complex phenotypes. Recent publications suggest that WES reveals the genetic cause in ~25% of these cases and is most successful when applied to children with neurological disease. The residual 75% of cases remain genetically elusive until more information becomes available in the literature or functional studies are pursued. WES performed on three families with presumed ciliopathy diagnoses, including orofaciodigital (OFD) syndrome, fetal encephalocele, or Joubert-related disorder, identified compound heterozygous variants in C2CD3. Biallelic variants in C2CD3 have previously been associated with ciliopathies, including OFD syndrome type 14 (OFD14; MIM: 615948). As three of the six identified variants were predicted to affect splicing, exon-skipping analysis using either RNA sequencing or PCR-based methods were completed to determine the pathogenicity of these variants, and showed that each of the splicing variants led to a frameshifted protein product. Using these studies in combination with the 2015 ACMG guidelines, each of the six identified variants were classified as either pathogenic or likely pathogenic, and are therefore likely responsible for our patients' phenotypes. Each of the families had a distinct clinical phenotype and severity of disease, extending from lethal to viable. These findings highlight that there is a broad phenotypic spectrum associated with C2CD3-mediated disease and not all patients present with the typical features of OFD14.

Development and Verification of an RNA Sequencing (RNA-Seq) Assay for the Detection of Gene Fusions in Tumors.

We assessed the performance characteristics of an RNA sequencing (RNA-Seq) assay designed to detect gene fusions in 571 genes to help manage patients with cancer. Polyadenylated RNA was converted to cDNA, which was then used to prepare next-generation sequencing libraries that were sequenced on an Illumina HiSeq 2500 instrument and analyzed with an in-house developed bioinformatic pipeline. The assay identified 38 of 41 gene fusions detected by another method, such as fluorescence in situ hybridization or RT-PCR, for a sensitivity of 93%. No false-positive gene fusions were identified in 15 normal tissue specimens and 10 tumor specimens that were negative for fusions by RNA sequencing or Mate Pair NGS (100% specificity). The assay also identified 22 fusions in 17 tumor specimens that had not been detected by other methods. Eighteen of the 22 fusions had not previously been described. Good intra-assay and interassay reproducibility was observed with complete concordance for the presence or absence of gene fusions in replicates. The analytical sensitivity of the assay was tested by diluting RNA isolated from gene fusion-positive cases with fusion-negative RNA. Gene fusions were generally detectable down to 12.5% dilutions for most fusions and as little as 3% for some fusions. This assay can help identify fusions in patients with cancer; these patients may in turn benefit from both US Food and Drug Administration-approved and investigational targeted therapies.

Validation of Osteogenic Properties of Cytochalasin D by High-Resolution RNA-Sequencing in Mesenchymal Stem Cells Derived from Bone Marrow and Adipose Tissues.

Differentiation of mesenchymal stromal/stem cells (MSCs) involves a series of molecular signals and gene transcription events required for attaining cell lineage commitment. Modulation of the actin cytoskeleton using cytochalasin D (CytoD) drives osteogenesis at early timepoints in bone marrow-derived MSCs and also initiates a robust osteogenic differentiation program in adipose tissue-derived MSCs. To understand the molecular basis for these pronounced effects on osteogenic differentiation, we investigated global changes in gene expression in CytoD-treated murine and human MSCs by high-resolution RNA-sequencing (RNA-seq) analysis. A three-way bioinformatic comparison between human adipose tissue-derived MSCs (hAMSCs), human bone marrow-derived MSCs (hBMSCs), and mouse bone marrow-derived MSCs (mBMSCs) revealed significant upregulation of genes linked to extracellular matrix organization, cell adhesion and bone metabolism. As anticipated, the activation of these differentiation-related genes is accompanied by a downregulation of nuclear and cell cycle-related genes presumably reflecting cytostatic effects of CytoD. We also identified eight novel CytoD activated genes-VGLL4, ARHGAP24, KLHL24, RCBTB2, BDH2, SCARF2, ACAD10, HEPH-which are commonly upregulated across the two species and tissue sources of our MSC samples. We selected the Hippo pathway-related VGLL4 gene, which encodes the transcriptional co-factor Vestigial-like 4, for further study because this pathway is linked to osteogenesis. VGLL4 small interfering RNA depletion reduces mineralization of hAMSCs during CytoD-induced osteogenic differentiation. Together, our RNA-seq analyses suggest that while the stimulatory effects of CytoD on osteogenesis are pleiotropic and depend on the biological state of the cell type, a small group of genes including VGLL4 may contribute to MSC commitment toward the bone lineage.

Network-directed cis-mediator analysis of normal prostate tissue expression profiles reveals downstream regulatory associations of prostate cancer susceptibility loci.

Large-scale genome-wide association studies have identified multiple single-nucleotide polymorphisms associated with risk of prostate cancer. Many of these genetic variants are presumed to be regulatory in nature; however, follow-up expression quantitative trait loci (eQTL) association studies have to-date been restricted largely to cis-acting associations due to study limitations. While trans-eQTL scans suffer from high testing dimensionality, recent evidence indicates most trans-eQTL associations are mediated by cis-regulated genes, such as transcription factors. Leveraging a data-driven gene co-expression network, we conducted a comprehensive cis-mediator analysis using RNA-Seq data from 471 normal prostate tissue samples to identify downstream regulatory associations of previously identified prostate cancer risk variants. We discovered multiple trans-eQTL associations that were significantly mediated by cis-regulated transcripts, four of which involved risk locus 17q12, proximal transcription factor HNF1B, and target trans-genes with known HNF response elements (MIA2, SRC, SEMA6A, KIF12). We additionally identified evidence of cis-acting down-regulation of MSMB via rs10993994 corresponding to reduced co-expression of NDRG1. The majority of these cis-mediator relationships demonstrated trans-eQTL replicability in 87 prostate tissue samples from the Gene-Tissue Expression Project. These findings provide further biological context to known risk loci and outline new hypotheses for investigation into the etiology of prostate cancer.

Gastroblastoma harbors a recurrent somatic MALAT1-GLI1 fusion gene.

Gastroblastoma is a rare distinctive biphasic tumor of the stomach. The molecular biology of gastroblastoma has not been studied, and no affirmative diagnostic markers have been developed. We retrieved two gastroblastomas from the consultation practices of the authors and performed transcriptome sequencing on formalin-fixed paraffin-embedded tissue. Recurrent predicted fusion genes were validated at genomic and RNA levels. The presence of the fusion gene was confirmed on two additional paraffin-embedded cases of gastroblastoma. Control cases of histologic mimics (biphasic synovial sarcoma, leiomyoma, leiomyosarcoma, desmoid-type fibromatosis, EWSR1-FLI1-positive Ewing sarcoma, Wilms' tumor, gastrointestinal stromal tumor, plexiform fibromyxoma, Sonic hedgehog-type medulloblastomas, and normal gastric mucosa and muscularis propria were also analyzed. The gastroblastomas affected two males and two females aged 9-56 years. Transcriptome sequencing identified recurrent somatic MALAT1-GLI1 fusion genes, which were predicted to retain the key domains of GLI1. The MALAT1-GLI1 fusion gene was validated by break-apart and dual-fusion FISH and RT-PCR. The additional two gastroblastomas were also positive for the MALAT1-GLI1 fusion gene. None of the other control cases harbored MALAT1-GLI1. Overexpression of GLI1 in the cases of gastroblastomas was confirmed at RNA and protein levels. Pathway analysis revealed activation of the Sonic hedgehog pathway in gastroblastoma and gene expression profiling showed that gastroblastomas grouped together and were most similar to Sonic hedgehog-type medulloblastomas. In summary, we have identified an oncogenic MALAT1-GLI1 fusion gene in all cases of gastroblastoma that may serve as a diagnostic biomarker. The fusion gene is predicted to encode a protein that includes the zinc finger domains of GLI1 and results in overexpression of GLI1 protein and activation of the Sonic hedgehog pathway.

The Role of the Histone Methyltransferase Enhancer of Zeste Homolog 2 (EZH2) in the Pathobiological Mechanisms Underlying Inflammatory Bowel Disease (IBD).

Regulatory T (Treg) cells expressing the transcription factor FOXP3 play a pivotal role in maintaining immunologic self-tolerance. We and others have shown previously that EZH2 is recruited to the FOXP3 promoter and its targets in Treg cells. To further address the role for EZH2 in Treg cellular function, we have now generated mice that lack EZH2 specifically in Treg cells (EZH2Δ/ΔFOXP3+). We find that EZH2 deficiency in FOXP3+ T cells results in lethal multiorgan autoimmunity. We further demonstrate that EZH2Δ/ΔFOXP3+ T cells lack a regulatory phenotype in vitro and secrete proinflammatory cytokines. Of special interest, EZH2Δ/ΔFOXP3+ mice develop spontaneous inflammatory bowel disease. Guided by these results, we assessed the FOXP3 and EZH2 gene networks by RNA sequencing in isolated intestinal CD4+ T cells from patients with Crohn's disease. Gene network analysis demonstrates that these CD4+ T cells display a Th1/Th17-like phenotype with an enrichment of gene targets shared by FOXP3 and EZH2. Combined, these results suggest that the inflammatory milieu found in Crohn's disease could lead to or result from deregulation of FOXP3/EZH2-enforced T cell gene networks contributing to the underlying intestinal inflammation.

Circular RNAs and their associations with breast cancer subtypes.

Circular RNAs (circRNAs) are highly stable forms of non-coding RNAs with diverse biological functions. They are implicated in modulation of gene expression thus affecting various cellular and disease processes. Based on existing bioinformatics approaches, we developed a comprehensive workflow called Circ-Seq to identify and report expressed circRNAs. Circ-Seq also provides informative genomic annotation along circRNA fused junctions thus allowing prioritization of circRNA candidates. We applied Circ-Seq first to RNA-sequence data from breast cancer cell lines and validated one of the large circRNAs identified. Circ-Seq was then applied to a larger cohort of breast cancer samples (n = 885) provided by The Cancer Genome Atlas (TCGA), including tumors and normal-adjacent tissue samples. Notably, circRNA results reveal that normal-adjacent tissues in estrogen receptor positive (ER+) subtype have relatively higher numbers of circRNAs than tumor samples in TCGA. Similar phenomenon of high circRNA numbers were observed in normal breast-mammary tissues from the Genotype-Tissue Expression (GTEx) project. Finally, we observed that number of circRNAs in normal-adjacent samples of ER+ subtype is inversely correlated to the risk-of-relapse proliferation (ROR-P) score for proliferating genes, suggesting that circRNA frequency may be a marker for cell proliferation in breast cancer. The Circ-Seq workflow will function for both single and multi-threaded compute environments. We believe that Circ-Seq will be a valuable tool to identify circRNAs useful in the diagnosis and treatment of other cancers and complex diseases.

Impact of RNA degradation on fusion detection by RNA-seq.

BACKGROUND: RNA-seq is a well-established method for studying the transcriptome. Popular methods for library preparation in RNA-seq such as Illumina TruSeq® RNA v2 kit use a poly-A pulldown strategy. Such methods can cause loss of coverage at the 5' end of genes, impacting the ability to detect fusions when used on degraded samples. The goal of this study was to quantify the effects RNA degradation has on fusion detection when using poly-A selected mRNA and to identify the variables involved in this process. RESULTS: Using both artificially and naturally degraded samples, we found that there is a reduced ability to detect fusions as the distance of the breakpoint from the 3' end of the gene increases. The median transcript coverage decreases exponentially as a function of the distance from the 3' end and there is a linear relationship between the coverage decay rate and the RNA integrity number (RIN). Based on these findings we developed plots that show the probability of detecting a gene fusion ("sensitivity") as a function of the distance of the fusion breakpoint from the 3' end. CONCLUSIONS: This study developed a strategy to assess the impact that RNA degradation has on the ability to detect gene fusions by RNA-seq.