Biochemistry Research International

Biochemistry Research International / 2020 / Article

Research Article | Open Access

Volume 2020 |Article ID 6285987 | https://doi.org/10.1155/2020/6285987

Ping Yang, Junjie Liu, Tianci Yang, Lei Zhang, Peiyou Gong, Boqing Li, Xiuzhi Zhou, "Construction and Investigation of MicroRNA-mRNA Regulatory Network of Gastric Cancer with Helicobacter pylori Infection", Biochemistry Research International, vol. 2020, Article ID 6285987, 10 pages, 2020. https://doi.org/10.1155/2020/6285987

Construction and Investigation of MicroRNA-mRNA Regulatory Network of Gastric Cancer with Helicobacter pylori Infection

Academic Editor: Bernardo Trigatti
Received23 Jan 2020
Revised04 Jun 2020
Accepted19 Jun 2020
Published25 Jul 2020

Abstract

Background. Helicobacter pylori (H. pylori) is a common human pathogen, which is closely correlated with gastric cancer (GC). However, the mechanism of H. pylori-related GC has not been elucidated. This study aimed to explore the role of H. pylori infection in GC and find biomarkers for early diagnosis of H. pylori-related GC. Methods. We identified differentially expressed microRNAs (DEMs) and genes (DEGs) from the Gene Expression Omnibus (GEO) dataset, constructed microRNA-(miRNA-)mRNA expression networks, analyzed the function and signal pathway of cross-genes, analyzed the relations between cross-genes and GC prognosis with the Cancer Genome Atlas (TCGA) data, and verified the expression of cross-genes in patients with H. pylori infection. Results. 22 DEMs and 68 DEGs were identified in GSE197694 and GSE27411 dataset. 16 miRNAs and 509 genes were involved in the expression network, while the cross-genes of the network were mainly enriched in MAP kinase (MAPK) signaling pathway and TGF-beta signaling pathway. Patients with higher expression of hsa-miR-196b-3p, CALML4, or SMAD6 or lower expression of PITX2 or TGFB2 had better outcomes than those with lower expression of hsa-miR-196b-3p, CALML4, or SMAD6 or higher expression of PITX2 or TGFB2 (). Patients with H. pylori infection had a higher expression of hsa-miR-196b-3p and CALML4 than those without H. pylori infection (). Conclusion. The study of miRNA-mRNA expression network would provide molecular support for early diagnosis and treatment of H. pylori-related GC.

1. Introduction

Gastric cancer (GC) is one of the most common cancers in China, with half a million deaths annually [1]. Because of population growth and life extension, the incidence and mortality rates of GC are increasing [2]. Therefore, it is urgent and important to identify the key genes in its pathogenesis. Helicobacter pylori (H. pylori) is acknowledged as a class I carcinogen [3], which colonizes in the gastric mucosa and causes chronic gastric, atrophic gastric, and GC [4]. However, how H. pylori infection is involved in the pathogenesis and development of GC is unknown [46]. Therefore, elucidating the molecular mechanism of H. pylori-related GC is of great importance for its early diagnosis and targeted therapy.

MicroRNAs (miRNAs) are 19∼25-nucleotide-long endogenous noncoding RNAs, which negatively regulate the expression of their targets at the posttranscription level and play significant roles in many biological processes, such as proliferation, differentiation, and apoptosis [79]. Evidence has been increasing that abnormal miRNAs expression is involved in the pathogenesis and development of many cancers, which suggests the promising biomarkers for early diagnosis and therapy of tumors [5, 10]. Now, a high-throughput platform combined with bioinformatics analysis has become a new way to identify biomarkers of disease [11, 12].

In this study, differentially expressed miRNAs (DEMs) and differentially expressed genes (DEGs) of gastric biopsy with H. pylori infection from the Gene Expression Omnibus (GEO) database [13, 14] were identified by R software. After predicting the potential targets of DEMs, we constructed the coexpression network of miRNA-mRNA to identify hub genes. The hub genes were identified using bioinformatics methods including gene ontology (GO) annotation [15] and Kyoto Encyclopedia of Genes and Genomes (KEGG) [16] signal pathway enrichment analysis, while the prognostic value of hub genes was analyzed from the Cancer Genome Atlas (TCGA) database, and the expression of hub genes was confirmed in 69 gastric specimens with or without H. pylori infection. We hope this study will provide new information for the molecular mechanism of H. pylori-related GC.

2. Materials and Methods

2.1. Microarray Data

The miRNA and gene expression profiles were obtained from the GEO dataset (https://www.ncbi.nlm.nih.gov/geo/). The screening criteria for GEO were as follows: (1) human gastric mucosa samples with or without H. pylori infection; (2) datasets were raw or standardized. The miRNA expression profile (GSE19769) [7] and the gene expression profile (GSE27411) [17] were included in this study. The GSE19769 profile was from the platform of GPL9081 including 10 cases of H. pylori-negative and 9 cases of H. pylori-positive gastritis specimen, while the GSE27411 from the platform of GPL6255 containing 6 cases of H. pylori-negative and 6 cases of H. pylori-positive atrophic gastritis.

2.2. Identification of DEMs and DEGs

The significance analysis of DEGs and DEMs in H. pylori-negative and -positive samples was performed using the R language software (version 3.6.0, https://www.r-project.org/) and a limma R package. The Benjamini and Hochberg false discovery rate (FDR) method was used to adjust the value to reduce the false-positive risk. The raw data of miRNAs and mRNAs expression were averaged and normalized, and the data of miRNAs were also log2-transformed, while the data with a median expression value of zero or less than zero were removed. value < 0.05 and |log2 fold change (FC)| > 1 were used as the filter threshold for identifying DEGs and DEMs. A pheatmap R package was used for hierarchical clustering analysis and for drawing heat maps of DEGs and DEMs, while a ggplot2 R package was used for drawing volcano plot [13, 14].

2.3. Interactive Analysis and Construction Expression Network of miRNA-mRNA

MiRNA-mRNA regulatory networks were constructed to predict the potential hub genes of DEMs and DEGs. First, TargetScan (http://www.targetscan.org/vert_72/), miRDB (http://mirdb.org/), PicTar (https://pictar.mdc-berlin.de/), and Miranda (http://miranda.org.uk/) were used to predict targets of DEMs, while different software programs with different algorithms and only genes predicted by at least 3 software programs were selected as the targets of DEMs. The overlapped genes of the targets of DEMs and differential expression genes of GSE27411 dataset were used to construct coexpression networks by Cytoscape software [18].

2.4. GO and KEGG Analysis of Cross-Genes

GO analysis is a common approach for gene annotation and gene classification from three aspects of cellular component (CC), biological process (BP), and molecular function (MF) [15]. KEGG is a comprehensive database resource with 17 main databases, which are used to understand advanced gene functions and practical biological systems [16]. The org.Hs.eg.db, clusterProfiler, enrichplot, and ggplot2 R packages were used for GO and KEGG analysis with a cut-off criterion of a value < 0.05 [19].

2.5. TCGA Data Processing

The TCGA database is a comprehensive and open database, which contains a variety of human cancer types [20], and was used for validating the relations between the hub genes of the network and the GC prognosis. The inclusion criteria were as follows: (1) the primary site is the stomach; (2) the disease type is the adenomas and adenocarcinomas; (3) the data category is transcriptome profiling; (4) gene expression quantification is used for gene data type, while miRNA expression quantification is used for miRNA data type. Eventually, the mature miRNA expression and clinical data of 452 GC cases (42 normal and 410 having tumors) were downloaded from the TCGA database, and the gene expression and clinical data of 374 GC cases (30 normal and 344 having tumors) were obtained. Survival analysis was performed with a survival R package, while the original data were standardized by the log2 (x + 1) method. The prognostic value of cross-genes was determined by Kaplan–Meier analysis, and was considered as a significant difference.

2.6. Patient Data

A total of 431 patients, who were confirmed as gastritis or primary GC from March 2019 to October 2019 at the Affiliated Hospital of Binzhou Medical University (Yantai, China) and the Yuhuangding Hospital (Yantai, China), were enrolled in this study. Blood urease was used to detect the H. pylori infection. However, 362 patients were excluded because of the lacking of accurate classification or detecting of the H. pylori infection. Finally, 69 patients, including 18 negative and 51 positive for H. pylori infection, were included in this study. The related clinical and pathological characteristics are listed in Table 1. The samples involved in this study have been approved by the ethics society of Binzhou Medical University and by the patients themselves or their families. None of the patients has received prior radiotherapy or chemotherapy.


CharacteristicsPatients (n)

Age (year)
 ≤6031
 >6038
Sex
 Male46
 Female23
Histopathological type
 Gastritis32
 Gastric adenocarcinoma37
H. pylori infection
 Negative18
 Positive51

2.7. Detecting DEMs by SYBR Green Real-Time PCR (RT-PCR)

Total RNA was extracted from paraffin-embedded tissue by miRNeasy FFPE Kit (no. 217504, Qiagen, German), and the RNA concentration and purity were measured by NanoDrop® ND-2000, while adjusting the A260/A280 ratio of RNA solution from 1.8 to 2.1. The cDNA synthesis was performed with the miRNA First Strand cDNA Synthesis (tailing reaction) (no. B532451, Sangon, China). The expression of hsa-miR-196b-3p and hsa-miR-196b-5p was detected in the ABI7500 quantitative PCR system (Applied Biosystems, USA) instrument by SYBR Green Premix Ex Taq II (Tli RNase H Plus, Takara, Japan). U6 small nuclear RNA was used as the internal controls. The 20 μl reaction mixture included 10.0 μl SYBR Green Master Mix, 4 μl cDNA template, 0.4 μl ROX Reference Dye, 1.0 μl primer pairs (10 μm), and 4.6 μl deionized water. PCR cycle was performed as follows: initial degeneration for 30 s at 95°C, followed by 40 cycles of 5 s at 95°C and 34 s at 60°C. The relative expression of hsa-miR-196b-3p and hsa-miR-196b-5p was calculated by comparing the cycle threshold (CT) method using the 2−△△ct method with U6 expression according to [21]. The primers of hsa-miR-196b-3p were 5′-TCGACAGCACGACACTGCCTTC-3′ (sense) and 5′- GACACGGACCCACAGACA-3′ (antisense), while the hsa-miR-196b-5p primers were 5′-GCACCAGCGTAGGTAGTTTCC-3′ (sense) and 5′-TATGCTTGTTCTCGTCTCTGTGTC-3′ (antisense).

2.8. Detecting Cross-Genes by Immunohistochemistry

Immunohistochemistry (IHC) was performed in paraffin-embedded sections on the basis of the standardized protocol. Briefly, paraffin-embedded sections (2–4 μm) were deparaffinized with xylene and rehydrated in series gradient ethanol (100%, 95%, and 85% for 1 min, respectively) at room temperature. Heat antigen repair was performed in an autoclave (121°C for 1.5 min) in a citrate sodium buffer (0.01 M), and then endogenous peroxidase was blocked using 3% hydrogen peroxide for 10 min. Sections were incubated with anti-human antibody SMAD6 (ab80049, Abcam), TGFB2 (ab53778, Abcam), PITX2 (ab98297, Abcam), CALML4 (A5086, Zhongshan Golden Bridge, Beijing), and NRP1 (ab25998, Abcam) for 1 h at 37°C and then incubated with a rabbit polyclonal antibody for 20 minutes. After dying with diaminobenzidine for 5–10 min at room temperature, sections were sealed with neutral balata, respectively. The sections were evaluated with semiquantitative method. Briefly, more than 400 cells were counted in each section, while some necrotic cells and peripheral-colored cells were elided. More than 10% of cells were nuclear; staining in all cells was defined as protein-positive expression, while less than 10% was protein-negative expression.

2.9. Statistical Analysis

All data were expressed as mean ± standard deviation (SD) of 3 independent experiments, and statistical analyses were performed with SPSS 17.0 (SPSS Ins., Chicago, IL, USA). The difference of miRNA and protein expression was analyzed with two-tail unpaired t-test with considered as significant difference.

3. Results

3.1. Identification of DEMs and DEGs

The dataset GSE19769 was selected to screen DEMs including 10 H. pylori-negative samples and 9 H. pylori-positive samples. 470 human miRNAs were analyzed in this dataset, and 22 DEMs have met the filtration criteria of logFC > 1 and value < 0.05, including 11 upregulated and 11 downregulated miRNAs (Table 2). Volcano map and heat map for the hierarchical clustering of the DEMs were carried out by a pheatmap R package (Figures 1(b) and 1(d)). 6 H. pylori-negative and 6 H. pylori-positive specimens of GSE27411 were analyzed in this study, while there are 18 samples in the dataset. A total of 18187 human mRNAs were expressed, and 68 genes reached the filtration criteria, among which 56 genes showed upregulated and 12 genes showed downregulated expression (Table S1). Volcano map and heat map for the hierarchical clustering of the DEGs were drawn by the pheatmap R package (Figures 1(a) and 1(c)).


miRNA IDLogFCtB valueAdj. valueExpression

hsa-miR-455−1.29961−13.288914.501581.32E − 106.20E − 08Down
hsa-miR-2232.0157257.9955677.0361152.91E − 075.02E − 05Up
hsa-miR-200a-5p−1.04839−7.298495.7924491.02E − 069.54E − 05Down
hsa-miR-146b1.8464937.0980555.4220751.47E − 060.000115Up
hsa-miR-200a-3p−1.18346−6.461584.2082094.97E − 060.000249Down
hsa-miR-1551.415176.2306433.7538077.83E − 060.000335Up
hsa-miR-411−1.54782−5.841292.9715891.71E − 050.00062Down
hsa-miR-551b−1.30788−5.449082.1644243.85E − 050.000994Down
hsa-miR-142-3p1.676885.4301552.125034.00E − 050.000994Up
hsa-miR-203−1.63669−5.129511.4940257.54E − 050.001476Down
hsa-miR-142-5p1.6779535.0514261.328668.90E − 050.001609Up
hsa-miR-153−1.47449−4.888010.9808640.0001260.002045Down
hsa-miR-204−1.55239−4.24223−0.410650.0005130.006184Down
hsa-miR-196b1.5942663.80117−1.366580.0013550.012736Up
hsa-miR-5091.0697633.765869−1.442880.0014650.013238Up
hsa-miR-326−1.20269−3.53852−1.932390.0024170.019915Down
hsa-miR-146a1.3244563.446357−2.12960.002960.022082Up
hsa-miR-299-5p−1.04861−3.32334−2.391320.0038770.027607Down
hsa-miR-520e1.322473.302749−2.434930.0040550.028448Up
hsa-miR-6501.7098053.294791−2.451770.0041270.028522Up
hsa-miR-1381.1149253.195115−2.66190.005130.033381Up
hsa-miR-368−1.02603−3.19023−2.672140.0051850.033381Down

3.2. Network Construction of the miRNA-mRNA

The miRNAs could bind to the 3′ UTR of their targets, resulting in the posttranscriptional suppression of these genes [22, 23]. The biological targets of DEMs were predicted by 4 different software programs, while gene targets of hsa-miR-455, hsa-miR-411, hsa-miR-551b, hsa-miR-509, and hsa-miR-520e were not predicted in the PicTar software. For network construction, the targets were selected in at least 3 databases. The numbers of targets of hsa-miR-455, hsa-miR-223, hsa-miR-200a-5p, hsa-miR-146b, hsa-miR-200a-3p, hsa-miR-155, hsa-miR-411, hsa-miR-551b, hsa-miR-142-3p, hsa-miR-203, hsa-miR-142-5p, hsa-miR-153, hsa-miR-204, hsa-miR-196b, hsa-miR-509, hsa-miR-326, hsa-miR-146a, hsa-miR-299-5p, hsa-miR-520e, and hsa-miR-138 were 88, 43, 9, 18, 131, 60, 36, 2, 65, 89, 65, 151, 69, 48, 189, 25, 18, 18, 20, and 61, respectively (Table S2 and Figures S1 and S2). Intersected genes were obtained with the target genes and all expressed mRNAs in the dataset GSE27411. And all the intersected genes were imported into the Cytoscape software to conduct the miRNAs-mRNA expression network (Figure 2). A total of 16 miRNAs and 509 genes were involved in the network.

3.3. GO and KEGG Analysis of Cross-Genes in the Network

To explore the biological functions of the cross-genes of the network, we used an enrichplot and a ggplot2 R package to analyze the GO categories and KEGG signal pathways of the cross-genes. In the BP, the cross-genes were concentrated in outflow tract septum morphogenesis, epithelial cell migration, and endothelial cell migration. In the CC, the cross-genes were enriched in transcription factor complex, synaptic membrane, and adherence junction, but RNA polymerase II proximal promoter sequence-specific DNA binding, proximal promoter sequence-specific DNA binding, and enhancer binding in the MF (Table 3 and Figure S3(A)). KEGG analysis demonstrated that the cross-genes were prominently enriched in the mitogen-activated protein kinase (MAPK) signaling pathway, Ras signaling pathway, and TGF-β signaling pathway (Table 4 and Figure S3(B)).


CategoryGene termsCount value

Biological processGO:0003148, outflow tract septum morphogenesis101.24E − 09
GO:0003151, outflow tract morphogenesis153.98E − 09
GO:0001655, urogenital system development306.76E − 09
GO:0003279, cardiac septum development171.22E − 08
GO:0060411, cardiac septum morphogenesis131.60E − 07
GO:0010632, regulation of epithelial cell migration261.74E − 07
GO:0010631, epithelial cell migration291.75E − 07
GO:0061564, axon development331.77E − 07
GO:0043542, endothelial cell migration251.77E − 07
GO:0090132, epithelium migration292.09E − 07

Cellular componentGO:0005667, transcription factor complex234.28E − 07
GO:0097060, synaptic membrane269.71E − 06
GO:0045211, postsynaptic membrane212.70E − 05
GO:0099699, integral component of synaptic membrane143.31E − 05
GO:0005912, adherens junction287.12E − 05
GO:0099240, intrinsic component of synaptic membrane147.62E − 05
GO:0033267, axon part220.000108619
GO:0000790, nuclear chromatin200.000128584
GO:0099055, integral component of postsynaptic membrane110.000162968
GO:0000151, ubiquitin ligase complex160.000198327

Molecular functionGO:0000978, RNA polymerase II proximal promoter sequence-specific DNA binding385.75E − 10
GO:0000987, proximal promoter sequence-specific DNA binding395.92E − 10
GO:0035326, enhancer binding191.24E − 09
GO:0000980, RNA polymerase II distal enhancer sequence-specific DNA binding151.89E − 08
GO:0001158 enhancer sequence-specific DNA binding165.18E − 08
GO:0003682, chromatin binding349.92E − 08
GO:0001228, DNA-binding transcription activator activity, RNA polymerase II-specific292.20E − 06
GO:0004714, transmembrane receptor protein tyrosine kinase activity97.94E − 06
GO:0003712, transcription coregulator activity321.19E − 05
GO:0019199, transmembrane receptor protein kinase activity101.30E − 05


KEGG termsCount value

MAPK signaling pathway253.70E − 06
Axon guidance173.89E − 05
Ras signaling pathway199.06E − 05
TGF-beta signaling pathway110.000134809
cGMP-PKG signaling pathway150.000184067
Human cytomegalovirus infection180.000186779
Adrenergic signaling in cardiomyocytes140.00018806
Hepatocellular carcinoma150.000196653
Rap1 signaling pathway170.000244921
ErbB signaling pathway100.000259955
Growth hormone synthesis, secretion, and action120.000281056

3.4. Survival Analysis Verification in TCGA

Survival analysis demonstrated that hsa-miR-196b-3p (, Figure 3(c)) was significantly related to the prognosis of GC patients, while hsa-miR-196b-5p was not (, Figure 3(f)), and the patients with lower hsa-miR-196b-3p expression had poor outcomes. CALML4, SMAD6, PITX2, and TGFB2 gene expression were also closely correlated with the GC patients’ prognosis. The survival time of GC patients with high expression of CALML4 or SMAD6 was significantly longer than that of patients with low mRNA expression of CALML4 (, Figure 3(a)) or SMAD6 (, Figure 3(d)). In addition, patients with high expression of PITX2 or TGFB2 had a significantly poor prognosis than those with low mRNA expression of PITX2 (, Figure 3(b)) or TGFB2 (, Figure 3(e)).

3.5. Validation of miRNA and Cross-Gene Expression

The expression of hsa-miR-196b-3p and hsa-miR-196b-5p was analyzed from the TCGA database to determine its prognostic value of GC. The results revealed that the expression of hsa-miR-196b-5p (log2FC = 3.269, ) and hsa-miR-196b-3p (log2FC = 4.674216894, ; Figure 4(c) and Table S3) in GC tissue was significantly higher than that in normal tissue. The expression of hsa-miR-196b-3p and hsa-miR-196b-5p was also detected by qPCR in 69 GC and gastritis patient specimens, which included 18 H. pylori negative and 51 positive. The results showed that hsa-miR-196b-3p and hsa-miR-196b-5p were overexpressed in the 69 patients with log2FC values of 2.01665 and 1.8458, respectively (Figure 4(d)). Further analysis showed that the expression level of hsa-miR-196b-3p in H. pylori-positive group was significantly higher than that in the negative group (, Figure 4(e)); however, there was no significant difference in hsa-miR-196b-5p expression between the H. pylori-negative and -positive groups (, Figure 4(f)). The protein expression of CALML4, SMAD6, PITX2, and TGFB2 was detected by immunohistochemistry in 69 specimens. The positive rate of CALML4 in H. pylori-positive samples (84.3%) was significantly higher than that in negative samples (55%, ), while there was no significant difference in SMAD6, PITX2, and TGFB2 between the two groups (, Figures 4(a) and 4(b)).

4. Discussion

Although H. pylori infection is closely related to GC, the pathogenesis of H. pylori-related GC has not been clarified [7, 2426]. Therefore, a comprehensive study of the molecular mechanism of H. pylori-related GC may be helpful to understand the disease and get better diagnosis and treatment methods. As an important part of bioinformatics, gene expression microarray, which has been widely used in tumor research [27, 28], can analyze the expression of thousands of genes. In this study, the miRNA and mRNA expression profile data of H. pylori-infected gastric tissue from GEO database were analyzed, DEMs and DEGs were analyzed with R software, DEMs targets were predicted, miRNA-mRNA expression network was constructed, prognostic value of hub genes for GC was verified, and hub genes expression was detected in the clinical sample. We screened 1 miRNA (hsa-miR-196b-3p) and 4 important nodal genes (CALML4, SMAD6, PITX2, and TGFB2) in H. pylori-related GC which can be used as biomarkers for GC prognosis.

miRNA participates in biological, pathological processes and infections by downregulating target expression [29]. The ectopic expression of miRNAs plays an important role in the pathogenesis of multiple cancers, including H. pylori-related GC [8, 30]. Several studies have shown that hsa-miR-196b is upregulated in GC and related to the outcomes of GC patients [31, 32], while other research showed that hsa-miR-196b expression in Epstein–Barr virus– (EBV-) positive GC was significantly lower than that in EBV negative [33]. However, the role of hsa-miR-196b in GC, especially in EBV or H. pylori infection, is still unknown.

Our study showed that the expression level of hsa-miR-196b in H. pylori-positive group was significantly higher than that in H. pylori-negative group. Analyzing TCGA data showed that the hsa-miR-196b-3p expression, rather than hsa-miR-196b-5p, could be used as a better biomarker for GC prognosis, which was consistent with the previous report [31]. And we also verified the hsa-miR-196b-3p and hsa-miR-196b-5p expression in clinical samples with or without H. pylori infection. The results showed that the expression of hsa-miR-196b-3p in H. pylori-positive group was significantly higher than that in the negative group (), while there was no significant difference in hsa-miR-196b-5p expression (). miRNA target analysis showed that hsa-miR-196b and hsa-miR-326 could regulate the expression of SMAD6 involved in many biological activities through phosphorylation of the TGF-β signaling pathway [34, 35]. Then, the results suggested that H. pylori may be involved in the pathogenesis of GC with hsa-miR-196b by regulating the expression of many genes and activating the infectious immune pathways. However, the molecular mechanism of hsa-miR-196b in infection and GC needs further experimental study.

Studies have shown that the hsa-miR-200 family, including hsa-miR-200a, participates in the negative feedback loop formed by ZEB1, ZEB2, and SIP1, in which hsa-miR-200 suppresses the expression of ZEB1, ZEB2, and, SIP1 and then downregulates their expression [36, 37]. ZEB2 and SIP1 have also been suggested to inhibit the transcription of cyclin D1 [38]. Some studies showed that, in H. pylori infection, the CagA gene can promote transformation from G1 into S/G2 in the cell cycle through activating AP-1 and cAMP, which suggested that hsa-miR-200 was involved in the transformation from gastric epithelial cells to EMT through ZEB loop [39]. Our study showed that hsa-miR-200a-3p showed low expression in H. pylori-related GC, and it could inhibit the expression of PITX2 and TGFB2, which participate in the TGF-β pathway. Analysis of the relationship between cross-genes and GC prognosis showed that PITX2 and TGFB2 were related to the prognosis of GC. The survival time of GC patients with overexpression of PITX2 and TGFB2 was remarkably shorter than those with underexpression, while hsa-miR-200a-3p had no prognostic value for GC.

Some results have shown that hsa-miR-223 was abnormally overexpressed in GC and was significantly upregulated in H. pylori-infected tissues, which could be involved in the pathogenesis of GC by targeting FBXW7 [40]. Some studies showed that hsa-miR-411 showed low expression in GC, and overexpression of hsa-miR-411 in GC cell led to decreasing proliferation and increasing apoptosis [40], while hsa-miR-411 could not be used as an independent predictor of GC prognosis [41]. However, the role of hsa-miR-411 in H. pylori infection has not been reported. Our results showed that the expression level of hsa-miR-223 was high in H. pylori-positive GC, while that of hsa-miR-411 was low, and both of them can regulate CALML4, which can activate the CGMP-PKG signaling pathway. Besides, patients with overexpression of CALML4 had better outcomes than those with underexpression. However, the mechanism of CALML4 regulated by hsa-miR-223 and hsa-miR-411 in the pathogenesis of H. pylori-related GC needs further experimental study.

In conclusion, we constructed a coexpression network of miRNA-mRNA and identified the key genes of hsa-miR-196b, CALML4, PITX2, TGFB2, and SMAD6 in H. pylori-related GC, which may provide a new way for the diagnosis and treatment of H. pylori-related GC.

Abbreviations

H. pylori:Helicobacter pylori
miRNAs:MicroRNAs
DEMs:Differentially expressed miRNAs
DEGs:Differentially expressed genes
GC:Gastric cancer
GEO:Gene expression omnibus
BP:Biological process
CC:Cellular component
MF:Molecular function
qPCR:Quantitative polymerase chain reaction
KEGG:Kyoto Encyclopedia of Genes and Genomes
GO:Gene ontology
TCGA:The Cancer Genome Atlas
EBV:Epstein–Barr virus
FC:Fold change
SD:Standard deviation
FDR:False discovery rate.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Ping Yang, Junjie Liu, and Tianci Yang contributed equally to this work.

Acknowledgments

This study was supported by the Science and Technology Program of Shandong Province (Grant no. J15LK02), Scientific Research Project of Yantai (Grant no. 2016ZH081), Youth Research Initiation of Yantai Yuhuangding Hospital (Grant no. 201801), and Science and Technology Innovation Plan for College Student (Grant no. 201910440022).

Supplementary Materials

Figure S1: Venn diagram of potential targets of DEMs predicted by 4 software programs. Venn for hsa-miR-455 (A), hsa-miR-223 (B), hsa-miR-200a-5p (C), hsa-miR-146b (D), hsa-miR-200a-3p (E), hsa-miR-155 (F), hsa-miR-411 (G), hsa-miR-551b (H), hsa-miR-142-3p (I), hsa-miR-203 (J), hsa-miR-142-5p (K), and hsa-miR-153 (L). Figure S2: Venn diagram of potential targets of DEMs predicted by 4 software programs. Venn for hsa-miR-204 (A), hsa-miR-196b (B), hsa-miR-509 (C), hsa-miR-326 (D), hsa-miR-146a (E), hsa-miR-299-5p (F), hsa-miR-520e (G), and hsa-miR-138 (H). Figure S3: GO and KEGG function analysis of the cross-genes. Circle diagram of (A) GO clusters and (B) KEGG pathway clusters. Table S1: differential expression genes of H. pylori-negative and -positive patients. Table S2: targets of DEMs in the network. Table S3: expression of hsa-miR-196b-3p and hsa-miR-196b-5p in TCGA. (Supplementary Materials)

References

  1. W. Chen, R. Zheng, P. D. Baade et al., “Cancer statistics in China, 2015,” CA: A Cancer Journal for Clinicians, vol. 66, no. 2, pp. 115–132, 2016. View at: Publisher Site | Google Scholar
  2. F. Bray, A. Jemal, N. Grey, J. Ferlay, and D. Forman, “Global cancer transitions according to the Human Development Index (2008-2030): a population-based study,” The Lancet Oncology, vol. 13, no. 8, pp. 790–801, 2012. View at: Publisher Site | Google Scholar
  3. B. Marshall and J. R. Warren, “Unidentified curved bacilli in the stomach of patients with gastritis and peptic ulceration,” The Lancet, vol. 323, no. 8390, pp. 1311–1315, 1984. View at: Publisher Site | Google Scholar
  4. M. Amieva and R. M. Peek, “Pathobiology of Helicobacter pylori-induced gastric cancer,” Gastroenterology, vol. 150, no. 1, pp. 64–78, 2016. View at: Publisher Site | Google Scholar
  5. K. Matsushima, H. Isomoto, N. Inoue et al., “MicroRNA signatures in Helicobacter pylori-infected gastric mucosa,” International Journal of Cancer, vol. 128, no. 2, pp. 361–370, 2011. View at: Publisher Site | Google Scholar
  6. Y. C. Lee and J. T. Lin, “Screening and treating Helicobacter pylori infection for gastric cancer prevention on the population level,” Journal of Gastroenterology and Hepatology, vol. 32, no. 6, pp. 1160–1169, 2017. View at: Publisher Site | Google Scholar
  7. K. Matsushima, H. Isomoto, N. Inoue et al., “MicroRNA signatures in Helicobacter pylori-infected gastric mucosa,” International Journal of Cancer, vol. 128, no. 2, pp. 361–370, 2011. View at: Publisher Site | Google Scholar
  8. J. H. Song and S. J. Meltzer, “MicroRNAs in pathogenesis, diagnosis, and treatment of gastroesophageal cancers,” Gastroenterology, vol. 143, no. 1, pp. 35–47, 2012. View at: Publisher Site | Google Scholar
  9. X. Li, Y. Zhang, Y. Zhang, J. Ding, K. Wu, and D. Fan, “Survival prediction of gastric cancer by a seven-microRNA signature,” Gut, vol. 59, no. 5, pp. 579–585, 2010. View at: Publisher Site | Google Scholar
  10. S. A. Melo and R. Kalluri, “Molecular pathways: microRNAs as cancer therapeutics,” Clinical Cancer Research, vol. 18, no. 16, pp. 4234–4239, 2012. View at: Publisher Site | Google Scholar
  11. K. Wu, C. Zhang, C. Zhang, and D. Dai, “A novel three-miRNA signature identified using bioinformatics predicts survival in esophageal carcinoma,” Biomed Research International, vol. 2020, Article ID 5973082, 2020. View at: Publisher Site | Google Scholar
  12. J. Pallarès-Albanell, M. T. Zomeño-Abellán, G. Escaramís et al., “A high-throughput screening identifies MicroRNA inhibitors that influence neuronal maintenance and/or response to oxidative stress,” Molecular Therapy - Nucleic Acids, vol. 17, pp. 374–387, 2019. View at: Publisher Site | Google Scholar
  13. R. Edgar, M. Domrachev, and A. E. Lash, “Gene expression omnibus: NCBI gene expression and hybridization array data repository,” Nucleic Acids Research, vol. 30, no. 1, pp. 207–210, 2002. View at: Publisher Site | Google Scholar
  14. T. Barrett, S. E. Wilhite, P. Ledoux et al., “NCBI GEO: archive for functional genomics data sets–update,” Nucleic Acids Research, vol. 41, no. D1, pp. D991–D995, 2012. View at: Publisher Site | Google Scholar
  15. M. Ashburner, C. A. Ball, J. A. Blake et al., “Gene ontology: tool for the unification of biology,” Nature Genetics, vol. 25, no. 1, pp. 25–29, 2000. View at: Publisher Site | Google Scholar
  16. M. Kanehisa and S. Goto, “KEGG: Kyoto encyclopedia of genes and genomes,” Nucleic Acids Research, vol. 28, no. 1, pp. 27–30, 2000. View at: Publisher Site | Google Scholar
  17. I. Nookaew, K. Thorell, K. Worah et al., “Transcriptome signatures in Helicobacter pylori-infected mucosa identifies acidic mammalian chitinase loss as a corpus atrophy marker,” BMC Medical Genomics, vol. 6, no. 1, p. 41, 2013. View at: Publisher Site | Google Scholar
  18. P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003. View at: Publisher Site | Google Scholar
  19. C. Zhang, Y. Zheng, X. Li, X. Hu, F. Qi, and J. Luo, “Genome-wide mutation profiling and related risk signature for prognosis of papillary renal cell carcinoma,” Annals of Translational Medicine, vol. 7, no. 18, p. 427, 2019. View at: Publisher Site | Google Scholar
  20. G. A. R. N. Cancer, J. N. Weinstein, E. A. Collisson et al., “The Cancer Genome Atlas Pan-Cancer analysis project,” Nature Genetics, vol. 45, no. 10, pp. 1113–1120, 2013. View at: Publisher Site | Google Scholar
  21. K. J. Livak and T. D. Schmittgen, “Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method,” Methods, vol. 25, no. 4, pp. 402–408, 2001. View at: Publisher Site | Google Scholar
  22. S. Xu, “microRNA expression in the eyes and their significance in relation to functions,” Progress in Retinal and Eye Research, vol. 28, no. 2, pp. 87–116, 2009. View at: Publisher Site | Google Scholar
  23. R. Garzon, G. Marcucci, and C. M. Croce, “Targeting microRNAs in cancer: rationale, strategies and challenges,” Nature Reviews Drug Discovery, vol. 9, no. 10, pp. 775–789, 2010. View at: Publisher Site | Google Scholar
  24. B. Marshall and J. R. Warren, “Unidentified curved bacilli in the stomach of patients with gastritis and peptic ulceration,” The Lancet, vol. 323, no. 8390, pp. 1311–1315, 1984. View at: Publisher Site | Google Scholar
  25. M. Amieva and R. M. Peek, “Pathobiology of Helicobacter pylori-induced gastric cancer,” Gastroenterology, vol. 150, no. 1, pp. 64–78, 2016. View at: Publisher Site | Google Scholar
  26. P. Rawla and A. Barsouk, “Epidemiology of gastric cancer: global trends, risk factors and prevention,” Gastroenterology Review, vol. 14, no. 1, pp. 26–38, 2019. View at: Publisher Site | Google Scholar
  27. H. Hackl, P. Charoentong, F. Finotello, and Z. Trajanoski, “Computational genomics tools for dissecting tumour-immune cell interactions,” Nature Reviews Genetics, vol. 17, no. 8, pp. 441–458, 2016. View at: Publisher Site | Google Scholar
  28. C. P. Bracken, H. S. Scott, and G. J. Goodall, “A network-biology perspective of microRNA function and dysfunction in cancer,” Nature Reviews Genetics, vol. 17, no. 12, pp. 719–732, 2016. View at: Publisher Site | Google Scholar
  29. J. Lu, G. Getz, E. A. Miska et al., “MicroRNA expression profiles classify human cancers,” Nature, vol. 435, no. 7043, pp. 834–838, 2005. View at: Publisher Site | Google Scholar
  30. S. A. Melo and R. Kalluri, “Molecular pathways: microRNAs as cancer therapeutics,” Clinical Cancer Research, vol. 18, no. 16, pp. 4234–4239, 2012. View at: Publisher Site | Google Scholar
  31. Y. Shindo, S. Hazama, Y. Nakamura et al., “miR-196b, miR-378a and miR-486 are predictive biomarkers for the efficacy of vaccine treatment in colorectal cancer,” Oncology Letters, vol. 14, no. 2, pp. 1355–1362, 2017. View at: Publisher Site | Google Scholar
  32. J. Ge, Z. Chen, R. Li, T. Lu, and G. Xiao, “Upregulation of microRNA-196a and microRNA-196b cooperatively correlate with aggressive progression and unfavorable prognosis in patients with colorectal cancer,” Cancer Cell International, vol. 14, no. 1, p. 128, 2014. View at: Publisher Site | Google Scholar
  33. A. L. Treece, D. L. Duncan, W. Tang et al., “Gastric adenocarcinoma microRNA profiles in fixed tissue and in plasma reveal cancer-associated and Epstein-Barr virus-related expression patterns,” Laboratory Investigation, vol. 96, no. 6, pp. 661–671, 2016. View at: Publisher Site | Google Scholar
  34. T. Imamura, M. Takase, A. Nishihara et al., “Smad6 inhibits signalling by the TGF-β superfamily,” Nature, vol. 389, no. 6651, pp. 622–626, 1997. View at: Publisher Site | Google Scholar
  35. M. J. Goumans and C. Mummery, “Functional analysis of the TGFbeta receptor/Smad pathway through gene ablation in mice,” The International Journal of Developmental Biology, vol. 44, no. 3, pp. 253–265, 2000. View at: Google Scholar
  36. U. Burk, J. Schubert, U. Wellner et al., “A reciprocal repression between ZEB1 and members of the miR-200 family promotes EMT and invasion in cancer cells,” EMBO Reports, vol. 9, no. 6, pp. 582–589, 2008. View at: Publisher Site | Google Scholar
  37. P. A. Gregory, A. G. Bert, E. L. Paterson et al., “The miR-200 family and miR-205 regulate epithelial to mesenchymal transition by targeting ZEB1 and SIP1,” Nature Cell Biology, vol. 10, no. 5, pp. 593–601, 2008. View at: Publisher Site | Google Scholar
  38. J. Mejlvang, M. Kriajevska, C. Vandewalle et al., “Direct repression of cyclin D1 by SIP1 attenuates cell cycle progression in cells undergoing an epithelial mesenchymal transition,” Molecular Biology of the Cell, vol. 18, no. 11, pp. 4615–4624, 2007. View at: Publisher Site | Google Scholar
  39. Y.-J. Chang, M.-S. Wu, J.-T. Lin, R. G. Pestell, M. J. Blaser, and C.-C. Chen, “Mechanisms for Helicobacter pylori CagA-induced cyclin D1 expression that affect cell cycle,” Cellular Microbiology, vol. 8, no. 11, pp. 1740–1752, 2006. View at: Publisher Site | Google Scholar
  40. X. Zhou, W. Jin, H. Jia, J. Yan, and G. Zhang, “MiR-223 promotes the cisplatin resistance of human gastric cancer cells via regulating cell cycle by targeting FBXW7,” Journal of Experimental & Clinical Cancer Research: CR, vol. 34, no. 1, p. 28, 2015. View at: Publisher Site | Google Scholar
  41. S. Azarbarzin, M. A. H. Feizi, R. R. Safaralizadeh et al., “Reduced expression of miR-411 in intestinal type of gastric adenocarcinoma,” Meta Gene, vol. 10, pp. 23–26, 2016. View at: Publisher Site | Google Scholar

Copyright © 2020 Ping Yang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views427
Downloads557
Citations

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.