Different effects of hydrogen

Scientific Reports volume 12, Article number: 7231 (2022) Cite this article

3797 Accesses

3 Citations

2 Altmetric

Metrics details

The potential for preventive and therapeutic applications of H2 have now been confirmed in various disease. However, the effects of H2 on health status have not been fully elucidated. Our previous study reported changes in the body weight and 13 serum biochemical parameters during the six-month hydrogen intervention. To obtain a more comprehensive understanding of the effects of long-term hydrogen consumption, the plasma metabolome and gut microbiota were investigated in this study. Compared with the control group, 14 and 10 differential metabolites (DMs) were identified in hydrogen-rich water (HRW) and hydrogen inhalation (HI) group, respectively. Pathway enrichment analysis showed that HRW intake mainly affected starch and sucrose metabolism, and DMs in HI group were mainly enriched in arginine biosynthesis. 16S rRNA gene sequencing showed that HRW intake induced significant changes in the structure of gut microbiota, while no marked bacterial community differences was observed in HI group. HRW intake mainly induced significant increase in the abundance of Lactobacillus, Ruminococcus, Clostridium XI, and decrease in Bacteroides. HI mainly induced decreased abundances of Blautia and Paraprevotella. The metabolic function was determined by metabolic cage analysis and showed that HI decreased the voluntary intake and excretions of rats, while HRW intake did not. The results of this study provide basic data for further research on hydrogen medicine. Determination of the effects of hydrogen intervention on microbiota profiles could also shed light on identification of mechanism underlying the biological effects of molecular hydrogen.

Hydrogen (H2) is the smallest and lightest gas molecule, which has been historically considered as a biologically inert molecule. Early in 1975, Dole et al. firstly reported the possible anti-cancer effect of hyperbaric treatment of 97.5% hydrogen gas in a mouse model of skin tumor1. However, medical researchers did not pay considerable attention to H2 until Ohsawa et al. reported that inhalation of 1–4% H2 gas significantly attenuates cerebral ischemia–reperfusion injury in rats by selectively neutralizing hydroxyl radicals and peroxynitrite2. The potential for preventive and therapeutic applications of H2 have now been confirmed in more than 170 different human and animal-disease models, including ischemia–reperfusion (I/R) injuries3,4, neurodegeneration5,6, cardiovascular diseases7,8, metabolic syndrome9,10, inflammation11,12, and cancer13,14. Several biological mechanisms have been proposed, including selectively reducing cytotoxic oxygen radicals2, anti-inflammatory effects15, recovering mitochondrial dysfunction16, regulation of endoplasmic reticulum stress17, but none of them can fully explain the multiple biological functions of H2.

In mammals, the gut microbiome forms a complex ecosystem consisting of a vast number of interacting bacteria, archaea, bacteriophages, eukaryotic virus and fungi, most of which are commensal or mutualistic microorganisms18. In the past decade the gut microbiota has been proved to play a profound role in the training of host immunity, digesting food, regulating gut endocrine function and neurological signalling, modifying drug action and metabolism, eliminating toxins and producing numerous compounds that influence the host19. At present, the research on the relationship between hydrogen consumption and gut microbiome is relatively limited. Most studies showed that hydrogen-rich water (HRW) could improve intestinal structural integrity and upregulation of butyrate-producing bacteria with ameliorated clinical features of gut microbiota disturbances20. However, these studies have been primarily focused on the modulatory effect of HRW consumption on intestinal flora in pathological conditions, whether administration of HRW regulates gut microbiome in healthy animals remains largely unknown. In addition, as another commonly used method of hydrogen consumption, whether hydrogen inhalation could also affect gut microbiome need to be further investigated.

One of the aims of the present study was to explore the possible regulatory effects of long term of HRW intake and hydrogen inhalation on gut microbiome. Our previous study showed that long term of HRW intake or hydrogen inhalation can influence some serum biochemical parameters of normal rats, indicating the potential modulatory effect of hydrogen consumption on metabolism in health status21. To further investigate the effect of hydrogen consumption on metabolism, LC–MS based pseudotargeted metabolomics analysis was performed to determine the changes in the levels of plasma metabolites. The relationship between differentially expressed plasma metabolites and altered gut microbiota was also investigated.

As reported in our previous study21, after 6-month hydrogen intervention, the body weights (BW) of rats in HI group were significantly decreased (455.20 ± 31.57 g vs. 525.00 ± 17.78 g, p = 0.007) compared with the controls, while no significant changes in BW were observed in HRW group. The BW of rats were lower in HI group than in HRW group (455.20 ± 31.57 g vs. 498.80 ± 18.50 g, p = 0.012) (Fig. 1A).

The effects of hydrogen intervention on the BW and metabolic function of rats. (A) BW changes at six month after hydrogen treatment; (B) BW changes during the first month after hydrogen treatment; (C–F) the effects of hydrogen treatment on food intake (C), water intake (D), defecation (E), and urinate (F).

To determine the effects of hydrogen intervention on metabolic function, we conducted another study using metabolic cages to measure 24-h voluntary intake and excretions of rats during the first month of hydrogen intervention. Compared with the control group, both HRW and HI had no significant effect on the BW of rats during the four weeks, however, the BW was significantly lower at the first (263.29 ± 19.98 g vs. 296.56 ± 19.69 g, p = 0.012) and the fourth week (383.21 ± 29.36 g vs. 428.81 ± 31.21 g, p = 0.024) after HI. Compared with HRW group, the BW of rats in HI group was lower with a marginal significance (p = 0.063) during the four weeks (Fig. 1B). As shown in Fig. 1C,D, compared with the control group, HI induced a significant decrease in both food intake (p < 0.0001) and water intake (p < 0.0001), while no reduction was observed in HRW group. During the four weeks, the defecation (p = 0.002, Fig. 1E) and urination (p = 0.006, Fig. 1F) was also markedly decreased in HI group, while no effect was observed in HRW group. Spearman’s correlation analysis revealed that the reduced defecation of rats in HI group was positively related to the reduction of food intake (r = 0.43, p = 0.006), and the decreased urinate was also positively related to the reduced water intake (r = 0.64, p < 0.0001).

To determine the effects of hydrogen intervention on plasma metabolites, LC–MS-based pseudotargeted metabolomics analysis was performed on fasting plasma samples. Eighty-six plasma metabolites were identified consisting of amino acids and their derivatives, intermediates in glycolysis and the citric acid cycle, lipid metabolites, nucleotide metabolites, urea cycle metabolites, carbohydrates, co-factors/vitamins, and hormones. To identified the difference in the levels of plasma metabolites between groups, Principal component analysis (PCA) and orthogonal partial least-squares discriminant analysis (OPLS-DA) were performed. PCA analysis showed that HRW and CON groups were clearly segregated by principle components (PC)1 (18.91%) and 2 (15.26%) (Fig. 2A). The segregation of HI and CON (PC1: 18.94%, PC2: 14.45%, Fig. 2B), as well as HRW and HI (PC1: 24.69%, PC2: 11.68%, Fig. 2C) was also observed. The OPLS-DA model from the plasma metabolic profile showed a good discrimination between HRW and CON (R2X = 0.262, R2Y = 0.991, Q2 = 0.645, Fig. 2D,G), HI and CON (R2X = 0.269, R2Y = 0.968, Q2 = 0.359, Fig. 2E,H), HRW and HI (R2X = 0.325, R2Y = 0.976, Q2 = 0.684, Fig. 2F,I).

The difference on plasma metabolite profile between groups was identified by both PCA and OPLS-DA analysis. (A) PCA plot of HRW vs. CON; (B) PCA plot of HI vs. CON; (C) PCA plot of HRW vs. HI; (D) OPLS-DA score plot of HRW vs. CON; (E) OPLS-DA score plot of HI vs. CON; (F) OPLS-DA score plot of HRW vs. HI; (G–I) Validation plot obtained from 200 permutation tests, respectively.

Differential metabolites (DMs) were selected by using the cutoff of OPLS-DA VIP score > 1.0 with a p value < 0.05 in the fold change of expression level between any two of the three groups. Thirty-five DMs were identified as shown in Supplementary Table 1. Compared with the control group, there are 14 and 10 DMs in HRW and HI group respectively. Twenty-two DMs were identified between HRW and HI group. Compared the control group, all the DMs were down-regulated in HRW group, while all the DMs were up-regulated in HI group. Compared HRW group, all the DMs were up-regulated in HI group. As shown in Fig. 3A, the dendrogram of hierarchical clustering showed the plasma samples in HRW group was clustered separately from the control group or HI group, however, the difference between HI group and the control group was much smaller. The pathway enrichment analysis based on metabolite quantitative alterations was performed by the MetaboAnalyst 5.0 (http://www.metaboanalyst.ca). The metabolic pathways with impact value > 0.1 and − log(p) > 2.0 are considered the most relevant pathways involved in the conditions under study. The results showed that the DMs between HRW and control group were mainly concentrated in starch and sucrose metabolism (Fig. 3B), the DMs between HI and control group were mainly involved in arginine biosynthesis (Fig. 3C), the DMs between HRW and HI group were mainly enriched in glycerolipid metabolism, inositol phosphate metabolism, starch and sucrose metabolism, glyoxylate and dicarboxylate metabolism, and ascorbate and aldarate metabolism (Fig. 3D).

Analysis of differential metabolites in plasma samples among three groups. (A) The hierarchical clustering results for differential metabolites. (B–D) The functional enrichment analysis for differential metabolites between two groups. NAAG N-Acetylaspartylglutamic acid, CMP Cytidine 5′-monophosphate.

To investigate the effects of hydrogen intervention on faecal microbiota structure, we analyzed the bacterial communities on the samples by targeted 16S rRNA gene (V3–V4 region) sequencing using the Illumina MiSeq. The high quality sequence count (CON: 35,896 ± 766, HRW: 36,063 ± 884, HI: 35,630 ± 845) and operational taxonomic units (OTUs) count (CON: 551 ± 76, HRW: 516 ± 11, HI: 559 ± 84) were similar in the three groups. Four measures of α-diversity [Chao1 (Fig. 4A), PD whole tree (Fig. 4B), Shannon (Fig. 4C), and Simpson (Fig. 4D)] all fail to report significant difference between the control group and HRW or HI group. However, the PD whole tree index was significantly higher in HI group than that in HRW group (Fig. 4B). Analysis of similarities (ANOSIM) was performed to quantitatively compare the bacterial community differences between different groups. As shown in Fig. 4E, ANOSIM revealed significant difference in the structure of gut microbiota between HRW and the control group (ANOSIM, r = 0.304, p = 0.005). Significant bacterial community differences was also observed between HRW and HI group (Fig. 4G, ANOSIM, r = 0.369, p = 0.003). However, no significant difference was observed between HI and the control group (Fig. 4F, ANOSIM, r = − 0.065, p = 0.65). The Wilcoxon rank-sum test was performed to identify differences in relative abundances of bacteria between the two groups. At the phyla levels, as shown in Fig. 4H, the dominant phyla of the three groups were Firmicutes, Bacteroidetes, Proteobacteria, and Actinobacteria. Compared with the controls, the relative percentage of Elusimicrobia was significantly increased in HRW group, while the proportion of Deferribacteres and Euryarchaeota were markedly decreased in HRW group. HI only induced markedly increase in the proportion of Elusimicrobia. Compared with HRW group, the proportion of Elusimicrobia was significantly lower and Spirochaetes and Euryarchaeota were markedly higher in HI group. At the family level, as shown in Fig. 4I, the dominant family of the faecal microbiota were Prevotellaceae, Ruminococcaceae, Lachnospiraceae, and Lactobacillaceae. The abundances of Lactobacillaceae, Peptostreptococcaceae, and Elusimicrobiaceae were increased, and the abundances of Bacteroidaceae and Desulfovibrionaceae were decreased in HRW group compared with the controls. The decreased abundance in Porphyromonadaceae and increased abundance in Elusimicrobiaceae were observed in HI group compared with the controls. Compared with HRW group, the increased abundance in Acidaminococcaceae and decreased abundances in some other phyla (e.g., Lactobacillaceae, Porphyromonadaceae, Peptostreptococcaceae, and Sphingomonadaceae) were observed in HI group. At the genus levels, as shown in Supplementary Table 2, compared with the control group, eleven and seven genera exhibited significantly different abundances in HRW and HI group respectively. Compared with HRW group, twenty-three genera showed significant differences in relative abundances of bacteria in HI group.

Analysis of bacterial community composition among three groups. (A–D) Alpha-diversity of bacterial community was measured by Chao1 (A), PD whole tree (B), Shannon (C), and Simpson (D); (E–G) ANOSIM was performed to quantitatively compare the bacterial community differences between different groups. (H,I) Community bar-plot analysis shows relative abundance of microbiota in each group at the phylum level (H) and family level (I), n = 6, in each group.

To determine whether the changes in gut microbial taxa would alter the gut microbiota function, the functional prediction at level-3 of KEGG pathway was performed by PICRUSt. As shown in Fig. 5, compared with the control group, the faecal microbiota of rats from HRW group had elevated pathways involved in ribosome biogenesis and replication, recombination and repair proteins. HI group had significantly enriched pathway involved in phenylalanine, tyrosine and tryptophan biosynthesis. Compared with HRW group, HI group had reduced pathways involved in ribosome biogenesis, chromosome, and recombination and repair proteins.

Differentially enriched KEGG pathways. The functional prediction at level-3 of KEGG pathway was performed by PICRUSt between HRW and Control group (A), HI and Control group (B), HRW and HI group (C).

To examine the possible relationship between differentially expressed plasma metabolites and altered gut microbiota at the genus level, multiomic analysis was performed by using Spearman’s correlation test. The results of Spearman’s correlation analysis between HRW and the control group are shown in Fig. 6A, the abundance of urocanic acid had the strongest correlation with Bacteroides (r = 0.944, p < 0.0001). l-lactic acid showed comparable moderate positive correlations with Desulfovibrio and Bacteroides. N-Acetyl-l-alanine and N-Acetyl-d-glucosamine also showed moderate positive correlation with Bacteroides and Elusimicrobium. d-Glucose 6-phosphate showed moderate negative correlations with Elusimicrobium, Clostridium XI, and Barnesiella. Figure 6B showed the Spearman’s correlation analysis between HI and the control group, malic acid had the strongest correlation with Blautia (r = 0.879, p = 0.0002). l-leucine showed moderate positive correlation with Blautia. N-Carbamoyl-l-aspartic acid, cis-Aconitic acid, and inosine also showed moderate positive correlations with Paraprevotella. N-Carbamoyl-l-aspartic acid showed moderate negative correlation with Bifidobacterium. Figure 6C showed the Spearman’s correlation analysis between HRW and HI group, l-phenylalanine showed high positive correlations with Sphingomonas (r = 0.866, p = 0.0003), Methylocacterium (r = 0.921, p < 0.0001), and Bosea (r = 0.848, p = 0.0005). Glyceric acid also showed high positive correlations with Lactobacillus (r = 0.937, p < 0.0001) and Mycoplasma (r = 0.881, p = 0.0002). N6-methyladenosine showed high positive correlation with Lactobacillus (r = 0.895, p < 0.0001), while had strong negative correlation with Methanosphaera (r = − 0.834, p = 0.0007).

Correlation analyses between significantly different genera and differential plasma metabolites. Spearman’s correlation test was performed between HRW and Control group (A), HI and Control group (B), HRW and HI group (C).

Our previous study reported that HI can induce significant decrease in BW of rats, while HRW did not21. Consistent with previous study, the results of this study also showed markedly decrease in BW in HI group, while no obvious effect was observed in HRW group. The metabolic cage experiments showed that HI could significantly decrease the food intake, water intake, defecation and urination. The decreased metabolic function may contributed to the reduction of BW in HI group. Our previous study also showed that both HRW intake and HI can induce significant changes in several serum biochemical parameters in normal rats21. To obtain a more comprehensive understanding of the metabolic alterations in response to HRW intake or HI, LC–MS based pseudotargeted metabolomics analysis was performed. The OPLS-DA models indicated clear separations between any two of the three groups based on their metabolomic responses. Compared with the control group, 14 and 10 DMs were identified in HRW and HI group, respectively. It is worth noting that all the DMs in HRW group were down-regulated, while all the DMs were up-regulated in HI group, 22 DMs were identified between HRW and HI groups, indicating that the modulatory effects of HRW intake on metabolism differ markedly from HI. Further functional enrichment analysis suggested the DMs were mainly involved in starch and sucrose metabolism and arginine biosynthesis in HRW and HI group, respectively.

Previous study showed that 4 weeks of HRW intake could significantly decrease the levels of blood glucose, lactate, and blood urea nitrogen (BUN) and exert antifatigue effects in chronic forced swimming mice22. It has also been reported that 3 months of HRW intake could markedly decrease the blood uric acid levels in male patients with hyperuricemia23. Consistent with these findings, our results also showed that 6 months of HRW intake could down-regulate plasma levels of d-glucose 6-phosphate, l-lactic acid, and uric acid. In addition, the plasma levels of nucleotides and their derivatives were also reduced after HRW intake, indicating the regulatory effects of HRW intake on nucleotides metabolism. Notably, all the changed nucleotides derivatives, including m6A, pseudouridine, and N2,N2-dimethylguanosine, were belong to RNA modifications. Among them, m6A is the most widespread epigenetic modification on mammalian mRNA and has been shown to act as a key regulator of numerous important biological processes in normal physiology and in disease, including cancer, heart failure, viral infection, and type 2 diabetes24,25. Pseudouridine was reported to undergo dynamic changes in response to serum starvation, hydrogen peroxide and heat shock in mammalian cells24. It has previously been reported that H2 administration could regulate expression of diverse genes26, our results suggest that molecular hydrogen may regulate gene expression by affecting epigenetic modifications. In addition, the plasma levels of N1-methyl-2-pyridone-5-carboxamide was significantly decreased by HRW intake. Nicotinamide adenine dinucleotide (NAD+) is an important coenzyme for redox reactions, making it central to energy metabolism27. Nicotinamide mononucleotide (NMN) is one of the substrates for NAD+ synthesis, which can be further metabolized to N1-methyl-2-pyridone-5-carboxamide. Thus, HRW intake may regulate energy metabolism via affecting NAD+ synthesis.

Unlike HRW intake, HI had an up-regulatory effects on plasma DMs. Six of those DMs, including l-citrulline, l-leucine, sarcosine, l-glutamic acid, N-carbamoyl-l-aspartic acid, and NAAG, were belong to amino acids and their derivatives. Among them, the plasma levels of l-citrulline showed the most considerable increase in HI group compared to the controls. It has been showed in rats that only the intestine produced circulating citrulline, and the increased citrulline levels may be caused by either increased production or decreased utilization28. The decreased utilization could arise from a decrease in clearance, i.e. renal failure. However, no significant change of creatinine plasma levels was observed in our study, indicating that the increased citrulline levels was not caused by the impairment of renal function. Previous study showed that the increase in plasma citrulline was correlated with protein absorption improvement in patients with short bowel syndrome (SBS) followed in the first year after resection29. We supposed that the increased citrulline plasma levels may be associated with the improvement of enterocyte function, although this need to be further investigated. HI also induced a significant increase in NAAG plasma levels. NAAG is the most prevalent and widely distributed dipeptide in the mammalian nervous system30. The levels of NAAG in plasma and cerebrospinal fluid (CSF) were much lower than those in brain tissues31. It is well established that the increase in NAAG is neuroprotective against N-methyl-d-aspartate (NMDA) receptor-mediated neurotoxicity, including ischemic brain injury30. HI was first reported to exert neuroprotective effects on ischemic stroke, and further studies also found its protective effects on other neurological impairment, including traumatic brain injury, subarachnoid hemorrhage, and neurodegenerative diseases32. Consistent with our results, the HI induced increase in NAAG was also observed in cortex tissues of mice with ischemic stroke33. Although the ability of molecular hydrogen to scavenge hydroxyl radicals may partly explain its neuroprotective effects, the regulatory effects on NAAG may also responsible for its protective benefits. In addition, the significant increase in plasma levels of two citric acid cycle intermediates, cis-aconitic acid and malic acid, were also observed in our study, indicating that HI may accelerate mitochondrial energy metabolism. N-carbamoyl-l-aspartic acid and inosine are intermediates of pyrimidine and purine metabolism, respectively. The increase in the two metabolites suggest that HI may have modulatory effects on nucleoside metabolism. The plasma levels of other DMs, including l-leucine, sarcosine, l-glutamic acid, and l-lactic acid, varied slightly by HI.

Recent studies have provided evidence that modulation of host gut microbiota may be one of the mechanisms contributing to the biological effects of exogenous hydrogen consumption. Qiu et al. showed that saturated hydrogen saline treatment could modulate the abundance of Bacteroides, Bifidobacteria, and Lactobacillus in feces, which may responsible for the improvement of lipid metabolism disorders in high-fat diet mice34. Jin et al. reported that sustained H2 release in the gut by hydrogen nanocapsule could increase the abundance of Akkermansia muciniphila and attenuate metabolic dysfunction-associated fatty liver disease35. In this study, HRW intake induced significant changes in the structure of gut microbiota, while no marked bacterial community differences was observed in HI group. Previous study showed that the peak of the hydrogen concentration in small intestine after oral intake of 5 ppm of HRW was approximately 20 times higher than that after inhalation of 4% hydrogen gas36. The significant difference in hydrogen concentration in intestine between HRW intake and HI may contribute to the different effects on microbiota composition.

In our study, HRW intake induced significant increase in the proportion of Lactobacillus, Ruminococcus, Clostridium XI, Elusimicrobium, Barnesiella, and Aquabacterium, and decrease in Bacteroides, Anaerotruncus, Desulfovibrio, Mucispirillum, and Bifidobacterium. Lactobacillus and Bifidobacterium are the most common probiotic bacteria with the reported beneficial effects including aid digestion, reduce constipation, resist infections, prevent traveler’s diarrhea and ameliorate intestinal bowel disease (IBD)37. In our study, HRW intake induced a significant increase in the abundance of Lactobacillus. The increased abundance of Lactobacillus induced by HRW intake may contribute to the beneficial effects of HRW. Although HRW also induced a marked decrease in Bifidobacterium, the relative abundances of Bifidobacterium is very low. A recent clinical study showed that drinking hydrogen-dissolved alkaline electrolyzed water (AEW) for two weeks induced an increase in Bifidobacterium in healthy volunteers38. The different pH values (HRW 7.5 vs. AEW 9.5) or duration of hydrogen treatment (HRW 6 months vs. AEW 2 weeks) may contribute to the different effects on the levels of Bifidobacterium. It has been shown that supplementation with Ruminococcus flavefaciens could attenuate the antidepressant effects of duloxetine on depressive-like behavior39, although the increased Ruminococcus has also been shown to be beneficial regarding antidepressant-induced constipation39. Previous study showed that 4 weeks of HRW intake could exert beneficial effects on depressive-like behavior in mice via suppression of the inflammasome activation40, the antidepressive effects of HRW may be diminished by long-term HRW intake induced increase in Ruminococcus according to our study. Clostridium XI belongs to class Clostridia, which have been reported to attenuate inflammation and allergic diseases41. It has also been demonstrated that Clostridium species can utilize indigestible polysaccharide and produce lots of short-chain fatty acids (SCFAs), which are now considered as key players in the interactions with the host that impact on health and disease, especially given recent evidence for their capacity to modify the epigenome and effects on tissues and organs beyond the gut42. The increase in Clostridium XI derived SCFAs may also contributed to the effects of HRW. In a study of 345 Chinese individuals, members of the genera Bacteroides has been shown to be more abundant in type II diabetic subjects compared to controls with normal glucose metabolism43. The improved glucose tolerance and hyperglycemia lowering effect of HRW intake have been previously reported9,44, which may be attributed by the HRW-induced decrease in Bacteroides levels. The results of the Spearman correlation analyses revealed a great number of significant correlations between the abundant of Bacteroides and DMs, including urocanic acid, l-lactic acid, N-acetyl-d-glucosamine, and N-Acetyl-l-tyrosine, however, the causal relationships between alterations in Bacteroides abundance and plasma DMs need to be further investigated. Although the levels of other genera, including Elusimicrobium, Barnesiella, Aquabacterium, Anaerotruncus, Desulfovibrio, Mucispirillum, and Bifidobacterium changed significantly, their relative abundance was very low.

Compared with HRW group, the changes in fecal microbiota were found to be much less in HI group. Among these changed genera, the abundances of Blautia and Paraprevotella were significantly increased. Blautia has been considered as a probiotic bacterium that occur widely in mammalian feces and intestines45. The increased abundances of Blautia may contributed to the effects of HI. The Paraprevotella genus has been found to be negatively correlated with the BMI index46. Consistently, the negative correlation was also observed between the abundance of Paraprevotella genus and the BW of rats. Spearman correlation analyses revealed a great number of significant correlations between the abundant of Blautia and DMs, Correlation analyses revealed significant negative correlations between the abundance of Blautia and DMs, including l-leucine and malic acid. The abundance of Paraprevotella was negatively correlated with DMs, including inosine, cis-aconitic acid, and N-carbamoyl-l-aspartic acid. The other changed genera, including Elusimicrobium, Propionibacterium, Porphyrobacter, Methanosphaera, and Bifidobacterium, all had relatively low abundance.

Although both HRW and HI have beneficial effects in various diseases, the underlying mechanisms of them may be different. For HRW intake, it has been reported that many kinds of diseases, especially gastrointestinal symptoms, such as constipation and diarrhea, could benefit from HRW. It is known that the gut microbiota plays an important role in health and disease, the modulatory effect of HRW on gut microbiota may significantly contribute to the improvement of these diseases. In our study, the results of metabolic cage experiments showed that HI induced significant changes in food/water intake of rats, which positively correlates with changes in defecation/urination, while no obvious effect was observed in HRW group. The reasons for the decline in food intake are multifactorial and involve both peripheral and central mechanisms47. Hyspler et al. used deuterium gas as a metabolic tracer to quantify hydrogen metabolism in the mammalian body and found that the deuterium can be oxidized to water48. In our study, the inhaled hydrogen gas may probably be oxidized to water, providing an endogenous supply of H2O, which might explain the decrease in the water intake of rats. Moreover, in our previous studies we have provided preliminary evidence that eukaryocytes have the ability of hydrogen metabolism, which may affect its metabolic activities49,50. It is, therefore, possible to hypothesize that HI may exert its biologic role through modulation of metabolic function of mitochondria, although this needs to be further investigated.

One limitation of the present study is that research on the effects of hydrogen intervention on faecal microbiota profiles was focused on the genus level, and did not conduct in-depth studies at the species or even strain levels. Another limitation is the lack of proven causal relationships between alterations in microbiota profiles and plasma metabolites. In addition, further study should also evaluate the effects of hydrogen intervention on the production of SFCAs, which play a key role in microbiota-host interactions.

Collectively, the results of this study could provide basic data for further research on hydrogen medicine. Our results also shed light on the effects of different routes of hydrogen intervention on microbiota profiles, which may significantly contributed to the therapeutic effects of hydrogen in various diseases.

The methods of animal experiment has been reported in our previous study as follows21: eighteen 3-week-old male Sprague–Dawley rats weighing 40–50 g were purchased from the Vital River Laboratory Animal Technology Co., Ltd (Beijing, China). Rats were housed under a constant temperature at 22–25 °C with a 12 h light–dark cycle and maintained on a normal diet. All procedures were approved by the Institutional Animal Experiment Committee of the Chinese PLA General Hospital and were conducted in compliance with the Regulations for the Administration of Affairs Concerning Experimental Animals (China) and the ARRIVE guidelines. Before the experiment, rats were adapted to laboratory conditions for one week. Rats were then randomly divided into three groups (6 in each group): (1) Control group: rats were maintained under normal conditions; (2) HRW group: rats were given HRW by oral intake for 1 h each time, two times per day; (3) HI group: rats were treated with HI (4%) for 1 h each time, two times per day. The experiment was last for six months. Fresh feces samples were collected via sterile operation before sacrifice, and were immediately stored at − 80 °C for microbiota analysis. After euthanasia, blood samples were immediately harvested from the portal vein. Plasma samples were obtained by centrifuging the whole blood at 1500g for 10 min at 4 °C, and stored at − 80 °C for pseudotargeted metabonomic analysis.

HRW (H2 concentration > 800 µM) was kindly provided by Shenzhen Kelieng Biomedical Co. Ltd. (Shenzhen, China) and stored under atmospheric pressure at 23 ± 2 °C in a stainless steel bucket (KLE-8). The hydrogen concentration was monitored using a hydrogen electrode (Unisense A/S, Aarhus, Denmark), ensuring that the hydrogen concentration of HRW for rats was maintained above 800 µM.

Rats were placed in a transparent closed box (72 × 53 × 45 cm, length × width × height) connected to the hydrogen gas generator which composed of an Oxy-Hydrogen Machine (SG-3000; Gang’an Health Management [Beijing] Co., Ltd., Beijing, China) and a gas mixer, allowed to spontaneous respiration (4% H2, 96% air containing 21% O2) for 1 h each time, and two times per day. The concentration of hydrogen and oxygen in the closed box was monitored by Thermal trace GC ultra-gas chromatography (Thermo Fisher, MA, USA).

Male Sprague–Dawley rats (220 ± 20 g) were purchased from the Vital River Laboratory Animal Technology Co., Ltd (Beijing, China) and randomly divided into three groups (n = 8 in each group): CON, HRW, and HI group. Metabolic function of rats were monitored during 4 weeks. Rats were housed individually in metabolic cages (Techniplast, UK) once a week. Food intake was calculated by the remaining weight of the food from the weight of food provided on the previous day. Water intake was calculated by subtracting the amount of water left from the measured amount of water provided on the previous day. The feces and urine during the 24 h was collected and measured. The BW was monitored during the four weeks. The data were analysed by two-way ANOVA repeated measures with Tukey’s post hoc test using GraphPad Prism 8.0.2 and presented as mean ± SD. Spearman correlation analysis was performed using GraphPad Prism 8.0.2.

To identify the DMs in plasma, the pseudo-targeted metabolomics approach was performed, which allows to identify and quantify 200 plasma metabolites. The plasma samples were thawed at 4 °C and thoroughly vortexed. For each sample, 100 µL was taken and transferred into a 1.5 mL Eppendorf tube. Then 400 µL ice-cold methanol/acetonitrile (1:1, v/v) was added. The samples were mixed by vigorous vortexed, sonicated in an ice-cold water bath sonicator for 20 min, and kept at − 20 °C for 1 h. Next, the samples were centrifuged at 14,000g for 20 min at 4 °C. The supernatant was transferred into a glass sampling vial and vacuum-dried at 4 °C. Then 100 µL of acetonitrile/water (1:1, v/v) was added before centrifugation at 14,000g for 20 min at 4 °C. The resultant supernatants were used for metabolomics analysis. The quality control (QC) samples were prepared by mixing equal volumes of all plasma samples with all other steps as described above.

The chromatographic separation of plasma was carried out on a Waters I-class liquid chromatography system with ACQUITY UPLC BEH Amide column (1.7 µm, 2.1 × 100 mm column, Waters). The mobile phase system consists of 25 mM ammonium acetate + 25 mM ammonia (pH 9.75) (A) and acetonitrile (B). Column temperature was maintained at 40 °C, and flow rate was 0.3 mL/min. Sample injection volume was 2 µL. The linear gradient for mobile phase B was as follows: 0–1 min 95% B; 1–14 min, 95% to 65%; 14–16 min, 65% to 40%; 16–18 min, 40% B; 18–18.1 min, 40% to 95%; 18.1–23 min, 95% B. The QC samples were injected at regular intervals (every 6 samples) throughout the analytical run.

Mass data acquisition was performed using a 5500 QqQ MS (AB SCIEX). Detection was performed through ESI positive and negative modes. The MS parameters were set as follows: sheath gas temperature, 350 °C; dry gas temperature, 350 °C; sheath gas flow, 11 L/min; dry gas flow, 10 L/min; capillary voltage, 4000 V or − 3500 V in positive or negative modes, respectively; nozzle voltage, 500 V; and nebulizer pressure, 30 psi. Multiple reaction monitoring (MRM) mode was used for detection. The dwell time of each ion pair was 3 ms, and the total cycle time was 1.263 s. The original MRM-based metabolomics data were analyzed using MRMAnalyzer (R) as previously reported34.

Unsupervised PCA analysis was first performed to visualize the general trends for all the samples. Then, a supervised OPLS-DA model was used to classify the samples and find the relevant variables related to the sample grouping. SIMCA-p software (version 14.1, https://www.sartorius.com) was used to perform PCA and OPLS-DA analysis. The model was validated with 200 random permutations to assess the predictive variation of the model. Variable importance in projections (VIP) scores obtained from the OPLS-DA model were used to assess the contribution of each variable to the established model. Metabolites that had VIP > 1 and p < 0.05 are identified as significantly changed metabolites. MetaboAnalyst 5.0 (https://www.metaboanalyst.ca/) was used for the functional enrichment analysis of the disturbed metabolites. All p values were corrected using Benjamini–Hochberg multiple test correction.

Total genomic DNA was isolated from feces using the QIAamp Fast DNA Stool Mini Kit (QIAGEN, Germany) following the manufacturer’s instructions. The quality and quantity of the extracted DNA were measured using a NanoDrop ND-2000 spectrophotometer (Thermo Fisher Scientific, USA). The V3-V4 region of the 16S rRNA gene were amplified with the forward primer 5′-ACTCCTACGGGAGGCAGCA-3′, and the reverse primer 5′-GGACTACHVGGGTWTCTAAT-3′. The sample-specific barcodes were added to the ends of the primers. The PCR program was set as follows: 98 °C 10 min, 25 cycles of 98 °C 15 s, 55 °C 30 s, 72 °C 30 s, and 72 °C 5 min. The products of amplification were purified and then sequenced on Illumina MiSeq platform with MiSeq Reagent Kit v3.

The sequencing data were analyzed using QIME package (version 1.9.1). High-quality paried-end reads (quality score ≥ 20) were assembled into tags by using FLASH (version 1.2.11). Tags with > 97% sequences identity were clustered into OTUs using USEARCH (version 10.0), and UCHIME (version 8.1) was utilized to identify and remove chimeric sequences. A representative sequence of each OTU was selected and subjected to BLAST to assign taxonomic classification using SILVA database (version 132).

The alpha diversity indices were calculated by QIIME (version 1.9.1). A one-way analysis of similarity (ANOSIM) was performed to determine the differences in bacterial communities among groups. The differences in Alpha diversity indexes and phyla, family, and genera relative abundances between groups were calculated by use of the Independent-sample t-test (for the normally distributed data) or Wilcoxon rank-sum test (for the non-normally distributed data). A p value < 0.05 was considered statistically significant. Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) was used to obtain relative Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway abundance information51,52,53,54.

Spearman correlation between the relative abundance of genera and the level of plasma metabolites was calculated in R software (version 3.2.1, https://www.r-project.org/) and visualized using ComplexHeatmap package (version 2.9.3, http://bioconductor.org/packages/release/bioc/html/ComplexHeatmap.html) in R software (version 3.2.1).

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Dole, M., Wilson, F. & Fife, W. Hyperbaric hydrogen therapy: A possible treatment for cancer. Science 190, 152–154. https://doi.org/10.1126/science.1166304 (1975).

Article ADS CAS PubMed Google Scholar

Ohsawa, I. et al. Hydrogen acts as a therapeutic antioxidant by selectively reducing cytotoxic oxygen radicals. Nat. Med. 13, 688–694. https://doi.org/10.1038/nm1577 (2007).

Article CAS PubMed Google Scholar

Li, S. et al. Molecular hydrogen protects against ischemia–reperfusion injury in a mouse fatty liver model via regulating HO-1 and Sirt1 expression. Sci. Rep. 9, 14019. https://doi.org/10.1038/s41598-018-32411-4 (2018).

Article ADS CAS Google Scholar

Huang, J., Liu, W. & Sun, X. Hydrogen inhalation improves mouse neurological outcomes after cerebral ischemia/reperfusion independent of anti-necroptosis. Med. Gas Res. 8, 1–5. https://doi.org/10.4103/2045-9912.229596 (2018).

Article CAS PubMed PubMed Central Google Scholar

Yoritaka, A. et al. Pilot study of H2 therapy in Parkinson’s disease: A randomized double-blind placebo-controlled trial. Mov. Disord. 28, 836–839. https://doi.org/10.1002/mds.25375 (2013).

Article CAS PubMed Google Scholar

Li, J. et al. Hydrogen-rich saline improves memory function in a rat model of amyloid-beta-induced Alzheimer’s disease by reduction of oxidative stress. Brain Res. 1328, 152–161. https://doi.org/10.1016/j.brainres.2010.02.046 (2010).

Article ADS CAS PubMed Google Scholar

Kato, R. et al. Hydrogen gas attenuates embryonic gene expression and prevents left ventricular remodeling induced by intermittent hypoxia in cardiomyopathic hamsters. Am. J. Physiol. Heart Circ. Physiol. 307, H1626–H1633. https://doi.org/10.1152/ajpheart.00228.2014 (2014).

Article CAS PubMed Google Scholar

Yu, Y. & Zheng, H. Chronic hydrogen-rich saline treatment reduces oxidative stress and attenuates left ventricular hypertrophy in spontaneous hypertensive rats. Mol. Cell. Biochem. 365, 233–242. https://doi.org/10.1007/s11010-012-1264-4 (2012).

Article CAS PubMed Google Scholar

Kajiyama, S. et al. Supplementation of hydrogen-rich water improves lipid and glucose metabolism in patients with type 2 diabetes or impaired glucose tolerance. Nutr. Res. 28, 137–143. https://doi.org/10.1016/j.nutres.2008.01.008 (2008).

Article CAS PubMed Google Scholar

Song, G. et al. Hydrogen-rich water decreases serum LDL-cholesterol levels and improves HDL function in patients with potential metabolic syndrome. J. Lipid Res. 54, 1884–1893. https://doi.org/10.1194/jlr.M036640 (2013).

Article CAS PubMed PubMed Central Google Scholar

Sada, H. et al. Peritoneal lavage with hydrogen-rich saline can be an effective and practical procedure for acute peritonitis. Surg. Today 51, 1860–1871. https://doi.org/10.1007/s00595-021-02271-z (2021).

Article CAS PubMed Google Scholar

Niu, Y. et al. Hydrogen attenuates allergic infammation by reversing energy metabolic pathway switch. Sci. Rep. 10, 1962. https://doi.org/10.1038/s41598-020-58999-0 (2020).

Article ADS CAS PubMed PubMed Central Google Scholar

Zhu, B., Cui, H. & Xu, W. Hydrogen inhibits the proliferation and migration of gastric cancer cells by modulating lncRNA MALAT1/miR-124-3p/EZH2 axis. Cancer Cell Int. 21, 70. https://doi.org/10.1186/s12935-020-01743-5 (2021).

Article CAS PubMed PubMed Central Google Scholar

Chen, J. et al. Hydrogen therapy can be used to control tumor progression and alleviate the adverse events of medications in patients with advanced non-small cell lung cancer. Med. Gas Res. 10, 75–80. https://doi.org/10.4103/2045-9912.285560 (2020).

Article CAS PubMed PubMed Central Google Scholar

Zhang, Y. et al. Anti-inflammatory effect of hydrogen-rich saline in a rat model of regional myocardial ischemia and reperfusion. Int. J. Cardiol. 148, 91–95. https://doi.org/10.1016/j.ijcard.2010.08.058 (2011).

Article PubMed Google Scholar

Cui, Y. et al. Hydrogen-rich saline attenuates neuronal ischemia–reperfusion injury by protecting mitochondrial function in rats. J. Surg. Res. 192, 564–572. https://doi.org/10.1016/j.jss.2014.05.060 (2014).

Article CAS PubMed Google Scholar

Li, H. et al. Hydrogen-rich saline protects against small-scale liver ischemia–reperfusion injury by inhibiting endoplasmic reticulum stress. Life Sci. 194, 7–14. https://doi.org/10.1016/j.lfs.2017.12.022 (2018).

Article CAS PubMed Google Scholar

Lynch, S. V. & Pedersen, O. The human intestinal microbiome in health and disease. N. Engl. J. Med. 375, 2369–2379. https://doi.org/10.1056/NEJMra1600266 (2016).

Article CAS PubMed Google Scholar

Fan, Y. & Pedersen, O. Gut microbiota in human metabolic health and disease. Nat. Rev. Microbiol. 19, 55–71. https://doi.org/10.1038/s41579-020-0433-9 (2021).

Article CAS PubMed Google Scholar

Ostojic, S. M. Hydrogen-rich water as a modulator of gut microbiota?. J. Funct. Foods 78, 104360. https://doi.org/10.1016/j.jff.2021.104360 (2021).

Article CAS Google Scholar

Xun, Z. et al. Effects of long-term hydrogen intervention on the physiological function of rats. Sci. Rep. 10, 18509. https://doi.org/10.1038/s41598-020-75492-w (2020).

Article CAS PubMed PubMed Central Google Scholar

Ara, J. et al. Hydrogen water drinking exerts antifatigue effects in chronic forced swimming mice via antioxidative and anti-inflammatory activities. Biomed. Res. Int. 2018, 2571269. https://doi.org/10.1155/2018/2571269 (2018).

Article CAS PubMed PubMed Central Google Scholar

Zheng, Y., Han, C., Chen, Z. & Zeng, Q. A study on the effect of hydrogen on the hyperuricemia in males. Chin. J. Clin. Healthc. 20, 118–121. https://doi.org/10.3969/J.issn.1672-6790.2017.02.002 (2017).

Article CAS Google Scholar

Wiener, D. & Schwartz, S. The epitranscriptome beyond m6A. Nat. Rev. Genet. 22, 119–131. https://doi.org/10.1038/s41576-020-00295-8 (2021).

Article CAS PubMed Google Scholar

He, P. C. & He, C. m6A RNA methylation: From mechanisms to therapeutic potential. EMBO J 40, e105977. https://doi.org/10.15252/embj.2020105977 (2021).

Article CAS PubMed PubMed Central Google Scholar

Ge, L., Yang, M., Yang, N., Yin, X. & Song, W. Molecular hydrogen: A preventive and therapeutic medical gas for various diseases. Oncotarget 8, 102653–102673. https://doi.org/10.18632/oncotarget.21130 (2017).

Article PubMed PubMed Central Google Scholar

Katsyuba, E., Romani, M., Hofer, D. & Auwerx, J. NAD+ homeostasis in health and disease. Nat. Metab. 2, 9–31. https://doi.org/10.1038/s42255-019-0161-5 (2020).

Article CAS PubMed Google Scholar

Crenn, P., Messing, B. & Cynober, L. Citrulline as a biomarker of intestinal failure due to enterocyte mass reduction. Clin. Nutr. 27, 328–339. https://doi.org/10.1016/j.clnu.2008.02.005 (2008).

Article CAS PubMed Google Scholar

Gong, J. et al. Serum citrulline is a simple quantitative marker for small intestinal enterocytes mass and absorption function in short bowel patients. J. Surg. Res. 127, 177–182. https://doi.org/10.1016/j.jss.2005.04.004 (2005).

Article CAS Google Scholar

Khacho, P., Wang, B. & Bergeron, R. The good and bad sides of NAAG. Adv. Pharmacol. 76, 311–349. https://doi.org/10.1016/bs.apha.2016.01.003 (2016).

Article CAS PubMed Google Scholar

Kinoshita, K., Arai, K., Kawaura, K., Hiyoshi, T. & Yamaguchi, J. Development, validation, and application of a surrogate analyte method for determining N-acetyl-l-aspartyl-l-glutamic acid levels in rat brain, plasma, and cerebrospinal fluid. J. Chromatogr. B Analyt. Technol. Biomed. Life Sci. 1003, 1–11. https://doi.org/10.1016/j.jchromb.2015.09.005 (2015).

Article CAS PubMed Google Scholar

Chen, W., Zhang, H. & Qin, S. Neuroprotective effects of molecular hydrogen: A critical review. Neurosci. Bull. 37, 389–404. https://doi.org/10.1007/s12264-020-00597-1 (2021).

Article CAS PubMed Google Scholar

Chen, L. et al. UPLC-QTOF/MS-based metabolomics reveals the protective mechanism of hydrogen on mice with ischemic stroke. Neurochem. Res. 44, 1950–1963. https://doi.org/10.1007/s11064-019-02829-x (2019).

Article CAS PubMed Google Scholar

Qiu, X. et al. Saturated hydrogen improves lipid metabolism disorders and dysbacteriosis induced by a high-fat diet. Exp. Biol. Med. (Maywood) 245, 512–521. https://doi.org/10.1177/1535370219898407 (2020).

Article CAS Google Scholar

Jin, Z. et al. Nanocapsule-mediated sustained H2 release in the gut ameliorates metabolic dysfunction-associated fatty liver disease. Biomaterials 276, 121030. https://doi.org/10.1016/j.biomaterials.2021.121030 (2021).

Article CAS PubMed Google Scholar

Liu, C. et al. Estimation of the hydrogen concentration in rat tissue using an airtight tube following the administration of hydrogen via various routes. Sci. Rep. 4, 5485. https://doi.org/10.1038/srep05485 (2014).

Article CAS PubMed PubMed Central Google Scholar

Danneskiold-Samsøe, N. B. et al. Interplay between food and gut microbiota in health and disease. Food Res. Int. 115, 23–31. https://doi.org/10.1016/j.foodres.2018.07.043 (2019).

Article CAS PubMed Google Scholar

Tanaka, Y., Kiuchi, M., Higashimura, Y., Naito, Y. & Koyama, K. The effects of ingestion of hydrogen-dissolved alkaline electrolyzed water on stool consistency and gut microbiota: A double-blind randomized trial. Med. Gas Res. 11, 138–144. https://doi.org/10.4103/2045-9912.318858 (2021).

Article CAS PubMed PubMed Central Google Scholar

Lukić, I. et al. Antidepressants affect gut microbiota and Ruminococcus flavefaciens is able to abolish their effects on depressive-like behavior. Transl. Psychiatry 9, 133. https://doi.org/10.1038/s41398-019-0466-x (2019).

Article PubMed PubMed Central Google Scholar

Zhang, Y. et al. Effects of hydrogen-rich water on depressive-like behavior in mice. Sci. Rep. 6, 23742. https://doi.org/10.1038/srep23742 (2016).

Article ADS CAS PubMed PubMed Central Google Scholar

Guo, P., Zhang, K., Ma, X. & He, P. Clostridium species as probiotics: Potentials and challenges. J. Anim. Sci. Biotechnol. 11, 24. https://doi.org/10.1186/s40104-019-0402-1 (2020).

Article PubMed PubMed Central Google Scholar

van der Hee, B. & Wells, J. M. Microbial regulation of host physiology by short-chain fatty acids. Trends Microbiol. 29, 700–712. https://doi.org/10.1016/j.tim.2021.02.001 (2021).

Article CAS PubMed Google Scholar

Qin, J. et al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature 490, 55–60. https://doi.org/10.1038/nature11450 (2012).

Article ADS CAS PubMed Google Scholar

Kamimura, N., Nishimaki, K., Ohsawa, I. & Ohta, S. Molecular hydrogen improves obesity and diabetes by inducing hepatic FGF21 and stimulating energy metabolism in db/db mice. Obesity 19, 1396–1403. https://doi.org/10.1038/oby.2011.6 (2011).

Article CAS PubMed Google Scholar

Liu, X. et al. Blautia—A new functional genus with potential probiotic properties?. Gut Microbes 13, 1875796. https://doi.org/10.1080/19490976.2021.1875796 (2021).

Article CAS PubMed Central Google Scholar

Yang, H. et al. Characteristics of the intestinal microorganisms in middle-aged and elderly patients: Effects of smoking. ACS Omega https://doi.org/10.1021/acsomega.1c02120 (2022).

Article PubMed PubMed Central Google Scholar

Morley, J. E. Decreased food intake with aging. J. Gerontol. A Biol. Sci. Med. Sci. 56, 81–88. https://doi.org/10.1093/gerona/56.suppl_2.81 (2001).

Article PubMed Google Scholar

Hyspler, R., Ticha, A., Schierbeek, H., Galkin, A. & Zadak, Z. The evaluation and quantification of dihydrogen metabolism using deuterium isotope in rats. PLoS One 10, e0130687. https://doi.org/10.1371/journal.pone.0130687 (2015).

Article CAS PubMed PubMed Central Google Scholar

Zhang, X. et al. Hydrogen evolution and absorption phenomena in the plasma membrane of Vigna radiata and Capsicum annuum. J. Plant Growth Regul. https://doi.org/10.1007/s00344-021-10544-7 (2022).

Article Google Scholar

Zhang X., et al. Mitochondria in higher plants possess H2 evolving activity which is closely related to complex I. arXiv:2001.02132 (2020).

Langille, M. et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31, 814–821. https://doi.org/10.1038/nbt.2676 (2013).

Article CAS PubMed PubMed Central Google Scholar

Kanehisa, M., Furumichi, M., Sato, Y., Ishiguro-Watanabe, M. & Tanabe, M. KEGG: Integrating viruses and cellular organisms. Nucleic Acids Res. 19, D545–D551. https://doi.org/10.1093/nar/gkaa970 (2021).

Article CAS Google Scholar

Kanehisa, M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 28, 1947–1951. https://doi.org/10.1002/pro.3715 (2019).

Article CAS PubMed PubMed Central Google Scholar

Kanehisa, M. & Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30. https://doi.org/10.1093/nar/28.1.27 (2000).

Article CAS PubMed PubMed Central Google Scholar

Download references

This work was supported by Military Logistics Key Open Research Projects (BHJ17L018).

These authors contributed equally: Fei Xie, Xue Jiang and Yang Yi.

Faculty of Environment and Life, Beijing University of Technology, No. 100, Pingleyuan, Chaoyang District, Beijing, 100124, China

Fei Xie, Xue Jiang, Yang Yi, Zi-Jia Liu, Chen Ma, Jin He, Zhi-ming Xun, Meng Wang, Meng-yu Liu, Yao Mawulikplimi Adzavon, Peng-xiang Zhao & Xue-mei Ma

Beijing Molecular Hydrogen Research Center, Beijing, 100124, China

Fei Xie, Xue Jiang, Yang Yi, Zi-Jia Liu, Chen Ma, Jin He, Zhi-ming Xun, Meng Wang, Meng-yu Liu, Yao Mawulikplimi Adzavon, Peng-xiang Zhao & Xue-mei Ma

You can also search for this author in PubMed Google Scholar

X.-M. and F.X. designed the overall project. F.X., X.J., Y.Y., Z.-j.L., and C.M. performed the animal experiments. J.H. and Z.-X. performed the LC–MS based pseudotargeted metabolomics analysis, M.W. and M.-L. performed the 16S rRNA gene sequencing, F.X. analyzed the data and wrote the manuscript, Y.-A., P.-Z., and X.-M. revised the manuscript. All authors reviewed the manuscript.

Correspondence to Xue-mei Ma.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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 http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Xie, F., Jiang, X., Yi, Y. et al. Different effects of hydrogen-rich water intake and hydrogen gas inhalation on gut microbiome and plasma metabolites of rats in health status. Sci Rep 12, 7231 (2022). https://doi.org/10.1038/s41598-022-11091-1

Download citation

Received: 01 August 2021

Accepted: 31 March 2022

Published: 04 May 2022

DOI: https://doi.org/10.1038/s41598-022-11091-1

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

News