Abstract
The human Periodontal Ligament Stem Cells (hPDLSCs) exhibit self-renewal capacity and clonogenicity potential. The Extracellular Vesicles (EVs) secreted by hPDLSCs are particles containing lipids, proteins, mRNAs, and non-coding RNAs, among which microRNAs, that are important in intercellular communication. The purpose of this study was the analysis of the non-coding RNAs contained in the EVs derived from hPDLSCs using Next Generation Sequencing. Moreover, our data were enriched using bioinformatic tools. The analysis highlighted the presence of non-coding RNAs and five microRNAs: MIR24-2, MIR142, MIR335, MIR490, and MIR296. Our results show that these miRNAs target the genes classified in two terms of the Gene Ontology: “Ras protein signal transduction” and “Actin/microtubule cytoskeleton organization.” Noteworthy, the in-deep analysis of our EVs highlights that the miRNAs could be implicated in the silencing of proto-oncogenes involved in 12 different types of tumors.
Introduction
The human Periodontal Ligament Stem Cells (hPDLSCs) are mesenchymal stem cells that can be easily harvested from periodontal tissue. Non-invasive surgery during standard dental scaling does not entail any additional risk to the donor (). The hPDLSCs showed self-renewal capacity, differentiation and immunomodulatory proprieties (). The Extracellular Vesicles (EVs) are particles delimited by a lipid bilayer capable of crossing biological barriers and being internalized in the target cells so that they play an essential role in intercellular communication (Valadi et al., 2007; Pizzicannella et al., 2019b). EVs contain different macromolecules such as lipids, proteins, DNA, mRNA and non-coding RNAs (ncRNAs) among which microRNAs (or miRNAs) (Zhang et al., 2019). The dynamic light-scattering analysis, previously demonstrated from our research group, highlighted the presence of two different populations of EVs derived from hPDLSCs. The average diameter of the EVs populations is between 90 ± 20 nm and 1,200 ± 400 nm and they have the same ζ-potential of −10.7 ± 0.9 mV. The tapping-mode topographic 3D atomic force-microscopy measurements showed that EVs derived from hPDLSCs have a globular shape with a central depression and a relatively smooth surface ().
The ncRNAs are circulating RNAs that represent the 99% of total RNAs and can be classified according with length (small 18–200 nucleotides and long >200), with function (ribosomal RNAs and transfer RNAs), or with regulation [miRNAs and long non-coding RNAs (lncRNAs)] (; ). The ncRNAs are involved in gene regulation by RNA interference, RNA modification or spliceosomal cycle. Specifically, miRNAs are small ncRNAs (18–25 nucleotides in length) that prevent the translation of mRNAs and consequently contribute to define the amount of proteins inside the cell (). miRNAs play an important role in physiological functions such as cell-cell communication, cell proliferation and vascularization. The ncRNAs and miRNAs are implicated in various diseases including cancer, in which they play an important biological role (Slaby et al., 2017).
Therefore, in this study, using Next Generation Sequencing (NGS), we want to analyze the ncRNAs and miRNAs contained in EVs derived from hPDLSCs in the early stages of stemness. Particularly, we investigated the possible implications of ncRNAs in the biological processes and in the proto-oncogenes modulation.
Materials and Methods
Culture and Extraction of the hPDLSCs
All the subjects that were involved in the study gave their informed consent. The protocol executed for this study was approved by the Medical Ethics Committee at the Medical School, “G. d’Annunzio” University of Chieti-Pescara, Chieti, Italy (n°266 17 April 2014) and it is in accordance with the Declaration of Helsinki. Five patients in healthy general conditions were chosen for tooth removal for orthodontic purposes. Subsequently, the cells were cultured in the Mesenchymal Stem Cell Growth Medium Chemically Defined (MSCGM-CD) Bulletin medium (Lonza, Basel, Switzerland). In order to minimize exposure to the non-human substances and facilitate the growth of human MSCs, the MSCGM-CD was changed twice a week. Cells isolation and characterization were performed as previously described (Rajan et al., 2016; Trubiani et al., 2016). The hPDLSCs were collected and washed with PBS (Lonza) several times. Finally, they were cultured at 37°C with 5% of CO2 with the MSCGM-CD medium.
Isolation of the hPDLSCs-Derived EVs
The hPDLSCs, isolated from each patient, were collected from conditioned medium (CM; 10 mL) at passage 2 after 48 h of incubation. For 15 min, the CM was centrifuged at 3,000 × g in order to remove debris and suspended cells. EVs were extracted using the commercial agglutinant ExoQuick TC (System Biosciences, Euroclone SpA, Milan, Italy). Specifically, 10 mL of CM of the hPDLSCs were mixed with 2 mL ExoQuick TC and the whole was incubated at 4°C overnight without rotation. Subsequently, in order to deposit the EVs, a centrifugation step was carried out for 30 min at 1,500 × g. Finally, the pellets were resuspended in 200 μL of PBS (Ca2+ and Mg2+). The EVs, obtained from hPDLSCs of each patient, were used for transcriptomic analyses.
RNA Extraction and Non-coding Analysis
The RNA extraction and processing were conducted as previously reported (Silvestro et al., 2020). Briefly, the Total Exosome RNA and Protein Isolation Kit (catalog #4478545; Thermo Fisher Scientific, Rockford, IL, United States) was used to isolate the RNA following the manufacturer’s protocol and 30 μL of RNA solution were collected from each sample (Pizzicannella et al., 2019a).
The manufacturer instructions and the TruSeq RNA Exome protocol (Illumina, San Diego, CA, United States) were followed in order to prepare the library. The normalized libraries were analyzed using the MiSeq instrument (Illumina) in a single read mode. For the NGS analysis, the RNA extracted from the EVs of each patient was repeated in triplicate.
The software fastQC (Babraham Institute, Cambridge, United Kingdom) was used to perform the quality check of the reads. Then, the adapters and the low-quality bases were identified and removed by the software Trimmomatic (Usadel Lab, Aachen, Germany) (). The “Homo Sapiens” reference genome, available from the University of California Santa Cruz website1, was used to align the reads taking advantage of the software Spliced Transcripts Alignment to a Reference RNA-seq aligner (STAR) (). Finally, the Cufflinks software (Trapnell Lab, Washington, DC, United States), version 2.2.1, was used to assign the gene symbol to each transcript (Trapnell et al., 2013). In order to study the transcripts, the whole set was classified with HUGO Gene Nomenclature Committee website2 (Stelzer et al., 2016) using the “non-coding RNAs” group (475). Specifically, the genes that are targeted by the identified miRNAs were analyzed with TargetScanHuman (version 7.2) (). Moreover, the target genes were enriched with the “Biological Process” terms included in the “Gene Ontology” using PANTHER. Finally, the Human microRNA Disease Database (HMDD) v3.2 was used in order to associate the genes target of our miRNAs to the disease in which their expression is deregulated. All the plots were depicted with the software R.
Results
Evaluation of Non-coding Transcripts of EVs Derived From hPDLSCs
The analysis of EVs derived from hPDLSCs reveals 955 non-coding transcripts (Supplementary Table S1). The HUGO website characterizes 212 transcripts distributed in 11 groups as showed in Figure 1. The most of the non-coding transcripts of EVs belong to the “Antisense RNAs” (109), to the “Long intergenic non-protein coding RNAs” (40) and to the “Long non-coding RNAs with non-systematic symbols” (25). The remaining transcripts are included in smaller categories. Our EVs contain also 5 miRNAs: MIR24-2, MIR142, MIR296, MIR335, MIR490. We studied the biological roles of our miRNAs using TargetScanHuman () that associates each miRNA to its target genes. We used the database PANTHER (Thomas et al., 2003) to visualize all the “Biological Process” terms in which the target genes are involved (Figure 2). The deepest terms that better characterize our transcripts are related to the cytoskeletal organization (“Actin/microtubule cytoskeleton organization”) and the signal transduction mediated by Ras protein (“Ras protein signal transduction”). PANTHER was also used to associate the genes with the pathway in which they are involved. In particular, the most characterized pathways include the genes CDC42, RAC1, RHOC, RHOA, RHOJ, SPRY3, PFN2, WASL, TUBB, TUBB4B, SSH3, SSH1, ARPC5L, PARVA, RHOQ. The first five genes are in common to “Actin/microtubule cytoskeleton organization” and “Ras protein signal transduction.” SPRY3 is found only in the pathway related to Ras protein while the other genes are exclusively related to the cytoskeletal organization. Moreover, we inspected the amount of genes that are simultaneously targeted by our miRNAs as represented in the Venn Diagram in Figure 3. There is no gene that is targeted by all the miRNAs. Nevertheless, we observed that one gene (GABRB1) is not targeted just by MIR142, one gene (ONECUT2) just by MIR24, one gene (CNNM2) just by MIR490. Furthermore, MIR-24 targets the most of the genes (1,500) immediately followed by MIR142 (1,250). MIR490, MIR335, and MIR296 target exclusively less than 200 genes each. Furthermore, using HMDD () we reported in Table 1 the associations found among our miRNAs and the target genes deregulated in specific diseases.
FIGURE 1
FIGURE 2
FIGURE 3
TABLE 1
| MiRNA | Target gene | Gene name | Disease | References |
| MIR24-2 | S100A8 | S100 calcium binding protein A8 | Laryngeal carcinoma | |
| MIR24-2 | AGPAT2 | 1-acylglycerol-3-phosphate O-acyltransferase 2 | Osteosarcoma | Song et al., 2013 |
| MIR142 | RAC1 | Rac family small GTPase 1 | Hepatocellular carcinoma | Wu et al., 2011 |
| MIR142 | PROM1 | Prominin 1 | Colon cancer | Shen et al., 2013 |
| MIR142 | ABCG2 | ATP binding cassette subfamily G member 2 | Colon cancer | Shen et al., 2013 |
| MIR142 | LGR5 | Leucine rich repeat containing G protein-coupled receptor 5 | Colon cancer | Shen et al., 2013 |
| MIR142 | HMGA1 | High mobility group AT-hook 1 | Osteosarcoma | Xu et al., 2014 |
| MIR296 | HMGA1 | High mobility group AT-hook 1 | Prostate cancer | Wei et al., 2011 |
| MIR296 | PIN1 | Peptidylprolyl cis/trans isomerase, NIMA-interacting 1 | Prostate cancer | |
| MIR296 | FGFR1 | Fibroblast growth factor receptor 1 | Hepatocellular carcinoma | Wang et al., 2016 |
| MIR296 | AKT2 | AKT serine/threonine kinase 2 | Pancreatic cancer | |
| MIR296 | PLK1 | Polo like kinase 1 | Lung cancer | Xu et al., 2016 |
| MIR296 | SP1 | Sp1 transcription factor | Cervical cancer | |
| MIR296 | CDX1 | Caudal type homeobox 1 | Gastric cancer | |
| MIR490 | CCND1 | Cyclin D1 | Lung cancer | |
| MIR490 | CDK1 | Cyclin dependent kinase 1 | Ovarian cancer | |
| MIR490 | FOS | Fos proto-oncogene, AP-1 transcription factor subunit | Bladder cancer | |
| MIR490 | HMGA2 | High mobility group AT-hook 2 | Osteosarcoma | |
| MIR335 | BCL2L2 | BCL2 like 2 | Gastric cancer | Xu et al., 2012 |
| Ovarian cancer | ||||
| MIR335 | SP1 | Sp1 transcription factor | Gastric cancer | Xu et al., 2012 |
| MIR335 | ROCK1 | Rho associated coiled-coil containing protein kinase 1 | Osteosarcoma | Wang et al., 2013 |
| Neuroblastoma | ||||
| MIR335 | MAPK1 | Mitogen-activated protein kinase 1 | Neuroblastoma | |
| MIR335 | LRG1 | Leucine rich alpha-2-glycoprotein 1 | Neuroblastoma |
MiRNAs found in hPDLSCs-derived EVs associated to specific cancer.
To each miRNA found in hPDLSC-derived EVs is associated a Target Gene and the kind of cancer in which it is found deregulated inside the Human MicroRNA Disease Database along with the Reference in which it is described.
Discussion
EVs are phospholipids membrane-enclosed organizations that carry proteins, lipids and nucleic acids including mRNAs and ncRNAs (Trubiani et al., 2019).
In this study, we analyzed the ncRNAs contained in EVs derived from hPDLSCs at early stage of stemness with a focus on miRNAs.
Our results show that the ncRNAs contained in EVs are classified in 11 families (Figure 1). The most represented class of ncRNAs is the “Antisense RNAs,” small molecules that regulate gene’s expression by binding the complementary mRNAs (Xu et al., 2018). Our results show also the presence of lncRNAs, an important class of transcripts that indirectly regulates the transcription recruiting transcription factors and affecting mRNA stability (). EVs transport lncRNAs from one cell to another and can induce epigenetic modifications in the receiving cells (). Other classes of ncRNAs represented in our EVs are the small nucleolar RNAs and miRNAs. The small nucleolar RNAs are involved in the chemical alteration of different RNAs and in the regulation of the alternative splicing in pre-mRNAs (; ; Sharma et al., 2016).
Among ncRNAs, miRNAs are the most studied since they can regulate the expression of 60% of the human genes. Specifically, miRNAs are post-transcriptional regulators that bind complementary mRNA and consequently reduce protein expression. miRNAs, released from the EVs, can determine the cell fate through their participation in differentiation and in the regulation of reprogramming processes (; ; ). In the last two decades, many studies contextualized the role of miRNAs in cancer development (Rupaimoole and Slack, 2017). Interestingly, in our EVs, we found five miRNAs: MIR24-2, MIR142, MIR335, MIR490, and MIR296.
Our bioinformatic analysis, represented in Figure 3, highlights the amount of target genes that are modulated by our miRNAs. Indeed, both MIR24-2 and MIR142 are able to influence more than 1,000 genes while MIR335, MIR490, and MIR296 less than 200 each. The representation with the Venn Diagram shows that each miRNA can regulate several genes and each gene can be regulated by more than one miRNA. This evidence suggests that MIR24-2 and MIR142 are ones most involved.
The MIR24-2, MIR142, MIR335, MIR490, and MIR296 found in our EVs target the genes mainly involved in “Ras protein signal transduction” and “Actin/microtubule cytoskeleton organization” (Figure 2).
The Ras proteins are GTPases that act as switches recruiting either the Map Kinases or the RHO/RAC family resulting in downstream signaling for transcription, cell cycle or cytoskeletal organization ().
The analysis performed using PANTHER database shows that our miRNAs target CDC42, RAC1, RHOA, RHOC, and RHOJ that simultaneously participate in the regulation of “Actin/microtubule cytoskeleton organization” and “Ras protein signal transduction.” These genes encode for smalls GTPases/Rho family that are active when bound to Guanosine Triphosphate (GTP) or inactive when bound to Guanosine Diphosphate (GDP). The CDC42, RAC1, RHOC, and RHOA genes play an important role in cytoskeletal rearrangement, morphology regulation, cell motility, cell adhesion and cell cycle regulation (; Thomas et al., 2019). RHOJ is a gene that encodes for an endothelial member of the Cdc42 subfamily involved in endothelial motility, tubulogenesis and microtubule lumen formation (). Indeed, miRNAs are very important regulators of these GTPase Rho family ().
Additionally, as regards the Gene Ontology term “Ras protein signal transduction,” the miRNAs inside our EVs target the SPRY3 gene, which encodes for the Protein sprouty homolog 3 (Spry3), an intracellular negative regulators of Receptor tyrosine kinase signaling and of the Fibroblast Growth Factor (FGF) pathways (). Furthermore, showed that Spry3 negatively regulates the release of calcium downstream of brain-derived neurotrophic factor signaling, and consequently cause the inhibition of the axons branching in cultured cortical neurons.
Moreover, the miRNAs present in our EVs regulate the “Actin/microtubule cytoskeleton organization” process by targeting different genes such as TUBB, TUBB4B, PFN2, WASL, SSH3, SSH1, ARPC5L, PARVA, and RHOQ. In detail, TUBB and TUBB4B encode respectively for a beta-tubulin protein and for the beta-4B tubulin chain, both representing structural components of the microtubules. RHOQ encodes for a small GTPase involved in the regulation of the cell shape through the organization of the actin cytoskeleton. SSH1 and SSH3 encode for protein phosphatases that play a role in the regulation of the dynamics of the actin filament. In addition, SSH1 and SSH3 activate the actin-binding protein family ADF/cofilin that promotes the disassembling of the actin filaments (). ARPC5L gene, part of Arp2/3 complex, encodes for the Actin Related Protein 2/3 Complex Subunit 5 Like. This gene regulates the actin polymerization in the cytoskeleton and generates an actin filament network (). PARVA encodes for the α-parvin protein that is involved in the regulation of the actin cytoskeleton dynamism (Sepulveda and Wu, 2006). PFN2 is an actin-binding protein that regulates the architecture of the synapse and the cytoskeleton (). WASL, which encodes for Neural Wiskott-Aldrich syndrome protein, appears to be involved in the polymerization of actin. In addition, PFN2 and WASL are highly expressed in neural tissues and participate to the axon guidance in the Slit/Robo pathway. These results show that our EVs contain miRNAs involved in axon guidance and consequently in neural differentiation.
Moreover, we evaluated (using HMDD) how much our miRNAs deregulate oncogenes involved in different kinds of cancers. Noteworthy, our miRNAs are involved at least in 12 different cancers (Table 1) and specifically in their development and progression ().
The MIR24-2 encodes miR-24 that plays an important role as a suppressor of the genes E2F2 and MYC. Moreover, by binding the complementary 3′-UTR mRNA of these genes, miR24-2 is involved in the regulation of the cell cycle (). also showed that the up-regulation of MIR24 leads to morphological changes, low cell proliferation and enhancement of cell invasion potential in laryngeal squamous cell carcinoma by inhibiting transcription of the S100A8 gene. Additionally, the overexpression of MIR24 inhibits osteosarcoma cell proliferation by blocking the transcription of AGPAT2 gene encoding for Lysophosphatidic Acid Acyltransferase β, an enzyme involved in osteosarcoma cell proliferation (Song et al., 2013).
The MIR142 appears to be a potential proto-oncogene suppressor. Indeed, through the transcriptional block of RAC1, a gene that encodes for a GTPase that regulates cell growth and migration and promotes the activation of protein kinases. Wu et al. (2011) showed that MIR142 has also a role in suppressing migration and invasion of hepatocellular carcinoma cells. In addition, MIR142 inhibits the proliferation of pancreatic cancer cells by decreasing the expression of heat shock protein 70 (). It inhibits the growth of cells responsible for colon cancer, hindering the transcription of PROM1 (CD133), ABCG2, and LRG5 genes (Shen et al., 2013). Moreover, by blocking the transcription of the HMGA1 gene, an important gene involved in promoting cancer and increasing invasiveness, Xu et al. (2014) reveals that MIR142 can act as a suppressor in osteosarcoma. Furthermore, as reported by the MIR142 is secreted by oral squamous cell carcinoma cells, promoting the growth of the tumor. The release of miR-142-3p also influences the tumor microenvironment, promoting angiogenesis (). Indeed, in oral squamous cell carcinoma, MIR142 was up-regulated, demonstrating how the differential expression of miRNAs could be useful in the diagnosis of oral cancer (Xu et al., 2019).
HMGA1 gene can be also targeted by MIR296, another miRNA presents in our EVs, supporting the knowledge by which different miRNAs can interact with the same mRNA. Negative regulation of this gene by MIR296 appears to reduce invasiveness in prostate cancer cells (Wei et al., 2011). In addition, MIR296 targeting PIN1, another gene involved in tumor development, suppresses cell proliferation and growth in prostate cancer cells (). In hepatocellular carcinoma, MIR296 inhibits the proliferative capacity and progression of the cell cycle and induces apoptosis blocking the transcription of the FGFR1 proto-oncogene (Wang et al., 2016). Other targets of MIR296 are the proto-oncogenes AKT2 in pancreatic cancer (), PLK1 in non-small cell lung cancer (Xu et al., 2016), SP1 in cervical cancer () and CDX1 in gastric cancer (). Through the negative regulation of the expression of these genes, MIR296 exerts its action as a tumor suppressor. Therefore, EVs derived from hPDLSCs, through the secretion of MIR296, could act as suppressors of the malignant progression by attenuating the transcription of oncogenic targets. The MIR296, as reported by Severino et al. (2015) was also detected in the in metastases of patients with oral squamous cell carcinoma.
In the transcriptome of our EVs, we found also MIR490 that seems to have the role of tumor suppressor in lung cancer cells A549 and ovarian cancer by blocking the transcription of CCND1 and CDK1 respectively, two important genes involved in cell cycle progression (; ). Additionally, as showed by , the MIR490 seems to be a suppressor of bladder cancer cell proliferation by blocking the transcription of FOS. Moreover, MIR490 inhibits the expression of HMGA2 and its downregulation appears to affect the development potential of the osteosarcoma ().
MIR335, present in our EVs, also appears to play the role of proto-oncogene suppressor in a wide variety of cancers. It is known that MIR335 negatively regulates the metastasis in gastric cancer by targeting BCL2L2 and SP1 (Xu et al., 2012). Likewise targeting the same BCL2L2 gene, it is capable of inhibiting invasiveness in ovarian cancer (). In neuroblastoma cells, MIR335 downregulates the ROCK1 and MAPK1 genes, involved in the non-canonical TGF-β pathways. This modulation implicates a reduced invasiveness of the neuroblastoma cells. In addition, MIR335 inhibits the mRNA transcribed by LRG1 gene reducing the migration of neuroblastoma cells (). Similarly, it negatively regulates at the post-transcriptional level the ROCK1 inhibiting tumor cell invasion and migration in osteosarcoma (Wang et al., 2013). Conversely, showed that MIR335 is up-regulated in fibroblasts associated with senescent cancer. These cells present in the tumor environment of oral neoplasms. Furthermore, the authors reported that MIR335 is also involved in the development of a secretory phenotype associated with senescence. In this way, MIR335 contributing to the progression of cancer ().
Interestingly, in our previous manuscript, we have already demonstrated in an analysis of the oral stem cell transcriptome that, following the in vitro expansion, many oncogenes were absent or poorly expressed ().
Conclusion
The EVs derived from hPDLSCs contain several ncRNAs among which the five miRNAs: MIR24-2, MIR142, MIR296, MIR335, and MIR490. Specifically, these miRNAs regulate genes that are involved in “Ras protein signal transduction” and “Actin/microtubule cytoskeleton organization” processes that regulate cell growth and differentiation during cytokinesis. Moreover, our results demonstrate that MIR24 and MIR142 are the most relevant miRNAs because they target more than 1,000 genes each. Furthermore, all the miRNAs detected in our EVs could have a potential role as proto-oncogenes suppressors. These findings indicated that the EVs obtained from hPDLSCs can be considered as anticancer therapeutic agents. This evidence supports the role of mesenchymal stem cell derivatives in the deregulatory functioning of human carcinogenesis in addition to their use in regenerative medicine.
Statements
Data availability statement
We uploaded our data to the repository Sequence Read Archive (accession number PRJNA630492).
Ethics statement
The study involving human participants was approved by the Medical Ethics Committee at the Medical School, “G. d’Annunzio” University of Chieti-Pescara, Chieti, Italy (n°266 17 April 2014) and it is in accordance with the Declaration of Helsinki. All the subjects that were involved in the study gave their informed consent.
Author contributions
JP and GM performed the hPDLSCs culturing and the isolation of the EVs. AG performed the RNA extraction and library preparation. LC performed the computational and statistical analysis. LC and SS wrote the manuscript. PB, OT, and EM designed the study and revised the manuscript.
Funding
This study was supported by a Current Research Fund 2019, Ministry of Health, Italy.
Acknowledgments
We would like to thank the Ministry of Health, Italy.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.00582/full#supplementary-material
Abbreviations
- EVs
Extracellular Vesicles
- FGF
Fibroblast Growth Factor
- GDP
Guanosine Diphosphate
- GTP
Guanosine Triphosphate
- hPDLSCs
human Periodontal Ligament Stem Cells
- lncRNAs
long non-coding RNAs
- miRNAs
microRNAs
- MSCGM-CD
Mesenchymal Stem Cell Growth Medium Chemically Defined
- ncRNAs
non-coding RNAs
- NGS
Next Generation Sequencing
- Spry3
protein sprout homolog 3
- STAR
Spliced Transcripts Alignment to a Reference RNA-seq aligner
- HMDD
Human microRNA Disease Database.
Footnotes
1.^http://labshare.cshl.edu/shares/gingeraslab/www-data/dobin/STAR/STARgenomes/Old/ENSEMBL/homo_sapiens/ENSEMBL.homo_sapiens.release-75/
References
1
AbellaJ. V.GalloniC.PernierJ.BarryD. J.KjaerS.CarlierM. F.et al (2016). Isoform diversity in the Arp2/3 complex determines actin filament dynamics.Nat. Cell Biol.1876–86. 10.1038/ncb3286
2
AgarwalV.BellG. W.NamJ. W.BartelD. P. (2015). Predicting effective microRNA target sites in mammalian mRNAs.eLife4:e05005.
3
AspenströmP. (2019). The intrinsic gdp/gtp exchange activities of Cdc42 and rac1 are critical determinants for their specific effects on mobilization of the actin filament system.Cells8:759. 10.3390/cells8070759
4
BartelD. P. (2004). MicroRNAs: genomics, biogenesis, mechanism, and function.Cell116281–297.
5
BolgerA. M.LohseM.UsadelB. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data.Bioinformatics302114–2120. 10.1093/bioinformatics/btu170
6
BusteloX. R.SauzeauV.BerenjenoI. M. (2007). GTP−binding proteins of the Rho/Rac family: regulation, effectors and functions in vivo.Bioessays29356–370. 10.1002/bies.20558
7
CalinG. A.CroceC. M. (2006). MicroRNA signatures in human cancers.Nat. Rev. Cancer6857–866. 10.1038/nrc1997
8
CaoJ.CaiJ.HuangD.HanQ.YangQ.LiT.et al (2013). miR-335 represents an invasion suppressor gene in ovarian cancer by targeting Bcl-w.Oncol. Rep.30701–706. 10.3892/or.2013.2482
9
ChenS.ChenX.XiuY.-L.SunK.-X.ZhaoY. (2015). MicroRNA-490-3P targets CDK1 and inhibits ovarian epithelial carcinoma tumorigenesis and progression.Cancer Lett.362122–130. 10.1016/j.canlet.2015.03.029
10
ChiricostaL.GugliandoloA.DiomedeF.PizzicannellaJ.TrubianiO.IoriR.et al (2019). Moringin pretreatment inhibits the expression of genes involved in mitophagy in the stem cell of the human periodontal ligament.Molecules24:3217. 10.3390/molecules24183217
11
DecaturW. A.FournierM. J. (2002). rRNA modifications and ribosome function.Trends Biochem. Sci.27344–351. 10.1016/s0968-0004(02)02109-6
12
DickmanC. T.LawsonJ.JabaleeJ.MaclellanS. A.LepardN. E.BennewithK. L.et al (2017). Selective extracellular vesicle exclusion of miR-142-3p by oral cancer cells promotes both internal and extracellular malignant phenotypes.Oncotarget815252–15266. 10.18632/oncotarget.14862
13
DiomedeF.D’auroraM.GugliandoloA.MerciaroI.EttorreV.BramantiA.et al (2018). A novel role in skeletal segment regeneration of extracellular vesicles released from periodontal-ligament stem cells.Int. J. Nanomed.13:3805. 10.2147/IJN.S162836
14
DobinA.DavisC. A.SchlesingerF.DrenkowJ.ZaleskiC.JhaS.et al (2013). STAR: ultrafast universal RNA-seq aligner.Bioinformatics2915–21. 10.1093/bioinformatics/bts635
15
DozmorovM. G.GilesC. B.KoelschK. A.WrenJ. D. (2013). Systematic classification of non-coding RNAs by epigenomic similarity.BMC Bioinformatics14(Suppl. 14):S2. 10.1186/1471-2105-14-S14-S2
16
EleuterioE.TrubianiO.SulpizioM.Di GiuseppeF.PierdomenicoL.MarchisioM.et al (2013). Proteome of human stem cells from periodontal ligament and dental pulp.PLoS One8:e71101. 10.1371/journal.pone.0071101
17
FalaleevaM.PagesA.MatuszekZ.HidmiS.Agranat-TamirL.KorotkovK.et al (2016). Dual function of C/D box small nucleolar RNAs in rRNA modification and alternative pre-mRNA splicing.Proc. Natl. Acad. Sci. U.S.A.113E1625–E1634. 10.1073/pnas.1519292113
18
FatimaF.NawazM. (2017). Vesiculated long non-coding RNAs: offshore packages deciphering trans-regulation between cells, cancer progression and resistance to therapies.Non Coding RNA3:10. 10.3390/ncrna3010010
19
FriedmanR. C.FarhK. K.BurgeC. B.BartelD. P. (2009). Most mammalian mRNAs are conserved targets of microRNAs.Genome Res.1992–105. 10.1101/gr.082701.108
20
GuH.YangT.FuS.ChenX.GuoL.NiY. (2014). MicroRNA-490-3p inhibits proliferation of A549 lung cancer cells by targeting CCND1.Biochem. Biophys. Res. Commun.444104–108. 10.1016/j.bbrc.2014.01.020
21
GugliandoloA.RajanT. S.SciontiD.DiomedeF.BramantiP.MazzonE.et al (2017). Reprogramming of oncogene expression in gingival mesenchymal stem cells following long-term culture in vitro.Cell Reprogram19159–170. 10.1089/cell.2016.0056
22
GuoY.FuW.ChenH.ShangC.ZhongM. (2012). miR-24 functions as a tumor suppressor in Hep2 laryngeal carcinoma cells partly through down-regulation of the S100A8 protein.Oncol. Rep.271097–1103. 10.3892/or.2011.1571
23
HacohenN.KramerS.SutherlandD.HiromiY.KrasnowM. A. (1998). sprouty encodes a novel antagonist of FGF signaling that patterns apical branching of the Drosophila airways.Cell92253–263. 10.1016/s0092-8674(00)80919-8
24
HuangZ.ShiJ.GaoY.CuiC.ZhangS.LiJ.et al (2019). HMDD v3.0: a database for experimentally supported human microRNA-disease associations.Nucleic Acids Res.47D1013–D1017. 10.1093/nar/gky1010
25
JeongD. H.ChoiY. N.SeoT. W.LeeJ. S.YooS. J. (2018). Ubiquitin-proteasome dependent regulation of Profilin2 (Pfn2) by a cellular inhibitor of apoptotic protein 1 (cIAP1).Biochem. Biophys. Res. Commun.506423–428. 10.1016/j.bbrc.2018.10.115
26
KabirT. D.LeighR. J.TasenaH.MelloneM.ColettaR. D.ParkinsonE. K.et al (2016). A miR-335/COX-2/PTEN axis regulates the secretory phenotype of senescent cancer-associated fibroblasts.Aging (Albany NY)81608–1635. 10.18632/aging.100987
27
KuritaS.GunjiE.OhashiK.MizunoK. (2007). Actin filaments-stabilizing and -bundling activities of cofilin-phosphatase Slingshot-1.Genes Cells12663–676. 10.1111/j.1365-2443.2007.01078.x
28
LalA.NavarroF.MaherC. A.MaliszewskiL. E.YanN.O’dayE.et al (2009). miR-24 Inhibits cell proliferation by targeting E2F2, MYC, and other cell-cycle genes via binding to “seedless” 3’ UTR microRNA recognition elements.Molecular Cell35610–625. 10.1016/j.molcel.2009.08.020
29
LeeK.-H.LinF.-C.HsuT.-I.LinJ.-T.GuoJ.-H.TsaiC.-H.et al (2014). MicroRNA-296-5p (miR-296-5p) functions as a tumor suppressor in prostate cancer by directly targeting Pin1.Biochim. Biophys. Acta Mol. Cell Res.18432055–2066. 10.1016/j.bbamcr.2014.06.001
30
LeszczynskaK.KaurS.WilsonE.BicknellR.HeathV. L. (2011). The role of RhoJ in endothelial cell biology and angiogenesis.Biochem. Soc. Trans.391606–1611. 10.1042/bst20110702
31
LiH.LiJ.ShiB.ChenF. (2017). MicroRNA-296 targets AKT2 in pancreatic cancer and functions as a potential tumor suppressor.Mol. Med. Rep.16466–472. 10.3892/mmr.2017.6602
32
LiS.XuX.XuX.HuZ.WuJ.ZhuY.et al (2013). MicroRNA-490-5p inhibits proliferation of bladder cancer by targeting c-Fos.Biochem. Biophys. Res. Commun.441976–981. 10.1016/j.bbrc.2013.11.006
33
LiT.LuY.ZhaoX.GuoH.LiuC.LiH.et al (2014). MicroRNA-296-5p increases proliferation in gastric cancer through repression of Caudal-related homeobox 1.Oncogene33783–793. 10.1038/onc.2012.637
34
LiuM.BiF.ZhouX.ZhengY. (2012). Rho GTPase regulation by miRNAs and covalent modifications.Trends Cell Biol.22365–373. 10.1016/j.tcb.2012.04.004
35
LiuW.XuG.LiuH.LiT. (2015). MicroRNA−490−3p regulates cell proliferation and apoptosis by targeting HMGA2 in osteosarcoma.FEBS Lett.5893148–3153. 10.1016/j.febslet.2015.08.034
36
LvL.WangX. (2018). MicroRNA-296 targets specificity protein 1 to suppress cell proliferation and invasion in cervical cancer.Oncol. Res. Featuring Preclin. Clin. Cancer Ther.26775–783. 10.3727/096504017X15132494420120
37
LynchJ.FayJ.MeehanM.BryanK.WattersK. M.MurphyD. M.et al (2012). MiRNA-335 suppresses neuroblastoma cell invasiveness by direct targeting of multiple genes from the non-canonical TGF-β signalling pathway.Carcinogenesis33976–985. 10.1093/carcin/bgs114
38
MackenzieT. N.MujumdarN.BanerjeeS.SangwanV.SarverA.VickersS.et al (2013). Triptolide induces the expression of miR-142-3p: a negative regulator of heat shock protein 70 and pancreatic cancer cell proliferation.Mol. Cancer Ther.121266–1275. 10.1158/1535-7163.MCT-12-1231
39
MensM. M. J.GhanbariM. (2018). Cell cycle regulation of stem cells by MicroRNAs.Stem Cell Rev. Rep.14309–322. 10.1007/s12015-018-9808-y
40
PalazzoA. F.LeeE. S. (2015). Non-coding RNA: what is functional and what is junk?Front. Genet.6:2. 10.3389/fgene.2015.00002
41
PanagiotakiN.Dajas-BailadorF.AmayaE.PapalopuluN.DoreyK. (2010). Characterisation of a new regulator of BDNF signalling, Sprouty3, involved in axonal morphogenesis in vivo.Development1374005–4015. 10.1242/dev.053173
42
PizzicannellaJ.DiomedeF.GugliandoloA.ChiricostaL.BramantiP.MerciaroI.et al (2019a). 3D Printing PLA/Gingival Stem Cells/EVs Upregulate miR-2861 and -210 during Osteoangiogenesis commitment.Int. J. Mol. Sci.20:3256. 10.3390/ijms20133256
43
PizzicannellaJ.GugliandoloA.OrsiniT.FontanaA.VentrellaA.MazzonE.et al (2019b). Engineered extracellular vesicles from human periodontal-ligament stem cells increase VEGF/VEGFR2 expression during bone regeneration.Front. Physiol.10:512. 10.3389/fphys.2019.00512
44
RajanT. S.GiacoppoS.DiomedeF.BalleriniP.PaolantonioM.MarchisioM.et al (2016). The secretome of periodontal ligament stem cells from MS patients protects against EAE.Sci. Rep.6:38743. 10.1038/srep38743
45
RupaimooleR.SlackF. J. (2017). MicroRNA therapeutics: towards a new era for the management of cancer and other diseases.Nat. Rev. Drug Discov.16203–222. 10.1038/nrd.2016.246
46
SepulvedaJ. L.WuC. (2006). The parvins.Cell. Mol. Life Sci.6325–35. 10.1007/s00018-005-5355-1
47
SeverinoP.OliveiraL. S.AndreghettoF. M.TorresN.CurioniO.CuryP. M.et al (2015). Small RNAs in metastatic and non-metastatic oral squamous cell carcinoma.BMC Med. Genomics8:31. 10.1186/s12920-015-0102-4
48
SharmaE.Sterne-WeilerT.O’hanlonD.BlencoweB. J. (2016). Global mapping of human RNA-RNA interactions.Mol. Cell62618–626. 10.1016/j.molcel.2016.04.030
49
ShenW.-W.ZengZ.ZhuW.-X.FuG.-H. (2013). MiR-142-3p functions as a tumor suppressor by targeting CD133, ABCG2, and Lgr5 in colon cancer cells.J. Mol. Med.91989–1000. 10.1007/s00109-013-1037-x
50
SilvestroS.ChiricostaL.GugliandoloA.PizzicannellaJ.DiomedeF.BramantiP.et al (2020). Extracellular vesicles derived from human gingival mesenchymal stem cells: a transcriptomic analysis.Genes (Basel)11:118. 10.3390/genes11020118
51
SlabyO.LagaR.SedlacekO. (2017). Therapeutic targeting of non-coding RNAs in cancer.Biochem. J.4744219–4251. 10.1042/BCJ20170079
52
SongL.YangJ.DuanP.XuJ.LuoX.LuoF.et al (2013). MicroRNA-24 inhibits osteosarcoma cell proliferation both in vitro and in vivo by targeting LPAATβ.Arch. Biochem. Biophys.535128–135. 10.1016/j.abb.2013.04.001
53
StelzerG.RosenN.PlaschkesI.ZimmermanS.TwikM.FishilevichS.et al (2016). The GeneCards suite: from gene data mining to disease genome sequence analyses.Curr. Protoc. Bioinformatics541.30.1–1.30.33. 10.1002/cpbi.5
54
ThomasP.PranatharthiA.RossC.SrivastavaS. (2019). RhoC: a fascinating journey from a cytoskeletal organizer to a cancer stem cell therapeutic target.J. Exp. Clin. Cancer Res.38:328. 10.1186/s13046-019-1327-4
55
ThomasP. D.CampbellM. J.KejariwalA.MiH.KarlakB.DavermanR.et al (2003). PANTHER: a library of protein families and subfamilies indexed by function.Genome Res.132129–2141. 10.1101/gr.772403
56
TrapnellC.HendricksonD. G.SauvageauM.GoffL.RinnJ. L.PachterL. (2013). Differential analysis of gene regulation at transcript resolution with RNA-seq.Nat. Biotechnol.31:46. 10.1038/nbt.2450
57
TrubianiO.GiacoppoS.BalleriniP.DiomedeF.PiattelliA.BramantiP.et al (2016). Alternative source of stem cells derived from human periodontal ligament: a new treatment for experimental autoimmune encephalomyelitis.Stem Cell Res. Ther.7:1. 10.1186/s13287-015-0253-4
58
TrubianiO.MarconiG. D.PierdomenicoS. D.PiattelliA.DiomedeF.PizzicannellaJ. (2019). Human oral stem cells, biomaterials and extracellular vesicles: a promising tool in bone tissue repair.Int. J. Mol. Sci.20:4987. 10.3390/ijms20204987
59
ValadiH.EkströmK.BossiosA.SjöstrandM.LeeJ. J.LötvallJ. O. (2007). Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells.Nat. Cell Biol.9654–659. 10.1038/ncb1596
60
WangL.BoX.ZhengQ.XiaoX.WuL.LiB. (2016). miR−296 inhibits proliferation and induces apoptosis by targeting FGFR 1 in human hepatocellular carcinoma.FEBS Lett.5904252–4262. 10.1002/1873-3468.12442
61
WangY.ZhaoW.FuQ. (2013). miR-335 suppresses migration and invasion by targeting ROCK1 in osteosarcoma cells.Mol. Cell. Biochem.384105–111. 10.1007/s11010-013-1786-4
62
WeiJ.-J.WuX.PengY.ShiG.OlcaB.YangX.et al (2011). Regulation of HMGA1 expression by microRNA-296 affects prostate cancer growth and invasion.Clin. Cancer Res.171297–1305. 10.1158/1078-0432.ccr-10-0993
63
WuL.CaiC.WangX.LiuM.LiX.TangH. (2011). MicroRNA−142−3p, a new regulator of RAC1, suppresses the migration and invasion of hepatocellular carcinoma cells.FEBS Lett.5851322–1330. 10.1016/j.febslet.2011.03.067
64
XuC.LiS.ChenT.HuH.DingC.XuZ.et al (2016). miR-296-5p suppresses cell viability by directly targeting PLK1 in non-small cell lung cancer.Oncol. Rep.35497–503. 10.3892/or.2015.4392
65
XuG.WangJ.JiaY.ShenF.HanW.KangY. (2014). MiR-142-3p functions as a potential tumor suppressor in human osteosarcoma by targeting HMGA1.Cell. Physiol. Biochem.331329–1339. 10.1159/000358700
66
XuG. Q.LiL. H.WeiJ. N.XiaoL. F.WangX. T.PangW. B.et al (2019). Identification and profiling of microRNAs expressed in oral buccal mucosa squamous cell carcinoma of Chinese hamster.Sci. Rep.9:15616. 10.1038/s41598-019-52197-3
67
XuJ.-Z.ZhangJ.-L.ZhangW.-G. (2018). Antisense RNA: the new favorite in genetic research.J. Zhejiang Univ. Sci. B19739–749. 10.1631/jzus.B1700594
68
XuY.ZhaoF.WangZ.SongY.LuoY.ZhangX.et al (2012). MicroRNA-335 acts as a metastasis suppressor in gastric cancer by targeting Bcl-w and specificity protein 1.Oncogene311398–1407. 10.1038/onc.2011.340
69
ZhangP.YeoJ. C.LimC. T. (2019). Advances in technologies for purification and enrichment of extracellular vesicles.SLAS Technol.24477–488. 10.1177/2472630319846877
Summary
Keywords
human periodontal ligament stem cells, extracellular vesicles, next generation sequencing, non-coding RNAs, microRNAs, cytokinesis
Citation
Chiricosta L, Silvestro S, Gugliandolo A, Marconi GD, Pizzicannella J, Bramanti P, Trubiani O and Mazzon E (2020) Extracellular Vesicles of Human Periodontal Ligament Stem Cells Contain MicroRNAs Associated to Proto-Oncogenes: Implications in Cytokinesis. Front. Genet. 11:582. doi: 10.3389/fgene.2020.00582
Received
05 March 2020
Accepted
12 May 2020
Published
04 June 2020
Volume
11 - 2020
Edited by
Saverio Bellusci, University of Giessen, Germany
Reviewed by
Federica Papaccio, INCLIVA Health Research Institute, Spain; Roberto Di Primio, Marche Polytechnic University, Italy
Updates
Copyright
© 2020 Chiricosta, Silvestro, Gugliandolo, Marconi, Pizzicannella, Bramanti, Trubiani and Mazzon.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Emanuela Mazzon, emanuela.mazzon@irccsme.it
†These authors have contributed equally to this work
This article was submitted to Stem Cell Research, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.