Identification of miRNAs differentially expressed during gastric malignant transformation process
To identify differentially expressed miRNAs through the progression from precancerous lesions to GC we analyzed the miRNomes available in ENA and TCGA databases. Fifteen miRNAs were differentially expressed through the gastric malignant transformation process (Table 1). Eight miRNAs (miR-873, miR-141-3p, miR-429-3p, miR-204-5p, miR-200a-3p, miR-592-5P, miR-146a-5p and miR-363-3p) were up-regulated in non-active chronic gastritis compared to normal gastric tissue; however, the expression of these miRNAs significantly decreased as the gastric carcinogenesis model progresses (Fig. 1a). In contrast, seven miRNAs (miR-122-5p, miR-196b-5p, miR-20b-5p, miR378a-3p, miR-99a-5p, miR-194-5p ad miR-101-2-3p) were down-regulated in non-active chronic gastritis compared to normal tissue, but as the degree of lesion progresses, the expression of this set of miRNAs increased significantly (Fig. 2a). These results show a set of miRNAs whose differential expression was associated with the progression of preneoplastic lesions towards tumor establishment. These miRNAs group could constitute biomarkers with diagnostic potential.
Differentially expressed miRNAs during gastric malignant transformation are associated to clinical outcome of GC patients.
To investigate the clinical significance of miRNAs on the outcome of GC patients, we evaluated OS of patients according to the expression levels of each miRNA by Kaplan-Meier assays. Results showed that low expression of miR-429-3p, miR-592-5p, miR-200a-3p, miR-146a-5p, and miR-141-3p whose expression decreases along the lesion grades was associated to poor OS in GC patients, while low levels of miR-363-3p was associated to improved OS (Fig. 1b). On the other hand, high levels of those miRNAs whose expression increased along the malignant progression (miR-194-5p, miR-196b-5p, and miR-378-3p) were associated with better OS (Fig. 2b). Contrariwise, higher expression of miR-99a-5p, miR20b-5p, and miR-122-5p were associated with poor OS of GC patients (Fig. 2b). miR-101-2-3p was not associated to OS in CG patients. Collectively, data showed that the expression levels of miRNAs could act as a predictive biomarker of the long-term outcome of GC patients. To evaluate the diagnostic power of miRNAs, ROC curves were constructed. Result indicated that combined AUC for the miRNAs significantly associated to OS was AUC = 0.68 to progressively downregulated miRNAs (Fig. 1c), while to those miRNAs progressively upregulated was AUC = 0.75 (Fig. 2c), suggesting that the combined expression of miRNAs might useful as potential tool for diagnosis of GC patients.
miRNAs differentially expressed during gastric malignant transformation regulate biological processes and signaling pathways associated to cancer
Gene ontology and biological pathways analysis showed that down-regulated miRNAs during progression of gastric lesions were significantly associated with biological pathways with relevance in cancer as adherent junction and metabolic processes including lysine degradation or fatty acid elongation, among others (Fig. S1a). Interestingly, several oncogenes as HDGF, TAZ, ZEB1/2, STAT4, KEAP1, Gli1, FSCN1, BCL2, ZEB1, L1CAM, MMP-9, ROCK1 and FBW7 have been reported as targets of this set of miRNAs in GC (Fig. S1b). On the other hand, up-regulated miRNAs were significantly associated to Hippo, p53 and TGF- β signaling pathways, ECM-receptor interaction and metabolic process related to proteoglycans in cancer as well as fatty acid biosynthesis (Fig. S1c). Reported target genes of these group of miRNAs include c-Myc, BCL2, PTEN, MAPK1, IGF1R, SOCS2 and FOXA1 (Fig. S1d). Taken together, these results suggest that progressive deregulation of these group of miRNAs could produce a deregulation of oncogenes or tumor suppressor genes associated to GC development. Remarkably, both up- and down-regulated miRNAs and their target genes have not been previously studied in preneoplastic lesions, neither their regulation by pathogenic strains of H. pylori plays through the carcinogenic process.
miRNAs expression profiles identify an overexpression of miR-18a-5p in gastric cancer
To determine whether miRNAs associated to GC multistep might also distinguish between GC and normal gastric tissue, we evaluate miRNAs expression profiles in GC tissues compared to normal tissues. We identified 159 down-regulated and 339 up-regulated miRNAs (Fig. 3a, Table S1). Among deregulated miRNAs, we detected the up-regulation of miR-196b, miR-200a, miR-141, miR-429 and miR-592, while miR-204 was down-regulated (Fig. 3b). The ten most deregulated miRNAs in GC were miR-139 (down-regulated) and miR-21, miR-196a-1, miR-196a-2, miR-196b, miR-135b, miR-146b, miR-501, and miR-18a (up-regulation) (Fig. 3c). Of these, up-regulation of miR-18a has previously been reported in GC patients positive to H. pylori-infection [24]; however, miR-18a has not been evaluated in gastric premalignant lesions. Hence, we compared the expression of miR-18a-5p in normal gastric tissues and GC tissues finding that miR-18a was overexpressed in GC (Fig. 3d). Moreover, higher expression of miR-18a was associated to better OS in GC patients (Fig. 3e) and ROC curves analysis showed AUC = 0.89 (Fig. 3f) demonstrated that the expression of miR-18a might useful as tool for diagnosis and prognosis of GC patients. These results indicate a potential role of miR-18a in gastric carcinogenic process, hence we aimed to evaluate miR-18a in gastric lesions.
miR-18a-5p is down-regulated in preneoplastic lesions infected with H. pylori
The expression of miR-18a was evaluated in premalignant gastric lesions collected from gastric antrum and body from 37 patients with gastric disease by endoscopic analysis (Fig. 4a-c). Seventy-four biopsies were examined (37 samples of gastric antrum and 37 samples gastric body). All samples were evaluated by hematoxylin and eosin (Fig. 4d-i) for histopathology diagnosis and to detect the presence of Helicobacter-like microorganisms. Genotyping of H. pylori was carried out in all samples (Table 2) by evaluating the presence of VacA (s1 and s2 alleles) gene, and cagPAI genes (cagA, cagL, virB10 and virD4), which distinguish pathogenic or non-pathogenic strains of H. pylori (Fig. S2a). The identification of cagPAI genes and the VacA s1 allele indicates an infection by pathogenic H. pylori strain (Fig. S2b), while the expression of VacA s2 allele and the absence of the expression of cagPAI genes indicated the infection by non-pathogenic H. pylori strain (Fig. S2 c, d). Samples with negative gene expression of bacterial genes were considered as negative to the infection by H. pylori.
Histopathological and genotyping examination of biopsies showed 82.4% of non-active gastritis (71.6% H. pylori positive), 12.2% follicular gastritis (77.7% H. pylori positive), 17.6% active gastritis (100 % H. pylori positive), 27.0% atrophic gastritis (85.0% H. pylori positive) and 5.4% intestinal metaplasia (75% H. pylori positive). Clinicopathologic characteristic of patients are shown in table 2. All samples positive to H. pylori were positive for pathogenic strains, except one sample from gastric body with diagnosis of non-active gastritis, that was positive to a non-pathogenic strain.
When evaluating the relative expression of miR-18a-5p in the different gastric lesions a down-regulation of miR-18a-5p was evidenced during malignancy progression (Fig. 5a). However, in metaplasia samples we observed an up-regulation of miR-18a. These results suggest a loss of expression of miR-18a at the initial premalignant gastric lesion, nevertheless, in advanced malignancy lesions the miR-18a expression might be reestablished to contribute with carcinogenesis. Afterwards, we evaluated the impact of the H. pylori infection on the expression of miR-18a hence we compared its expression among H. pylori-positive and -negative premalignant lesions (Fig. 5b). We observed diminished expression of miR-18a in non-active and follicular gastritis samples positive to H. pylori, however in atrophic gastritis the expression of miR-18b was unchanged. Interestingly, in H. pylori positive metaplasia, the expression of miR-18a was lower than those samples negative for bacterial presence. Moreover, miR-18a was down-regulated in infected tissues compared with non-active gastric tissues uninfected (Fig. 5c) however, the change of expression in infected tissues was unclear (Fig. 5d). These results suggested that miR-18a may differentially regulated in absence or presence of H. pylori.
H. pylori is up-regulating miR-18a in infected AGS cells
To explore the role of H. pylori on the expression of miR-18a-5p we performed co-cultures of H. pylori-AGS cells (Fig. 6a). Infection of AGS cells were evaluated at 6 and 12 hours after co-culturing by quantifying VapD gene expression. We demonstrate that VapD gene was overexpressed in co-cultures indicating that infection was established (Fig. 6b). Surprisingly, the expression of miR-18a-5p was significantly elevated in AGS cells at 12 hours of infection (Fig. 6c). These results suggesting that H. pylori favor the expression of miR-18a in cancer cells.