Characterization and selective incorporation of small non-coding RNAs in non-small cell lung cancer extracellular vesicles
- Chuang Li†1,
- Fang Qin†1,
- Fen Hu2,
- Hui Xu2,
- Guihong Sun3,
- Guang Han4, 5,
- Tao Wang2Email author and
- Mingxiong Guo1Email authorView ORCID ID profile
© The Author(s) 2018
Received: 20 November 2017
Accepted: 4 January 2018
Published: 10 January 2018
Extracellular vesicles (EVs) play important roles in intercellular communication through the delivery of their cargoes, which include proteins, lipids, and RNAs. Increasingly, multiple studies have reported the association between EV small non-coding RNAs and cancer, due to their regulatory functions in gene expression. Hence, analysis of the features of small non-coding RNA expression and their incorporation into EVs is important for cancer research.
We performed deep sequencing to investigate the expression of small RNAs in plasma EVs from lung adenocarcinoma (ADC) patients, lung squamous cell carcinoma (SQCC) patients, and healthy controls. Then, eighteen differently expressed miRNAs in plasma EVs was validated by QRT-PCR. The small RNA expression profiles of plasma EVs were different among lung ADC, SQCC patients, and healthy controls. And many small RNAs, including 5′ YRNA hY4-derived fragments, miR-451a, miR-122-5p, miR-20a-5p, miR-20b-5p, miR-30b-5p, and miR-665, were significantly upregulated in non-small cell lung cancer (NSCLC) EVs. And the cell viability assays indicated that hY4-derived fragments inhibited the proliferation of lung cancer cell A549. By comparing the cellular and EV expression levels of six miRNAs in NSCLC cells, we found that miR-451a and miR-122-5p were significantly downregulated in NSCLC cell lysates, while significantly upregulated in NSCLC EVs.
The differently expressed EV small RNAs may serve as potential circulating biomarkers for the diagnosis of NSCLC. Particularly, YRNA hY4-derived fragments can serve as a novel class of biomarkers, which function as tumor suppressors in NSCLC. Additionally, miR-451a and miR-122-5p may be sorted into NSCLC EVs in a selective manner.
Extracellular vesicles (EVs) are membrane surrounded structures released by cells, and contain complex cargoes including nucleic acid, protein, and lipids . EVs are heterogeneous in size and intracellular origins, and the main two subtypes of EVs are exosomes and microvesicles. Exosomes are small EVs (< 150 nm) released through endosomal pathway, and microvesicles are 100–1000 nm vesicles budding from plasma membrane. Importantly, EVs, especially exosomes, are secreted by almost all types of cells in culture, and are abundant in various body fluids including blood, saliva, urine, and breast milk . It is widely confirmed that exosomes are originally generated by endocytosis and persist in late endosomes; invagination of the limiting late endosomal membrane results in the formation multivesicular bodies (MVBs), from which exosomes are finally released into extracellular space upon fusing of MVBs and the cell membrane . Exosomes play important roles in intercellular communication, through direct interaction of exosomal proteins with signaling receptors of target cells , or by delivery of their contents into the recipient cell cytosol following fusion with recipient cell membranes through endocytic pathways . Other studies have demonstrated that exosomes are involved in cancer initiation, progression, and metastasis . Among the many types of molecules carried in exosomes, small RNAs, especially miRNAs, have attracted the most attention due to their roles in gene expression regulation.
Small RNAs are small (< 200 nucleotides) non-coding RNA molecules that include microRNA (miRNA), small nuclear RNA (snRNA), small nucleolar RNA (snoRNA), piwi-interacting RNA (piRNA), transfer RNA (tRNA), small ribosome RNA (rRNA), and small cytoplasmic RNA (Y RNA). These specialized RNAs play important roles in many different biological processes, such as RNA silencing, transcriptional regulation, chromosome replication, RNA processing, and regulation of protein trafficking and degradation . It has increasingly become clear that non-coding RNAs can be sorted into exosomes, protected from RNase degradation, and internalized by neighboring or distant cells, where they subsequently modulate cellular processes .
Multiple EV small RNAs play important roles in tumorigenesis; however, the function of small RNAs in tumor-derived EV is complicated. For example, tumor-derived exosomes can transfer oncogenic small RNAs to recipient cells, and promote proliferation, migration, angiogenesis, and extracellular matrix remodeling, thus promoting the formation of a premetastatic niche that contributes to tumor metastasis . In contrast, cancer cell-derived exosomes can mediate secretion of tumor suppressor small RNAs during tumor progression as a mechanism of coordinately suppressing tumor metastasis [10, 11].
Based on their cancer-specific expression profiles and the ability to detect them non-invasively in patient blood, small non-coding RNA molecules are promising biomarkers for cancer diagnosis [12, 13]. For example, circulating exosomal snRNA RNU6-1 alone or in combination with miR-320/miR-574-3p/RNU6-1 is a potential biomarker for glioblastoma multiforme , and exosomal miR-21 is a diagnostic biomarker for breast cancer, hepatocellular carcinoma, and esophageal squamous cell carcinoma [15–17]. In pancreatic ductal adenocarcinoma, the exosomal miRNA signature is more diagnostically relevant than exosomal protein glypican-1 levels , illustrating the potential advantages of exosomal small RNAs as biomarkers for cancer diagnosis.
Lung cancer is the leading worldwide cause of cancer-related death for both men and women . The two main pathological types of lung cancer are NSCLC and small cell lung cancer (SCLC). Approximately 85–90% of all lung cancers are categorized as NSCLC , and the most common pathological types of NSCLC are adenocarcinoma (ADC; 30–50%) and squamous cell carcinoma (SQCC; 30%) . In China in 2015, there were an estimated 733,300 new cases and 610,200 deaths due to lung cancer . Similarly, in the USA, the average 5-year survival rate for lung cancer is only 18%, and more than one-half of lung cancer patients are diagnosed at a late stage, for which the 5-year survival rate is as low as 4% due to lack of effective treatments [19, 23]. Hence, investigating the genetic differences between NSCLC and healthy patients is important for understanding the pathogenesis of NSCLC and discovery of novel diagnostic biomarkers and therapeutic targets.
Adenocarcinoma and SQCC have clinical and pathological differences , and in order to identify specific diagnostic biomarkers or treatment targets for NSCLC subtypes, researchers have focused on genetic alterations or gene expression differences between lung ADC and SQCC [25, 26]. Importantly, although miRNA expression profiles are available, there is little information about the global expression differences of small non-coding RNAs in lung ADC and SQCC.
In recent years, it has become clear that exosomal miRNAs are closely associated with lung cancer. For example, miR-21, miR-155, miR-200b, and miR-379, are abnormally expressed in lung cancer-derived exosomes [27, 28], and tumor-derived exosomal miRNAs have shown promise as biomarkers for lung cancer subtypes, including ADC-specific miRNAs miR-181-5p, miR-30a-3p, miR-30e-3p, and miR-361-5p, as well as SQCC-specific miRNAs miR-10b-5p, miR-15b-5p and miR-320b . Importantly, however, a comprehensive analysis and extensive verification of circulating exosomal small RNAs in NSCLC subtypes is still necessary for further understanding of NSCLC pathogenesis.
In the present study, we performed high-throughput sequencing of small RNAs in plasma EVs from lung ADC patients, SQCC patients, and healthy controls (CTRL). We found that a variety of small non-coding RNA species are present in plasma EVs, and EV small RNA fragments are likely to be associated with NSCLC pathogenesis. YRNA hY4-derived fragments were significantly upregulated in plasma EVs from both ADC and SQCC patients, and may be promising biomarkers for NSCLC diagnosis. Furthermore, we identified that the RNY4P7 gene on chromosome 2 may not be a pseudogene, because the corresponding transcript detected in EVs and cells. Additionally, analysis of miRNAs demonstrated clear EV miRNA expression profile differences between ADC, SQCC, and CTRL groups, and suggested that certain miRNAs might be selectively sorted into NSCLC EVs.
Characterization and properties of plasma EVs
Clinical characteristics of plasma EV samples from NSCLC patients and healthy controls
Age of diagnosis (year, mean ± SD)
60.04 ± 7.89
63.08 ± 7.19
62.25 ± 10.21
60.15 ± 8.11
62.30 ± 7.68
53.38 ± 18.09
Characterization of EV particle diameter
Average of particle diameter (nm)
Polydispersity index (PDI)a
Major peak of particle diameter (nm)
Percentage of 20–200 nm diameter (%)
Small RNA expression profiles are different between NSCLC plasma EVs and controls
The data filtering statistics of RNA-Seq data from NSCLC and control EVs
Removed inferior quality
Smaller than 17 nt
All sample groups had a major tRNA fragment peak at 29–34 nt, but the abundance of tRNA fragments in the SQCC group was about threefold higher than in ADC and CTRL groups (Fig. 3a). The reads for rRNA fragments in all sample groups were mainly distributed at 20–32 nt, and the RPM of the CTRL group at 20–26 nt was higher than ADC and SQCC (Fig. 3b). The snRNA fragment reads were mainly distributed at 19–24 nt, and the RPM of the CTRL group was higher than in the ADC and SQCC groups at the longer snRNA fragment (> 24 nt) (Fig. 3c). We found no significant difference in snoRNA fragments at the major peak (about 27 nt) or minor peak (18–24 nt), but the RPM of the CTRL group was slightly higher than ADC and SQCC at 30–42 nt (Fig. 3d). In addition, the RPM for piRNA fragments in the CTRL group was significantly higher than in the ADC and SQCC groups at the major peak of 20–22 nt (Fig. 3e). All groups had a main YRNA peak at 30–33 nt, and the YRNA fragment abundance in ADC and SQCC groups at the main peak was significantly higher than CTRL (Fig. 3f). Together, our initial data indicates that the expression signature of EV small RNA species is different in NSCLC and healthy samples.
hY4 RNA fragments are upregulated in NSCLC EVs and inhibit NSCLC cell proliferation
We confirmed that a small RNA on chromosome 8 was miR-320a (TCGCCCTCTCAACCCAGCTTTT), which is consistent with EV miRNA expression profiles (Additional file 1: Table S1). Importantly, we found that miR-320a was the most abundant EV miRNA in our samples, and was significantly downregulated in SQCC EVs. The small RNAs that mapped to chromosomes 2, 3, 6, and 7 were confirmed as the hY4 fragment (hY4F, AGTTCTGATAACCCACTACCATCGGACCAGCC) and hY4 reverse fragment (hY4rF, GGCTGGTCCGATGGTAGTGGGTTATCAGAACT), which aligned to the 5′-end of the hY4 RNA sequence. hY4F mapped to the Y_RNA gene (predicted YRNA from the Rfam database) on chromosome 3, the RNY4P10 pseudogene on chromosome 6, and the RNY4 gene on chromosome 7. Additionally, hY4rF mapped to the RNY4P7 pseudogene on chromosome 2 (Fig. 4d–g).
It is generally acknowledged that YRNA genes are clustered together at a single locus on human chromosome 7q36 , however, 966 pseudogenes derived from YRNAs have been characterized in the human genome . Recent reports demonstrate that some hY4 sequences previously annotated as pseudogenes are transcribed, processed, and secreted, suggesting their pseudogene annotation may not be accurate . Our results are consistent with this, as most of the upregulated YRNA fragments in NSCLC EVs were identified as hY4 fragments. In addition, our results suggest that RNY4P7 may not be a hY4 RNA pseudogene, as hY4rF, which is fully aligned to RNY4P7, was detected in both NSCLC cells and normal cells (Fig. 4h).
Then, we detected expression of hY4 RNA fragments in NSCLC cells, and found that both hY4F and hY4rF were significantly down-regulated in all the NSCLC cells we tested (Fig. 4h). Furthermore, we found that overexpression of hY4F and hY4rF could significantly suppress proliferation of NSCLC cell A549 (Fig. 4i, j). The above results indicate that hY4 RNA fragments, including hY4F and hY4rF may function as tumor suppressors in NSCLC.
miRNAs are differentially expressed in NSCLC plasma EVs
We next compared the expression of all miRNAs identified by small RNA sequencing to determine whether there were any correlations between different EV sample types. Overall, we found a strong positive correlation between ADC and CTRL (r = 0.9490), and a weaker correlation between ADC and SQCC (r = 0.8039), and SQCC and CTRL (r = 0.7085) (P values < 0.0001) (Fig. 5b). This data collectively suggests that while there are many differentially expressed miRNAs between the three groups, the majority of EV miRNAs have similar expression.
The specific miRNAs that were expressed in plasma EVSs among ADC, SQCC patients, and CTRL group is shown in Additional file 1: Table S1, and the significantly changed miRNAs are listed in Additional file 1: Tables S2–S5. For ADC vs. CTRL, ADC EVs contained 37 miRNAs that were significantly upregulated, and 11 miRNAs that were significantly downregulated. For SQCC vs. CTRL, 98 miRNAs were significantly upregulated in SQCC, while 47 miRNAs were significantly downregulated. For NSCLC (average of ADC and SQCC) vs. CTRL, 78 miRNAs were significantly upregulated, while 27 miRNAs were significantly downregulated. For ADC vs. SQCC, 33 miRNAs were significantly upregulated in ADC EVS s, while 56 miRNAs were significantly downregulated. In summary, our data demonstrates that there are many distinct EV miRNA expression signatures in ADC and SQCC patients compared to CTRL group, which suggests that there may be certain miRNA candidates that may serve as specific biomarkers for either ADC and SQCC.
qRT-PCR verification of miRNAs in NSCLC plasma EVs
RNA-Seq results of EV miRNAs selected for qRT-PCR detection
Our above results indicate that expression of miR-20a-5p, miR-20b-5p, miR-30b-5p, miR-451a, and miR-665 are significantly upregulated in NSCLC plasma EVs, which is consistent with the RNA-Seq results. However, the expression of other 13 miRNAs tested by qRT-PCR was not significantly different among the three EV sample groups (Additional file 1: Figure S3), which was in contrast to our RNA-Seq results. These differences may be explained by individual variation of EV miRNAs due to biological heterogeneity  of lung tumors [38, 39].
Enrichment of certain microRNAs in NSCLC EVs
Next, we compared the expression of these miRNAs in EVs from each cell line (Fig. 8b). We found that miR-451a, miR-122-5p, miR-20a-5p, miR-20b-5p, miR-30b-5p, and miR-665 were significantly upregulated in EVs from most of the tested NSCLC cells, consistent with our RNA-Seq and qRT-PCR results. Finally, we compared the cellular expression of these miRNAs to the matched EV expression for each cell line. As shown in Fig. 8c and d, miR-451a and miR-122-5p were significantly higher in EVs from all cell types compared to matched cells, suggesting that these miRNAs may be selectively sorted into NSCLC EVs. We also found that the EV vs. cellular ratio fold change expression of miR-451a in all NSCLC cell lines, except for H1975, was higher than in control cells (Fig. 8d), further indicated its potential to be selectively sorted into EVs.
GO and KEGG analysis of differentially expressed miRNA
To investigate the potential biological functions and pathways in which differentially expressed NSCLC plasma EV miRNAs are involved, we performed GO (Gene Ontology) biological process enrichment analysis and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis with the predicted miRNA targets. The GO analysis results (Additional file 1: Figure S4) showed that the predicted targets of miR-20a-5p, miR-20b-5p, miR-30b-5p, and miR-122-5p were similarly enriched for GO biological process terms, including regulation of macromolecule metabolic processes, regulation of nucleobase, nucleoside, nucleotide and nucleic acid metabolic processes, regulation of nitrogen compound metabolic processes, and regulation of macromolecule biosynthetic processes. However, we found that the GO enrichment results of predicted miR-451a and miR-665 targets were extremely different. The predicted targets of miR-451a were enriched in regulation of catalytic activity, regulation of cell proliferation, regulation of phosphate metabolism, and regulation of phosphorylation. In contrast, the predicted miR-665 targets had major enrichment in development-related biological processes, such as respiratory system development, nervous system development, lung development, respiratory tube development, blood vessel morphogenesis, and embryonic skeletal system development. Together, these results demonstrate that these EV miRNAs may have both common and specific functions, suggesting that complicated regulation networks exist among potential cancer-related EV miRNAs.
The KEGG results (Additional file 1: Figure S5) show that the predicted targets of some miRNAs were significantly enriched in common pathways. For example, the predicted targets of all analyzed miRNAs, except miR-665, were enriched for endocytosis (kegg: hsa04144), the predicted targets of miR-122-5p, miR-20a-5p, miR-20b-5p, miR-30b-5p, and miR-665 were enriched for cancer pathways (kegg: hsa05205), and miR-122-5p, miR-20a-5p, miR-20b-5p, and miR-665 were enriched for cancer proteoglycans (kegg: hsa05205). Additionally, the KEGG results showed that predicted targets of these miRNAs are enriched for multiple pathways involved in NSCLC, such as MAPK signaling (kegg: hsa04010), Wnt signaling (kegg: hsa04310), and FoxO signaling (kegg: hsa04068), suggesting that these miRNAs are likely to regulate signaling pathways involved in NSCLC pathogenesis.
Our study has demonstrated that miR-451a and miR-122-5p are selectively sorted into NSCLC EVs, and given the importance of the endocytic pathway in the selective sorting and uptake of EVs [5, 43], we analyzed the KEGG enrichment of predicted miR-451a and miR-122-5p targets in the endocytic pathway. As shown in Additional file 1: Figure S6, predicted targets of miR-451a and miR-122-5p are enriched for both clathrin-dependent and clathrin-independent endocytosis. In addition, some enriched targets of either miRNA are the same genes, or belong to the same gene families, including CAV1/CAV2, KIF5B/KIF5C, ARPC2/ARPC3, and GIT1/GIT2. Together, the above results indicate that miR-122-5p and miR-451a are likely to regulate the endocytic pathway by similar mechanisms.
Since the presence of miRNAs in EVs was published about 10 years ago , multiple types of small non-coding RNAs have been found in EVs , and many studies have indicated that some miRNAs are differentially expressed in cancer exosomes, and are closely related to tumorigenesis [8, 9]. Circulating EV small RNAs, especially miRNAs, are also potential non-invasive biomarkers for cancer diagnosis [13–15]. Overall, we have demonstrated that the small RNA expression profiles of plasma EVs are different among lung ADC, SQCC patients, and healthy controls. In particular, YRNA fragments are the most abundant small RNA subtype, and are more abundant in NSCLC plasma EVs than healthy controls. Our analyses suggest that most small RNA types, including YRNA, tRNA, piRNA, rRNA, and snRNA fragments, have different expression and distribution signatures in EVs from NSCLC patients and healthy controls, and are likely to be involved in NSCLC pathogenesis. Further analysis of YRNA fragments indicates that the majority of upregulated YRNA fragments in NSCLC plasma EVs are hY4 RNA-derived fragments. Interestingly, we found that the RNY4P7 gene on chromosome 2, corresponding to hY4rF, may not be a pseudogene. Besides, we found that hY4-derived fragments, hY4F and hY4rF, could inhibit NSCLC cell proliferation. Our analysis of miRNAs indicated that some miRNAs expressed in plasma EVs differ among ADC, SQCC, and CTRL groups. We further found that miR-451a, miR-122-5p, miR-20a-5p, miR-20b-5p, miR-30b-5p, and miR-665 are significantly upregulated in NSCLC EVs, suggesting that these miRNAs may serve as potential circulating biomarkers for the diagnosis of NSCLC. We also found that expression of miR-451a and miR-122-5p was significantly downregulated in NSCLC cell lysates, while significantly upregulated in NSCLC EVs, indicating that miR-451a and miR-122-5p may be selectively sorted into NSCLC EVs.
Differentially expressed YRNA fragments in NSCLC EVs
We found that EV YRNA expression was significantly higher in NSCLC than in healthy individuals, suggesting that EV YRNA may be a potential biomarker for NSCLC diagnosis. YRNAs are 84–113 nt small non-coding RNAs with stem-loop structures that bind Ro60 and La proteins to form the Ro ribonucleoprotein complex . There are four YRNAs in human: hY1, hY3, hY4, and hY5. YRNAs are upregulated in some human tumor tissues, and hY1 and hY3 are required for cell proliferation , suggesting that YRNAs are closely associated with carcinogenesis. YRNA fragments have also been found in EVs from both normal and cancer cells , as well as in human serum and plasma , suggesting a potential function in EV-mediated intercellular communication. Other studies have identified the presence of YRNA fragments in EVs from both breast cancer and normal breast cells, however there was no significant difference in YRNA expression between the two types of EVs . Although it is generally acknowledged that YRNAs are clustered together at a single locus on human chromosome 7q36 , some studies show that YRNAs also have pseudogenes on other chromosomes in the human genome [35, 36]. A recent report showed that over 85% of YRNA-derived fragments in EVs from human breast cancer cells and non-tumor breast cells corresponded to hY4 RNA . Similarly, we found that most of the YRNA fragments upregulated in plasma EVs of NSCLC patients were hY4 RNA fragments, and the RNY4P7 gene on chromosome 2 may not be a hY4 RNA pseudogene. Additionally, we found that hY4-derived fragments (hY4F and hY4rF) were significantly down-regulated in NSCLC cells, and inhibited the proliferation of NSCLC cell A549, suggesting that hY4-derived fragments play important roles in NSCLC genesis and progress.
Expression of other small RNAs in NSCLC EVs
We also found that other types of small RNAs in plasma EVs have different expression profiles among ADC, SQCC, and CTRL groups. For example, SQCC EVs had higher levels of tRNA fragments, and NSCLC EVs had lower levels of rRNA, piRNA, and snRNA fragments than in EVs from healthy controls. tRNA fragments are the cleavage forms of the mature tRNAs or tRNA precursors, and are heterogeneous in size, nucleotide composition, biogenesis, and function . tRNA fragments are associated with breast cancer  and prostate cancer  progression, and may act as miRNA-like molecules and as post-transcriptional regulators . Studies have reported that tRNA-derived fragments are present in EVs from human semen , plasma , and breast cancer cells , and the tRF-Leu-CAG tRNA fragment, which is upregulated in NSCLC and NSCLC cell lines, can promote cell proliferation . Similarly, our studies suggest that EV tRNA fragments may be associated with NSCLC pathogenesis, particularly in SQCC.
In addition to tRNAs, our study also demonstrated that rRNA, snRNA, and piRNA have differential expression signatures in EVs from NSCLC patients compared to healthy controls, and therefore, may serve as potential biomarkers for NSCLC. snRNAs are a class of small non-coding RNAs (about 150 nt), that form the main component of the spliceosome complex . Of interest, circulating U2 snRNA fragments can be a diagnostic biomarker for ovarian cancer , pancreatic cancer, and colorectal cancer . piRNAs are 26–31 nt small non-coding RNA molecules that form RNA–protein complexes through interactions with piwi proteins ; they are involved in gene silencing  and epigenetics , and are biomarkers in various cancer types . Multiple studies have identified piRNA in exosomes from human saliva  and plasma . Although our study has identified small RNA species that are candidate biomarkers for NSCLC, further research is required to determine the specific sequences of these EV RNA species, and more importantly, to evaluate their functions and mechanisms in NSCLC pathogenesis.
Differently expressed miRNAs in NSCLC EVs
Our miRNA expression analysis showed that the majority of EV miRNAs from the CTRL group have a lower expression than those from ADC and SQCC groups, and there is less diversity of EV miRNAs between ADC and CTRL than between SQCC and CTRL. The EV miRNA profiles also showed that the majority of miRNAs were not significantly differentially expressed, and about 10% of miRNAs were significantly upregulated in NSCLC plasma EVs, while only about 3% of miRNAs were significantly downregulated in NSCLC plasma EVs.
Our sequencing results also indicate that there is bigger difference of EV miRNA between SQCC and CTRL groups than between ADC and CTRL groups, which may represent a distinction between ADC and SQCC. Multiple studies have shown many differences between ADC and SQCC clinical and pathological characteristics , as well as genomic alteration  and gene expression profiles . For example, recent studies on exome sequencing and copy number profiles of 1144 lung cancers demonstrated that somatic gene alterations, including gene mutations and recurrent somatic copy number alterations, are largely distinct in lung ADC and SQCC . It has also been widely reported that miRNA expression profiles are very different between lung ADC and SQCC tumors . Furthermore, recent studies on tumor-derived exosomal miRNAs have also demonstrated the exosomal miRNA expression variety between lung ADC and SQCC, and identified ADC- or SQCC-specific miRNAs for the diagnosis of NSCLC subtypes . Overall, our EV miRNA profiles illustrate differences among lung ADC, SQCC, and healthy individuals, including distinctions between lung ADC and SQCC, consistent with previous reports.
Our qRT-PCR results showed that expression of miR-20a-5p, miR-20b-5p, miR-30b-5p, miR-451a, and miR-665 are significantly upregulated in NSCLC plasma EVs, consistent with our RNA-Seq results. The recent reports on deep sequencing of tumor-derived exosomal miRNA from NSCLC (16 lung ADC and 10 lung SQCC) clinical stage I patients have shown that miR-451a is significantly upregulated and miR-665 is significantly downregulated in lung SQCC exosomes, and expression of miR-20a-5p, miR-20b-5p, and miR-30b-5p is not significantly different between ADC and SQCC . Our findings regarding the expression pattern of miR-20a-5p, miR-20b-5p, miR-30b-5p, and miR-665 differ from these previous reports, perhaps because samples in our study were from patients of various clinical stages, while the previous studies contained only samples from stage I patients. In fact, other studies have shown that miRNA expression profiles from different clinical stages are distinct [65, 66]; for example, expressions of miR-21 and miR-4257 are significantly upregulated in plasma exosomes from NSCLC patients of more advanced stage .
Selective sorting of small RNAs into NSCLC EVs
Although most small RNAs are released into extracellular circulation in a non-selective manner , some specific small RNAs are selectively sorted into exosomes and other EVs [42, 69]. There are two possible routes for miRNA export via exosomes, an Ago-associated pathway and a RNA-binding protein-dependent pathway. For example, some proteins, such as SYNCRIP and hnRNPA2B1 , can guide the sorting of certain miRNAs into exosomes by binding to specific motifs. Other proteins, such as the RISC (RNA induced silence complex [RISC]) protein AGO2, may control incorporation of exosomal miRNAs in a miRNA sequence-independent manner . Sheckman and his co-workers found that RNA-binding protein YBX1 binds to miR-223 and is required for its selective sorting into exosomes, however they failed to discover the recognition motif sequences, suggesting the possibility that the recognition motif for sorting into exosomes may be based on secondary rather than primary RNA structures . Of interest, we found that the expression of miR-451a and miR-122-5p were significantly downregulated in NSCLC cell lysates, but significantly upregulated EVs of the same cells, indicating that miR-451a and miR-122-5p are likely to be selectively sorted into NSCLC EVs.
Both miR-122-5p [73–75] and miR-451a [76–79] are tumor suppressor miRNAs in various tumors ; in lung cancer, miR-122-5p inhibits metastasis and epithelial mesenchymal transition of NSCLC cells , while miR-451a inhibits proliferation of NSCLC cells, enhances NSCLC cell apoptosis, and suppresses lung tumor growth . Similarly, it was reported that major vault protein mediates selective sorting of the tumor suppressor miR-193a into colon cancer cell exosomes, which paradoxically promotes cancer cell proliferation and tumor progression . Therefore, the sorting of tumor suppressor miR-122-5p and miR-451a by NSCLC cells into EVs may indirectly play a role in promoting NSCLC progression.
Additionally, another recent report demonstrated that exosomal delivery of miR-122-5p suppresses glucose uptake by niche cells by downregulating the glycolytic enzyme pyruvate kinase, thus promoting tumor metastasis by increasing nutrient availability in the premetastatic niche . And miR-451a has also been reported as a negative regulator of glucose production and glucose homeostasis by targeting glycerol kinase-mediated gluconeogenesis in normal hepatic cells . Our KEGG pathway analysis showed that predicted targets of miR-122-5p and miR-451a are significantly enriched for the endocytic pathway, that is, miR-122-5p and miR-451a are very likely to regulate the endocytic pathway. Thus, NSCLC EVs delivery of miR-122-5p and miR-451a may also function to reprogram metabolism processes or interfere with endocytosis of a non-tumor cell, which conversely may provide a more appropriate environment for tumor progression and metastasis.
Overall, we have illustrated the different expression signatures of various small RNA species in EVs from lung ADC, SQCC, and healthy controls. We first identified EV YRNA fragments as a novel candidate class of circulating biomarkers for NSCLC diagnosis, and found that hY4-derived fragments functioned as tumor suppressors in NSCLC. Additionally, we demonstrated that miR-451a and miR-122-5p may be selectively sorted into NSCLC EVs. Although our studies were properly designed, including the use of an appropriate sample number and multi-aspect validation with clinical samples and cultured cell lines, some deficiencies and limitations exist. For instance, the robustness of our small RNA sequencing results was not high enough, and many differentially expressed miRNAs identified by RNA-Seq were not consistent with the qRT-PCR validation results. Importantly, our studies have indicated a wider range of biological effects in NSCLC that may be mediated by EV small RNAs than previously known. Nevertheless, further thorough investigations of the function and mechanism of these EV small RNAs in tumorigenesis are urgently needed.
The human lung adenocarcinoma cell lines (A549, H157, H1975, H2228, H3122, H2030, and H1299), and human large cell lung cancer cell line H460 were cultured in RPMI-1640 medium (GIBCO); the human squamous cell lung carcinoma cells (SK-MES-1) were cultured in MEM medium (GIBCO); and the human normal embryonic lung fibroblast cell line IMR-90 was cultured in MEM medium supplemented with 10% FBS (GIBCO), GlutaMAX (GIBCO), MEM non-essential amino acids solution (GIBCO), sodium pyruvate (GIBCO), and penicillin–streptomycin solution. All cells were maintained in a humidified chamber at 37 °C in 5% CO2.
Patient plasma samples
Plasma samples from lung ADC patients, SQCC patients, and healthy individuals were provided by Tongji Hospital and Hubei Cancer Hospital. The diagnosis of lung cancer subtypes was referred to 2015 World Health Organization classification of lung tumors. And healthy individuals without lung cancer or any other diseases were selected as the control group. Blood was collected and dispensed into an EDTA anticoagulation tube, then mixed gently to ensure exposure to the EDTA-coated walls. Plasma was separated by centrifugation at 5000 rpm for 5 min at 4 °C. The clear supernatant was transferred to a labeled tube, and stored at − 80 °C. This study was approved by ethics committees of Tongji Hospital and Hubei Cancer Hospital, and all methods were carried out in accordance with relevant guidelines and regulations. All donors gave permission and written informed consent was obtained from all participants.
Preparation of EV-free FBS
Sixty milliliters of fetal bovine serum (FBS; GIBCO) was added to 11 ml ultracentrifuge tubes with an adaptor, placed in a Type P40ST-2079 rotor (Hitachi), and spun in a Hitachi CP80WX ultracentrifuge at 100,000g at 4 °C for 16 h. The supernatant of each tube was then transferred to a 50 ml tube and stored at − 20 °C for future use.
EV isolation from plasma
Extracellular vesicles from plasma were isolated using Ribo™ Exosome Isolation Reagent (CAT. NO. C10110-2, Ribobio, China) according to the manufacturer’s protocol. First, plasma samples were centrifuged at 2000g at room temperature for 20 min to remove residual cells and debris, and then the supernatant was transferred to a new tube and centrifuged at 10,000g at room temperature for 20 min. The supernatant was transferred to a new tube and mixed adequately with 1/3 volume isolation reagent, and the mixture was placed at 4 °C for 30 min and centrifuged at 15,000g at 4 °C for 2 min. Finally, the supernatant was removed and the EV pellet was recovered by re-suspending in phosphate-buffered saline (PBS).
EV isolation from cell culture medium
After incubating cells with medium containing EV-free FBS for 48 h, EVs from cell culture media were isolated using Ribo™ Exosome Isolation Reagent (CAT. NO. C10130-2, Ribobio, China) according to the manufacturer’s protocol. First, cell culture medium was centrifuged at 2000g at room temperature for 30 min to remove residual cells and debris, and then the supernatant was transferred to a new tube and mixed with 1/3 volume isolation reagent. Next, the mixture was placed at 4 °C overnight, and then centrifuged at 1500g at 4 °C for 30 min. Finally, the supernatant was removed and the EV pellet was recovered by re-suspending in PBS.
Western blot analysis
Extracellular vesicle protein was prepared by re-suspending EV pellets in RIPA buffer supplemented with protease inhibitor (Roche). Twenty micrograms of protein from each sample was separated on a 10% SDS-PAGE gel in parallel with a protein marker (Thermo). Proteins were transferred to a polyvinylidene fluoride (PVDF) membrane, and blocked with 5% skim milk dissolved in Tween/Tris buffered saline (TTBS) for 1 h at room temperature. Membranes were then incubated with primary antibodies against TSG101 (ProteinTech) or CD9 (ProteinTech) at a dilution of 1:5000 with 3% skim milk diluted in TTBS for 1 h at room temperature. After washing 6 times with TTBS, membranes were incubated with horseradish peroxidase conjugated goat anti-rabbit IgG secondary antibodies (Santa Cruz Biotechnology) at a dilution of 1:10,000 in 3% skim milk dissolved in TTBS for 1 h at room temperature. After washing 6 times with TTBS, membranes were covered with ECL western blot substrate solution and visualized by exposing to film and developing in a film processor.
Transmission electron microscopy (TEM)
The EV pellets were re-suspended in 100 µl of 2% paraformaldehyde, and 5 µl was placed on Formvar–carbon-coated electron microscopy grids and adsorbed by the membrane for 20 min in a dry environment. After washing the samples with several drops of PBS on a sheet of Parafilm, samples were incubated on drops of buffered 1% glutaraldehyde for 5 min and then washed several times on drops of distilled water. Next, samples were negatively stained on drops of 4% uranyl acetate for 10 min on ice. Then, the grids were removed and excess fluid was blotted with filter paper, so that a thin film was left behind over the EV side of the grid. The samples were allowed to dry, and observed using a JEM-1400plus transmission electron microscope (JEOL) at an accelerating voltage of 80 kV.
Flow cytometry analysis of exosomal marker proteins
Isolated EV pellets were diluted in PBS, then incubated with FITC-conjugated anti-CD63 or anti-CD81 antibodies (Abcam, US) diluted in PBS/0.5% BSA for 1 h at 37 °C. The EV samples were analyzed with a Cytoflex flow cytometer (Beckman, US), and IgG incubated samples were used as a negative control.
Nanoparticle tracking analysis
Isolated EV pellets were diluted in PBS and analyzed using the Nanosight NS300 system (Malvern Instruments, UK). The particle size distribution was calculated using Nanosight Tracking Analysis software. Size distribution profiles were averaged across three replicates or each sample to derive the representative size distribution profiles.
To assess cell proliferation, cells were transiently transfected with miRNAs mimics (Ribobio, China) for 24 h using lipofectamine 2000 (Thermo, US), and then seeded in 96-well culture plates at a concentration of 6000 cells per well, with each group comprised of six replicates. The cells were incubated with 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide (MTT) (Sigma, US) for 4 h at 37 °C, and then the supernatant was removed and dimethyl sulfoxide was added. Absorbance was measured at 570 nm with microplate reader (BioTek, US), and above tests were performed every 24 h for consecutive days.
Colony formation assay
After transfected for 24 h, cells are seeded into 6-well culture plates in dilutions of 200 cells per well to form colonies in 9 days. Colonies are fixed with 4% (v/v) glutaraldehyde for 15 min, stained with 1% (w/v) crystal violet for 20 min, and counted using a stereomicroscope.
Isolation of total RNA from cells or EVs was performed with Trizol reagent (Invitrogen, Life Technologies) and the Direct-zol™ RNA MiniPrep kit (Zymo Research Corp.) according to the manufacturer’s protocol. RNA concentration was quantitated using a NanoDrop 2000 Spectrophotometer (Thermo Scientific).
Quantitative real-time reverse transcription PCR (qRT-PCR)
Spike-in control cel-miR-39 mimic (miRNeasy Serum/Plasma Spike-In Control, Qiagen) was added to total RNA samples as internal control, and the mixed RNA was reverse-transcribed to cDNA using a reverse transcriptase kit (miScript II RT Kit, Qiagen). SYBR-based qRT-PCR was performed on an ABI7500 real-time PCR amplifier (Applied Biosystems) using specific primers (miScript primer Assay, Qiagen) and a universal PCR kit (miScript SYBR Green, Qiagen) according to the manufacturer’s protocol. For hY4 RNA fragments detection, U6 snRNA was used as internal control. Amplifying conditions were as follows: 95 °C for 15 min, followed by 40 cycles of 94 °C for 15 s, 55 °C for 30 s, and 70 °C for 34 s [84, 85]. Gene expression data were normalized to internal control expression, and the relative expression was determined as 2−Δ Ct, where Δ Ct = Ct (target gene) − Ct (internal control).
cDNA library generation and high-throughput sequencing
Sequencing libraries were generated as follows: following isolation of total RNA from plasma EVs, RNA samples were analyzed with an Agilent 2100 Bioanalyzer (Agilent) using a total RNA nanochip, then the small RNA fractions (< 60 nt) were separated by electrophoresis. Next, small RNA fragments were ligated to 5′-adapter and 3′-adapters, and then reverse-transcribed and PCR amplified, with indexes included in the PCR reaction. PCR products were run on a 6% polyacrylamide gel and selected to include RNAs with an insert size < 51 nt. The DNA-containing index was eluted from the gel, precipitated, and dissolved in nuclease-free water. Small RNA inserts were sequenced using a HiSeq 2500 high-throughput sequencing system (Illumina).
Sequencing data processing and analysis
To get clean reads, the adaptor sequences, contaminated reads, low quality reads, less than 17 nt reads, and poly-A/T/C/G/N repeat sequences were removed. Next, the clean reads (17–45 nt) were mapped to the human genome database using Burrows–Wheeler Aligner software . Uniquely aligned reads were categorized and annotated to known non-coding small RNA transcripts using Rfam version 11.0  for YRNA, rRNA, snRNA, snoRNA, and tRNA. miRBase version 21  was used for miRNA, and piRNABank  was used for piRNA.
Next, the quantity of small RNA in each sample was obtained and normalized as RPM (number of mapped reads per million clean reads). The clean reads of miRNA, containing the entire reference sequence from the miRBase database and a length difference within 4 nt, were confirmed, as well as miRNA generated from the same end of the same precursor and a length discrepancy below 3 nt. Considering sequencing error, reads with less than ten copies were removed. The sequence length distribution and percentage of different small RNA species were also analyzed with the statistical data of reads. The distribution of small RNA across the whole human genome was performed using a R language package RCircos . The cumulative frequency distribution was performed using GraphPad Prism 5, and the relative frequency is shown as a fraction. The correlation analysis was also determined using GraphPad Prism 5 with a Gaussian distribution.
Significant differences in miRNA expression were determined using R package edgeR , with thresholds of a P value < 0.05, a false discovery rate (FDR) < 0.05, and a |log2 (Fold Change)| ≥ 1. Scatter plots and MA plots of miRNA expression differences were generated using R package edgeR. The heatmap of miRNA expression differences was plotted using function heatmap.2; in R package gplots, and miRNAs with extremely low abundance (less than 1 RPM) in all groups were removed. Considering that RNA-Seq produces a wide range of read counts per gene, and genes with a low coverage of reads can produce artificially high fold-change values , RPM was added with a moderate cutoff of 0.1 and transformed as log2 for drawing heatmaps.
Gene ontology (GO) and KEGG enrichment analysis of predicted miRNA targets
miRNA targets were predicted with default parameters using starBase version 2.0, a database of predicted miRNA targets gathered from targetScan, picTar, RNA22, PITA, and miRanda databases . The GO enrichment analysis of miRNA targets was performed with the BinGO plugin  for Cytoscape software  version 3.4.0 using whole annotation as a reference set; hypergeometric test, false discovery rate (FDR) correction, P < 0.05 significance level and GO BP (Biological Process) ontology files were selected. For KEGG pathway enrichment analysis, the DAVID online tool (https://david.ncifcrf.gov)  was used with Fisher’s exact test, FDR correction, and P < 0.05 significance level.
MG, TW and GS designed the study. CL and FQ performed experiments. TW, FH, HX, CL and GH collected patient plasma samples. CL and MG performed bioinformatics analyses. MG and CL analyzed and interpreted the data, prepared the manuscript, wrote the manuscript with contributions from all authors. All authors read and approved the final manuscript.
This work was supported by National Natural Science Foundation of China Grants (31271511 to M. G., 81470252 to T. W., 31370187 and 81572447 to G. S.). The funders had no role in study design, collection, analysis and interpretation of data, in the writing of the report and in the decision to submit this article for publication.
The authors declare that they have no competing interests.
High through-put sequencing data presented in this article have been deposited in the NCBI Gene Expression Omnibus and are accessible under the Accession Number GSE103149.
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided 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 Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Maas SL, Breakefield XO, Weaver AM. Extracellular vesicles: unique intercellular delivery vehicles. Trends Cell Biol. 2017;27(3):172–88.PubMedView ArticleGoogle Scholar
- Vlassov AV, Magdaleno S, Setterquist R, Conrad R. Exosomes: current knowledge of their composition, biological functions, and diagnostic and therapeutic potentials. Biochem Biophys Acta. 2012;1820(7):940–8.PubMedView ArticleGoogle Scholar
- Thery C, Zitvogel L, Amigorena S. Exosomes: composition, biogenesis and function. Nat Rev Immunol. 2002;2(8):569–79.PubMedView ArticleGoogle Scholar
- Munich S, Sobo-Vujanovic A, Buchser WJ, Beer-Stolz D, Vujanovic NL. Dendritic cell exosomes directly kill tumor cells and activate natural killer cells via TNF superfamily ligands. Oncoimmunology. 2012;1(7):1074–83.PubMedPubMed CentralView ArticleGoogle Scholar
- Mulcahy LA, Pink RC, Carter DR. Routes and mechanisms of extracellular vesicle uptake. J Extracell Vesicles. 2014;3:24641.View ArticleGoogle Scholar
- Tkach M, Thery C. Communication by extracellular vesicles: where we are and where we need to go. Cell. 2016;164(6):1226–32.PubMedView ArticleGoogle Scholar
- Storz G. An expanding universe of noncoding RNAs. Science. 2002;296(5571):1260–3.PubMedView ArticleGoogle Scholar
- Zhang J, Li S, Li L, Li M, Guo C, Yao J, Mi S. Exosome and exosomal microRNA: trafficking, sorting, and function. Genom Proteom Bioinform. 2015;13(1):17–24.View ArticleGoogle Scholar
- Hannafon BN, Ding WQ. Intercellular communication by exosome-derived microRNAs in cancer. Int J Mol Sci. 2013;14(7):14240–69.PubMedPubMed CentralView ArticleGoogle Scholar
- Teng Y, Ren Y, Hu X, Mu J, Samykutty A, Zhuang X, Deng Z, Kumar A, Zhang L, Merchant ML, et al. MVP-mediated exosomal sorting of miR-193a promotes colon cancer progression. Nat Commun. 2017;8:14448.PubMedPubMed CentralView ArticleGoogle Scholar
- Ostenfeld MS, Jeppesen DK, Laurberg JR, Boysen AT, Bramsen JB, Primdal-Bengtson B, Hendrix A, Lamy P, Dagnaes-Hansen F, Rasmussen MH, et al. Cellular disposal of miR23b by RAB27-dependent exosome release is linked to acquisition of metastatic properties. Can Res. 2014;74(20):5758–71.View ArticleGoogle Scholar
- Yang H, Fu H, Xu W, Zhang X. Exosomal non-coding RNAs: a promising cancer biomarker. Clin Chem Lab Med. 2016;54(12):1871–9.PubMedView ArticleGoogle Scholar
- Thind A, Wilson C. Exosomal miRNAs as cancer biomarkers and therapeutic targets. J Extracell Vesicles. 2016;5:31292.PubMedView ArticleGoogle Scholar
- Manterola L, Guruceaga E, Perez-Larraya JG, Gonzalez-Huarriz M, Jauregui P, Tejada S, Diez-Valle R, Segura V, Sampron N, Barrena C, et al. A small noncoding RNA signature found in exosomes of GBM patient serum as a diagnostic tool. Neuro-oncology. 2014;16(4):520–7.PubMedPubMed CentralView ArticleGoogle Scholar
- Hannafon BN, Trigoso YD, Calloway CL, Zhao YD, Lum DH, Welm AL, Zhao ZJ, Blick KE, Dooley WC, Ding WQ. Plasma exosome microRNAs are indicative of breast cancer. Breast Cancer Res BCR. 2016;18(1):90.PubMedView ArticleGoogle Scholar
- Wang H, Hou L, Li A, Duan Y, Gao H, Song X. Expression of serum exosomal microRNA-21 in human hepatocellular carcinoma. Biomed Res Int. 2014;2014:864894.PubMedPubMed CentralGoogle Scholar
- Tanaka Y, Kamohara H, Kinoshita K, Kurashige J, Ishimoto T, Iwatsuki M, Watanabe M, Baba H. Clinical impact of serum exosomal microRNA-21 as a clinical biomarker in human esophageal squamous cell carcinoma. Cancer. 2013;119(6):1159–67.PubMedView ArticleGoogle Scholar
- Lai X, Wang M, McElyea SD, Sherman S, House M, Korc M. A microRNA signature in circulating exosomes is superior to exosomal glypican-1 levels for diagnosing pancreatic cancer. Cancer Lett. 2017;393:86–93.PubMedPubMed CentralView ArticleGoogle Scholar
- Siegel RL, Miller KD, Jemal A. Cancer statistics, 2017. CA Cancer J Clin. 2017;67(1):7–30.PubMedView ArticleGoogle Scholar
- Wakelee H, Kelly K, Edelman MJ. 50 years of progress in the systemic therapy of non-small cell lung cancer. In: American Society of Clinical Oncology educational book American Society of Clinical Oncology Meeting. 2014. p. 177–89.Google Scholar
- Alamgeer M, Peacock CD, Matsui W, Ganju V, Watkins DN. Cancer stem cells in lung cancer: evidence and controversies. Respirology. 2013;18(5):757–64.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, Jemal A, Yu XQ, He J. Cancer statistics in China, 2015. CA Cancer J Clin. 2016;66(2):115–32.PubMedView ArticleGoogle Scholar
- Heuvers ME, Wisnivesky J, Stricker BH, Aerts JG. Generalizability of results from the national lung screening trial. Eur J Epidemiol. 2012;27(9):669–72.PubMedView ArticleGoogle Scholar
- Kawase A, Yoshida J, Ishii G, Nakao M, Aokage K, Hishida T, Nishimura M, Nagai K. Differences between squamous cell carcinoma and adenocarcinoma of the lung: are adenocarcinoma and squamous cell carcinoma prognostically equal? Jpn J Clin Oncol. 2012;42(3):189–95.PubMedView ArticleGoogle Scholar
- Lu C, Chen H, Shan Z, Yang L. Identification of differentially expressed genes between lung adenocarcinoma and lung squamous cell carcinoma by gene expression profiling. Mol Med Rep. 2016;14(2):1483–90.PubMedPubMed CentralView ArticleGoogle Scholar
- Campbell JD, Alexandrov A, Kim J, Wala J, Berger AH, Pedamallu CS, Shukla SA, Guo G, Brooks AN, Murray BA, et al. Distinct patterns of somatic genome alterations in lung adenocarcinomas and squamous cell carcinomas. Nat Genet. 2016;48(6):607–16.PubMedPubMed CentralView ArticleGoogle Scholar
- Munagala R, Aqil F, Gupta RC. Exosomal miRNAs as biomarkers of recurrent lung cancer. Tumour Biol. 2016;37(8):10703–14.PubMedView ArticleGoogle Scholar
- Cazzoli R, Buttitta F, Di Nicola M, Malatesta S, Marchetti A, Rom WN, Pass HI. microRNAs derived from circulating exosomes as noninvasive biomarkers for screening and diagnosing lung cancer. J Thorac Oncol. 2013;8(9):1156–62.PubMedPubMed CentralView ArticleGoogle Scholar
- Jin X, Chen Y, Chen H, Fei S, Chen D, Cai X, Liu L, Lin B, Su H, Zhao L, et al. Evaluation of tumor-derived exosomal miRNA as potential diagnostic biomarkers for early stage non-small-cell lung cancer using next-generation sequencing. Clin Cancer Res. 2017;23:5311–9.PubMedView ArticleGoogle Scholar
- Kamerkar S, LeBleu VS, Sugimoto H, Yang S, Ruivo CF, Melo SA, Lee JJ, Kalluri R. Exosomes facilitate therapeutic targeting of oncogenic KRAS in pancreatic cancer. Nature. 2017;546(7659):498–503.PubMedPubMed CentralGoogle Scholar
- Melo SA, Luecke LB, Kahlert C, Fernandez AF, Gammon ST, Kaye J, LeBleu VS, Mittendorf EA, Weitz J, Rahbari N, et al. Glypican-1 identifies cancer exosomes and detects early pancreatic cancer. Nature. 2015;523(7559):177–82.PubMedPubMed CentralView ArticleGoogle Scholar
- Huang X, Yuan T, Tschannen M, Sun Z, Jacob H, Du M, Liang M, Dittmar RL, Liu Y, Liang M, et al. Characterization of human plasma-derived exosomal RNAs by deep sequencing. BMC Genom. 2013;14:319.View ArticleGoogle Scholar
- Lasser C, Eldh M, Lotvall J. Isolation and characterization of RNA-containing exosomes. J Vis Exp JoVE. 2012;59:e3037.Google Scholar
- Maraia RJ, Sasaki-Tozawa N, Driscoll CT, Green ED, Darlington GJ. The human Y4 small cytoplasmic RNA gene is controlled by upstream elements and resides on chromosome 7 with all other hY scRNA genes. Nucleic Acids Res. 1994;22(15):3045–52.PubMedPubMed CentralView ArticleGoogle Scholar
- Perreault J, Noel JF, Briere F, Cousineau B, Lucier JF, Perreault JP, Boire G. Retropseudogenes derived from the human Ro/SS-A autoantigen-associated hY RNAs. Nucleic Acids Res. 2005;33(6):2032–41.PubMedPubMed CentralView ArticleGoogle Scholar
- Dhahbi JM, Spindler SR, Atamna H, Boffelli D, Mote P, Martin DI. 5′-YRNA fragments derived by processing of transcripts from specific YRNA genes and pseudogenes are abundant in human serum and plasma. Physiol Genom. 2013;45(21):990–8.View ArticleGoogle Scholar
- Ferguson SW, Nguyen J. Exosomes as therapeutics: the implications of molecular composition and exosomal heterogeneity. J Control Release. 2016;228:179–90.PubMedView ArticleGoogle Scholar
- Errico A. Lung cancer: heterogeneity in space and time. Nat Rev Clin Oncol. 2014;11(12):684.Google Scholar
- de Bruin EC, McGranahan N, Mitter R, Salm M, Wedge DC, Yates L, Jamal-Hanjani M, Shafi S, Murugaesu N, Rowan AJ, et al. Spatial and temporal diversity in genomic instability processes defines lung cancer evolution. Science. 2014;346(6206):251–6.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang R, Wang ZX, Yang JS, Pan X, De W, Chen LB. MicroRNA-451 functions as a tumor suppressor in human non-small cell lung cancer by targeting ras-related protein 14 (RAB14). Oncogene. 2011;30(23):2644–58.PubMedView ArticleGoogle Scholar
- Ma D, Jia H, Qin M, Dai W, Wang T, Liang E, Dong G, Wang Z, Zhang Z, Feng F. MiR-122 induces radiosensitization in non-small cell lung cancer cell line. Int J Mol Sci. 2015;16(9):22137–50.PubMedPubMed CentralView ArticleGoogle Scholar
- Nolte-’t Hoen EN, Buermans HP, Waasdorp M, Stoorvogel W, Wauben MH, t Hoen PA. Deep sequencing of RNA from immune cell-derived vesicles uncovers the selective incorporation of small non-coding RNA biotypes with potential regulatory functions. Nucleic Acids Res. 2012;40(18):9272–85.PubMedPubMed CentralView ArticleGoogle Scholar
- Hazan-Halevy I, Rosenblum D, Weinstein S, Bairey O, Raanani P, Peer D. Cell-specific uptake of mantle cell lymphoma-derived exosomes by malignant and non-malignant B-lymphocytes. Cancer Lett. 2015;364(1):59–69.PubMedPubMed CentralView ArticleGoogle Scholar
- Valadi H, Ekstrom K, Bossios A, Sjostrand M, Lee JJ, Lotvall JO. Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat Cell Biol. 2007;9(6):654–9.PubMedView ArticleGoogle Scholar
- Tosar JP, Gambaro F, Sanguinetti J, Bonilla B, Witwer KW, Cayota A. Assessment of small RNA sorting into different extracellular fractions revealed by high-throughput sequencing of breast cell lines. Nucleic Acids Res. 2015;43(11):5601–16.PubMedPubMed CentralView ArticleGoogle Scholar
- Lerner MR, Boyle JA, Hardin JA, Steitz JA. Two novel classes of small ribonucleoproteins detected by antibodies associated with lupus erythematosus. Science. 1981;211(4480):400–2.PubMedView ArticleGoogle Scholar
- Christov CP, Trivier E, Krude T. Noncoding human Y RNAs are overexpressed in tumours and required for cell proliferation. Br J Cancer. 2008;98(5):981–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Dhahbi JM, Spindler SR, Atamna H, Boffelli D, Martin DI. Deep sequencing of serum small RNAs identifies patterns of 5′ tRNA half and YRNA fragment expression associated with breast cancer. Biomark Cancer. 2014;6:37–47.PubMedPubMed CentralView ArticleGoogle Scholar
- Kowalski MP, Krude T. Functional roles of non-coding Y RNAs. Int J Biochem Cell Biol. 2015;66:20–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Gebetsberger J, Polacek N. Slicing tRNAs to boost functional ncRNA diversity. RNA Biol. 2013;10(12):1798–806.PubMedPubMed CentralView ArticleGoogle Scholar
- Guzman N, Agarwal K, Asthagiri D, Yu L, Saji M, Ringel MD, Paulaitis ME. Breast cancer-specific miR signature unique to extracellular vesicles includes “microRNA-like” tRNA fragments. Mol Cancer Res MCR. 2015;13(5):891–901.PubMedView ArticleGoogle Scholar
- Olvedy M, Scaravilli M, Hoogstrate Y, Visakorpi T, Jenster G, Martens-Uzunova ES. A comprehensive repertoire of tRNA-derived fragments in prostate cancer. Oncotarget. 2016;7(17):24766–77.PubMedPubMed CentralView ArticleGoogle Scholar
- Mleczko AM, Celichowski P, Bakowska-Zywicka K. Ex-translational function of tRNAs and their fragments in cancer. Acta Biochim Pol. 2014;61(2):211–6.PubMedGoogle Scholar
- Vojtech L, Woo S, Hughes S, Levy C, Ballweber L, Sauteraud RP, Strobl J, Westerberg K, Gottardo R, Tewari M, et al. Exosomes in human semen carry a distinctive repertoire of small non-coding RNAs with potential regulatory functions. Nucleic Acids Res. 2014;42(11):7290–304.PubMedPubMed CentralView ArticleGoogle Scholar
- Shao Y, Sun Q, Liu X, Wang P, Wu R, Ma Z. tRF-Leu-CAG promotes cell proliferation and cell cycle in non-small cell lung cancer. Chem Biol Drug Des. 2017;90:730–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Matera AG, Terns RM, Terns MP. Non-coding RNAs: lessons from the small nuclear and small nucleolar RNAs. Nat Rev Mol Cell Biol. 2007;8(3):209–20.PubMedView ArticleGoogle Scholar
- Kuhlmann JD, Baraniskin A, Hahn SA, Mosel F, Bredemeier M, Wimberger P, Kimmig R, Kasimir-Bauer S. Circulating U2 small nuclear RNA fragments as a novel diagnostic tool for patients with epithelial ovarian cancer. Clin Chem. 2014;60(1):206–13.PubMedView ArticleGoogle Scholar
- Baraniskin A, Nopel-Dunnebacke S, Ahrens M, Jensen SG, Zollner H, Maghnouj A, Wos A, Mayerle J, Munding J, Kost D, et al. Circulating U2 small nuclear RNA fragments as a novel diagnostic biomarker for pancreatic and colorectal adenocarcinoma. Int J Cancer. 2013;132(2):E48–57.PubMedView ArticleGoogle Scholar
- Seto AG, Kingston RE, Lau NC. The coming of age for piwi proteins. Mol Cell. 2007;26(5):603–9.PubMedView ArticleGoogle Scholar
- Weick EM, Miska EA. piRNAs: from biogenesis to function. Development. 2014;141(18):3458–71.PubMedView ArticleGoogle Scholar
- Brennecke J, Malone CD, Aravin AA, Sachidanandam R, Stark A, Hannon GJ. An epigenetic role for maternally inherited piRNAs in transposon silencing. Science. 2008;322(5906):1387–92.PubMedPubMed CentralView ArticleGoogle Scholar
- Assumpcao CB, Calcagno DQ, Araujo TM, Santos SE, Santos AK, Riggins GJ, Burbano RR, Assumpcao PP. The role of piRNA and its potential clinical implications in cancer. Epigenomics. 2015;7(6):975–84.PubMedPubMed CentralView ArticleGoogle Scholar
- Ogawa Y, Taketomi Y, Murakami M, Tsujimoto M, Yanoshita R. Small RNA transcriptomes of two types of exosomes in human whole saliva determined by next generation sequencing. Biol Pharm Bull. 2013;36(1):66–75.PubMedView ArticleGoogle Scholar
- Hamamoto J, Soejima K, Yoda S, Naoki K, Nakayama S, Satomi R, Terai H, Ikemura S, Sato T, Yasuda H, et al. Identification of microRNAs differentially expressed between lung squamous cell carcinoma and lung adenocarcinoma. Mol Med Rep. 2013;8(2):456–62.PubMedView ArticleGoogle Scholar
- Slattery ML, Herrick JS, Mullany LE, Valeri N, Stevens J, Caan BJ, Samowitz W, Wolff RK. An evaluation and replication of miRNAs with disease stage and colorectal cancer-specific mortality. Int J Cancer. 2015;137(2):428–38.PubMedView ArticleGoogle Scholar
- Hart M, Nolte E, Wach S, Szczyrba J, Taubert H, Rau TT, Hartmann A, Grasser FA, Wullich B. Comparative microRNA profiling of prostate carcinomas with increasing tumor stage by deep sequencing. Mol Cancer Res MCR. 2014;12(2):250–63.PubMedView ArticleGoogle Scholar
- Dejima H, Iinuma H, Kanaoka R, Matsutani N, Kawamura M. Exosomal microRNA in plasma as a non-invasive biomarker for the recurrence of non-small cell lung cancer. Oncol Lett. 2017;13(3):1256–63.PubMedPubMed CentralView ArticleGoogle Scholar
- Turchinovich A, Weiz L, Langheinz A, Burwinkel B. Characterization of extracellular circulating microRNA. Nucleic Acids Res. 2011;39(16):7223–33.PubMedPubMed CentralView ArticleGoogle Scholar
- van Balkom BW, Eisele AS, Pegtel DM, Bervoets S, Verhaar MC. Quantitative and qualitative analysis of small RNAs in human endothelial cells and exosomes provides insights into localized RNA processing, degradation and sorting. J Extracell Vesicles. 2015;4:26760.PubMedView ArticleGoogle Scholar
- Villarroya-Beltri C, Gutierrez-Vazquez C, Sanchez-Cabo F, Perez-Hernandez D, Vazquez J, Martin-Cofreces N, Martinez-Herrera DJ, Pascual-Montano A, Mittelbrunn M, Sanchez-Madrid F. Sumoylated hnRNPA2B1 controls the sorting of miRNAs into exosomes through binding to specific motifs. Nat Commun. 2013;4:2980.PubMedPubMed CentralView ArticleGoogle Scholar
- Guduric-Fuchs J, O’Connor A, Camp B, O’Neill CL, Medina RJ, Simpson DA. Selective extracellular vesicle-mediated export of an overlapping set of microRNAs from multiple cell types. BMC Genom. 2012;13:357.View ArticleGoogle Scholar
- Shurtleff MJ, Temoche-Diaz MM, Karfilis KV, Ri S, Schekman R. Y-box protein 1 is required to sort microRNAs into exosomes in cells and in a cell-free reaction. Elife. 2016;5:e19276.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang Y, Xing QF, Liu XQ, Guo ZJ, Li CY, Sun G. MiR-122 targets VEGFC in bladder cancer to inhibit tumor growth and angiogenesis. Am J Transl Res. 2016;8(7):3056–66.PubMedPubMed CentralGoogle Scholar
- Lu W, Zhang Y, Zhou L, Wang X, Mu J, Jiang L, Hu Y, Dong P, Liu Y. miR-122 inhibits cancer cell malignancy by targeting PKM2 in gallbladder carcinoma. Tumour Biol. 2015;37:15615–25.View ArticleGoogle Scholar
- Hsu SH, Wang B, Kota J, Yu J, Costinean S, Kutay H, Yu L, Bai S, La Perle K, Chivukula RR, et al. Essential metabolic, anti-inflammatory, and anti-tumorigenic functions of miR-122 in liver. J Clin Investig. 2012;122(8):2871–83.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang J, Zhao X, Shi J, Pan Y, Chen Q, Leng P, Wang Y. miR-451 suppresses bladder cancer cell migration and invasion via directly targeting c-Myc. Oncol Rep. 2016;36(4):2049–58.PubMedView ArticleGoogle Scholar
- Zang WQ, Yang X, Wang T, Wang YY, Du YW, Chen XN, Li M, Zhao GQ. MiR-451 inhibits proliferation of esophageal carcinoma cell line EC9706 by targeting CDKN2D and MAP3K1. World J Gastroenterol. 2015;21(19):5867–76.PubMedPubMed CentralView ArticleGoogle Scholar
- Li Y, Wang J, Dai X, Zhou Z, Liu J, Zhang Y, Li Y, Hou Y, Pang L, Wang X, et al. miR-451 regulates FoxO3 nuclear accumulation through Ywhaz in human colorectal cancer. Am J Transl Res. 2015;7(12):2775–85.PubMedPubMed CentralGoogle Scholar
- Huang JY, Zhang K, Chen DQ, Chen J, Feng B, Song H, Chen Y, Zhu Z, Lu L, De W, et al. MicroRNA-451: epithelial-mesenchymal transition inhibitor and prognostic biomarker of hepatocellular carcinoma. Oncotarget. 2015;6(21):18613–30.PubMedPubMed CentralView ArticleGoogle Scholar
- Pan X, Wang R, Wang ZX. The potential role of miR-451 in cancer diagnosis, prognosis, and therapy. Mol Cancer Ther. 2013;12(7):1153–62.PubMedView ArticleGoogle Scholar
- Qin H, Sha J, Jiang C, Gao X, Qu L, Yan H, Xu T, Jiang Q, Gao H. miR-122 inhibits metastasis and epithelial-mesenchymal transition of non-small-cell lung cancer cells. OncoTargets Ther. 2015;8:3175–84.Google Scholar
- Fong MY, Zhou W, Liu L, Alontaga AY, Chandra M, Ashby J, Chow A, O’Connor ST, Li S, Chin AR, et al. Breast-cancer-secreted miR-122 reprograms glucose metabolism in premetastatic niche to promote metastasis. Nat Cell Biol. 2015;17(2):183–94.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhuo S, Yang M, Zhao Y, Chen X, Zhang F, Li N, Yao P, Zhu T, Mei H, Wang S, et al. MicroRNA-451 negatively regulates hepatic glucose production and glucose homeostasis by targeting glycerol kinase-mediated gluconeogenesis. Diabetes. 2016;65(11):3276–88.PubMedView ArticleGoogle Scholar
- Zeng Y, Liu H, Kang K, Wang Z, Hui G, Zhang X, Zhong J, Peng W, Ramchandran R, Raj JU, et al. Hypoxia inducible factor-1 mediates expression of miR-322: potential role in proliferation and migration of pulmonary arterial smooth muscle cells. Sci Rep. 2015;5:12098.PubMedPubMed CentralView ArticleGoogle Scholar
- Niu Y, Zhang L, Qiu H, Wu Y, Wang Z, Zai Y, Liu L, Qu J, Kang K, Gou D. An improved method for detecting circulating microRNAs with S-Poly(T) Plus real-time PCR. Sci Rep. 2015;5:15100.PubMedPubMed CentralView ArticleGoogle Scholar
- Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25(14):1754–60.PubMedPubMed CentralView ArticleGoogle Scholar
- Daub J, Eberhardt RY, Tate JG, Burge SW. Rfam: annotating families of non-coding RNA sequences. Methods Mol Biol. 2015;1269:349–63.PubMedView ArticleGoogle Scholar
- Kozomara A, Griffiths-Jones S. miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014;42(Database issue):D68–73.PubMedView ArticleGoogle Scholar
- Sai Lakshmi S, Agrawal S. piRNABank: a web resource on classified and clustered piwi-interacting RNAs. Nucleic acids Res. 2008;36(Database issue):D173–7.PubMedView ArticleGoogle Scholar
- Zhang H, Meltzer P, Davis S. RCircos: an R package for Circos 2D track plots. BMC Bioinform. 2013;14:244.View ArticleGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.PubMedView ArticleGoogle Scholar
- Yuan Y-C, Wu X, Warden C. Optimal calculation of RNA-Seq fold-change values. Int J Comput Bioinform Silico Model. 2013;2(6):285–92.Google Scholar
- Li JH, Liu S, Zhou H, Qu LH, Yang JH. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2014;42(Database issue):D92–7.PubMedView ArticleGoogle Scholar
- Maere S, Heymans K, Kuiper M. BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics. 2005;21(16):3448–9.PubMedView ArticleGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.PubMedPubMed CentralView ArticleGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.View ArticleGoogle Scholar