ABSTRACT Metabolites and their interactions with microbiota may be involved in Helicobacter pylori-associated gastric lesion development. This study aimed to explore metabolite alterations upon H. pylori eradication and possible roles of microbiota-metabolite interactions in progression of precancerous lesions. Targeted metabolomics assays and 16S rRNA gene sequencing were conducted to investigate metabolic and microbial alterations of paired gastric biopsy specimens in 58 subjects with successful and 57 subjects with failed anti-H. pylori treatment. Integrative analyses were performed by combining the metabolomics and microbiome profiles from the same intervention participants. A total of 81 metabolites were significantly altered after successful eradication compared to failed treatment, including acylcarnitines, ceramides, triacylglycerol, cholesterol esters, fatty acid, sphingolipids, glycerophospholipids, and glycosylceramides, with P values of <0.05 for all. The differential metabolites showed significant correlations with microbiota in baseline biopsy specimens, such as negative correlations between Helicobacter and glycerophospholipids, glycosylceramide, and triacylglycerol (P < 0.05 for all), which were altered by eradication. The characteristic negative correlations between glycosylceramides and Fusobacterium, Streptococcus, and Gemella in H. pylori-positive baseline biopsy specimens were further noticed in active gastritis and intestinal metaplasia (P < 0.05 for all). A panel including differential metabolites, genera, and their interactions may help to discriminate high-risk subjects who progressed from mild to advanced precancerous lesions in short-term and long-term follow-up periods with areas under the curve (AUC) of 0.914 and 0.801, respectively. Therefore, our findings provide new insights into the metabolites and microbiota interactions in H. pylori-associated gastric lesion progression. IMPORTANCE In this study, a panel was established including differential metabolites, genera, and their interactions, which may help to discriminate high-risk subjects for progression from mild lesions to advanced precancerous lesions in short-term and long-term follow-up. KEYWORDS: Helicobacter pylori, gastric metabolites, gastric microbiota, interactions, precancerous lesions INTRODUCTION Helicobacter pylori is an important risk factor for various gastric disorders, including chronic gastritis, glandular atrophy, intestinal metaplasia (IM), epithelial dysplasia, and even gastric cancer (GC) ([54]1). Eradication treatment has shown great benefits for GC prevention and regression of precancerous lesions ([55]2[56]–[57]4). Pathogenic mechanisms of gastric carcinogenesis involve complex interactions among H. pylori, other gastric microbiota, and associated metabolites, although systematic studies are still needed. The advances in sequencing technology have revealed altered microbial diversity and different bacterial interactions in GC and precancerous lesions ([58]5, [59]6). Our previous intervention study confirmed that H. pylori may induce gastric microbial dysbiosis, which can be restored by successful eradication ([60]7). Our subsequent follow-up study further showed that the differential bacteria in the progression-to-dysplasia/GC subjects were enriched in protein and adipose metabolism pathways by microbial functional-capacity prediction ([61]8). However, the exact microbiota-metabolite interactions still need integrative microbiome and metabolomics confirmation. Recent studies have reported that interplays between the gut microbiota and metabolites may regulate inflammation and the immune system in colorectal carcinogenesis, such as between Bifidobacterium, Lactobacillus, and short-chain fatty acids ([62]9[63]–[64]11). An integrative microbiome and metabolomics study in GC tissues found significant correlations between Helicobacter, Lactobacillus, and differential metabolites ([65]12). Nevertheless, the microbiota-metabolite interactions in the early stage of gastric lesion progression during H. pylori infection, which are important for a deeper understanding of GC etiology and prevention, still remain unclear. In the present study, we compared the metabolomics profiles of paired gastric biopsy specimens at baseline and follow-up time points from subjects with successful or failed anti-H. pylori treatment based on a prospective population-based cohort. The differential metabolites by H. pylori eradication were integrated with the differential microbiota in the same intervention subjects. This study provided us a unique opportunity to unravel the possible interactions between the gastric microbiota and metabolites in the early stage of precancerous lesion progression. RESULTS The baseline characteristics of successful eradication and failed treatment groups are presented in Table S1 in the supplemental material. No significant differences in baseline age, delta-over-baseline (DOB) value in [^13C]urea breath test (^13C-UBT), body mass index (BMI), smoking habits, alcohol consumption, or presence of active gastritis or gastric lesions were found between successful H. pylori eradication and failed treatment groups (P > 0.05 for all). There was a higher frequency of male subjects in the failed treatment group than in the successful eradication group (70.2% versus 51.7%; P = 0.043). Alterations of gastric metabolites by H. pylori eradication. A total of 267 metabolites were quantified using the Biocrates P500 platform in 230 gastric biopsy specimens, including 58 pairs before and after successful eradication and 57 pairs before and after failed treatment. Orthogonal partial least-squares discrimination analysis (OPLS-DA) found significant differences in gastric metabolomics profiles before and after successful eradication (R^2 = 0.84, Q^2 = 0.68, and P = 0.038 [[66]Fig. 1A and [67]B]). However, in the failed treatment group, gastric metabolomics profiles were not significantly changed by medical therapy (R^2 = 0.71, Q^2 = 0.53, and P = 0.167 [[68]Fig. 1C and [69]D]). FIG 1. [70]FIG 1 [71]Open in a new tab Metabolomics profiles in gastric biopsy specimens before and after anti-H. pylori treatment. Orthogonal partial least-squares discrimination analysis (OPLS-DA) found significant differences in gastric metabolomics profiles before and after successful H. pylori eradication (A); however, no significant changes were found for participants who failed to clear H. pylori by medical therapy (C). The OPLS-DA models in panels A and C were validated in panels B and D, respectively. A total of 81 metabolites showed significant alterations after successful eradication compared to failed treatment, with fold changes of >1.5 and P values of <0.05 adjusted for multiple comparison by the false-discovery rate (FDR) (Table S2). The differential metabolites include 4 acylcarnitines (Cx), 3 ceramides (Cer), 3 cholesterol esters (CE), 1 fatty acid (FA), 10 sphingolipids (SM), 6 triacylglycerols (TG), 44 glycerophospholipids, and 10 glycosylceramides. The glycerophospholipids and glycosylceramides can be further subdivided as 1 lyso-phosphatidylcholine (lysoPC), 43 phosphatidylcholines (PC), 5 hexosylceramides (HexCer), 3 dihexosylceramides (Hex2Cer), and 2 trihexosylceramides (Hex3Cer). We investigated the potentially relevant influence factors for gastric metabolite alterations after H. pylori eradication. To represent the overall metabolic status in each gastric biopsy specimen, we calculated a comprehensive metabolic index using the 81 differential metabolites. Potentially relevant factors include DOB value in ^13C-UBT representing H. pylori infection status, microbial Shannon and Richness indexes representing gastric microbial diversity, and gastric juice pH value. [72]Figure 2A shows significant increases in microbial diversity indexes (P < 0.001 for both) and decreases in pH values (P = 0.001) as well as in metabolic indexes (P < 0.001) accompanying the dramatically decreasing trend of DOB values (P < 0.001) after successful eradication. However, except for the fact that pH values were increased (P = 0.020) in the failed treatment group, the factors showed no significant alterations (P < 0.05 for all). FIG 2. [73]FIG 2 [74]Open in a new tab Potentially relevant influence factors and pathway enrichment analyses for gastric metabolite alterations after H. pylori eradication. (A) Paired plots showed significant decreasing trends of DOB values (P < 0.001), pH values (P = 0.001), and metabolic indexes (P < 0.001) and significant increasing trends of microbial diversity indexes (P < 0.001 for both) after successful H. pylori eradication. However, in the failed treatment group, except that pH values were increased significantly (P = 0.020), the factors showed no significant alterations (P > 0.05 for all). a, metabolic index was calculated using the 81 differential metabolites after H. pylori eradication to represent the overall metabolic status in each gastric biopsy specimen by the logistic regression equation in Materials and Methods. DOB, delta over baseline. (B) A total of 26 pathways were enriched by both the differential metabolites and microbiota after successful eradication. Among them, 5 pathways were significantly changed by eradication in both metabolomics and microbiome analyses, including pathways for leishmaniasis, adipocytokine signaling, lipid and atherosclerosis, choline metabolism in cancer, and cholesterol metabolism (P values <0.05 for all). *, P values < 0.05; **, P values < 0.01. The putative functions of the 81 differential metabolites after H. pylori eradication were found to be enriched in 29 pathways using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. We compared the 29 metabolic pathways with the previously predicted microbiota functional-capacity changes ([75]7) in the same intervention participants. A total of 26 pathways were enriched both by the differential metabolites and microbiota, covering 5 main classes, such as cellular processes, environmental information processing, human diseases, metabolism, and organismal systems ([76]Fig. 2B). Among them, 5 pathways were significantly changed after eradication in both metabolomics and microbiome analyses, including leishmaniasis, adipocytokine signaling, lipid and atherosclerosis, choline metabolism in cancer, and cholesterol metabolism pathways (P < 0.05 for all). Interactions between differential metabolites and microbiota before and after H. pylori eradication. In our microbiome analysis, 65 bacterial taxa were significantly changed after successful H. pylori eradication ([77]7). We used Spearman’s correlation analysis to assess the potential interplays between the 81 differential metabolites and the 65 differential taxa in the 91 participants with both microbiome and metabolomics results. We found 71 negative (−0.29 < r < −0.21; P < 0.05 for all) and 177 positive (0.21< r < 0.28; P < 0.05 for all) significant correlations between 34 metabolites and 58 taxa in H. pylori-positive baseline biopsy specimens ([78]Fig. 3A). After successful eradication, stronger correlations were found between 30 metabolites and 41 taxa, including 49 negative (−0.44 < r < −0.34; P < 0.05 for all) and 44 positive (0.34 < r < 0.52; P < 0.05 for all) correlations ([79]Fig. 3B). After failed treatment ([80]Fig. 3C), 20 metabolites and 63 taxa were significantly correlated, with 163 negative associations (−0.52 < r <−0.26; P < 0.05 for all) and 60 positive associations (0.26 < r < 0.48; P < 0.05 for all). FIG 3. [81]FIG 3 [82]Open in a new tab Interactions between differential metabolites and the microbiota before and after anti-H. pylori treatment. The heat maps show correlations between 81 differential metabolites and 65 previously identified differential bacterial taxa in biopsy specimens of subjects at H. pylori-positive baseline (A), after successful eradication (B), or after failed treatment (C). *, P values <0.05. The correlation networks showed the significant associations between differential metabolites and previously identified differential genera in biopsy specimens of subjects at H. pylori-positive baseline (D), after successful eradication (E), or after failed treatment (F). Red lines represent positive genus-metabolite correlations. Blue lines represent negative genus-metabolite correlations. Green nodes represent differential metabolites. Orange nodes represent gastric microbiota. Node radii are based on the number of significant genus-metabolite correlations. Cx, acylcarnitines; Cer, ceramides; CE, cholesterol esters; FA, fatty acids; HexCer, hexosylceramides; Hex2Cer, dihexosylceramides; Hex3Cer, trihexosylceramides; PC, phosphatidylcholines; lysoPC, lyso-phosphatidylcholine; SM, sphingomyelins; TG, triacylglycerols. Of the 65 differential taxa, we further focused on the 18 differential genera after H. pylori eradication ([83]7). Significant interactions between differential genera and metabolites were visualized by correlation network construction. Helicobacter was negatively correlated with glycerophospholipids, glycosylceramide, and triacylglycerol (−0.24 < r < −0.21; P < 0.05 for all) in baseline biopsy specimens ([84]Fig. 3D), while it was positively correlated with glycerophospholipids and sphingolipids (0.37 < r < 0.51; P < 0.05 for all) after successful eradication ([85]Fig. 3E). Positive correlations between triacylglycerols and Prevotella (0.21