Skip to main content

Mir-421 and mir-550a-1 are potential prognostic markers in esophageal adenocarcinoma

Abstract

Objective

To identify the prognostic indicators of esophageal adenocarcinoma (EAC) for future EAC diagnosis and treatment.

Methods

The EAC dataset from The Cancer Genome Atlas was screened for differentially expressed microRNAs (miRNAs) and mRNAs associated with EAC. Weighted gene coexpression network analysis was performed to cluster miRNAs or mRNA with similar expression patterns to identify the miRNAs or mRNA that are highly associated with EAC. Prognostic miRNAs for overall survival (OS) were identified using Cox proportional-hazards regression analysis and least absolute shrinkage and selection operator based on survival duration and status. Two types of miRNAs were selected to develop a prognostic signature model for EAC using multiple Cox regression analysis. Furthermore, the signature was validated using internal validation sets 1 and 2. The receiver operating characteristic curve and concordance index were used to evaluate the accuracy of the signature and validation sets. The expression of miR-421, miR-550a-3p, and miR-550a-5p was assessed using quantitative polymerase chain reaction (qPCR). The proliferation, invasion, and migration of EAC cells were assessed using CCK8 and transwell assays. The OS of target mRNAs was assessed using Kaplan–Meier analysis. Functional enrichment analysis of the target mRNAs was performed using Metascape.

Results

The prognostic signature and validation sets comprising mir-421 and mir-550a-1 had favorable predictive power in OS. Compared with the patients with EAC in the high-expression group, those assigned to the low-expression group displayed increased OS according to survival analysis. Differential and qPCR analysis showed that miR-421, miR-550a-3p, and miR-550a-5p were highly expressed in the EAC tissues and cell lines. Moreover, the downregulation of miR-421 and miR-550a-3p with inhibitor markedly suppressed the proliferation, invasion, and migration in OE33 cells compared with the negative control. A total of 20 target mRNAs of three miRNAs were predicted, among which seven target mRNAs—ASAP3, BCL2L2, LMF1, PPM1L, PTPN21, SLC18A2, and NR3C2—had prognostic value; PRKACB, PDCD4, RPS6KA5, and BCL2L2 were enriched in the miRNA cancer pathway.

Conclusion

Prognostic indicators of EAC may be useful in future EAC diagnosis and treatment.

Introduction

The Global Cancer Statistics 2018 report revealed that esophageal carcinoma (ESCA), which is one of the most common malignant cancers worldwide, ranked seventh in terms of cancer incidence and sixth in terms of cancer-associated mortality globally [1]. ESCA can be divided into esophageal adenocarcinoma (EAC) and esophageal squamous carcinoma [2]. The incidence of EAC has increased substantially over the recent decades, and it remains a highly fatal disease with a poor prognosis [3, 4]. The prognostic outcomes for EAC is unfavorable, and the 5-year overall survival (OS) is < 15% despite the advancements in diagnosis and treatment [5, 6]. Exploring the mechanism that leads to the development and progression of EAC would be significant in understanding the poor outcomes of EAC. Thus, there is an urgent need to identify more novel and potential indicators that can help understand EAC and improve decision-making in clinical practice to effectively predict the occurrence, progression, and prognosis of this disease [7].

MicroRNAs (miRNAs) are noncoding RNAs that play an important role in cell differentiation, proliferation, and survival by binding to complementary target mRNAs, leading to the inhibition or degradation of mRNA translation [8]. Increasing evidence has indicated that miRNAs can be used in the diagnosis and prognosis of various tumors and diseases, including colorectal cancer, cutaneous melanoma, acute coronary syndrome, and acute myeloid leukemia; however, these are not part of standard clinical practice [9,10,11,12].

Studies reporting the miRNA biomarkers that improve the diagnosis and treatment of EAC are limited. Therefore, identification of new diagnostic and prognostic biomarkers is urgently needed because of the biological complexity and unfavorable prognosis of EAC.

In the past decade, novel biomarkers and therapeutic targets in various cancers have been explored using microarray and next-generation sequencing technologies [13,14,15]. However, the analyses of related makers are inadequate and even contradictory because of the different statistical data processing methods used [16]. To solve this problem, integrated bioinformatics methods such as weighted gene coexpression network analysis (WGCNA) have been used in various cancer studies; these were initially used to identify miRNA modules that are highly associated with EAC [17,18,19].

In this study, a prognostic model of a signature comprising two miRNAs was constructed based on a miRNA dataset from The Cancer Genome Atlas (TCGA); two prognostic miRNAs were screened for the prognosis of EAC. Moreover, the influence of two prognostic miRNAs on the proliferation, invasion, and migration of EAC cells was assessed. The Metascape database was explored for the possible cellular functions and pathways of target mRNAs related to this signature. Furthermore, seven prognostic target mRNAs were identified through Kaplan–Meier analysis. The entire analysis flow chart is shown in Fig. 1.

Fig. 1
figure 1

The entire analysis flow chart. QC Quality control, DEA Differentially expressed analysis, DEMs Differentially expressed miRNAs, UCox Univariate Cox, MCox Multivariate Cox

Materials and methods

Data analysis using the EAC cohort of TCGA

The miRNA-Seq data of 88 tumor tissues and 9 adjacent nontumorous tissues (Additional file 1: 1) and clinical information for EAC (Additional file 1: 2) were derived from TCGA (https://tcga-data.nci.nih.gov/tcga/) at the miRNA gene level rather than at mature miRNA level on October 19, 2020. For clinical data, clinical information from normal patients was deleted. Patient miRNAs with a mean value of expression = 0 were removed.

Screening of differentially expressed miRNAs and genes

To identify differentially expressed miRNAs (DEMs) and differentially expressed genes (DEGs) in EAC and normal samples, edgeR was performed to normalize the gene read counts to the trimmed mean of M values [20]. P–values were adjusted to minimize the false discovery rate (FDR) using the Benjamini & Hochberg method [21]. A |Fold change| of > 1 and adjusted P-value < 0.05 were considered significant [22].

To identify DEMs and DEGs associated with EAC, WGCNA [23]—an algorithm used to identify gene coexpression network via high-throughput expression profiles with different traits—was performed to cluster miRNAs with similar expression patterns for obtaining the most related miRNAs from patients with EAC. WGCNA with a softpower of 2 and module size cutoff of 15 were set as the threshold for miRNAs, whereas a softpower of 10 and module size cutoff of 100 were set as the threshold on the basis of the expression values of all mRNAs.

The VennDiagram package [24] was used to intersect the results of WGCNA and DEMs identified to screen for DEMs associated with EAC.

Construction of a multimicroRNA-based classifier

Univariate Cox proportional-hazards regression was used to identify the relationship between the DEM expression and OS of 88 patients. To prevent data overfitting, miRNAs with a P-value of < 0.05 (candidate variables) were rescreened via Lasso Cox regression analysis and the sample with a survival time of = 0 was removed [25]. Thereafter, stepwise multivariate Cox proportional-hazards regression analysis was performed to evaluate the contribution of the selected miRNAs. The regression methods were used to construct a classifier. Using a multivariate Cox regression model to calculate the risk scores for all patients and considering the median risk score as the cutoff value, the patients were divided into high- and low-risk groups. The analysis was performed using the survival R package (version 2.38, https://CRAN.R-project.org/package=survival) and rms (version 6.1.1, https://CRAN.R-project.org/package=rms). The R package “survivalROC” (version 1.0.3, https://CRAN.R-project.org/package=survivalROC), which performs receiver operating characteristic curve (ROC) and concordance index (CI) analyses using 1-, 2-, 3-, 4-, 5- OS data, was used to assess the performance of the prognostic signature model. The R package “caret” was used to randomly assign 43 of the patients to the internal training set and 44 to the internal testing set. Furthermore, the model of prognostic signature was validated using internal validation sets 1 and 2.

Kaplan–Meier survival analysis

The Kaplan–Meier plotter, a database that integrates the resources of TCGA and Gene Expression Omnibus, evaluates the effects of mRNAs, miRNAs, and proteins on patient survival in 21 cancer types (http://kmplot.com/). The Kaplan–Meier plotter database was used to analyze the effects of miRNAs and their target genes on patient survival in EAC.

Expressions in the GC subgroup and target gene prediction for microRNAs

The miRbase was used to obtain mature miRNAs of precursor miRNAs (premiRNA). The miRDB (http://www.mirdb.org/), R and Hiplot (https://hiplot.com.cn) were performed to reshape and visual data of expressions in the GC subgroup for mature miRNAs, respectively. The miRWalk (http://mirwalk.umm.uni-heidelberg.de/), miRTarBase (http://mirtarbase.mbc.nctu.edu.tw/), and TargetScan (http://www.targetscan.org/) databases were used to explore target mRNAs. The Cytoscape software (Version 3.7.2; http://www.cytoscape.org/) was used to visualize and construct the miRNA–mRNA network.

Functional enrichment analysis

To determine the biological functions and potential signaling pathways of overlapping target mRNAs, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed using the Metascape database (http://metascape.org/gp/index.html#/main/step1).

Cell culture and real-time quantitative polymerase chain reaction

A normal human esophageal epithelium cell line (Heec) purchased from the Chinese National Infrastructure of Cell Line Resource, and two EAC cell lines (OE33 and SEG-1) were purchased from the Chinese Binsui Bio for in vitro validation.

OE33 cells were cultured in complete medium comprising 90% Dulbecco's modified Eagle’s medium (Invitrogen, Carlsbad, USA) and 10% fetal bovine serum (HyClone, USA). SEG-1 and Heec cells were cultured in Roswell Park Memorial Institute Medium-1640 (Invitrogen, Carlsbad, CA, USA).

TRIzol (Solarbio, China) was used to isolate total RNA from whole-cell lysates using the PrimeScript™ RT Reagent Kit (TaKaRa, USA). Bulge-LoopTM miRNA qRT-PCR primers were used to synthesize cDNA according to the manufacturer’s instructions. Real-time quantitative polymerase chain reaction (qRT-PCR) was performed using Takara Bio Green Premix Ex Taq™ II (Tli RNaseH Plus, TaKaRa, USA) following the manufacturer’s instructions. U6 was used as the endogenous control. All primers were obtained from RiboBio Co., Ltd., China.

miRNA transfection

miR-421 and miR-550a-3p inhibitors and corresponding controls were synthesized and purchased from RiboBio. Cells were cultured overnight in six-well plates (2 × 105 cells/well) and then transfected with miR-421 inhibitor (concentration, 100 nM), miR-550a-3p (100 nM), and corresponding negative control (100 nM) when the cell confluence reached 30–40%. The transfection reagent LipofectamineTM3000 (Invitrogen, USA) was used.

CCK8 assay

Approximately 2 × 103 cells per chamber were cultured in 24-well plates, and then the CCK8 assay was performed using the CCK8 kit (Solarbio, China). A microplate reader (Bio-Rad, USA) was used to measure the absorbance at an optical density of 450 nm.

Transwell invasion and migration assay

The mobility of transfected cells was detected using transwell chambers (Corning, China); for migration, 5 × 104 cells per chamber were cultured in transwell chambers without basement membrane coating. For invasion, 1 × 105 cells/chamber were cultured in transwell chambers containing a basement membrane coating. Afterward, the cells were incubated for 24 h in an incubator at 37 °C with 5% CO2; the cells of the inner chamber were wiped using a swab and cells on the outer chamber were fixed using methanol for 15 min and stained using 0.5% crystal violet (Solarbio, China) for 30 min. Then, the cells were washed three or four times with phosphate buffered saline and observed under an inverted phase microscope (ZEISS, Germany) to obtain field images.

Statistical analysis

The pairwise Pearson correlation was used to assess the weighted coexpression relationship among all dataset subjects in an adjacency matrix. Logrank P-value was used to analyze the differences between the low- and high-risk/expression groups in survival analyses. Differences between the variables were considered significant if |log2FC|≥ 1 and adjusted P-value < 0.05 in difference expression analysis or A fold change of > 1 and P-value < 0.05 in qRT-PCR analysis. The R software (version 4.0.1) was used.

Results

Identification of DEMs associated with EAC

DEMs obtained from TCGA included 52 upregulated and 56 downregulated miRNAs (Fig. 2A–B). Using WGCNA, 16 coexpression modules were identified, among which a brown module, which included 60 miRNAs, showed a high association with EAC (Fig. 2C). To obtain the DEMs associated with EAC, 28 miRNAs were obtained using data between DEMs and the brown module, including 17 upregulated and 11 downregulated miRNAs identified using the VennDiagram package in R (Fig. 2D).

Fig. 2
figure 2

Identification of differentially expressed microRNAs (DEMs) associated with esophageal adenocarcinoma (EAC). A Volcano plots of DEMs. B Heatmap of DEMs in both EAC and adjacent normal samples. C Relationship of DEMs between the modules of the normal samples and those of EAC samples. D Venn plots of microRNAs between DEMs and microRNAs in the brown module

Identification of prognostic microRNA signatures

A total of 4 out of 28 significant miRNAs associated with the expression levels in patients with EAC were found to be significantly related to the OS on univariate Cox regression analysis (Fig. 3A). Following initial filtration, Lasso penal Cox analysis with penalty parameters adjusted by tenfold cross-validation was performed to further narrow the mRNAs by selecting those that occurred > 900 times out of a total of 1000 replicates (Fig. 3B–C). Thus, three miRNAs were identified. Multivariate Cox proportional-hazards regression model analysis revealed that the prognostic signature contained the following two genes: has-mir-550a-1 and has-mir-421 (Table 1).

Fig. 3
figure 3

Identification of prognostic microRNAs (miRNAs). A Univariate Cox regression analysis of four differentially expressed miRNAs (DEMs). B, C Lasso-penalized Cox regression analysis of four DEMs. D Multivariate Cox regression analysis of three DEMs

Table 1 Prognostic value of the two microRNAs in patients with esophageal adenocarcinoma in The Cancer Genome Atlas cohort

Interestingly, all miRNAs had a hazard ratio of > 1 and were considered prognostic risk factors (Fig. 3D). The results of prognostic risk score, survival status, and expression of the two OS-related miRNAs for patients with EAC were displayed in (Additional file 3: Fig. 1A–C). Kaplan–Meier survival analysis indicated that the OS of patients in the low-risk groups was remarkably high compared with that of those in the high-risk groups, implying that the high-risk groups were associated with a poor prognosis (Fig. 4A–C). To evaluate the prognostic model, the 1-, 2-, 3-, 4- and 5-year ROC curves and CI values of entirety as well as internal validation sets 1 and 2 were plotted. The area under the curve (AUC) values of 1-, 2-, 3-, 4-, and 5-year entirety were 0.71, 0.72, 0.71, 0.72, and 0.77, respectively (Fig. 4D). The AUC of internal validation sets 1 and 2 were > 0.7. Furthermore, the CI values of entirety as well as validation sets 1 and 2 were 0.69, 0.68, and 0.67, respectively. (Fig. 4E). These results implied that the prognostic gene signature performed well in survival prediction.

Fig. 4
figure 4

Kaplan–Meier curves and prediction of signature model and validation sets for the low- and high-risk groups

OS of prognostic miRNAs

Kaplan–Meier analysis based on optimal cut-off value of miRNAs expression was used to explore the prognostic information of mir-421 and mir-550a-1 in order to identify OS. The OS of two miRNAs was poorer in the high-expression group than in the low-expression group (P < 0.05; Fig. 5A–B). To further ensure the accuracy of the results, the median was designed to divide expression grouping of EAC patients. The results showed that patients with high expression had a lower overall survival than those with low expression (Additional file 3: Fig. 2). Previous results have shown that mir-421 and mir-550a-1 show a higher expression in EAC samples than in paracancerous samples (Fig. 5C–D). Therefore, these data indicate that mir-421 and mir-550a-1 have prognostic value for patients with EAC.

Fig. 5
figure 5

Overall survival (OS) and expression analyses. A, B OS of mir-421 and mir-550a-1, respectively. C, D Expression value of mir-421 and mir-550a-1, respectively

Verification of the expression of miR-421, miR-550a-3p, and miR-550a-5p

The mature miRNAs of premiRNAs, miR-421, miR-550a-3p, and miR-550a-5p were identified using miRbase database, and the expression of Mature mir-421 and mir-550a-1 in the EAC subgroup was analyzed. As a result, they expressions are correlated with the vast majority of tumor grade, stage, M, N, T and gender compared with normal samples (Additional file 3: Fig. 3 and Additional file 3: Fig. 4). Importantly, the expression of miR-421, miR-550a-3p, and miR-550a-5p was verified using qRT-PCR analysis of Heec, OE33, and SEG-1 cells. The expression of these miRNAs was higher in OE33 and SEG-1 cells than in Heec cells (Fig. 6). Because the miR-550a-3p expression was the highest in miR550a-3p and miR-550a-5p, miR-550a-3p and miR-421 were selected for the subsequent experiments.

Fig. 6
figure 6

Expression of prognostic microRNAs indicators in the cell lines

Suppression of miR-421 and miR550a-3p inhibited the proliferation, invasion, and migration of EAC cells

miR-421 and miR-550a-3p inhibitors were transfected into OE33 cells. The results showed that the downregulation of miR-421 and miR-550a-3p with inhibitor markedly suppressed the proliferation (Fig. 7A–D), migration (Figs. 7E and F), and invasion (Figs. 7G and H) of OE33 cells compared with the negative control.

Fig. 7
figure 7

Effects of the suppression of miR-421 and miR-550a-3p on the proliferation, invasion, and migration of OE33 cells. Paired t test *P < 0.05, **P < 0.01, and ***P < 0.001

Identification of DEGs associated with EAC

Differential analysis showed that 3,390 mRNAs were DEGs according to the critical value (FDR < 0.05 and |log2FC|> 1). Of these, 1,155 were upregulated and 2,235 were downregulated (Additional file 3: Fig. 5A). Moreover, the distribution of the significant DEGs was assessed via heatmapping (Additional file 3: Fig. 5B).

Five mRNA modules were identified using WGCNA (Additional file 3: Fig. 5C). Among these five mRNA modules, the blue module had a strongly-connected gene coexpression network with EAC (Additional file 3: Fig. 5C). A total of 610 DEGs associated with EAC were obtained from the genes between DEGs and the blue module, including 432 upregulated and 178 downregulated mRNAs (Additional file 3: Fig. 5D).

Construction of the miRNA–mRNA network and functional enrichment analysis of downregulated target mRNAs

Upregulated miRNAs are known to promote or inhibit the occurrence and progression of tumors by downregulated target mRNAs. Therefore, we used these downregulated target genes for subsequent analysis. The downstream target mRNAs with reference to the three mature miRNAs are identified using miRDB, miRTarBase, miRWalk, and TargetScan. Moreover, potential mRNAs, which were only shared by the downregulated DEGs associated with EAC and at least two databases, were selected to enhance the accuracy of the prediction (Fig. 8A–C). The results revealed an association between the miRNAs and their corresponding 20 target mRNAs (Fig. 8D).

Fig. 8
figure 8

Prediction of target mRNAs and establishment of a microRNA (miRNA)–mRNA axis and Gene ontology enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of the 20 target mRNAs. A, B, and C overlapped target mRNAs were analyzed using the predicted target mRNAs, differentially expressed genes (DEGs) associated with esophageal adenocarcinoma, and significantly-downregulated mRNAs, respectively. D miRNA–mRNA network was constructed using the three miRNAs for EAC prognosis and 20 overlapped target mRNAs. E Biological process. F Molecular function. G KEGG pathways

GO enrichment analysis of biological processes indicated that the target mRNAs of three maturing miRNAs were mainly enriched in behavior, cellular response to lipid, cellular response to hormone stimulus, dephosphorylation, and positive regulation of defense response (Fig. 8E). Figure 8E shows the enriched molecular function GO terms. Unfortunately, cellular components were not enriched when the P-value was < 0.01. The KEGG pathway analysis revealed that the target mRNAs participate in miRNAs in cancer pathways (Fig. 8G). Therefore, these target mRNAs may be associated with cellular response, transmembrane transporter activities, and miRNAs in cancer pathways.

Survival analysis of downregulated target mRNAs

To explore the prognostic value of downregulated target mRNAs, their OS was assessed using the Kaplan–Meier plotter database. Patients in the low-expression group had a reduced OS compared with those in the high-expression group when the P-value was < 0.05 (Fig. 9A–G and Table 2), except for ELAVL2 (Fig. 9H and Table 2), where the seven target mRNAs in the EAC and adjacent samples and those of they all were poorly expressed (Fig. 9I). Regarding the other target mRNAs, these genes had no prognostic value (P < 0.05; Additional file 3: Table 1). Therefore, these data suggested that the signature of the following seven downregulated target mRNAs—ASAP3, BCL2L2, LMF1, PPM1L, PTPN21, SLC18A2, and NR3C2—have a prognostic value in EAC.

Fig. 9
figure 9

Overall survival of the target genes and their expression in EAC and adjacent tissues. Wilcoxon test *P < 0.05, **P < 0.01, and ***P < 0.001

Table 2 Prognostic value of the eight mRNAs in patients with esophageal adenocarcinoma

Discussion

In this study, we used bioinformatics, particularly WGCNA, to identify 28 DEMs associated with EAC; two of these miRNAs were then selected to establish the prognostic signature model of EAC via multiple Cox regression analysis.

This study discovered two crucial miRNAs associated with a poor EAC prognosis and constructed a prognostic model comprising these miRNAs that had upregulated expression in the EAC tissues. For the two prognostic miRNAs (mir-550a-1 and mir-421), the mature miRNAs of mir-550a-1, mir-421, miR-421, miR-550a-3p, and miR-550a-5p were identified using the miRbase database. The suppression of miR-421 and miR-550a-3p inhibited the proliferation, invasion, and migration of EAC cells. The unfavorable prognosis of miR-421 in EAC has been reported previously, and upregulated miR-421 expression has been identified as an unfavorable prognostic marker in EAC [26]. The overexpression of miR-421 facilitates cell proliferation in non-small-cell lung cancer (NSCLC) [27]. In addition, the association of miR-421 with various tumors has garnered growing attention. For example, the upregulated miR-421 expression has been shown to inhibit the proliferation and metastasis of colorectal cancer by targeting MTA1 [28]. Moreover, circAHNAK1 suppresses the proliferation and metastasis of triple-negative breast cancer cells by modulating the miR-421 expression (29). Similarly, circSETD3 acts as a sponge for miR-421 and inhibits the growth of hepatocellular carcinoma [30]. N-Myc has been shown to promote the development of neuroendocrine therapy resistance in prostate cancer through the differential regulation of the miR-421 pathway [31]. Although the prognostic value of miR-421 in EAC and the mechanism of miR-421 in other cancers have been validated, the mechanism through which this operates in EAC remains unclear. The present study shows that the high expression of miR-421 in EAC cells may be a useful prognostic indicator in EAC.

Regarding the two mature miRNAs of mir-550a-1 (miR-550a-3p and miR-550a-5p), to the best of our knowledge, the association of miR-550a-3p and miR-550a-5p with EAC has not been studied previously. miR-550a-3p reportedly shows high expression and is associated with significantly reduced OS in melanoma [32]. The upregulated expression of miR-550a-3p significantly promotes the cellular proliferation, invasion, and migration in NSCLC, whereas the knockdown of the miR-550a-3p expression inhibits cancer growth and metastasis [33]. Recent studies have shown that high expression of miR‐550a‐3p reverses the inhibitory effect of the increased LINC00261 expression, and the reduced SDPR expression reverses the growth-promoting effect of miR‐550a‐3p in breast cancer stem cells [34]. Regarding miRNA-550a-5p, a previous study demonstrated that patients with lung adenocarcinoma in the miR-550a-5p high-expression group had a shorter OS than those in the low-expression group [35]. Transfection of an miR-550a-5p mimic promotes cell migration in colorectal cancer [36]. Recent evidence suggests that miR-550a-5p shows high expression and promotes tumor proliferation by binding to LIMD1 in lung adenocarcinoma [37]. Furthermore, the in vitro expression of miR-550a-3p and miR-550a-5p was higher in EAC cells than in esophageal epithelial cells; however, the expression of miR-550a-3p was more significant, suggesting that miR-550a-3p plays an important role in the prognosis, occurrence, and progression of EAC. Nonetheless, these findings must be validated in future studies.

Metascape database analysis revealed that the target mRNAs PRKACB and PDCD4 of miR-550a-3p as well as RPS6KA5 and BCL2L2 of miR-550a-5p were associated with miRNA cancer pathway, providing important clues for subsequent studies on the mechanism underlying EAC development. Unfortunately, no target gene for miR-421 was enriched in miRNA cancer pathway. Nevertheless, these findings should be further explored in subsequent studies. Furthermore, seven prognostic target mRNAs—ASAP3, BCL2L2, LMF1, PPM1L, PTPN21, SLC18A2, and NR3C2—were obtained through Kaplan–Meier analysis. Among these mRNAs, the prognostic values and mechanisms of LMF1, PPM1L, and PTPN21 have rarely been reported in cancers. To our best knowledge, the correlation between the four other prognostic target mRNAs and EAC has not been studied previously. miRNA-143-3p reportedly inhibits the metastases of colorectal cancer by targeting ASAP3 [38]. The inhibition of BCL2L2 results in the high expression of miR-335-5p, which increases cisplatin sensitivity in ovarian cancer [39]. The low expression of SLC18A2 reduces OS in prostate cancer [40]. Downregulated NR3C2 expression is associated with a poor prognosis and aggressive characteristics in nondistant metastatic clear-cell renal cell carcinoma [41].

Overall, the consistency of our findings with previous studies confirms the high reliability of the analysis methods used in this study. However, the detailed mechanisms of miR-421, miR-550a-3p, and miR-550a-5p that influence the tumorigenesis and development of EAC by targeting mRNAs require urgent attention.

In conclusion, we analyzed the EAC data obtained from TCGA using combined bioinformatic approaches and cell experiments and identified prognostic indicators of two miRNAs associated with seven prognostic target mRNAs as the possible indicators for future EAC diagnosis and treatment.

Availability of supporting data

The original contributions presented in the study are included in the article/Additional files Material, further inquiries can be directed to the corresponding author/s.

References

  1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.

    Article  PubMed  Google Scholar 

  2. Talukdar FR, di Pietro M, Secrier M, Moehler M, Goepfert K, Lima SSC, et al. Molecular landscape of esophageal cancer: implications for early detection and personalized therapy. Ann N Y Acad Sci. 2018;1434(1):342–59.

    Article  CAS  PubMed  Google Scholar 

  3. Yoon HH, Shi Q, Sukov WR, Lewis MA, Sattler CA, Wiktor AE, et al. Adverse prognostic impact of intratumor heterogeneous HER2 gene amplification in patients with esophageal adenocarcinoma. J Clin Oncol. 2012;30(32):3932–8.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Ong CA, Shapiro J, Nason KS, Davison JM, Liu X, Ross-Innes C, et al. Three-gene immunohistochemical panel adds to clinical staging algorithms to predict prognosis for patients with esophageal adenocarcinoma. J Clin Oncol. 2013;31(12):1576–82.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Li X, Francies HE, Secrier M, Perner J, Miremadi A, Galeano-Dalmau N, et al. Organoid cultures recapitulate esophageal adenocarcinoma heterogeneity providing a model for clonality studies and precision therapeutics. Nat Commun. 2018;9(1):2983.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Rubenstein JH, Shaheen NJ. Epidemiology, diagnosis, and management of esophageal adenocarcinoma. Gastroenterology. 2015;149(2):302-17.e1.

    Article  PubMed  Google Scholar 

  7. Lan T, Lu Y, Xiao Z, Xu H, He J, Hu Z, et al. A six-microRNA signature can better predict overall survival of patients with esophagus adenocarcinoma. PeerJ. 2019;7: e7353.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Rupaimoole R, Slack FJ. MicroRNA therapeutics: towards a new era for the management of cancer and other diseases. Nat Rev Drug Discov. 2017;16(3):203–22.

    Article  CAS  PubMed  Google Scholar 

  9. Hibner G, Kimsa-Furdzik M, Francuz T. Relevance of MicroRNAs as potential diagnostic and prognostic markers in colorectal cancer. Int J Mol Sci. 2018;19(10):2944.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Ross CL, Kaushik S, Valdes-Rodriguez R, Anvekar R. MicroRNAs in cutaneous melanoma: role as diagnostic and prognostic biomarkers. J Cell Physiol. 2018;233(7):5133–41.

    Article  CAS  PubMed  Google Scholar 

  11. Barraclough JY, Joan M, Joglekar MV, Hardikar AA, Patel S. MicroRNAs as prognostic markers in acute coronary syndrome patients—a systematic review. Cells. 2019;8(12):1572.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Bhayadia R, Krowiorz K, Haetscher N, Jammal R, Emmrich S, Obulkasim A, et al. Endogenous tumor suppressor microRNA-193b: therapeutic and prognostic value in acute myeloid Leukemia. J Clin Oncol. 2018;36(10):1007–16.

    Article  CAS  PubMed  Google Scholar 

  13. Wang Z, Yang B, Zhang M, Guo W, Wu Z, Wang Y, et al. lncRNA epigenetic landscape analysis identifies EPIC1 as an oncogenic lncRNA that Interacts with MYC and promotes cell-cycle progression in cancer. Cancer Cell. 2018;33(4):706-20.e9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Sihag S, Nussenzweig SC, Walch HS, Hsu M, Tan KS, Sanchez-Vega F, et al. Next-generation sequencing of 487 esophageal adenocarcinomas reveals independently prognostic genomic driver alterations and pathways. Clin Cancer Res. 2021;27(12):3491–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Hoppe S, Jonas C, Wenzel MC, Velazquez Camacho O, Arolt C, Zhao Y, et al. Genomic and transcriptomic characteristics of esophageal adenocarcinoma. Cancers (Basel). 2021;13(17):4300.

    Article  CAS  PubMed  Google Scholar 

  16. Liu X, Wu J, Zhang D, Bing Z, Tian J, Ni M, et al. Identification of potential key genes associated with the pathogenesis and prognosis of gastric cancer based on integrated bioinformatics analysis. Front Genet. 2018;9:265.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Niemira M, Collin F, Szalkowska A, Bielska A, Chwialkowska K, Reszec J, et al. Molecular signature of subtypes of non-small-cell lung cancer by large-scale transcriptional profiling: identification of key modules and genes by weighted gene co-expression network analysis (WGCNA). Cancers (Basel). 2019;12(1):37.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Wang M, Wang L, Pu L, Li K, Feng T, Zheng P, et al. LncRNAs related key pathways and genes in ischemic stroke by weighted gene co-expression network analysis (WGCNA). Genomics. 2020;112(3):2302–8.

    Article  CAS  PubMed  Google Scholar 

  19. Long J, Huang S, Bai Y, Mao J, Wang A, Lin Y, et al. Transcriptional landscape of cholangiocarcinoma revealed by weighted gene coexpression network analysis. Brief Bioinform. 2020. https://doi.org/10.1093/bib/bbaa224.

    Article  PubMed  Google Scholar 

  20. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.

    Article  CAS  PubMed  Google Scholar 

  21. Kovacs B, Lumayag S, Cowan C, Xu S. MicroRNAs in early diabetic retinopathy in streptozotocin-induced diabetic rats. Invest Ophthalmol Vis Sci. 2011;52(7):4402–9.

    Article  CAS  PubMed  Google Scholar 

  22. Lai Y. A statistical method for the conservative adjustment of false discovery rate (q-value). BMC Bioinform. 2017;18(Suppl 3):69.

    Article  Google Scholar 

  23. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinform. 2008;9:559.

    Article  Google Scholar 

  24. Chen H, Boutros PC. VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinform. 2011;12:35.

    Article  Google Scholar 

  25. Friedman J, Hastie T, Tibshirani R. Regularization paths for generalized linear models via coordinate descent. J Stat Softw. 2010;33(1):1–22.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Lin XF, Zhang CQ, Dong BR. MiR-421 expression independently predicts unfavorable overall survival in patients with esophageal adenocarcinoma. Eur Rev Med Pharmacol Sci. 2019;23(9):3790–8.

    PubMed  Google Scholar 

  27. Li X, Chen SH, Zeng JW. MiR-421 Is overexpressed and promotes cell proliferation in non-small cell lung cancer. Med Princ Pract. 2020;29(1):80–9.

    Article  PubMed  Google Scholar 

  28. Xue L, Yang D. MiR-421 inhibited proliferation and metastasis of colorectal cancer by targeting MTA1. J buon. 2018;23(6):1633–9.

    PubMed  Google Scholar 

  29. Xiao W, Zheng S, Zou Y, Yang A, Xie X, Tang H, et al. CircAHNAK1 inhibits proliferation and metastasis of triple-negative breast cancer by modulating miR-421 and RASA1. Aging (Albany NY). 2019;11(24):12043–56.

    Article  CAS  PubMed  Google Scholar 

  30. Xu L, Feng X, Hao X, Wang P, Zhang Y, Zheng X, et al. CircSETD3 (Hsa_circ_0000567) acts as a sponge for microRNA-421 inhibiting hepatocellular carcinoma growth. J Exp Clin Cancer Res. 2019;38(1):98.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Yin Y, Xu L, Chang Y, Zeng T, Chen X, Wang A, et al. Correction to: N-Myc promotes therapeutic resistance development of neuroendocrine prostate cancer by differentially regulating miR-421/ ATM pathway. Mol Cancer. 2019;18(1):107.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Xiong J, Xue Y, Xia Y, Zhao J, Wang Y. Identification of key microRNAs of plasma extracellular vesicles and their diagnostic and prognostic significance in melanoma. Open Med (Wars). 2020;15(1):464–82.

    Article  CAS  PubMed  Google Scholar 

  33. Yang JZ, Bian L, Hou JG, Wang HY. MiR-550a-3p promotes non-small cell lung cancer cell proliferation and metastasis through down-regulating TIMP2. Eur Rev Med Pharmacol Sci. 2018;22(13):4156–65.

    PubMed  Google Scholar 

  34. Li Y, Wu C. LINC00261/microRNA-550a-3p/SDPR axis affects the biological characteristics of breast cancer stem cells. IUBMB Life. 2021;73(1):188–201.

    Article  CAS  PubMed  Google Scholar 

  35. Lin Y, Lv Y, Liang R, Yuan C, Zhang J, He D, et al. Four-miRNA signature as a prognostic tool for lung adenocarcinoma. Onco Targets Ther. 2018;11:29–36.

    Article  PubMed  Google Scholar 

  36. Wang G, Fu Y, Yang X, Luo X, Wang J, Gong J, et al. Brg-1 targeting of novel miR550a-5p/RNF43/Wnt signaling axis regulates colorectal cancer metastasis. Oncogene. 2016;35(5):651–61.

    Article  PubMed  Google Scholar 

  37. Guo ZZ, Ma ZJ, He YZ, Jiang W, Xia Y, Pan CF, et al. miR-550a-5p functions as a tumor promoter by targeting LIMD1 in Lung adenocarcinoma. Front Oncol. 2020;10: 570733.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Guo L, Fu J, Sun S, Zhu M, Zhang L, Niu H, et al. MicroRNA-143-3p inhibits colorectal cancer metastases by targeting ITGA6 and ASAP3. Cancer Sci. 2019;110(2):805–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Liu R, Guo H, Lu S. MiR-335-5p restores cisplatin sensitivity in ovarian cancer cells through targeting BCL2L2. Cancer Med. 2018;7(9):4598–609.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Haldrup C, Lynnerup AS, Storebjerg TM, Vang S, Wild P, Visakorpi T, et al. Large-scale evaluation of SLC18A2 in prostate cancer reveals diagnostic and prognostic biomarker potential at three molecular levels. Mol Oncol. 2016;10(6):825–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Zhao Z, Zhang M, Duan X, Deng T, Qiu H, Zeng G. Low NR3C2 levels correlate with aggressive features and poor prognosis in non-distant metastatic clear-cell renal cell carcinoma. J Cell Physiol. 2018;233(10):6825–38.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the support from Henan Key Laboratory for Pharmacology of Liver Diseases and the financial support from the National Natural Science Foundation of China (3160067), Project of Basic Research Fund of Henan Institute of Medical and Pharmacological Sciences (2021BP0101), Key Scientific Research Projects of Colleges and Universities in Henan Province (21A416011), and Medical Science and Technique Project of Henan Province (SBGJ202002104).

Funding

This work was supported by the National Natural Science Foundation of China (3160067), Project of Basic Research Fund of Henan Institute of Medical and Pharmacological Sciences (2021BP0101), Key Scientific Research Projects of Colleges and Universities in Henan Province (21A416011), and Medical Science and Technique Project of Henan Province (SBGJ202002104).

Author information

Authors and Affiliations

Authors

Contributions

Yun Ji and Liguo Zhang conceived and planned the study. Yun Ji performed majority of bioinformatic analysis. Yun Ji and Lulu Wang prepared the figures and wrote the first draft of the manuscript. Lulu Wang and Guanglei Chang performed majority of experiments, Juan Yan and Jingjing Liu participated in the identification of differentially expressed genes associated with patients with esophageal adenocarcinoma. Liping Dai, Zhenyu Ji, Meixia He and Hongliang Xu help Kaplan–Meier survival analysis and Functional enrichment analysis. Yun Ji, Lulu Wang, and Liguo Zhang interpreted the results. All authors contributed to the article and approved the submitted version. All authors read and approved by the final manuscript.

Corresponding author

Correspondence to Liguo Zhang.

Ethics declarations

Ethical approval and consent to participate

Not applicable.

Consent for publication

All the authors have consented to publication.

Competing interests

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Human ethics

Not applicable.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Esophageal adenocarcinoma prognostic markers-1.

Additional file 2.

Esophageal adenocarcinoma prognostic markers-2.

Additional file 3.

Supplementary Material.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ji, Y., Wang, L., Chang, G. et al. Mir-421 and mir-550a-1 are potential prognostic markers in esophageal adenocarcinoma. Biol Direct 18, 5 (2023). https://doi.org/10.1186/s13062-022-00352-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13062-022-00352-8

Keywords