Microarray analysis of miRNA based on the regional lymph node metastasis status of esophageal squamous cell carcinoma.

BACKGROUND: Tendency toward extensive regional lymph node metastasis (LNM) is an important clinical characteristic of esophageal squamous cell carcinoma (ESCC) and differs greatly between patients. MicroRNAs (miRNAs) are involved in the invasion and metastasis of ESCC. We performed a microarray analysis of miRNAs based on LNM status to identify tumor-specific miRNAs for the prediction of LNM in ESCC.
METHODS: Four pairs of ESCC tumor tissues with or without LNM were selected for microarray analysis to identify differentially expressed miRNAs, then 50 tumor tissue samples were selected to verify the differences of the screened miRNAs with quantitative reverse transcription-polymerase chain reaction (qRT-PCR). The relations between candidate miRNAs and clinicopathologic features were analyzed to confirm tumor specificity in the prediction of LNM in ESCC. Target gene prediction using miRWalk2.0 was used to analyze the potential mechanisms of the tumor-specific miRNAs.
RESULTS: The present microarray analysis identified significant differential expression in 62 miRNAs in the ESCC samples with LNM, compared to those without LNM. Of these, 19 miRNAs were selected for qRT-PCR verification, and three miRNAs were significantly upregulated in ESCC samples with LNM compared to those without LNM. The three miRNAs were not significantly associated with any other clinicopathologic features except for the TNM stage and could be regarded as tumor-specific miRNAs capable of predicting LNM in ESCC. Finally, 858 mRNAs were significantly associated with tumor-specific miRNAs and possibly involved in the regulation of LNM in ESCC.
CONCLUSIONS: The present microarray analysis based on LNM status identified three tumor-specific miRNAs for predicting regional LNM in ESCC, which provides valuable clues for potential mechanism research and guarantees further investigation. 2021 Translational Cancer Research. All rights reserved.

Introduction

Esophageal cancer (EC) is one of the most common malignant tumors in the world, ranking eighth in cancer-related morbidity and sixth in cancer-related mortality worldwide (1). Esophageal squamous cell carcinoma (ESCC) is the dominant histological subtype in EC, accounting for almost 90% of EC patients, especially in China (2). For decades, although researchers have endeavored to find breakthroughs in the treatment of ESCC, the outcome has been disappointing and the survival rate remains dismal (3,4). One of the most important reasons for the failure of ESCC treatment is regional lymph node metastasis (LNM); the regional LNM rate in ESCC is reported to reach 40–70% (5,6). In ESCC, the lymphatic drainage regions include the lower neck, supraclavicular, mediastinum, and even the upper abdomen. It is very difficult to control such a wide range of lymphatic drainage regions for surgery and radiotherapy, which are the two main treatments of ESCC, because the control of such a wide range of lymphatic drainage regions will greatly increase the difficulty and complication of treatment, even might lead to the death of patients (7-9). Therefore, the control of a selective range of lymphatic drainage regions will be conducive to reduce the difficulty and complication of treatment, and will be worthy of be studied.

In fact, although in clinical practice there is a relatively high and wide incidence of LNM in ESCC, there is a great difference among patients (5,6). In current clinical practice, in order to distinguish this difference, lymph node dissection or irradiation is usually carried out selectively according to the location of ESCC (10). However, due to the characteristics of lymph node jumping in the metastasis of ESCC (11), these disposals often lead to the omission and increase the risk of LNM after treatment. With the progression of modern molecular biology, researchers are gradually realizing that intrinsic molecular biological characteristics may lead to the great difference in LNM among different patients (12). Therefore, the search for an appropriate molecular biomarker to predict regional LNM to guide the selective lymphatic dissection or irradiation of ESCC will undoubtedly reduce the difficulty and complication of treatment, decrease the incidence of treatment failure for LNM, increase survival rate of ESCC, and bring great progress in the treatment of ESCC.

MicroRNAs (miRNAs) are a class of non-coding small RNAs, which can modulate the target gene expression at the post-transcriptional level. MiRNAs participate in a variety of physiological and pathological processes, including cell proliferation, cell death, stress response, apoptosis, fat metabolism, cell differentiation, and so on (13). Accumulating evidence has shown that miRNAs are closely associated with malignant tumors, and many miRNAs act as oncogenes or tumor suppressors in several malignant tumors (14). In recent years, some studies have found that some miRNAs participate in the regulation of cancer invasion and metastasis, and can be used for the prediction of cancer invasion and metastasis (15,16). In ESCC, previous studies have found that miRNA-145, miRNA-200b, and miRNA-208 are involved in regulating the invasion and metastasis of ESCC, and miRNA-17a, miRNA-19b, miRNA-302b, miRNA-100, miRNA-199a, and miRNA-92a were associated with LNM in ESCC (17,18). However, these conclusions are mostly based on experiments in vitro or with a small sample of clinical cases. Due to the lack of characteristic studies based on the regional LNM status of ESCC, there are no miRNAs that can be used for the prediction of regional LNM in ESCC at present.

Therefore, in this study, mature miRNA microarray analysis based on regional LNM status was used to screen for miRNAs with significant differential expression in ESCC. Simultaneously, the verification of a large sample of ESCC tissues with or without regional LNM was performed to obtain tumor-specific miRNAs for the prediction of regional LNM in ESCC to guide the clinical treatment of ESCC. We present the following article in accordance with the MDAR checklist (available at http://dx.doi.org/10.21037/tcr-20-1969).

Methods

Samples

Primary tumor tissues were taken from EC patients and immediately frozen in liquid nitrogen after surgical resection and prepared for the experiments. All the samples selected were in accord with the following criteria: patients underwent radical EC resection and regional lymphatic dissection; samples were pathologically confirmed as ESCC postoperatively and regional LNM status was pathologically evaluated postoperatively; none of preoperative chemotherapy or radiotherapy was taken; all the patients did not accompanied with any other serious diseases and the expected survive was more than 6 months. Eligible patients were registered from the Shandong Cancer Hospital and the Second Hospital of Shandong University from September 2017 to March 2018. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). All the patients signed the informed consent, and the study was approved by the Ethics Committees of the two hospitals [registration number: KYLL-2017(GJ)P-0036].

RNA isolation and microarray analysis

According to the manufacturer’s protocol, total RNA was respectively isolated from four pairs of ESCC tissues with and without LNM using a TRIzol reagent (Invitrogen, Carlsbad, CA, USA), and then was purified using a mirVana miRNA Isolation Kit (Ambion, Austin, TX, USA). To ensure that the total RNA isolated met the experimental requirements, a NanoDrop ND-1000 Spectrophotometer (Thermo Fisher Scientific) was used to quantify the purity and concentration of total RNA, and an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) was used to evaluate the integrity of total RNA.

After verifying that the quality of total RNA extraction from four pairs of ESCC tissues reached the experimental requirements, miRNA expression profiles of ESCC tissue with and without LNM were analyzed with the human miRNA array version 21.0 (Agilent). During the miRNA array analysis, a miRNA Complete Labeling and Hyb Kit (Agilent) was firstly used to label miRNAs according to the manufacturer’s protocol. Then the Cy3-labeled samples were hybridized according to the guidelines for the Agilent miRNA array. Through the whole miRNA array analysis, the labeling spike-in RNA and the Hyb spike-in RNA from Agilent were used for the quality control of the experimental process. After the hybridization, the slides were washed and dried according to the manufacturer’s manuals, and then were scanned using the Agilent microarray scanner.

Following the scan, the images were imported into the Agilent Feature Extraction software (v10.7) (Agilent Instruments) for further grid alignment and data extraction. After that, Agilent GeneSpring software was used for the normalization of the raw signal data extracted. And then, the differentially expressed miRNAs between two groups were identified using a standard fold change (abs) >2 and P<0.05. The entire miRNAs datasets mentioned here are available from the Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo/)

miRNA target gene prediction analysis

After microarray analysis, miRNA target gene prediction analysis was performed for the differentially expressed miRNAs. MiRNA target gene prediction analysis using miRWalk2.0 (http://zmf.umm.uni-heidelberg.de/mirwalk2/) generated the prediction results of miRNA target genes from 12 existing miRNA-target prediction programs, including DIANA-microTv4.0, DIANA-microT-CDS, miRanda-rel2010, mirBridge, miRDB4.0, miRmap, miRNAMap, PicTar2, PITA, RNA22v2, RNAhybrid2.1, and Targetscan6.2. A gene would be regarded as a target of a differentially expressed miRNA if it was predicted in at least seven of the 12 miRNA-target prediction programs.

Gene ontology enrichment and pathway analysis

After miRNA target gene prediction analysis, gene ontology (GO) enrichment and pathway analyses of the target mRNAs for the differentially expressed miRNAs were performed. The GO project provides the annotations for the genes and gene products of any organism (http://www.geneontology.org), which covers the roles of mRNAs in the three domains including biological process, cellular component, and molecular function. The Fisher’s exact method was used to test whether the overlap between the differentially expressed list and the GO annotated list was greater than that expected by chance. While pathway analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis was used to explore the function of target genes. The Fisher’s exact test was used to identify a significant pathway, and the P value was used to defined the threshold of significance; the lower P value represented the more significant pathway.

MiRNA screen and verification

Based on the results of the microarray analysis, the differentially expressed miRNAs were screened for verification according to the criteria as follow: (I) fold change >2.0 and P<0.05, and (II) the raw signal intensity of each miRNA >1 and with a “detected” Feature Extraction Flag. After the miRNA screen, 50 ESCC samples, including 25 without LNM and 25 with LNM were selected for miRNA verification using quantitative reverse transcription-polymerase chain reaction (qRT-PCR). The significant miRNAs were regarded as candidate miRNAs capable of predicting regional LNM in ESCC.

Correlation analysis between candidate miRNAs and clinicopathologic features

To verify the specificity of the candidate miRNAs capable of predicting regional LNM in ESCC, the relations between verified miRNAs and clinicopathologic features were also analyzed. The clinicopathologic features analyzed included sex, age, tumor size, tumor location, T stage, TNM stage, tumor clinical subtype and tumor differentiation degree. Finally, tumor-specific miRNAs capable of predicting regional LNM in ESCC were obtained. We defined the postoperative stage according to the 2017 version TNM staging classification system of the American Joint Committee on Cancer (AJCC).

Potential mechanism analysis of the specific miRNAs

According to the results of the miRNA target gene prediction analysis, potential target genes of the tumor-specific miRNAs and their enriched pathways were also analyzed, and the potential mechanisms of the tumor-specific miRNAs regulating regional LNM in ESCC were predicted.

Statistical analysis

Statistical analysis was carried out using SPSS 20.0 software (IBM SPSS Statistics for Windows, Armonk, NY, USA) except the analysis methods indicated above. The expression levels of miRNAs were displayed as the mean ± standard deviation (SD). The expression differences of miRNAs between the two groups were compared with the independent-samples t-test. The differences in clinicopathologic features between the two groups were compared using the chi-square test. The correlations between specific miRNAs and tumor size were tested with Pearson correlation analysis. The value of P<0.05 was considered statistically significant.

Results

Sample selection

According to the inclusion criteria of samples, except LNM status, the clinicopathologic features of the four pairs of ESCC tissues selected for microarray analysis were matched, and those of the 50 tumor samples including 25 without LNM and 25 with LNM used to verify the effectiveness of the screened miRNAs were similar between the two groups (all P values >0.05). The clinicopathologic features of the samples selected in this study are presented in .

Table 1

The clinical information of samples for microarray analysis

Sample	Age (years)	Gender	Primary location	Pathology grade	Clinical subtype	Tumor size (cm2)	Stages†
NLNM‡1	59	Male	Lower thoracic esophagus	Moderate differentiation	Ulcerative type	5.8×2.5	pT3N0M0
NLNM2	62	Male	Midthoracic esophagus	Moderate differentiation	Ulcerative type	6.0×2.3	pT2N0M0
NLNM3	68	Male	Lower thoracic esophagus	Poor differentiation	Ulcerative type	5.0×1.6	pT3N0M0
NLNM4	67	Female	Lower thoracic esophagus	Moderate differentiation	Ulcerative type	6.5×3.0	pT3N0M0
LNM§1	58	Male	Lower thoracic esophagus	Moderate differentiation	Ulcerative type	5.0×3.0	pT3N3M0
LNM2	73	Male	Midthoracic esophagus	Moderate differentiation	Ulcerative type	4.0×3.5	pT2N1M0
LNM3	75	Male	Lower thoracic esophagus	Poor differentiation	Ulcerative type	3.0×2.5	pT3N2M0
LNM4	56	Female	Lower thoracic esophagus	Moderate differentiation	Ulcerative type	5.0×4.0	pT3N1M0

†, postoperative pathological stage; ‡, NLNM: represent sample without lymph node metastasis; §, LNM: represent sample with lymph node metastasis.

Table 2

The clinicopathological characteristics for verified samples

Main features	Group 1†	Group 2‡	P value
Gender (n)			1.000
Male	20	20
Female	5	5
Median age (year)	62 (46 to 76)	63 (45 to 82)
Age (n)			0.564
≥60 years	14	16
<60 years	11	9
Primary location (n)			0.258
Midthoracic esophagus	14	10
Lower thoracic esophagus	11	15
Pathology grade (n)			0.248
Poor	8	12
Moderate and well	17	13
Clinical subtype (n)			0.395
Ulcerative type	10	13
Medullary and other type	15	12
Tumor size (mean ± SD) (cm2)	8.14±5.55	9.21±6.59	0.539
T stage (n)			0.480
T1–2	6	4
T3	19	21
N stage (n)			–
N0	25	0
N1	0	18
N2	0	6
N3	0	1
TNM stage§ (n)			–
IB	1	0
IIA	14	0
IIB	10	1
IIIA	0	3
IIIB	0	20
IVA	0	1

†, Group 1: expression levels of miRNAs in cancer tissues without lymph node metastasis; ‡, Group 2: expression levels of miRNAs in cancer tissues with lymph node metastasis; §, postoperative pathological stage.

Overview of the microarray analysis

To guarantee the reliability of the microarray analysis results, strict quality control was performed in the sample selection, storage to RNA extraction, and microarray analysis. The final results of the microarray analysis reached our expectations.

According to the results of the comparative microarray analysis, 62 miRNAs showed obvious differential expression (fold change ≥2.0, P<0.05) in ESCC samples with LNM compared to those without LNM. Among them, 19 miRNAs were upregulated and 43 miRNAs were downregulated in ESCC samples with LNM compared to those without LNM. The results are expressed as a cluster plot, scatter plot, volcano plot, and Circos plot in .

Figure 1

Plots for the microarray analysis of the samples. (A) Hierarchical cluster analysis plots for the differentially expressed miRNAs in ESCC tissues with and without LNM show good similarities among cancer tissues of different groups. Yellow indicates relatively high expression and blue indicates relatively low expression of miRNAs. Each RNA is represented by a single row of colored boxes and each sample by a single column. (B) The Circos plot for the locations of differentially expressed miRNAs in ESCC tissues with and without LNM. Red represents significantly upregulated miRNAs, green represents significantly downregulated miRNAs, the length of the pillar represents a multiple of the difference. (C,D) The scatter plots and the volcano plots illustrate the distribution of miRNAs. Red dots represent miRNAs that are significantly upregulated, green dots represent miRNAs that are significantly downregulated, and black dots represent miRNAs without significant differences. LNM, lymph node metastasis; ESCC, esophageal squamous cell carcinoma.

MiRNA target gene prediction analysis

Further miRNA target gene prediction analysis identified 50 miRNAs, including 19 upregulated miRNAs and 31 downregulated miRNAs, among the 62 differentially expressed miRNAs with related target genes (target mRNA) according to the 12 existing miRNA-target prediction programs based on miRWalk2.0. The number of predicted target genes reached 8,637 and the potentially related pairs of miRNA-mRNA reached 37,231 according to the miRNA target gene prediction criteria.

GO and pathway analyses

To understand the functions of the underlying target mRNAs, GO analysis and pathway analysis were performed. The results of GO analysis showed that the majority of underlying target mRNAs was involved in biological process, one fraction of them was defined as cellular component, and a few of them participated in molecular function (these are shown in ). Among the roles of mRNAs, cellular processes, single-organism processes, metabolic processes, biological regulation, and regulation of biological processes were the top five biological processes involved in the greatest number of potential target mRNA. Cell, cell parts, and organelles were the top three cellular components involved in the greatest number of potential target mRNA. While binding and catalytic activity were the top two molecular functions involving the greatest number of potential target mRNA. The top 30 significantly enriched GO terms are also shown in .

Figure 2

GO enrichment analysis of miRNA target genes. The x-axis indicates GO terms, including biological process (blue), cell component (green), and molecular function (red). The y-axis includes three data series, percent of genes presented on the left, numbers of DEGs (black number), and background genes (gray number, representing all gene analyzed) on the right. The solid histogram represents DEGs and the virtual histogram represents background genes. GO, gene ontology; DEGs, differentially expressed genes.

Figure 3

The significantly enriched top 30 GO terms of biological process, cellular component, and molecular function according to P value and q-value (corrected P value). GO, gene ontology.

In the KEGG pathway analysis, the top 30 significantly enriched KEGG pathways are shown, and among them, pathways in cancer, proteoglycans in cancer, Rap1 signaling pathway, Ras signaling pathway, MAPK signaling pathway, Wnt signaling pathway, PI3K-Akt signaling pathway, ErbB signaling pathway, and TNF signaling pathway were the significant signaling pathways associated with the pathogenesis of cancer and involved in potential target mRNAs. The results are shown in .

Figure 4

The top 30 significant enrichment KEGG pathways of mRNAs according to P value and q-value (corrected P value). KEGG, Kyoto Encyclopedia of Genes and Genomes.

miRNA screen and qRT-PCR verification

According to the screening criteria for miRNAs and the attempt to verify more miRNAs, 19 miRNAs, including 11 upregulated and eight downregulated miRNAs, were selected among the differentially expressed miRNAs for further qRT-PCR verification. In the verification, 50 ESCC samples, including 25 without LNM and 25 with LNM were selected to verify the effectiveness of the screened miRNAs. The results of qRT-PCR verification showed that the expression tendency of 14 miRNAs was consistent with that of the microarray data, and the consistency was 73.68%. Among the 19 miRNAs verified, five miRNAs were significantly upregulated, no miRNAs were significantly downregulated in ESCC tissues with LNM compared with those without LNM. Among the five significantly upregulated miRNAs, the expression of miR-4440 and miR-6127 were inconsistent between the qRT-PCR verification data and microarray data, and they were therefore excluded. Finally, three obviously upregulated miRNAs, including miR-19a-3p, miR-17-3p, and miR-18b-5p were identified as potentially involved in the LNM of ESCC and were regarded as candidate miRNAs capable of predicting regional LNM in ESCC. The detailed data of the 19 miRNAs verified are showed in .

Table 3

Expression data of verified miRNAs in verified samples

Systematic name	Group 1† (mean ± SD)	Group 2§ (mean ± SD)	Expression tendency	P value	Consistency with microarray data
miR-19a-3p	2.7755±1.9353	4.3625±2.4464	Up	0.0142	Consistent
miR-17-3p	11.5405±7.1147	20.8601±11.6442	Up	0.0015	Consistent
miR-18b-5p	3.8370±3.5448	7.4825±4.7258	Up	0.0034	Consistent
miR-4440	15.0672±9.7576	28.5884±12.8204	Up	0.0001	Inconsistent
miR-6127	1.8697±2.4215	7.0070±10.2709	Up	0.0219	Inconsistent
miR-18a-5p	11.8070±17.2393	15.4644±12.1872	Up	0.3907	Consistent
miR-4728-5p	12.1957±8.2501	10.1191±5.7729	Down	0.3077	Consistent
miR-301a-3p	4.0238±5.3562	6.0518±4.9430	Up	0.1706	Consistent
miR-1587	2.8291±4.4520	2.0499±3.4587	Down	0.4929	Consistent
miR-6779-5p	11.1996±15.6308	6.5792±9.1691	Down	0.2100	Consistent
miR-1249-5p	11.5868±25.2595	1.9651±2.0145	Down	0.0696	Consistent
miR-25-3p	5.8477±7.8516	5.9716±15.2436	Up	0.9713	Consistent
miR-590-5p	1.7158±1.9584	1.9663±3.1066	Up	0.7346	Consistent
miR-106b-5p	3.1245±3.8459	5.2519±5.1641	Up	0.1056	Consistent
miR-1290	1.4002±1.2459	1.1800±1.3934	Down	0.5588	Consistent
miR-93-5p	2.4981±1.8751	4.1980±3.9970	Up	0.0626	Consistent
miR-197-5p	3.0791±4.4242	7.5768±11.3776	Up	0.0750	Inconsistent
miR-17-5p	3.8699±6.0198	3.1352±2.3148	Down	0.5716	Inconsistent
miR-20a-5p	2.3972±4.5983	1.5096±1.1455	Down	0.3537	Inconsistent

†, Group 1: expression levels of miRNAs in cancer tissues without lymph node metastasis; §, Group 2: expression levels of miRNAs in cancer tissues with lymph node metastasis.

Relations between candidate miRNAs and clinicopathologic features

To further verify the specificity of the three candidate miRNAs capable of predicting regional LNM in ESCC, the relations between candidate miRNAs and clinicopathologic features were also analyzed. The results showed that the three candidate miRNAs were not significantly associated with all other clinicopathologic features, except TNM stage, which was determined by the N stage, so the three candidate miRNAs could be regarded as tumor-specific miRNAs capable of predicting regional LNM of ESCC. All data are detailed in .

Table 4

The relations between candidate miRNAs and clinicopathological characteristics

Main features	n	hsa-miR-19a-3p			hsa-miR-17-3p			hsa-miR-18b-5p
		Mean ± SD	P value		Mean ± SD	P value		Mean ± SD	P value
Gender			0.784			0.612			0.797
Male	40	3.6147±2.4346			15.8129±11.0116			5.5763±4.5577
Female	10	3.3863±1.9273			17.7501±9.4024			5.9937±4.6209
Age (year)			0.156			0.613			0.954
<60	20	2.9944±2.0373			15.2544±10.2401			5.7061±4.8784
≥60	30	3.9521±2.4576			16.8309±11.0367			5.6289±4.3607
Primary location			0.457			0.93			0.382
Lower thoracic esophagus	26	3.8069±2.0788			16.3295±10.8015			6.2036±4.1978
Midthoracic esophagus	24	3.3113±2.5869			16.0603±10.7074			5.0706±4.8773
Pathology grade			0.401			0.593			0.756
Moderate and well	30	3.7973±2.4582			15.5333±10.5382			5.4949±3.9463
Poor	20	3.2266±2.1258			17.2009±11.0041			5.9071±5.3804
Clinical subtype			0.806			0.07			0.662
Ulcerative type	23	3.6581±1.9898			13.2506±8.8091			5.3519±4.0881
Medullary and other type	27	3.4932±2.6133			18.7131±11.5639			5.9219±4.9303
T stage			0.303			0.911			0.726
T1–2	10	2.8846±2.0499			16.5407±10.3503			5.2051±4.4572
T3	40	3.7401±2.3823			16.1152±10.8475			5.7735±4.5919
TNM stage†			0.008			0.002			0.002
I–II	26	2.7521±1.8999			11.7499±7.0523			3.8157±3.4749
III–IV	24	4.4539±2.4551			21.0215±11.8661			7.6575±4.7441
Tumor size (cm2)	–	−0.047§	0.745		0.102§	0.479		0.034§	0.815

†, postoperative pathological stage; §, Pearson correlation coefficient.

Potential mechanism of the tumor-specific miRNAs

For the three tumor-specific miRNAs capable of predicting the regional LNM of ESCC, we could not obtain the potential target genes of miR-19a-3p and miR-18b-5p, according to the above miRNA target gene prediction criteria, but 858 potential target mRNAs were obtained for miR-17-3p. The 858 target genes of miR-17-3p were enriched in 180 KEGG pathways, which included some pathways associated with cancer such as the Wnt signaling pathway, PI3K-Akt signaling pathway, AMPK signaling pathway, and MAPK signaling pathway, etc. All these were possibly involved in the miR-17-3p potential mechanism to regulate LNM in ESCC.

Discussion

The tendency toward extensive regional LNM is an important clinical characteristic of ESCC and also a major cause of failure in surgery and radiotherapy (5,6), which results in a disappointing 5-year overall survival rate (3,4). Therefore, the search for appropriate molecular biomarkers to predict regional LNM of ESCC has become the focus at present. Li et al. reviewed potential biomarkers for predicting regional LNM in ESCC in some studies in recent years (19). However, these biomarkers are still not able to guide the treatment of regional LNM in ESCC due to the lack of pertinence and specificity on LNM status. Based on the cognition that the molecular biological properties of cancer determine the biological behaviors of cancer and the tumor specificity of miRNA in tumorigenesis, development, invasion, and metastasis (12,15,16), this microarray analysis of miRNA based on LNM status obtained three tumor-specific miRNAs for predicting regional LNM in ESCC, which would guide the selective lymphatic dissection or irradiation of ESCC, change the treatment modalities of ESCC, and be conducive to bring great progress in the treatment of ESCC.

In this study, to assure the reliabilities of the experimental results, strict quality control measures were taken for sample selection, pathological sampling, sample storage, RNA isolation, and microarray analysis to eliminate any possible interference with the experimental results. Finally, we obtained the expected results, which were intuitively presented with cluster plot, Circos plot, volcano plot, and scatter plot as shown in .

A comparative microarray analysis identified 62 miRNAs that were differentially expressed in ESCC tissues with LNM compared to those without LNM. In combination with the target prediction analysis based on miRWalk2.0, 50 miRNAs and 8,637 target genes were considered to have potential regulatory functions in the regional LNM of ESCC. Subsequent GO and KEGG pathway analyses showed the functions of these potential target mRNAs and their enriched pathways. Several pathways were significantly involved in the pathogenesis of cancer. These results, on the one hand, exemplify the need for this experiment, and on the other hand, have provided valuable research evidence of a select number of miRNAs and their target genes for further studies.

Among the 62 miRNAs that had underlying regulatory function in the regional LNM of ESCC according to the screening criteria, 19 miRNAs, which accounted for almost one-third of differentially expressed miRNAs identified, were selected for further qRT-PCR verification. In the verification process, an additional 50 samples were used, therefore, the results of this confirmatory experiment not only verified the consistency between the qRT-PCR verification data and microarray data, but also verified the expression levels of these screened miRNAs in ESCC tissue. Finally, the results of this confirmatory experiment showed good consistency between the qRT-PCR verification data and the microarray data.

After excluding the inconsistent miRNAs between the qRT-PCR verification data and microarray data, miRNA-17-3p, miRNA-19a-3p, and miRNA-18b-5p were identified as significantly upregulated in ESCC. Further univariate analysis among the three miRNAs and clinicopathologic features showed that the three upregulated miRNAs were not significantly associated with all other clinicopathologic features except the TNM stage, which was determined by N status. These results indicated that miRNA-17-3p, miRNA-19a-3p, and miRNA-18b-5p were involved in the regulation of regional LNM of ESCC, and could be identified as tumor-specific miRNAs used for the prediction of regional LNM in ESCC.

Among the three tumor-specific miRNAs, we could not obtain the target genes of miRNA-19a-3p and miRNA-18b-5p from target prediction analysis based on miRWalk2.0. However, several studies (20-23) revealed that miRNA-19a-3p was involved in regulating the biological processes of a few benign diseases and malignant tumors via different mechanisms. Wang et al. reported that miRNA-19a-3p might contribute to the pathogenesis of preeclampsia during pregnancy via the lncRNA-PSG10P/miRNA-19a-3p/IL1RAP pathway (20). Chen et al. reported that miRNA-19a-3p promoted rheumatoid arthritis fibroblast-like synoviocytes via targeting SOCS3 (21). Bai et al. found that miRNA-19a-3p could enhance the malignancy grade of human ovarian cancer cells via inhibiting the expression of IGFBP-3 (22). Li et al. found miRNA-19a-3p could regulate cell growth and apoptosis in gastric cancer via the miRNA-19a-3p/PTEN/PI3K/AKT signaling pathway (23). Although there were many such studies, the function of miRNA-19a-3p in ESCC was not reported. With respect to miRNA-18b-5p, the number of studies is rare and only a few researchers have published reports (24-26). Zhang et al. found miR-18a-5p were significantly upregulated in colorectal cancer tissues and significantly elevated in colorectal cancer plasma exosomes, but its function and mechanism were unclear (24). Xue et al. confirmed that miRNA-18b-5p could regulate lung adenocarcinoma cell proliferation via the miR-18b-5p/VMA21 axis (25). Cochetti et al. found that miR-18b-5p owned the highest sensitivity and specificity for the prediction of prostate cancer, and decreased in expression concomitantly with an increase in malignancy, the mechanism was not studied (26). Finally, for tumor-specific miRNA-17-3p, 858 potential target mRNAs were obtained based on the present study criteria, and these target genes enriched many tumor-related pathways. Overall, this indicates that the function and regulatory mechanism of miRNA-17-3p are intricate and extensive in the pathogenesis of cancer. However, over the past ten years, studies associated with miRNA-17-3p focused more on benign diseases such as inhibiting angiogenesis by downregulating Flk-1 (27), inhibiting mouse cardiac fibroblast senescence by targeting Par4 (28), and promoting keratinocyte cell growth and metastasis in cutaneous wound healing by targeting MYOT and regulating the Notch1/NF-κB pathways (29), and so on. Moreover, only a few studies report results associated with miRNA-17-3p in cancer including prostate cancer (30), colon cancer (31), and glioblastoma (32). However, among the studies associated with miRNA-17-3p, there was no research related to ESCC. All the studies above indicated that the potential mechanisms of the three tumor-specific miRNAs regulating regional LNM of ESCC were worthy of being studied and provides evidence for further studies. Moreover, further function and mechanism studies will be conducive to confirm the application value of the three tumor-specific miRNAs and provide theoretical foundation for their clinical application.

Conclusions

In summary, the present microarray analysis of miRNAs based on the LNM status of ESCC identified three tumor-specific miRNAs for the prediction of regional LNM in ESCC and provides significant evidence for future studies of the potential mechanism, which warrants further investigation.