Skip to content

Advertisement

  • Letter to the Editor
  • Open Access

Differential distribution and enrichment of non-coding RNAs in exosomes from normal and Cancer-associated fibroblasts in colorectal cancer

  • 1, 2,
  • 3,
  • 4, 5,
  • 1,
  • 6,
  • 1, 7,
  • 8,
  • 9,
  • 10,
  • 11,
  • 11,
  • 10,
  • 12,
  • 11,
  • 1, 10Email author and
  • 1, 11Email authorView ORCID ID profile
Molecular Cancer201817:114

https://doi.org/10.1186/s12943-018-0863-4

  • Received: 24 May 2018
  • Accepted: 24 July 2018
  • Published:

Abstract

Exosome production from cancer-associated fibroblasts seems to be an important driver of tumor progression. We report the first in-depth biotype characterization of ncRNAs, analyzed by Next Generation Sequencing and Bioinformatics, expressed in established primary human normal and cancer-associated fibroblasts (CAFs) from cancer and normal mucosa tissues from 9 colorectal cancer patients, and/or packaged in their derived exosomes. Differential representation and enrichment analyses based on these ncRNAs revealed a significant number of differences between the ncRNA content of exosomes and the expression patterns of the normal and cancer-associated fibroblast cells. ncRNA regulatory elements are specifically packaged in CAF-derived exosomes, supporting a specific cross-talk between CAFs and colon cancer cells and/or other stromal cells, mediated by exosomes. These sncRNAs are potential biomarkers present in cancer-associated fibroblast-derived exosomes, which should thereby contribute to developing new non-invasive diagnostic, prognostic and predictive methods for clinical applications in management of cancer patients.

Keywords

  • Colon Cancer
  • Liquid biopsy
  • Tumor microenvironment
  • Exosomes
  • Non-coding RNAs
  • Next generation sequencing

Tumor progression is deeply influenced by the local microenvironment. Fibroblasts usually named as cancer-associated fibroblasts (CAFs), are one of the most abundant and active cell types of the tumor microenvironment (TME). CAFs seem to regulate many aspects of tumorogenesis involving interactions between the malignant cells and other cells of the TME [1].

Exosomes, released by cells are important mediators of intercellular communication. In the tumor context, exosomes released by cancer cells transmit signals to cancer cells and also to stromal cells generating an active TME which promotes tumor progression [2]. Exosomes are also released by both cancer cells and stromal cells not only into the cancer microenvironment, but also into circulation [2]. Significant amounts of miRNAs packaged into exosomes were detected in many types of liquid biopsy samples [2], and certain types of circulating miRNAs strongly correlates with progression of different cancer types [3], however, few information about others ncRNAs contained in exosomes is unknown.

In the past decade large-scale analyses have focused on the comprehensive identification of non-coding RNAs (ncRNAs) and of new ones such us the long ncRNAs (lncRNAs). lncRNAs are poorly conserved and regulate gene expression by diverse mechanisms that are not yet fully understood. Recent data suggest that they may serve as master drivers of carcinogenesis in various types of tumors, such as breast, prostate and colon ones [4].

Results and discussion

Identification and quantification of ncRNA biotypes in cell and exosomal fractions from human primary normal and cancer-associated fibroblasts

Primary NFs and CAFs were propagated from 9 colorectal cancer patients (Additional file 1). RNA samples from CAF and NF derived exosomes were combined to get 3 pools of NF-derived exosomes consisting of 3 samples each (NF-EXO samples) vs 3 pools of CAF-derived exosomes consisting of 3 samples each (CAF-EXO samples). RNA samples from cells were combined in the same way (NF-CELL versus CAF-CELL samples).

The expression patterns of 1668 lncRNA and 1761 sncRNA species related to 16 biotypes were analyzed by Next Generation Sequencing and differential expression analyses. The fractions of cell samples (NF-CELL and CAF-CELL) were more similar to each other than to the samples of the exosomal fraction (NF-EXO and CAF-EXO) and vice-versa (Fig. 1). Thus, biotypes such as sense intronic RNAs, lincRNAs and miRNAs appear to be higher in the cell fraction while YRNAs and piRNAs are more enriched in the exosomal fraction. For an absolute and average relative counts of reads mapped to each ncRNA biotype and samples details about the dispersion and biological variation, see Additional files 2, 3 and Additional file 4: Figure S1.
Fig. 1
Fig. 1

Above, 3D line chart, where the most prominent lines based on the average relative counts of reads mapped to each ncRNA biotype in NF-CELL, CAF-CELL, NF-EXO and CAF-EXO samples, are emphasized. Below, two heatmaps with double clustering based on the same information (to the left, heatmap based on lncRNA frequencies; to the right, the one for sncRNAs). The following ranges (0, 0.001), (0.0011, 0.05), (0.1, 0.9) were defined to color the breaks in yellow, orange and red, respectively. Dendrograms were inferred using the complete linkage with the Euclidean distance measure as clustering method

Differential content and enrichment of ncRNAs in fibroblast-derived exosomes compared with their respective NFs or CAFs

Differentially distributed ncRNAs between the exosomal and cellular contents in both NF and CAF data are shown in Fig. 2a and b and Additional files 5 and 6. Samples of the cellular fraction are significantly more enriched in “snoRNAs” than the exosomal fraction, which in turn resulted significantly enriched in YRNAs (Fig. 2c and d). The functional role/s of YRNAs is a current subject of debate. Specifically on cancer research, [5] recently reported enrichment of YRNAs in extracellular vesicles released from the human colon cancer LIM1863 cell line. Our results are consistent with that study, supporting the idea of a potential role for encapsulated YRNAs in exosomes.
Fig. 2
Fig. 2

a Volcano plots from the differential expression analyses of lncRNAs and sncRNAs between NF-CELL versus NF-EXO samples. Red dots mean differentially distributed ncRNAs supported by FDR < 0.05. Below each plot we summarize the results of each analysis. b Volcano plots from the differential expression analyses of lncRNAs and sncRNAs between CAF-CELL and CAF-EXO samples. c Bar-plot showing the results of the differential enrichment analyses of biotypes performed between NF-CELL and NF-EXO samples accompanied by the p values for each biotype. d Results of the differential enrichment analysis of biotypes in CAF-CELL and CAF-EXO. e Venn diagram showing the relationships between significant and non-significant ncRNAs of both analyses, NF-CELL versus NF-EXO and CAF-CELL versus CAF-EXO. f Venn diagram showing the relationships among all over-represented ncRNAs of both comparisons

According to Fig. 2e and f, 2 ncRNAs, respectively expressed by 2 related but different snRNA loci named U1 and RNVU1–17, are over-represented in the cell compartment when the analysis is based on NF cells and their exosomes, while they are over-represented in CAF-derived exosomes when the analysis is performed between CAF and their exosomes, which points to these 2 sncRNAs as potential biomarkers worth investigating. This suggests that the expression patterns and exosomal distribution of differentially expressed ncRNAs are apparently orchestrated by the normal or cancer-associated condition of fibroblast cells (See also Additional files 7 and 8).

Deregulation and selective distribution of ncRNAs in CAFs and derived exosomes in comparison to NFs and their exosomes

The analyses performed between NF-CELL versus CAF-CELL samples showed 2 lncRNAs under-expressed in CAF-CELL samples. Regarding sncRNAs, 2 miRNAs (MIR3935, MIR7853), 1 piRNA (piR-46,427) and 1 snRNA (RNU7-130P) were significantly under-expressed in CAF-CELL samples, while 1 miRNA (miR-21-related) was over-expressed (Fig. 3a and Additional file 9).
Fig. 3
Fig. 3

a Volcano plots from the differential expression analyses of lncRNAs and sncRNAs in NF-CELL and CAF-CELL samples. Red dots are significant ncRNAs supported by FDR < 0.05 and are 2× increased to let the reader see them clearly. Names for significant ncRNAs are also included. b Volcano plots from the differential expression analyses of lncRNAs and sncRNAs in NF-EXO and CAF-EXO samples created using similar criteria for the volcano plot previously shown in plot A. For more information about the 42 significant sncRNAs (red dots) see Additional file 11: Table S1

Fifty-two ncRNAs were differentially distributed between NF- and CAF-derived exosomes (Fig. 3b). 2 were lncRNAs (LINC00326 and WEE2-AS1), which were over represented in CAF-EXO samples (Additional file 10). The remaining 50 are summarized in Additional file 11: Table S1 and correspond to 8 sncRNAs that were under-represented and 42 sncRNAs over-represented in CAF-EXO samples (Additional file 12).

About the 42 sncRNAs over-represented in CAF-EXO, twelve of them correspond to sequences with no homologs, while the remaining 30 split into 4 families, designated as Families 1, 2, 3 and 4, each of which consisted of paralog sequences mapping in different loci (Additional file 13, Additional file 11: Table S1).

Among the over-represented sncRNAs of CAF-EXO samples, 26 snRNAs split into 2 and 23 paralogs, designated as Family 1 and 4, respectively, and a single sequence (RNU11) with no other homologs in our analyses. snRNAs are involved in splicing regulation, which are frequently altered in cancer [6]. Nevertheless, there is not much information about splicing alteration depending on snRNA regulation. Our data show over-distribution of many snRNAs in exosomes derived from CAFs. Remarkably, we also showed in this study that 2 of these sncRNAs, U1 and RNVU1–17, were more represented in CAF-derived exosomes than in their parental cells, but less than in NF-derived exosomes. This suggests that the tumor-associated condition plays a role in the expression patterns of these snRNAs.

We found up to 7 piRNAs over-represented in CAF-EXO samples. Three of them are paralogs constituting Family 2, while the remaining 4 are sequences with no homologs detected in our analyses. piRNAs interact with PIWI proteins and play a role in gene regulation and silencing by methylation or euchromatic histone modifications, which are mechanism clearly linked to the hallmarks of cancer [7], but they are also associated with exosomes and other microvesicles. The number of studies on piRNAs in cancer has increased during the last decade although their specific tumorigenic functions remain unclear [8].

Of the 9 miRNAs observed, 2 of them, mir6087 and its paralog AL353644.1, constitute Family 3. It is worth noting that miR-6087 has been related to colorectal cancer survival [9].

Target genes for over-represented sncRNAs in CAF exosomes and associated gene ontologies

Prediction of target genes for CAF-EXO over-distributed sncRNAs (supported by FDR < 1E-04) is summarized in Additional file 14: Table S2. The piRNAs belonging to Family 2 were predicted to interact with 65 target genes, and similarly piR-57,251, which had 62 target genes. The Family 4 were predicted to interact with 15 target genes; and the RNU11 snRNA, with 7 target genes (remaining predictions are provided as Additional files 10 and 12).

The biological process identified by GO and metabolic pathway annotations (Additional file 15: Figure. S2) represents functions or processes of immediate relevance to tumor progression and microenvironment. Presumably, the effect of CAF-derived exosomes could affect not only tumor cells but also other cells from the tumor microenvironment, prompting tumor progression.

In summary, the ncRNA species act in the regulation of diverse cell functions [10]. Our data support a regulatory role of the deregulated ncRNAs in CAFs, possibly involved with CAFs’ activation process or specifically with the cross-talk between stromal and cancer cells. The levels of some of these regulatory elements are differentially represented in CAF-derived exosomes pointing to the clinical relevance of our study, since the identification of these markers in exosomes could be translated into two kinds of new clinical tools. Firstly, these ncRNA could be potentially used as new biomarkers for patient prognoses and prediction of response to therapy in cancer. But secondarily and more important, cancer exosomes are delivered to circulating body fluids, which gives an opportunity to develop new tools for oncology patients management by minimally invasive blood analyses. Moreover, the low cost together with the feasibility of a non-invasive method offer the chance to define intervals during treatment and the course of the disease to monitor cancer progression and patients’ evolution.

Abbreviations

CAF: 

Cancer associated fibroblast

ECM: 

Extracellular matrix

lncRNA: 

Long non-coding RNA

logFC: 

Log Fold Change

MDS: 

Multidimensional Scalin

miRNA: 

miRNA

ncRNA: 

Non-Coding RNA

NF: 

Normal fibroblast

piRNA: 

PIWI Interacting RNA

RNA-Seq: 

RNA Sequencing

sncRNA: 

Short non-coding RNA

snoRNA: 

Small nucleolar RNA

snRNA: 

Small nuclear RNA

TME: 

Tumor microenvironment

Declarations

Acknowledgements

M. Eaude helped with the English text. We thank lab members for help and advice throughout this research.

Funding

This research is supported by PI12/02037, PI12/01635, PI15/02101, PI17/01847 and RD12/0036/0041 from the Instituto de Salud Carlos III (Plan Estatal de I + D + i 2013–2016) and cofinanced by the European Development Regional Fund “A way to achieve Europe” (ERDF); by “CIBER de Cáncer”, CB16/12/00273, CB16/12/00301, CB16/12/00446 and CB16/12/00291, from the Instituto de Salud Carlos III-FEDER; by the Fundación Científica AECC (a multifaceted approach to target pancreatic cancer); by SAF2010–20750, and SAF 2016-76377-R from the Ministerio de Economía y Competitividad of Spain-FEDER; by S2010/BMD-2344 from the Comunidad de Madrid; and by the Fundación Banco Santander. Cristina Peña is a recipient of a Miguel Servet Contract from the Instituto de Salud Carlos III.

Availability of data and materials

Raw data have been deposited at the NCBI SRA archive with BioProject record PRJNA432148, SRA accession SRP131727 and BioSample records (SAMN08435722, SAMN08435723, SAMN08435724, SAMN08435725).

Authors’ contributions

Conceptualization, MH, FB, VG and CP; Methodology, MH, MR, AH, RR, BG and JE; Software, formal analysis and data curation, CL, MH; Investigation, MH, MR and MJL; Resources, AC, PG and MRG; Writing and visualization, MH, CL, VG and CP; Supervision, VG and CP; Funding acquisition, FB, TC, AC, VG and CP. All authors read and approved the final manuscript.

Ethics approval and consent to participate

Fresh tissue from 9 patients operated for colorectal primary tumors at the Puerta de Hierro-Majadahonda University Hospital was used for the propagation of primary CAFs and NFs. Ethics approval is shown as Supplemental Information. Informed written consent was obtained from all participants after an explanation of the nature of the study, as approved by the Research Ethics Board of Puerta de Hierro Majadahonda University Hospital.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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.

Authors’ Affiliations

(1)
Department of Medical Oncology, Hospital Universitario Puerta de Hierro de Majadahonda, Majadahonda, Madrid, Spain
(2)
Department of Oncology-Pathology, Karolinska Institutet, Stockholm, Sweden
(3)
Biotechvana, Scientific Park, University of Valencia, Valencia, Spain
(4)
Department of Molecular Cell Biology, Institute for Cancer Research, University Hospital-The Norwegian Radium Hospital, and Centre for Cancer Biomedicine, Faculty of Medicine, University of Oslo, Oslo, Norway
(5)
Pathology Department, Fundación Instituto de Investigación Jiménez Díaz, CIBERONC, Madrid, Spain
(6)
Unidad de Genómica, Campus de Cantoblanco, Scientific Park of Madrid, Madrid, Spain
(7)
Grupo de Investigación de oncología traslacional, Departamento de tumores digestivos, Hospital doce de Octubre, Madrid, Spain
(8)
Pathology Department, Hospital Universitario de Guadalajara, Guadalajara, Spain
(9)
Instituto de Investigaciones Biomédicas Alberto Sols, Consejo Superior de Investigaciones Científicas - Universidad Autónoma de Madrid, CIBERONC, Madrid, Spain
(10)
Laboratorio de Oncología Molecular, Hospital Clínico San Carlos. Centro de Investigación Biomédica en Red de Cancer (CIBERONC), Madrid, Spain
(11)
Medical Oncology Department, Ramón y Cajal University Hospital, IRYCIS, CIBERONC, Alcalá University, Madrid, Spain
(12)
Centro de Estudios Biosanitarios, Madrid, Spain

References

  1. Östman A, Augsten M. Cancer-associated fibroblasts and tumor growth – bystanders turning into key players. Curr Opin Genet Dev. 2009;19(1):67–73.View ArticlePubMedGoogle Scholar
  2. Ge R, Tan E, Sharghi-Namini S, Asada HH. Exosomes in Cancer microenvironment and beyond: have we overlooked these extracellular messengers? Cancer Microenviron. 2012;5(3):323–32.View ArticlePubMedPubMed CentralGoogle Scholar
  3. Cheng G. Circulating miRNAs: roles in cancer diagnosis, prognosis and therapy. Adv Drug Deliv Rev. 2015;81:75–93.View ArticlePubMedGoogle Scholar
  4. Prensner JR, Chinnaiyan AM. The emergence of lncRNAs in Cancer biology. Cancer Discov. 2011;1(5):391–407.View ArticlePubMedPubMed CentralGoogle Scholar
  5. Chen M, Xu R, Ji H, Greening DW, Rai A, Izumikawa K, et al. Transcriptome and long noncoding RNA sequencing of three extracellular vesicle subtypes released from the human colon cancer LIM1863 cell line. Sci Rep. 2016;6(1):38397.View ArticlePubMedPubMed CentralGoogle Scholar
  6. Lee SC-W, Abdel-Wahab O. Therapeutic targeting of splicing in cancer. Nat Med. 2016;22(9):976–86.View ArticlePubMedPubMed CentralGoogle Scholar
  7. Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646–74.View ArticlePubMedPubMed CentralGoogle Scholar
  8. Ng KW, Anderson C, Marshall EA, Minatel BC, Enfield KS, Saprunoff HL, et al. Piwi-interacting RNAs in cancer: emerging functions and clinical utility. Mol Cancer. 2016;15(1):5.View ArticlePubMedPubMed CentralGoogle Scholar
  9. Slattery ML, Herrick JS, Pellatt DF, Mullany LE, Stevens JR, Wolff E, et al. Site-specific associations between miRNA expression and survival in colorectal cancer cases. Oncotarget. 2016;7(37):60193–205.View ArticlePubMedPubMed CentralGoogle Scholar
  10. Szymański M, Barciszewska MZ, Zywicki M, Barciszewski J. Noncoding RNA transcripts. J Appl Genet. 2003;44(1):1–19.PubMedGoogle Scholar

Copyright

© The Author(s). 2018

Advertisement