Integrated analysis of the differential cellular and EBV miRNA expression profiles in microdissected nasopharyngeal carcinoma and non-cancerous nasopharyngeal tissues

  • Authors:
    • Xun-Xun Wan
    • Hong Yi
    • Jia-Quan Qu
    • Qiu-Yan He
    • Zhi-Qiang Xiao
  • View Affiliations

  • Published online on: September 1, 2015     https://doi.org/10.3892/or.2015.4237
  • Pages: 2585-2601
Metrics: Total Views: 0 (Spandidos Publications: | PMC Statistics: )
Total PDF Downloads: 0 (Spandidos Publications: | PMC Statistics: )


Abstract

Nasopharyngeal carcinoma (NPC) is commonly diagnosed in southern Asia. MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression post-transcriptionally. Increasing evidence suggests that the dysregulation of miRNAs promotes NPC tumorigenesis. Epstein-Barr virus (EBV) infection and EBV-encoded miRNAs are also associated with the development of NPC. However, it is unclear how cellular and EBV miRNAs jointly regulate target genes and signaling pathways in NPC. In the present study, we analyzed the differential cellular and EBV miRNA expression profiles in 20 pooled NPC tissues using microarrays. We found that 19 cellular miRNAs and 9 EBV miRNAs were upregulated and 31 cellular miRNAs were downregulated in NPC tissues. Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated that the 19 upregulated miRNAs target mainly the p53 signaling pathway in cancer, whereas the downregulated miRNAs regulate pathways related to cancer, focal adhesion and Erb, and MAPK signaling. In contrast, the upregulated EBV miRNAs target primarily the TGF-β and Wnt signaling pathways. Data also suggested that cellular miR-34b, miR-34c, miR-18a, miR‑200a/b, miR-449a, miR-31 and let-7 may be dysregulated in NPCs, and that the aberrant activation of their target genes in the p53 pathway and cell cycle enhance NPC cell survival and proliferation. In addition, EBV-miRNAs such as BART3 and BART5 target genes in the p53, TGF-β and Wnt signaling pathways to modulate NPC apoptosis and transformation. To better elucidate the interaction between miRNAs and target genes, we constructed an anti-correlated cellular and EBV miRNA/target gene regulatory network. The current findings may help dissect the roles played by cellular and EBV miRNAs during NPC tumorigenesis, and also provide useful biomarkers for the diagnosis and treatment of NPCs.

Introduction

Nasopharyngeal carcinoma (NPC) is a human malignancy derived from the epithelium of the nasopharyngeal recess. It has a particularly high incidence in Southeast Asia, and is associated with genetic and epigenetic events, as well as Epstein-Barr virus (EBV) infection. It is important to elucidate the molecular mechanisms of NPC carcinogenesis.

MicroRNAs (miRNAs) are short non-coding RNA molecules and are involved in post-transcriptional (13). The dysregulation of miRNAs appears to play a crucial role in cancer, where they exert effects as oncogenes or tumor suppressors (4). Approximately 50% of miRNAs are localized in cancer-associated genomic regions, and their expression can be altered by genomic amplification, loss of heterozygosity, viral integration and genomic rearrangements (5). Large-scale profiling studies revealed a global alteration of the miRNA expression patterns in human types of cancers (68), and distinct miRNA expression signatures have been proposed as diagnostic and prognostic markers for various types of human cancer (9,10).

A growing number of cellular miRNAs are now known to be involved in NPC via targeting specific mRNAs. Most of these, including miR-29c, miR-218, miR-26a, miR-663, let-7, miR-141, miR-216b and miR-200a, play an important role in the onset and progression of NPC (1118). A unique feature of NPC is its strong association with EBV latent infection. The full-length EBV genome is detected in all cases of NPC, suggesting that products of the EBV genome are involved in the pathogenesis of this malignancy (19,20). As well as EBV-encoded protein-coding genes, latently infected NPC cells and tissues also express high levels of viral miRNAs (21). EBV is the first human virus known to express miRNAs. To date, >44 EBV miRNAs have been discovered in EBV-positive cell lines, and they exhibit different expression patterns depending on the sample and stage of infection (21). Since EBV severely dysregulates the mRNA profiles of host cells (22), EBV miRNAs may have oncogenic properties. Functional studies revealed that several EBV-encoded miRNAs play roles in NPC by targeting a variety of cellular (2326) and viral (2730) mRNAs. Therefore, both cellular and EBV miRNAs play a crucial role in NPC (31).

Various genome-wide miRNA expression profiling studies using microarray-based approaches yielded abundant information regarding the phenotypes of cancers. Distinct patterns of cellular miRNA expression and specific miRNA signatures in NPC have been identified, which are associated with clinicopathological characteristics as well as prognosis, suggesting that miRNAs play key roles in the development, invasion and metastasis of NPC (3236). Several studies have performed EBV miRNA expression profiling in NPC cells and tissues (21,3739). However, there has been no systematic analysis of both differentially expressed cellular and EBV miRNA expression profiles in NPC tissues.

In the present study, we used microarray analyses to compare the differences in the cellular and EBV miRNA expression profiles in pooled laser capture microdissected (LCM) NPC tissues and normal nasopharyngeal epithelial tissues (NNETs). We identified the target genes of differential cellular and EBV miRNAs using database predictions, and comparing the resulting information with the gene expression in NPC reported in the Gene Expression Omnibus (GEO) (GSE34573 and GSE12452). In addition, we analyzed the functions and pathways modulated by the cellular and EBV miRNA target genes. Finally, we constructed an anti-correlated cellular and EBV miRNA/target gene regulatory network. The resulting data may be helpful for elucidating the roles of cellular and EBV miRNA deregulation in NPC carcinogenesis, and also revealing the global trends in cellular and EBV miRNAs and gene interactions in NPC.

Materials and methods

Sample collection, laser capture microdissection and total RNA extraction

Paired NPC tissues and adjacent NNETs from 20 patients with confirmed NPC were obtained from the Tumor Hospital of Hunan Province (China) at the time of diagnosis before any treatment, and were used for microarrays and quantitative reverse transcription-PCR (qRT-PCR). Tissue specimens were obtained using fiber optic nasopharyngoscopy directly at the tumor growth site, and the adjacent side with an observed normal mucosal morphology was used for the corresponding non-tumor pair. All patients signed an informed consent form for participation in the present study, which had been previously reviewed by the Institutional Review Board. All tissue samples were verified by histopathology before microdissection. Laser capture microdissection (LCM) was used to purify the target cells from the fresh NPC tissues and NNETs using a Leica AS LMD system (Leica, Mannheim, Germany), as previously described (40). The cells obtained from LCM were collected in the cap of a 0.5-ml Eppendorf (EP) tube containing TRIzol reagent (Invitrogen, Life Technologies, Shanghai, China). Total RNA was extracted using TRIzol according to the manufacturer's instructions. The concentration and integrity of the resulting RNA were evaluated using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). Only total RNA samples with an RNA integrity number (RIN) ≥7 were used for miRNA microarray analyses and qRT-PCR. To diminish the effect of biological sample variation on the results of the microarray analysis, equal amounts of total RNA from the microdissected cells of 10 different individuals were pooled to generate one common sample for each tissue type.

miRNA microarrays

RNA samples were labeled using a miRCURY Hy3/Hy5 labeling kit and hybridized on the miRCURY locked nucleic acid (LNA) array (version 10.0; Exiqon, KangChen Biotech, Shanghai, China) according to the manufacturer's instructions. The array contains probes for 850 human miRNAs and 39 EBV miRNAs, and uses a one-color technique. The arrays were scanned using a laser scanner (GenePix 4000B; Molecular Devices, LLC, Sunnyvale, CA, USA), and the raw data were normalized and analyzed using image analysis software (GenePix Pro v6.0) to generate a median value for each miRNA from the repeated probes. The miRNA signal intensities were log2 transformed, and differentially expressed miRNAs were analyzed using the significance analysis of microarrays (SAM, version 3.01). P-values were calculated using t-tests, and differentially detected miRNAs with a ≥2.0- or ≤0.50 fold-change and a P-value <0.01 were considered to be statistically significant.

Quantitative RT-PCR

To validate the miRNA expression levels determined using microarrays, stem-looped qRT-PCR was performed to detect the expression of 15 cellular miRNAs (hsa-miR-200b, hsa-miR-34c-5p, hsa-miR-31, hsa-miR-34b*, hsa-miR-656, hsa-miR-374a, hsa-miR-150, hsa-miR-200a, hsa-miR-365-1, hsa-miR-15a, hsa-miR-20a, hsa-miR-17, hsa-miR-18a, hsa-miR-106a and hsa-miR-412) and five EBV miRNAs (ebv-miR-BART5, ebv-miR-BART6-3p, ebv-miR-BART3, ebv-miR-BART1-5p and ebv-miR-BART7) in an independent set of samples consisting of 10 pairs of LCM-purified NPCs and NNETs. Briefly, 2 µg of total RNA was reverse transcribed into cDNA using a reverse transcription kit according to the manufacturer's protocol (Promega, Madison, WI, USA). miRNA-specific primers (Bulge-Loop™ miRNA qPCR primers) for the differentially expressed miRNAs were synthesized by GeneChem Co. (Shanghai, China) and are summarized in Table I. The RT products were amplified using real-time PCR with the miScript SYBR-Green PCR kit (Qiagen, Shanghai, China) according to the manufacturer's instructions; U6 was used as the internal control. The following PCR program was used: 95°C for 5 min, followed by 40 cycles of 95°C for 10 sec, 60°C for 20 sec and 78°C for 20 sec. The real-time PCR primers used to amplify the miRNAs were designed based on the miRNA sequences provided by the Sanger Center miRNA Registry and synthesized by GeneChem; their sequences are provided in Table II. All qRT-PCR reactions were performed in triplicate and products were quantified using the 2−ΔΔCt method normalized against U6. All qRT-PCR reactions were performed on an ABI GeneAmp PCR System 9700 (ABI, Carlsbad, CA, USA).

Table I

The specific stem-looped RT primers to representatives of differentially expressed miRNAs.

Table I

The specific stem-looped RT primers to representatives of differentially expressed miRNAs.

miRNART primers
hsa-miR-200b 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACTCATCAT-3′
hsa-miR-34c-5p 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACGCAATC-3′
hsa-miR-31 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACCAGCTAT-3′
hsa-miR-34b* 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACCAATCA-3′
hsa-miR-656 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACAGAGGTT-3′
hsa-miR-374a 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACCACTTAT-3′
hsa-miR-200a 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACACATCGT-3′
hsa-miR-150 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACCACTGGT-3′
hsa-miR-365-1 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACATAAGGA-3′
hsa-miR-15a 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACCACAAAC-3′
hsa-miR-412 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACACGGCT-3′
hsa-miR-20a 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACCTACCTG-3′
hsa-miR-18a 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACCTATCTG-3′
hsa-miR-17 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACCTACCTG-3′
hsa-miR-106a 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACCTACCTG-3′
ebv-miR-BART7 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACCCCTGG-3′
ebv-miR-BART1-5p 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACCACAGCA-3′
ebv-miR-BART3 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACACACCT-3′
ebv-miR-BART6-3p 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACTCTAAGG-3′
ebv-miR-BART5 5′-GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACCGATGGG-3′
U6 5′-CGCTTCACGAATTTGCGTGTCAT-3′

[i] RT, reverse transcription.

Table II

The real-time PCR primers for amplification of the mature miRNAs.

Table II

The real-time PCR primers for amplification of the mature miRNAs.

miRNABidirectional primersAnnealing temperature (°C)Product length (bp)
hsa-miR-200bGSP 5′-GGGGTAATACTGCCTGGT-3′
R 5′-TGCGTGTCGTGGAGTC-3′6063
hsa-miR-34c-5pGSP 5′-GGGAGGCAGTGTAGTTAGC-3′
R 5′-CAGTGCGTGTCGTGGAGT-3′6066
hsa-miR-31GSP 5′-GGAGGCAAGATGCTGGC-3′
R 5′-CAGTGCGTGTCGTGGAGT-3′6064
hsa-miR-34b*GSP 5′-GGGTAGGCAGTGTCATTAGC-3′
R 5′-CAGTGCGTGTCGTGGAGT-3′6066
hsa-miR-656GSP 5′-GGCGGAATATTATACAGTCAA-3′
R 5′-CAGTGCGTGTCGTGGAGT-3′6066
hsa-miR-374aGSP 5′-GGCACCTTATAATACAACCTG-3′
R 5′-TGCGTGTCGTGGAGTC-3′6065
hsa-miR-200aGSP 5′-GGGGTAACACTGTCTGGTAG-3′
R 5′-TGCGTGTCGTGGAGTC-3′6063
hsa-miR-150GSP 5′-GCTCTCCCAACCCTTGT-3′
R 5′-TGCGTGTCGTGGAGTC-3′6061
hsa-miR-365-1GSP 5′-GCACTTCACCTGGTCCACT-3′
R 5′-CAGTGCGTGTCGTGGAGT-3′6065
hsa-miR-15aGSP 5′-GGGTAGCAGCACATAATGG-3′
R 5′-CAGTGCGTGTCGTGGAGT-3′6067
hsa-miR-412GSP 5′-GCACTTCACCTGGTCCACT-3′
R 5′-CAGTGCGTGTCGTGGAGT-3′6065
hsa-miR-20aGSP 5′-GGGTAAAGTGCTTATAGTGC-3′
R 5′-TGCGTGTCGTGGAGTC-3′6063
hsa-miR-18aGSP 5′-GGGTAAGGTGCATCTAGTGC-3′
R 5′-TGCGTGTCGTGGAGTC-3′6064
hsa-miR-17GSP 5′-GGGCAAAGTGCTTACAGTGC-3′
R 5′-TGCGTGTCGTGGAGTC-3′6065
hsa-miR-106aGSP 5′-GGGGAAAAGTGCTTACAGTG-3′
R 5′-CAGTGCGTGTCGTGGAGT-3′6067
ebv-miR-BART7GSP 5′-GGGCATCATAGTCCAGTGT-3′
R 5′-CAGTGCGTGTCGTGGAGT-3′6065
ebv-miR-BART1-5pGSP 5′-GGCTCTTAGTGGAAGTGACG-3′
R 5′-CAGTGCGTGTCGTGGAGT-3′6067
ebv-miR-BART3GSP 5′-GGCGCACCACTAGTCACC-3′
R 5′-CAGTGCGTGTCGTGGAGT-3′6064
ebv-miR-BART6-3pGSP 5′-AAGCGGGGATCGGACTA-3′
R 5′-CAGTGCGTGTCGTGGAGT-3′6065
ebv-miR-BART5GSP 5′-CGGCAAGGTGAATATAGC-3′
R 5′-TGCGTGTCGTGGAGTC-3′6064
U6F 5′-GCTTCGGCAGCACATATACTAAAAT-3′
R 5′-CGCTTCACGAATTTGCGTGTCAT-3′6089

[i] GSP are specific primers corresponding to the miRNA. R are primters matching with RT primers.

Identification of cellular and EBV miRNA target genes

The putative target genes of the cellular miRNAs expressed differentially in miRNA microarrays were predicted using three databases (miRanda, PITA and TargetScan, 2013). The cellular miRNA target genes in any of the three databases were selected, and their potential differential gene expression profiles (GSE34573 and GSE12452) in NPC biopsies from the National Center for Biotechnology Information (NCBI)-Gene Expression Omnibus (GEO) were compared. Potential target genes, which were differentially expressed in the gene expression profiles in NPC and also had a negative correlation with the cellular miRNA expression in the microarrays, were identified and further analyzed.

The putative host target genes of the differentially expressed EBV miRNAs in the microarrays were predicted using three databases (RepTar (41), DIANA-microT v3.0 (42) and starBase (43), 2013). The host target genes recorded in any of these databases were selected and their potential differential gene expression profiles (GSE34573 and GSE12452) in NPC biopsies from the NCBI-GEO were compared. The predicted host target genes that were correlated negatively with the EBV miRNA expression profiles were identified and further analyzed. The putative viral target genes of the differentially expressed EBV miRNAs in the microarrays were predicted using the vHoT database (2013) (44), which predicts the viral target genes of virus-derived miRNAs using popular miRNA target prediction algorithms including TargetScan, miRanda and PITA. The targets of the differentially expressed EBV miRNAs in the viral genome were identified according to the match sites of these miRNAs in the EBV genome and the following criteria: a total context score 0 (TargetScan), a minimum free energy (MFE) ≤0 kcal/mol (miRanda) and a ddG threshold ≤10 (PITA).

GO and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of the cellular and EBV miRNA target genes

To understand the functions of the predicted host and virus target genes of the cellular and EBV miRNAs, we performed ontology classification and KEGG pathway analysis using DAVID Bioinformatics Resources 6.7 (http://david.abcc.ncifcrf.gov/home.jsp). The target genes were assigned to functional groups based on biological processes. Correlation was defined as significantly enriched when P≤0.05.

Integrative analysis of the cellular miRNA and gene expression profiles to construct an miRNA/gene regulatory network

The post-transcriptional regulatory network of an miRNA and gene is defined as a directed, bipartite graph in which the expression of the miRNA/target gene interacting pairs are anti-correlated. The differentially expressed miRNAs and genes of interest were determined using pathways extracted from KEGG as the primary nodal networks. The networks were drawn using Cytoscape 2.8.3. The resulting post-transcriptional network identified the genes regulated by the same miRNAs, as well as miRNAs that regulate the specific genes of functional relevance.

Establishing an intercross regulatory network of cellular and EBV miRNAs in NPC

To understand the interaction between cellular and EBV miRNAs, we identified the target genes regulated by both cellular and viral miRNAs that were anti-correlated with the miRNA expression observed in the microarrays. We then established a gene regulatory network that was co-regulated by both cellular and viral miRNAs in NPC. The networks were drawn using Cytoscape 2.8.3.

Results

Differentially expressed cellular and EBV miRNAs in NPC tissues and NNETs

miRNA expression profiling was performed using pooled LCM-purified NPC tissues and paired NNETs (Fig. 1), and miRCURY LNA arrays (version 10.0) containing 850 human and 39 EBV miRNAs. miRNAs were scored as positive using SAM; differential expression with a ≥2 fold-change was regarded as significant. These analyses identified 50 cellular miRNAs and 9 differential EBV miRNAs that were expressed differentially between NPC and normal tissues. Among these, 19 cellular miRNAs and 9 EBV miRNAs were upregulated, whereas 31 cellular miRNAs were downregulated in the NPC tissues compared with the NNETs (Table III). Importantly, miR-106a, miR-142-3p, miR-155, miR-17, miR-18a, let-7, miR-200a, miR-200b, miR-31, miR-34b, miR-34c and miR449a were previously reported to be differentially regulated (11,32,33,36,39).

Table III

Microarray-based detection of the differentially expressed miRNAs in the NPC tissues and the NNETs.

Table III

Microarray-based detection of the differentially expressed miRNAs in the NPC tissues and the NNETs.

miRNAChromosome locationNormalized value
Fold-change NPC/NNET
NNETNPC
Upregulated cellular miRNAs
 hsa-miR-106aXq26.20.080.404.77
 hsa-miR-142-3p17q220.200.402.00
 hsa-miR-142-5p17q220.110.222.00
 hsa-miR-1435q320.050.153.22
 hsa-miR-15521q21.30.150.352.35
 hsa-miR-1613q14.20.120.282.30
 hsa-miR-1713q31.30.110.363.15
 hsa-miR-18a13q31.30.030.154.33
 hsa-miR-193a-3p17q11.20.030.134.13
 hsa-miR-2066p12.20.230.552.35
 hsa-miR-20a13q31.30.961.962.04
 hsa-miR-2117q23.10.280.582.04
 hsa-miR-221*Xp11.30.631.432.27
 hsa-miR-30b*8q24.225.3214.142.66
 hsa-miR-30c-2*6q137.3317.962.45
 hsa-miR-338-5p17q25.30.180.925.17
 hsa-miR-41214q32.310.361.143.16
 hsa-miR-5507p14.319.2345.502.37
 hsa-miR-63719p13.30.170.362.06
Downregulated cellular miRNAs
 hsa-let-7b*22q13.311.260.040.04
 hsa-miR-120q13.332.001.030.51
 hsa-miR-105*Xq280.190.090.50
 hsa-miR-133b6p12.23.690.030.01
 hsa-miR-138-1*3p21.326.913.130.45
 hsa-miR-15019q13.332.200.650.30
 hsa-miR-15a13q14.21.930.350.18
 hsa-miR-15b3q25.331.980.160.08
 hsa-miR-19517p13.10.240.010.04
 hsa-miR-200a1p36.330.350.170.47
 hsa-miR-200b1p36.330.420.160.38
 hsa-miR-200c*12p13.311.330.420.31
 hsa-miR-319p21.30.120.010.12
 hsa-miR-32611q13.40.870.140.16
 hsa-miR-32816q22.10.290.140.49
 hsa-miR-3357q32.20.150.040.23
 hsa-miR-34b*11q23.10.630.000.01
 hsa-miR-34c-5p11q23.11.810.100.06
 hsa-miR-365 16p13.12/17q11.21.510.250.16
 hsa-miR-374aXq13.23.270.050.02
 hsa-miR-449a5q11.20.630.020.04
 hsa-miR-450b-3pXq26.30.130.030.25
 hsa-miR-486-5p8p11.211.900.790.42
 hsa-miR-506Xq27.30.360.140.38
 hsa-miR-513-5pXq27.312.494.950.40
 hsa-miR-518a-3p19q13.425.232.680.51
 hsa-miR-574-5p4p142.190.610.28
 hsa-miR-64720q13.330.860.370.43
 hsa-miR-65614q32.311.590.090.06
 hsa-miR-886-3p5q31.10.410.180.44
 hsa-miR-886-5p5q31.10.700.300.42
Upregulated EBV miRNA
 ebv-miR-BART8*0.080.283.37
 ebv-miR-BART40.040.123.47
 ebv-miR-BART50.120.524.45
 ebv-miR-BART1-5p0.060.254.19
 ebv-miR-BART18-3p0.140.282.02
 ebv-miR-BART100.000.22110.83
 ebv-miR-BART6-3p0.020.177.22
 ebv-miR-BART30.0040.64160.19
 ebv-miR-BART70.0040.44109.9

[i] EBV, Epstein-Barr virus; NNET normal nasopharyngeal epithelial tissue; NPC, nasopharyngeal carcinoma.

Compared with all previous published studies, several novel miRNAs that were upregulated in NPC were identified in the present study: hsa-miR-142-5p, hsa-miR-193a-3p, hsa-miR-221*, hsa-miR-30b*, hsa-miR-30c-2*, hsa-miR-338-5p, hsa-miR-412, hsa-miR-550 and hsa-miR-637. In addition several novel miRNAs that were downregulated in NPC were identified: hsa-miR-105*, hsa-miR-133b, hsa-miR-138-1*, hsa-miR-200c*, hsa-miR-326, hsa-miR-328, hsa-miR-365, hsa-miR-374a, hsa-miR-450b-3p, hsa-miR-506, hsa-miR-513-5p, hsa-miR-518a-3p, hsa-miR-574-5p, hsa-miR-647, hsa-miR-656, hsa-miR-886-3p and hsa-miR-886-5p.

Validation of the miRNA microarray results using qRT-PCR

To validate the microarray data, qRT-PCR was performed to quantify the expression of 15 cellular miRNAs and 5 EBV miRNAs in an independent set of 10 LCM-purified NPCs and paired NNETs. Data revealed that the expression of hsa-miR-200b, hsa-miR-34c-5p, hsa-miR-31, hsa-miR-34b*, hsa-miR-656, hsa-miR-374a, hsa-miR-150, hsa-miR-200a, hsa-miR-365-1 and hsa-miR-15a was significantly decreased, whereas hsa-miR-20a, hsa-miR-17, hsa-miR-18a, hsa-miR-106a, hsa-miR-412, ebv-miR-BART5, ebv-miR-BART6-3p, ebv-miR-BART3, ebv-miR-BART1-5p and ebv-miR-BART7 were significantly increased in the NPC tissues when compared with the NNETs. These observations were consistent with the results of the miRNA microarrays (Fig. 2).

Prediction of the miRNA target genes

The putative target genes of 50 differentially expressed cellular miRNAs were predicted using the three databases as described in Materials and methods. A total of 69,869 target genes were recorded by the databases (data not shown). To identify more genuine target genes, all the predicted cellular miRNA targets were compared with the differential gene expression profiles in NPC tissues within NCBI-GEO. This identified 555 cellular miRNA target genes that were expressed differentially in NPC and negatively correlated with the expression of the differentially expressed miRNAs. Among these target genes, 218 and 337 were regulated by the upregulated and downregulated miRNAs, respectively (data not shown).

The putative host target genes of the nine upregulated EBV miRNAs were predicted using the three databases as described in the Materials and methods. This identified 18,993 host target genes (data not shown), among which 636 were regulated by 4 or more EBV miRNAs (data not shown).

The putative viral target genes of the 9 differentially expressed EBV miRNAs were predicted using vHoT database. This revealed 68 target genes within the EBV genome for 8 differentially expressed EBV miRNAs using any of the three algorithms. Among these, 37 viral target genes were predicted at least by two algorithms (Table IV). Notably, there were no predicted results for ebv-miR-BART18-3p.

Table IV

The 37 putative viral target genes of the 8 upregulated EBV miRNAs predicted by at least two of the algorithms TargetScan, miRanda and PITA from the vHoT database.

Table IV

The 37 putative viral target genes of the 8 upregulated EBV miRNAs predicted by at least two of the algorithms TargetScan, miRanda and PITA from the vHoT database.

CategoryTarget geneEBV miRNAsGene description
CapsidBBRF1ebv-miR-BART10, ebv-miR-BART1-5p, ebv-miR-BART4, ebv-miR-BART8*Transcription factor, portal protein
CapsidBcLF1ebv-miR-BART1-5p, ebv-miR-BART3, ebv-miR-BART4MCP, VCA, viral lytic gene (late)
CapsidBDLF1ebv-miR-BART7, ebv-miR-BART3, ebv-miR-BART5, ebv-miR-BART1-5pMCP
CapsidBTRF1ebv-miR-BART10, ebv-miR-BART7, ebv-miR-BART6-3pCapsid
CapsidBVRF2ebv-miR-BART6-3p, ebv-miR-BART3, ebv-miR-BART8*Maturational protease
CapsidBdRF1ebv-miR-BART3, ebv-miR-BART4, ebv-miR-BART5Packaging, terminase small subunit, internal scaffold protein
GlycoproteinBALF4ebv-miR-BART3, ebv-miR-BART1-5p, ebv-miR-BART7 Glycoprotein-encoding genes, product gB
GlycoproteinBDLF3ebv-miR-BART10, ebv-miR-BART1-5p, ebv-miR-BART3 Glycoprotein-encoding genes, product gp150
GlycoproteinBILF1ebv-miR-BART4, ebv-miR-BART10, ebv-miR-BART1-5p, ebv-miR-BART5, ebv-miR-BART7Product gp64, undetected in mature virus
GlycoproteinBILF2ebv-miR-BART3, ebv-miR-BART5, ebv-miR-BART6-3p, ebv-miR-BART7Product gp78
GlycoproteinBLLF1ebv-miR-BART3, ebv-miR-BART1-5p, ebv-miR-BART7 Glycoprotein-encoding genes, product gp350
GlycoproteinBXLF2ebv-miR-BART1-5p, ebv-miR-BART4, ebv-miR-BART7Product gH
LatentEBNA-1ebv-miR-BART1-5p, ebv-miR-BART4, ebv-miR-BART3, ebv-miR-BART6-3p, ebv-miR-BART10, ebv-miR-BART7, ebv-miR-BART8*EBNA-1
LatentEBNA-LPebv-miR-BART5, ebv-miR-BART6-3p, ebv-miR-BART1-5p
LatentLMP-1ebv-miR-BART10, ebv-miR-BART8*LMP-1
LatentLMP-2Aebv-miR-BART1-5p, ebv-miR-BART4, ebv-miR-BART10, ebv-miR-BART5LMP-2A
LatentLMP-2B ebv-miR-BART8*, ebv-miR-BART1-5p, ebv-miR-BART7
LatentRPMS1ebv-miR-BART1-5p, ebv-miR-BART3, ebv-miR-BART4, ebv-miR-BART6-3p, ebv-miR-BART7, ebv-miR-BART8*, ebv-miR-BART10, ebv-miR-BART5
ReplicationBALF5ebv-miR-BART10, ebv-miR-BART1-5p, ebv-miR-BART4, ebv-miR-BART6-3p, ebv-miR-BART7Polymerase
ReplicationBBLF2/BBLF3 ebv-miR-BART6-3pPrimase-associated factor
TegumentBALF2ebv-miR-BART6-3p, ebv-miR-BART1-5p, ebv-miR-BART3, ebv-miR-BART4, ebv-miR-BART5, ebv-miR-BART8*Replication, viral lytic gene (EA), ssDNABP
TegumentBFRF1Aebv-miR-BART7Form the terminase complex, cleavage of DNA concatemers
TegumentBGLF1ebv-miR-BART3, ebv-miR-BART5, ebv-miR-BART7Capsid-assoc, packaging
TegumentBGLF3ebv-miR-BART5, ebv-miR-BART10, ebv-miR-BART4Undetected in mature virus
TegumentBNRF1ebv-miR-BART3, ebv-miR-BART5, ebv-miR-BART7, ebv-miR-BART4MTP
TegumentBOLF1ebv-miR-BART7, ebv-miR-BART3, ebv-miR-BART5, ebv-miR-BART6-3pLTPBP
TegumentBORF2ebv-miR-BART7, ebv-miR-BART3, ebv-miR-BART5, ebv-miR-BART8*Ribonucleotide reductase large subunit, RNR-L
TegumentBPLF1ebv-miR-BART3, ebv-miR-BART6-3p, ebv-miR-BART7, ebv-miR-BART8*, ebv-miR-BART10, ebv-miR-BART1-5p, ebv-miR-BART4, ebv-miR-BART5LTP
TegumentBRLF1ebv-miR-BART5, ebv-miR-BART7Product Rta
TegumentBRRF2ebv-miR-BART5, ebv-miR-BART6-3pLikely undergoes post-translational modification, including phosphorylation
TegumentBSLF1ebv-miR-BART7, ebv-miR-BART8*Virally encoded replication proteins, primase
TegumentBVRF1ebv-miR-BART6-3p, ebv-miR-BART1-5p, ebv-miR-BART4Capsid-assoc, portal plug
Transcription factorBHRF1ebv-miR-BART3, ebv-miR-BART8*Bcl-2 homolog
Transcription factorBZLF1ebv-miR-BART10, ebv-miR-BART3, ebv-miR-BART4Zta, viral lytic genes (immediate early, ZEBRA) site-specific DNA binding protein, transcriptional activator, origin binding protein, A BZLF1 mutation affects global expression of EBV lytic genes
?BDLF4ebv-miR-BART4, ebv-miR-BART8*
?BLLF2ebv-miR-BART3
?BWRF1 ebv-miR-BART6-3p

[i] EA, early antigen; EBV, Epstein-Barr virus; LTP, large tegument protein; LTPBP, LTP-binding protein; MCP, major capsid protein; MTP, major tegument protein; ssDNABP, single-stranded DNA-binding protein; VCA, viral capsid antigen.

Gene Ontology (GO) function and KEGG pathway analyses of the miRNA target genes

A total of 218 and 337 target genes whose expression was regulated by upregulated and downregulated cellular miRNAs, respectively, and negatively correlated with the expression of miRNAs in NPC were uploaded into DAVID for functional and pathway enrichment analyses. The results showed that 118 and 447 GO functions were annotated for upregulated and downregulated cellular miRNA targets (P≤0.05), respectively. Of these, 2 and 35 pathways were enriched statistically (P≤0.05) among the upregulated and the downregulated cellular miRNA targets, respectively.

In addition, the 636 targets that were regulated by 4 or more EBV miRNAs were also uploaded into DAVID for functional and pathway enrichment analyses. Data revealed that 101 GO functions and 12 pathways were statistically enriched (P≤0.05). The most enriched GO terms and KEGG pathways of the target genes are shown in Table V, most of which are involved in carcinogenesis. Notably, the novel differential miRNAs were all enriched in major cellular functions associated with tumorigenesis, including the regulation of apoptosis, cell proliferation, the cell cycle and cell differentiation (data not shown).

Table V

GO functional annotation and KEGG pathway analysis of the target genes regulated by the upregulated and downregulated cellular miRNAs and EBV miRNAs in the NPC biopsies.

Table V

GO functional annotation and KEGG pathway analysis of the target genes regulated by the upregulated and downregulated cellular miRNAs and EBV miRNAs in the NPC biopsies.

TermCount%P-valueFDR
GO functional annotation of upregulated cellular miRNA target genes
 GO:0043065-positive regulation of apoptosis200.742.78E-064.62E−03
 GO:0006917-induction of apoptosis160.591.64E-052.74E−02
 GO:0008285-negative regulation of cell proliferation170.631.71E-052.84E−02
 GO:0006915-apoptosis220.823.08E-055.13E−02
 GO:0042981-regulation of apoptosis260.963.56E-055.92E−02
 GO:0007242-intracellular signaling cascade331.221.29E-042.15E−01
 GO:0042127-regulation of cell proliferation230.854.85E-048.04E−01
 GO:0019220-regulation of phosphate metabolic process170.635.26E-048.72E−01
 GO:0006959-humoral immune response70.265.31E-048.80E−01
 GO:0042744-hydrogen peroxide catabolic process40.151.25E-032.05E+00
 GO:0045859-regulation of protein kinase activity130.481.70E-032.78E+00
 GO:0002252-immune effector process80.301.72E-032.83E+00
 GO:0007264-small GTPase mediated signal transduction120.441.98E-033.25E+00
 GO:0043066-negative regulation of apoptosis130.482.10E-033.44E+00
 GO:0034599-cellular response to oxidative stress50.192.21E-033.61E+00
GO functional annotation of downregulated cellular miRNA target genes
 GO:0042127-regulation of cell proliferation641.461.53E-192.69E−16
 GO:0008284-positive regulation of cell proliferation380.874.54E-137.99E−10
 GO:0043067-regulation of programmed cell death541.239.94E-131.75E−09
 GO:0051726-regulation of cell cycle320.731.23E-112.16E−08
 GO:0045449-regulation of transcription1042.383.94E-106.92E−07
 GO:0007242-intracellular signaling cascade631.441.28E-092.26E−06
 GO:0010604-positive regulation of macromolecule metabolic process491.122.53E-094.46E−06
 GO:0006357-regulation of transcription from441.013.83E-096.73E−06
 RNA polymerase II promoter
 GO:0007264-small GTPase mediated signal transduction270.624.49E-097.90E−06
 GO:0043065-positive regulation of apoptosis320.738.02E-091.41E−05
 GO:0007049-cell cycle451.038.86E-091.56E−05
 GO:0045596-negative regulation of cell differentiation220.501.49E-082.62E−05
 GO:0006355-regulation of transcription, DNA-dependent761.741.62E-082.84E−05
 GO:0009890-negative regulation of biosynthetic process370.851.76E-083.09E−05
 GO:0045786-negative regulation of cell cycle140.322.45E-084.30E−05
GO functional annotation of targets regulated by four or more EBV miRNAs
 GO:0010557-positive regulation of macromolecule biosynthetic process440.731.43E-052.47E−02
 GO:0010628-positive regulation of gene expression400.662.30E-053.99E−02
 GO:0045941-positive regulation of transcription390.642.71E-054.70E−02
 GO:0051254-positive regulation of RNA metabolic process350.582.82E-054.90E−02
 GO:0006350-transcription1031.702.98E-055.16E−02
 GO:0048598-embryonic morphogenesis260.433.39E-055.87E−02
 GO:0045449-regulation of transcription1212.005.08E-058.80E−02
 GO:0051173-positive regulation of nitrogen compound metabolic process400.662.05E-043.56E−01
 GO:0010604-positive regulation of macromolecule metabolic process480.794.63E-048.00E−01
 GO:0010828-positive regulation of glucose transport60.106.50E-041.12E+00
 GO:0010906-regulation of glucose metabolic process70.129.26E-041.59E+00
 GO:0035295-tube development180.301.08E-031.86E+00
 GO:0007268-synaptic transmission210.352.30E-033.91E+00
 GO:0007167-enzyme linked receptor protein signaling pathway230.382.36E-034.02E+00
 GO:0007169-transmembrane receptor protein tyrosine kinase signaling pathway170.283.32E-035.61E+00
KEGG pathways of upregulated cellular miRNA target genes
 hsa04115: p53 signaling pathway70.268.20E-049.03E−01
 hsa05200: pathways in cancer130.487.10E-037.58E+00
KEGG pathways of downregulated cellular miRNA target genes
 hsa05200: pathways in cancer420.969.06E-151.04E−11
 hsa04510: focal adhesion240.557.78E-088.91E−05
 hsa04012: ErbB signaling pathway130.301.84E-052.11E−02
 hsa04010: MAPK signaling pathway230.533.76E-054.30E−02
 hsa04722: neurotrophin signaling pathway150.343.77E-054.31E−02
 hsa04512: ECM-receptor interaction120.276.78E-057.76E−02
 hsa04115: p53 signaling pathway100.232.91E-043.33E−01
 hsa04530:tight junction140.323.38E-043.86E−01
 hsa04110: cell cycle130.306.31E-047.20E−01
 hsa04660: T cell receptor signaling pathway120.276.45E-047.36E−01
 hsa04062: chemokine signaling pathway160.379.24E-041.05E+00
 hsa04620: Toll-like receptor signaling pathway110.251.41E-031.61E+00
 hsa04662: B cell receptor signaling pathway90.212.65E-032.99E+00
 hsa04664: Fc ε RI signaling pathway90.213.39E-033.81E+00
 hsa04060: cytokine-cytokine receptor interaction180.414.29E-034.80E+00
 hsa04630:Jak-STAT signaling pathway120.271.10E-021.19E+01
 hsa04910: insulin signaling pathway110.251.14E-021.22E+01
 hsa04210: apoptosis80.182.15E-022.20E+01
 hsa04150: mTOR signaling pathway60.142.52E-022.53E+01
 hsa04370: VEGF signaling pathway70.163.31E-023.20E+01
KEGG pathways of the targets regulated by four or more EBV miRNAs
 hsa04350: TGF-β signaling pathway100.171.87E-032.14E+00
 hsa04920: adipocytokine signaling pathway80.135.64E-036.32E+00
 hsa00534: heparan sulfate biosynthesis50.089.15E-031.01E+01
 hsa04530: tight junction110.181.12E-021.22E+01
 hsa04910: insulin signaling pathway110.181.18E-021.28E+01
 hsa04520: adherens junction80.131.19E-021.29E+01
 hsa04310: Wnt signaling pathway110.182.42E-022.46E+01
 hsa04514: CAMs100.172.69E-022.70E+01
 hsa05200: pathways in cancer180.303.60E-023.45E+01

[i] CAM, cell adhesion molecules; EBV, Epstein-Barr virus; GO, Gene Ontology; NPC, nasopharyngeal carcinoma; KEGG, Kyoto Encyclopedia of Genes and Genomes; FDR, false discovery rate.

The targets of the upregulated miRNAs were significantly enriched in two pathway annotations (P<0.05): the p53 signaling and the cancer pathways. Among the targets of the downregulated miRNAs, the most statistically significant pathway annotations (P<0.05) were cancer pathways, focal adhesion, the ErbB signaling and the MAPK signaling pathways. Notably, target genes of both the upregulated and downregulated cellular miRNAs are involved in the p53 signaling pathway. The target genes of EBV miRNAs are involved in both the TGF-β and Wnt signaling pathways. This suggests that these pathways may play important roles in the development of NPC.

Construction of a cellular miRNA/gene regulatory network

An miRNA/gene regulatory network was constructed using miRNA-target gene pairs. In this network, all target genes were co-regulated by 4 or more miRNAs, enriched in signaling pathways and their expression was negatively correlated with cellular miRNA expression. Among the upregulated miRNAs, 15 miRNAs and 13 target genes were included in a canonical miRNA/gene regulatory network (Fig. 3A); for the down-regulated miRNAs, 27 miRNAs and 36 target genes were included (Fig. 3B).

Establishing an intercross regulatory network of cellular and EBV miRNAs in NPC

Analyses identified 25 target genes that were co-regulated by cellular and viral miRNAs and whose expression was negatively correlated with that of the cellular and viral miRNAs. As a result, we established a cellular miRNA/viral miRNA/gene regulatory network in NPC tissues (Fig. 3C). This network combined with the GO and pathway analyses suggests that cellular and EBV miRNAs may inhibit cell apoptosis and promote cell proliferation by jointly regulating the expression of ARHGEF12, SH3GLB1, DHCR24, IGFBP5, ITGA2, FOXO1 and FZD4.

Discussion

Although the biological effects of many miRNAs have been individually characterized, the effect of the dysregulation of multiple miRNAs on cellular functions and in tumors remains largely unknown. Furthermore, although each miRNA regulates hundreds of genes, their effect on an individual gene is moderate at best (45,46). Recent studies suggested that multiple miRNAs works in concert to regulate targets in a common pathway (47,48). Therefore, pathway analysis rather than individual target gene characterization, provides a better solution for evaluating the biological consequences of global miRNA dysregulation.

Extensive miRNA profiling has become an important tool for elucidating their roles in cancer as well as for identifying biomarkers. To link miRNA profiling data with biological functions a strategy was developed that analyzes the functions and pathways that are regulated collectively by co-expressed miRNAs using computationally predicted targets (49). However, since all target prediction algorithms generate a certain number of false positives, these results significantly reduce the reliability of the data when the targets predicted for multiple miRNAs are combined to analyze functional pathways. Therefore, the results of the pathway enrichment analysis according to co-regulated targets provide a good insight into the functional role of dysregulated miRNAs.

In the present study, we analyzed both the expression profiles of both cellular and EBV miRNAs in pooled laser capture microdissected NPC tissues using microarrays. The analysis of pooled NPC tissues avoids the large variation that occurs between samples to help identify bona fide miRNAs and target genes that are dysregulated in NPCs. In addition, few studies have performed a systematic analysis of the cellular and EBV miRNA expression profiles in 20 NPC tissues. Our analyses identified 50 differentially expressed cellular miRNAs and 9 differentially expressed EBV miRNAs, which regulate a total of 555 and 636 target genes, respectively. Furthermore, DAVID analysis revealed 37 and 12 pathways that were modulated by cellular and EBV miRNAs in NPC, respectively. Among these, the cellular miRNAs mainly regulate pathways involved in apoptosis and cell proliferation, whereas the EBV miRNAs mainly target genes in the TGF-β and Wnt signaling pathways. Furthermore, we constructed a cellular and EBV miRNA/target gene regulatory network (Fig. 3). This global miRNA and target gene analysis in NPC tissues is important to further understand the mechanism of NPC carcinogenesis and provide guidance for the diagnosis and treatment.

Previously, different groups identified a number of miRNAs that were dysregulated in NPC tissues using microarrays or RT-PCR, yet there was poor inter-platform consistency (11,32,33,36). Sengupta et al (11) analyzed 31 tumor and 10 normal surrounding tissues using microarrays, and identified 8 cellular miRNAs with ≥5-fold differential expression between tumor and normal tissues. Chen et al (32) quantified the expression of 270 human miRNAs in 13 NPCs and 9 adjacent normal tissues using RT-PCR, and identified 11 and 24 miRNAs that were significantly upregulated and down-regulated in NPC samples. Li et al (36) screened the miRNA expression changes in 8 NPCs and 4 normal nasopharyngeal tissues using microarrays, and identified 34 miRNAs that were expressed differentially in NPC tissues. Luo et al (33) performed miRNA microarrays using 20 different stage NPCs and 6 normal samples, and identified 11 miRNAs that may mediate the development of NPC. Among these four studies, only miR-34b and miR-34c were commonly identified as downregulated in NPC tissues. This poor consistency may be partially attributable to variations in NPC samples and the different platforms used for the measurements. Notably, miR-34b and miR-34c were also significantly downregulated in our pooled NPC tissues (Table III). miR-34 is a downstream target of p53, and promotes apoptosis (50). The downregulation of miR-34b and miR-34c strongly suggest that NPC exhibits enhanced cell survival due to inhibition of apoptosis pathways.

As previously mentioned, Sengupta et al (11) analyzed only miRNAs with ≥5-fold differential expression between tumor and normal tissues. They may have excluded some important miRNAs from their dataset. Therefore, we searched for other miRNAs that were commonly dysregulated in the remaining three datasets. This revealed that miR-18a was commonly overexpressed in NPC in these datasets as well as the present study. miR-18a belongs to the miR-17-92 cluster, and enhances tumorigenesis by inducing c-Myc expression (51). Similarly, miR-449a was commonly downregulated in NPC in the current dataset as well as by Chen et al (32) and Li et al (36). In addition, let-7b was downregulated significantly in NPC in the current dataset, which was consistent with the findings of Li et al (36). Since let-7 targets key cell cycle proto-oncogenes such as RAS, CDC25a, CDK6, and cyclin D to repress cell proliferation (52), we speculate that the downregulation of Let-7 in NPC may play a major role in tumor growth. Furthermore, miR-200a, miR-200b, and miR-31 were all downregulated in the present study and in the report by Chen et al (32). The miR-200 family inhibits NPC cell growth, migration and invasion (18). miR-31 is a tumor suppressor gene in breast cancer, and its reduced expression in NPC contributes to tumorigenesis and metastasis (53). Therefore, miR-18a, miR-449a, Let-7, the miR-200 family and miR-31 may be commonly dysregulated in NPC, and play roles in its tumorigenesis, proliferation and metastasis.

GO analysis found that the cellular miRNAs upregulated in NPC mainly regulate apoptosis-related genes, whereas the downregulated miRNAs primarily modulate genes related to cell proliferation. KEGG pathway analysis demonstrated that the upregulated miRNAs mainly target the p53 signaling pathway in cancer, whereas the downregulated miRNAs modulate pathways related to cancer, focal adhesion, ErbB and MAPK signaling (Table V). Chen et al (32) also observed that miRNAs modulated apoptosis and survival pathways and enhanced cell cycle progression and signaling pathways in NPC. Therefore, it is possible that the upregulated miRNAs such as miR-34b and -34c in NPC inhibit p53-mediated apoptosis to enhance cell survival; conversely, the downregulated miRNAs promote cell growth. This may be an important mechanism for the development of NPCs. Furthermore, we constructed a cellular miRNA/gene regulatory network, which provides vital information for the role of miRNAs in the tumorigenesis of NPCs.

EBV miRNAs are very abundant in NPC, and they may play important regulatory roles in tumor progression (31). Previously, Wong et al (39) pooled 5 NPC samples and identified 29 upregulated EBV-encoded miRNAs in NPC tissues. In the present study, the 9 EBV miRNAs that were upregulated in 20 NPC tissues (BART1-5p, BART3, BART4, BART5, BART6, BART7, BART8, BART10 and BART18-3p) were also identified by Wong et al (39). In addition, Zhu et al (38) reported that EBV-miR-BART4, BRAT6, and BRAT7 were expressed significantly in NPC samples using northern blotting. Some of these EBV miRNAs modulate the transformation, growth and survival of NPC. For example, miR-BART3 targets the DICE1 tumor suppressor to promote cellular growth and transformation in NPC (54). In addition, miR-BART5 regulates a cellular protein named p53-upregulated modulator of apoptosis (PUMA) to modulate host cell survival (25). miR-BART7 was proposed as a novel serological biomarker for the diagnosis of NPC (55). In the present study, KEGG pathway analysis showed that the upregulated EBV miRNAs mainly target genes in the TGF-β and Wnt signaling pathways. Both the TGF-β and Wnt signaling pathways were previously shown to be dysregulated in NPC (32,56). In addition, the TGF-β signaling pathway plays a critical role in the progression of human cancer (57). Activation of the Wnt signaling pathway in NPC is also associated with NPC development (58). Therefore, the dysregulation of EBV miRNAs such as BART3 and BART5 in NPC plays key roles in the transformation and growth of NPC.

To globally understand the genesis of NPC, we also established a cellular miRNA/viral miRNA/gene regulatory network in NPC tissues (Fig. 3C). This network may help explain how the EBV works with cellular miRNAs to regulate gene expression and promote the transformation of NPC.

In conclusion, we performed a systematic characterization of the cellular and EBV miRNA expression profiles in NPC tissues. Our data suggest that cellular miRNAs such as miRNA-34b, miRNA-34c, miR-18a and Let-7 mainly target genes involved in apoptosis and cell proliferation. In addition, we found that the upregulated EBV miRNAs in NPCs mainly target the TGF-β and Wnt signaling pathways. These identified cellular and EBV miRNAs may be potential therapeutic targets for the treatment of NPC. In addition, we identified 25 target genes that were co-regulated by cellular and viral miRNAs. Further analysis of this cellular miRNA/viral miRNA/gene regulatory network may unveil the mechanism of NPC carcinogenesis.

Acknowledgments

The present study was supported by the National Natural Science Foundation of China (nos. 81230053, 81172559 and 81272959), the Hunan Provincial Natural Science Foundation of China (no. 10JJ3045), and the Scientific Research Fund of Hunan Provincial Education Department (no. 08C566).

Abbreviations:

NPC

nasopharyngeal carcinoma

NNET

normal nasopharyngeal epithelial tissue

EBV

Epstein-Barr virus

GO

Gene Ontology

KEGG

Kyoto Encyclopedia of Genes and Genomes

LCM

laser capture microdissected

NCBI-GEO

National Center for Biotechnology Information-Gene Expression Omnibus

qRT-PCR

quantitative reverse transcription-PCR

References

1 

Ambros V: The functions of animal microRNAs. Nature. 431:350–355. 2004. View Article : Google Scholar : PubMed/NCBI

2 

Croce CM and Calin GA: miRNAs, cancer, and stem cell division. Cell. 122:6–7. 2005. View Article : Google Scholar : PubMed/NCBI

3 

Bushati N and Cohen SM: microRNA functions. Annu Rev Cell Dev Biol. 23:175–205. 2007. View Article : Google Scholar : PubMed/NCBI

4 

Chen CZ: MicroRNAs as oncogenes and tumor suppressors. N Engl J Med. 353:1768–1771. 2005. View Article : Google Scholar : PubMed/NCBI

5 

Calin GA, Sevignani C, Dumitru CD, Hyslop T, Noch E, Yendamuri S, Shimizu M, Rattan S, Bullrich F, Negrini M, et al: Human microRNA genes are frequently located at fragile sites and genomic regions involved in cancers. Proc Natl Acad Sci USA. 101:2999–3004. 2004. View Article : Google Scholar : PubMed/NCBI

6 

Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, et al: MicroRNA expression profiles classify human cancers. Nature. 435:834–838. 2005. View Article : Google Scholar : PubMed/NCBI

7 

Murakami Y, Yasuda T, Saigo K, Urashima T, Toyoda H, Okanoue T and Shimotohno K: Comprehensive analysis of microRNA expression patterns in hepatocellular carcinoma and non-tumorous tissues. Oncogene. 25:2537–2545. 2006. View Article : Google Scholar

8 

Kida Y and Han YP: MicroRNA expression in colon adenocarcinoma. JAMA. 299:2628–2629. 2008. View Article : Google Scholar : PubMed/NCBI

9 

Yanaihara N, Caplen N, Bowman E, Seike M, Kumamoto K, Yi M, Stephens RM, Okamoto A, Yokota J, Tanaka T, et al: Unique microRNA molecular profiles in lung cancer diagnosis and prognosis. Cancer Cell. 9:189–198. 2006. View Article : Google Scholar : PubMed/NCBI

10 

Schetter AJ, Leung SY, Sohn JJ, Zanetti KA, Bowman ED, Yanaihara N, Yuen ST, Chan TL, Kwong DLW, Au GKH, et al: MicroRNA expression profiles associated with prognosis and therapeutic outcome in colon adenocarcinoma. JAMA. 299:425–436. 2008. View Article : Google Scholar : PubMed/NCBI

11 

Sengupta S, den Boon JA, Chen IH, Newton MA, Stanhope SA, Cheng YJ, Chen CJ, Hildesheim A, Sugden B and Ahlquist P: MicroRNA 29c is down-regulated in nasopharyngeal carcinomas, up-regulating mRNAs encoding extracellular matrix proteins. Proc Natl Acad Sci USA. 105:5874–5878. 2008. View Article : Google Scholar : PubMed/NCBI

12 

Alajez NM, Lenarduzzi M, Ito E, Hui AB, Shi W, Bruce J, Yue S, Huang SH, Xu W, Waldron J, et al: miR-218 suppresses nasopharyngeal cancer progression through downregulation of survivin and the SLIT2-ROBO1 pathway. Cancer Res. 71:2381–2391. 2011. View Article : Google Scholar : PubMed/NCBI

13 

Lu J, He ML, Wang L, Chen Y, Liu X, Dong Q, Chen YC, Peng Y, Yao KT, Kung HF, et al: MiR-26a inhibits cell growth and tumorigenesis of nasopharyngeal carcinoma through repression of EZH2. Cancer Res. 71:225–233. 2011. View Article : Google Scholar : PubMed/NCBI

14 

Yi C, Wang Q, Wang L, Huang Y, Li L, Liu L, Zhou X, Xie G, Kang T, Wang H, et al: miR-663, a microRNA targeting p21WAF1/CIP1, promotes the proliferation and tumorigenesis of nasopharyngeal carcinoma. Oncogene. 31:4421–4433. 2012. View Article : Google Scholar : PubMed/NCBI

15 

Wong TS, Man OY, Tsang CM, Tsao SW, Tsang RK, Chan JY, Ho WK, Wei WI and To VS: MicroRNA let-7 suppresses naso-pharyngeal carcinoma cells proliferation through downregulating c-Myc expression. J Cancer Res Clin Oncol. 137:415–422. 2011. View Article : Google Scholar :

16 

Zhang L, Deng T, Li X, Liu H, Zhou H, Ma J, Wu M, Zhou M, Shen S, Li X, et al: microRNA-141 is involved in a nasopharyngeal carcinoma-related genes network. Carcinogenesis. 31:559–566. 2010. View Article : Google Scholar : PubMed/NCBI

17 

Deng M, Tang H, Zhou Y, Zhou M, Xiong W, Zheng Y, Ye Q, Zeng X, Liao Q, Guo X, et al: miR-216b suppresses tumor growth and invasion by targeting KRAS in nasopharyngeal carcinoma. J Cell Sci. 124:2997–3005. 2011. View Article : Google Scholar : PubMed/NCBI

18 

Xia H, Ng SS, Jiang S, Cheung WK, Sze J, Bian XW, Kung HF and Lin MC: miR-200a-mediated downregulation of ZEB2 and CTNNB1 differentially inhibits nasopharyngeal carcinoma cell growth, migration and invasion. Biochem Biophys Res Commun. 391:535–541. 2010. View Article : Google Scholar

19 

Busson P, Keryer C, Ooka T and Corbex M: EBV-associated nasopharyngeal carcinomas: From epidemiology to virus-targeting strategies. Trends Microbiol. 12:356–360. 2004. View Article : Google Scholar : PubMed/NCBI

20 

Wei WI and Sham JS: Nasopharyngeal carcinoma. Lancet. 365:2041–2054. 2005. View Article : Google Scholar : PubMed/NCBI

21 

Qiu J, Cosmopoulos K, Pegtel M, Hopmans E, Murray P, Middeldorp J, Shapiro M and Thorley-Lawson DA: A novel persistence associated EBV miRNA expression profile is disrupted in neoplasia. PLoS Pathog. 7:e10021932011. View Article : Google Scholar : PubMed/NCBI

22 

Barth S, Meister G and Grässer FA: EBV-encoded miRNAs. Biochim Biophys Acta. 1809:631–640. 2011. View Article : Google Scholar : PubMed/NCBI

23 

Xia T, O'Hara A, Araujo I, Barreto J, Carvalho E, Sapucaia JB, Ramos JC, Luz E, Pedroso C, Manrique M, et al: EBV microRNAs in primary lymphomas and targeting of CXCL-11 by ebv-mir-BHRF1-3. Cancer Res. 68:1436–1442. 2008. View Article : Google Scholar : PubMed/NCBI

24 

Nachmani D, Stern-Ginossar N, Sarid R and Mandelboim O: Diverse herpesvirus microRNAs target the stress-induced immune ligand MICB to escape recognition by natural killer cells. Cell Host Microbe. 5:376–385. 2009. View Article : Google Scholar : PubMed/NCBI

25 

Choy EY, Siu KL, Kok KH, Lung RW, Tsang CM, To KF, Kwong DL, Tsao SW and Jin DY: An Epstein-Barr virus-encoded microRNA targets PUMA to promote host cell survival. J Exp Med. 205:2551–2560. 2008. View Article : Google Scholar : PubMed/NCBI

26 

Li Z, Chen X, Li L, Liu S, Yang L, Ma X, Tang M, Bode AM, Dong Z, Sun L, et al: EBV encoded miR-BHRF1-1 potentiates viral lytic replication by downregulating host p53 in nasopharyngeal carcinoma. Int J Biochem Cell Biol. 44:275–279. 2012. View Article : Google Scholar

27 

Lo AK, To KF, Lo KW, Lung RW, Hui JW, Liao G and Hayward SD: Modulation of LMP1 protein expression by EBV-encoded microRNAs. Proc Natl Acad Sci USA. 104:16164–16169. 2007. View Article : Google Scholar : PubMed/NCBI

28 

Barth S, Pfuhl T, Mamiani A, Ehses C, Roemer K, Kremmer E, Jäker C, Höck J, Meister G and Grässer FA: Epstein-Barr virus-encoded microRNA miR-BART2 down-regulates the viral DNA polymerase BALF5. Nucleic Acids Res. 36:666–675. 2008. View Article : Google Scholar :

29 

Lung RWM, Tong JHM, Sung YM, Leung PS, Ng DCH, Chau SL, Chan AWH, Ng EKO, Lo KW and To KF: Modulation of LMP2A expression by a newly identified Epstein-Barr virus-encoded microRNA miR-BART22. Neoplasia. 11:1174–1184. 2009. View Article : Google Scholar : PubMed/NCBI

30 

Lo AK, Dawson CW, Jin DY and Lo KW: The pathological roles of BART miRNAs in nasopharyngeal carcinoma. J Pathol. 227:392–403. 2012. View Article : Google Scholar : PubMed/NCBI

31 

Marquitz AR and Raab-Traub N: The role of miRNAs and EBV BARTs in NPC. Semin Cancer Biol. 22:166–172. 2012. View Article : Google Scholar :

32 

Chen HC, Chen GH, Chen YH, Liao WL, Liu CY, Chang KP, Chang YS and Chen SJ: MicroRNA deregulation and pathway alterations in nasopharyngeal carcinoma. Br J Cancer. 100:1002–1011. 2009. View Article : Google Scholar : PubMed/NCBI

33 

Luo Z, Zhang L, Li Z, Li X, Li G, Yu H, Jiang C, Dai Y, Guo X, Xiang J, et al: An in silico analysis of dynamic changes in microRNA expression profiles in stepwise development of nasopharyngeal carcinoma. BMC Med Genomics. 5:32012. View Article : Google Scholar : PubMed/NCBI

34 

Liu N, Chen NY, Cui RX, Li WF, Li Y, Wei RR, Zhang MY, Sun Y, Huang BJ, Chen M, et al: Prognostic value of a microRNA signature in nasopharyngeal carcinoma: A microRNA expression analysis. Lancet Oncol. 13:633–641. 2012. View Article : Google Scholar : PubMed/NCBI

35 

Zeng X, Xiang J, Wu M, Xiong W, Tang H, Deng M, Li X, Liao Q, Su B, Luo Z, et al: Circulating miR-17, miR-20a, miR-29c, and miR-223 combined as non-invasive biomarkers in nasopharyngeal carcinoma. PLoS One. 7:e463672012. View Article : Google Scholar : PubMed/NCBI

36 

Li T, Chen JX, Fu XP, Yang S, Zhang Z, Chen KhH and Li Y: microRNA expression profiling of nasopharyngeal carcinoma. Oncol Rep. 25:1353–1363. 2011.PubMed/NCBI

37 

Chen SJ, Chen GH, Chen YH, Liu CY, Chang KP, Chang YS and Chen HC: Characterization of Epstein-Barr virus miRNAome in nasopharyngeal carcinoma by deep sequencing. PLoS One. 5:e127452010. View Article : Google Scholar : PubMed/NCBI

38 

Zhu JY, Pfuhl T, Motsch N, Barth S, Nicholls J, Grässer F and Meister G: Identification of novel Epstein-Barr virus microRNA genes from nasopharyngeal carcinomas. J Virol. 83:3333–3341. 2009. View Article : Google Scholar : PubMed/NCBI

39 

Wong AM, Kong KL, Tsang JW, Kwong DL and Guan XY: Profiling of Epstein-Barr virus-encoded microRNAs in nasopharyngeal carcinoma reveals potential biomarkers and oncomirs. Cancer. 118:698–710. 2012. View Article : Google Scholar

40 

Cheng AL, Huang WG, Chen ZC, Peng F, Zhang PF, Li MY, Li F, Li JL, Li C, Yi H, et al: Identification of novel biomarkers for differentiation and prognosis of nasopharyngeal carcinoma by laser capture microdissection and proteomic analysis. Clin Cancer Res. 14:435–445. 2008. View Article : Google Scholar : PubMed/NCBI

41 

Elefant N, Berger A, Shein H, Hofree M, Margalit H and Altuvia Y: RepTar: A database of predicted cellular targets of host and viral miRNAs. Nucleic Acids Res. 39:D188–D194. 2011. View Article : Google Scholar :

42 

Maragkakis M, Reczko M, Simossis VA, Alexiou P, Papadopoulos GL, Dalamagas T, Giannopoulos G, Goumas G, Koukis E, Kourtis K, et al: DIANA-microT web server: Elucidating microRNA functions through target prediction. Nucleic Acids Res. 37:W273–W276. 2009. View Article : Google Scholar : PubMed/NCBI

43 

Yang JH, Li JH, Shao P, Zhou H, Chen YQ and Qu LH: starBase: A database for exploring microRNA-mRNA interaction maps from Argonaute CLIP-Seq and Degradome-Seq data. Nucleic Acids Res. 39:D202–D209. 2011. View Article : Google Scholar

44 

Kim H, Park S, Min H and Yoon S: vHoT: A database for predicting interspecies interactions between viral microRNA and host genomes. Arch Virol. 157:497–501. 2012. View Article : Google Scholar

45 

Baek D, Villén J, Shin C, Camargo FD, Gygi SP and Bartel DP: The impact of microRNAs on protein output. Nature. 455:64–71. 2008. View Article : Google Scholar : PubMed/NCBI

46 

Selbach M, Schwanhäusser B, Thierfelder N, Fang Z, Khanin R and Rajewsky N: Widespread changes in protein synthesis induced by microRNAs. Nature. 455:58–63. 2008. View Article : Google Scholar : PubMed/NCBI

47 

Cloonan N, Brown MK, Steptoe AL, Wani S, Chan WL, Forrest AR, Kolle G, Gabrielli B and Grimmond SM: The miR-17-5p microRNA is a key regulator of the G1/S phase cell cycle transition. Genome Biol. 9:R1272008. View Article : Google Scholar : PubMed/NCBI

48 

Liu Q, Fu H, Sun F, Zhang H, Tie Y, Zhu J, Xing R, Sun Z and Zheng X: miR-16 family induces cell cycle arrest by regulating multiple cell cycle genes. Nucleic Acids Res. 36:5391–5404. 2008. View Article : Google Scholar : PubMed/NCBI

49 

Gusev Y, Schmittgen TD, Lerner M, Postier R and Brackett D: Computational analysis of biological functions and pathways collectively targeted by co-expressed microRNAs in cancer. BMC Bioinformatics. 8(Suppl 7): S162007. View Article : Google Scholar : PubMed/NCBI

50 

Chang TC, Wentzel EA, Kent OA, Ramachandran K, Mullendore M, Lee KH, Feldmann G, Yamakuchi M, Ferlito M, Lowenstein CJ, et al: Transactivation of miR-34a by p53 broadly influences gene expression and promotes apoptosis. Mol Cell. 26:745–752. 2007. View Article : Google Scholar : PubMed/NCBI

51 

Dews M, Homayouni A, Yu D, Murphy D, Sevignani C, Wentzel E, Furth EE, Lee WM, Enders GH, Mendell JT, et al: Augmentation of tumor angiogenesis by a Myc-activated microRNA cluster. Nat Genet. 38:1060–1065. 2006. View Article : Google Scholar : PubMed/NCBI

52 

Johnson CD, Esquela-Kerscher A, Stefani G, Byrom M, Kelnar K, Ovcharenko D, Wilson M, Wang X, Shelton J, Shingara J, et al: The let-7 microRNA represses cell proliferation pathways in human cells. Cancer Res. 67:7713–7722. 2007. View Article : Google Scholar : PubMed/NCBI

53 

O'Day E and Lal A: MicroRNAs and their target gene networks in breast cancer. Breast Cancer Res. 12:2012010. View Article : Google Scholar : PubMed/NCBI

54 

Lei T, Yuen KS, Xu R, Tsao SW, Chen H, Li M, Kok KH and Jin DY: Targeting of DICE1 tumor suppressor by Epstein-Barr virus-encoded miR-BART3* microRNA in nasopharyngeal carcinoma. Int J Cancer. 133:79–87. 2013. View Article : Google Scholar : PubMed/NCBI

55 

Zhang G, Zong J, Lin S, Verhoeven RJ, Tong S, Chen Y, Ji M, Cheng W, Tsao SW, Lung M, et al: Circulating Epstein-Barr virus microRNAs miR-BART7 and miR-BART13 as biomarkers for nasopharyngeal carcinoma diagnosis and treatment. Int J Cancer. 136:E301–E312. 2015. View Article : Google Scholar

56 

Zeng ZY, Zhou YH, Zhang WL, Xiong W, Fan SQ, Li XL, Luo XM, Wu MH, Yang YX, Huang C, et al: Gene expression profiling of nasopharyngeal carcinoma reveals the abnormally regulated Wnt signaling pathway. Hum Pathol. 38:120–133. 2007. View Article : Google Scholar

57 

Meulmeester E and Ten Dijke P: The dynamic roles of TGF-β in cancer. J Pathol. 223:205–218. 2011. View Article : Google Scholar

58 

Chou J, Lin YC, Kim J, You L, Xu Z, He B and Jablons DM: Nasopharyngeal carcinoma - review of the molecular mechanisms of tumorigenesis. Head Neck. 30:946–963. 2008. View Article : Google Scholar : PubMed/NCBI

Related Articles

Journal Cover

November-2015
Volume 34 Issue 5

Print ISSN: 1021-335X
Online ISSN:1791-2431

Sign up for eToc alerts

Recommend to Library

Copy and paste a formatted citation
x
Spandidos Publications style
Wan X, Yi H, Qu J, He Q and Xiao Z: Integrated analysis of the differential cellular and EBV miRNA expression profiles in microdissected nasopharyngeal carcinoma and non-cancerous nasopharyngeal tissues. Oncol Rep 34: 2585-2601, 2015.
APA
Wan, X., Yi, H., Qu, J., He, Q., & Xiao, Z. (2015). Integrated analysis of the differential cellular and EBV miRNA expression profiles in microdissected nasopharyngeal carcinoma and non-cancerous nasopharyngeal tissues. Oncology Reports, 34, 2585-2601. https://doi.org/10.3892/or.2015.4237
MLA
Wan, X., Yi, H., Qu, J., He, Q., Xiao, Z."Integrated analysis of the differential cellular and EBV miRNA expression profiles in microdissected nasopharyngeal carcinoma and non-cancerous nasopharyngeal tissues". Oncology Reports 34.5 (2015): 2585-2601.
Chicago
Wan, X., Yi, H., Qu, J., He, Q., Xiao, Z."Integrated analysis of the differential cellular and EBV miRNA expression profiles in microdissected nasopharyngeal carcinoma and non-cancerous nasopharyngeal tissues". Oncology Reports 34, no. 5 (2015): 2585-2601. https://doi.org/10.3892/or.2015.4237