- Primary research
- Open Access
Identification of potential plasma biomarkers in early-stage nasopharyngeal carcinoma-derived exosomes based on RNA sequencing
Cancer Cell International volume 21, Article number: 185 (2021)
Early diagnosis of nasopharyngeal carcinoma (NPC) is vital to improve the prognosis of these patients. However, early diagnosis of NPC is typically challenging. Therefore, we explored the pathogenetic roles and associated mechanisms of exosomes in plasma of patients with early-stage NPC.
Exosomes in plasma were extracted by ultra-high-speed centrifugation. Western blot and transmission electron microscopy (TEM) were used to verify the purity of exosomes. The sequencing data (6 plasma samples from healthy volunteers vs. 6 NPC plasma samples) were analyzed by principal component analysis (PCA), DESeq2, gene ontology (GO), Kyoto encyclopedia of genes and genomes (KEGG), and TargetScan. The differentially expressed miRNAs (DEmiRNAs) were obtained from the dataset (GSE118720) downloaded from the Gene Expression Omnibus (GEO) repository. Additionally, the datasets downloaded from the GEO database (GSE12452, GSE13597, GSE53819, GSE64634) were used to predict the target genes and functions of hsa-miR-1301-3p. qPCR was applied to verify the differences in the expressions of hsa-miR-1301-3p between 10 normal plasma and 10 NPC plasma samples.
Western blot, TEM, and Nanoparticle Tracking Analysis showed adequate purity of the extracted exosomes. RNA-seq analysis revealed 21 upregulated miRNAs, and 10 downregulated miRNAs in plasma exosomes of early-stage NPC patients. GO analysis showed that the target genes of DEmiRNAs were mainly enriched in DNA synthesis and transcription regulation. KEGG analysis revealed that DEmiRNAs were mainly enriched in PI3K-Akt and MAPK signaling pathways. Moreover, the expression of hsa-mir-1301-3p was verified to be significantly upregulated in enlarged samples of plasma exosomes.
We identified several DEmiRNAs extracted from tumor-derived exosomes between normal plasma and early-stage NPC plasma. Bioinformatics analyses indicated that these DEmiRNAs may be related to NPC development. Our study may provide novel insights into underlying biomarkers and mechanisms of plasma exosomes in early-stage NPC.
Nasopharyngeal carcinoma (NPC) is a malignant cancer that shows wide variability in geographical distribution. The condition is highly prevalent in Southeast Asia, especially in the Guangdong Province of China . In 2018, an estimated 129,079 new cases of NPC were diagnosed across the world, accounting for 0.7% of all carcinomas . The reported 5-year overall survival (OS) rate of patients with early-stage NPC is 92.8–98.1% while the 5-year OS rate of patients with advanced NPC is 68.9–78.6% [3, 4]. Unfortunately, more than 70% of NPC patients are diagnosed at an advanced stage [5, 6]. Therefore, early diagnosis of NPC is vital to improve the prognosis of these patients. Currently, nasopharyngeal biopsy and imaging examination are the main modalities for the diagnosis of NPC. However, patients with early-stage NPC are typically asymptomatic and are unlikely to be subjected to these investigations; moreover, biopsy is an invasive procedure. Consequently, use of routine screening methods may not be very effective in facilitating early diagnosis of NPC in clinical practice . Therefore, unraveling the mechanisms of the occurrence and development of early-stage NPC and identification of potential molecular biomarkers for early diagnosis of NPC are key research imperatives.
Several risk factors have been implicated in the pathogenesis and progression of NPC, such as Epstein-Barr virus (EBV) infection, genetic susceptibility, and environmental factors [8, 9]. Several studies have indicated that EBV DNA in plasma could be regarded as a biomarker for the early diagnosis of NPC [10, 11]. However, there is no robust evidence supporting its role as a biomarker. Moreover, use of different polymerase chain reaction methods may yield inconsistent results. Different standards for detection tend to limit the clinical application of plasma EBV-DNA . Exosomes contain proteins, nucleic acids (DNA, mRNA, and miRNA), lipids and metabolites . Exosomes have been shown to regulate tumor-induced immune response and promote tumor progression as well as metastasis [14, 15]. Exosomes extracted from breast cancer and prostate cancer cells were shown to induce tumors by transferring miRNAs [16, 17]. miR-200 extracted from exosomes in metastatic breast carcinoma cells promoted epithelial-mesenchymal transformation (EMT) and metastasis of breast cancer cells .
In the contemporary literature, an increasing number of studies have investigated the role of exosomes in NPC diagnosis. David et al. found that exosomes extracted from NPC may regulate the tumor micro-environment and affect the growth of adjacent cells through the intercellular transfer of EBV latent membrane protein 1 (LMP1), signal molecules, and viral miRNAs . miR-34c in exosomes was found to inhibit the malignant behavior and radiation resistance of NPC by blocking the EMT process . However, just a few studies have focused on the early-stage NPC, and no reliable biomarkers for the early diagnosis of NPC have been identified. To the best of our knowledge, this is the first study to explore the underlying plasma biomarkers in early-stage NPC exosomes.
The purpose of our study was to explore the underlying mechanism of early-stage NPC and to identify potential markers for early diagnosis of NPC. In the present study, we identified the differentially expressed genes (DEGs) in the exosomes of early-stage NPC. Using the Gene Expression Omnibus (GEO) and the Cancer Genome Atlas (TCGA) databases, bioinformatics analyses were performed to identify the differentially expressed miRNAs (DEmiRNAs) and uncover the potential mechanisms of miRNAs in NPC. Finally, the DEmiRNAs were verified by plasma samples collected at our hospital. Our study may provide novel insights into the early diagnosis of NPC.
Materials and methods
Blood samples were collected from 16 patients with early-stage NPC who were treated at our hospital between November 2018 and September 2019. Plasma samples from 16 age- and sex-matched healthy donors (8 men and 8 women; median age: 50 years) were also collected. In line with the eighth edition of the International Union Against Cancer/American Joint Committee on Cancer guidelines for NPC , 7 patients had stage I NPC and 9 patients had stage II NPC. The ethics committee of our hospital approved the use of human tissue specimens in this work (project ethics number: SQ2019-018-01). Six pairs of blood samples were used for the extraction of plasma exosomes, and the other ten pairs of blood samples were used for further sample verification.
Isolation and identification of exosomes
Melt the samples at 37 °C at medium speed, and move the samples to a new centrifugal tube with centrifugation at 2000g at 4 °C for 30 min. Afterwards, carefully move the supernatant to a new centrifugal tube with centrifugation again at 10,000g at 4 °C for 45 min to remove the large vesicles. The supernatant was filtrated by a 0.45 μm filter membrane and the filtrate was collected. The filtrate was moved to a new centrifugal tube, and an overspeed rotor was selected for centrifugation at 4 °C and 100,000g for 70 min. The supernatant was removed and resuspended with 1 × PBS precooled by 10 mL. The overspeed rotor was selected and centrifuged again at 4 °C for 100,000g for 70 min. The supernatant was removed, and 100 μL precooled 1 × PBS was used for resuspend.
Dilute 20 μL exosomes to 90 μL, and add 20 μL fluorescently labeled antibodies (CD9, CD81, IgG) to 30 μL diluted exosomes, respectively. Then mix well, and incubate for 30 min at 37 °C without light. 1 mL precooled PBS was then added, an overspeed rotor was selected, and the rotor was centrifuged at 4 °C for 110,000g for 70 min. The procedure was repeated twice. The supernatant was carefully removed and resuspended with 1 × PBS precooled by 50 μL. Finally, the protein index results were detected by nanoflow cytometry measurement(NFCM).
Exosomes were mixed with an equal volume of 4% paraformaldehyde. Next, 10 µL specimens were moved onto Formvar/carbon-coated nickel grids. Subsequently, the membranes were adsorbed for 20 min at room temperature (RT). After washing and fixating using 50 µL 1% w/v glutaraldehyde, specimens were contrasted with 4% w/v uranyl acetate (UA) for 5 min and embedded on ice in a mixture of 4% UA and 2% methylcellulose (100:900) solution for 10 min. The excess liquid was sucked dry and the grill was dried by air. The shape and size of exosomes were checked by transmission electron microscopy (TEM) (JEM-2100 JEOL, Tokyo, Japan).
RNA extraction and qualification
Total RNA was isolated from the exosomes of early-stage NPC plasma and normal plasma using TRlzol reagent (Invitrogen, South San Francisco, CA, USA) and treated with RNase-free DNase I to remove genomic DNA. RNA degradation and contamination were assessed on 1% agarose gels. The purity, concentration, and integrity of the RNAs were checked using the NanoPhotometer® spectrophotometer (Implen, CA, USA), Qubit®. RNA Assay Kit in Qubit® 2.0 Fluorometer (Life Technologies, CA, USA),degradation and RNA Nano 6000 Assay Kit with the Agilent Bioanalyzer 2100 system (Agilent Technologies, CA, USA). All procedures were conducted according to the manufacturer’s instructions.
Library Preparation and Sequencing
The entire 3 mg total RNA of each sample was divided into aliquots and a small RNA (sRNA) library prepared. The NEBNext Multiplex sRNA Library Prep Set for Illumina (New England Biolabs, Ipswich, MA, USA) was applied to generate the sequencing library. The index code was added to assign a sequence to every specimen. TruSeq SR Cluster Kit v3-cBot-HS (Illumina, San Diego, CA, USA) was applied to cluster index-coded specimens. The libraries were sequenced on an Illumina Hiseq 2000 machine to generate 50 bp single-end reads. The entire process of sequencing was carried out at Shanghai Biotechnology Co. Ltd.
First, the adaptors and low-quality reads (Q30 < 80%) and low complexity reads were removed from the raw reads (unfiltered total reads obtained from the Hiseq 2500) using Fastp (version 0.20.1) and in-house perl scripts; the remaining reads were regarded as clean reads. Contaminants, the reads < 5 nt and trimming the adapters were filtered from these clean reads to obtain the valid reads. The effective reads were used to calculate length distribution and base preference. Moreover, the effective reads were aligned against a reference human genome (GRCh38) using Bowtie. Subsequently, the effective reads were aligned against Rfam (version 14.2) using Bowtie tools to annotate the rRNA, snRNA, snoRNA, sRNA, cis-reg, miRNA, and tRNA species. In the end, the unannotated reads were applied to predict new miRNAs using the miRDeep* software (Version v38). A new GTF file was constructed using the information of the novel miRNAs and the gff3 file in the miRBase. Subsequently, HTSeq v0.6.1 was applied to measure the reads mapped to every miRNA.
The expressions of miRNAs were measured as transcripts per million reads (TPM) by the criteria: Normalized expression = miRNA read count × 1,000,000/lib size (libsize is the sum of miRNA read count). The DESeq2 R package (1.26.0) was used to identify the DEmiRNAs between NPC and normal group. P-values were adjusted for multiple testing by DESeq2 default adjustments. The criteria of DEmiRNAs were P-values less than 0.05 and fold change more than 2.
Target prediction, GO, and KEGG
The potential target genes of DEmiRNAs were obtained from TargetScan (http://www.targetscan.org). TargetScan predicts targets of miRNAs by looking for sites that match the seed region of every miRNA. GO and KEGG analyses were applied to investigate the functions of target genes for the DEmiRNAs by the Database for Annotation, Visualization and Integrated Discovery (DAVID, version 6.8) (http://david.abcc.ncifcrf.gov/) and KOBAS software, respectively. GO and KEGG term with a P-value less than 0.05 were considered as significantly enriched by target genes.
The mRNA-miRNA network construction
Based on the predicted binding relationship between miRNAs and mRNAs, the mRNA-miRNA networks were constructed using Cytoscape v3.7.0 software (http://www.cytoscape.org/).
Quantitative polymerase chain reaction assay
The hsa-miR-1301-3p expression was quantified by qPCR. In brief, total RNAs were extracted from plasma exosomes obtained from early-stage NPC patients and normal donors using TRIzol reagent (Invitrogen; Carlsbad, CA, USA) and transcribed into cDNA using the RevertAid First Strand cDNA Synthesis Kit (Thermo; Waltham, MA, USA). qPCR assay was conducted by the SYBR Green qPCR Master (Roche). The cel-miR-39 was regarded as external control for quantification of the expressions of miRNAs. All determinations were performed in triplicate. Statistical analysis was performed by the SPSS 22.0 statistical software (SPSS, Inc.; Chicago, IL, USA). The Graph Pad Prism 8 drawing software (GraphPad; La Jolla, CA, USA) was used for generating graphs. Between-group differences were assessed using the Student’s t-test. P-values less than 0.05 were considered indicative of statistical significance.
GEO data processing
The Limma package  was used to conduct the normalization and logarithmic transformation with base 2; the DEGs between NPC and normal specimens were identified by Limma. The genetic integration of DEGs screened from the six sets of data was performed using an R package “RobustRankAggreg”  according to a robust rank aggregation (RRA) method,which is based on the null hypothesis of irrelevant inputs and identifies genes that consistently rank better than expected.The criterion was as follows: P-value < 0.01.
Isolation and characterization of exosomes
Exosomes are microvesicles (length: 30–120 nm) that can transport miRNAs, mRNAs and proteins, which are essential for the cell to cell communication. Exosomes in plasma were extracted by ultra-high-speed centrifugation. Several methods were applied to confirm that the substances isolated were exosomes as recommended in the minimal information for studies of extracellular vesicles 2018 (MISEV2018).
Western blot was used to determine the expressions of two exosome markers CD81 and CD63. We found abundant expressions of both CD81 and CD63 (Fig. 1a). These results indicated the high efficiency of ultra-high-speed centrifugation in extracting exosomes. Subsequently, TEM analysis revealed round particles with a diameter of approximately 100 nm, which corresponded to the unique characteristics of exosomes (Fig. 1b). Besides, the size of exosomes was checked using NFCM. We noticed that the plasma exosomes were round vesicles ranging from 50 to 150 nm (Fig. 1c, d). The concentration of these particles was further verified by NFCM (Additional file 1: Fig. S1). The results showed that the plasma exosomes were pure enough for further exploration.
Deep sequencing data
To verify the miRNA expression profile during the progression of NPC, 12 small RNA libraries (6 control and 6 NPC) were constructed for deep sequencing. After removal of poor-quality raw reads, total clean reads were obtained for subsequent bioinformatics analyses. Additional file 2: Table S1 in the supplemental material shows that the sequence reads were produced from all samples, with the single base sequencing error rate below 1%. Both Q20 and Q30 of each sample were > 90%. No GC bias was found. These clean reads were mapped to Rfam databases (Additional file 3: Table S2). The majority of reads were mapped to miRNAs and rRNAs.
We performed principal component analysis (PCA) of the miRNA expression profiles across 12 specimens to verify the clustering or outliers in the specimen set. Tight clustering was observed between normal samples, but not the tumor ones. Moreover, there was a partial overlap between NPC and normal samples, with NPC samples showing a significantly higher PC1 value than normal. Overall, there was considerable variability with respect to the relative expression patterns across the NPC samples obtained from patients with early-stage NPC (Fig. 2a).
Profiling the small RNA between control and sample
The types and contents of small RNAs in each sample of NPC and normal were recorded, respectively. The percentage content of small RNAs was calculated by groups (Fig. 2f, g). There was a significant difference in the content of small RNAs of the same kind among the groups, and the contents of miRNAs and rRNAs were the most in both groups.
The abundance of most miRNAs was relatively low in the samples. HTSeq was used to assess the expression of miRNAs in our study. The expression of identified miRNAs was standardized using TPM (transcripts per million; miRNAs were normalized by TPM). The detectable miRNAs identified in at least two samples were retained for subsequent analysis. Consequently, 1669 and 1395 miRNAs were screened; of these, 1356 were common between NPC and normal groups, and 313 and 39 were unique to NPC and normal groups (Fig. 2b). On further analysis of length distribution, the lengths of reads in all samples mostly ranged from 5 to 35 bp. The distribution of the length of small RNAs is depicted in the cluster shaped bar chart (Fig. 2d, e). The annotated sRNA size distribution peaked at 21 nt and 22 nt; 21-22 nt long sequences are regarded as the length of miRNAs in animals. Similar to the data shown in Fig. 2F, G, most small RNAs are miRNAs and rRNAs.
It is well known that the miRNA content is generally small; therefore, high expression of miRNA typically has a special function, which can be used as the focus of further research. We performed statistical analysis of the miRNA expression and drew a histogram. Read counts of screened miRNAs ranged from 0 to 355,139, the 20 most abundant miRNAs accounted for approximately 70% of all mapped reads. The remaining miRNAs accounted for less than 30% (Fig. 2c, Additional file 4: Table S3). Since the total amount of miRNA detected in the two groups is quite different, the ranking is based on the sum of the percentage of each miRNA in the two groups.
Next, we investigated the miRNAs associated with NPC. RNA-seq analysis showed the expression levels of 1708 miRNAs in all samples. Differential expression analysis using DEseq revealed 21 (15 known and 6 novel) down-regulated and 10 (9 known and 1 novel) up-regulated miRNAs (fold change ≥ 2, adjusted P-value < 0.05) (Fig. 3a). A heatmap of DEmiRNAs between NPC and normal groups is shown in Fig. 3b.
GO, and KEGG of predicted target genes
To investigate the underlying mechanisms of these DEmiRNAs, the target sequences of DEmiRNAs among the groups were predicted using Targetscan. We predicted the target genes for the upregulated and downregulated miRNAs. Next, the target genes were analyzed by GO and KEGG analyses.
GO analysis showed that most of the functional terms were related to the regulation of transcription from RNA polymerase II promoter and MAPK cascade. Figure 3c shows the top 10 downregulated biological processes, while Fig. 3d shows the top 10 upregulated biological processes.
KEGG analysis consistently revealed that most genes targeted by the down-regulated miRNAs were enriched in PI3K-Akt pathway, pathway in cancer, and focal adhesion (Fig. 3e), while the genes targeted by the up-regulated miRNAs were enriched in PI3K-Akt pathway, p53 pathway, and Hippo pathway (Fig. 3f). Among these, the PI3K-Akt pathway has been previously implicated in NPC progression [24, 25]. Our findings are in line with those of previous studies.
Network analysis of miRNA-mRNA pathway
To investigate the functions of miRNAs and the related target genes in NPC, we analyzed the interplay between the miRNAs and the relevant target genes which may help characterize the key regulatory mechanisms of miRNAs. All 1876 genes were predictably targeted by the down-regulated miRNAs, and 697 genes were predictably targeted by the up-regulated miRNAs. Target genes of the DEmiRNAs that were significantly enriched in the KEGG terms were annotated in the network (Fig. 4).
Important miRNAs related to NPC
Among the DEmiRNAs, we sought to identify the important miRNAs that are related to the progression of NPC. Consequently, we downloaded GSE118720 from the GEO database to study the miRNAs that are differentially expressed between NPC and normal samples. After analyzing the data (Fig. 5a), and then comparing with our data, we finally selected miRNA-1301-3p for further study. There was a significant difference between the two sets of data with respect to the expression of miRNA-1301-3p. We also compared the expression of miR-1301-3p between cancer and paracancerous samples in the context of different tumors in TCGA and drew a box diagram (Fig. 5b). The results showed a significant difference in the expressions of miR-1301-3p between cancer and paracancerous tissues in the context of many tumors.
Target prediction and functional analysis of hsa-miR-1301-3p
Microarray data (GSE12452, GSE13597, GSE53819, GSE64634) were obtained from the GEO database. RRA was applied to select significantly DEGs in the four datasets (Fig. 5c). Finally, 1406 DEGs were selected. Subsequently, the target genes of hsa-miR-1301-3p were predicted using TargetScan (http://www.targetscan.org/). To improve the reliability of bioinformatics analysis, the overlapping genes were selected for subsequent research (Fig. 5d).
GO and KEGG analyses of the overlapping genes were also performed. The significant GO category involved DNA synthesis in DNA repair, apoptotic process, meiotic nuclear division, and cell proliferation (Fig. 5e). The pathways of target genes that were enriched included MAPK pathway, cell cycle, and PI3K-Akt pathway (Fig. 5f).
RT-qPCR was conducted to measure the content of hsa-miR-1301-3p in plasma exosomes of 10 patients with early-stage NPC and normal controls; the results showed significant differences between the two groups (P < 0.05) (Fig. 6).
Early diagnosis of NPC is a key strategy for improving the prognosis of patients [25, 26]. However, accurate early diagnosis of NPC is typically challenging . Therefore, identification of novel diagnostic biomarkers of NPC is a key research imperative. The role of exosomes in carcinogenesis has evoked considerable interest in the context of several tumors such as hepatocellular carcinoma , breast cancer , and ovarian cancer . The present study is the first to explore potential plasma biomarkers in early-stage NPC exosomes. Exosomes were extracted from the plasma of patients with early-stage NPC at our hospital and healthy donors. We identified 21 downregulated miRNAs and 10 up-regulated miRNAs. The hsa-miR-1301-3p was regarded as a vital (upregulated) DEmiRNAs in early-stage NPC. Moreover, GO analysis revealed that the target genes of DEmiRNAs were mainly enriched in DNA synthesis and transcription regulation. KEGG analysis also showed that the target mRNAs of DEmiRNAs were primarily enriched in the PI3k-Akt and MAPK pathways. Our findings may offer novel insights into the role of exosomes in early-stage NPC.
To obtain the RNAs present in exosomes, exosomes from plasma of early-stage NPC were isolated and characterized. Several methods were used to demonstrate that the exosomes isolated from plasma were indeed exosomes and were pure enough for subsequent experiments. Tumor-derived exosomes represent a liquid biopsy tool with diagnostic and prognostic relevance; in addition, miRNAs present in exosomes may serve as potential biomarkers. Therefore, exosomes may provide novel insights for the early diagnosis of NPC.
We determined the types and contents of small RNAs in plasma samples obtained from patients with early NPC and healthy donors to identify the most abundant small RNAs in both groups. The results showed that miRNAs accounted for most small RNAs in these samples. The 20 most abundant miRNAs accounted for about 70% of all mapped reads. In recent years, the role of miRNAs in carcinogenesis and tumor progression has evoked increasing interest in the context of several cancers, such as lung carcinoma , colorectal carcinoma , and breast carcinoma . In a previous study, miR-99a was found to inhibit metastasis of NPC by targeting the HOXA1 gene . Jiang et al. found that EBV miRNA BART2-5p may promote the metastasis of NPC by suppressing the RND3 gene . Hence, miRNAs represent a focus area for further research.
Identification of the DEmiRNAs is a critical step for investigating the mechanisms by which miRNAs contribute to the initiation of NPC for exploring novel diagnostic biomarkers. In the present study, RNA-seq identified 21 down-regulated and 10 up-regulated miRNAs in NPC. Many of these are associated with cancer; for example, decreased expressions of has-miR-143 and has-miR-145 have also been identified in chronic lymphocytic leukemia , lung cancer [37, 38], and colorectal cancer , respectively. Daniela et al. demonstrated a positive correlation of miR-183 with cervical cancer . The details of the DEmiRNAs are provided in Fig. 3a and b. Accordingly, we selected these DEmiRNAs for further investigation.
Target genes of these DEmiRNAs were analyzed by GO and KEGG analyses. KEGG analysis showed that the target genes were significantly enriched in the PI3k-Akt and MAPK pathway. The enriched pathways discovered in the present study are mainly involved in cancer progression and metastasis. Several studies have shown that the PI3K-Akt pathway is crucial for the proliferation, migration, and invasion of NPC cells [41, 42]. Downregulation of the SPEN gene leads to the inactivation of PI3k-Akt and c-JUN signals, thus inhibiting the migration and invasion of NPC . Zhu et al. revealed that inhibition of autophagy by the ANXA1 gene may promote the migration and metastasis of NPC by activation of the PI3K-Akt pathway .
MAPK plays a vital role in cell proliferation, invasion, and metastasis. Overexpression of miR-4500 may downregulate the RRM2 gene and inhibit the activation of the MAPK pathway, thereby reducing breast carcinoma cell proliferation, migration, and angiogenesis . In a study by Wan et al., Aur-A was found to promote EMT and invasive growth of NPC through the downstream MAPK pathway . GO and KEGG analyses of target genes of has-miR-1301-3p were largely consistent with the above-mentioned results. Our results indicate that the PI3k-Akt and MAPK pathway play a crucial role in the genesis of NPC. However, the specific underlying mechanisms of the PI3k-Akt pathway and MAPK pathway need to be investigated in further study.
Studies have shown that miRNAs can match their target mRNAs and act on them to inhibit protein translation by degradation of mRNAs . We predicted the target mRNAs of DEmiRNAs by TargetScan. Subsequently, the target mRNAs that were significantly enriched in the KEGG terms were annotated in the network. Furthermore, miRNA-mRNA-pathway networks were constructed via an interplay analysis. Visual networks of NPC were constructed to explore the direct interactions between miRNAs and the corresponding target mRNAs and pathways in early-stage NPC, as well as to investigate the potential regulatory functions of the DEmiRNAs.
To identify important DEmiRNAs related to the oncogenesis of NPC for subsequent research, we analyzed our sequencing data and the GSE118720 dataset downloaded from the GEO database. The findings suggested significant upregulation of the hsa-miR-1301-3p in the plasma exosomes. The results are consistent with high expression of has-miR-1301-3p in the TCGA database. The association of has-miR-1301-3p with proliferation, invasion, and migration of cells has been found in the context of many tumors. In a study by Fang et al., has-miR-1301 was found to inhibit the migration and invasion of HepG2 cells and promote cell apoptosis . However, no studies have investigated the target genes and mechanism of hsa-miR-1301-3p in the context of early-stage NPC. Our study explored the underlying mechanism of hsa-miR-1301-3p using GO and KEGG analyses, which may facilitate a better understanding of the development and progression of early-stage NPC.
Some limitations of our study should be considered. First, the sample size of early-stage NPC patients and normal donors was relatively small. Therefore, larger studies are required to provide more robust evidence. In addition our study was based on bioinformatics analyses. Further experiments are required to provide direct evidence of the mechanisms in early-stage NPC. Furthermore, we did not analyze data pertaining to clinical prognosis of patients with early-stage NPC. Due to the lack of appropriate datasets, we could not investigate the correlation between the DEmiRNAs and the prognosis.
In conclusion, we identified 21 downregulated and 10 upregulated DEmiRNAs in early-stage NPC using bioinformatics analyses. The PI3k-Akt and MAPK pathways might play a vital role in the pathogenesis of NPC. Our findings may provide novel insights into the pathogenetic mechanisms of early-stage NPC and help identify underlying diagnostic biomarkers in plasma exosomes.
We identified several DEmiRNAs and analyzed the biological processes and signaling pathways of their predicted target genes in early-stage NPC using bioinformatics analyses. Our study may offer novel insights into the potential biomarkers and mechanisms of early-stage NPC. However, further functional experiments are required for in-depth exploration of the molecular mechanisms in early-stage NPC.
Availability of data and materials
Transmission electron microscopy
Principal component analysis
Kyoto encyclopedia of genes and genomes
Differentially expressed miRNAs
Gene Expression Omnibus
Latent membrane protein 1
The Cancer Genome Atlas
Nanoflow cytometry measurement
Phosphate buffer saline
Transcripts per million reads
Robust rank aggregation
Mitogen-activated protein kinase
Phosphatidylinositol-3 OH kinase/protein kinase B
Shield KD, Ferlay J, Jemal A, Sankaranarayanan R, Chaturvedi AK, Bray F, Soerjomataram I. The global incidence of lip, oral cavity, and pharyngeal cancers by subsite in 2012. CA Cancer J Clin. 2017;67:51–64. https://doi.org/10.3322/caac.21384.
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: 431 GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. Cancer J Clin. 2018;68(6):394–424. https://doi.org/10.3322/caac.21492.
Su SF, Han F, Zhao C, et al. Intensity-modulated radiotherapy prolongs the survival of patients with nasopharyngeal carcinoma compared with conventional two-dimensional radiotherapy: a 10-year experience with a large cohort and long follow-up. Eur J Cancer. 2015;51:2587–95. https://doi.org/10.1016/j.ejca.2015.08.006.
Zhang MX, Li J, Shen GP, et al. The seventh edition of the UICC/AJCC staging system for nasopharyngeal carcinoma is prognostically useful for patients treated with intensity-modulated radiotherapy from an endemic area in China. Radiother Oncol. 2012;104:331–7. https://doi.org/10.1016/j.radonc.2011.10.009.
Mao YP, Xie FY, Liu LZ, Sun Y, Li L, Tang LL, Liao XB, Xu HY, Chen L, Lai SZ, Lin AH, Liu MZ, Ma J. Re-evaluation of 6th edition of AJCC staging system for nasopharyngeal carcinoma and proposed improvement based on magnetic resonance imaging. Int J Radiat Oncol Biol Phys. 2009;73:1326–34. https://doi.org/10.1016/j.ijrobp.2008.07.062.
Pan JJ, Ng WT, Zong JF, Chan LL, O’Sullivan B, Lin SJ, Sze HC, Chen YB, Choi HC, Guo QJ, Kan WK, Xiao YP, Wei X, Le QT, Glastonbury CM, Colevas AD, Weber RS, Shah JP, Lee AW. Proposal for the 8th edition of the AJCC/UICC staging system for nasopharyngeal cancer in the era of intensity-modulated radiotherapy. Cancer. 2016;122:546–58. https://doi.org/10.1002/cncr.29795.
Peng H, Dong D, Fang MJ, Li L, Tang LL, Chen L, Li WF, Mao YP, Fan W, Liu LZ, Tian L, Lin AH, Sun Y, Tian J, Ma J. Prognostic value of deep learning PET/CT-based radiomics: potential role for future individual induction chemotherapy in advanced nasopharyngeal carcinoma. Clin Cancer Res. 2019. https://doi.org/10.1158/1078-0432.CCR-18-3065.
Lo K, Chung G, To K. Deciphering the molecular genetic basis of NPC through molecular, cytogenetic, and epigenetic approaches. Semin Cancer Biol. 2012;22(2):79–86. https://doi.org/10.1016/j.semcancer.2011.12.011.
Xu FH, Xiong D, Xu YF, Cao SM, Xue WQ, Qin HD, Liu WS, Cao JY, Zhang Y, Feng QS, Chen LZ, Li MZ, Liu ZW, Liu Q, Hong MH, Shugart YY, Zeng YX, Zeng MS, Jia WH. An epidemiological and molecular study of the relationship between smoking, risk of nasopharyngeal carcinoma, and Epstein-Barr virus activation. J Natl Cancer Inst. 2012;104:1396–410. https://doi.org/10.1093/jnci/djs320.
Guo SS, Tang LQ, Chen QY, Zhang L, Liu LT, Guo L, Mo HY, Luo DH, Huang PY, Xiang YQ, Sun R, Chen MY, Wang L, Lv X, Zhao C, Guo X, Cao KJ, Qian CN, Zeng MS, Bei JX, Hong MH, Shao JY, Sun Y, Ma J, Mai HQ. Inductionchemotherapy followed by concurrent chemoradiotherapy versus concurrent chemoradiotherapy alone in stage III-IVb nasopharyngeal carcinoma patients with Epstein-Barr virus DNA >/=4000 copies/ml: a matched study. Oncotarget. 2016;7:29739–48. https://doi.org/10.18632/oncotarget.8828.
Peng H, Chen L, Li WF, Guo R, Zhang Y, Zhang F, et al. Prognostic value of neoadjuvant chemotherapy in locoregionally advanced nasopharyngeal carcinoma with low pre-treatment Epstein-Barr virus DNA: a propensity-matched analysis. J Cancer. 2016;7:1465–71.
Le QT, Jones CD, Yau TK, Shirazi HA, Wong PH, Thomas EN, Patterson BK, Lee AW, Zehnder JL. A comparison study of different PCR assays in measuring circulating plasma epstein-barr virus DNA levels in patients with nasopharyngeal carcinoma. Clin Cancer Res. 2005;11:5700–7. https://doi.org/10.1158/1078-0432.CCR-05-0648.
Kowal J, Tkach M, Thery C. Biogenesis and secretion of exosomes. Curr Opin Cell Biol. 2014;29:116–25. https://doi.org/10.1016/j.ceb.2014.05.004.
Bobrie A, Thery C. Exosomes and communication between tumours and the immune system: are all exosomes equal? Biochem Soc T rans. 2013;41:263–7. https://doi.org/10.1042/BST20120245.
Nogues L, Benito-Martin A, Hergueta-Redondo M, Peinado H. The influence of tumour-derived extracellular vesicles on local and distal metastatic dissemination. Mol Asp Med. 2018;60:15–26. https://doi.org/10.1016/j.mam.2017.11.012.
Melo SA, Sugimoto H, O’Connell JT, Kato N, Villanueva A, Vidal A, Qiu L, Vitkin E, Perelman LT, Melo CA, Lucci A, Ivan C, Calin GA, Kalluri R. Cancer exosomes perform cell-independent microRNA biogenesis and promote tumorigenesis. Cancer Cell. 2014;26:707–21. https://doi.org/10.1016/j.ccell.2014.09.005.
Abd Elmageed ZY, Yang Y, Thomas R, Ranjan M, Mondal D, Moroz K, Fang Z, Rezk BM, Moparty K, Sikka SC, Sartor O, Abdel-Mageed AB. Neoplastic reprogramming of patient-derived adipose stem cells by prostate cancer cell-associated exosomes. Stem Cells. 2014;32:983–97. https://doi.org/10.1002/stem.1619.
Le MT, Hamar P, Guo C, Basar E, Perdigão-Henriques R, Balaj L, Lieberman J. miR-200-containing extracellular vesicles promote breast cancer cell metastasis. J Clin Invest. 2014;124:5109–28. https://doi.org/10.1172/JCI75695.
Meckes DG, Shair KH, Marquitz AR, Kung CP, Edwards RH, Raab-Traub N. Human tumor virus utilizes exosomes for intercellular communication. Proc Natl Acad Sci USA. 2010. https://doi.org/10.1073/pnas.1014194107.
Wan FZ, Chen KH, Sun YC, Chen XC, Liang RB, Chen L, Zhu XD. Exosomes overexpressing miR-34c inhibit malignant behavior and reverse the radioresistance of nasopharyngeal carcinoma. J Transl Med. 2020. https://doi.org/10.1186/s12967-019-02203-z.
Amin MBES, Greene FL, et al., editors. AJCC Cancer Staging Manual. 8th ed. New York: Springer; 2017.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015. https://doi.org/10.1093/nar/gkv007.
Kolde R, Laur S, Adler P, Vilo J. Robust rank aggregation for gene list integration and meta-analysis. Bioinformatics. 2012;28(4):573–80. https://doi.org/10.1093/bioinformatics/btr709.
Xuefang Z, Ruinian Z, Liji J, Chun Z, Qiaolan Z, Jun J, Yuming C, Junrong H. miR-331-3p inhibits proliferation and promotes apoptosis of nasopharyngeal carcinoma cells by targeting elf4B-PI3K-AKT pathway. Technol Cancer Res Treat. 2020. https://doi.org/10.1177/1533033819892251.
Huang Q, Hou S, Zhu X, Liu S. MicroRNA-192 promotes the development of nasopharyngeal carcinoma through targeting RB1 and activating PI3K/AKT pathway. World J Surg Oncol. 2020. https://doi.org/10.1186/s12957-020-1798-y.
Zhang MX, Li J, Shen GP, Zou X, Xu JJ, Jiang R, You R, Hua YJ, Sun Y, Ma J, Hong MH, Chen MY. Intensity-modulated radiotherapy prolongs the survival of patients with nasopharyngeal carcinoma compared with conventional two-dimensional radiotherapy: a 10-year experience with a large cohort and long follow-up. Eur J Cancer. 2015;51(17):2587–95. https://doi.org/10.1016/j.ejca.2015.08.006.
Chen L, Mao YP, Xie FY, Liu LZ, Sun Y, Tian L, Tang LL, Lin AH, Li L, Ma J. The seventh edition of the UICC/AJCC staging system for nasopharyngeal carcinoma is prognostically useful for patients treated with intensity-modulated radiotherapy from an endemic area in China. Radiother Oncol. 2012;104(3):331–7. https://doi.org/10.1016/j.radonc.2011.10.009.
Zhuo X, Zhou W, Li D, Chang A, Wang Y, Wu Y, Zhou Q. Plasma microRNA expression signature involving miR-548q, miR-630 and miR-940 as biomarkers for nasopharyngeal carcinoma detection. Cancer Biomark. 2018;23(4):579–87. https://doi.org/10.3233/CBM-181852.
Li X, Li C, Zhang L, Wu M, Cao K, Jiang F, Chen D, Li N, Li W. The significance of exosomes in the development and treatment of hepatocellular carcinoma. Mol Cancer. 2020;19(1):1. https://doi.org/10.1186/s12943-019-1085-0.
Jiang M, Zhang W, Zhang R, Liu P, Ye Y, Yu W, Guo X, Yu J. Cancer exosome-derived miR-9 and miR-181a promote the development of early-stage MDSCs via interfering with SOCS3 and PIAS3 respectively in breast cancer. Oncogene. 2020;39(24):4681–94. https://doi.org/10.1038/s41388-020-1322-4.
Hayashita Y, Osada H, Tatematsu Y, Yamada H, Yanagisawa K, Tomida S, Yatabe Y, Kawahara K, Sekido Y, Takahashi T. A polycistronic microRNA cluster, miR-17-92, is overexpressed in human lung cancers and enhances cell proliferation. Cancer Res. 2005;65(21):9628–32.
Nakajima G, Hayashi K, Xi Y, Kudo K, Uchida K, Takasaki K, Yamamoto M, Ju J. Non-coding MicroRNAs hsa-let-7 g and hsa-miR-181b are associated with chemoresponse to S-1 in colon cancer. Cancer Genom Proteom. 2006;3(5):317–24.
Xin F, Li M, Balch C, Thomson M, Fan M, Liu Y, Hammond SM, Kim S, Nephew KP. Computational analysis of microRNA profiles and their target genes suggests significant involvement in breast cancer antiestrogen resistance. Bioinformatics. 2009;25(4):430–4. https://doi.org/10.1093/bioinformatics/btn646.
Wang JG, Tang WP, Liao MC, Liu YP, Ai XH. MiR-99a suppresses cell invasion and metastasis in nasopharyngeal carcinoma through targeting HOXA1. Onco Targets Ther. 2017;10:753–61. https://doi.org/10.2147/OTT.S126781.
Jiang C, Li L, Xiang YQ, Lung ML, Zeng T, Lu J, Tsao SW, Zeng MS, Yun JP, Kwong DLW, Guan XY. Epstein-Barr virus miRNA BART2-5p promotes metastasis of nasopharyngeal carcinoma by suppressing RND3. Cancer Res. 2020;80(10):1957–69. https://doi.org/10.1158/0008-5472.CAN-19-0334.
Calin GA, Liu CG, Sevignani C, Ferracin M, Felli N, Dumitru CD, Shimizu M, Cimmino A, Zupo S, Dono M, Dell’Aquila ML, Alder H, Rassenti L, Kipps TJ, Bullrich F, Negrini M, Croce CM. MicroRNA profiling reveals distinct signatures in B cell chronic lymphocytic leukemias. Proc Natl Acad Sci USA. 2004;101(32):11755–60. https://doi.org/10.1073/pnas.0404432101.
Johnson SM, Grosshans H, Shingara J, Byrom M, Jarvis R, Cheng A, Labourier E, Reinert KL, Brown D, Slack FJ. RAS is regulated by the let-7 microRNA family. Cell. 2005;120(5):635–47. https://doi.org/10.1016/j.cell.2005.01.014.
Takamizawa J, Konishi H, Yanagisawa K, Tomida S, Osada H, Endoh H, Harano T, Yatabe Y, Nagino M, Nimura Y, Mitsudomi T, Takahashi T. Reduced expression of the let-7 microRNAs in human lung cancers in association with shortened postoperative survival. Cancer Res. 2004;64(11):3753–6. https://doi.org/10.1158/0008-5472.CAN-04-0637.
Michael MZ, O’Connor SM, van Holst Pellekaan NG, Young GP, James RJ. Reduced accumulation of specific microRNAs in colorectal neoplasia. Mol Cancer Res. 2003;1(12):882–91.
Witten D, Tibshirani R, Gu SG, Fire A, Lui WO. Ultra-high throughput sequencing-based small RNA discovery and discrete statistical biomarker analysis in a collection of cervical tumours and matched controls. BMC Biol. 2010;8:58. https://doi.org/10.1186/1741-7007-8-58.
Xu M, Tian GL, Hao CC, Shi M, Zha DJ, Liang K. MicroRNA-29 targets FGF2 and inhibits the proliferation, migration and invasion of nasopharyngeal carcinoma cells via PI3K/AKT signaling pathway. Eur Rev Med Pharmacol Sci. 2020;24(13):7199. https://doi.org/10.26355/eurrev_201906_18186.
Xie X, Xiong G, Chen W, Fu H, Li M, Cui X. FOXD3 inhibits cell proliferation, migration and invasion in nasopharyngeal carcinoma through regulation of the PI3K/Akt pathway. Biochem Cell Biol. 2020. https://doi.org/10.1139/bcb-2020-0011.
Li Y, Lv Y, Cheng C, Huang Y, Yang L, He J, Tao X, Hu Y, Ma Y, Su Y, Wu L, Yu G, Jiang Q, Liu S, Liu X, Liu Z. SPEN induces miR-4652-3p to target HIPK2 in nasopharyngeal carcinoma. Cell Death Dis. 2020;11(7):509. https://doi.org/10.1038/s41419-020-2699-2.
Zhu JF, Huang W, Yi HM, Xiao T, Li JY, Feng J, Yi H, Lu SS, Li XH, Lu RH, He QY, Xiao ZQ. Annexin A1-suppressed autophagy promotes nasopharyngeal carcinoma cell invasion and metastasis by PI3K/AKT signaling activation. Cell Death Dis. 2018;9(12):1154. https://doi.org/10.1038/s41419-018-1204-7.
Li S, Mai H, Zhu Y, Li G, Sun J, Li G, Liang B, Chen S. MicroRNA-4500 inhibits migration, invasion, and angiogenesis of breast cancer cells via RRM2-dependent MAPK signaling pathway. Mol Ther Nucleic Acids. 2020;21:278–89. https://doi.org/10.1016/j.omtn.2020.04.018.
Wan XB, Long ZJ, Yan M, Xu J, Xia LP, Liu L, Zhao Y, Huang XF, Wang XR, Zhu XF, Hong MH, Liu Q. Inhibition of Aurora-A suppresses epithelial-mesenchymal transition and invasion by downregulating MAPK in nasopharyngeal carcinoma cells. Carcinogenesis. 2008;29(10):1930–7. https://doi.org/10.1093/carcin/bgn176.
Shenouda SK, Alahari SK. MicroRNA function in cancer: oncogene or a tumor suppressor? Cancer Metastasis Rev. 2009;28(3–4):369–78. https://doi.org/10.1007/s10555-009-9188-5.
Fang L, Yang N, Ma J, Fu Y, Yang GS. microRNA-1301-mediated inhibition of tumorigenesis. Oncol Rep. 2012;27(4):929–34. https://doi.org/10.3892/or.2011.1589.
The Gene Expression Omnibus Database. http://www.ncbi.nlm.nih.gov/geo. Accessed 20 Dec 2019.
We would like to thank Medjaden Bioscience Limited (MJD) for the paper revision support.
This project was funded by Joints Funds for the Innovation of Science and Technology, Fujian Province, China (Grant numbers: 2018Y9107, 2019Y9034), and Science and Technology Program of Fujian Province, China (Grant number: 2018Y2003).
Ethics approval and consent to participate
The ethics committee of Fujian Cancer Hospital approved the human tissue samples related to this work (project ethics number: SQ2019-018-01).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Concentration of exosome. A, B: Concentration of exosomes as detected by Nanoparticle Tracking Analysis (NTA). Fig.S1A represents three NPC samples, while Fig.S1B represents three normal samples
The sequence reads produced from all samples.
The sequence reads of 12 samples were mapped to Rfam database.
: Table S3. The read counts of the 20 most abundant miRNAs.
About this article
Cite this article
Zheng, W., Ye, W., Wu, Z. et al. Identification of potential plasma biomarkers in early-stage nasopharyngeal carcinoma-derived exosomes based on RNA sequencing. Cancer Cell Int 21, 185 (2021). https://doi.org/10.1186/s12935-021-01881-4
- Nasopharyngeal carcinoma (NPC)
- Gene Expression Omnibus (GEO)
- Bioinformatics analysis
- Kyoto encyclopedia of genes and genomes (KEGG)