Next Article in Journal
An Expanded Landscape of Unusually Short RNAs in 11 Samples from Six Eukaryotic Organisms
Previous Article in Journal
Comprehensive Transcriptional Profiling and Mouse Phenotyping Reveals Dispensable Role for Adipose Tissue Selective Long Noncoding RNA Gm15551
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Two Exosomal miRNAs in Circulating Blood of Cancer Patients by Using Integrative Transcriptome and Network Analysis

by
Andrés Rincón-Riveros
1,
Josefa Antonia Rodríguez
2,
Victoria E. Villegas
3,* and
Liliana López-Kleine
4,*
1
Bioinformatics and Systems Biology Group, Universidad Nacional de Colombia, Bogotá 111221, Colombia
2
Cancer Biology Research Group, Instituto Nacional de Cancerología, Bogotá 111221, Colombia
3
Centro de Investigaciones en Microbiología y Biotecnología-UR (CIMBIUR), Facultad de Ciencias Naturales, Universidad del Rosario, Bogotá 111221, Colombia
4
Department of Statistics, Faculty of Science, Universidad Nacional de Colombia, Bogotá 111221, Colombia
*
Authors to whom correspondence should be addressed.
Non-Coding RNA 2022, 8(3), 33; https://doi.org/10.3390/ncrna8030033
Submission received: 21 February 2022 / Revised: 20 April 2022 / Accepted: 22 April 2022 / Published: 12 May 2022
(This article belongs to the Topic Intelligent Computing Unlocks the Molecular Code of Complex Diseases)
(This article belongs to the Section Small Non-Coding RNA)

Abstract

:
Exosomes carry molecules of great biological and clinical interest, such as miRNAs. The contents of exosomes vary between healthy controls and cancer patients. Therefore, miRNAs and other molecules transported in exosomes are considered a potential source of diagnostic and prognostic biomarkers in cancer. Many miRNAs have been detected in recent years. Consequently, a substantial amount of miRNA-related data comparing patients and healthy individuals is available, which contributes to a better understanding of the initiation, development, malignancy, and metastasis of cancer using non-invasive sampling procedures. However, a re-analysis of available ncRNA data is rare. This study used available data about miRNAs in exosomes comparing healthy individuals and cancer patients to identify possible global changes related to the presence of cancer. A robust transcriptomic analysis identified two common miRNAs (miR-495-3p and miR-543) deregulated in five cancer datasets. They had already been implicated in different cancers but not reported in exosomes circulating in blood. The study also examined their target genes and the implications of these genes for functional processes.

1. Introduction

In the past three decades, the accumulation of gene expression data in databases al-lowed an extensive study to understand the molecular processes involved in diseases and their diagnosis using transcriptomic biomarkers. It also helped refine treatments and identified pharmacologically interesting genes [1]. High-throughput sequencing has recently allowed identifying and quantifying small RNAs (sRNAs) in eukaryotes, revealing their central role in diseases and their potential use as biomarkers. Since they can be found both in disease-affected tissue and extracellular vesicles circulating in the blood, they are far more useful for this purpose than transcriptomic gene biomarkers [2].
Extracellular vesicles are small structures with a lipid bilayer, produced by the endocytic machinery of the cell and secreted by most nucleated cells [3]. Exosomes and other types of extracellular vesicles circulating in the body fluids of cancer patients could be a source of diagnostic and prognostic cancer biomarkers. During the past decade, there was a dramatic increase in studies demonstrating that extracellular vesicles carry molecules of great biological and clinical interest. For example, exosomes contain non-coding RNAs, especially miRNAs, which have been strongly suggested to be associated with the development of cancer [4,5,6,7,8]. MicroRNAs are small non-coding RNAs of approximately 20–25 nucleotides (nt) in length, with the ability to modulate gene expression by binding to specific regions at the 3′ UTR end of genes, inducing the repression or complete degradation of mRNAs [9,10].
The first step towards a possible use of sRNAs (especially miRNAs) contained in exosomes is to explore available information related to their presence, relationship with disease, and potential use as biomarkers for prognosis, diagnosis, and treatment. Many sRNAs have been detected in recent years; consequently, a substantial amount of sRNA-related data comparing patients and healthy individuals is available at present, especially in cancer. Although there exists some specific knowledge on the role of miRNAs in isolated types of cancer, the use of these datasets for an overall and integrative analysis of vesicle content has been very limited. It can, however, provide significant insight about the potential use of non-coding RNAs contained in them and promote the use of liquid biopsies. Combining this information with evidence on the roles of sRNAs mainly in regulatory processes should deepen our understanding of the initiation, development, malignancy, and metastasis of cancer using noninvasive sampling procedures. This study partly fulfills this need.
Methods for analyzing gene expression data need to be adapted and integrated properly with information in the literature and databases to clarify and develop biological hypotheses about the roles of miRNAs in cancer [11,12]. In this study, we analyzed available data about miRNAs contained in exosomes comparing healthy individuals and cancer patients to identify possible global changes related to the presence of cancer. Using robust transcriptome analysis, we identified two common miRNAs (miR-495-3p and miR-543) deregulated in all five cancer datasets analyzed in this work. They had already been implicated in different cancer types but not reported in exosomes circulating in blood. These common miRNAs can serve as a starting point for monitoring tumorigenesis by studying their target genes and the involvement of these genes in functional processes. We, thus, investigated the molecular processes possibly affected by them by identifying their target genes using a network-based approach. We also discussed how the direction of deregulation in each cancer type and the molecular functions affected by their target genes make these two miRNAs potentially useful biomarkers in cancer.

2. Results

2.1. miRNA Expression in Different Types of Cancer

The established criteria for identifying miRNAs with differential expression in different types of cancer included a fold change >1 or ≤1 and a p-value adjusted by the FDR method <0.2. We chose this soft threshold given that our main goal was to globally identify potential central miRNAs common to several and very different types of cancer (Table 1) in the sense of a meta-analysis. Table S1 summarizes the results of the differential expression analysis.
To determine a possible miRNA master regulator between different tumor types, we inspected the differential expression datasets of these seven cancer types using a Venn diagram (Figure S1). In these datasets, we found that at least one miRNA with altered expression in common was found among five types of cancer. In gastric cancer and colorectal cancer, miR-495-3p and miR-543 were overexpressed, while in multiple myeloma, only miR-495-3p was overexpressed. In contrast, in prostate cancer, miR-495-3p and miR-543 were found to be deregulated; in glioblastoma, only miR-543 was deregulated (Table 2).

2.2. External Validation of Differentially Expressed miRNAs

We used the UALCAN tool [13] (http://ualcan.path.uab.edu/analysis.html, accessed on 20 March 2022) to test the expression of common miRNAs in the studied cancer types (except multiple myeloma, which is not available) (Table 2). Next, the findings were validated, including the direction of the deregulation of miRNAs present in extracellular vesicles.
Subsequently, a network was built with miRNAs in exosomes that overlapped among the five types of cancer in the miRNet web application (https://www.mirnet.ca/ accessed on 20 March 2022), to search for the target genes of these miRNAs.

2.3. miRNA–mRNA Bipartite Interaction Network

The constructed miRNA–mRNA network had 1121 nodes and 1255 edges based on information from the literature in the miRNet database. For miR-495-3p, we found connections with respect to 981 genes; on the other hand, for miR-543, there were connections with respect to 183 genes. Notably, two exosomal miRNAs’ shared connections to 45 genes (Figure S2) were as follows: HMGA2, MIEF1, FNDC3B, ZNF703, KLHL15, SIGLEC14, AMD1, BACH1, CLTC, DDX3X, NOTCH2, NRAS, UBL3, PTP4A1, STX7, MATR3, WDFY3, ARL6IP1, C16orf72, LIN7C, GPCPD1, FBXO11, ADAM22, TMEM33, COL12A1, ANKRD13C, YPEL5, ACSL3, ZBTB1, BTBD3, DYNC1LI2, EGLN1, SERPINE2, RPLP1, ARHGAP29, SMARCA5, GNAQ, ZNF281, IRGQ, KIF5C, CFL1, NAP1L1, F1, PGAP1, and GFPT1.
Of these common target genes of miR-495-3p and miR-543, NRAS, HMGA2, and EGLN1 present in KEGG cancer metabolic pathways stand out (p-value < 0.05) (Figure S3). The top three target genes with the highest degree and betweenness in the bipartite network were AKT1, JUN, and GSK3B, which were involved in the miR-495-3p community. Interestingly, these hubs were not present in the large community detected with the Louvain algorithm. However, GSK3B (Glycogen Synthase Kinase 3 Beta) acts as a negative regulator in the phosphorylation of key cancer-related genes, such as APC, JUN, and CTNNB1/beta-catenin. For its part, AKT1 is a regulator of GSK3B [14,15].

2.4. Community Detection and Functional Enrichment Analysis

We used the Louvain algorithm to detect communities in the bipartite network. The community associated with miRNA-495-3p was the largest one, and it was further employed to perform functional enrichment analysis. For analysis, we selected overrepresented KEGG pathways and GO terms identified as significant based on Fisher’s exact test with a p-value adjusted by the Benjamini–Hochberg method of less than 0.05. Table 3 lists the top six metabolic pathways (Figure S4). The GO terms in the category of biological process (BP) with the highest representation among the genes of this community were synapse organization, axo-dendritic transport, RNA splicing, regulation of RNA splicing, and dendrite development. A network was built with the most representative terms to assess the principal biological process represented in miRNA–mRNA targets. It showed how these genes connect with other BPs of interest in transcriptomic and RNA regulation processes, which can be affected in diverse cancer types (Figure 1).
For the enrichment of KEGG metabolic pathways, we identified genes annotated to more than one metabolic pathway were also found, demonstrating the complexity of the biological processes. Among the top five metabolic pathways with a several members and low associated adjusted p-values, there were interconnections between potential target genes of miR-495-3p and ncRNA participating in several metabolic functions (Figure 2).
The community with the second largest number of members was the one associated with miRNA-543. This community had 70 members, on which functional enrichment was carried out with enrichKEGG from the clusterProfiler package. For this community, we only identified one metabolic pathway was identified with adjusted p-value < 0.05, namely, SNARE interactions in vesicular transport (SNAP23/STX4/STX7/STX1B). This pathway is composed of a family of proteins that participates in the biogenesis and secretion of extracellular vesicles, which has been studied in various fields such as autophagy and diseases such as Parkinson’s [16,17,18].

3. Discussion

Epigenetic modifications are closely associated with the development of cancer. Specifically, histone modifications, DNA methylation, and the regulation of gene expression by ncRNAs have been revealed to be related to this disease [19,20]. An miRNA can have hundreds or thousands of mRNAs targets; in turn, a single gene can be modulated by several miRNAs [21]. Therefore, it is necessary to bear in mind that the molecular signature of each type of tumor is unique. Cancer is considered a heterogeneous disease, due to which it is crucial to start focusing efforts on deciphering common patterns of different tumors [22,23].
Furthermore, miRNAs are closely related to other species of ncRNAs such as lncRNAs and circRNAs, generating a close network of interactions and competition for molecular targets. This network forms the basis for Salmena’s competitive endogenous RNA (ceRNA) hypothesis proposed in 2011 [24,25].
Given the complexity of molecular interactions, focusing only on a single type of biomolecule for the understanding of a complex disease such as cancer creates difficulties in the understanding of biological processes and disease development. Thus, it is vital to apply systems biology and network analysis, including protein–protein interaction networks, metabolites, and transcription factor genes [26,27,28]. Recent advances in systems biology and bioinformatics have also increased our understanding of heterogeneous interaction networks in cancer [29] and have highlighted the importance of miRNAs in cancer. Therefore, identifying two common miRNAs in five cancer types is an important starting point to develop biomarkers for diagnosis and treatment tuning.
Notably, we found that two miRNAs circulating in exosomes were upregulated in each cancer studied here in comparison to levels in healthy patients: hsa-miR-543 and hsa-miR-495-3p. External validation with TCGA data using the UALCAN tool highlighted that screening deregulated miRNAs in available databases generates valuable results. Although the literature has reported that hsa-miR-543 and hsa-miR-495-3p are involved in different types of cancer, there are no mentions of their detection in circulating exosomes as potential oncogenic markers. miR-543 has already been shown to suppress breast tumor cell viability, proliferation, and progression by repressing VCAN [30]. miR-543 also plays an oncogenic role in prostate cancer cells by suppressing Numb and promoting tumor growth, metastasis, and the acquisition of stem cell-like traits [31]. The knowledge that this miRNA is deregulated in different cancer types but is also related to other genes not directly involved in cancer renders this non-coding RNA a promising candidate for monitoring tumor-related processes; therefore, it should be a subject of further investigation [32]. CircTLK1 sequesters miR-495-3p by sponging, thus contributing to tumor growth and metastasis in renal cell carcinoma [33]. In contrast, the hypermethylation of the miR-495-3p promoter increases the expression of at least 10 epigenetically modified oncogenes that are overexpressed in gastric cancer [34]. Moreover, studies have demonstrated the tumor-suppressing potential of miR-495-3p in stomach, melanoma, and prostate tumors, which raises many questions about the role of this miRNA in all cancers analyzed here [35,36,37].
Our findings showed that miR-495-3p and miR-543, recognized modulators of carcinogenesis, were downregulated in prostate cancer. This result was possibly due to competition with other ncRNAs, such as NORAD and MCM3AP-AS1 lncRNAs that silence these miRNAs by increasing the expression of TRIP13 and the SLC39A10/PTEN/Akt axis, accelerating tumor progression [37,38]. We also found deregulated miR-543 in glioblastoma, which has been reported as a preoperative and classification biomarker in glioma [39]. In contrast, our results indicated that miR-495-3p and miR-543 are upregulated in gastric and colorectal cancer, in line with some published reports. In colorectal cancer, the overexpression of miR-543 was reported to increase chemoresistance by blocking tumor suppressor PTEN [40,41]. This miRNA also suppresses histone deacetylase SIRT1 in gastric tumor cells, increasing proliferation and tumor progression [42]. As for multiple myeloma, recent studies have indicated the importance of using non-invasive biomarkers to study this disease [43]. We found that miR-495-3p was deregulated, which has been proposed as a tumor suppressor in some hematological malignancies, such as acute myeloid leukemia and mixed lineage leukemia [44]. Previous studies suggested that miR495-3p modulates the expression of aquaporin-1 (AQP1), a protein with a potential role in osteosarcoma and multiple myeloma development [45,46,47].
To further investigate the role of these two common miRNAs, we identified their target genes using the miRNet tool. Experimental evidence allowed us to build a network where nodes were miRNAs and their target genes and edges represented interactions between them. Based on this network, network theory metrics were applied to find miRNA-mRNA regulatory patterns. The application of modularity algorithms for detecting communities allowed the identification of groups with a high density of connections. These groups of clustered genes associated with hsa-miR-543 and hsa-miR-495-3p indicated a closer community of relevant target genes, which helped clarify the potential global role of these target genes in cancer and physiological processes.
The metabolic pathways associated with the miRNA-495-3p community have already been described in cancer development. Gamma-aminobutyric acid (GABA) synapse is one of the main metabolic pathways of neurotransmitters in mammals, with implications in regulating inflammation and immune response. It has been found to promote cell proliferation due to the overexpression of the GABAA receptor activating MAPK signaling in brain, gastric, breast, and prostate cancers [48,49,50].

4. Material and Methods

In this study, we reanalyzed openly available data on miRNAs contained in exosomes from recognized genomic databases. We examined datasets of cancer-associated exo-somes circulating in human plasma or serum, from tumor samples and normal controls, to detect differentially expressed miRNAs. We proceeded to predict the target genes of the most significant miRNAs, using an miRNA-gene bipartite network to illustrate the most important findings regarding the roles of miRNAs in biological functions and metabolic pathways (Figure 3).

4.1. High-Throughput Gene Expression Data Retrieval

For the present study, we retrieved datasets from the Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo, accessed on 1 June 2021), Expression Atlas (https://www.ebi.ac.uk/gxa, accessed on 5 June 2021), and Genomic Expression Archive (GEA, https://www.ddbj.nig.ac.jp/gea, accessed on 5 June 2021). We also applied specialized search tools, such as Google Dataset Search (https://datasesearch.reearch.google.com, accessed on 6 June 2021), OmicsDI (http://www.omicsdi.org, accessed on 6 June 2021), and DataMed (https://datamed.org, accessed on 7 June 2021), using the following keywords: “cancer”, “exosomes”, “RNAseq”, and “exosomal RNA”. Table 4 provides detailed information on the six RNAseq datasets selected for this study, with seven representative cancer types. All selected datasets had respective healthy or normal controls to perform differential expression analysis.

4.2. Data Processing

Each selected database was subject to quality control before processing, which involved filtering out miRNAs with low counts and outliers and independent normalization. To identify potentially deregulated miRNAs, we used the R DESeq2 and siggenes packages [55,56]. For this first filter of potentially deregulated miRNAs, we established a False Discovery Rate (FDR) of 20% and a fold-change threshold of >1.0.

4.3. Bipartite miRNA–Gene Network

Using miRNAs that passed the first filter in the five types of cancer studied to predict possible target genes, we applied the miRNet web tool (http://www.mirnet.ca/, accessed on 20 August 2021). In association with 20 databases, this tool forms a unique database of information on exosomes in various species [57]. The miRNA–mRNA interactions found in miRNet were processed with R package igraph [58]. Subsequently, we built a network with the miRNAs of exosomes overlapping among the five types of cancer in the miRNet web application, looking for the target genes of these miRNAs. As a source of the miRNAs, we selected “exosomes”, which is an advantage of miRNet and is achieved by linking ExoCarta to the analysis [59]. The versions used were miRTarBase v8.0 [60] and TarBase v8.0 [61] for gene search.

4.4. Network Community Detection

The basic topological properties of the network, such as degree (representing the level of connectedness), were analyzed using the R package igraph. For the detection of communities, we evaluated three algorithms, Louvain [62], label propagation [63], and Walktrap [64], in the same R package [58]. The communities found in the constructed interaction network were identified using the modularity index with respect to their members [65].
Q = i = 1 k ( e i i a i 2 )  

4.5. Enrichment Analysis

We analyzed metabolic pathway overrepresentation (KEGG) and Gene Ontology (GO) using the enrichKEGG (KEGG metabolic pathways), enrichMKEGG (KEGG modules), and enrichGO functions of the R clusterProfiler package [66], after collecting terms and pathways with an adjusted p-value < 0.05, a minimum count of 10, and an enrichment factor > 1.5 (ratio between the observed counts and the counts expected by chance) for analysis.

5. Conclusions

Our study provides an integrative analysis of miRNA targets expressed in extracellular vesicles in different types of cancer, their possible relationships with vital physiological functions, and their potential roles in cancer. miRNAs and their relationships with target genes of diverse functional groups are valuable for diagnostic and therapeutic purposes. Beyond acting as potential transcriptomic biomarkers, they serve for monitoring disease status and treatment tuning because they are polyfunctional molecules: in certain types of tumors, they can play an oncogenic role (oncomiR) while acting as tumor suppressors in others. Therefore, there is much interest in further investigating the content of vesicles, including more data sets, which has not happened very often, and wet lab experiments. The use of miRNAs as biomarkers of cancer onset or progression and in monitoring therapeutic outcomes is promising. Nevertheless, considering their multiple roles and partners, their use needs to be preceded by large-scale prospective studies on various types of cancer to identify the molecular signatures of each disease.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ncrna8030033/s1: Table S1. DE across seven types of cancer; Figure S1. A Venn diagram was elaborated as a tool to find miRNAs in common between different datasets Figure S2. The complete miRNA–mRNA network based on interaction data from miRNet; Figure S3. Three miRNAs in common to miR-495-3p and miR-543 in KEGG cancer pathways; Figure S4. KEGG metabolic pathways analysis of the target genes; Supplementary file S2. miRNA–mRNA bipartite interaction network.

Author Contributions

Conceptualization: A.R.-R., J.A.R., V.E.V. and L.L.-K.; methodology: A.R.-R., V.E.V. and L.L.-K.; formal analysis: A.R.-R. and L.L.-K.; writing and drafting of the manuscript: A.R.-R., J.A.R., V.E.V. and L.L.-K.; writing—review and scientific editing: A.R.-R. and L.L.-K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Instituto Nacional de Cancerología Grant Number XRPM:C19010300-454 (to J.A.R.), a Colombian Ministry of Sciences scholarship (to A.R.-R.), and Grant Number IV-FPC007 from the Universidad del Rosario, Bogotá (to V.E.V.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data used are publicly available from the URLs included in the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

3′ UTR3′ Untranslated region
AMPKAMP-activated protein kinase
ceRNAcompetitive endogenous RNA
circRNACircular RNA
CRCColorectal cancer
FDRFalse Discovery Rate
GABAGamma-aminobutyric acid
GEAGenomic Expression Archive
GEOGene Expression Omnibus
GOGene Ontology
KEGGKyoto Encyclopedia of Genes and Genomes pathways
lncRNALong non-coding RNA
MAPKMitogen-activated protein kinase
miRmiRNA
ncRNAnon-coding RNA
RNARibonucleic acid
SNARESoluble N-ethylmaleimide-sensitive-factor attachment protein receptor
sRNAsmall RNA

References

  1. Tsimberidou, A.M.; Fountzilas, E.; Bleris, L.; Kurzrock, R. Transcriptomics and solid tumors: The next frontier in precision cancer medicine. Semin. Cancer Biol. 2020, 17. [Google Scholar] [CrossRef] [PubMed]
  2. Beckers, M.; Mohorianu, I.; Stocks, M.; Applegate, C.; Dalmay, T.; Moulton, V. Comprehensive processing of high-throughput small RNA sequencing data including quality checking, normalization, and differential expression analysis using the UEA sRNA Workbench. RNA 2017, 23, 823–835. [Google Scholar] [CrossRef] [Green Version]
  3. Kowal, J.; Tkach, M.; Théry, C. Biogenesis and secretion of exosomes. Curr. Opin. Cell Biol. 2014, 29, 116–125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Vykoukal, J.; Sun, N.; Aguilar-Bonavides, C.; Katayama, H.; Tanaka, I.; Fahrmann, J.F.; Capello, M.; Fujimoto, J.; Aguilar, M.; Wistuba, I.I.; et al. Plasma-derived extracellular vesicle proteins as a source of biomarkers for lung adenocarcinoma. Oncotarget 2017, 8, 95466–95480. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Semina, E.V.; Rysenkova, K.D.; Troyanovskiy, K.E.; Shmakova, A.A.; Rubina, K.A. MicroRNAs in Cancer: From Gene Expression Regulation to the Metastatic Niche Reprogramming. Biochemistry 2021, 86, 785–799. [Google Scholar] [CrossRef] [PubMed]
  6. Lima, L.G.; Ham, S.; Shin, H.; Chai, E.P.Z.; Lek, E.S.H.; Lobb, R.J.; Müller, A.F.; Mathivanan, S.; Yeo, B.; Choi, Y.; et al. Tumor microenvironmental cytokines bound to cancer exosomes determine uptake by cytokine receptor-expressing cells and biodistribution. Nat. Commun. 2021, 12, 3543. [Google Scholar] [CrossRef] [PubMed]
  7. Muluhngwi, P.; Klinge, C.M. Identification and Roles of miR-29b-1-3p and miR29a-3p-Regulated and Non-Regulated lncRNAs in Endocrine-Sensitive and Resistant Breast Cancer Cells. Cancers 2021, 13, 3530. [Google Scholar] [CrossRef]
  8. Romano, R.; Picca, A.; Eusebi, L.H.; Marzetti, E.; Calvani, R.; Moro, L.; Bucci, C.; Guerra, F. Extracellular Vesicles and Pancreatic Cancer: Insights on the Roles of miRNA, lncRNA, and Protein Cargos in Cancer Progression. Cells 2021, 10, 1361. [Google Scholar] [CrossRef]
  9. Angius, A.; Cossu-Rocca, P.; Arru, C.; Muroni, M.R.; Rallo, V.; Carru, C.; Uva, P.; Pira, G.; Orrù, S.; De Miglio, M.R. Modulatory Role of microRNAs in Triple Negative Breast Cancer with Basal-Like Phenotype. Cancers 2020, 12, 3298. [Google Scholar] [CrossRef]
  10. Torii, C.; Maishi, N.; Kawamoto, T.; Morimoto, M.; Akiyama, K.; Yoshioka, Y.; Minami, T.; Tsumita, T.; Alam, M.T.; Ochiya, T.; et al. miRNA-1246 in extracellular vesicles secreted from metastatic tumor induces drug resistance in tumor endothelial cells. Sci. Rep. 2021, 11, 13502. [Google Scholar] [CrossRef]
  11. Stocks, M.B.; Mohorianu, I.; Beckers, M.; Paicu, C.; Moxon, S.; Thody, J.; Dalmay, T.; Moulton, V. The UEA sRNA Workbench (version 4.4): A comprehensive suite of tools for analyzing miRNAs and sRNAs. Bioinformatics 2018, 34, 3382–3384. [Google Scholar] [CrossRef] [PubMed]
  12. Fallmann, J.; Videm, P.; Bagnacani, A.; Batut, B.; Doyle, M.A.; Klingstrom, T.; Eggenhofer, F.; Stadler, P.F.; Backofen, R.; Grüning, B. The RNA workbench 2.0: Next generation RNA data analysis. Nucleic Acids Res. 2019, 47, W511–W515. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Chandrashekar, D.S.; Bashel, B.; Balasubramanya, S.A.H.; Creighton, C.J.; Ponce-Rodriguez, I.; Chakravarthi, B.V.S.K.; Varambally, S. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia 2017, 19, 649–658. [Google Scholar] [CrossRef] [PubMed]
  14. Racaud-Sultan, C.; Vergnolle, N. GSK3β, a Master Kinase in the Regulation of Adult Stem Cell Behavior. Cells 2021, 10, 225. [Google Scholar] [CrossRef] [PubMed]
  15. Qi, X.-J.; Wildey, G.M.; Howe, P.H. Evidence That Ser87 of BimEL Is Phosphorylated by Akt and Regulates BimEL Apoptotic Function *. J. Biol. Chem. 2006, 281, 813–823. [Google Scholar] [CrossRef] [Green Version]
  16. Kimura, T.; Jia, J.; Kumar, S.; Choi, S.W.; Gu, Y.; Mudd, M.; Dupont, N.; Jiang, S.; Peters, R.; Farzam, F.; et al. Dedicated SNAREs and specialized TRIM cargo receptors mediate secretory autophagy. EMBO J. 2017, 36, 42–60. [Google Scholar] [CrossRef]
  17. Yu, H.; Sun, T.; An, J.; Wen, L.; Liu, F.; Bu, Z.; Cui, Y.; Feng, J. Potential Roles of Exosomes in Parkinson’s Disease: From Pathogenesis, Diagnosis, and Treatment to Prognosis. Front. Cell Dev. Biol. 2020, 8, 86. [Google Scholar] [CrossRef] [Green Version]
  18. Yu, Z.; Shi, M.; Stewart, T.; Fernagut, P.-O.; Huang, Y.; Tian, C.; Dehay, B.; Atik, A.; Yang, D.; De Giorgi, F.; et al. Reduced oligodendrocyte exosome secretion in multiple system atrophy involves SNARE dysfunction. Brain 2020, 143, 1780–1797. [Google Scholar] [CrossRef]
  19. Jung, G.; Hernández-Illán, E.; Moreira, L.; Balaguer, F.; Goel, A. Epigenetics of colorectal cancer: Biomarker and therapeutic potential. Nat. Rev. Gastroenterol. Hepatol. 2020, 17, 111–130. [Google Scholar] [CrossRef]
  20. Lee, J.E.; Kim, M.-Y. Cancer epigenetics: Past, present and future. Semin. Cancer Biol. 2021, 31, S1044-579X(21)00075-4. [Google Scholar] [CrossRef]
  21. Dragomir, M.; Mafra, A.C.P.; Dias, S.M.G.; Vasilescu, C.; Calin, G.A. Using microRNA Networks to Understand Cancer. Int. J. Mol. Sci. 2018, 19, 1871. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Zhang, C.; Yang, J.; Chen, Y.; Jiang, F.; Liao, H.; Liu, X.; Wang, Y.; Kong, G.; Zhang, X.; Li, J.; et al. miRNAs derived from plasma small extracellular vesicles predict organo-tropic metastasis of gastric cancer. Gastric Cancer 2022, 25, 360–374. [Google Scholar] [CrossRef]
  23. Pisarska, J.; Baldy-Chudzik, K. MicroRNA-Based Fingerprinting of Cervical Lesions and Cancer. J. Clin. Med. 2020, 9, 3668. [Google Scholar] [CrossRef]
  24. Salmena, L.; Poliseno, L.; Tay, Y.; Kats, L.; Pandolfi, P.P. A ceRNA Hypothesis: The Rosetta Stone of a Hidden RNA Language? Cell 2011, 146, 353–358. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Zhou, R.-S.; Zhang, E.-X.; Sun, Q.-F.; Ye, Z.-J.; Liu, J.-W.; Zhou, D.-H.; Tang, Y. Integrated analysis of lncRNA-miRNA-mRNA ceRNA network in squamous cell carcinoma of tongue. BMC Cancer 2019, 19, 779. [Google Scholar] [CrossRef] [PubMed]
  26. Barabási, A.-L.; Oltvai, Z.N. Network biology: Understanding the cell’s functional organization. Nat. Rev. Genet. 2004, 5, 101–113. [Google Scholar] [CrossRef]
  27. Masoomy, H.; Askari, B.; Tajik, S.; Rizi, A.K.; Jafari, G.R. Topological analysis of interaction patterns in cancer-specific gene regulatory network: Persistent homology approach. Sci. Rep. 2021, 11, 16414. [Google Scholar] [CrossRef]
  28. Otálora-Otálora, B.A.; Florez, M.; López-Kleine, L.; Canas Arboleda, A.; Grajales Urrego, D.M.; Rojas, A. Joint Transcriptomic Analysis of Lung Cancer and Other Lung Diseases. Front. Genet. 2019, 10, 1260. [Google Scholar] [CrossRef] [Green Version]
  29. Capobianco, E. Next Generation Networks: Featuring the Potential Role of Emerging Applications in Translational Oncology. J. Clin. Med. 2019, 8, 664. [Google Scholar] [CrossRef] [Green Version]
  30. Li, R.; Hou, S.; Zou, M.; Ye, K.; Xiang, L. miR-543 impairs cell proliferation, migration, and invasion in breast cancer by suppressing VCAN. Biochem. Biophys. Res. Commun. 2021, 570, 191–198. [Google Scholar] [CrossRef]
  31. Wang, X.; Yang, J.-Y.; Cai, J.; Zhang, D.-J.; Zhao, L.; Luo, L.-H.; Xiong, Y.; Zhang, T.; Jin, M. MiR-543/Numb promotes proliferation, metastasis, and stem-like cell traits of prostate cancer cells. Am. J. Transl. Res. 2021, 13, 617–631. [Google Scholar] [PubMed]
  32. Zhou, C.; Zhao, X.; Duan, S. The role of miR-543 in human cancerous and noncancerous diseases. J. Cell. Physiol. 2021, 236, 15–26. [Google Scholar] [CrossRef] [PubMed]
  33. Lei, X.; Yang, M.; Xiao, Z.; Zhang, H.; Tan, S. circTLK1 facilitates the proliferation and metastasis of renal cell carcinoma by regulating miR-495-3p/CBL axis. Open Life Sci. 2021, 16, 362–374. [Google Scholar] [CrossRef] [PubMed]
  34. Yuan, X.; Yang, T.; Xu, Y.; Ou, S.; Shi, P.; Cao, M.; Zuo, X.; Liu, Q.; Yao, J. SNHG10 Promotes Cell Proliferation and Migration in Gastric Cancer by Targeting miR-495-3p/CTNNB1 Axis. Dig. Dis. Sci. 2021, 66, 2627–2636. [Google Scholar] [CrossRef]
  35. Eun, J.W.; Kim, H.S.; Shen, Q.; Yang, H.D.; Kim, S.Y.; Yoon, J.H.; Park, W.S.; Lee, J.Y.; Nam, S.W. MicroRNA-495-3p functions as a tumor suppressor by regulating multiple epigenetic modifiers in gastric carcinogenesis. J. Pathol. 2018, 244, 107–119. [Google Scholar] [CrossRef] [PubMed]
  36. Xia, Y.; Zhou, Y.; Han, H.; Li, P.; Wei, W.; Lin, N. lncRNA NEAT1 facilitates melanoma cell proliferation, migration, and invasion via regulating miR-495-3p and E2F3. J. Cell. Physiol. 2019, 234, 19592–19601. [Google Scholar] [CrossRef]
  37. Chen, F.; Liu, L.; Wang, S. Long non-coding RNA NORAD exhaustion represses prostate cancer progression through inhibiting TRIP13 expression via competitively binding to miR-495-3p. Cancer Cell Int. 2020, 20, 323. [Google Scholar] [CrossRef]
  38. Jia, Z.; Li, W.; Bian, P.; Liu, H.; Pan, D.; Dou, Z. LncRNA MCM3AP-AS1 Promotes Cell Proliferation and Invasion Through Regulating miR-543-3p/SLC39A10/PTEN Axis in Prostate Cancer. OncoTargets Ther. 2020, 13, 9365–9376. [Google Scholar] [CrossRef]
  39. Ebrahimkhani, S.; Vafaee, F.; Hallal, S.; Wei, H.; Lee, M.Y.T.; Young, P.E.; Satgunaseelan, L.; Beadnall, H.; Barnett, M.H.; Shivalingam, B.; et al. Deep sequencing of circulating exosomal microRNA allows non-invasive glioblastoma diagnosis. NPJ Precis. Oncol. 2018, 2, 28. [Google Scholar] [CrossRef] [Green Version]
  40. Liu, G.; Zhou, J.; Dong, M. Down-regulation of miR-543 expression increases the sensitivity of colorectal cancer cells to 5-Fluorouracil through the PTEN/PI3K/AKT pathway. Biosci. Rep. 2019, 39, BSR20190249. [Google Scholar] [CrossRef] [Green Version]
  41. Sun, J.; Zhou, J.; Dong, M.; Sheng, W. Dysregulation of MicroRNA-543 expression in colorectal cancer promotes tumor migration and invasion. Mol. Carcinog. 2017, 56, 250–257. [Google Scholar] [CrossRef] [PubMed]
  42. Li, J.; Dong, G.; Wang, B.; Gao, W.; Yang, Q. miR-543 promotes gastric cancer cell proliferation by targeting SIRT1. Biochem. Biophys. Res. Commun. 2016, 469, 15–21. [Google Scholar] [CrossRef] [PubMed]
  43. Manier, S.; Liu, C.-J.; Avet-Loiseau, H.; Park, J.; Shi, J.; Campigotto, F.; Salem, K.Z.; Huynh, D.; Glavey, S.V.; Rivotto, B.; et al. Prognostic role of circulating exosomal miRNAs in multiple myeloma. Blood 2017, 129, 2429–2436. [Google Scholar] [CrossRef] [PubMed]
  44. Jiang, X.; Huang, H.; Li, Z.; He, C.; Li, Y.; Chen, P.; Gurbuxani, S.; Arnovitz, S.; Hong, G.-M.; Price, C.; et al. miR-495 is a tumor-suppressor microRNA down-regulated in MLL-rearranged leukemia. Proc. Natl. Acad. Sci. USA 2012, 109, 19397–19402. [Google Scholar] [CrossRef] [Green Version]
  45. Pelagalli, A.; Nardelli, A.; Fontanella, R.; Zannetti, A. Inhibition of AQP1 Hampers Osteosarcoma and Hepatocellular Carcinoma Progression Mediated by Bone Marrow-Derived Mesenchymal Stem Cells. Int. J. Mol. Sci. 2016, 17, 1102. [Google Scholar] [CrossRef] [Green Version]
  46. Vacca, A.; Frigeri, A.; Ribatti, D.; Nicchia, G.P.; Nico, B.; Ria, R.; Svelto, M.; Dammacco, F. Microvessel overexpression of aquaporin 1 parallels bone marrow angiogenesis in patients with active multiple myeloma. Br. J. Haematol. 2001, 113, 415–421. [Google Scholar] [CrossRef] [PubMed]
  47. Zhu, L.; Lin, Z.-W.; Wang, G.; Zhang, H.; Liu, B.; Xu, Q.-J. MicroRNA-495 downregulates AQP1 and facilitates proliferation and differentiation of osteoblasts in mice with tibial fracture through activation of p38 MAPK signaling pathway. Sci. Rep. 2019, 9, 16171. [Google Scholar] [CrossRef] [Green Version]
  48. Takehara, A.; Hosokawa, M.; Eguchi, H.; Ohigashi, H.; Ishikawa, O.; Nakamura, Y.; Nakagawa, H. γ-Aminobutyric Acid (GABA) Stimulates Pancreatic Cancer Growth through Overexpressing GABA-A Receptor π Subunit. Cancer Res. 2007, 67, 9704–9712. [Google Scholar] [CrossRef] [Green Version]
  49. Xia, S.; He, C.; Zhu, Y.; Wang, S.; Li, H.; Zhang, Z.; Jiang, X.; Liu, J. GABA(B)R-Induced EGFR Transactivation Promotes Migration of Human Prostate Cancer Cells. Mol. Pharmacol. 2017, 92, 265–277. [Google Scholar] [CrossRef] [Green Version]
  50. Jiang, S.-H.; Hu, L.-P.; Wang, X.; Li, J.; Zhang, Z.-G. Neurotransmitters: Emerging targets in cancer. Oncogene 2020, 39, 503–515. [Google Scholar] [CrossRef]
  51. Tang, S.; Cheng, J.; Yao, Y.; Lou, C.; Wang, L.; Huang, X.; Zhang, Y. Combination of Four Serum Exosomal MiRNAs as Novel Diagnostic Biomarkers for Early-Stage Gastric Cancer. Front. Genet. 2020, 11, 237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Yao, B.; Qu, S.; Hu, R.; Gao, W.; Jin, S.; Liu, M.; Zhao, Q. A panel of miRNAs derived from plasma extracellular vesicles as novel diagnostic biomarkers of lung adenocarcinoma. FEBS Open Bio 2019, 9, 2149–2158. [Google Scholar] [CrossRef] [PubMed]
  53. Tan, C.; Cao, J.; Chen, L.; Xi, X.; Wang, S.; Zhu, Y.; Yang, L.; Ma, L.; Wang, D.; Yin, J.; et al. Noncoding RNAs Serve as Diagnosis and Prognosis Biomarkers for Hepatocellular Carcinoma. Clin. Chem. 2019, 65, 905–915. [Google Scholar] [CrossRef] [PubMed]
  54. Yuan, T.; Huang, X.; Woodcock, M.; Du, M.; Dittmar, R.; Wang, Y.; Tsai, S.; Kohli, M.; Boardman, L.; Patel, T.; et al. Plasma extracellular RNA profiles in healthy and cancer patients. Sci. Rep. 2016, 6, 19413. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Tusher, V.G.; Tibshirani, R.; Chu, G. Significance analysis of microarrays applied to the ionizing radiation response. Proc. Natl. Acad. Sci. USA 2001, 98, 5116–5121. [Google Scholar] [CrossRef] [Green Version]
  57. Chang, L.; Zhou, G.; Soufan, O.; Xia, J. miRNet 2.0: Network-based visual analytics for miRNA functional analysis and systems biology. Nucleic Acids Res. 2020, 48, W244–W251. [Google Scholar] [CrossRef]
  58. Csardi, G.; Nepusz, T. The igraph software package for complex network research. InterJournal Complex Syst. 2006, 1695, 1–9. [Google Scholar]
  59. Keerthikumar, S.; Chisanga, D.; Ariyaratne, D.; Al Saffar, H.; Anand, S.; Zhao, K.; Samuel, M.; Pathan, M.; Jois, M.; Chilamkurti, N.; et al. ExoCarta: A Web-Based Compendium of Exosomal Cargo. J. Mol. Biol. 2016, 428, 688–692. [Google Scholar] [CrossRef] [Green Version]
  60. Huang, H.-Y.; Lin, Y.-C.-D.; Li, J.; Huang, K.-Y.; Shrestha, S.; Hong, H.-C.; Tang, Y.; Chen, Y.-G.; Jin, C.-N.; Yu, Y.; et al. miRTarBase 2020: Updates to the experimentally validated microRNA–target interaction database. Nucleic Acids Res. 2020, 48, D148–D154. [Google Scholar] [CrossRef] [Green Version]
  61. Karagkouni, D.; Paraskevopoulou, M.D.; Chatzopoulos, S.; Vlachos, I.S.; Tastsoglou, S.; Kanellos, I.; Papadimitriou, D.; Kavakiotis, I.; Maniou, S.; Skoufos, G.; et al. DIANA-TarBase v8: A decade-long collection of experimentally supported miRNA-gene interactions. Nucleic Acids Res. 2018, 46, D239–D245. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Blondel, V.D.; Guillaume, J.-L.; Lambiotte, R.; Lefebvre, E. Fast unfolding of communities in large networks. J. Stat. Mech. Theory Exp. 2008, 2008, P10008. [Google Scholar] [CrossRef] [Green Version]
  63. Raghavan, U.N.; Albert, R.; Kumara, S. Near linear time algorithm to detect community structures in large-scale networks. Phys. Rev. E 2007, 76, 36106. [Google Scholar] [CrossRef] [Green Version]
  64. Pons, P.; Latapy, M. Computing Communities in Large Networks Using Random Walks BT—Computer and Information Sciences—ISCIS 2005; Yolum, P., Güngör, T., Gürgen, F., Özturan, C., Eds.; Springer: Berlin/Heidelberg, Germany, 2005; pp. 284–293. [Google Scholar]
  65. Clauset, A.; Newman, M.E.J.; Moore, C. Finding community structure in very large networks. Phys. Rev. E 2004, 70, 66111. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Yu, G.; Wang, L.-G.; Han, Y.; He, Q.-Y. clusterProfiler: An R Package for Comparing Biological Themes among Gene Clusters. OMICS A J. Integr. Biol. 2012, 16, 284–287. [Google Scholar] [CrossRef]
Figure 1. Enrichment network for results obtained from the hypergeometric enrichment test: linkages of genes, GO biological process, and enriched terms are connecting overlapping gene sets.
Figure 1. Enrichment network for results obtained from the hypergeometric enrichment test: linkages of genes, GO biological process, and enriched terms are connecting overlapping gene sets.
Ncrna 08 00033 g001
Figure 2. Gene-concept networks from the KEGG enrichment analysis; the beige circles represent pathways while the gray ovals indicate genes.
Figure 2. Gene-concept networks from the KEGG enrichment analysis; the beige circles represent pathways while the gray ovals indicate genes.
Ncrna 08 00033 g002
Figure 3. Flowchart of the pipeline employed for the multi-cancer exosomal miRNA study.
Figure 3. Flowchart of the pipeline employed for the multi-cancer exosomal miRNA study.
Ncrna 08 00033 g003
Table 1. Results summary of exosomal miRNA expression in seven cancer datasets (patients vs. controls).
Table 1. Results summary of exosomal miRNA expression in seven cancer datasets (patients vs. controls).
CancerNumber of Upregulated miRNAsNumber of Downregulated miRNAs
Prostate cancer625
Gastric cancer7-
Colon cancer39-
Glioblastoma 21
Multiple myeloma10997
Lung cancer310
Liver cancer1946
Table 2. Identification of common miRNAs (patients vs. controls).
Table 2. Identification of common miRNAs (patients vs. controls).
CancerUpregulated miRNA (Fold Change)Downregulated miRNA (Fold Change)
Prostate cancermiR-543 + (1.16 × 10−6)miR-543 (−5.07), miR-495-3p (−4.25)
Gastric cancermiR-495-3p (1.86), miR-543 (1.47)
miR-543 + (1.35 × 10−10)
Colon cancermiR-495-3p (4.62), miR-543 (5.7)
miR-495-3p + (3.72 × 10−12)
Glioblastoma miR-543 (−4.64)
miR-495-3p–(0.0120)
Multiple myelomamiR-495-3p (1.32)
Table 3. Kyoto Encyclopedia of Genes and Genomes pathways enriched for mRNAs involved in the miR-495-3p community.
Table 3. Kyoto Encyclopedia of Genes and Genomes pathways enriched for mRNAs involved in the miR-495-3p community.
DescriptionAdjusted p-ValueGene ID
GABAergic synapse0.001GABRA1/NSF/GABRA4/GLS/SLC12A5/GAD1/SLC38A2/GAD2/GABBR2/GABRB2/CACNA1D/GABRG1/GABRB3/PRKCB/GABARAP
AMPK signaling pathway0.002SREBF1/RAB10/EEF2K/ELAVL1/PPP2R2C/PPP2R5C/LIPE/SCD/RHEB/PRKAG2/AKT3/RAB14/PFKFB2/MAP3K7/SCD5/PIK3R1/IRS2
MAPK signaling pathway0.020HSPA1B/LDLR/TGFBR2/RNF41/RAB10/RAB31/TFRC/EPS15/SNX5/EEA1/RAB11A/RAB5B/ACAP2/SNX4/STAM2/CHMP3/CAPZA1/DNAJC6/RAB22A/HSPA2/PSD3/KIF5B/DNM3/CAPZA2
Endocytosis0.0204017HSPA1B/LDLR/TGFBR2/RNF41/RAB10/RAB31/TFRC/EPS15/SNX5/EEA1/RAB11A/RAB5B/ACAP2/SNX4/STAM2/CHMP3/CAPZA1/DNAJC6/RAB22A/HSPA2/PSD3/KIF5B/DNM3/CAPZA2
Oxytocin signaling pathway0.0204017EEF2K/CAMK2A/PTGS2/MEF2C/PRKAG2/ITPR2/CAMK2G/GUCY1A2/CACNA2D1/CACNA1D/RCAN1/CALM3/PRKCB/CAMK1D/ACTG1/CALM1/PPP1CB
Table 4. Summary of miRNA datasets among different cancer types for this study.
Table 4. Summary of miRNA datasets among different cancer types for this study.
Accession NumberSample TypePatient FeaturesMethodological Analysis
GSE130654
[51]
Extracellular vesicles derived from gastric cancer patients36 non-cardia adenocarcinoma patients (stages I and II) and 12 healthy individualsDEseq2
GSE111803
[52]
Extracellular vesicles derived from lung cancer patients5 patients with lung adenocarcinoma and 5 healthy controlsDEseq2
GSE94564
[43]
Extracellular vesicles derived from multiple myeloma patients10 patients newly diagnosed with MM and 5 healthy individualsDEseq2
GSE123972
[53]
Extracellular vesicles derived from Hepatocellular carcinoma patients10 individuals with HCC pooled into 2 libraries and 10 healthy donors pooled into 2 librariesDEseq2
GSE71008
[54]
Extracellular vesicles derived from colon and
prostate cancer patients
100 colon cancer patients and 36 prostate cancer patients, and 50 healthy controlsSAM
GSE122488
[39]
Extracellular vesicles derived from glioblastoma patients12 patients with glioblastoma,
10 with glioma stages II–III,
and 16 healthy controls
SAM
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rincón-Riveros, A.; Rodríguez, J.A.; Villegas, V.E.; López-Kleine, L. Identification of Two Exosomal miRNAs in Circulating Blood of Cancer Patients by Using Integrative Transcriptome and Network Analysis. Non-Coding RNA 2022, 8, 33. https://doi.org/10.3390/ncrna8030033

AMA Style

Rincón-Riveros A, Rodríguez JA, Villegas VE, López-Kleine L. Identification of Two Exosomal miRNAs in Circulating Blood of Cancer Patients by Using Integrative Transcriptome and Network Analysis. Non-Coding RNA. 2022; 8(3):33. https://doi.org/10.3390/ncrna8030033

Chicago/Turabian Style

Rincón-Riveros, Andrés, Josefa Antonia Rodríguez, Victoria E. Villegas, and Liliana López-Kleine. 2022. "Identification of Two Exosomal miRNAs in Circulating Blood of Cancer Patients by Using Integrative Transcriptome and Network Analysis" Non-Coding RNA 8, no. 3: 33. https://doi.org/10.3390/ncrna8030033

APA Style

Rincón-Riveros, A., Rodríguez, J. A., Villegas, V. E., & López-Kleine, L. (2022). Identification of Two Exosomal miRNAs in Circulating Blood of Cancer Patients by Using Integrative Transcriptome and Network Analysis. Non-Coding RNA, 8(3), 33. https://doi.org/10.3390/ncrna8030033

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop