Quanzheng Zhang1†
Rui Mao- 1Department of Critical Care Medicine, The Third People’s Hospital of Chengdu, The Affiliated Hospital of Southwest Jiaotong University, Chengdu, Sichuan, China
- 2West China Fourth Hospital, Sichuan University, Chengdu, China
- 3Department of Dermatology, Xiangya Hospital, Central South University, Changsha, China
Background: Observational studies have provided evidence of a close association between gut microbiota and the progression of chronic hepatitis B (CHB). However, establishing a causal relationship between gut microbiota and CHB remains a subject of investigation.
Methods: Genome-wide association study (GWAS) summary data of gut microbiota came from the MiBioGen consortium, while the GWAS summary data of CHB came from the Medical Research Council Integrative Epidemiology Unit (IEU) Open GWAS project. Based on the maximum likelihood (ML), Mendelian randomization (MR)-Egger regression, inverse variance weighted (IVW), MR Pleiotropy RESidual Sum and Outlier (MR-PRESSO), and weighted-mode and weighted-median methods, we conducted a bidirectional, two-sample, MR analysis to explore the causal relationship between the gut microbiota and CHB. Additionally, we evaluated the genetic associations between individual gut microbes and CHB using the Linkage disequilibrium score regression (LDSC) program.
Results: According to the IVW method estimates, genetically predicted class Alphaproteobacteria (odds ratio [OR] = 0.57; 95% confidence interval [CI], 0.34–0.96; false discovery rate [FDR] = 0.046), genus Family XIII AD3011 group (OR = 0.60; 95% CI, 0.39–0.91; FDR = 0.026), genus Prevotella 7 (OR = 0.73; 95% CI, 0.56–0.94; FDR = 0.022) exhibited a protective effect against CHB. On the other hand, family Family XIII (OR = 1.79; 95% CI, 1.03–3.12; FDR = 0.061), genus Eggerthella group (OR = 1.34; 95% CI, 1.04–1.74; FDR = 0.043), genus Eubacterium ventriosum group (OR = 1.59; 95% CI, 1.01–2.51; FDR = 0.056), genus Holdemania (OR = 1.35; 95% CI, 1.00–1.82; FDR = 0.049), and genus Ruminococcus gauvreauii group (OR = 1.69; 95% CI, 1.10–2.61; FDR = 0.076) were associated with an increased risk of CHB. The results from LDSC also indicated a significant genetic correlation between most of the aforementioned gut microbiota and CHB. Our reverse MR analysis demonstrated no causal relationship between genetically predicted CHB and gut microbiota, and we observed no significant horizontal pleiotropy or heterogeneity of instrumental variables (IVs).
Conclusion: In this study, we identified three types of gut microbiota with a protective effect on CHB and five types with an adverse impact on CHB. We postulate that this information will facilitate the clinical prevention and treatment of CHB through fecal microbiota transplantation.
Introduction
Nearly two billion people worldwide are affected by hepatitis B virus (HBV) infection, with 10–20% of them eventually developing decompensated liver cirrhosis or liver cancer, posing significant health risks (Trépo et al., 2014). The persistent nature of chronic hepatitis B (CHB) stems from a complex interplay between host immunity, the virus itself, environmental factors, and other variables, with impaired host immune function being a primary factor.
The gut microbiota refers to the complex community of microorganisms, including bacteria, viruses, fungi, and archaea, that reside in the digestive tracts of humans. This intricate ecosystem is comprised of trillions of microbes, with thousands of different species cohabitating and interacting within the host. The human gut microbiome is not only involved in essential metabolic processes but also plays a crucial role in immune regulation and maintenance of gut integrity (Chen et al., 2021; Ahmad et al., 2022; Rasheed et al., 2022). Current research suggests a potential link between gut microbiota and host immune function. Gut microbiota can directly or indirectly influence host immunity through the gut-liver axis (Compare et al., 2012; Milosevic et al., 2019). Under normal conditions, the interaction between gut microbiota and Toll-like receptors (TLR) on intestinal epithelial cells and immune cells helps maintain immune system homeostasis. In addition, gut microbiota maintains a dynamic balance with the liver by regulating metabolism and immune response and mediates immune and inflammatory responses by participating in bile acid metabolism (Yiu et al., 2017; Wang et al., 2020). Imbalances in gut microbiota may occur at all stages of CHB and may contribute to the chronicity of infection and the progression of the disease. Conversely, by regulating gut microbes, it may provide new avenues for the treatment of hepatitis B. Some studies are exploring how to use probiotics, prebiotics and other methods to improve the clinical outcome of hepatitis (Chen et al., 2011; Khoruts and Sadowsky, 2016; Morrison and Preston, 2016; Plaza-Diaz et al., 2019).
However, the causal relationship between gut microbiota and CHB remains uncertain. The comprehensive analysis of 86 independent studies by Sun et al. (2021). revealed that the changes in intestinal microbial metabolites abrogated the intestinal microbiota-bile acid-host axis, influencing the progression of CHB. Yet, several cross-sectional observations studies suggest that the occurrence and development of CHB can lead to alterations in gut microbiota (Chen et al., 2020; Zeng et al., 2020; Joo et al., 2021; Shen et al., 2023). Most earlier analyses were designed as case–control studies, making it challenging to ascertain exposure times and definitive outcomes. Additionally, in observational studies, the association between the gut microbiota and CHB is easily affected by confounding factors such as dietary patterns, environment, age, and lifestyle; and it is challenging to control these factors via observational studies effectively (Rinninella et al., 2019). These conditions limit the inference of a causal relationship between the gut microbiota and CHB.
Research design with MR follows the Mendelian inheritance law where “parental alleles are randomly assigned to offspring.” If genotype determines phenotype, the genotype is associated with a particular disease through the phenotype. Thus the association between CHB and gut microbiota can be inferred using the genotype as an instrumental variable. As an extension of the MR approach, bidirectional MR can be used to determine the direction of causality between two related phenotypes.
In the present study, we conducted a bidirectional, two-sample MR analysis adopting summary statistics of large GWAS from the MiBioGen consortium and the IEU Open GWAS project to assess the causal association between gut microbiota and CHB.
Methods
Data sources
We obtained the genetic variation in gut microbiotas from the largest genome-wide meta-analysis of gut microbiota to date conducted by the MiBioGen consortium (Kurilshikov et al., 2021; MiBioGen consortium, 2021). This study comprised 18,340 individuals from 24 cohorts (including the United States, Canada, Israel, South Korea, Germany, Denmark, the Netherlands, Belgium, Sweden, Finland, and the United Kingdom), targeting the variable regions V4, V3-V4, and V1-V2 of the 16S rRNA gene to analyse the microbial composition and classify the specific microbes using direct classification. We analysed microbiota quantitative trait loci (mbQTL) of microflora to identify the host genetic variation of genetic loci related to the abundance level of bacterial taxa in gut microflora. The genus was the lowest taxonomic level described in our study, and 211 taxa (131 genera, 35 families, 20 orders, 16 classes, and nine phyla) with an average abundance of more than 1% were then identified. The GWAS summary data (UK biobank, FinnGen consortium, and Biobank Japan) for CHB came from the data released by Sakaue et al. (2021), and these were then sorted and stored by the GWAS catalog (GWAS ID: GCST90018804, Trait name: Chronic hepatitis B) (IEU OpenGWAS project, 2021). The GWAS data included 19,082,360 single-nucleotide polymorphisms (SNPs) and 523,707 samples, including 2,379 in the case group and 521,328 in the control group. Age, sex, top 10 principal components, and genotyping batches were corrected during analysis (Sakaue et al., 2021).
Instrumental variables
We used the following four criteria to select instrumental variables. First, the SNP-phenotype association level was required to reach the locus-wide significance threshold (p < 1 × 10−8). Unfortunately, only a small number of intestinal microbiota were selected as IV after we selected SNP. In order to explore the relationship between CHB and intestinal microflora to obtain more comprehensive results, we used the second threshold to identify a lower level of SNP- significance than the locus (p < 1 × 10−5) and selected them as the second IV set to find more potential causal associations. Second, SNPs with secondary (minor) allelic frequencies (MAFs) ≤ 0.01 were removed. Third, we only retained the lowest p-value among the SNPs with an r2 < 0.001 (with a clumping window size = 10,000 kb). When performing LD clumping on SNP data, we used EUR as a reference panel. Finally, we used the allelic frequency information to infer the forward strand alleles when palindromic SNPs were uncovered.
Statistical analysis
The design and flow of the entire study are presented in Figure 1. In the present study, we evaluated the causal relationship between gut microflora and CHB using various methods, including maximum likelihood (ML), inverse variance weighted (IVW), MR Pleiotropy RESidual Sum and Outlier (MR-PRESSO), weighted-mode, MR-Egger regression, and weighted-median methods. IVW is characterized by ignoring the existence of an intercept term in regression and using the reciprocal of the outcome variance (se2) as the weight used for fitting. In the absence of heterogeneity and horizontal pleiotropy, IVWs constitute the most reliable estimates (Burgess et al., 2016). The ML method is similar to the IVW method, and under the assumption that there was no heterogeneity and that horizontal multiplicity was established, the evaluation result for ML showed almost no deviation, and its standard error was less than the IVW (Pierce and Burgess, 2013). MR-Egger regression is based on the assumption that instrument strength is independent of direct effect (InSIDE), which makes it possible to evaluate the existence of poly-efficacity with intercept terms. An intercept term of zero indicated that horizontal pleiotropy did not exist, and the MR-Egger regression result was consistent with IVW (Bowden et al., 2015). In addition, MR-Egger regression was also used to ascertain whether or not there was horizontal poly-tropism. When more than 50% of the instrumental variables are invalid, the weighted-median method can estimate causality more correctly. If the InSIDE hypothesis is violated, the weighted-model estimation method exhibits greater detection ability with regard to causal effect, less deviation, and a lower class-I error rate than with the MR-Egger regression (Hartwig et al., 2017). In addition, significant outliers were detected using MR-PRESSO (Verbanck et al., 2018) tests and MR-Egger regression, and horizontal pleiotropic effects were corrected by removing outliers. The global test also evaluates whether horizontal pleiotropy among all instruments is present (Verbanck et al., 2018). In addition, we further examined the heterogeneity among all SNPs using Cochran’s Q-test statistics. We identified potentially heterogeneous SNPs by conducting a “leave-one-out” analysis and excluding each instrumental SNP in turn. Finally, we also performed reverse MR analysis of CHB and gut microflora. The methods and settings we adopted were consistent with those of forward MR.
Figure 1. The design and workflow of the entire study.
The F statistic was calculated using the formula , where R2 was the proportion of the trait variance explained by the SNP, k was the number of IVs, and N was the sample size of the GWAS regarding the SNP for the particular trait; this was employed to quantify the strength of the instrument and a value >10 was considered sufficient. We used the TwoSampleMR R package to calculate the coefficient of determination (R2) regarding exposure to genetic variants. The R2 value was estimated using the formula , where EAF was the effect allele frequency (EAF) of the SNP, SD was the standard deviation, and β was the estimated effect size of the SNP on the trait. Based on the online MR-power calculation tool1 (Burgess, 2014), we calculated the power of the MR estimates. Furthermore, we calculated the genetic association between each gut microbe and CHB using the LDSC program.
False discovery rate (FDR) correction was conducted by applying the fdrtools procedure, with a false discovery rate of q-value < 0.1 (Strimmer, 2008). All analyses were performed using R (version 4.2.2), and MR analysis was based on TwoSampleMR (Hemani et al., 2017) (version 0.5.6), meta (version 6.2.1) (Balduzzi et al., 2019), and MR-PRESSO (version 1.0) packages. The heatmap was drawn with the ImageGP tool (Chen et al., 2022).
Results
SNP selection
Based on the screening criteria of the IVs, we initially identified 1, 1, 3, 6, and 19 SNPs associated with gut microbiota at the phylum, class, order, family, and genus levels, respectively, at a significance level of p < 5 × 10−8. Details of each SNP, including specific chromosome, EAF value, and R2, are provided in Supplementary Table S1. However, the number of IVs obtained was not sufficient for a robust MR analysis. Therefore, under the screening of another threshold p < 1 × 10−5, we get 125, 224, 280, 503, and 1,667 SNPs associated with gut microbiota at the phylum, class, order, family, and genus levels. The details of SNP in the exposure and four outcome variables are shown in Supplementary Table S2. After excluding duplicated SNPs across different microorganisms, we got 21 SNP at a threshold of p < 5 × 10−8 and 1,676 SNP at a threshold of p < 1 × 10−5.
Causal effects of gut microbiota on the development of CHB
Based on the IVs set filtered by a threshold of p < 5 × 10−8, We used gut microbiota as exposure and CHB as outcomes for MR analysis. However, due to the lack of IVs, the results of MR analysis can not be obtained. Therefore, we use the IVs filtered by the threshold of p < 1 × 10−5 to perform the following MR analysis (Kurilshikov et al., 2021; Li et al., 2022; Long et al., 2023).
Four gut microbiotas were excluded because the available IVs were less than 3. The causal effects of the remaining 207 gut microbiotas on CHB are shown in Supplementary Table S3. Figure 2 illustrates 11 gut microbiotas that exert significant causal effects on CHB. Briefly, the estimate of IVW suggested that genetically predicted class Alphaproteobacteria (OR = 0.57; 95% CI, 0.34–0.96; p = 0.033; FDR = 0.046), genus Family XIII AD3011 group (OR = 0.60; 95% CI, 0.39–0.91; p = 0.017; FDR = 0.026), genus Prevotella 7 (OR = 0.73; 95% CI, 0.56–0.94; p = 0.015; FDR = 0.022), unknown genus id.2001 (OR = 0.71; 95% CI, 0.52–0.97; p = 0.010; FDR = 0.029), unknown genus id.2755 (OR = 0.70; 95% CI, 0.53–0.93; p = 0.013; FDR = 0.023), and unknown genus id.959 (OR = 0.76; 95% CI, 0.58–0.99; p = 0.047; FDR = 0.065) all exhibited a protective causal relationship with CHB. Conversely, genetically predicted family Family XIII (OR = 1.79; 95% CI, 1.03–3.12; p = 0.039; FDR = 0.061), genus Eggerthella (OR = 1.34; 95% CI, 1.04–1.74; p = 0.023; FDR = 0.043), genus Eubacterium ventriosum group (OR = 1.59; 95% CI, 1.01–2.51; p = 0.046; FDR = 0.056), genus Holdemania (OR = 1.35; 95% CI, 1.00–1.82; p = 0.040; FDR = 0.049), and genus Ruminococcus gauvreauii group (OR = 1.69; 95% CI, 1.10–2.61; p = 0.017; FDR = 0.076) exhibited a risk effect on CHB.
Figure 2. MR analysis of the causal relationship between genetically predicted gut microbiotas and CHB.
As shown in Supplementary Table S3, we noted a 0.35–8.32% contribution to the total variation (R2 values) of the 207 gut microbiotas, and the F values ranged from 10.81 to 160.49—excluding the possibility of weak genetic tool variables. Cochran’s IVW Q test revealed no significant heterogeneity among these IVs (Supplementary Table S4). Additionally, except for the genus Eubacterium ventriosum group (p = 0.034), MR-Egger regression intercept analysis showed no significant horizontal pleiotropy among the other 206 gut microbiotas (Supplementary Table S5). Further MR-PRESSO analysis did not reveal horizontal pleiotropy between the genus Eubacterium ventriosum group (Supplementary Table S6, global test p-value = 0.991) and CHB. Notably, in the resulting diagram of the leave-one-out method analysis, we observed abnormal IVs in family Family XIII, genus Eubacterium ventriosum group, genus Holdemania, and unknown genus id.959 (Figure 3). However, except for the genus Ruminococcus gauvreauii group, MR-PRESSO analysis did not reveal any significant outliers (global test p > 0.05, Supplementary Table S6). In addition, when rs431418 was excluded, outlier-corrected MR-PRESSO analysis did not disclose any horizontal pleiotropy in the genus Ruminococcus gauvreauii group (global test p = 0.122, Supplementary Table S6). Therefore, there was insufficient evidence to demonstrate horizontal pleiotropy between these gut microbiotas and CHB.
Figure 3. “Leave-one-out” plots for the causal association between genetically predicted gut microbiotas and CHB. These gut microbiotas included class Alphaproteobacteria (A), family Family XIII (B), genus Eggerthella (C), genus Eubacterium ventriosum group (D), genus Prevotella 7 (E), genus Family XIII AD3011 group (F), Ruminococcus gauvreauii group (G), genus Holdemania (H), unknown genus id.2001 (I), unknown genus id.2755 (J), and unknown genus id.959 (K).
In reverse MR analysis, we identified 11 SNPs (p < 1 × 10−5; R2 < 0.001) significantly associated with CHB, explaining 45.5% of the total variation (summaries and details for each SNP are presented in Supplementary Table S7). The results of reverse MR analysis did not reveal any significant causal relationship between CHB and the 11 gut microbiotas mentioned above (Supplementary Table S8). Cochran’s Q statistics indicated no notable heterogeneity (Supplementary Table S9, p > 0.05), and the MR-Egger intercept test (Supplementary Table S10) and MR-PRESSO global test (Supplementary Table S11) suggested no evidence of pleiotropy.
The genetic correlation between each gut microbe and CHB is presented in Supplementary Table S12. The results suggest that there is a significant genetic correlation between gut microbiota such as genus Family XIII AD3011, genus Holdemania, class Alphaproteobacteria, genus Eubacterium ventriosum, genus Eggerthella, genus Ruminococcus gauvreauii, family Family XIII and chronic hepatitis B.
Finally, we employed ImageGP to illustrate the relative abundance of bacteria in individuals exhibiting a non-zero abundance of the 11 gut microbiota significantly associated with chronic hepatitis B, along with the average read count of these 11 gut microbiota across the study participants (see Supplementary Figure S1). Evidently, a substantial majority of these 11 gut microbiota are observed in over 50% of the population. The fecal transplantation of any of these 11 gut microbiota could thus ensure its applicability to a broad segment of the population.
Discussion
Despite the significant progress made in reducing annual hepatitis B virus (HBV) infections through the use of vaccines and antiviral drugs, HBV remains a substantial economic burden and health threat to society (Lok et al., 2016; Tan et al., 2021). To effectively treat patients with chronic liver disease and gain further insights into its pathogenesis, it is crucial to explore appropriate interventions and treatment strategies.
The relationship between gut microbiota and CHB represents an intricate intersection of microbial ecology, host metabolism, and hepatic immunology. Emerging evidence points to alterations in the gut microbiome, known as dysbiosis, playing a significant role in the progression of CHB. The gut-liver axis serves as a bidirectional communication channel where microbial metabolites and pathogen-associated molecular patterns (PAMPs) can affect liver inflammation and fibrosis. In CHB patients, dysbiosis may lead to an increased permeability of the gut barrier, allowing the translocation of bacterial products into the liver and subsequently triggering inflammatory responses. This can contribute to liver damage and exacerbate the disease progression. Moreover, the complex interplay of gut microbiota with antiviral immune responses may influence the efficacy of antiviral treatments in CHB. Some researchers are investigating therapeutic interventions, such as the use of probiotics or prebiotics, to modulate the gut microbiome and potentially improve clinical outcomes for CHB patients. However, the mechanistic understanding of how gut microbiota influences CHB is still in its infancy, and further rigorous investigations are warranted to translate these findings into clinical practice (Chen et al., 2011; Lu et al., 2011; Wang et al., 2017).
In this study, we conducted the first-ever investigation into the causal relationship between gut microbiota and CHB using the summary statistics of the largest genome-wide meta-analysis of CHB and gut microbiotas conducted by the MiBioGen consortium. MR is less prone to confounding factors because germline genetic variation is randomly assigned during meiosis and therefore reflects exposure without being affected by reverse causality. The results of our bidirectional, two-sample MR study suggested that genetic predisposition to family Family XIII (OR = 1.79), genus Eggerthella (OR = 1.34), genus Eubacterium ventriosum group (OR = 1.59), genus Holdemania (OR = 1.35), and genus Ruminococcus gauvreauii group (OR = 1.69) was associated with a 79, 34, 59, 35, and 69% increased risk of CHB, respectively. Genetically predicted class Alphaproteobacteria (OR = 0.57), genus Family XIII AD3011 group (OR = 0.60), and genus Prevotella 7 (OR = 0.73) were also associated with a 43, 40, and 27% decreased risk of CHB, respectively. However, reverse MR analysis displayed no significant causal relationship between CHB and gut microbiota. In addition, various validation methods were used to ensure the reliability and robustness of our findings.
The genus Ruminococcus gauvreauii group belongs to the family Ruminococcaceae, which is an integral part of the Firmicutes phylum. Members of this genus are widely recognized as critical constituents of the gut microbiota in both humans and other animals (Chen et al., 2021). It has been hypothesized that these bacteria may influence liver inflammation and fibrosis through the modulation of short-chain fatty acids production, bile acid metabolism, and interaction with Toll-like receptors, thus impacting the overall gut-liver axis. Komiyama et al. (2021) and Shen et al. (2023) demonstrated that the intestinal microbial abundance of the genus Ruminococcus gauvreauii group increased significantly in patients with hepatitis B and C, and this bacterial group was also used as a biomarker of liver cancer caused by these two pathogens. This result was similar to our findings in that the Ruminococcus gauvreauii group increased the risk of CHB.
Several observational studies and a meta-analysis showed that genus Eggerthella was closely related to colitis, cancer, and mental health disorders (Nikolova et al., 2021; Alexander et al., 2022; Bai et al., 2022); and that it was able to induce Th17 activation by relieving the inhibition of Th17 transcription factor Ror γ t, a cell-and-antigen-independent mechanism (Alexander et al., 2022). Several mechanistic studies have shown that Th-17 cells and their secreted cytokine IL-17 promote the progression and long-term survival of CHB (Ge et al., 2010; Shen et al., 2023). In the current study, Eggerthella significantly increased the risk of CHB, most likely also through the Th-17 pathway.
The genus Eubacterium ventriosum group and genus Holdemania are both part of the Firmicutes phylum and are found within the human gut microbiome. The Eubacterium ventriosum group is recognized for its anaerobic properties and ability to ferment carbohydrates, while Holdemania is a less characterized genus known for its presence in various environments, including the gastrointestinal tract. In the context of hepatitis, particularly chronic liver diseases, these genera might be involved in modulating hepatic inflammation. Many observational studies have also revealed increased intestinal abundances of the genus Eubacterium ventriosum group and genus Holdemania in patients with various inflammatory diseases (Jang et al., 2020; Moreno-Arrones et al., 2020; Radjabzadeh et al., 2022; Ye et al., 2022). Alterations in their abundance could impact gut barrier integrity, allowing the translocation of microbial products into the liver, thereby triggering inflammatory responses. Though the exact mechanisms remain incompletely understood, some studies have indicated that shifts in the composition of gut microbiota, including these genera, might be correlated with liver disease progression and could influence the host’s immune responses to hepatic infections (Bajaj et al., 2014; Ren et al., 2019). These studies further confirmed our results of the Eubacterium ventriosum group and Holdemania as significantly increasing the risk of CHB.
The class Alphaproteobacteria includes Gram-negative, free-living, symbiotic and obligate intracellular bacteria and important plant, animal, and human pathogens (Hallez et al., 2017). Zhang et al. (2019) discerned that class Alphaproteobacteria was significantly lower in patients with chronic and acute liver failure associated with hepatitis B than in the control group. Similarly, we ascertained that class Alphaproteobacteria reflected a protective causal relationship with respect to CHB.
The genus Prevotella 7, part of the Bacteroidetes phylum, is a group of gram-negative bacteria that are widely recognized as essential constituents of the human gut microbiota. They are involved in the fermentation of polysaccharides, thereby contributing to producing short-chain fatty acids (SCFAs) such as butyrate, which plays a crucial role in maintaining gut barrier integrity and modulating the host’s immune response. In addition, it has been reported that Prevotella 7 participates in various inflammatory pathologic processes by producing redox proteins or increasing host resistance (Hofer, 2014). Observational studies by Juan et al. also showed that Prevotella 7 was associated with a reduced risk of gastric cancer (Wu et al., 2018). In the context of hepatitis, Prevotella 7 might have protective effects. This protective effect is speculated to arise from its potential role in enhancing intestinal barrier function, which may limit the translocation of microbial metabolites and prevent their pro-inflammatory effects in the liver. Moreover, the immunomodulatory effect exerted by SCFAs can potentially attenuate hepatic inflammation, thereby preventing the progression of hepatitis (Shen et al., 2017). In our research, Prevotella 7 was the protective factor against CHB and may thus reduce the risk of CHB by regulating the inflammatory and immune responses of the host. While promising, these observations require further investigation to fully elucidate the role of Prevotella 7 in liver health.
The genus Family XIII AD3011 group represents a less well-characterized classification within the Firmicutes phylum. Research by Lin et al. (2022) showed that the intestinal abundance of the genus Family XIII AD3011 group in psoriatic arthritis was significantly lower than that in the normal group. However, there is a paucity of research on the relationship between the genus Family XIII AD3011 group and CHB. Importantly, our study revealed for the first time that the genus Family XIII AD3011 group exhibited a protective causal relationship concerning CHB.
Strengths and limitations
This study possessed several advantages. Genetic variation in gut microbiota and CHB were obtained from the most extensive available GWAS meta-analysis of the MiBioGen consortium and GWAS catalog, ensuring the strength of the instruments applied in the MR analyses. The GWAS data we used encompassed many races worldwide, including those in Europe, Asia, and Australia, ensuring the universal applicability of our findings.
Some limitations to our investigation should also be acknowledged. First, the data from gut microbiota and CHB included both Europeans and Asians, where a population stratification bias may exist. However, at present, there is no GWAS summary data of large gut microbiota of the same race, and most of the GWAS data of CHB are from East Asian people (Trépo et al., 2014) only the GWAS summary data of CHB used in this study includes both European and East Asian races. Thus, although this partial bias was inevitable, we employed as many analytical methods as possible to eliminate this potential error. Second, there was no assessment of the existence of survival bias. In the future, separate MR analyses should be conducted for different age groups, and the consistency of the results should be assessed. The SNPs used in the analysis did not reach the traditional GWAS significance threshold (i.e., p < 5 × 10−8); however, we need to include additional genetic variation as IV for sensitivity analysis and the accurate detection of horizontal polymorphism. When the p value was set at 5 × 10−8 or 1 × 10−6, only one SNP or none were available for each microbe. For this, we used FDR correction to restrict the possibility of false positives. There are many high-level studies screened this way (Kurilshikov et al., 2021; Li et al., 2022; Long et al., 2023). Finally, the bacterial groups were analysed solely at the order or family level. Implementing additional advanced macro genome sequencing analysis in the future will thus enable the results to be more specific and accurate.
Conclusion
Through this two-sample MR study, we established a causal link between gut microbiota and CHB, identifying specific gut microbiotas that either protect or exacerbate CHB risk. This information could potentially contribute to the clinical prevention and treatment of CHB through fecal microbiota transplantation.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Ethics statement
The data we used were from published studies approved by the appropriate ethics committees, and no further ethical approval was required for this study.
Author contributions
QZ: conceptualization, methodology, software, investigation, visualization, and writing of the revised draft. XZ and RM: data provision and methodology. JZ: conceptualization, writing-revised draft, and funding acquisition. CZ: conceptualization, writing-original draft, funding acquisition, and project supervision. All authors contributed to the article and approved the submitted version.
Acknowledgments
The authors wish to express our gratitude to all researchers involved herein for providing their GWAS data. The authors also acknowledge and thank the participants and researchers of the IEU open GWAS project and the MiBioGen consortium. The authors thank LetPub (www.letpub.com) for its linguistic assistance while preparing this manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1243811/full#supplementary-material
Abbreviations
MR, Mendelian randomization; ML, Maximum likelihood; IVW, Inverse variance weighted; MR-PRESSO, MR Pleiotropy RESidual Sum and Outlier; GWAS, Genome-wide association study; CHB, Chronic hepatitis B; SNP, Single nucleotide polymorphism.
Footnotes
References
Ahmad, S., Yousaf, M. S., Tahir, S. K., Mustafa, G., Majeed, K. A., Rashid, M. A., et al. (2022). Effects of co-supplementation of β-Galacto-oligosaccharides and methionine on production performance, blood metabolites, and gut Histomorphometry in broilers. Pak. Vet. J. 42, 173–178. doi: 10.29261/pakvetj/2022.031
Alexander, M., Ang, Q. Y., Nayak, R. R., Bustion, A. E., Sandy, M., Zhang, B., et al. (2022). Human gut bacterial metabolism drives Th17 activation and colitis. Cell Host Microbe 30, 17–30.e9. doi: 10.1016/j.chom.2021.11.001
Bai, X., Wei, H., Liu, W., Coker, O. O., Gou, H., Liu, C., et al. (2022). Cigarette smoke promotes colorectal cancer through modulation of gut microbiota and related metabolites. Gut 71, 2439–2450. doi: 10.1136/gutjnl-2021-325021
Bajaj, J. S., Heuman, D. M., Hylemon, P. B., Sanyal, A. J., White, M. B., Monteith, P., et al. (2014). Altered profile of human gut microbiome is associated with cirrhosis and its complications. J. Hepatol. 60, 940–947. doi: 10.1016/j.jhep.2013.12.019
Balduzzi, S., Rücker, G., and Schwarzer, G. (2019). How to perform a meta-analysis with R: a practical tutorial. Evid. Based Ment. Health 22, 153–160. doi: 10.1136/ebmental-2019-300117
Pierce, B. L., and Burgess, S. (2013). Efficient design for Mendelian randomization studies: subsample and 2-sample instrumental variable estimators. Am. J. Epidemiol. 178, 1177–1184. doi: 10.1093/aje/kwt084
Bowden, J., Davey Smith, G., and Burgess, S. (2015). Mendelian randomization with invalid instruments: effect estimation and bias detection through egger regression. Int. J. Epidemiol. 44, 512–525. doi: 10.1093/ije/dyv080
Burgess, S. (2014). Sample size and power calculations in Mendelian randomization with a single instrumental variable and a binary outcome. Int. J. Epidemiol. 43, 922–929. doi: 10.1093/ije/dyu005
Burgess, S., Dudbridge, F., and Thompson, S. G. (2016). Combining information on multiple instrumental variables in Mendelian randomization: comparison of allele score and summarized data methods. Stat. Med. 35, 1880–1906. doi: 10.1002/sim.6835
Chen, T., Liu, Y.-X., and Huang, L. (2022). ImageGP: an easy-to-use data visualization web server for scientific researchers. iMeta 1:e5. doi: 10.1002/imt2.5
Chen, Z., Xie, Y., Zhou, F., Zhang, B., Wu, J., Yang, L., et al. (2020). Featured gut microbiomes associated with the progression of chronic hepatitis B disease. Front. Microbiol. 11:383. doi: 10.3389/fmicb.2020.00383
Chen, Y., Yang, F., Lu, H., Wang, B., Chen, Y., Lei, D., et al. (2011). Characterization of fecal microbial communities in patients with liver cirrhosis. Hepatology 54, 562–572. doi: 10.1002/hep.24423
Chen, D., Zhong, G., Su, H., Rahman, M. U., Chen, K., Tang, J., et al. (2021). Physiological variation in ruminal microbiota under altered energy levels in starter ration of suckling angus calves. Pak. Vet. J. 41, 409–413. doi: 10.1371/journal.pone.0284509
Compare, D., Coccoli, P., Rocco, A., Nardone, O. M., De Maria, S., Cartenì, M., et al. (2012). Gut--liver axis: the impact of gut microbiota on non alcoholic fatty liver disease. Nutr. Metab. Cardiovasc. Dis. 22, 471–476. doi: 10.1016/j.numecd.2012.02.007
Ge, J., Wang, K., Meng, Q.-H., Qi, Z.-X., Meng, F.-L., and Fan, Y.-C. (2010). Implication of Th17 and Th1 cells in patients with chronic active hepatitis B. J. Clin. Immunol. 30, 60–67. doi: 10.1007/s10875-009-9328-2
Hallez, R., Delaby, M., Sanselicio, S., and Viollier, P. H. (2017). Hit the right spots: cell cycle control by phosphorylated guanosines in alphaproteobacteria. Nat. Rev. Microbiol. 15, 137–148. doi: 10.1038/nrmicro.2016.183
Hartwig, F. P., Davey Smith, G., and Bowden, J. (2017). Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int. J. Epidemiol. 46, 1985–1998. doi: 10.1093/ije/dyx102
Hemani, G., Tilling, K., and Davey, S. G. (2017). Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. 13:e1007081. doi: 10.1371/journal.pgen.1007081
Hofer, U. (2014). Microbiome: pro-inflammatory Prevotella? Nat. Rev. Microbiol. 12:5. doi: 10.1038/nrmicro3180
IEU OpenGWAS project. (2021). Trait: Chronic hepatitis B - IEU OpenGWAS project. Available at: https://gwas.mrcieu.ac.uk/datasets/bbj-a-99/ (Accessed March 5, 2023).
Jang, J.-H., Yeom, M.-J., Ahn, S., Oh, J.-Y., Ji, S., Kim, T.-H., et al. (2020). Acupuncture inhibits neuroinflammation and gut microbial dysbiosis in a mouse model of Parkinson’s disease. Brain Behav. Immun. 89, 641–655. doi: 10.1016/j.bbi.2020.08.015
Joo, E.-J., Cheong, H. S., Kwon, M.-J., Sohn, W., Kim, H.-N., and Cho, Y. K. (2021). Relationship between gut microbiome diversity and hepatitis B viral load in patients with chronic hepatitis B. Gut Pathog. 13:65. doi: 10.1186/s13099-021-00461-1
Khoruts, A., and Sadowsky, M. J. (2016). Understanding the mechanisms of faecal microbiota transplantation. Nat. Rev. Gastroenterol. Hepatol. 13, 508–516. doi: 10.1038/nrgastro.2016.98
Komiyama, S., Yamada, T., Takemura, N., Kokudo, N., Hase, K., and Kawamura, Y. I. (2021). Profiling of tumour-associated microbiota in human hepatocellular carcinoma. Sci. Rep. 11:10589. doi: 10.1038/s41598-021-89963-1
Kurilshikov, A., Medina-Gomez, C., Bacigalupe, R., Radjabzadeh, D., Wang, J., Demirkan, A., et al. (2021). Large-scale association analyses identify host factors influencing human gut microbiome composition. Nat. Genet. 53, 156–165. doi: 10.1038/s41588-020-00763-1
Li, P., Wang, H., Guo, L., Gou, X., Chen, G., Lin, D., et al. (2022). Association between gut microbiota and preeclampsia-eclampsia: a two-sample Mendelian randomization study. BMC Med. 20:443. doi: 10.1186/s12916-022-02657-x
Lin, C.-Y., Hsu, C.-Y., He, H.-R., Chiang, W.-Y., Lin, S.-H., Huang, Y.-L., et al. (2022). Gut microbiota differences between psoriatic arthritis and other undifferentiated arthritis: a pilot study. Medicine (Baltimore) 101:e29870. doi: 10.1097/MD.0000000000029870
Lok, A. S. F., McMahon, B. J., Brown, R. S., Wong, J. B., Ahmed, A. T., Farah, W., et al. (2016). Antiviral therapy for chronic hepatitis B viral infection in adults: a systematic review and meta-analysis. Hepatology 63, 284–306. doi: 10.1002/hep.28280
Long, Y., Tang, L., Zhou, Y., Zhao, S., and Zhu, H. (2023). Causal relationship between gut microbiota and cancers: a two-sample Mendelian randomisation study. BMC Med. 21:66. doi: 10.1186/s12916-023-02761-6
Lu, H., Wu, Z., Xu, W., Yang, J., Chen, Y., and Li, L. (2011). Intestinal microbiota was assessed in cirrhotic patients with hepatitis B virus infection. Intestinal microbiota of HBV cirrhotic patients. Microb. Ecol. 61, 693–703. doi: 10.1007/s00248-010-9801-8
MiBioGen consortium. (2021). MiBioGen. Available at: https://mibiogen.gcc.rug.nl/ (Accessed March 5, 2023).
Milosevic, I., Vujovic, A., Barac, A., Djelic, M., Korac, M., Radovanovic Spurnic, A., et al. (2019). Gut-liver Axis, gut microbiota, and its modulation in the Management of Liver Diseases: a review of the literature. Int. J. Mol. Sci. 20:395. doi: 10.3390/ijms20020395
Moreno-Arrones, O. M., Serrano-Villar, S., Perez-Brocal, V., Saceda-Corralo, D., Morales-Raya, C., Rodrigues-Barata, R., et al. (2020). Analysis of the gut microbiota in alopecia areata: identification of bacterial biomarkers. J. Eur. Acad. Dermatol. Venereol. 34, 400–405. doi: 10.1111/jdv.15885
Morrison, D. J., and Preston, T. (2016). Formation of short chain fatty acids by the gut microbiota and their impact on human metabolism. Gut Microbes 7, 189–200. doi: 10.1080/19490976.2015.1134082
Nikolova, V. L., Smith, M. R. B., Hall, L. J., Cleare, A. J., Stone, J. M., and Young, A. H. (2021). Perturbations in gut microbiota composition in psychiatric disorders: a review and Meta-analysis. JAMA Psychiat. 78, 1343–1354. doi: 10.1001/jamapsychiatry.2021.2573
Plaza-Diaz, J., Ruiz-Ojeda, F. J., Gil-Campos, M., and Gil, A. (2019). Mechanisms of action of probiotics. Adv. Nutr. 10, S49–S66. doi: 10.1093/advances/nmy063
Radjabzadeh, D., Bosch, J. A., Uitterlinden, A. G., Zwinderman, A. H., Ikram, M. A., van Meurs, J. B. J., et al. (2022). Gut microbiome-wide association study of depressive symptoms. Nat. Commun. 13:7128. doi: 10.1038/s41467-022-34502-3
Rasheed, N., Saeed, N., Dyary, H., Faqe, S., Sulaiman, R., Muhammad, P., et al. (2022). Outbreaks of inclusion body hepatitis caused by fowl adenovirus in commercial broiler farms in the Kurdistan region, North Iraq from 2013 to 2021. Pak. Vet. J. 42, 201–207. doi: 10.29261/pakvetj/2021.080
Ren, Z., Li, A., Jiang, J., Zhou, L., Yu, Z., Lu, H., et al. (2019). Gut microbiome analysis as a tool towards targeted non-invasive biomarkers for early hepatocellular carcinoma. Gut 68, 1014–1023. doi: 10.1136/gutjnl-2017-315084
Rinninella, E., Raoul, P., Cintoni, M., Franceschi, F., Miggiano, G. A. D., Gasbarrini, A., et al. (2019). What is the healthy gut microbiota composition? A changing ecosystem across age, environment, diet, and diseases. Microorganisms 7:14. doi: 10.3390/microorganisms7010014
Sakaue, S., Kanai, M., Tanigawa, Y., Karjalainen, J., Kurki, M., Koshiba, S., et al. (2021). A cross-population atlas of genetic associations for 220 human phenotypes. Nat. Genet. 53, 1415–1424. doi: 10.1038/s41588-021-00931-x
Shen, Y., Wu, S.-D., Chen, Y., Li, X.-Y., Zhu, Q., Nakayama, K., et al. (2023). Alterations in gut microbiome and metabolomics in chronic hepatitis B infection-associated liver disease and their impact on peripheral immune response. Gut Microbes 15:2155018. doi: 10.1080/19490976.2022.2155018
Shen, F., Zheng, R.-D., Sun, X.-Q., Ding, W.-J., Wang, X.-Y., and Fan, J.-G. (2017). Gut microbiota dysbiosis in patients with non-alcoholic fatty liver disease. Hepatobiliary Pancreat. Dis. Int. 16, 375–381. doi: 10.1016/S1499-3872(17)60019-5
Strimmer, K. (2008). Fdrtool: a versatile R package for estimating local and tail area-based false discovery rates. Bioinformatics 24, 1461–1462. doi: 10.1093/bioinformatics/btn209
Sun, X., Pan, C. Q., and Xing, H. (2021). Effect of microbiota metabolites on the progression of chronic hepatitis B virus infection. Hepatol. Int. 15, 1053–1067. doi: 10.1007/s12072-021-10230-6
Tan, M., Bhadoria, A. S., Cui, F., Tan, A., Van Holten, J., Easterbrook, P., et al. (2021). Estimating the proportion of people with chronic hepatitis B virus infection eligible for hepatitis B antiviral treatment worldwide: a systematic review and meta-analysis. Lancet Gastroenterol. Hepatol. 6, 106–119. doi: 10.1016/S2468-1253(20)30307-1
Trépo, C., Chan, H. L. Y., and Lok, A. (2014). Hepatitis B virus infection. Lancet 384, 2053–2063. doi: 10.1016/S0140-6736(14)60220-8
Verbanck, M., Chen, C.-Y., Neale, B., and Do, R. (2018). Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat. Genet. 50, 693–698. doi: 10.1038/s41588-018-0099-7
Wang, X., Chen, L., Wang, H., Cai, W., and Xie, Q. (2020). Modulation of bile acid profile by gut microbiota in chronic hepatitis B. J. Cell. Mol. Med. 24, 2573–2581. doi: 10.1111/jcmm.14951
Wang, J., Wang, Y., Zhang, X., Liu, J., Zhang, Q., Zhao, Y., et al. (2017). Gut microbial Dysbiosis is associated with altered hepatic functions and serum metabolites in chronic hepatitis B patients. Front. Microbiol. 8:2222. doi: 10.3389/fmicb.2017.02222
Wu, J., Xu, S., Xiang, C., Cao, Q., Li, Q., Huang, J., et al. (2018). Tongue coating microbiota community and risk effect on gastric Cancer. J. Cancer 9, 4039–4048. doi: 10.7150/jca.25280
Ye, S., Wang, L., Li, S., Ding, Q., Wang, Y., Wan, X., et al. (2022). The correlation between dysfunctional intestinal flora and pathology feature of patients with pulmonary tuberculosis. Front. Cell. Infect. Microbiol. 12:1090889. doi: 10.3389/fcimb.2022.1090889
Yiu, J. H. C., Dorweiler, B., and Woo, C. W. (2017). Interaction between gut microbiota and toll-like receptor: from immunity to metabolism. J. Mol. Med. (Berl) 95, 13–20. doi: 10.1007/s00109-016-1474-4
Zeng, Y., Chen, S., Fu, Y., Wu, W., Chen, T., Chen, J., et al. (2020). Gut microbiota dysbiosis in patients with hepatitis B virus-induced chronic liver disease covering chronic hepatitis, liver cirrhosis and hepatocellular carcinoma. J. Viral Hepat. 27, 143–155. doi: 10.1111/jvh.13216
Keywords: gut microbiota, chronic hepatitis B, Mendelian randomization, genome-wide association study, linkage disequilibrium score regression
Citation: Zhang Q, Zhou J, Zhang X, Mao R and Zhang C (2023) Mendelian randomization supports causality between gut microbiota and chronic hepatitis B. Front. Microbiol. 14:1243811. doi: 10.3389/fmicb.2023.1243811
Edited by:
Khalid Mehmood, Islamia University of Bahawalpur, PakistanReviewed by:
Yong-Xin Liu, Chinese Academy of Agricultural Sciences, ChinaReza Rastmanesh, American Physical Society, United States
Copyright © 2023 Zhang, Zhou, Zhang, Mao and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chuan Zhang, zhangchuan0905@sina.com
†These authors have contributed equally to this work and share first authorship