Oncogenic Role of THOR, a Conserved Cancer/Testis Long Non-coding RNA.

Large-scale transcriptome sequencing efforts have vastly expanded the catalog of long non-coding RNAs (lncRNAs) with varying evolutionary conservation, lineage expression, and cancer specificity. Here, we functionally characterize a novel ultraconserved lncRNA, THOR (ENSG00000226856), which exhibits expression exclusively in testis and a broad range of human cancers. THOR knockdown and overexpression in multiple cell lines and animal models alters cell or tumor growth supporting an oncogenic role. We discovered a conserved interaction of THOR with IGF2BP1 and show that THOR contributes to the mRNA stabilization activities of IGF2BP1. Notably, transgenic THOR knockout produced fertilization defects in zebrafish and also conferred a resistance to melanoma onset. Likewise, ectopic expression of human THOR in zebrafish accelerated the onset of melanoma. THOR represents a novel class of functionally important cancer/testis lncRNAs whose structure and function have undergone positive evolutionary selection.

Expressed pseudogenes in the transcriptional landscape of human cancers.

Pseudogene transcripts can provide a novel tier of gene regulation through generation of endogenous siRNAs or miRNA-binding sites. Characterization of pseudogene expression, however, has remained confined to anecdotal observations due to analytical challenges posed by the extremely close sequence similarity with their counterpart coding genes. Here, we describe a systematic analysis of pseudogene "transcription" from an RNA-Seq resource of 293 samples, representing 13 cancer and normal tissue types, and observe a surprisingly prevalent, genome-wide expression of pseudogenes that could be categorized as ubiquitously expressed or lineage and/or cancer specific. Further, we explore disease subtype specificity and functions of selected expressed pseudogenes. Taken together, we provide evidence that transcribed pseudogenes are a significant contributor to the transcriptional landscape of cells and are positioned to play significant roles in cellular differentiation and cancer progression, especially in light of the recently described ceRNA networks. Our work provides a transcriptome resource that enables high-throughput analyses of pseudogene expression.

Comprehensive Assessment of the Intrinsic Pancreatic Microbiome.

OBJECTIVE: We sought to comprehensively profile tissue and cyst fluid in patients with benign, precancerous, and cancerous conditions of the pancreas to characterize the intrinsic pancreatic microbiome. SUMMARY BACKGROUND DATA: Small studies in pancreatic ductal adenocarcinoma (PDAC) and intraductal papillary mucinous neoplasm (IPMN) have suggested that intra-pancreatic microbial dysbiosis may drive malignant transformation. METHODS: Pancreatic samples were collected at the time of resection from 109 patients. Samples included tumor tissue (control, n=20; IPMN, n=20; PDAC, n=19) and pancreatic cyst fluid (IPMN, n=30; SCA, n=10; MCN, n=10). Assessment of bacterial DNA by quantitative PCR and 16S ribosomal RNA gene sequencing was performed. Downstream analyses determined the relative abundances of individual taxa between groups and compared intergroup diversity. Whole-genome sequencing data from 140 patients with PDAC in the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium (CPTAC) were analyzed to validate findings. RESULTS: Sequencing of pancreatic tissue yielded few microbial reads regardless of diagnosis, and analysis of pancreatic tissue showed no difference in the abundance and composition of bacterial taxa between normal pancreas, IPMN, or PDAC groups. Low-grade dysplasia (LGD) and high-grade dysplasia (HGD) IPMN were characterized by low bacterial abundances with no difference in tissue composition and a slight increase in Pseudomonas and Sediminibacterium in HGD cyst fluid. Decontamination analysis using the CPTAC database confirmed a low-biomass, low-diversity intrinsic pancreatic microbiome that did not differ by pathology. CONCLUSIONS: Our analysis of the pancreatic microbiome demonstrated very low intrinsic biomass that is relatively conserved across diverse neoplastic conditions and thus unlikely to drive malignant transformation.

Timing of Adjuvant Immunotherapy in Stage III Melanoma, Does it Matter?

BACKGROUND: The optimal time to initiate adjuvant immune checkpoint inhibitors (ICI) following resection remains undefined. Herein, we investigated the impact of time to adjuvant ICI on survival in patients with stage III melanoma. METHODS: Patients with resected stage III melanoma receiving adjuvant immune therapy were identified within a multi-institutional retrospective cohort. Patients were stratified by time to adjuvant ICI: within 6 weeks, 6-12 weeks, and greater than 12 weeks from surgery. Recurrence-free survival (RFS) was compared among time strata with Kaplan-Meier and Cox proportional hazards methods in the multi-institutional cohort. RESULTS: Altogether, 626 patients were identified within the multi-institutional cohort: 39% of patients initiated adjuvant ICI within 6 weeks, 42.2% within 6-12 weeks, and 18.8% greater than 12 weeks from surgery. In a multivariate Cox model, adjusting for histology, nodal tumor burden, and pathologic stage, we found that increased time to adjuvant ICI was associated with improved RFS. Patients who initiated adjuvant ICI within 6 weeks of surgery had worse RFS. These findings were preserved in a conditional landmark analysis and separate subgroups of patients with (1) new melanoma diagnoses, (2) occult stage III disease, and (3) those receiving anti-PD-1 monotherapy. CONCLUSIONS: Outcomes for patients with stage III melanoma are not compromised when adjuvant ICI is initiated beyond 6 weeks from resection. Additional work is needed to better understand the underlying mechanisms and implications of timing of adjuvant ICI on long-term outcomes.

TACO produces robust multisample transcriptome assemblies from RNA-seq.

Accurate transcript structure and abundance inference from RNA sequencing (RNA-seq) data is foundational for molecular discovery. Here we present TACO, a computational method to reconstruct a consensus transcriptome from multiple RNA-seq data sets. TACO employs novel change-point detection to demarcate transcript start and end sites, leading to improved reconstruction accuracy compared with other tools in its class. The tool is available at http://tacorna.github.io and can be readily incorporated into RNA-seq analysis workflows.

RNA biomarkers associated with metastatic progression in prostate cancer: a multi-institutional high-throughput analysis of SChLAP1.

BACKGROUND: Improved clinical predictors for disease progression are needed for localised prostate cancer, since only a subset of patients develop recurrent or refractory disease after first-line treatment. Therefore, we undertook an unbiased analysis to identify RNA biomarkers associated with metastatic progression after prostatectomy. METHODS: Prostate cancer samples from patients treated with radical prostatectomy at three academic institutions were analysed for gene expression by a high-density Affymetrix GeneChip platform, encompassing more than 1 million genomic loci. In a discovery cohort, all protein-coding genes and known long non-coding RNAs were ranked by fold change in expression between tumours that subsequently metastasised versus those that did not. The top ranked gene was then validated for its prognostic value for metastatic progression in three additional independent cohorts. 95% of the gene expression assays were done in a Clinical Laboratory Improvements Amendments certified laboratory facility. All genes were assessed for their ability to predict metastatic progression by receiver-operating-curve area-under-the-curve analyses. Multivariate analyses were done for the primary endpoint of metastatic progression, with variables including Gleason score, preoperative prostate-specific antigen concentration, seminal vesicle invasion, surgical margin status, extracapsular extension, lymph node invasion, and expression of the highest ranked gene. FINDINGS: 1008 patients were included in the study: 545 in the discovery cohort and 463 in the validation cohorts. The long non-coding RNA SChLAP1 was identified as the highest-ranked overexpressed gene in cancers with metastatic progression. Validation in three independent cohorts confirmed the prognostic value of SChLAP1 for metastatic progression. On multivariate modelling, SChLAP1 expression (high vs low) independently predicted metastasis within 10 years (odds ratio [OR] 2·45, 95% CI 1·70-3·53; p<0·0001). The only other variable that independently predicted metastasis within 10 years was Gleason score (8-10 vs 5-7; OR 2·14, 95% CI 1·77-2·58; p<0·0001). INTERPRETATION: We identified and validated high SChLAP1 expression as significantly prognostic for metastatic disease progression of prostate cancer. Our findings suggest that further development of SChLAP1 as a potential biomarker, for treatment intensification in aggressive prostate cancer, warrants future study. FUNDING: Prostate Cancer Foundation, National Institutes of Health, Department of Defense, Early Detection Research Network, Doris Duke Charitable Foundation, and Howard Hughes Medical Institute.

RNASTAR: an RNA STructural Alignment Repository that provides insight into the evolution of natural and artificial RNAs.

Automated RNA alignment algorithms often fail to recapture the essential conserved sites that are critical for function. To assist in the refinement of these algorithms, we manually curated a set of 148 alignments with a total of 9600 unique sequences, in which each alignment was backed by at least one crystal or NMR structure. These alignments included both naturally and artificially selected molecules. We used principles of isostericity to improve the alignments from an average of 83%-94% isosteric base pairs. We expect that this alignment collection will assist in a wide range of benchmarking efforts and provide new insight into evolutionary principles governing change in RNA structural motifs. The improved alignments have been contributed to the Rfam database.

Digital spatial profiling of intraductal papillary mucinous neoplasms: Toward a molecular framework for risk stratification.

The histopathologic heterogeneity of intraductal papillary mucinous neoplasms (IPMN) complicates the prediction of pancreatic ductal adenocarcinoma (PDAC) risk. Intratumoral regions of pancreaticobiliary (PB), intestinal (INT), and gastric foveolar (GF) epithelium may occur with either low-grade dysplasia (LGD) or high-grade dysplasia (HGD). We used digital spatial RNA profiling of dysplastic epithelium (83 regions) from surgically resected IPMN tissues (12 patients) to differentiate subtypes and predict genes associated with malignancy. The expression patterns of PB and GF lesions diverged from INT, suggesting that PB and GF arise from a common lineage. Transcriptional dysregulation within PB lesions mirrored that of PDAC, whereas INT and GF foci did not. Tumor necrosis factor/nuclear factor κB (TNF-NFκB) and cell cycle (cycling S and cycling G2-M) programs occurred with relative prominence in PB and INT subtypes, respectively. Together, this study delineates markers of high-risk IPMN and insights into malignant progression.

Oculus: faster sequence alignment by streaming read compression.

BACKGROUND: Despite significant advancement in alignment algorithms, the exponential growth of nucleotide sequencing throughput threatens to outpace bioinformatic analysis. Computation may become the bottleneck of genome analysis if growing alignment costs are not mitigated by further improvement in algorithms. Much gain has been gleaned from indexing and compressing alignment databases, but many widely used alignment tools process input reads sequentially and are oblivious to any underlying redundancy in the reads themselves. RESULTS: Here we present Oculus, a software package that attaches to standard aligners and exploits read redundancy by performing streaming compression, alignment, and decompression of input sequences. This nearly lossless process (> 99.9%) led to alignment speedups of up to 270% across a variety of data sets, while requiring a modest amount of memory. We expect that streaming read compressors such as Oculus could become a standard addition to existing RNA-Seq and ChIP-Seq alignment pipelines, and potentially other applications in the future as throughput increases. CONCLUSIONS: Oculus efficiently condenses redundant input reads and wraps existing aligners to provide nearly identical SAM output in a fraction of the aligner runtime. It includes a number of useful features, such as tunable performance and fidelity options, compatibility with FASTA or FASTQ files, and adherence to the SAM format. The platform-independent C++ source code is freely available online, at http://code.google.com/p/oculus-bio.

ChimeraScan: a tool for identifying chimeric transcription in sequencing data.

SUMMARY: Next generation sequencing (NGS) technologies have enabled de novo gene fusion discovery that could reveal candidates with therapeutic significance in cancer. Here we present an open-source software package, ChimeraScan, for the discovery of chimeric transcription between two independent transcripts in high-throughput transcriptome sequencing data. AVAILABILITY: http://chimerascan.googlecode.com CONTACT: cmaher@dom.wustl.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Analysis of the androgen receptor-regulated lncRNA landscape identifies a role for ARLNC1 in prostate cancer progression.

The androgen receptor (AR) plays a critical role in the development of the normal prostate as well as prostate cancer. Using an integrative transcriptomic analysis of prostate cancer cell lines and tissues, we identified ARLNC1 (AR-regulated long noncoding RNA 1) as an important long noncoding RNA that is strongly associated with AR signaling in prostate cancer progression. Not only was ARLNC1 induced by the AR protein, but ARLNC1 stabilized the AR transcript via RNA-RNA interaction. ARLNC1 knockdown suppressed AR expression, global AR signaling and prostate cancer growth in vitro and in vivo. Taken together, these data support a role for ARLNC1 in maintaining a positive feedback loop that potentiates AR signaling during prostate cancer progression and identify ARLNC1 as a novel therapeutic target.

The lncRNA landscape of breast cancer reveals a role for DSCAM-AS1 in breast cancer progression.

Molecular classification of cancers into subtypes has resulted in an advance in our understanding of tumour biology and treatment response across multiple tumour types. However, to date, cancer profiling has largely focused on protein-coding genes, which comprise <1% of the genome. Here we leverage a compendium of 58,648 long noncoding RNAs (lncRNAs) to subtype 947 breast cancer samples. We show that lncRNA-based profiling categorizes breast tumours by their known molecular subtypes in breast cancer. We identify a cohort of breast cancer-associated and oestrogen-regulated lncRNAs, and investigate the role of the top prioritized oestrogen receptor (ER)-regulated lncRNA, DSCAM-AS1. We demonstrate that DSCAM-AS1 mediates tumour progression and tamoxifen resistance and identify hnRNPL as an interacting protein involved in the mechanism of DSCAM-AS1 action. By highlighting the role of DSCAM-AS1 in breast cancer biology and treatment resistance, this study provides insight into the potential clinical implications of lncRNAs in breast cancer.

Analysis of long non-coding RNAs highlights tissue-specific expression patterns and epigenetic profiles in normal and psoriatic skin.

BACKGROUND: Although analysis pipelines have been developed to use RNA-seq to identify long non-coding RNAs (lncRNAs), inference of their biological and pathological relevance remains a challenge. As a result, most transcriptome studies of autoimmune disease have only assessed protein-coding transcripts. RESULTS: We used RNA-seq data from 99 lesional psoriatic, 27 uninvolved psoriatic, and 90 normal skin biopsies, and applied computational approaches to identify and characterize expressed lncRNAs. We detect 2,942 previously annotated and 1,080 novel lncRNAs which are expected to be skin specific. Notably, over 40% of the novel lncRNAs are differentially expressed and the proportions of differentially expressed transcripts among protein-coding mRNAs and previously-annotated lncRNAs are lower in psoriasis lesions versus uninvolved or normal skin. We find that many lncRNAs, in particular those that are differentially expressed, are co-expressed with genes involved in immune related functions, and that novel lncRNAs are enriched for localization in the epidermal differentiation complex. We also identify distinct tissue-specific expression patterns and epigenetic profiles for novel lncRNAs, some of which are shown to be regulated by cytokine treatment in cultured human keratinocytes. CONCLUSIONS: Together, our results implicate many lncRNAs in the immunopathogenesis of psoriasis, and our results provide a resource for lncRNA studies in other autoimmune diseases.

The landscape of long noncoding RNAs in the human transcriptome.

Long noncoding RNAs (lncRNAs) are emerging as important regulators of tissue physiology and disease processes including cancer. To delineate genome-wide lncRNA expression, we curated 7,256 RNA sequencing (RNA-seq) libraries from tumors, normal tissues and cell lines comprising over 43 Tb of sequence from 25 independent studies. We applied ab initio assembly methodology to this data set, yielding a consensus human transcriptome of 91,013 expressed genes. Over 68% (58,648) of genes were classified as lncRNAs, of which 79% were previously unannotated. About 1% (597) of the lncRNAs harbored ultraconserved elements, and 7% (3,900) overlapped disease-associated SNPs. To prioritize lineage-specific, disease-associated lncRNA expression, we employed non-parametric differential expression testing and nominated 7,942 lineage- or cancer-associated lncRNA genes. The lncRNA landscape characterized here may shed light on normal biology and cancer pathogenesis and may be valuable for future biomarker development.

Targeting the MLL complex in castration-resistant prostate cancer.

Resistance to androgen deprivation therapies and increased androgen receptor (AR) activity are major drivers of castration-resistant prostate cancer (CRPC). Although prior work has focused on targeting AR directly, co-activators of AR signaling, which may represent new therapeutic targets, are relatively underexplored. Here we demonstrate that the mixed-lineage leukemia protein (MLL) complex, a well-known driver of MLL fusion-positive leukemia, acts as a co-activator of AR signaling. AR directly interacts with the MLL complex via the menin-MLL subunit. Menin expression is higher in CRPC than in both hormone-naive prostate cancer and benign prostate tissue, and high menin expression correlates with poor overall survival of individuals diagnosed with prostate cancer. Treatment with a small-molecule inhibitor of menin-MLL interaction blocks AR signaling and inhibits the growth of castration-resistant tumors in vivo in mice. Taken together, this work identifies the MLL complex as a crucial co-activator of AR and a potential therapeutic target in advanced prostate cancer.

A novel RNA in situ hybridization assay for the long noncoding RNA SChLAP1 predicts poor clinical outcome after radical prostatectomy in clinically localized prostate cancer.

Long noncoding RNAs (lncRNAs) are an emerging class of oncogenic molecules implicated in a diverse range of human malignancies. We recently identified SChLAP1 as a novel lncRNA that demonstrates outlier expression in a subset of prostate cancers, promotes tumor cell invasion and metastasis, and associates with lethal disease. Based on these findings, we sought to develop an RNA in situ hybridization (ISH) assay for SChLAP1 to 1) investigate the spectrum of SChLAP1 expression from benign prostatic tissue to metastatic castration-resistant prostate cancer and 2) to determine whether SChLAP1 expression by ISH is associated with outcome after radical prostatectomy in patients with clinically localized disease. The results from our current study demonstrate that SChLAP1 expression increases with prostate cancer progression, and high SChLAP1 expression by ISH is associated with poor outcome after radical prostatectomy in patients with clinically localized prostate cancer by both univariate (hazard ratio = 2.343, P = .005) and multivariate (hazard ratio = 1.99, P = .032) Cox regression analyses. This study highlights a potential clinical utility for SChLAP1 ISH as a novel tissue-based biomarker assay for outcome prognostication after radical prostatectomy.

The long non-coding RNA PCAT-1 promotes prostate cancer cell proliferation through cMyc.

Long non-coding RNAs (lncRNAs) represent an emerging layer of cancer biology, contributing to tumor proliferation, invasion, and metastasis. Here, we describe a role for the oncogenic lncRNA PCAT-1 in prostate cancer proliferation through cMyc. We find that PCAT-1-mediated proliferation is dependent on cMyc protein stabilization, and using expression profiling, we observed that cMyc is required for a subset of PCAT-1-induced expression changes. The PCAT-1-cMyc relationship is mediated through the post-transcriptional activity of the MYC 3' untranslated region, and we characterize a role for PCAT-1 in the disruption of MYC-targeting microRNAs. To further elucidate a role for post-transcriptional regulation, we demonstrate that targeting PCAT-1 with miR-3667-3p, which does not target MYC, is able to reverse the stabilization of cMyc by PCAT-1. This work establishes a basis for the oncogenic role of PCAT-1 in cancer cell proliferation and is the first study to implicate lncRNAs in the regulation of cMyc in prostate cancer.

The IncRNAs PCGEM1 and PRNCR1 are not implicated in castration resistant prostate cancer.

Long noncoding RNAs (IncRNAs) are increasingly implicated in cancer biology, contributing to essential cancer cell functions such as proliferation, invasion, and metastasis. In prostate cancer, several lncRNAs have been nominated as critical actors in disease pathogenesis. Among these, expression of PCGEM1 and PRNCR1 has been identified as a possible component in disease progression through the coordination of androgen receptor (AR) signaling (Yang et al., Nature 2013, see ref. [1]). However, concerns regarding the robustness of these findings have been suggested. Here, we sought to evaluate whether PCGEM1 and PRNCR1 are associated with prostate cancer. Through a comprehensive analysis of RNA-sequencing data (RNA-seq), we find evidence that PCGEM1 but not PRNCR1 is associated with prostate cancer. We employ a large cohort of >230 high-risk prostate cancer patients with long-term outcomes data to show that, in contrast to prior reports, neither gene is associated with poor patient outcomes. We further observe no evidence that PCGEM1 nor PRNCR1 interact with AR, and neither gene is a component of AR signaling. Thus, we conclusively demonstrate that PCGEM1 and PRNCR1 are not prognostic lncRNAs in prostate cancer and we refute suggestions that these lncRNAs interact in AR signaling.