Background: Transcriptome analysis of cancer tissues has been instrumental in defining tumor subtypes, diagnostic signatures and cancer regulatory networks. Cancer transcriptomes are still predominantly analyzed at the level of gene expression. Few studies have addressed transcript-level variations, and most of these only looked at splice variants. Previously we introduced a k-mer based, reference-free method, DE-kupl, that performs differential analysis of RNA-seq data at the k-mer level, which enables distinguishing RNAs differing by a single nucleotide. Here we evaluate the significance of differential events discovered by this method in two independent lung adenocarcinoma RNA-seq datasets (N=583 and N=154). Results: Focusing on differential events in a tumor vs normal setting, we found events in endogenous repeats, alternative splicing and polyadenylation sites, long non-coding RNAs, retained introns and unmapped RNAs. Replicability was highly significant for most event classes (assessed by comparing to events shared between unrelated tumors). Overall about 160,000 differential k-mer contigs were shared between datasets, including a large set of sequences from hypervariable genes such as immunoglobulins, SFTP and mucin genes. Most interestingly, we identified a set of novel tumor-specific long non-coding RNAs in intergenic and intronic regions. We found that expressed endogenous transposons defined two major groups of patients (high/low repeat expression) with distinct clinical characteristic. A number of repeats, intronic RNAs and lincRNA achieved strong patient stratification in univariate or multivariate survival models. Finally, using antigen presentation prediction, we identified 55 contigs predicted to produce recurrent tumor-specific antigens. Conclusions: K-mer based RNA-seq analysis enables description of cancer transcriptomes at nucleotide precision, independently of prior transcript annotation. Application to lung cancer data uncovered events stemming from a wide variety of transcriptional and postranscriptional mechanisms. Among those events, a significant subset was replicable between cohorts, thus constituting novel RNA hallmarks of cancer. The code is available at: https://github.com/Transipedia/dekupl-lung-cancer-inter-cohort.

/pdf/the-contribution-of-uncharted-rna-sequences-to-tumor-3pjeea0cet.pdf

Nonparametric simulation to evaluate the performance of statistical methods for differential gene expression analysis of long non-coding RNA-seq data in cancer studies

The contribution of uncharted RNA sequences to tumor identity in lung adenocarcinoma

https://biblio.ugent.be/publication/8599291/file/8599326.pdf

Nonparametric simulation to evaluate the performance of statistical methods for differential gene expression analysis of long non-coding RNA-seq data in cancer studies

Citations

The contribution of uncharted RNA sequences to tumor identity in lung adenocarcinoma

Related Papers (5)

Comparative analysis of differential gene expression tools for RNA sequencing time course data.

A comparison of statistical methods for detecting differentially expressed genes from RNA-seq data.

NPEBseq: nonparametric empirical bayesian-based procedure for differential expression analysis of RNA-seq data

Discrete distributional differential expression (D3E)--a tool for gene expression analysis of single-cell RNA-seq data

Modeling and analysis of RNA-seq data: a review from a statistical perspective