1. Introduction
Oral squamous cell carcinoma (OSCC), with a 5-year survival rate of 50-60%, is the major cause of oral cancer morbidity and mortality and has become a major global health challenge [
1,
2]. In the past 20 years, despite therapeutic improvements, the prognosis of OSCC has not improved significantly, mainly due to late diagnosis, frequent locoregional recurrence and metastasis [
3]. Angiogenesis has been reported to be associated with tumor metastasis, disease progression, and poor prognosis in OSCC [
4]. It is the complex biological process of new blood vessels formation involving interactions between vascular cells and the extracellular environment [
5]. Tumor cells can induce angiogenesis by secreting various growth factors, cytokines, proteases, and even exosomes to the surrounding environment, and then provide oxygen and nutrients through new blood vessels to maintain tumor growth [
5]. However, the detailed mechanisms between tumorigenicity and angiogenesis in OSCC are not fully understood.
Heparan sulfate proteoglycans (HSPGs) are a group of distinct glycoproteins, found on the cell surface and in the extracellular matrix (ECM), composed of a core protein with one or more covalently linked heparan sulfate (HS), a type of glycosaminoglycan (GAG) [
6,
7]. Both HS GAGs and HSPGs bind various ligands and play important roles in normal physiology and pathological effects, including embryonic development, homeostasis, inflammatory responses, microbial infection, angiogenesis, and tumor initiation and progression [
7,
8]. However, HS chains appear to determine the binding affinity to various ligands such as growth factors, cytokines, chemokines and ECM components [
7,
9]. HS are unbranched polysaccharides composed of uronic acid (glucuronic acid, GlcA, or iduronic acid, IdoA) and N-acetyl-glucosamine (GlcNAc) disaccharide repeats with various sulfation modifications [
10]. Sulfation of HS is an important post-translational modification catalyzed by 2-O, 3-O and 6-O sulfotransferases. 2-O sulfation occurs on GlcA and IdoA residues with a strong preference for IdoA, while 3-O and 6-O sulfation occur on GlcNAc residues [
11,
12]. 2-O-sulfation within HS is catalyzed by HS 2-O sulfotransferase (HS2ST1) and is critical for HS to participate in various ligand-binding capabilities and signaling transductions [
13]. Abnormal expression and function of HS2ST1 have been detected in breast cancer, where upregulation of HS2ST1 is associated with reduced invasive behavior and cancer stem cells phenotype [
14,
15]. Furthermore, it has been shown that upregulation of HS2ST1 lead to conformational change in HS and reduces basic fibroblast growth factor binding, resulting in reduced p44/42 MAPK and p38 MAPK activity, which correlates with a reduced invasive phenotype. [
14]. These results highlight the tumor suppressor role of HS2ST1 in cancer. HSPGs are also present on the surface of vascular endothelial cells and contribute to the regulation of vascular permeability, immune cell trafficking, ligands binding and blood vessels formation [
16,
17]. Although HS2ST1 plays a key role in cancer progression, its regulation and mechanisms involved in vascular endothelial cells have not been clearly elucidated.
It is widely recognized that tumor-derived exosomes possess the capacity to regulate cell-cell communication between the tumor and surrounding tumor microenvironment (TME) by transferring cargoes such as proteins, nucleic acids, and lipids [
18]. One of the most important cargoes in exosomes are microRNAs (miRNAs), which belong to the family of single-stranded non-coding RNAs, around 18-25 nucleotides, involved in the gene regulation at the post-transcriptional level [
19]. Accumulating evidence suggests that tumor-derived exosomal miRNAs play an key role in stimulating endothelial cell migration and blood vessel formation that in turn support tumor progression and metastasis [
20]. For example, in OSCC derived exosomes, miR-221 has been reported to downregulate phosphoinositide-3-kinase regulatory subunit 1 in endothelial cells, thereby promoting cellular migration and tube formation [
21]. In addition, Wang et al. also reported that OSCC exosomal miR-210-3p targeting EFNA3 promotes angiogenesis in oral cancer via PI3K/AKT pathway [
22]. In this study, we found that miR-421 is abundantly present in OSCC-derived exosomes and was able to significantly induce tube formation in endothelial cell. MiR-421 has been shown to exert oncogenic and inhibitory effects on programmed cell death protein 4 in OSCC cells [
23]. However, there are currently no studies in the literature on the related role of OSCC-derived exosomal miR-421 in human endothelial cells.
Here, we focused on exosomal miR-421 to study its effect on OSCC angiogenesis. We demonstrated that miR-421 could be transferred to human umbilical vein vascular endothelial cells (HUVECs) via exosomes secreted by OSCC cells. Furthermore, the internalized exosome-miR-421 targets HS2ST1 through VEGF-mediated activation of ERK and AKT, thus promoting HUVEC migration, invasion, and tube formation. Maternally expressed gene 3 (MEG3), a long non-coding RNA (lncRNA), has been reported to function as a tumor suppressor gene in various solid cancers [
24,
25]. We further found that downregulation of lncRNA MEG3 can lead to upregulation of miR-421 in OSCC cells. In addition, treatment the exosomes derived from MEG3-overexpressing OSCC cells significantly increased the expression of HS2ST1 protein and, conversely, decreased the phosphorylation of VEGFR2, ERK, AKT and tube formation of HUVEC cells. Taken together, our results suggest that the lncRNA MEG3-mediated exosomal miR-421/HS2ST1 axis may be a new therapeutic target for OSCC angiogenesis.
3. Discussion
Angiogenesis is a crucial process in solid cancer tumorigenesis, promoting tumor growth, recurrence and metastasis [
5]. Despite large evidence suggested that targeting angiogenesis can provide clinical advantage in the treatment of solid tumors [
40]. However, there are no FDA-approved anti-angiogenic drugs for OSCC, although this cancer occurs in highly vascularized tissues [
41]. In this study, we showed that OSCC-derived exosomal miR-421 can target endothelial cells to increase angiogenesis by targeting HS2ST1 via VEGFR2 and downstream signaling pathways. Furthermore, we also found that miR-421 was negatively regulated by lncRNA MEG3, which acts as a ceRNAs in OSCC (
Figure 8). Therefore, our study not only suggests that blocking the transmission of pro-angiogenic miRNAs in exosomes might be an effective antiangiogenic therapeutic approach for OSCC, but also provides insights about the regulation and functionality of lncRNA MEG3.
LncRNA MEG3 is located in the imprinted DLK1–MEG3 locus on human chromosome 14q32 region [
42]. An increasing number of studies have shown that lncRNA MEG3 is hypermethylated and downregulated in various cancer tissues and cancer cell lines, playing the role of a tumor suppressor gene, and is a good biomarker for cancer diagnosis and prognosis [
43]. In terms of the role of ceRNA, lncRNA MEG3 can bind different miRNAs like a “sponge” and contribute to the regulatory network of miRNAs. For example, the neuroprotective effect of MEG3 silencing on PC12 cells is through the regulation of PI3K/AKT and NF-kB pathways by sponging miR-21 [
44]. Likewise, another study showed that down-expression of MEG3 could promote the proliferation of colonic stem cells in early colorectal cancer through binding miR-708 [
45]. Furthermore, previous studies have shown that MEG3 regulated EMT by sponging miR-421 targeting E-cadherin in HNSCC [
46]. Ye et al. found that anisomycin inhibited the downstream signaling of the Notch1 pathway by attenuating the molecular sponge effect of the lncRNA Meg3 with miR-421 in ovarian cancer cells [
47]. Consistent with these findings, we also showed that MEG3 exerted its regulatory role in oral cancer by interacting with miR-421, and includes several novel observations: (1) MEG3 restoration not only reduced the levels of miR-421 in OSCC cells but also resulted in decreased miR-421 levels in OSCC-deriver exosomes. (2) OSCC-deriver exosomes from MEG3 restoration resulted in reduced VEGFR2, ERK, and AKT activities, leading to reduced tube formation in HUVECs. Our study shows that down-regulated MEG3 can increase the expression of miR-421, which is then packaged in OSCC cell-derived exosomes and secreted into the extracellular microenvironment to target vascular endothelial cells and induce angiogenesis. To our knowledge, this is the first study to show that lncRNA MEG3 can affect exosomal miRNA secretion from OSCC and cause angiogenesis in the tumor microenvironment.
Recently, exosomes have been receiving enormous attention in both basic and translational medicine [
48]. This is mainly because exosomes derived from cancer cells can change the behavior of target cells in tumor microenvironment, and more importantly, exosomes have shown utility in the diagnosis and treatment of various cancers [
49]. Zhang et al. reported that treatment of esophageal squamous cell carcinoma cells with the exosome inhibitor SW4869 largely attenuated the migration and tube formation of HUVECs [
50]. Collectively, these results demonstrate that cancer cells-derived exosomes exerts an angiogenic effect and suggest that exosomes are indeed potential targets for cancer therapy, whether by blocking the generation of exosomes in cancer cells or by inhibiting the exosome internalization in HUVECs. Many in vitro studies and pre-clinical in vivo studies indicate that lots of compounds have the ability to block the formation and release of exosomes, which may prove useful for interfering the progression of certain types of cancers [
51]. For example, imipramine (acid sphingomyelinase inhibitor), simvastatin (HMG-CoA reductase inhibitor) and DPTIP (neutral sphingomyelinase 2 inhibitor) are designed for blocking exosomes biogenesis. Calpeptin (an inhibitor of the Ca
2+-dependent protease), and manumycin A (Ras farnesyltransferase inhibitor) are designed for inhibiting exosomes release. On the other hand, methyl-β-cyclodextrin inhibits caveolin-mediated endocytosis, which blocks exosomes uptake in recipient cells. These inhibitors affect different stages of exosome biogenesis and trafficking, and disruption of cell-to-cell communication from cancer cell-derived exosomes appears to provide an exciting new approach to cancer therapy [
52]. Currently, based on research results and clinical trials, exosomes are considered as immunotherapeutic vaccines for various cancers [
48]. Dendritic cells (DCs)-derived exosomes have been found to activate naturel killer cells and antigen-specific cytotoxic T lymphocytes in the treatment of patients [
53]. DCs-derived exosomes not only contain various bioactive molecules responsible for antigen presentation, but also have a smaller size to escape clearance by the mononuclear phagocyte system [
48]. Thus, using natural sources of exosomes from immune cells to suppress tumor growth represent an important strategy for cancer therapy.
Emerging evidence suggests that exosomal cargo, especially miRNA, is valuable for predicting disease progression and prognosis [
54]. In addition to avoiding RNase-mediated degradation, another advantage of packaging miRNAs into cellular exosomes is that they can be easily taken up by neighboring cells or carried to distant sites, thus potentially inducing phenotypic changes in recipient cells. Since oral cancer cells-derived exosomes can enhance endothelial cell migration and tube formation, the purpose of this study was to identify and characterize selective exosomal miRNA profiles and speculate on their potential targets involved in angiogenesis. To this end, we performed miRNA array analysis on oral cancer cell-derived exosomes and tissue samples from OSCC patients. Our study showed that high expression of miR-421 in oral cancer cell and OSCC patients is clinically associated with poor patient survival, suggesting that miR-421 may be a useful biomarker for OSCC prognosis. Additionally, we found that HS2ST1 expression was lower in HUVECs treated with exosomal miR-421. Restoration of HS2ST1 reversed exosomal miR-421-induced HUVECs tube formation, suggesting that HS2ST1 contributes to miR-421-mediated angiogenesis. Aberrant expression of HE2ST1 is frequently observed in tumors, suggesting that HS2ST1 plays an important role in tumor progression. For example, upregulation of HS2ST1 was associated with reduced invasive ability by attenuating FGF-2-induced MAPK activation in breast cancer cells [
14]. In addition, overexpression of HS2ST1 can decrease the stemness phenotype in triple negative MDA-MB-231 breast cancer cell line, possibly by regulating notch and Wnt-signaling pathways [
55]. However, little literature discusses the role of HS2ST1 in vascular endothelial cells. In this report, we demonstrate a novel role for HS2ST1 in regulating tube formation in vascular endothelial cells through VEGFR2-mediated MAPK signaling. Importantly, HS2ST1 overexpression reduced VEGF-induced phosphorylation of VEGFR2, ERK, and AKT, suggesting that HS 2-O-sulfation by HS2ST1 interferes the binding of VEGF to its receptor, VEGFR2. The mechanism by which HS2ST1 inhibits VEGF signaling is not well understood, however, previous studies have shown that 2-O-sulfated HS oligosaccharides are highly negatively-charge, which interacts with positively charged ligand proteins, thereby blocking ligand-receptor binding [
37]. Another possible speculation is that the sulfated HS by HS2ST1 may induce a conformational change in HS that hinders the binding of VEGF and VEGFR2. Therefore, whether it is the charge change or the configuration change that blocks the ligand engagement needs further investigation.
4. Materials and Methods
4.1. Cell Culture
OSCC cell lines (OEC-M1 and TW2.6) were cultured as previously described [
26]. Transformed normal human keratinocyte (OKF4/hTERT) was obtained from Rheinwald lab, and cultured in oral keratinocyte medium (OKM; ScienCell Research Laboratories, Carlsbad, CA, USA) according to the manufacturer’s instructions. HUVEC cells were cultured in M199 medium (Gibco, Grand Island, NY, USA) supplemented with 500 units/mL of heparin (Sigma-Aldrich, St. Louis, MO, USA) and 0.03 mg/mL of endothelial cell growth supplement (Sigma-Aldrich). All cells were definite by cell morphology and growth features, and cultured at 37 °C in 10% fetal bovine serum (FBS, Kibbutz BeitHaemek, Israel) within 2 months of resuscitation from the frozen stock.
4.2. Exosome Purification, Characterization and Quantification
Exosomes from the culture medium were purified by ultracentrifugation method. Briefly, cell culture medium was replaced with serum-free medium for 48 h. The collected culture medium was centrifuged at 300xg for 10 min to remove cellular debris. The supernatant was then centrifuged at 1,200xg for 20 min and 10,000xg for 30 min to remove large vesicles. Finally, the supernatant was ultracentrifuged at 110,000×g for 2 h to collect exosome pellet using Beckman Optima XL-90 ultracentrifuge (GMI Inc., Ramsey, MN, USA). The obtained exosome pellets were lastly washed and resuspended in PBS.
The morphology of exosome pellets was examined under a Hitachi H-7650 transmission electron microscope (TEM)( Hitachi Ltd., Tokyo, Japan). Briefly, 5 μL of exosomes in PBS drops were loaded onto carbon-coated copper grids and left at room temperature for 10 min. The exosomes on the copper mesh were washed with sterile distilled water, and then incubated with 1% phosphotungstic acid for 5 min. Let the samples dry under the lamp for 10 min before viewing on the TEM. For examining the particle size and quantity of exosomes, dynamic light scattering (DLS) was performed using a Nano Zetasizer (Malvern Instruments Inc., Westborough, MA, USA) according to the operating instructions.
The BCA kit was used to measure the protein concentration of exosomes. 25 μg of exosomes were used for each experiment. Exosome-specific surface proteins, CD9 (GTX100912, GeneTex, Irvine, CA, USA) and CD81 (GTX101766, GeneTex), were detected by western blot.
4.3. Tube Formation Assay
For tube formation assay, a precooled 96-well plate was coated with 10 μl of Matrigel (BD Biosciences) in each well and incubated at room temperature for 1 h. A total of 3,500 HUVEC cells of cell suspension was added to each well and incubated for 6 h, and the tubes that formed were imaged and counted with image J software package. Mean number of meshes was statistically analyzed and used as angiogenesis index.
4.4. Cell Migration and Invasion Assays
The 24-well Fluoro-Blok inserts (#351152, BD Biosciences, Franklin Lakes, NJ) were utilized to measure the invasion and migration ability of HUVECs. 3x104 cells were seeded in the upper chamber coated with (for invasion assay) or without (for migration assay) Matrigel (BD Biosciences) in culture medium containing 10% NuSerum. Following an 8-hour incubation for migration assays and 24-hour incubation for invasion assays, cells that had crossed the Fluoro-Blok membrane were fixed with 95% methanol and stained with propidium iodide. Fluorescent images were then counted using the Analytical Imaging Station software e (Imaging Research, Ontario, Canada).
4.5. RNA Extraction, Reverse Transcription and Quantitative-PCR (q-PCR)
TRIzol reagent (Life Technologies, Gaithersburg, MD, USA) was used to purify total RNA and exosomes based on the manufacturer’s protocol. RNA concentration was determined by NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA). For mRNA and lncRNA analysis, the cDNA was synthesized using random hexamer primers and SuperScript III reverse transcriptase (Invitrogen, Carlsbad, CA, USA). For miRNA analysis, the cDNA was synthesized using specific stem-loop RT primers and TaqMan MicroRNA Reverse Transcription Kit (Applied Biosystems, Carlsbad, CA, USA). q-PCR analysis was performed using Omics Green EvaGreen q-PCR Master Mix (OMICS Biotechnology, New Taipei City, Taiwan) to detect gene expression levels, and using QuantiTect SYBR Green PCR System (Qiagen, Hilden, Germany) to detect miRNA expression levels, respectively, according to the manufacturer’s instructions on the ABI StepOnePlus Real-time PCR system (Applied Biosystems). Bdi-miR-159-5p was used as the spike-in control for exosomal miRNA analysis, while, GAPDH and RUN44 were used as the internal control for the detection of mRNA and miRNA in cells. All primers used for this study are listed in
Supplementary Table S4.
4.6. Exosomal miRNA Profiling
Exosomal miRNA isolated from OKF4/hTERT, OEC-M1 and TW2.6 cells were extracted using TRIzol reagent (Life Technologies). The miRNA RNA purity, concentration and integrity were determined using a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA) and Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA), respectively. Exosomal miRNA profiling was performed using Affymetrix, GeneChipTM miRNA 3.0 array (Thermo Fisher Scientific) according to the manufacturer’s protocol.
4.7. In Vivo Matrigel Plug Assay
A total of 400 μl of growth factor-reduced Matrigel HC (#354262, Corning Inc., Corning, NY, USA), containing 15U of heparin (#H3393, Sigma-Aldrich) and 100 μl of miR-421 mimics (PM421) or scramble control (SC) was injected into the right flank of 6-week-old male nude mice (BALB/cAmM.CgFoxn1nu/CrlNarl), 4 mice per group. The plug was removed after 12 days, half of the plug was formalin-fixed and sliced for immunohistochemistry [
27], and the other half of the plug was used for hemoglobin analysis. The presence of the blood endothelial cell-specific marker was assessed with CD31 antibody (#3528, Cell signaling, Danvers, MA, USA). Hemoglobin concentration was measured using a hemoglobin kit (#EIAHGBC, Invitrogen, Carlsbad, CA) following the manufacturer’s protocol to indicate as an angiogenic index.
4.8. Clinical Samples and Microarray Profiling
A total of 40 paired OSCC tissue, including tumor and their adjacent non-tumorous tissues, were obtained from the National Cheng Kung University Hospital (Tainan, Taiwan) as previously described [
26]. The research protocol was reviewed and approved by the Research Ethics Committee of National Health Research Institutes (No: EC1040409-E) and Institutional Human Experiment and Ethic Committee of National Cheng Kung University Hospital (HR-97-100). These matched pairs of oral tumor/adjacent normal (T/N) tissues were prepared for microarray analysis using the whole genome DASL HumanRef-8-v3 chip and the Human-2v MicroRNA Expression BeadChips (Illumina, Inc.). These microarray profiling data are deposited in Gene Expression Omnibus (GEO) under accession number GSE37991 for gene expression and GSE45238 for miRNA expression.
4.9. Plasmids, Transfection and Virus Infections Assays
To construct the dual-luciferase reporter vectors, the entire 3’-untranslated region (3’-UTR) of HS2ST1 fragment and the entire MEG3 gene, containing target sequences of miR-421, were PCR amplified and cloned into the pmirGLO vector (Promega, Madison, WI, USA) according to the manufacturer’s instructions. The miR-421 binding site mutation vectors were also constructed by using Site-Directed Mutagenesis Kit (Stratagene, La Jolla, CA), and all the constructs were verified by DNA sequencing.
For transfection of the plasmids, cells were transfected with 2 μg of plasmids using Lipofectamine 2000 (Invitrogen) according to the manufacturer’s protocol. The miRNA mimics (PM) and scramble control (SC) were obtained from Thermo Fisher Scientific, the small interfering RNAs (siRNAs) of HS2ST1 (si-HS2ST1#1, #2) and si-non-targeting control (si-NT) were purchased from Horizon Discovery (Cambridge, UK). The nucleotide transfection was performed using Lipofectamine RNAiMAX (Thermo Fisher Scientific) according to the manufacturer’s instructions.
The lentiviral vectors expressing HS2ST1 or the control (NT) were transfected into the packaging cell line 293FT, along with pMD. G and pCMV△R8.91 plasmid, using the Polyjet transfection reagent (SignaGen Lab, Ijamsville, MA, USA), and the lentivirus was collected after 48 h to infect HUVEC cells.
4.10. Luciferase Reporter Assay
Cells were cotransfected with reporter vector containing the 3’-UTR of HS2ST1 or MEG3 sequence and miR-421 mimics (PM421, 20 nM) or scramble control (SC, 20 nM). After 48 h, firefly and Renilla luciferase activities were detected using the Dual Luciferase Reporter Assay System (Promega, Madison, WI, USA) on the Orion L luminometer (Berthold, GmbH, Pforzheim, Germany) according to the manufacturer’s protocol.
4.11. Protein Extraction and Western Blotting Analysis
Protein extraction and western blotting were performed as previously described [
28]. The primary antibodies used in the present study included HS2ST1 (sc-376530, Santa Cruz Biotechnology, Santa Cruz, CA, USA), CK18 (sc-6259, Santa Cruz Biotechnology), ERK (#9102, Cell Signaling), phosphor-ERK (p-ERK, #9101, Cell Signaling), AKT (#9272, Cell Signaling), phosphor-AKT (p-AKT, #9271, Cell Signaling), VEGFR2 (#24792, Cell Signaling), phosphor-VEGFR2 (p-VEGFR2, #2478, Cell Signaling). GAPDH (GTX100118, GeneTex) was used as an internal control. Signals from HRP-conjugated secondary antibodies were visualized using the enhanced chemiluminescence (ECL) detection system (PerkinElmer, Waltham, MA) and the chemiluminescence was exposed onto Kodak X-Omat film (Kodak, Chalon/Paris, France).
4.12. RNA-Immunoprecipitation (RNA-IP)
RNA-IP was performed and modified as previously described [
29]. Briefly, cells were infected with virus particles which contained pCDH-HA-2xMCP-SBP, pCDH-24xMS2 or pCDH-MEG3-24xMS2 plasmids for 48 h. Then, cells were lysed using RIPA extraction buffer containing RNase inhibitor and protease inhibitor cocktail (Sigma-Aldrich, Inc. Louis, MO, USA). The protein concentration was then determined using the BCA assay kit (Thermo Fisher Scientific) with bovine serum albumin as standard. Before immunoprecipitation, Streptavidin-Sepharose beads slurry (GE Healthcare, Chicago, IL) were pre-swollen in NT2 buffer (50 mM Tris-HCl, pH7.4, 150 mM NaCl, 1mM MgCl
2, 0.05% NP-40) supplemented with 20 mM DTT, 20 mM EDTA and 200U RNase inhibitor. Equal amounts of protein lysates were subjected to immunoprecipitation overnight mixed with Streptavidin-Sepharose beads at 4℃. Pellet beads were then washed with ice-cold NT2 buffer plus 0.5 M Urea (Mallinckrodt) for 4-5 times. To isolate RNA from immunoprecipitated pellets, phenol-chloroform-isoamyl alcohol with glycogen was added to the beads at -20℃ for overnight. Precipitate RNA was resuspend in a volume suitable for subsequent assay.
4.13. Fluorescence Analysis and Exosome Uptake Assay
For exosome uptake experiments, HUVECs were pretreated with endocytosis inhibitor (dynasore) or dimethyl sulfoxide (DMSO) solvent for 6 h. And then, HUVECs were washed with PBS and incubated with PKH26 (Sigma-Aldrich)-labeled exosomes (25 μg) isolated from OEC-M1 cells at 37 °C for 2 h. As a negative control, HUVECs were incubated with DMSO only. Uptake of exosomes in HUVECs was photographed under a TCS SP5 II Confocal Microscopy (Leica, Wetzlar, Germany). Exosome uptake was also assessed by measuring PKH26 fluorescence intensity using a FACSCalibur cytometer equipped with CellQuest™ Pro software (BD Biosciences).
4.14. Statistical Analysis
All statistical analyzes were performed using GraphPad Prism 5 software Ver.5.01 (San Diego, CA, USA), and results are expressed as mean ± standard deviation (SD). Differences between groups were analyzed by two-tailed Student’s t-test. Correlations between data pairs were performed by parametric Pearson Spearman correlation analysis. Values of p < 0.05 were considered statistically significant.
Author Contributions
S.-G.S., and C.-Y.H., conceived and designed the experiments. C.-Y.H., and S.-G.S., performed the data analysis and interpretation. C.-Y.H., S.-T.C., G.-H.W., Y.-M.H., and W.-J.C., performed the experiments. S.-T.C., performed the bioinformatics analysis. J.-R.H., contributed clinical materials. C.-Y.H., H.-D.W., and S.-G.S., were involved in the manuscript preparation. S.-G.S., contributed in the coordination of the study. All authors have read and agreed to the published version of the manuscript.
Figure 1.
Exosomal miR-421 derived from OSCC cells promotes HUVEC migration, invasion and tube formation. (A-C) Exosome isolated from normal oral keratinocyte (OKF4/hTERT) or OSCC cells (OEC-M1 and TW2.6) was incubated with HUVEC cells for migration assay (A), invasion assay (B), and tube formation assay (C). (D) Heat map of differential expressed exosomal miRNAs. Two biological replicates (for OKF4/hTERT) and three biological replicates (for OEC-M1 and TW2.6) were used for Affymetrix miRNA array. (E) Venn diagram of 12 differentially expressed miRNAs based on the intersection of OSCC-derived exosomal miRNAs with OSCC patient miRNA database (GSE37991). (F) Tube formation assay of HUVECs treated with 12 differentially expressed miRNA mimics. (G) Transwell migration and invasion assay of HUVEC cells treated with miR-421 mimics (PM421). (H) Relative expression level of miR-421 in exosomes isolated from miR-421 (PM421) or scramble control (SC)-transfected OEC-M1 cells using quantitative-RT-PCR (left). Exosome isolated from miR-421 (PM421) or scramble control (SC)-transfected OEC-M1 cells was incubated with HUVEC cells for migration, invasion, and tube formation assay (middle and right). Data are represented as mean ± SD; *p < 0.05; **p < 0.01; ***p < 0.001.
Figure 1.
Exosomal miR-421 derived from OSCC cells promotes HUVEC migration, invasion and tube formation. (A-C) Exosome isolated from normal oral keratinocyte (OKF4/hTERT) or OSCC cells (OEC-M1 and TW2.6) was incubated with HUVEC cells for migration assay (A), invasion assay (B), and tube formation assay (C). (D) Heat map of differential expressed exosomal miRNAs. Two biological replicates (for OKF4/hTERT) and three biological replicates (for OEC-M1 and TW2.6) were used for Affymetrix miRNA array. (E) Venn diagram of 12 differentially expressed miRNAs based on the intersection of OSCC-derived exosomal miRNAs with OSCC patient miRNA database (GSE37991). (F) Tube formation assay of HUVECs treated with 12 differentially expressed miRNA mimics. (G) Transwell migration and invasion assay of HUVEC cells treated with miR-421 mimics (PM421). (H) Relative expression level of miR-421 in exosomes isolated from miR-421 (PM421) or scramble control (SC)-transfected OEC-M1 cells using quantitative-RT-PCR (left). Exosome isolated from miR-421 (PM421) or scramble control (SC)-transfected OEC-M1 cells was incubated with HUVEC cells for migration, invasion, and tube formation assay (middle and right). Data are represented as mean ± SD; *p < 0.05; **p < 0.01; ***p < 0.001.

Figure 2.
miR-421 induces neovascularization in vivo. (A) Image of Matrigel plugs from mice for 12 days after injection of Matrigel mixed with scrambled control (SC) or miR-421 mimic (PM421)(upper). Hemoglobin concentration of plugs was determined and used as angiogenesis index (lower). (B) Representative images of IHC staining for H&E and CD31 from Matrigel plugs. Selective case number was showed from (A) (scale bar: 100 μm). Arrowheads show the location of microvascular structures with CD31 positive staining. (C) Clinical expression profile of miR-421 in Taiwanese OSCC cohort (GSE37991, left) and TCGA cohort (GDC Head and Neck Cancer, right).(D) Kaplan–Meier analysis of overall survival according to miR-421 expression from 239 OSCC patients in TCGA’s cohort (GDC Head and Neck Cancer). Data are represented as mean ± SD; *p < 0.05; ***p < 0.001.
Figure 2.
miR-421 induces neovascularization in vivo. (A) Image of Matrigel plugs from mice for 12 days after injection of Matrigel mixed with scrambled control (SC) or miR-421 mimic (PM421)(upper). Hemoglobin concentration of plugs was determined and used as angiogenesis index (lower). (B) Representative images of IHC staining for H&E and CD31 from Matrigel plugs. Selective case number was showed from (A) (scale bar: 100 μm). Arrowheads show the location of microvascular structures with CD31 positive staining. (C) Clinical expression profile of miR-421 in Taiwanese OSCC cohort (GSE37991, left) and TCGA cohort (GDC Head and Neck Cancer, right).(D) Kaplan–Meier analysis of overall survival according to miR-421 expression from 239 OSCC patients in TCGA’s cohort (GDC Head and Neck Cancer). Data are represented as mean ± SD; *p < 0.05; ***p < 0.001.
Figure 3.
HS2ST1 is a target of miR-421 in HUVECs. (A) A heatmap of 371 differentially expressed genes was obtained in HUVECs transfected with miR-421 (PM421, 20 nM) compared to control cells (SC) with three biological replicates using Affymetrix microarrays (ClariomTM S Assay). (B) Venn diagram of 27 miR-421-targeted genes based on the intersection of Affymetrix microarray with miRWalk database. (C) Heatmap of 27 miR-421-targeted genes. (D) Validation of the expression level of 27 genes in miR-421-transfected HUVEC cells using quantitative real time-PCR. (E) Schematic representation of the putative miR-421 binding sequence in the 3’-UTR of HS2ST1 with wild-type form (HS2ST1 3’-UTR WT) and deleted form (HS2ST1 3’-UTR Mut). The matched nucleotides are labeled with green and red color. (F) The effect of miR-421 (PM421, 20 nM) on the luciferase activities of the constructs containing the wild-type or mutant-type 3’-UTR in HUVEC cells. (G) Western blot analysis of HS2ST1 in HUVEC cells following miR-421 (PM4421, 20 nM) or scramble control (SC) transfection for 48 h. GAPDH was used as an internal control. (H) Western blot analysis of HS2ST1 protein in HUVEC cells following transfection of si-HS2ST1 (si-HS2ST1 #1 and #2) or si-non-targeting control (si-NT) for 48 h. (I-J) Tube formation assay of HUVECs treated with the indicated transfection. Data are represented as mean ± SD; **p < 0.01; ***p < 0.001.
Figure 3.
HS2ST1 is a target of miR-421 in HUVECs. (A) A heatmap of 371 differentially expressed genes was obtained in HUVECs transfected with miR-421 (PM421, 20 nM) compared to control cells (SC) with three biological replicates using Affymetrix microarrays (ClariomTM S Assay). (B) Venn diagram of 27 miR-421-targeted genes based on the intersection of Affymetrix microarray with miRWalk database. (C) Heatmap of 27 miR-421-targeted genes. (D) Validation of the expression level of 27 genes in miR-421-transfected HUVEC cells using quantitative real time-PCR. (E) Schematic representation of the putative miR-421 binding sequence in the 3’-UTR of HS2ST1 with wild-type form (HS2ST1 3’-UTR WT) and deleted form (HS2ST1 3’-UTR Mut). The matched nucleotides are labeled with green and red color. (F) The effect of miR-421 (PM421, 20 nM) on the luciferase activities of the constructs containing the wild-type or mutant-type 3’-UTR in HUVEC cells. (G) Western blot analysis of HS2ST1 in HUVEC cells following miR-421 (PM4421, 20 nM) or scramble control (SC) transfection for 48 h. GAPDH was used as an internal control. (H) Western blot analysis of HS2ST1 protein in HUVEC cells following transfection of si-HS2ST1 (si-HS2ST1 #1 and #2) or si-non-targeting control (si-NT) for 48 h. (I-J) Tube formation assay of HUVECs treated with the indicated transfection. Data are represented as mean ± SD; **p < 0.01; ***p < 0.001.

Figure 4.
HS2ST1 blocks VEGF-induced tube formation. (A) Western blot analysis of ERK, phosphor-ERK (p-ERK), AKT, phosphor-AKT (p-AKT) after transfecting the si-HS2ST1 (si-HS2ST1 #1 and #2) or si-non-targeting control (si-NT) for 48 h. GAPDH was used as an internal control. (B) Transwell migration and invasion assay of HUVEC cells following transfection of siRNA for 24h. (C) Tube formation assay of HUVEC cells after transfection with HS2ST1 expression vector or control vector (Vec). Western blot analysis showed the HS2ST1 protein level (right upper). Mean number of meshes was statistically analyzed (right bottom). (D) Transwell migration and invasion assay of HUVEC cells following transfection of HS2ST1 expression vector or control vector (Vec) for 24h. (E-H) Western blot analysis of ERK, phosphor-ERK (p-ERK), AKT, phosphor-AKT (p-AKT), HS2ST1, VEGFR2, phosphor-VEGFR2 (p-VEGFR2), N-Cadherin, snall, slug, and VE-Cadherin in HUVEC cells treated with the indicated treatment and transfection. GAPDH was used as an internal control. Data are represented as mean ± SD; *p < 0.05; **p < 0.01.
Figure 4.
HS2ST1 blocks VEGF-induced tube formation. (A) Western blot analysis of ERK, phosphor-ERK (p-ERK), AKT, phosphor-AKT (p-AKT) after transfecting the si-HS2ST1 (si-HS2ST1 #1 and #2) or si-non-targeting control (si-NT) for 48 h. GAPDH was used as an internal control. (B) Transwell migration and invasion assay of HUVEC cells following transfection of siRNA for 24h. (C) Tube formation assay of HUVEC cells after transfection with HS2ST1 expression vector or control vector (Vec). Western blot analysis showed the HS2ST1 protein level (right upper). Mean number of meshes was statistically analyzed (right bottom). (D) Transwell migration and invasion assay of HUVEC cells following transfection of HS2ST1 expression vector or control vector (Vec) for 24h. (E-H) Western blot analysis of ERK, phosphor-ERK (p-ERK), AKT, phosphor-AKT (p-AKT), HS2ST1, VEGFR2, phosphor-VEGFR2 (p-VEGFR2), N-Cadherin, snall, slug, and VE-Cadherin in HUVEC cells treated with the indicated treatment and transfection. GAPDH was used as an internal control. Data are represented as mean ± SD; *p < 0.05; **p < 0.01.

Figure 5.
Exosomal miR-421 promotes tube formation by targeting HS2ST1in HUVEC cells. (A and D) Images of tube formation assay of HUVECs treated with the indicated transfection. (B and E) Transwell migration, invasion and tube formation assay of HUVECs treated with the indicated transfection. Data are represented as mean ± SD; *p < 0.05; **p < 0.01; ***p < 0.001. (C and F) Western blot analysis of HS2ST1, VEGFR2, phosphor-VEGFR2 (p-VEGFR2), ERK, phosphor-ERK (p-ERK), AKT and phosphor-AKT (p-AKT) in HUVEC cells treated with the indicated transfection. GAPDH was used as an internal control.
Figure 5.
Exosomal miR-421 promotes tube formation by targeting HS2ST1in HUVEC cells. (A and D) Images of tube formation assay of HUVECs treated with the indicated transfection. (B and E) Transwell migration, invasion and tube formation assay of HUVECs treated with the indicated transfection. Data are represented as mean ± SD; *p < 0.05; **p < 0.01; ***p < 0.001. (C and F) Western blot analysis of HS2ST1, VEGFR2, phosphor-VEGFR2 (p-VEGFR2), ERK, phosphor-ERK (p-ERK), AKT and phosphor-AKT (p-AKT) in HUVEC cells treated with the indicated transfection. GAPDH was used as an internal control.
Figure 6.
LncRNA MEG3 inhibits tube formation via targeting miR-421. (A) Following MS2 RNA-immunoprecipitation in OEC-M1 cells, RNA was extracted by Trizol-chloroform and precipitated by isopropanol plus glycogen at -20 oC. The level of miR-421 was measure by qRT-PCR in pCDH-MEG3-MS2 group (MEG3) compared with pCDH-MS2 control vector group (Vec). (B) Upper, schematic representation of the putative miR-421 binding sequence in the LncRNA MEG3 with wild-type form (MEG3 WT) and mutant form (MEG3 Mut). The matched nucleotides are labeled with green and red color. The mutant nucleotides are labeled with blue color. Lower, the effect of miR-421 (PM421, 20 nM) on the luciferase activities of the constructs containing the wild-type or mutant-type MEG3 in OEC-M1 cells. (C & D) Relative expression of MEG3 and miR-421 in OEC-M1 and TW2.6 cells after transfection of MEG3 expression vector or control vector (Vec). (E-G) Relative miR-421 levels, western blot analysis and tube formation assay in HUVECs after incubation with exosomes derived from OEC-M1 and TW2.6 cells which were transfected with MEG3 expression vector or control vector (Vec). (H) MEG3 or miR-421 expression level by qRT-PCR in 40 of OSCC tumors (T) compared with their own adjacent normal tissues (N). Expression levels are expressed as the log2 ratios. (I) Correlation analysis of miR-421 and MEG3 in human OSCC patients. Data are represented as mean ± SD; **p < 0.01; ***p < 0.001.
Figure 6.
LncRNA MEG3 inhibits tube formation via targeting miR-421. (A) Following MS2 RNA-immunoprecipitation in OEC-M1 cells, RNA was extracted by Trizol-chloroform and precipitated by isopropanol plus glycogen at -20 oC. The level of miR-421 was measure by qRT-PCR in pCDH-MEG3-MS2 group (MEG3) compared with pCDH-MS2 control vector group (Vec). (B) Upper, schematic representation of the putative miR-421 binding sequence in the LncRNA MEG3 with wild-type form (MEG3 WT) and mutant form (MEG3 Mut). The matched nucleotides are labeled with green and red color. The mutant nucleotides are labeled with blue color. Lower, the effect of miR-421 (PM421, 20 nM) on the luciferase activities of the constructs containing the wild-type or mutant-type MEG3 in OEC-M1 cells. (C & D) Relative expression of MEG3 and miR-421 in OEC-M1 and TW2.6 cells after transfection of MEG3 expression vector or control vector (Vec). (E-G) Relative miR-421 levels, western blot analysis and tube formation assay in HUVECs after incubation with exosomes derived from OEC-M1 and TW2.6 cells which were transfected with MEG3 expression vector or control vector (Vec). (H) MEG3 or miR-421 expression level by qRT-PCR in 40 of OSCC tumors (T) compared with their own adjacent normal tissues (N). Expression levels are expressed as the log2 ratios. (I) Correlation analysis of miR-421 and MEG3 in human OSCC patients. Data are represented as mean ± SD; **p < 0.01; ***p < 0.001.

Figure 7.
Effects of exosomes uptake in HUVECs. (A) Images of exosome uptake at 2 h after dynasore treatment in HUVECs under a confocal microscope. Exosomes isolated from OEC-M1 were stained with PKH26 and added into the culture medium of dynasore-treated (20 or 40 μM) or untreated (Mock) HUVECs. As a negative control (NC), HUVECs were incubated with DMSO only. Scale bar: 50µm. (B) Exosome uptake was assessed by measuring PKH26 fluorescence intensity using a FACSCalibur cytometer. (C) Effect of exosome uptake on miR-421 level after dynasore treatment in HUVECs using qRT-PCR. (D) Effect of exosome uptake on protein level after dynasore treatment in HUVECs using western blot analysis. GAPDH was used as an internal control. (E) Effect of exosome uptake on tube formation at 6 h after dynasore treatment in HUVECs. Data are represented as mean ± SD; *p < 0.05; **p < 0.01; ***p < 0.001.
Figure 7.
Effects of exosomes uptake in HUVECs. (A) Images of exosome uptake at 2 h after dynasore treatment in HUVECs under a confocal microscope. Exosomes isolated from OEC-M1 were stained with PKH26 and added into the culture medium of dynasore-treated (20 or 40 μM) or untreated (Mock) HUVECs. As a negative control (NC), HUVECs were incubated with DMSO only. Scale bar: 50µm. (B) Exosome uptake was assessed by measuring PKH26 fluorescence intensity using a FACSCalibur cytometer. (C) Effect of exosome uptake on miR-421 level after dynasore treatment in HUVECs using qRT-PCR. (D) Effect of exosome uptake on protein level after dynasore treatment in HUVECs using western blot analysis. GAPDH was used as an internal control. (E) Effect of exosome uptake on tube formation at 6 h after dynasore treatment in HUVECs. Data are represented as mean ± SD; *p < 0.05; **p < 0.01; ***p < 0.001.
Figure 8.
Proposed model illustrating the exosomal-miR-421 modulates angiogenesis in OSCC. Downregulated lncRNA MEG3 increases the amount of miR-421, which is then packaged in OSCC-derived exosomes and secreted into the tumor microenvironment. MiR-421-enriched exosomes are taken up by HUVEC cells to induce angiogenesis through targeting HS2ST1. Downregulation of HS2ST1 leads to decreased 2-O-sulfation of HSPG, thereby facilitating VEGF binding to VEGFR2, which then promotes EMT, migration, invasion and angiogenic potential of HUVEC cells.
Figure 8.
Proposed model illustrating the exosomal-miR-421 modulates angiogenesis in OSCC. Downregulated lncRNA MEG3 increases the amount of miR-421, which is then packaged in OSCC-derived exosomes and secreted into the tumor microenvironment. MiR-421-enriched exosomes are taken up by HUVEC cells to induce angiogenesis through targeting HS2ST1. Downregulation of HS2ST1 leads to decreased 2-O-sulfation of HSPG, thereby facilitating VEGF binding to VEGFR2, which then promotes EMT, migration, invasion and angiogenic potential of HUVEC cells.