A 2023 study published in Scientific Data (Nature Portfolio) by Liu and colleagues used comprehensive RNA sequencing of plasma exosomes to simultaneously profile mRNA, lncRNA, and circRNA in osteosarcoma patients and healthy controls — demonstrating the power of multi-RNA analysis from a single library preparation.
Osteosarcoma is the most common primary malignant bone tumor in children and adolescents. Plasma exosomal RNAs represent a promising non-invasive biomarker source, but the full spectrum of coding and non-coding RNAs packaged into osteosarcoma-derived exosomes had not been comprehensively characterized prior to this study.
Figure 1. Study overview and experimental design.
Plasma exosome isolation and characterization from osteosarcoma patients and healthy controls, RNA extraction, library construction, and sequencing workflow. From Liu et al. 2023 (CC BY 4.0).
Methodology: Plasma exosomes were isolated from 10 osteosarcoma patients and 5 healthy donors by ultracentrifugation. Total exosomal RNA was extracted, and RNA-seq libraries were constructed using Ribo-depletion (rRNA removal) without polyA selection — enabling simultaneous capture of mRNA, lncRNA, and circRNA. Sequencing was performed on an Illumina NovaSeq 6000 (PE150).
Reads were aligned to the reference genome (hg19) using STAR, and transcript quantification was performed using featureCounts. Differentially expressed RNAs were identified by edgeR with FDR < 0.05 and |log₂FC| > 1. circRNAs were detected and quantified using find_circ and CIRI2 with a minimum of 2 unique back-spliced reads.
Results: The study identified 1,754 lincRNAs, 7,096 known circRNAs, and 1,935 novel circRNAs in plasma exosomes. Comparative analysis between osteosarcoma and healthy controls revealed:
- 331 differentially expressed mRNAs — genes involved in extracellular matrix organization, cell adhesion, and cancer signaling pathways
- 132 differentially expressed lincRNAs — including known cancer-associated lncRNAs and novel candidates
- 489 differentially expressed circRNAs — demonstrating that circular RNAs form a substantial component of the exosomal transcriptome
These results demonstrate that comprehensive multi-RNA profiling from a single library can simultaneously capture dysregulated coding and non-coding transcripts relevant to cancer biology.

Figure 2. Differential expression of mRNAs, lincRNAs, and circRNAs in osteosarcoma exosomes.
Volcano plots and heatmaps showing 331 differentially expressed mRNAs, 132 lincRNAs, and 489 circRNAs identified in plasma exosomes of osteosarcoma patients vs. healthy controls. From Liu et al. 2023 (CC BY 4.0).



RNA species composition detected by OmniRNA-seq — 8 RNA types identified from a single library.
Multi-RNA expression profiling — simultaneous detection of mRNA, lncRNA, and small RNA dynamics across experimental conditions.
Representative OmniRNA-seq data outputs — size distribution, group separation, differential expression, and pathway enrichment.