Skip to main content

A multi-omics study to investigate the progression of the Correa pathway in gastric mucosa in the context of cirrhosis



Patients with liver cirrhosis (LC) are prone to gastric mucosa damage. We investigated the alterations of gastric mucosa in LC patients and their possible mechanisms through multi-omics.


We observed significant gastric mucosa microbial dysbiosis in LC subjects. Gastric mucosal microbiomes of LC patients contained a higher relative abundance of Streptococcus, Neisseria, Prevotella, Veillonella, and Porphyromonas, as well as a decreased abundance in Helicobacter and Achromobacter, than control subjects. The LC patients had higher levels of bile acids (BAs) and long-chain acylcarnitines (long-chain ACs) in serum. The gastric mucosal microbiomes were associated with serum levels of BAs and long-chain ACs. Transcriptome analyses of gastric mucosa revealed an upregulation of endothelial cell specific molecule 1, serpin family E member 1, mucin 2, caudal type homeobox 2, retinol binding protein 2, and defensin alpha 5 in LC group. Besides, the bile secretion signaling pathway was significantly upregulated in the LC group.


The alterations in the gastric mucosal microbiome and transcriptome of LC patients were identified. The impaired energy metabolism in gastric mucosal cells and bile acids might aggravate the inflammation of gastric mucosa and even exacerbate the Correa’s cascade process. The gastric mucosal cells might reduce bile acid toxicity by bile acid efflux and detoxification.

Trial registration: ChiCTR2100051070.


Liver cirrhosis (LC) is a common chronic progressive liver disease, the 11th most common cause of death worldwide. About 1 million deaths annually are attributable to liver cirrhosis [1]. The liver interacts directly with the gut through the hepatic portal and bile secretion systems [2]. Increasing studies indicate that patients with cirrhosis have altered gut microbiota, which is one of the significant factors in promoting the progression of liver cirrhosis [3].

Changes in microbiota with cirrhosis have been demonstrated in the stool, intestinal mucosa, serum, ascites, and saliva samples, pointing toward a global mucosal immune impairment in patients with liver cirrhosis [4]. Among published studies, some results indicated that patients with liver cirrhosis and portal hypertension had mucosal abnormalities. They include mucosal inflammatory-like abnormalities (edema, erythema, granularity, and friability), erosive gastritis, gastric ulcer, and intestinal metaplasia [5,6,7]. Enrichment of certain pathogenic bacteria in the stomach can promote further progression of gastric inflammation or gastric cancer [8]. However, the characterization of gastric microbiota in cirrhosis and their relationship with gastric mucosal abnormalities is unclear.

In addition, patients with cirrhosis also have alterations in serum metabolite [9], which could collaborate with mucosal microbiota to irritate lesions. Technological advances in liquid chromatography coupled with high-resolution mass spectrometry are beginning to revolutionize our understanding of the cause of liver cirrhosis. Recent advances indicate that a general metabolic alteration is a hallmark of the intricate road to decompensated cirrhosis [10]. This metabolic alteration is characterized by decreased β-oxidation of fatty acids in the mitochondria and a concomitant increase in extramitochondrial glucose utilization via glycolysis [10, 11]. Nevertheless, the interaction between gastric flora and blood metabolism in patients with liver cirrhosis is still unknown. It is also unclear whether gastric mucosal dysbiosis in cirrhosis patients has any underlying relationship with serum metabolite.

We hypothesized that there are significant alterations in the gastric mucosal microbiota in cirrhosis, like the gut, which may exacerbate mucosal abnormalities in the stomach. We also postulated that perturbations of the gastric mucosal microbiota have significant correlations with serum metabolite in patients with liver cirrhosis, and poses significant effects on the development of the disease. To answer these questions, we performed 16S rRNA high-throughput sequencing to characterize the gastric mucosal microbial communities and liquid chromatography–mass spectrometry (LC–MS) to assess serum metabolism in the same group of patients with liver cirrhosis.


Subjects and sample collection

We performed an observational study of 64 patients with liver cirrhosis and 59 non-liver cirrhosis patients. Patients were recruited at the Qilu Hospital of Shandong University from September 2021 to April 2022. According to international guidelines, LC is diagnosed based on clinical symptoms, physical signs, radiologic examination, laboratory tests, medical history, and cirrhosis-associated complications in chronic liver disease [12]. Participants aged above 18 years who signed informed consent were included into the study. The exclusion criteria were as follows: (1) cardiovascular disease, diabetes, inflammatory bowel disease, and mental disease, (2) individuals who received proton pump inhibitors, antibiotics, hormones, immunosuppressants, or chemotherapy drugs within one month of enrollment, (3) patients with positive rapid urease test. Detailed demographic, clinical, and etiology of cirrhosis model for end-stage liver disease scores were collected for all patients at the time of inclusion (Tables 1, 2). The study was registered in Chinese Clinical (, ChiCTR2100051070). The study was approved by the Ethics Committee of Qilu Hospital of Shandong University and complied with the Declaration of Helsinki.

Table 1 Demographic characteristics of the groups
Table 2 Etiology of cirrhotic patients in each group

We performed high-definition gastroscopy (Pentax EG29-i10, Pentax, Tokyo, Japan) for all enrolled participants. Each piece of mucosa from the gastric antrum to body of the stomach was examined by gastroscopy, and a rapid urease test was carried out. Two pieces of gastric tissue fragments with negative rapid urease test were collected using 2.2 mm sterile biopsy forceps, and all were stored immediately at − 80 °C for the subsequent procedure. One piece was used for 16S rRNA sequencing, and the other was stored in RNALater according to the manufacturer’s instructions and stored at − 80 ℃ for subsequent transcriptomics analysis. All examinations were performed by two experienced gastroenterologists with at least five years’ experience (YYL, YJG).

The blood samples were drawn by experienced nurses after 8h fasting on the 2nd day morning after admission, and were collected in 10 mL EDTA tubes. Samples were centrifuged at 13,000g for 5 min at 4 ℃ within one hour after collection, and the plasma was frozen at − 80 °C freezer. Standard laboratory parameters were evaluated by the central laboratory of the Qilu Hospital of Shandong University (Additional file 1: Table S1). Blood samples, and gastric mucosa were transported under dry ice to laboratory, where they were thawed at the time of analysis. The flow chart of the study population is shown in Additional file 2: Figure S1A.

Microbiome analysis

We performed microbiota analysis using alpha and beta diversity analyses. Samples with < 1% Helicobacter pylori relative abundance were grouped as H. pylori-negative, while samples with > 1% H. pylori relative abundance were grouped as H. pylori-positive [13, 14]. Differentially abundant bacterial taxa are identified by the linear discriminant analysis (LDA) effect size (LEfSe) method. We used microbial network analysis to identify clusters of microbial taxa that were highly correlated (correlation coefficients < -0.5 or > 0.5, q < 0.05).

Metabolome analysis

Unsupervised principal component analysis and partial least squares discriminate analysis (PLS-DA) was used to assess the global metabolic alterations between groups. The partial least squares discriminate analysis model was used with the first principal component of variable importance in projection values (VIP > 1.0) combined with Student’s t-test (P < 0.001) to determine the significantly different metabolites between the LC group and control group. The differentially accumulated metabolites were mapped to the Kyoto Encyclopedia of Genes and Genomes database (KEGG) for descriptive annotation. Spearman method was used to analyze the correlation coefficients for the metabolome and microbiome data integration.

Transcriptome analysis

We identified genes with false discovery rate (FDR) < 0.05 and |log2FC|> 1 in a comparison as significant differential expression gene (DEG). DEGs were considered significantly enriched in a KEGG pathway at q ≤ 0.05 compared with the whole transcriptome background.

Statistical analysis

Statistical analysis was performed using the Mann–Whitney U test and one-way ANOVA when appropriate. Statistical significance was taken as P < 0.05. Data were analyzed using SPSS software version 25.0. All authors had access to the study data and had reviewed and approved the final manuscript.


Increased alpha-diversity and altered overall microbial composition in LC

Total of 69 subjects, including 30 control subjects (S-C) and 39 patients with cirrhosis (S-LC), with similar demographics, were included in this study and the information were shown in Tables 1 and 2. We compared the alpha diversity of gastric microbiota between the S-LC and S-C group, and the Sobs (P < 0.05), Shannon (P < 0.001), Ace (P < 0.05), and Chao1 (P < 0.05) values were significantly higher in S-LC (Additional file 1: Table S2, Fig. 1A). Meanwhile, beta diversity analysis showed separate clusters for S-LC and S-C (P = 0.001, Fig. 1B). The gastric microbiota was dominated by eight phyla: Proteobacteria, Campilobacterota, Firmicutes, Bacteroidetes, Actinobacteria, unclassified_k__norank_d__Bacteria, Fusobacteria, and Cyanobacteria (Fig. 1C), although the two groups presented in different order of relative abundance at the phylum level. The gastric microbiota in liver cirrhosis had an over-representation of Firmicutes, Bacteroidetes, Actinobacteria, and Fusobacteriota (P < 0.001; Additional file 3: Figure S2). At the genus level, several genera, including Streptococcus, Prevotella, Neisseria, Fusobacterium, Haemophilus, Veillonella, Porphyromonas, Actinomyces, Gemella, Alloprevotella, Rothia, Granulicatella, and Peptostreptococcus, significantly increased in relative abundance between the S-C and S-LC (Fig. 1D). The Helicobacter and Achromobacter decreased in LC (Fig. 1D).

Fig. 1
figure 1

Gastric mucosal microbiome dysbiosis between S-C and S-LC. S-C represents gastric mucosa derived from the control group, and S-LC represents gastric mucosa derived from the liver cirrhosis group. In A, B, DF, red is the S-C group, and blue is the S-LC group, *p < 0.05, **p < 0.01, ***p < 0.001. A Increased microbial richness, estimated by Shannon index. B Gastric mucosal microbiota showed relative clustering between control subjects compared with all patients with cirrhosis. C Relative abundance of microbial species at the phylum level. D The difference in species composition at the genus level. E The LDA value distribution histogram shows the species between S-C and S-LC at the genus level (LDA > 4). F The relative abundance of H. pylori was higher in S-C than in S-LC. Statistical significance was determined by the Mann–Whitney U test

Bacteria differentially abundant in LC versus controls

We performed linear discriminant analysis effect size at the genus level to further identify the gastric-specific species signatures. Seven bacterial taxa showed distinct relative abundances between the two groups. Increased abundance in bacteria, including Streptococcus, Neisseria, Prevotella, Veillonella, and Porphyromonas, and decreased abundance in Helicobacter and Achromobacter were observed in S-LC (LDA score > 4, P < 0.05; Fig. 1E). Furthermore, we found that the H. pylori infection was significantly lower in S-LC than in S-C (Fig. 1F).

To understand the potential interplay among differentially abundant bacteria in S-LC and control, we performed the network topology analysis (Additional file 4: Figure S3) at the genus level. In S-C group, co-occurrence interactions were observed, reflecting synergistic interkingdom interactions’ contribution to gastric microbiota homeostasis (Additional file 4: Figure S3A). Much fewer co-occurrence interactions were observed in S-LC (Additional file 4: Figure S3B). In addition, H. pylori had a co-exclusive association with other gastric microbes in S-C predominantly. However, these correlations were relatively infrequent in S-LC. Collectively, the above microbial analysis indicated a state of dysbiosis in the mucosal microbiome of LC patients.

Blood metabolism changes with LC

We recruited 22 cirrhotic patients (age 53.04 ± 11.08 years, 15 men) and 44 control patients (age 48.32 ± 14.37 years, 25 men) who agreed to give blood serum (Tables 1, 2). We used liquid chromatography-mass spectrometry; to analyze the serum samples, and the abundance profiles were obtained for 1540 annotated serum metabolites. It was found that 492 of 1540 metabolites had significantly different abundances (Fig. 2A). We subsequently performed the partial least squares discriminate analysis and the results revealed visual separation between these groups without overfitting (Fig. 2B). In addition, PLSDA-VIP table (Fig. 2C, Additional file 1: Table S3) showed top 30 metabolites with VIP > 1.0 and P < 0.001. Of them, 9 were increased in the control group, and 21 were increased in the LC group. Notably, the VIP value of taurochenodeoxycholate-3-sulfate (VIP = 3.72), taurodeoxycholic acid (VIP = 2.93), and cis-5-tetradecenoylcarnitine (VIP = 2.89) were more than 2, indicating their significant contribution to the disease. Furthermore, we carried out an enrichment analysis for all differential metabolites using the KEGG. The results indicated that the metabolites in LC were mainly associated with sphingolipid metabolism, glycerophospholipid metabolism, cysteine and methionine metabolism (Additional file 5: Figure S4).

Fig. 2
figure 2

The liver cirrhosis group is associated with altered serum metabolites. A Differential metabolites by volcano diagram. B PLS-DA of the metabolites across the two groups. C VIP scores with the corresponding expression heatmap. On the left side is the metabolite heatmap. On the right side is the metabolite VIP bar graph. The bar length indicates the contribution of the metabolite to the difference between the two groups. The higher value means the metabolite is more difference between the two groups. The bar color indicates the P value of the metabolite between the two groups, *p < 0.05, **p < 0.01, and ***p < 0.001

Association of microbes and metabolites with LC

To further determine the relationships between gastric microbiota and metabolic changes, we subsequently perform the Spearman’s correlation analysis of serum differential metabolite and microbiota in stomach. Interestingly, there was a strong correlation between a large number of stomach microbiota and altered metabolites (Fig. 3A). The correlation was considered statistically significant with correlation coefficient |r|> 0.7 and p < 0.05 in this study, and results were visualized as heatmaps.

Fig. 3
figure 3

Correlation analysis between the microbiota and metabolites. A Heatmap of gastric mucosal flora and serum metabolite correlation. Asterisk(***) denotes a P value less than 0.05 and |r|> 0.7. B Heatmap of gastric mucosal flora and serum long-chain ACs correlation. Asterisk(***) denotes a P value less than 0.05 and |r|> 0.6. C Heatmap of gastric mucosal flora and serum bile acids correlation. Asterisk(***) denotes a P value less than 0.05 and |r|> 0.7

It was found that genus Streptococcus, enhanced in the gastric microbiota community in LC patients, was positively correlated with taurodeoxycholic acid (r = 0.71) and (E)-2-Penten-1-ol (r = 0.82). In addition, His Phe Trp and phosphatidylserine were negatively correlated with almost all the bacteria except Helicobacter (r = 0.65; r = 0.69). The relative concentrations of some long-chain acylcarnitines (palmitoyl-l-carnitine, octadecenoylcarnitine, 3-hydroxyhexadecanoyl carnitine, gamma-linolenic carnitine, and cis-5-tetradecenoylcarnitine) were generally seen to increase in the serum metabolites from cirrhotic patients. As shown in Fig. 3B, increased long-chain ACs were positively associated with almost all the bacteria, such as Actinomyces and cis-5-tetradecenoylcarnitine (r = 0.66). Besides, taurodeoxycholic acid, deoxycholic acid glycine conjugate, glycocholic acid, taurocholic acid, cholic acid, taurochenodeoxycholate-3-sulfate, taurochenodeoxycholate-7-sulfate, chenodeoxycholic acid 3-sulfate, and their derivatives were positively correlated with almost all the bacteria, except 7a,12a-dihydroxy-3-oxo-4-cholenoic acid (Fig. 3C).

Gastric mucosal transcriptome changes in LC

To obtain a comprehensive view of gastric mucosa influenced by microbial colonization and metabolic alterations, we further investigated the gastric mucosal transcriptome in two groups comprising 10 control and 10 LC patients. It was found that 181 and 124 genes were differentially down and upregulated, respectively. Of them, the endothelial cell specific molecule 1, serpin family E member 1, mucin 2, caudal type homeobox 2, and retinol binding protein 2 were significantly upregulated genes (Fig. 4A, B) and were associated with intestinal metastasis and gastric carcinogenesis [15,16,17,18]. We also found that defensin alpha 5 (DEFA5), a gene encoding an antimicrobial peptide, is significantly upregulated (P < 0.001). Then, we mapped all DEGs to KEGG pathways, and the top 11 specific pathways were represented in a bubble chart (Fig. 4C). Furthermore, genes associated with the most significantly enriched pathways (q < 0.05) were shown in Additional file 6: Figure S5. Among these 11 pathways, “neuroactive ligand-receptor interaction” was the most represented pathway. Significantly, it was found that the bile secretion signaling pathway was significantly upregulated in the LC group. In the bile secretion signaling pathway, solute carrier family 51 subunit alpha (SLC51A), solute carrier family 51 subunit beta(SLC51B), cytochrome P450 3A4(CYP3A4) were up-regulated.

Fig. 4
figure 4

S-C represents gastric mucosa derived from the control group, and S-LC represents gastric mucosa derived from the treatment group. A Volcano map of DEGs. B Venn graph of DEGs. C DEGs enriched in the KEGG pathway. The X-axis represents the rich factor, indicating the ratio of enriched genes to total genes in this pathway. A more prominent rich factor indicates more significant enrichment. ESM1 endothelial cell specific molecule 1, SERPINE1 serpin family E member 1, MUC2 mucin 2, CDX2 caudal type homeobox 2, RBP2 retinol binding protein 2, DEFA5 defensin alpha 5, DEG differential expression gene


The major conclusions derived from the current studies are that altered gastric flora and elevated bile acids might aggravate the injury to the gastric mucosa and even exacerbate Correa’s cascade process. The gastric mucosal cells might reduce bile acid toxicity by bile acid efflux and detoxification. We base this conclusion on three lines of evidence. First, we demonstrated that opportunistic pathogenic bacteria colonized the stomach mucosa of patients with cirrhosis. Next, we provided metabolic evidence to support the conclusion that gastric mucosal cells had an impaired energy metabolism. Finally, the bile secretion pathway, including genes involved bile acid efflux and detoxification, was found to upregulate in the gastric mucosal cell by transcriptome analysis.

The microbiota alteration in the skin, intestinal mucosa, ascites fluid, serum, and the oral cavity have been studied before, except for the gastric mucosa. Our data demonstrate that higher bacterial diversity and increased relative abundance of multiple bacterial genera characterize S-LC microbial dysbiosis. In our study, only Helicobacter pylori (Hp) was detected in the Helicobacter; therefore, we consider the Helicobacter to be Hp in the following. There was a greater relative abundance of Streptococcus sp. and Prevotella_melaninogenica, Neisseria spp., and Fusobacterium_periodonticum, with lower Hp. These significantly increased bacteria are pathogenic oral bacteria [19,20,21] that have the potential to elicit an inflammatory response in epithelial cells.

Interestingly, emerging findings suggest the magnitude of roles of specific oral and gastric microbiota correlated with inflammation in the development of early-stage gastric adenocarcinoma [13]. Furthermore, as the diseases develop into more severe stages, such as atrophic gastritis, intestinal metaplasia, and gastric adenocarcinoma, the dominance of Hp begins to be displaced by other bacteria, including Streptococcus, Prevotella, and other bacteria [13, 22]. Taken together, we hypothesize that bacteria in the stomach of patients with cirrhosis may originate from the oral cavity, which may induce gastric mucosal abnormalities similar to Correa’s cascade process.

Previous studies have demonstrated that BAs modulate intestinal immunity, inflammation, and tumorigenesis [23]. We found that serum BAs were predominantly conjugated, and primary BAs were elevated. BAs can alter membrane lipid composition, and increased BAs concentrations can solubilize membranes and dissociate integral membrane proteins [24]. Of note, we also found that long-chain ACs were increased and positively correlated with most gastric mucosal flora, except Hp. In the process of β-oxidation, acylcarnitine transports acyl groups (organic acids and fatty acids) from the cytoplasm to the mitochondria so that they can be broken down to produce energy for cell activities [25]. According to the Human Metabolome Database, the primary function of most long-chain acylcarnitine is to ensure long-chain fatty acid transport into the mitochondria. Blood accumulation of long-chain ACs is a marker for incomplete fatty acid oxidation. Moreau R. et al. found that acute-on-chronic liver failure was characterized by extra-mitochondrial glucose metabolism through glycolysis and depressed mitochondrial ATP-producing fatty acid β-oxidation, which may contribute to the development of organ failures [10]. The evidence above indicates impaired energy utilization in the microcirculation of cirrhotic patients who did not have other organ complications at enrollment.

We performed a small sample transcriptome analysis to validate further the alterations in gastric mucosa—upregulation of genes associated with gastric mucosa malignancy of cirrhotic patients identified in differential gene analysis. We also found that the gene encoding DEFA5 upregulates, contributing to direct antimicrobial, mucosal host defense, and immunomodulatory properties [26]. Furthermore, the over-expression of defensins in multiple types of cancer, such as colon cancer, lung cancer, and renal cell carcinomas, suggests a potential involvement of defensins in cancer development [26,27,28]. Another study indicates that DEFA5 produced from metaplastic Paneth cells may accelerate the initiation of Barrett’s esophagus, which is thought to be a precancerous lesion of esophageal adenocarcinoma [29]. An ex vivo animal study shows the down-regulation of the DEFA5 gene in gastric cancer cells and that DEFA5 inhibits the growth of gastric cancer cells [30]. The underlying mechanisms of DEFA5 in the initiation and progression of gastric cancer await further studies.

Based on the pathway enrichment analysis of DEGs from gastric mucosa transcriptome, we found that the bile secretion signaling pathway was significantly upregulated in the LC group. We also found that the genes encoding the organic solute transporterα–β (Ostα–Ostβ) upregulate. Ostα–Ostβ responsible for transporting bile acids across the enterocyte basolateral membrane into the portal circulation for subsequent renal excretion [31]. CYP3A4 is the major enzyme that catalyzes the hydroxylation of bile acids at various positions, converting bile acids into more hydrophilic and less toxic molecules [32]. Upregulation of these genes is an adaptive change in gastric mucosal cells in response to bile acid.

There are some limitations of the study that are important to the trial. First, this is a single-center study with a limited sample size. Second, the etiology of LC is complex and diverse, which may yield different flora microenvironments and serum metabolites. Meanwhile, this difference can reflect the generality of the study’s conclusions. Other factors, such as drinking and smoking, were not adjusted, which may impact gastric flora and damage to the gastric mucosa, but relevant antibiotic exposures have been excluded.

In conclusion, this study is the first work on integrated metabolomic, transcriptomic, and microbial analyses to identify critical metabolites and provide insight into the molecular and metabolic mechanisms underlying the alteration in gastric mucosa. Importantly, we showed that members of gastric pathogenic taxa were accumulated, which might originate from the oral cavity. The over-represented bacteria and serum bile acids co-exacerbated the damage of gastric mucosa and even accelerated the Correa’s cascade process. Our study highlights the enrichment of the bile secretion pathway in gastric mucosal cells in the context of LC, which may potentially serve as a protective cellular mechanism for preventing gastric lesions.


The major conclusions are that altered gastric flora and elevated bile acids might aggravate the injury to the gastric mucosa and even exacerbate Correa’s cascade process. We also provided transcriptomic evidence to support a protective mechanism for preventing gastric mucosal cells from bile acid toxicity. Understanding how microbiome and metabolite operate in the gastric mucosa should guide the development of new therapeutics for gastric abnormalities in patients with liver cirrhosis.

Availability of data and materials

All data from this work are public and have been included in this article in the main manuscript or as additional files. Additional datasets and information are available from the corresponding author under reasonable request.



Liver cirrhosis


Bile acids

long-chain ACs:

Long-chain acylcarnitines


Liquid chromatography–mass spectrometry


Linear discriminant analysis


Linear discriminant analysis effect size


Partial least squares discriminate analysis


Variable important in projection

KEGG database:

Kyoto Encyclopedia of Genes and Genomes database


False discovery rate


Fold change


Differential expression genes


Gastric mucosa in the control group


Gastric mucosa in patients with cirrhosis


Defensin alpha 5


Solute carrier family 51 subunit alpha


Solute carrier family 51 subunit beta


Cytochrome P450 3A4


Helicobacter pylori


Organic solute transporterα–β


  1. Sepanlou SG, Safiri S, Bisignano C, Ikuta KS, Merat S, Saberifiroozi M, Poustchi H, Tsoi D, Colombara DV, Abdoli A, Adedoyin RA. The global, regional, and national burden of cirrhosis by cause in 195 countries and territories, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet Gastroenterol Hepatol. 2020;5:245–66.

    Article  Google Scholar 

  2. Albillos A, de Gottardi A, Rescigno M. The gut-liver axis in liver disease: pathophysiological basis for therapy. J Hepatol. 2020;72:558–77.

    Article  CAS  PubMed  Google Scholar 

  3. Qin N, Yang F, Li A, Prifti E, Chen Y, Shao L, et al. Alterations of the human gut microbiome in liver cirrhosis. Nature. 2014;513:59–64.

    Article  CAS  PubMed  Google Scholar 

  4. Acharya C, Sahingur SE, Bajaj JS. Microbiota, cirrhosis, and the emerging oral-gut-liver axis. JCI Insight. 2017;2: e94416.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Norman K, Pirlich M. Gastrointestinal tract in liver disease: which organ is sick? Curr Opin Clin Nutr Metab Care. 2008;11:613–9.

    Article  PubMed  Google Scholar 

  6. Ibrisim D, Cevikbas U, Akyuz F, Poturoglu S, Ahishali E, Gulluoglu M, et al. Intestinal metaplasia in portal hypertensive gastropathy: a frequent pathology. Eur J Gastroenterol Hepatol. 2008;20:874–80.

    Article  PubMed  Google Scholar 

  7. Mohamed AA, Elshimy AA, El Sadik AO, Ezzat E, Nasar M, Elshaer SSM, et al. Association between severity of liver disease, frequency of Helicobacter pylori infection, and degree of gastric lesion in egyptian patients with hepatitis B virus infection. Am J Trop Med Hyg. 2018;98:221–6.

    Article  PubMed  Google Scholar 

  8. Sung JJY, Coker OO, Chu E, Szeto CH, Luk STY, Lau HCH, et al. Gastric microbes associated with gastric inflammation, atrophy and intestinal metaplasia 1 year after Helicobacter pylori eradication. Gut. 2020;69:1572–80.

    Article  CAS  PubMed  Google Scholar 

  9. Usami M, Miyoshi M, Yamashita H. Gut microbiota and host metabolism in liver cirrhosis. World J Gastroenterol. 2015;21:11597–608.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Moreau R, Clària J, Aguilar F, Fenaille F, Lozano JJ, Junot C, et al. Blood metabolomics uncovers inflammation-associated mitochondrial dysfunction as a potential mechanism underlying ACLF. J Hepatol. 2020;72:688–701.

    Article  CAS  PubMed  Google Scholar 

  11. Engelmann C, Clària J, Szabo G, Bosch J, Bernardi M. Pathophysiology of decompensated cirrhosis: portal hypertension, circulatory dysfunction, inflammation, metabolism and mitochondrial dysfunction. J Hepatol. 2021;75(Suppl 1):S49-s66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Tsochatzis EA, Bosch J, Burroughs AK. Liver cirrhosis. Lancet. 2014;383:1749–61.

    Article  PubMed  Google Scholar 

  13. Coker OO, Dai Z, Nie Y, Zhao G, Cao L, Nakatsu G, et al. Mucosal microbiome dysbiosis in gastric carcinogenesis. Gut. 2018;67:1024–32.

    Article  CAS  PubMed  Google Scholar 

  14. Kim J, Kim N, Jo HJ, Park JH, Nam RH, Seok YJ, et al. An appropriate cutoff value for determining the colonization of Helicobacter pylori by the pyrosequencing method: comparison with conventional methods. Helicobacter. 2015;20:370–80.

    Article  CAS  PubMed  Google Scholar 

  15. Yu JH, Zheng JB, Qi J, Yang K, Wu YH, Wang K, et al. Bile acids promote gastric intestinal metaplasia by upregulating CDX2 and MUC2 expression via the FXR/NF-κB signalling pathway. Int J Oncol. 2019;54:879–92.

    CAS  PubMed  PubMed Central  Google Scholar 

  16. Li L, Wang L, Song P, Geng X, Liang X, Zhou M, et al. Critical role of histone demethylase RBP2 in human gastric cancer angiogenesis. Mol Cancer. 2014;13:81.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Liu N, Zhang LH, Du H, Hu Y, Zhang GG, Wang XH, et al. Overexpression of endothelial cell specific molecule-1 (ESM-1) in gastric cancer. Ann Surg Oncol. 2010;17:2628–39.

    Article  PubMed  Google Scholar 

  18. Pei JP, Zhang CD, Yusupu M, Zhang C, Dai DQ. Screening and validation of the hypoxia-related signature of evaluating tumor immune microenvironment and predicting prognosis in gastric cancer. Front Immunol. 2021;12: 705511.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Li WH, Yu YC, Lee WS. Vertebral osteomyelitis caused by Fusobacterium periodonticum. J Infect. 2021;82:e42–4.

    Article  PubMed  Google Scholar 

  20. Iljazovic A, Roy U, Gálvez EJC, Lesker TR, Zhao B, Gronow A, et al. Perturbation of the gut microbiome by Prevotella spp. enhances host susceptibility to mucosal inflammation. Mucosal Immunol. 2021;14:113–24.

    Article  CAS  PubMed  Google Scholar 

  21. Larsen JM. The immune response to Prevotella bacteria in chronic inflammatory disease. Immunology. 2017;151:363–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Liu X, Shao L, Liu X, Ji F, Mei Y, Cheng Y, et al. Alterations of gastric mucosal microbiota across different stomach microhabitats in a cohort of 276 patients with gastric cancer. EBioMedicine. 2019;40:336–48.

    Article  PubMed  Google Scholar 

  23. Cai J, Sun L, Gonzalez FJ. Gut microbiota-derived bile acids in intestinal immunity, inflammation, and tumorigenesis. Cell Host Microbe. 2022;30:289–300.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Coleman R, Lowe PJ, Billington D. Membrane lipid composition and susceptibility to bile salt damage. Biochim Biophys Acta. 1980;599:294–300.

    Article  CAS  PubMed  Google Scholar 

  25. Tarasenko TN, Cusmano-Ozog K, McGuire PJ. Tissue acylcarnitine status in a mouse model of mitochondrial β-oxidation deficiency during metabolic decompensation due to influenza virus infection. Mol Genet Metab. 2018;125:144–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Ganz T. Defensins: antimicrobial peptides of innate immunity. Nat Rev Immunol. 2003;3:710–20.

    Article  CAS  PubMed  Google Scholar 

  27. Gobert AP, Latour YL, Asim M, Barry DP, Allaman MM, Finley JL, et al. Protective role of spermidine in colitis and colon carcinogenesis. Gastroenterology. 2022;162:813-827.e818.

    Article  CAS  PubMed  Google Scholar 

  28. Ghosh SK, McCormick TS, Weinberg A. Human beta defensins and cancer: contradictions and common ground. Front Oncol. 2019;9:341.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Nomura Y, Tanabe H, Moriichi K, Igawa S, Ando K, Ueno N, et al. Reduction of E-cadherin by human defensin-5 in esophageal squamous cells. Biochem Biophys Res Commun. 2013;439:71–7.

    Article  CAS  PubMed  Google Scholar 

  30. Wu Z, Ding Z, Cheng B, Cui Z. The inhibitory effect of human DEFA5 in growth of gastric cancer by targeting BMI1. Cancer Sci. 2021;112:1075–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Dawson PA, Hubbert M, Haywood J, Craddock AL, Zerangue N, Christian WV, et al. The heteromeric organic solute transporter alpha-beta, Ostalpha-Ostbeta, is an ileal basolateral bile acid transporter. J Biol Chem. 2005;280:6960–8.

    Article  CAS  PubMed  Google Scholar 

  32. Chen J, Zhao KN, Chen C. The role of CYP3A4 in the biotransformation of bile acids and therapeutic implication for cholestasis. Ann Transl Med. 2014;2:7.

    PubMed  PubMed Central  Google Scholar 

Download references


The authors thank all of the participants who recruited patients in this study. Thanks to Majorbio Bio-Pharm Technology Co. Ltd. (Shanghai, China) for all the help they have given us.


This project was supported by National Natural Science Foundation of China (Grant No. 82261160396), Shandong Provincial Key Research and Development Program (Major Scientific and Technological Innovation Project) (Grant NO.2021CXGC010506), and Clinical Research Center of Shandong University (Grant No. 2020SDUCRCC022). The authors declare that no funds, grants, or other support were received during the preparation of this manuscript.

Author information

Authors and Affiliations



RM and QL drafted the manuscript. LL, GY, JW, RM, and QL performed bioinformatics analyses. YL and YG performed gastroscopies. ZL, QL, RM, XX, YZ, and MD collected human samples. LL and ZL designed and supervised the study and revised the paper. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Lixiang Li or Zhen Li.

Ethics declarations

Ethics approval and consent to participate

All procedures followed were in accordance with the ethical standards of the responsible committee on human experimentation (Ethics Committee of Qilu Hospital of Shandong University, China) and with the Helsinki Declaration of 1975, as revised in 2008. The study was registered in Chinese Clinical (, ChiCTR2100051070). Informed consent was obtained from all patients for being included in the study. This article does not contain any studies with animal subjects.

Consent for publication

All authors consent to publication.

Competing interests

The authors disclose no financial or non-financial conflicts, and conflict of interest statements have been signed by each author.

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: Methods details.

Table S1. Biochemical indices of cirrhotic patients in each group. Table S2. Gastric mucosa microbiome alpha diversity index. Table S3. PLSDA-VIP table between patients with cirrhosis and controls.

Additional file 2: Figure S1.

S-C represents gastric mucosa derived from the control group, and S-LC represents gastric mucosa derived from the liver cirrhosis group. (A) Fig A is the flow chart of the study. (B) Validation of the PLS-DA model using permutation testing. (C) The difference in species composition at the specie level. *p < 0.05, **p < 0.01, and ***p < 0.001.

Additional file 3: Figure S2.

S-C represents gastric mucosa derived from the control group, and S-LC represents gastric mucosa derived from the treatment group. Wilcoxon rank-sum test bar plot on phylum level. **p < 0.01, and ***p < 0.001. S-C, gastric mucosa in the control group; S-LC, gastric mucosa in patients with cirrhosis.

Additional file 4: Figure S3.

Network analysis applied to the gastric mucosal microbiota genera network analysis reveals significant interactions. The size of the node is proportional to the species' abundance. Node color corresponds to phylum taxonoic classification. The edge color represents positive (green) and negative (red) correlations, and the edge thickness is equivalent to the correlation values. All Spearman correlations have |r|> 0.5, p < 0.05. Fig A and Fig B represent the control and liver cirrhosis groups, respectively.

Additional file 5: Figure S4.

Metabolite KEGG pathway enrichment map. Each bubble in the graph represents a KEGG Pathway. The X-axis represents the relative importance of metabolites in the pathway in terms of impact value; the Y-axis represents the enrichment significance of metabolite involvement in the pathway—log10(P value). The size of the bubble represents the Impact value. The larger the bubble, the greater the importance of the pathway. KEGG, Kyoto Encyclopedia of Genes and Genomes.

Additional file 6: Figure S5.

Cnetplot of top 11 enriched KEGG terms between control and liver cirrhosis group.

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 The Creative Commons Public Domain Dedication waiver ( 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

Ma, R., Li, Q., Yu, G. et al. A multi-omics study to investigate the progression of the Correa pathway in gastric mucosa in the context of cirrhosis. Gut Pathog 15, 45 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: