Characterization of cell-type specific circular RNAs associated with colorectal cancer metastasis.

Colorectal cancer (CRC) is the most common gastrointestinal malignancy and a leading cause of cancer deaths in the United States. More than half of CRC patients develop metastatic disease (mCRC) with an average 5-year survival rate of 13%. Circular RNAs (circRNAs) have recently emerged as important tumorigenesis regulators; however, their role in mCRC progression remains poorly characterized. Further, little is known about their cell-type specificity to elucidate their functions in the tumor microenvironment (TME). To address this, we performed total RNA sequencing (RNA-seq) on 30 matched normal, primary and metastatic samples from 14 mCRC patients. Additionally, five CRC cell lines were sequenced to construct a circRNA catalog in CRC. We detected 47 869 circRNAs, with 51% previously unannotated in CRC and 14% novel candidates when compared to existing circRNA databases. We identified 362 circRNAs differentially expressed in primary and/or metastatic tissues, termed circular RNAs associated with metastasis (CRAMS). We performed cell-type deconvolution using published single-cell RNA-seq datasets and applied a non-negative least squares statistical model to estimate cell-type specific circRNA expression. This predicted 667 circRNAs as exclusively expressed in a single cell type. Collectively, this serves as a valuable resource, TMECircDB (accessible at https://www.maherlab.com/tmecircdb-overview), for functional characterization of circRNAs in mCRC, specifically in the TME.

Integrated analysis of genomic and transcriptomic data for the discovery of splice-associated variants in cancer.

Somatic mutations within non-coding regions and even exons may have unidentified regulatory consequences that are often overlooked in analysis workflows. Here we present RegTools ( www.regtools.org ), a computationally efficient, free, and open-source software package designed to integrate somatic variants from genomic data with splice junctions from bulk or single cell transcriptomic data to identify variants that may cause aberrant splicing. We apply RegTools to over 9000 tumor samples with both tumor DNA and RNA sequence data. RegTools discovers 235,778 events where a splice-associated variant significantly increases the splicing of a particular junction, across 158,200 unique variants and 131,212 unique junctions. To characterize these somatic variants and their associated splice isoforms, we annotate them with the Variant Effect Predictor, SpliceAI, and Genotype-Tissue Expression junction counts and compare our results to other tools that integrate genomic and transcriptomic data. While many events are corroborated by the aforementioned tools, the flexibility of RegTools also allows us to identify splice-associated variants in known cancer drivers, such as TP53, CDKN2A, and B2M, and other genes.

Pan-Cancer Analysis Reveals Recurrent BCAR4 Gene Fusions across Solid Tumors.

Chromosomal rearrangements often result in active regulatory regions juxtaposed upstream of an oncogene to generate an expressed gene fusion. Repeated activation of a common downstream partner-with differing upstream regions across a patient cohort-suggests a conserved oncogenic role. Analysis of 9,638 patients across 32 solid tumor types revealed an annotated long noncoding RNA (lncRNA), Breast Cancer Anti-Estrogen Resistance 4 (BCAR4), was the most prevalent, uncharacterized, downstream gene fusion partner occurring in 11 cancers. Its oncogenic role was confirmed using multiple cell lines with endogenous BCAR4 gene fusions. Furthermore, overexpressing clinically prevalent BCAR4 gene fusions in untransformed cell lines was sufficient to induce an oncogenic phenotype. We show that the minimum common region to all gene fusions harbors an open reading frame that is necessary to drive proliferation. IMPLICATIONS: BCAR4 gene fusions represent an underappreciated class of gene fusions that may have biological and clinical implications across solid tumors.

LINC00355 regulates p27KIP expression by binding to MENIN to induce proliferation in late-stage relapse breast cancer.

Late-stage relapse (LSR) in patients with breast cancer (BC) occurs more than five years and up to 10 years after initial treatment and has less than 30% 5-year relative survival rate. Long non-coding RNAs (lncRNAs) play important roles in BC yet have not been studied in LSR BC. Here, we identify 1127 lncRNAs differentially expressed in LSR BC via transcriptome sequencing and analysis of 72 early-stage and 24 LSR BC patient tumors. Decreasing expression of the most up-regulated lncRNA, LINC00355, in BC and MCF7 long-term estrogen deprived cell lines decreases cellular invasion and proliferation. Subsequent mechanistic studies show that LINC00355 binds to MENIN and changes occupancy at the CDKN1B promoter to decrease p27Kip. In summary, this is a key study discovering lncRNAs in LSR BC and LINC00355 association with epigenetic regulation and proliferation in BC.

DANSR: A Tool for the Detection of Annotated and Novel Small RNAs.

Existing small noncoding RNA analysis tools are optimized for processing short sequencing reads (17-35 nucleotides) to monitor microRNA expression. However, these strategies under-represent many biologically relevant classes of small noncoding RNAs in the 36-200 nucleotides length range (tRNAs, snoRNAs, etc.). To address this, we developed DANSR, a tool for the detection of annotated and novel small RNAs using sequencing reads with variable lengths (ranging from 17-200 nt). While DANSR is broadly applicable to any small RNA dataset, we applied it to a cohort of matched normal, primary, and distant metastatic colorectal cancer specimens to demonstrate its ability to quantify annotated small RNAs, discover novel genes, and calculate differential expression. DANSR is available as an open source tool.

Long non-coding RNA RAMS11 promotes metastatic colorectal cancer progression.

Colorectal cancer (CRC) is the most common gastrointestinal malignancy in the U.S.A. and approximately 50% of patients develop metastatic disease (mCRC). Despite our understanding of long non-coding RNAs (lncRNAs) in primary colon cancer, their role in mCRC and treatment resistance remains poorly characterized. Therefore, through transcriptome sequencing of normal, primary, and distant mCRC tissues we find 148 differentially expressed RNAs Associated with Metastasis (RAMS). We prioritize RAMS11 due to its association with poor disease-free survival and promotion of aggressive phenotypes in vitro and in vivo. A FDA-approved drug high-throughput viability assay shows that elevated RAMS11 expression increases resistance to topoisomerase inhibitors. Subsequent experiments demonstrate RAMS11-dependent recruitment of Chromobox protein 4 (CBX4) transcriptionally activates Topoisomerase II alpha (TOP2α). Overall, recent clinical trials using topoisomerase inhibitors coupled with our findings of RAMS11-dependent regulation of TOP2α supports the potential use of RAMS11 as a biomarker and therapeutic target for mCRC.

Long non-coding RNA LCAL62 / LINC00261 is associated with lung adenocarcinoma prognosis.

BACKGROUND: More than half of non-small cell lung cancer (NSCLC) patients present with metastatic disease at initial diagnosis with an estimated five-year survival rate of ~5%. Despite advances in understanding primary lung cancer oncogenesis metastatic disease remains poorly characterized. Recent studies demonstrate important roles of long non-coding RNAs (lncRNAs) in tumor physiology and as prognostic markers. Therefore, we present the first transcriptome analysis to identify lncRNAs altered in metastatic lung adenocarcinoma leading to the discovery and characterization of the lncRNA LCAL62 as a prognostic biomarker. PATIENTS AND METHODS: RNA-Seq, microarray, nanoString expression, and clinical data from 1,116 LUAD patients across six independent cohorts and 83 LUAD cell lines were used to discover and evaluate the survival association of metastasis associated lncRNAs. Coexpression and gene set enrichment analyses were used to establish gene regulatory networks and implicate metastasis associated lncRNAs in specific biological processes. RESULTS: Our integrative analysis discovered LCAL62 as the most down-regulated lncRNA in metastasis. Further low LCAL62 expression promoted aggressive phenotypes and regulated genes associated with metastasis (such as metastasis repressor FOXA2). Low LCAL62 expression corresponded to poor overall patient survival across five independent lung adenocarcinoma cohorts (n = 881) including our own nanoString validation cohort. CONCLUSION: We discovered that LCAL62 was down-regulated in lung cancer progression to promote invasive phenotypes, and lower expression was significantly associated with poor patient outcome and aggressive lung adenocarcinoma.

Multi-institutional Analysis Shows that Low PCAT-14 Expression Associates with Poor Outcomes in Prostate Cancer.

BACKGROUND: Long noncoding RNAs (lncRNAs) are an emerging class of relatively underexplored oncogenic molecules with biological and clinical significance. Current inadequacies for stratifying patients with aggressive disease presents a strong rationale to systematically identify lncRNAs as clinical predictors in localized prostate cancer. OBJECTIVE: To identify RNA biomarkers associated with aggressive prostate cancer. DESIGN, SETTING, AND PARTICIPANTS: Radical prostatectomy microarray and clinical data was obtained from 910 patients in three published institutional cohorts: Mayo Clinic I (N=545, median follow-up 13.8 yr), Mayo Clinic II (N=235, median follow-up 6.7 yr), and Thomas Jefferson University (N=130, median follow-up 9.6 yr). OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: The primary clinical endpoint was distant metastasis-free survival. Secondary endpoints include prostate cancer-specific survival and overall survival. Univariate and multivariate Cox regression were used to evaluate the association of lncRNA expression and these endpoints. RESULTS AND LIMITATIONS: An integrative analysis revealed Prostate Cancer Associated Transcript-14 (PCAT-14) as the most prevalent lncRNA that is aberrantly expressed in prostate cancer patients. Down-regulation of PCAT-14 expression significantly associated with Gleason score and a greater probability of metastatic progression, overall survival, and prostate cancer-specific mortality across multiple independent datasets and ethnicities. Low PCAT-14 expression was implicated with genes involved in biological processes promoting aggressive disease. In-vitro analysis confirmed that low PCAT-14 expression increased migration while overexpressing PCAT-14 reduced cellular growth, migration, and invasion. CONCLUSIONS: We discovered that androgen-regulated PCAT-14 is overexpressed in prostate cancer, suppresses invasive phenotypes, and lower expression is significantly prognostic for multiple clinical endpoints supporting its significance for predicting metastatic disease that could be used to improve patient management. PATIENT SUMMARY: We discovered that aberrant prostate cancer associated transcript-14 expression during prostate cancer progression is prevalent across cancer patients. Prostate cancer associated transcript-14 is also prognostic for metastatic disease and survival highlighting its importance for stratifying patients that could benefit from treatment intensification.