Profiling of repetitive RNA sequences in the blood plasma of patients with cancer

Reggiardo, Roman E.; Maroli, Sreelakshmi Velandi; Peddu, Vikas; Davidson, Andrew E.; Hill, Alexander; LaMontagne, Erin; Aaraj, Yassmin Al; Jain, Miten; Chan, Stephen Y.; Kim, Daniel H.

doi:10.1038/s41551-023-01081-7

Download PDF

Article
Open access
Published: 31 August 2023

Profiling of repetitive RNA sequences in the blood plasma of patients with cancer

Nature Biomedical Engineering volume 7, pages 1627–1635 (2023)Cite this article

17k Accesses
2 Citations
176 Altmetric
Metrics details

Subjects

Abstract

Liquid biopsies provide a means for the profiling of cell-free RNAs secreted by cells throughout the body. Although well-annotated coding and non-coding transcripts in blood are readily detectable and can serve as biomarkers of disease, the overall diagnostic utility of the cell-free transcriptome remains unclear. Here we show that RNAs derived from transposable elements and other repeat elements are enriched in the cell-free transcriptome of patients with cancer, and that they serve as signatures for the accurate classification of the disease. We used repeat-element-aware liquid-biopsy technology and single-molecule nanopore sequencing to profile the cell-free transcriptome in plasma from patients with cancer and to examine millions of genomic features comprising all annotated genes and repeat elements throughout the genome. By aggregating individual repeat elements to the subfamily level, we found that samples with pancreatic cancer are enriched with specific Alu subfamilies, whereas other cancers have their own characteristic cell-free RNA profile. Our findings show that repetitive RNA sequences are abundant in blood and can be used as disease-specific diagnostic biomarkers.

Terminal modifications independent cell-free RNA sequencing enables sensitive early cancer detection and classification

Article Open access 02 January 2024

Non-coding RNAs as liquid biopsy biomarkers in cancer

Article 10 January 2022

A comprehensive characterization of the cell-free transcriptome reveals tissue- and subtype-specific biomarkers for cancer detection

Article Open access 21 April 2021

Main

Of the 3 billion base pairs in the human genome, approximately 75% are transcribed into RNA¹. The vast majority of these RNAs are not translated into proteins and are thus considered non-coding RNAs. Although non-coding RNAs such as microRNAs^2,3 and long non-coding RNAs^4,5 (lncRNAs) are well annotated, many other non-coding RNAs are generated throughout the genome, including RNAs transcribed from repeat elements such as transposable elements (TEs)⁶. There are over 5 million repeat element insertions in the human genome, with repeat sequences comprising roughly half the genomic sequence content⁷. TE RNAs in particular are aberrantly expressed in diseases such as cancer^{8,9,10,11,12,13}, highlighting their potential as abundant and specific biomarkers of disease^14,15.

Cell-free RNAs¹⁶ are released from cells that comprise the various tissues and organ systems throughout the human body^17,18. The diagnostic and prognostic potential of cell-free RNA is evidenced by the prediction of pre-eclampsia in pregnancy^19,20,21, and cell-free RNAs serve as biomarkers of diseases such as cancer^{22,23,24,25,26} and Alzheimer’s disease^27,28,29. Cell-free RNAs have been profiled predominantly via whole-exome RNA sequencing (RNA-seq), which precludes the detection of repeat-derived and other non-coding RNA, or ribosomal-RNA-depleted total RNA-seq³⁰. Studies using total RNA-seq for cell-free RNA have identified many well-annotated non-coding RNAs in human plasma, and a small fraction of repeat-derived cell-free RNAs (1–2%) in healthy individuals³¹. However, the diagnostic potential of the repeat-derived cell-free RNA transcriptome in the context of disease remains unknown.

Here we report that repeat-aware profiling of the cell-free RNA transcriptome (COMPLETE-seq) enables the in-depth characterization of disease-specific, repeat-derived cell-free RNAs, and the accurate classification of patients with cancer by leveraging the rich feature space of the repeat-derived cell-free RNA transcriptome. In marked contrast to the cell-free RNA transcriptomes of healthy individuals, patients with cancer show strong enrichment of TE- and other repeat-derived cell-free RNAs in their blood, with patients with pancreatic cancer showing high levels of short interspersed nuclear element (SINE)-derived cell-free RNAs from various Alu subfamily elements. We further show the generalizability of COMPLETE-seq by showing that repeat-aware classification of liver, lung, oesophageal, colorectal and stomach cancer cell-free RNA data³² shows improved performance compared with repeat-naive classifiers. Taken together, our results show that repeat-aware COMPLETE-seq profiling of the cell-free RNA transcriptome identifies a robust and dynamic repeat-element-derived RNA signature for the diagnosis of diseases such as cancer.

Results

COMPLETE-seq enables repeat-aware profiling of the cell-free RNA transcriptome

We developed the COMPLETE-seq technology to enable repeat-aware characterization of the cell-free RNA transcriptome. To generate cell-free RNA-seq data from human plasma, we leveraged a highly sensitive RNA-seq protocol that robustly detects both coding and non-coding RNAs⁵. Given that the human genome contains millions of repeat element insertions that have not been examined in the context of cell-free RNA, we created a custom transcriptome annotation for cell-free RNA quantification that incorporates both well-annotated coding and non-coding RNAs (that is, GENCODE) and over 5 million repeat element insertions found in the human genome (that is, RepeatMasker). We then aggregated the RNA signal from individual repeat element insertions to the subfamily element level³³, reducing the number of repeat features from over 5 million to approximately 15,000 repeat features for disease classification and other downstream analyses (Fig. 1a).

**Fig. 1: Cell-free RNA transcriptome profiling using repeat-aware COMPLETE-seq.**

Compared with using only well-annotated GENCODE coding and non-coding genes (repeat naive) for cell-free RNA quantification, the application of our COMPLETE-seq technology significantly enhanced the percentage of mapped reads in our cell-free RNA data from patients with pancreatic cancer (Fig. 1b and Extended Data Fig. 1). For healthy control cell-free RNA, however, the mapping rate difference between repeat-naive versus our repeat-aware approach was negligible (Fig. 1b). Notably, there were no significant differences between the total number of repeat subfamilies that were represented in the cell-free RNA of patients with pancreatic cancer and healthy individuals (Fig. 1c).

Both repeat-naive and repeat-aware quantification of the cell-free RNA transcriptomes of patients with pancreatic cancer and healthy individuals enabled robust, unsupervised disease identification in low-dimensional space via principal component analysis (Extended Data Fig. 1e,f). Compared with using only well-annotated GENCODE coding and non-coding genes for cell-free RNA quantification, however, the application of our repeat-aware technology increased the sample-to-sample correlation (Pearson) of cell-free RNA data from patients with pancreatic cancer (Extended Data Fig. 1c,d), indicating greater overall similarity when examining a more robust annotation of their cell-free RNA transcriptomes.

TE RNAs and other repeat element RNAs are enriched in pancreatic cancer cell-free RNA

To determine the repeat composition of cell-free RNA, we first examined repeat content at the superfamily level. In healthy individuals, less than 10% of the DESeq2-normalized cell-free RNA counts consistently corresponded to repeats. However, we found substantially larger fractions of repeat-derived cell-free RNAs across almost all patients with pancreatic cancer, with most of these cell-free RNAs being derived from SINE elements (Fig. 2a). We also found significant differences in the information content, as quantified by Shannon entropy³⁴, of protein-coding RNAs, lncRNAs, and long terminal repeat (LTR), SINE and simple repeat superfamilies (Fig. 2b). These differences in biotype or repeat superfamily transcriptome diversity suggest dynamic changes in both the abundance and identity of cell-free RNAs in the context of diseases such as cancer.

**Fig. 2: Disease-specific enrichment of repeat-derived cell-free RNA.**

To further characterize the features of repeat-derived cell-free RNAs, we used full-length complementary DNA to perform single-molecule sequencing using nanopore technology³⁵ on cell-free RNA samples from patients with pancreatic cancer that were also sequenced using Illumina technology. We examined the size distributions of protein-coding RNA, lncRNA, and long interspersed nuclear element (LINE)-, SINE- and LTR-derived cell-free RNAs from patients with pancreatic cancer, which revealed cell-free RNA transcripts up to 1,337 nt in length for protein-coding RNA (median 456 nt), 970 nt for lncRNA (median 303 nt), 1,002 nt for LINE (median 258 nt), 368 nt for SINE (median 185 nt) and 477 nt for LTR (median 167 nt) (Extended Data Fig. 2a). For SINE-derived cell-free RNA, we observed a bimodal length distribution, reflecting both full-length, ~300-nt-long Alu-derived RNA, along with a shorter species of Alu-derived RNA (Extended Data Fig. 2a). We then compared the expected length of SINE-derived RNA based on genomic alignment with the observed length of cell-free SINE RNA and found both full-length transcripts of expected size, and half-length cell-free SINE RNAs (Extended Data Fig. 2b,c). In addition, we compared the cell-free RNA abundances of all COMPLETE-seq-annotated genes and repeat subfamilies in matched nanopore and Illumina libraries, which revealed strong concordance between well-annotated coding and non-coding genes, with a bias towards Illumina for TE RNAs and towards nanopore for some simple repeat RNAs (Extended Data Fig. 3).

We next performed differential expression analysis and found that Alu subfamily elements were the most enriched TE signal in cell-free RNA from patients with pancreatic cancer, with AluY, AluSc, AluSg7, AluSc8, AluSx3 and AluSg subfamily elements among the most significantly enriched (P < 0.01) in patients with pancreatic cancer compared with healthy individuals (Fig. 2c). Further analysis of the significantly enriched repeat subfamilies showed that the upregulated Alu elements in pancreatic cancer cell-free RNA were highly abundant (Extended Data Fig. 1g), despite the lack of increase in overall SINE complexity in the pancreatic cancer cell-free RNA transcriptome (Fig. 2b). Capturing both the robust enrichment of Alu elements and the significant increase in simple repeat entropy in the pancreatic cancer cell-free RNA transcriptome, hierarchical clustering achieved perfect clustering of patients with pancreatic cancer (Fig. 2d). These repeat-aware analyses further contextualized the differences in the respective repeat superfamilies, where SINE-derived cell-free RNAs were uniform in their enrichment in our cohort of patients with pancreatic cancer, whereas simple repeat-derived cell-free RNAs were far more divergent, with some simple repeat-derived cell-free RNAs being enriched in healthy individuals (Fig. 2d).

COMPLETE-seq reveals cancer-specific repeat element RNA signatures in cell-free RNA

To show the generalizability and applicability of COMPLETE-seq technology in enabling RNA liquid biopsy for cancer diagnosis, we used COMPLETE-seq quantification to analyse lung, liver, oesophageal, colorectal and stomach cancer cell-free RNA-seq data, along with their corresponding healthy controls³². We observed repeat superfamily variability in both healthy and cancer cell-free RNA transcriptomes (Extended Data Fig. 4 and Extended Data Fig. 5) and a significant (P < 0.05) increase in mapping rate by using our repeat-aware COMPLETE-seq annotation for analysing oesophageal, liver and stomach cancer cell-free RNA transcriptomes (Extended Data Fig. 4b).

Performing pairwise comparisons between five different cancers and healthy individuals captured robust and significant (P < 0.01) differential expression of repeat-derived cell-free RNAs that were characteristic to each cancer type (Fig. 3a–e). Additional analyses of the significantly differentially expressed repeat subfamilies showed that these repeat-derived cell-free RNAs were also highly abundant (Extended Data Fig. 6), with significant changes to biotype or repeat superfamily entropy (P < 0.05) (Extended Data Fig. 6f). By comparing the significantly differentially expressed repeat RNA signal across cancer types, we identified cancer-specific TE- and other repeat-derived cell-free RNA enrichment or depletion across all repeat superfamilies (Fig. 3f,g).

**Fig. 3: Disease-specific repeat-derived cell-free RNA signatures.**

COMPLETE-seq features improve classification performance of diagnostic models

To show proof of concept for diagnostic modelling using repeat-aware COMPLETE-seq analysis of cell-free RNA-seq data³², we trained logistic regression classifiers with tenfold cross-validation on training sets created for each cancer and healthy comparison (Fig. 4). To determine the utility of repeat-aware COMPLETE-seq features for disease classification, we trained models on repeat-naive and repeat-aware feature sets comprising DESeq2-normalized counts or biotype/repeat superfamily entropy for each cancer type. This resulted in eight feature sets comprising repeat-aware and repeat-naive counts, entropy and counts filtered by training set differential expression. Optimized repeat-aware models were compared with their repeat-naive counterparts, revealing repeat-driven increases in both area under the curve (AUC) (Fig. 4a,e,i,m,q) and training sensitivity for liver cancer (86% sensitivity) (Fig. 4b,c), oesophageal cancer (56% sensitivity) (Fig. 4f,g), colorectal cancer (91% sensitivity) (Fig. 4j,k), stomach cancer (86% sensitivity) (Fig. 4n,o) and lung cancer (93% sensitivity) (Fig. 4r,s) at 90% specificity.

**Fig. 4: COMPLETE-seq features improve performance of diagnostic models.**

Liver cancer, which showed a large repeat fraction (Extended Data Fig. 2a), had a corresponding dependence on repeat-aware features for classification, which resulted in a significant (P < 0.05) improvement in training sensitivity (Fig. 4a). Classification performance in the respective testing cohorts largely reflected the improvements seen in training, suggesting that our models have the potential to generalize to unseen data (Fig. 4c,g,k,o,s). Notably, we observed cancer-specific differences in repeat-aware feature dependence for disease classification. Stomach (Fig. 4p) and colorectal (Fig. 4l) cancer models each used one repeat-aware feature, liver (Fig. 4d) and oesophageal (Fig. 4h) cancer models used five and ten repeat-aware features, respectively, and the lung (Fig. 4t) cancer model used many repeat-aware features and the most overall features. For all five cancer models, repeat-aware features enhanced disease classification, highlighting the potential of COMPLETE-seq for highly sensitive and specific disease diagnosis.

Discussion

Our study reveals the value and utility of broadly characterizing the cell-free RNA transcriptome using our repeat-aware COMPLETE-seq technology for RNA liquid biopsies. Although other studies have provided valuable insights into protein-coding cell-free RNA dynamics, we find that the vast non-coding and repeat-derived cell-free RNA transcriptome is a rich source of abundant and disease-specific RNA biomarkers. We show that repeat-derived cell-free RNAs, including simple repeat RNAs and TE RNAs transcribed from LINE, SINE and LTR elements, are cancer-specific RNAs that are normally present at low or undetectable levels in healthy individuals. By creating a custom, repeat-aware transcriptome annotation for cell-free RNA quantification that incorporates over 5 million repeat element insertions throughout the human genome, we show that repeat-derived cell-free RNAs are highly enriched in the plasma of patients with cancer, with each cancer type showing its own characteristic repeat-derived cell-free RNA signature. COMPLETE-seq also greatly reduces the number of features used for downstream analysis and disease classification from over 5 million repeat element insertions to ~15,000 aggregated repeat elements at the subfamily level. Notably, our repeat-aware approach achieves highly accurate disease classification by incorporating both protein-coding RNAs and non-coding RNAs, such as lncRNAs and repeat-derived RNAs.

Although cell-free RNA studies so far have focused on short-read RNA-seq technologies, we show that long-read RNA-seq technologies, such as single-molecule nanopore sequencing, provide additional information regarding the full length of cell-free RNAs. We see differences in cell-free RNA length (that is, bimodal SINE-derived cell-free RNAs) that may serve as additional disease-specific features to further improve disease classification via RNA liquid biopsy (that is, RNA fragmentomics). Moreover, we also show that COMPLETE-seq robustly characterizes repeat-derived cell-free RNAs in both poly(A)-selected and total RNA library preparation protocols. In both cell-free RNA-seq contexts, COMPLETE-seq analysis increases mapping rate significantly and provides a richer feature space that leverages highly abundant and disease-specific repeat-derived cell-free RNAs to improve classification performance.

COMPLETE-seq also provides systemic insights into disease pathogenesis, and opportunities to discover drug targets for diseases such as cancer. In addition, our RNA liquid-biopsy technology enables non-invasive, systemic monitoring of protein-coding and repeat-derived cell-free RNA responses to targeted therapies, such as KRAS inhibitors³⁶, which induce treated cancer cells to preferentially release TE-derived cell-free RNAs in extracellular vesicles⁹. Given the preferential upregulation and secretion of TE-derived cell-free RNAs in response to mutant KRAS^8,9, companion diagnostic tests developed using repeat-aware RNA liquid biopsy would enable the robust detection of repeat-derived cell-free RNA signatures of oncogenic RAS signalling.

To move towards clinical implementation of COMPLETE-seq, future studies will require the generation of larger and more diverse cell-free RNA transcriptomic datasets across additional early-stage cancer types to further improve diagnostic performance and to accurately deconvolve the cell-free RNA transcriptome to determine cancer tissue of origin. Moreover, multi-cancer early detection using COMPLETE-seq will also require larger prospective studies to evaluate repeat-aware classification performance in an asymptomatic population. These studies may enable the application of repeat-aware RNA liquid-biopsy technology for the early detection of cancer and, more generally, for precision health.

Methods

Cell-free RNA isolation from blood plasma

The ExoRNeasy kit (Qiagen) was used to isolate cell-free RNA from blood plasma of de-identified healthy controls (blood collected in K2EDTA tubes; Discovery Life Sciences) and patients with pancreatic cancer (blood collected in K2EDTA tubes; BioIVT). Samples were initially filtered through a 0.8 µm filter to remove any contaminants, such as buffy coat. Filtered plasma was then processed using the ExoRNeasy kit to isolate cell-free RNA according to manufacturer instructions.

Library preparation for cell-free RNA-seq

Full-length cDNA was synthesized from cell-free RNA from pancreatic cancer and healthy control samples (Takara SMART-Seq HT kit). Size distribution of cell-free RNA and resulting cDNA were evaluated using an Agilent bioanalyser. Final libraries were made using the Illumina Nextera XT DNA Prep kit. These libraries were then sequenced (PE150) on an Illumina NextSeq 500.

Nanopore full-length cDNA libraries were prepared as described above, followed by manufacturer instructions for the Oxford Nanopore ligation kit LSK109. Sequencing for each library was performed on an Oxford Nanopore MinION device with R9.4 flow cells.

RNA-seq quality control, alignment and quantification

RNA-seq reads (FASTQ) were trimmed with Trimmomatic³⁷ (v.0.39), quality assessed using FastQC³⁸ (v.0.11.9) and visualized using MultiQC³⁹ (v.1.11). Quantification was performed using Salmon⁴⁰ (v.1.9.0) with two separate transcriptome annotations:

(1)
The GENCODE consortium Hg38 reference annotation (v.39): 61,488 genes
(2)
A concatenation of the above reference and the Hg38 repeat element track available at the University of California Santa Cruz genome browser mySQL server genome-mysql.cse.ucsc.edu: ~5 million insertions

The following optional arguments were used:

–validateMappings, –gcBias, –seqBias, –recoverOrphans, –rangeFactorizationBins 4

to enable selective alignment, reduce sequence biases, rescue reads with unmapped pairs and improved quantification, respectively.

The repeat-aware annotation aggregation was performed much as the transcript-to-gene aggregation is performed for alignment-free quantification approaches such as Salmon. Briefly, we assembled a reference transcriptome that includes all annotated repeats in Hg38 and associated each individual repeat instance ‘transcript’ (n ≈ 5 million) with its subfamily ‘gene’ (n ≈ 15,000). Repeat instances were summated to the subfamily level.

Nanopore signals were base-called with guppy (v.5.0.7) and alignments were quantified using StringTie2 (ref. ⁴¹) with the ‘-L’ argument for long reads using the repeat-aware reference described above (number 2).

Differential expression analysis

Salmon quantifications were loaded into R using tximeta (v.1.12.4)⁴² and converted to DESeq (v.1.34.0)^43,44 objects where counts were aggregated to the gene level from transcripts. Count normalization was performed with DESeq2 and pairwise comparisons were calculated with the following models:

Internal cohort: ~age + gender + input_volume + condition

External cohort: ~age + gender + condition

Significant differential expression was considered at adjusted P values (P_adj) of <0.01.

Unsupervised analysis

Using the gene-level, variance-stabilized counts computed above, principal component analysis was performed with the prcomp from the R stats (v.4.1.3) package on a count matrix filtered to include only genes with non-zero s.d. across the samples using the following optional arguments: centre = T scale. = T rank. = 50 to calculate the 50 principal components from centred, scaled counts.

Pearson correlation was calculated using the cor function from the R stats (v.4.1.3) package. K-means clustering was calculated and plotted with the ComplexHeatmap package in R using, where appropriate, Pearson correlation values (described above) or expression Z scores calculated with the scale function in R with ‘centre = T’ to centre our scaled values at zero.

Cell-free RNA length analysis

lncRNA, protein-coding and TE reference sequences were retrieved as described above. Sequences were extracted from Hg38 to create the biotype reference genomes used in the length analyses. Nanopore reads were aligned to Hg38 using Minimap2. To determine alignments in genomic regions with overlapping annotations, the length of the aligned fragment was compared with the lengths of the overlapping repeats. The annotation with the closest length to that of the fragment was chosen as the correct alignment. Fragment length was extracted using the PySam (v.0.15.4) template length.

Modelling and statistical analysis

Feature engineering

DESeq2-normalized counts of features with a non-zero s.d. were used directly in all cases except where they were used to calculate Shannon entropy. Entropy (H) was calculated on a per-sample basis for each biotype or subfamily as:

$${H}_{{\mathrm{biotype}}}=-\mathop{\sum }\limits_{i}^{n}{p}_{i}{{\rm{log }}}_{2}({p}_{i})$$

where p_i represents the fractional contribution of a given feature i (total of n) belonging to the biotype of interest to the total biotype abundance.

For classification as described below, eight feature sets belonging to three categories were used as input matrices for model training:

(1)
Total: repeat-naive, repeat-aware and repeat-alone features
(2)
Differential expression: repeat-naive, repeat-aware and repeat-alone differentially expressed features (calculated on training set, excluding test set)
(3)
Entropy: TE clade entropy and TE clade entropy plus repeat-naive features

Classification and performance evaluation

Each disease cohort was paired with healthy samples and split into stratified training (80%) and testing (20%) subsets. Training splits were used to optimize logistic regression models by performing tenfold cross-validated classification using elastic net penalty values (ɑ) from zero (lasso) to one (ridge) to optimize feature selection via the regularization parameter (λ), producing a final model trained on the entire training split with selected features. Top-performing models were identified by training sensitivity at 90% specificity, which was determined by calculating the probability threshold that achieved ~90% specificity and AUC.

When feature sets including differentially expressed genes were used, differential expression was calculated using only the training split and excluding testing samples. Model performance was finally evaluated on held-out test data by generating prediction probabilities on the test split samples and classifying based on the 90% specificity probability threshold defined in training. Features with non-zero coefficients (β) in the final models were identified to determine total feature size. Confidence intervals for sensitivity were estimated as binomial confidence intervals based on the successful observations and the total training/testing cohort. All modelling was performed with the cv.glmnet function from the glmnet package and custom code written in R.

Statistical analysis

Unless otherwise stated, comparison of means was performed with a two-sided, unpaired Wilcoxon rank-sum test. Where paired tests are used, lines are drawn to connect the dependent observations. When represented using symbolic ranks (*), statistical significance is defined as follows: non-significant P > 0.05, *P ≤ 0.05, **P ≤ 0.01, ***P ≤ 0.001 and ****P ≤ 0.0001.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Data availability

The main data supporting the results in this study are available within the article and its Supplementary Information. RNA-seq data are available at the NCBI Gene Expression Omnibus repository, under accession number GSE136651. Publicly available data used in this study are available at the NCBI Gene Expression Omnibus repository, under accession number GSE174302. All data generated in this study, including source data for the figures, are available from the corresponding author on reasonable request.

Code availability

Custom code used in this study is available on GitHub at https://github.com/rreggiar/exRNA_disease_biomarkers.

References

Djebali, S. et al. Landscape of transcription in human cells. Nature 489, 101–108 (2012).
Article PubMed PubMed Central CAS Google Scholar
Bartel, D. P. MicroRNAs: target recognition and regulatory functions. Cell 136, 215–233 (2009).
Article PubMed PubMed Central CAS Google Scholar
Kim, D. H., Saetrom, P., Snove, O. Jr & Rossi, J. J. MicroRNA-directed transcriptional gene silencing in mammalian cells. Proc. Natl Acad. Sci. USA 105, 16230–16235 (2008).
Article PubMed PubMed Central CAS Google Scholar
Rinn, J. L. & Chang, H. Y. Genome regulation by long noncoding RNAs. Annu. Rev. Biochem. 81, 145–166 (2012).
Article PubMed CAS Google Scholar
Kim, D. H. et al. Single-cell transcriptome analysis reveals dynamic changes in lncRNA expression during reprogramming. Cell Stem Cell 16, 88–101 (2015).
Article PubMed PubMed Central CAS Google Scholar
Elbarbary, R. A., Lucas, B. A. & Maquat, L. E. Retrotransposons as regulators of gene expression. Science 351, aac7247 (2016).
Article PubMed PubMed Central Google Scholar
Fernandes, J. D. et al. The UCSC repeat browser allows discovery and visualization of evolutionary conflict across repeat families. Mob. DNA 11, 13 (2020).
Article PubMed PubMed Central Google Scholar
Reggiardo, R. E. et al. Epigenomic reprogramming of repetitive noncoding RNAs and IFN-stimulated genes by mutant KRAS. Preprint at bioRxiv https://doi.org/10.1101/2020.11.04.367771 (2020).
Khojah, R. et al. Extracellular RNA signatures of mutant KRAS(G12C) lung adenocarcinoma cells. Preprint at bioRxiv https://doi.org/10.1101/2022.02.23.481574 (2022).
Burns, K. H. Transposable elements in cancer. Nat. Rev. Cancer 17, 415–424 (2017).
Article PubMed CAS Google Scholar
Kong, Y. et al. Transposable element expression in tumors is associated with immune infiltration and increased antigenicity. Nat. Commun. 10, 5228 (2019).
Article PubMed PubMed Central Google Scholar
Reggiardo, R. E. et al. Mutant KRAS regulates transposable element RNA and innate immunity via KRAB zinc-finger genes. Cell Rep. 40, 111104 (2022).
Article PubMed PubMed Central CAS Google Scholar
Carrillo, D. et al. Transposable element RNA dysregulation in mutant KRAS(G12C) 3D lung cancer spheroids. Preprint at bioRxiv https://doi.org/10.1101/2023.02.27.530369 (2023).
Reggiardo, R. E., Maroli, S. V. & Kim, D. H. lncRNA biomarkers of inflammation and cancer. Adv. Exp. Med. Biol. 1363, 121–145 (2022).
Article PubMed CAS Google Scholar
Wang, J., Ma, P., Kim, D. H., Liu, B. F. & Demirci, U. Towards microfluidic-based exosome isolation and detection for tumor therapy. Nano Today https://doi.org/10.1016/j.nantod.2020.101066 (2021).
Lo, Y. M. & Chiu, R. W. The biology and diagnostic applications of plasma RNA. Ann. N. Y. Acad. Sci. 1022, 135–139 (2004).
Article PubMed CAS Google Scholar
Vorperian, S. K., Moufarrej, M. N. & Quake, S. R. Cell types of origin of the cell-free transcriptome. Nat. Biotechnol. 40, 855–861 (2022).
Article PubMed PubMed Central CAS Google Scholar
Anfossi, S., Babayan, A., Pantel, K. & Calin, G. A. Clinical utility of circulating non-coding RNAs—an update. Nat. Rev. Clin. Oncol. 15, 541–563 (2018).
Article PubMed Google Scholar
Moufarrej, M. N. et al. Early prediction of preeclampsia in pregnancy with cell-free RNA. Nature 602, 689–694 (2022).
Article PubMed PubMed Central CAS Google Scholar
Rasmussen, M. et al. RNA profiles reveal signatures of future health and disease in pregnancy. Nature 601, 422–427 (2022).
Article PubMed PubMed Central CAS Google Scholar
Ng, E. K. et al. The concentration of circulating corticotropin-releasing hormone mRNA in maternal plasma is increased in preeclampsia. Clin. Chem. 49, 727–731 (2003).
Article PubMed CAS Google Scholar
Larson, M. H. et al. A comprehensive characterization of the cell-free transcriptome reveals tissue- and subtype-specific biomarkers for cancer detection. Nat. Commun. 12, 2357 (2021).
Article PubMed PubMed Central CAS Google Scholar
Hulstaert, E. et al. Charting extracellular transcriptomes in the Human Biofluid RNA Atlas. Cell Rep. 33, 108552 (2020).
Article PubMed CAS Google Scholar
Roskams-Hieter, B. et al. Plasma cell-free RNA profiling distinguishes cancers from pre-malignant conditions in solid and hematologic malignancies. npj Precis. Oncol. 6, 28 (2022).
Article PubMed PubMed Central CAS Google Scholar
Lo, K.-W. et al. Analysis of cell-free Epstein–Barr virus-associated RNA in the plasma of patients with nasopharyngeal carcinoma. Clin. Chem. 45, 1292–1294 (1999).
Article PubMed CAS Google Scholar
Kopreski, M. S., Benko, F. A., Kwak, L. W. & Gocke, C. D. Detection of tumor messenger RNA in the serum of patients with malignant melanoma. Clin. Cancer Res. 5, 1961–1965 (1999).
PubMed CAS Google Scholar
Koh, W. et al. Noninvasive in vivo monitoring of tissue-specific global gene expression in humans. Proc. Natl Acad. Sci. USA 111, 7361–7366 (2014).
Article PubMed PubMed Central CAS Google Scholar
Toden, S. et al. Noninvasive characterization of Alzheimer’s disease by circulating, cell-free messenger RNA next-generation sequencing. Sci. Adv. https://doi.org/10.1126/sciadv.abb1654 (2020).
Yan, Z. et al. Presymptomatic increase of an extracellular RNA in blood plasma associates with the development of Alzheimer’s disease. Curr. Biol. 30, 1771–1782 (2020).
Article PubMed CAS Google Scholar
Moufarrej, M. N., Wong, R. J., Shaw, G. M., Stevenson, D. K. & Quake, S. R. Investigating pregnancy and its complications using circulating cell-free RNA in women’s blood during gestation. Front. Pediatr. 8, 605219 (2020).
Article PubMed PubMed Central Google Scholar
Yao, J., Wu, D. C., Nottingham, R. M. & Lambowitz, A. M. Identification of protein-protected mRNA fragments and structured excised intron RNAs in human plasma by TGIRT-seq peak calling. eLife 9, e60743 (2020).
Article PubMed PubMed Central CAS Google Scholar
Chen, S. et al. Cancer type classification using plasma cell-free RNAs derived from human and microbes. eLife https://doi.org/10.7554/eLife.75181 (2022).
Bourque, G. et al. Ten things you should know about transposable elements. Genome Biol. 19, 199 (2018).
Article PubMed PubMed Central CAS Google Scholar
Garcia-Nieto, P. E., Wang, B. & Fraser, H. B. Transcriptome diversity is a systematic source of variation in RNA-sequencing data. PLoS Comput. Biol. 18, e1009939 (2022).
Article PubMed PubMed Central CAS Google Scholar
Deamer, D., Akeson, M. & Branton, D. Three decades of nanopore sequencing. Nat. Biotechnol. 34, 518–524 (2016).
Article PubMed PubMed Central CAS Google Scholar
Moore, A. R., Rosenberg, S. C., McCormick, F. & Malek, S. RAS-targeted therapies: is the undruggable drugged? Nat. Rev. Drug Discov. 19, 533–552 (2020).
Article PubMed PubMed Central CAS Google Scholar
Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).
Article PubMed PubMed Central CAS Google Scholar
Andrews, S. FastQC v.0.11.9 (Babraham Bioinformatics, Babraham Institute, 2010).
Ewels, P., Magnusson, M., Lundin, S. & Käller, M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics 32, 3047–3048 (2016).
Article PubMed PubMed Central CAS Google Scholar
Patro, R., Duggal, G., Love, M. I., Irizarry, R. A. & Kingsford, C. Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 14, 417–419 (2017).
Article PubMed PubMed Central CAS Google Scholar
Kovaka, S. et al. Transcriptome assembly from long-read RNA-seq alignments with StringTie2. Genome Biol. 20, 278 (2019).
Article PubMed PubMed Central CAS Google Scholar
Love, M. I. et al. Tximeta: reference sequence checksums for provenance identification in RNA-seq. PLoS Comput. Biol. 16, e1007664 (2020).
Article PubMed PubMed Central CAS Google Scholar
Zhu, A., Ibrahim, J. G. & Love, M. I. Heavy-tailed prior distributions for sequence count data: removing the noise and preserving large differences. Bioinformatics 35, 2084–2092 (2019).
Article PubMed CAS Google Scholar
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).
Article PubMed PubMed Central Google Scholar

Download references

Acknowledgements

We thank members of the Kim Lab for helpful discussions. This work was supported by funds from the American Cancer Society (https://doi.org/10.53354/pc.gr.158353; D.H.K.) and the Baskin School of Engineering at UC Santa Cruz (D.H.K.). R.E.R. was supported by the National Institutes of Health (DK131504), S.V.M. was supported by the California Institute for Regenerative Medicine (EDUC4-12759), V.P. was supported by the Tobacco-Related Disease Research Program (T32DT4904), S.Y.C. was supported by the American Heart Association Established Investigator Award (18EIA33900027) and the National Institutes of Health (HL122596 and HL124021), and D.H.K. was supported by a Research Scholar Grant (RSG-22-099-01-CDP) from the American Cancer Society.

Author information

Miten Jain
Present address: Department of Bioengineering, Northeastern University, Boston, MA, USA
Miten Jain
Present address: Department of Physics, Northeastern University, Boston, MA, USA

Authors and Affiliations

Department of Biomolecular Engineering, University of California Santa Cruz, Santa Cruz, CA, USA
Roman E. Reggiardo, Vikas Peddu, Andrew E. Davidson, Alexander Hill, Erin LaMontagne, Miten Jain & Daniel H. Kim
Department of Molecular, Cell and Developmental Biology, University of California Santa Cruz, Santa Cruz, CA, USA
Sreelakshmi Velandi Maroli
Center for Pulmonary Vascular Biology and Medicine, Pittsburgh Heart, Lung, Blood Vascular Medicine Institute, Division of Cardiology, Department of Medicine, University of Pittsburgh School of Medicine, Pittsburgh, PA, USA
Yassmin Al Aaraj & Stephen Y. Chan
Institute for the Biology of Stem Cells, University of California Santa Cruz, Santa Cruz, CA, USA
Daniel H. Kim
Genomics Institute, University of California Santa Cruz, Santa Cruz, CA, USA
Daniel H. Kim
Center for Molecular Biology of RNA, University of California Santa Cruz, Santa Cruz, CA, USA
Daniel H. Kim
Canary Center at Stanford for Cancer Early Detection, Stanford University School of Medicine, Palo Alto, CA, USA
Daniel H. Kim

Authors

Roman E. Reggiardo
View author publications
You can also search for this author in PubMed Google Scholar
Sreelakshmi Velandi Maroli
View author publications
You can also search for this author in PubMed Google Scholar
Vikas Peddu
View author publications
You can also search for this author in PubMed Google Scholar
Andrew E. Davidson
View author publications
You can also search for this author in PubMed Google Scholar
Alexander Hill
View author publications
You can also search for this author in PubMed Google Scholar
Erin LaMontagne
View author publications
You can also search for this author in PubMed Google Scholar
Yassmin Al Aaraj
View author publications
You can also search for this author in PubMed Google Scholar
Miten Jain
View author publications
You can also search for this author in PubMed Google Scholar
Stephen Y. Chan
View author publications
You can also search for this author in PubMed Google Scholar
Daniel H. Kim
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

D.H.K. conceived of and designed the study. R.E.R. performed experiments, devised and performed computational analyses, and generated figures. S.V.M. and E.L. performed Illumina sequencing experiments. A.E.D. and A.H. performed RNA-seq analyses. V.P. and M.J. performed nanopore-sequencing experiments and generated data. Y.A.A. and S.Y.C. provided materials. D.H.K. supervised research. D.H.K. and R.E.R. wrote the article with input from all co-authors.

Corresponding author

Correspondence to Daniel H. Kim.

Ethics declarations

Competing interests

D.H.K. and R.E.R. are inventors on patent applications covering the methods and compositions to detect cancer using cell-free RNA submitted by the Regents of the University of California. D.H.K. and R.E.R. are founders and shareholders and D.H.K. is a board member of LincRNA Bio. S.Y.C. has served as a consultant to United Therapeutics and Acceleron Pharma. S.Y.C. has held research grants from Actelion, Bayer and Pfizer. S.Y.C. is a director, officer and shareholder of Synhale Therapeutics. S.Y.C. has submitted patent applications regarding metabolism in pulmonary hypertension.

Peer review

Peer review information

Nature Biomedical Engineering thanks the anonymous reviewer(s) for their contribution to the peer review of this work.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data

Extended Data Fig. 1 Performance overview of COMPLETE-seq on the internal cohort.

a, Comparison of age distributions between cohorts (Wilcoxon, two-sided, ns: p > 0.05) enter line, median; box limits, upper and lower quartiles; whiskers, 1.5x interquartile range. b, Number of samples, stratified by gender, in each cohort. c, Heatmap (K-means) of Pearson correlation between panc samples using Repeat-naïve quantification. d, Heatmap (K-means) of Pearson correlation between panc samples using Repeat-aware quantification. e, PCA dimensions 1 & 2 calculated using variance-stabilized, Repeat-naive quantifications for normal and panc samples. f, PCA dimensions 1 & 2 calculated using variance-stabilized, Repeat-aware quantifications for normal and panc samples. g, MA plot of log2FoldChange between panc and normal samples compared to log-scale baseMean derived from DESeq2. Significantly DE genes/subfamilies are full opacity and colour.

Extended Data Fig. 2 Nanopore sequencing of cell-free RNA reveals biotype-specific fragment-size patterns.

a, Distribution of cell-free RNA lengths in base pairs (bp) for GENCODE biotypes or - repeat superfamily elements in pancreatic (panc 6, panc 7) cancer patients. b, Density plots depicting the relationship between expected (genomic SINE locus length) and observed SINE cell-free RNA length in pancreatic (panc 6, panc 7) cancer patients. c, Cumulative distribution function plot of SINE cell-free RNA length empirically calculated in pancreatic (panc 6, panc 7) cancer patients.

Extended Data Fig. 3 Nanopore and Illumina show agreement in the quantification of most GENCODE-annotated genes.

a, Scatter plot depicting transcripts-per-million abundance for transcripts detected in matched nanopore and Illumina libraries from sample panc 7. Linear fit described. b, Scatter plot depicting transcripts-per-million abundance for transcripts detected in matched nanopore and Illumina libraries from sample panc 6. Linear fit described.

Extended Data Fig. 4 Repeat-aware analysis of cell-free RNA from 5 different cancers.

a, Distribution of biotype representation (by DESeq2 normalized count) in cell-free RNA-seq quantifications for each cancer type, coloured by GENCODE biotype or repeat subfamily, and facetted by stage. b, Comparison of mapping rates between use of a Repeat-naïve (GENCODE v39) reference annotation and Repeat-aware reference annotation (Wilcoxon, paired, two-sided). center line, median; box limits, upper and lower quartiles; whiskers, 1.5x interquartile range. ns: p > 0.05, *: p <= 0.05, **: p <= 0.01, ***: p <= 0.001.

Extended Data Fig. 5 Repeat-specific diversity across internal and external cohorts.

a, Distribution of repeat representation (by DESeq2 normalized count) in cell-free RNA-seq quantifications for pancreatic cancer, coloured by repeat subfamily, and facetted by stage. b, Distribution of repeat representation (by DESeq2 normalized count) in cell-free RNA-seq quantifications for each cancer type, colored by repeat subfamily.

Extended Data Fig. 6 Differential expression and variability of repeat-elements in 5 different cancers.

a-e, MA plots of log2FoldChange between labeled cancer type and healthy donor samples compared to log-scale baseMean derived from DESeq2. Significantly DE genes/repeat subfamilies are full opacity and color. f, Comparison of significantly different (Wilcoxon, two-sided) Shannon Entropy distributions for GENCODE biotypes and repeat subfamilies. center line, median; box limits, upper and lower quartiles; whiskers, 1.5x interquartile range. ns: p > 0.05, *: p <= 0.05, **: p <= 0.01, ***: p <= 0.001, ****: p <= 0.0001.

Supplementary information

Reporting Summary

Supplementary Information

Metadata describing critical clinical features of the 32-member cohort generated for this study, and available at https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE138651.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Reggiardo, R.E., Maroli, S.V., Peddu, V. et al. Profiling of repetitive RNA sequences in the blood plasma of patients with cancer. Nat. Biomed. Eng 7, 1627–1635 (2023). https://doi.org/10.1038/s41551-023-01081-7

Download citation

Received: 22 August 2022
Accepted: 26 July 2023
Published: 31 August 2023
Issue Date: December 2023
DOI: https://doi.org/10.1038/s41551-023-01081-7