Atrophic gastritis and gastric cancer tissue miRNome analysis reveals hsa-miR-129-1 and hsa-miR-196a as potential early diagnostic biomarkers.

BACKGROUND: Gastric cancer (GC) is one of the most frequently diagnosed tumor globally. In most cases, GC develops in a stepwise manner from chronic gastritis or atrophic gastritis (AG) to cancer. One of the major issues in clinical settings of GC is diagnosis at advanced disease stages resulting in poor prognosis. MicroRNAs (miRNAs) are small noncoding molecules that play an essential role in a variety of fundamental biological processes. However, clinical potential of miRNA profiling in the gastric cancerogenesis, especially in premalignant GC cases, remains unclear. AIM: To evaluate the AG and GC tissue miRNomes and identify specific miRNAs' potential for clinical applications (e.g., non-invasive diagnostics).
METHODS: Study included a total of 125 subjects: Controls (CON), AG, and GC patients. All study subjects were recruited at the Departments of Surgery or Gastroenterology, Hospital of Lithuanian University of Health Sciences and divided into the profiling (n = 60) and validation (n = 65) cohorts. Total RNA isolated from tissue samples was used for preparation of small RNA sequencing libraries and profiled using next-generation sequencing (NGS). Based on NGS data, deregulated miRNAs hsa-miR-129-1-3p and hsa-miR-196a-5p were analyzed in plasma samples of independent cohort consisting of CON, AG, and GC patients. Expression level of hsa-miR-129-1-3p and hsa-miR-196a-5p was determined using the quantitative real-time polymerase chain reaction and 2-ΔΔCt method.
RESULTS: Results of tissue analysis revealed 20 differentially expressed miRNAs in AG group compared to CON group, 129 deregulated miRNAs in GC compared to CON, and 99 altered miRNAs comparing GC and AG groups. Only 2 miRNAs (hsa-miR-129-1-3p and hsa-miR-196a-5p) were identified to be step-wise deregulated in healthy-premalignant-malignant sequence. Area under the curve (AUC)-receiver operating characteristic analysis revealed that expression level of hsa-miR-196a-5p is significant for discrimination of CON vs AG, CON vs GC and AG vs GC and resulted in AUCs: 88.0%, 93.1% and 66.3%, respectively. Compar-ing results in tissue and plasma samples, hsa-miR-129-1-3p was significantly down-regulated in GC compared to AG (P = 0.0021 and P = 0.024, tissue and plasma, respectively). Moreover, analysis revealed that hsa-miR-215-3p/5p and hsa-miR-934 were significantly deregulated in GC based on Helicobacter pylori (H. pylori) infection status [log2 fold change (FC) = -4.52, P-adjusted = 0.02; log2FC = -4.00, P-adjusted = 0.02; log2FC = 6.09, P-adjusted = 0.02, respectively].
CONCLUSION: Comprehensive miRNome study provides evidence for gradual deregulation of hsa-miR-196a-5p and hsa-miR-129-1-3p in gastric carcinogenesis and found hsa-miR-215-3p/5p and hsa-miR-934 to be significantly deregulated in H. pylori carrying GC patients. ©The Author(s) 2022. Published by Baishideng Publishing Group Inc. All rights reserved.

Core Tip: In this research we aimed to evaluate microRNAs profiles of premalignant and malignant stages of gastric cancer (GC). To date this is the first study analyzing atrophic gastritis (AG) and GC tissue miRNomes in the subjects of European origin using next-generation sequencing approach. We showed that hsa-miR-196a-5p expression in tissue is significant for discrimination between controls and AG or GC, while hsa-miR-129-1-3p is potential candidate for non-invasive GC diagnostic. This study provides novel insights into complex GC pathogenesis cascade and might be highly significant for future studies of new AG or GC associated epigenetic markers or even diagnostic targets.

INTRODUCTION

Gastric cancer (GC) is the one of the most common malignancy and the fourth leading cause of cancer-related death worldwide[1]. Studies show, that in most cases GC development is a stepwise process: Chronic gastric mucosa inflammation progresses to atrophic gastritis (AG) or intestinal metaplasia (IM), which eventually may become predisposition to GC. This complex cascade involves many factors: Helicobacter pylori (H. pylori) infection, lifestyle, dietary habits, and genetic or epigenetic alterations, including miRNA expression changes[2,3]. One of the major concerns in diagnostics of GC is poor survival rate and prognosis, while this tumor is usually diagnosed at late stages. Therefore, investigation of the molecular mechanisms that are critical in the complex GC pathological cascade may help to identify novel therapeutic targets and consequently improve the disease prognosis. MicroRNAs (miRNAs) are small (approx 22 nt) non-coding RNA molecules that regulate gene expression by binding to the specific sites within 3’ untranslated regions of target mRNAs[4,5]. MiRNAs play very important role in many physiological and pathological processes as well as tumorigenesis and may function as either tumor-suppressors or as oncogenic miRNAs[6-8]. Studies have reported numerous differentially expressed miRNAs in malignant gastric tissues including members of miR-20, miR-451, miR-148, miR-223 families[9-11]. Despite the previous efforts and conducted miRNA studies in GC, the miRNome characterization of premalignant gastric condition - AG - remains largely unknown.

In this study, we aimed to investigate miRNome profile through GC tumorigenesis cascade including precancerous lesions, such as AG. Also, expression of two miRNAs (hsa-miR-129-1 and hsa-miR-196a) was analyzed in plasma samples of the independent cohort of AG and GC patients. Tissue miRNome analysis results revealed distinct miRNA profiles comparing controls (CON), AG, and GC groups. Also, our study findings show that two miRNAs: Hsa-miR-129-1 and hsa-miR-196a may be a relevant biomarker for GC diagnostics.

MATERIALS AND METHODS

Study population

The study included a total of 125 CON and patients diagnosed with AG and GC, who were divided into the profiling cohort of 60 subjects and validation cohort of 65 subjects. Tissue samples of profiling cohort were collected during the years 2007-2015, while plasma of participants in validation cohort was collected from years 2011-2019 at the Departments of Surgery and Gastroenterology, Hospital of Lithuanian University of Health Sciences (Kaunas, Lithuania). Clinical and phenotypic characteristics of subjects investigated in profiling and validation cohorts are presented in Table 1. H. pylori status was assessed using indirect ELISA to detect serum-specific IgG antigen (Virion/Serion GmbH, Germany). Control group consisted of subjects, who had no signs of atrophy or IM according to Operative Link on Gastritis Assessment (OLGA) staging system (stage 0)[12]. AG group consisted of individuals that had stage I-IV atrophy score in gastric mucosa by OLGA classification. Gastric adenocarcinoma in GC patients was verified by histology and classified according to the American Joint Committee on Cancer TNM Staging Classification and Lauren Classification[13,14]. Adjacent GC (GCaj) samples were biopsy samples obtained from endoscopically healthy appearing gastric mucosa at least 2 cm away from the primary tumor.

Table 1

Demographic characteristics of profiling and validation cohorts

		Profiling cohort (n = 60)			Validation cohort (n = 65)
		CON (n = 21)	AG (n = 19)	GC (n = 20)	CON (n = 11)	AG (n = 30)	GC (n = 24)
Age	Mean ± SD	58.29 ± 15.52	69.21 ± 8.78	64.95 ± 10.89	42.27 ± 12.89	68.01 ± 11.81	68.33 ± 11.27
Gender (n)	Male	5	3	15	5	9	18
	Female	16	16	5	6	21	6
Helicobacter pylori infection (n)	Negative	12	10	8	-	17	9
	Positive	9	9	9	-	10	4
	Unknown	-	-	3	11	3	11
Differentiation grade (n)	G1	-	-	4	-	-	-
	G2	-	-	4	-	-	12
	G3	-	-	12	-	-	12
Lauren classification (n)	Diffuse	-	-	10	-	-	8
	Intestinal	-	-	10	-	-	13
	Mixed	-	-	-	-	-	2
	Unknown	-	-	-	-	-	1
T (n)	T1	-	-	6	-	-	3
	T2	-	-	2	-	-	5
	T3	-	-	8	-	-	9
	T4	-	-	4	-	-	6
	Unknown	-	-	-	-	-	1
N (n)	N0	-	-	10	-	-	6
	N1	-	-	2	-	-	5
	N2	-	-	3	-	-	4
	N3	-	-	5	-	-	8
	Unknown	-	-	-	-	-	1
M (n)	M0	-	-	7	-	-	14
	M1	-	-	2	-	-	9
	Unknown	-	-	11	-	-	1

SD: Standard deviation; CON: Control; AG: Atrophic gastritis; GC: Gastric cancer.

The study was approved by the Kaunas Regional Biomedical Research Ethics Committee (approval No BE-2-10 and BE-2-31) and performed in accordance with the Declaration of Helsinki. All study participants provided written informed consent before enrollment.

Total RNA extraction

Total RNA, including small RNA fraction, was isolated from CON, AG and GC tissues using miRNeasy Mini Kit (Qiagen, Germany) according to the manufacturer’s instructions. Quantification of RNA was performed using Nanodrop2000 spectrophotometer (Thermo Fisher Scientific, United States) and quality of RNA samples was evaluated by Agilent 2100 Bioanalyzer (Agilent Technologies, United States). Circulating nucleic acids, including circulating miRNA fraction, was isolated using QIAamp Circulating Nucleic Acid Kit (Qiagen, Germany) according to manufacturer’s instructions. All isolated samples were stored at -80 °C prior to further analysis.

Small RNA-seq library preparation and next-generation sequencing

Small RNA libraries were prepared using Illumina TruSeq Small RNA Sample Preparation Kit (Illumina, United States) according to the manufacturer’s protocol with 1 μg RNA input per sample followed by RNA 3’ adapter ligation, RNA 5’ adapter ligation, cDNA synthesis, polymerase chain reaction (PCR) amplification using unique barcode sequences for each sample and gel size-selection of small RNA library. The yield and quality of sequencing libraries were assessed using the Agilent 2100 Bioanalyzer (Agilent Technologies, United States). The small RNA libraries were randomized, pooled 24 samples per lane and sequenced using Illumina HiSeq 2500 (1 × 50 bp single-end reads).

Bioinformatics analysis of small RNA-seq data

Analysis of raw small RNA-seq data was performed by nf-core/smrnaseq pipeline v.1.0.0 including Nextflow v.20.07.1[15], Java v.11.0.7, and Docker v.19.03.12. In brief, all steps consisted of read quality control using FastQC v.0.11.9, removing 3’ adapter sequences with TrimGalore! v.0.6.5, mapping to mature and hairpin miRNAs (miRBase v.22.1[16]), and GRCh37 human reference genome with Bowtie v.1.3.0[17]. After alignment and trimming sorted BAM files were used for further analysis with edgeR v.3.32.1[18] and mirtop v.0.4.23. MiRNA quality was assessed and summarized using MultiQC v.1.9[19]. Normalized counts were generated using isomiRs package and differential expression analysis was carried out using the DESeq2 Bioconductor package v.1.26.0[20]. The threshold for significant differential expression was Bonferroni[21] adjusted P-value < 0.05 and absolute value of log2 fold change (FC) |log2FC| > 1.

Validation of miRNA expression in plasma by reverse transcription quantitative real-time PCR

To validate differentially expressed miRNAs in plasma samples, isolated plasma circulating microRNA was reverse transcribed to cDNA using the TaqMan™ MicroRNA Reverse Transcription Kit (Thermo Fisher Scientific, United States). The material was preamplified using the TaqMan PreAmp Master Mix (Applied Biosystems, United States) according to the manufacturer’s protocol. Quantitative real-time PCR (RT-PCR) was performed using the TaqMan MicroRNA Assays: Hsa-miR-129* (Assay ID: 002298), hsa-miR-196a (Assay ID: 241070_mat) on 7500 Fast Real-Time PCR System (Applied Biosystems, United States). All RT-qPCR reactions were run in duplicate in a 20 μL reaction and the relative fold change in miRNA expression was estimated using the 2-ΔΔCt method[22]. Ct values were normalized to the RNU6B (Assay ID: 001093, Thermo Fisher Scientific, United States) endogenous control.

Statistical analysis

Statistical analysis was performed using RStudio software (R v.3.6.3). Shapiro-Wilk normality test was used to test the normal distribution of data. For normally distributed data, statistical significance was assessed by Student's t-test. If the data did not pass normality tests was performed non-parametric Wilcoxon rank-sum test. A P < 0.05 was considered statistically significant. Area under the receiver operating characteristic curve (AUC-ROC) analysis was performed using pROC R package.

RESULTS

Small RNA sequencing reveals distinct miRNomes of healthy, premalignant, and malignant stages of GC

Small RNA sequencing of CON, AG, and paired GC (cancerous and adjacent) tissues in total identified 1037 miRNAs annotated in the miRBase v22.1. Sequencing yielded approx 250 M raw sequencing reads (from 359 K to 16 M reads per sample). After quality control steps 396 low-abundant and non-variable miRNAs and 5 outlying samples were removed resulting in 641 miRNAs and 75 samples which were used for further analysis (Supplementary Figures 1 and 2). The number of deregulated miRNAs corresponded to pathological cascade of GC development. The highest number of deregulated miRNAs were determined when comparing GC and CON groups (129 differentially expressed miRNAs, 82 up-regulated and 47 down-regulated; Supplementary Table 1). Next, 99 differentially expressed miRNAs were identified analyzing GC compared to AG (67 up-regulated and 32 down-regulated; Supplementary Table 2). The lowest number, 20 miRNAs, were found to be deregulated comparing AG and CON (6 up-regulated and 14 down-regulated; Supplementary Table 3). Differential expression results comparing GC vs GCaj, AG vs GCaj, and CON vs GCaj are presented in Supplementary Tables 4, 5 and 6 respectively.

Differential expression results and top five deregulated miRNAs in each case are represented in Figure 1A. Multidimensional scaling analysis of normalized expression values, assessing the similarity structure of miRNomes (Spearman’s correlation distance), revealed 4 clusters, corresponding to the CON, AG, GC cancerous and adjacent tissues (Figure 1B). The AG cluster was intermediate between GC and CON, whereas GCaj was overlapping with AG and CON groups.

Figure 1

Results of microRNA differential expression analysis. A: Differentially expressed gastric tissue microRNAs among different conditions. P-adjusted < 0.05 and |log2 fold change| > 1; B: Multidimensional scaling plot based on normalized data showing a clustering corresponding to control, atrophic gastritis, gastric cancerous and adjacent tissues. The density plots show distributions of the first and second dimensions. CON: Control; AG: Atrophic gastritis; GC: Gastric cancerous; GCaj: Gastric adjacent tissue; MDS: Multidimensional scaling.

Hsa-miR-129-1-3p and hsa-miR-196a-5p may be employed for discrimination of healthy, premalignant, and malignant GC cases

To further study miRNome profiles, altered expression of miRNAs was analyzed in three main comparison groups: AG vs CON, GC vs CON and AG vs GC according to clinical significance. Analyzing uniquely deregulated miRNAs, 40 differentially expressed miRNAs were found when compared GC to CON (25.8% of all deregulated miRNAs), 18 (11.6%) - AG compared to GC, and 6 (3.9%) - AG compared to CON (Figure 2). Most of the deregulated miRNAs (n = 79, 68.7%) were similar between GC vs CON and GC vs AG comparison groups. 12 miRNAs (7.7%) were deregulated in both AG and GC groups when compared to CON. Four miRNAs (2.6%) were similarly deregulated between AG vs CON and AG vs GC groups. Finally, only 2 miRNAs (hsa-miR-129-1-3p and hsa-miR-196a-5p) (1.29%) were identified as deregulated between all comparison groups. AUC-ROC analysis revealed that expression level of hsa-miR-129-1-3p in tissues resulted in AUCs: 68.1%; 86.3%, and 78.1%, CON vs AG, CON vs GC, and AG vs GC, respectively (Figures 3A, 3B and 3C). In addition to this, expression level of hsa-miR-196a-5p could be significant for discrimination of CON vs AG, CON vs GC and AG vs GC and resulted in AUCs: 88.0%, 93.1% and 66.3% (Figures 3D, 3E and 3F).

Figure 2

Venn diagram representing the number of commonly and uniquely differentially expressed microRNAs in three different comparison groups. P-adjusted < 0.05 and |log2 fold change| > 1. CON: Control; AG: Atrophic gastritis; GC: Gastric cancer.

Figure 3

Receiver operating characteristic curves showing prediction performances of expression levels. A-C: Hsa-miR-129-1-3p; D-F: Hsa-miR-196a-5p in tissue samples between different comparison groups: Control vs atrophic gastritis; control vs gastric cancer; and atrophic gastritis vs gastric cancer. AUC: Area under the curve; CON: Control; AG: Atrophic gastritis; GC: Gastric cancer.

Hsa-miR-129-1-3p and hsa-miR-196a-5p expression in the plasma follows the expression pattern of CON, AG, and GC tissues

Differential expression analysis of NGS data in tissue samples revealed that hsa-miR-129-1-3p was significantly down-regulated and hsa-miR-196a-5p was up-regulated in AG and GC tissues compared to CON (P = 0.002 and P = 0.00018; P = 1.2 × 10-5 and P = 3.1 × 10-5, respectively). Moreover, hsa-miR-129-1-3p was significantly down-regulated in the case of AG compared to GC (P = 0.0021) and reflected a stepwise process of a pathology (Figure 4A). Therefore, to identify whether the expression changes of these two miRNAs can be detected noninvasively in the body fluids of the patients, hsa-miR-129-1-3p and hsa-miR-196a-5p were selected for RT-qPCR analysis in plasma samples of independent cohort. The analysis showed similar expression patterns in the case of hsa-miR-129-1-3p, which was significantly down-regulated when comparing AG and GC groups (P = 0.024). There were no other significant findings between the groups (Figure 4B).

Figure 4

Hsa-miR-129-1-3p and hsa-miR-196a-5p expression levels in study comparison groups. A: Atrophic gastritis and gastric cancer tissue samples compared to controls; B: Atrophic gastritis and gastric cancer plasma samples compared to controls. Box plot graphs; boxes correspond to the median value and interquartile range. CON: Control; AG: Atrophic gastritis; GC: Gastric cancerous; GCaj: Gastric adjacent tissue.

Hsa-miR-215-3p/5p and hsa-miR-934 may be associated with H. pylori-induced GC

To investigate role of miRNAs in AG atrophy progression (OLGA classification) and H. pylori-induced GC, differential miRNAs profile analysis in the subgroups of the study was performed. The analysis revealed a minor clustering in AG tissues corresponding to OLGA stages (Supplementary Figures 3A and 3H). H. pylori status in GC tissues (Supplementary Figure 3B). However, no significantly deregulated miRNAs were determined comparing I-II OLGA stages vs III-IV OLGA stages (AG tissue samples). On the other hand, analyzing GC group based on H. pylori infection status [H. pylori (neg.) vs H. pylori (pos.)], three miRNAs were shown to be significantly deregulated: Hsa-miR-215-3p (log2FC = -4.52, P-adjusted = 0.02), hsa-miR-215-5p (log2FC = -4.00, P-adjusted = 0.02), and hsa-miR-934 (log2FC = 6.09, P-adjusted = 0.02).

DISCUSSION

This study represents comprehensive miRNome profiling of premalignant and malignant GC cases by implementing high throughput technologies such as NGS. Although there are several studies reporting profiles of GC tissue miRNAs[23,24], analysis of the association between miRNA expression and AG is very scarce reporting only individual miRNAs[25]. Moreover, based on small RNA-seq findings, two miRNAs were analyzed in subjects’ plasma samples to investigate potential non-invasive markers. To our best knowledge this is the first study analyzing AG and GC tissue miRNomes in the subjects of European origin.

First, our study showed different profiles of deregulated miRNAs between tissue samples of studied groups. In total, 20 differentially expressed miRNAs were identified in AG and 129 - in GC comparing to CON; also 99 deregulated miRNAs - comparing GC and AG groups. MiRNAs such as hsa-miR-3131, hsa-miR-483, hsa-miR-150, hsa-miR-200a-3p, hsa-miR-873-5p were previously reported by the GC profiling studies of Pereira et al[23] and Assumpção et al[24]. Yet, we were able to identify number of novel miRNAs (of which hsa-miR-548ba, hsa-miR-4521, hsa-miR-549a were the most deregulated). There are no data showing the role of these novel miRNAs in inflammatory or tumorous processes of gastric tissue. However, recent studies have shown that hsa-miR-548ba was associated with bladder cancer, hsa-miR-549a with the metastasis of renal cancer, and hsa-miR-4521 with H. pylori infection in esophageal epithelial cells[26-28]. Taking into consideration miRNome of AG, hsa-miR-3591-3p, hsa-miR-122-3p and hsa-miR-122-5p, hsa-miR-451a miRNAs were already reported by Liu et al[29], while the most deregulated miRNAs including hsa-miR-215, hsa-miR-4497, and hsa-miR-1251 were reported for the first time in our study. Previous research showed that hsa-miR-215-5p was deregulated in different lesions of the gastrointes-tinal tract (Barrett’s esophagus, intraepithelial neoplastic lesions, ulcerative colitis)[30-32]. However, hsa-miR-4497 and hsa-miR-452 were not previously associated with AG but were reported to play an important role in GC development[33,34].

Next, we identified hsa-miR-215-3p and hsa-miR-215-5p to be down-regulated while hsa-miR-934 - up-regulated in GC group comparing negative and positive H. pylori infection status. Studies revealed the altered expression of various miRNAs in H. pylori-induced GC tissue samples, including miR-934, miR-146a, miR-375, miR-204[35-37]. Although, hsa-miR-215 deregulation was previously associated with GC[38-40], there is no data showing its link with H. pylori infection.

In addition to this, we showed that two miRNAs (hsa-miR-129-1-3p and hsa-miR-196a-5p) were gradually deregulated comparing all three study groups (CON, AG, and GC) which also corresponds to pathological cascade of GC. In concordance to our results, it has already been shown that hsa-miR-129-1-3p was down-regulated in GC tissues, function as a tumor suppressor in GC and even corresponds to the same expression pattern in gastric juice[41,42]. There is no data regarding the hsa-miR-196a expression in AG tissue, however, investigators have revealed that hsa-miR-196a is overexpressed in GC tissue, plasma, commercial cell lines and promotes cell proliferation[43,44]. ROC-AUC analysis suggests great potential of hsa-miR-196a-5p expression in tissue for discrimination of AG and GC in contrast to CON (AUC = 89.5% and AUC = 89.5%, respectively). Therefore, further studies are needed to confirm this finding.

Finally, selected miRNAs were analyzed in independent cohort of CON, AG, and GC plasma samples by using RT-qPCR. Results showed similar deregulation direction in plasma samples as in the tissue samples. However, significant differences were only determined comparing the expression of hsa-miR-129-1-3p between AG and GC suggesting its potential role in non-invasive diagnostics of malignant cases. No significant expression changes were observed between study groups and hsa-miR-196a-5p. Other studies have shown controversial results: Tsai et al[45] reported that miR-196a/b was up-regulated in both the plasma and tissue of metastatic GC patients, while miRNome profiling study revealed that miR-196a-5p was found to be down-regulated in plasma of patients with precursor lesions of GC compared to non-active gastritis[46].

In our study, using NGS and RT-qPCR techniques we have shown the distinct miRNome profiles of CON, AG, GC, GCaj tissues, and potential of specific miRNAs as non-invasive biomarkers. In addition to this, novel miRNAs not previously reported as AG or GC associated epigenetic markers were identified. We have shown that hsa-miR-196a-5p expression in tissue could be significant for discrimination between CON and AG or GC, confirmed hsa-miR-129-1-3p as non-invasive biomarker in disease progression monitoring, and showed that miRNAs could be a great candidate for future research of new diagnostic approaches.

CONCLUSION

In conclusion, we showed gradual deregulation of hsa-miR-196a-5p and hsa-miR-129-1-3p in the gastric carcinogenesis pathway and confirmed hsa-miR-129-1-3p as a possible non-invasive biomarker. We also found hsa-miR-215-3p/5p and hsa-miR-934 to be significantly deregulated in GC based on H. pylori infection status. These data provide novel insights into complex GC pathogenesis cascade which could be highly significant for future studies of new diagnostic GC targets.

ARTICLE HIGHLIGHTS

Research background

Gastric cancer (GC) is a complex disease arising from the interaction of environmental (e.g., diet, smoking, etc.) and host-associated factors [e.g., Helicobacter pylori (H. pylori) infection, genetics, etc.]. Due to its silent course, it is also one of the most lethal cancers worldwide as it is usually diagnosed at the advanced stages.

Research motivation

Novel biomarkers that would help to improve GC patients’ diagnosis and prognosis are highly needed. Studies show that microRNAs (miRNAs) play an important role in many cancers and could be a promising biomarker or even therapeutic target.

Research objectives

The objectives of the study were to analyze whole miRNome profiles of control, premalignant and malignant gastric tissues, and select the potential miRNA markers that could have a potential for minimally invasive GC diagnostics.

Research methods

Total RNA from gastric tissue samples was subjected for small RNA sequencing (smRNA-seq). Plasma total circulating nucleic acids were used for the expression analysis of the most tissue deregulated miRNAs by real-time quantitative polymerase chain reaction. Statistical analysis involved the differential expression and discrimination analyses.

Research results

The abundance of altered expression miRNAs corresponded to a pathological cascade of GC development. Hsa-miR-129-1-3p and has-miR-196a-5p were shown to be deregulated in healthy-premalignant-malignant sequence. In addition to this, we showed that down-regulation of hsa-miR-129-1-3p could also be detected non-invasively in GC patients’ plasma samples. Finally, results indicated that hsa-miR-215-3p/5p and hsa-miR-934 were significantly deregulated based on H. pylori infection status for GC patients.

Research conclusions

Gastric tissue miRNome study provides extensive profiling of control, premalignant and malignant cases. Based on smRNA-seq results several miRNAs were shown as potential gastric carcinogenesis (hsa-miR-196a-5p and hsa-miR-129-1-3p); and H. Pylori-related (hsa-miR-215-3p/5p and hsa-miR-934) biomarkers.

Research perspectives

This study provides novel insights into complex GC pathogenesis cascade and could serve as a reference for future research to support our findings.