ORIGINAL RESEARCH article

Front. Cell. Infect. Microbiol., 05 June 2023

Sec. Intestinal Microbiome

Volume 13 - 2023 | https://doi.org/10.3389/fcimb.2023.1119875

Alterations of the gut microbiota associated with the occurrence and progression of viral hepatitis

  • 1. School of Public Health, Guangxi Medical University, Nanning, China

  • 2. Guangxi Key Laboratory of AIDS Prevention and Treatment, Guangxi Universities Key Laboratory of Prevention and Control of Highly Prevalent Disease, Guangxi Medical University, Nanning, China

  • 3. Life Science Institute, Guangxi Medical University, Nanning, China

  • 4. Collaborative Innovation Centre of Regenerative Medicine and Medical BioResource Development and Application Co-constructed by the Province and Ministry, Guangxi Medical University, Nanning, China

  • 5. Guangxi Colleges and Universities Key Laboratory of Prevention and Control of Highly Prevalent Diseases, Guangxi Medical University, Nanning, China

Abstract

Background:

Gut microbiota is the largest population of microorganisms and is closely related to health. Many studies have explored changes in gut microbiota in viral hepatitis. However, the correlation between gut microbiota and the occurrence and progression of viral hepatitis has not been fully clarified.

Methods:

PubMed and BioProject databases were searched for studies about viral hepatitis disease and 16S rRNA gene sequencing of gut microbiota up to January 2023. With bioinformatics analyses, we explored changes in microbial diversity of viral hepatitis, screened out crucial bacteria and microbial functions related to viral hepatitis, and identified the potential microbial markers for predicting risks for the occurrence and progression of viral hepatitis based on ROC analysis.

Results:

Of the 1389 records identified, 13 studies met the inclusion criteria, with 950 individuals including 656 patient samples (HBV, n = 546; HCV, n = 86; HEV, n = 24) and 294 healthy controls. Gut microbial diversity is significantly decreased as the infection and progression of viral hepatitis. Alpha diversity and microbiota including Butyricimonas, Escherichia-Shigella, Lactobacillus, and Veillonella were identified as the potential microbial markers for predicting the risk of development of viral hepatitis (AUC>0.7). Microbial functions including tryptophan metabolism, fatty acid biosynthesis, lipopolysaccharide biosynthesis, and lipid metabolism related to the microbial community increased significantly as the development of viral hepatitis.

Conclusions:

This study demonstrated comprehensively the gut microbiota characteristics in viral hepatitis, screened out crucial microbial functions related to viral hepatitis, and identified the potential microbial markers for predicting the risk of viral hepatitis.

1 Introduction

Viral hepatitis is a significant healthcare burden worldwide, resulting in around 300 million infections and more than 1.3 million deaths each year (Policy, 2017). It is caused by five types of viruses: hepatitis A virus (HAV), hepatitis B virus (HBV), hepatitis C virus (HCV), hepatitis D virus (HDV), and hepatitis E virus (HEV). HBV and HCV are responsible for about 90% of the mortality from viral hepatitis (Jefferies et al., 2018), and long-term infection could lead to other dangerous consequences such as liver cirrhosis (LC), hepatocellular carcinoma (HCC), and even death (Winer et al., 2017; Datfar et al., 2021; Su et al., 2022b). HDV is an incomplete virus that requires HBV for its replication and globally affects around 5% of people with HBV infection(WHO). HAV and HEV are the common infectious etiologies of acute hepatitis throughout the world (Lin et al., 2017; Aslan and Balaban, 2020). Given the great threat of viral hepatitis to public health, it is meaningful to explore how to reverse hepatitis virus infection and to prevent the development of end-stage liver disease.

Researchers have proved that the immune system is involved in one way or another in patients with viral hepatitis (Barathan et al., 2018; Su et al., 2023). It could control the viral infection; however, the long-lasting inflammation could lead to hepatic damage (Hou et al., 2018). Thus, the immune responses are critical for determining the outcome of hepatitis virus infection (Su et al., 2022a; Chi et al., 2023). The gut microbiota is the largest population of microorganisms in the human body. It resides in the intestine and is made up of trillions of bacteria, fungi, and other microbes. It is critical for the host defense and immune homeostasis in humans (Rinninella et al., 2019; Xue et al., 2020). Alterations in the gut microbiota are associated with immune homeostasis disturbances (Yang and Cong, 2021). Recent studies have raised the possibility that gut microbial dysbiosis is involved in the occurrence and progression of viral hepatitis (Chen et al., 2022; Ali et al., 2023).

Hepatitis virus infection alters the diversity of the gut microbiome. One study showed increased diversity of microbiota in the treatment-naive HCV group (Sultan et al., 2021), while another study showed the opposite result, with lower alpha diversity among HCV-infected patients (Aly et al., 2016). Other researchers have found significant changes in the gut microbiome among HBV-induced disease groups when compared with healthy individuals. For example, the abundance of Veilonella, Streptococcus, and Enterococcus increased significantly in the HBV-related acute-on-chronic liver failure (HBV-ACLF) group, and Bacteroides, Lachnospiracea incertae sedis, and Clostridium cluster XIVa were enriched in patients with HBV-related HCC patients (HBV-HCC) (Yang et al., 2018; Huang et al., 2020; Yao et al., 2021). Nevertheless, it is not completely clear whether there is a correlation between the microbiome change in different types of viral hepatitis. How the gut microbiome is altered during the process from hepatitis virus infection to chronic hepatic disorders needs to be studied in greater detail. Besides, how alterations in the gut microbiota affect the development of viral hepatitis is not well understood.

Therefore, we conducted this study using available data of 16S ribosomal RNA (rRNA) gene amplicon sequencing of gut microbiota-related viral hepatitis disease studies. We reanalyzed these data by using rigorous bioinformatics methods to elucidate shifts in the microbial community structure and diversity along with the development of viral hepatitis and to reveal crucial microbial taxonomy and functions related to viral hepatitis, aiming to provide recommendations for delaying disease progression during hepatic exacerbation.

2 Materials and methods

2.1 Search strategy and selection criteria

The PubMed (https://pubmed.ncbi.nlm.nih.gov/) and BioProject (https://www.ncbi.nlm.nih.gov/bioproject) databases were searched for articles and data about viral hepatitis and the gut microbiota published up to January 2023. The search strategy was developed based on keywords, medical subject headings (MeSH) terms, and synonyms (Table S1). Studies were included according to the following criteria: (1) it included a viral hepatitis or its related hepatic disease group and a control group that can be distinguished, (2) it used stool or rectal swab samples used for sequencing, and (3) it employed 16S ribosomal RNA (16S rRNA) gene sequencing and the data available. Raw 16S rRNA gene sequence and metadata information were downloaded from publicly available databases or obtained from the authors. Animal experiments or in vitro studies, reviews, meta-analyses, comments, letters, poster abstracts, and studies with less than three individuals (either case or control groups) were excluded.

2.2 Processing of raw data

The Sequence Read Archive (SRA) files were downloaded from the National Center for Biotechnology Information (NCBI) SRA database and converted to FASTQ format by the SRA Toolkit (https://github.com/ncbi/sratoolkit, Version 3.0.0). These files were then imported into Quantitative Insights in Microbial Ecology (QIIME) version 2 for microbiome bioinformatics analysis(Estaki et al., 2020). The QIIME2 DADA2 plug-in, which takes input dereplicated amplicon sequencing reads as an input, removed the primers and low-quality reads from the sequences, and finally returned the inferred composition of the samples. This information was used to implement sequence quality control for single-end and paired-end 16S rRNA gene reads(Callahan et al., 2016; Estaki et al., 2020). The outputs of the DADA2 pipeline were representative sequences, some statistics on the procedure, and a feature table of amplicon sequence variants (ASVs) with ≥ 97% identity matches to the expected sequences in the extreme dataset(Nearing et al., 2018).

2.3 Filtering data

After generating an ASV table and representative sequences, a q2-feature-table plug-in in QIIME2 was used to filter the low-quality features to improve the quality of the downstream statistical analysis. The samples containing the minimum count of sequences that might be caused by sequencing errors or a low level of the microbiome leading to low uptake of deoxyribonucleic acid (DNA) were filtered. Meanwhile, the low abundance features were filtered based on the interquartile range (IQR). Then the generated filtered sequences were classified with the q2-feature-classifier plug-in based on the SILVA database (https://www.arb-silva.de/, version 138) released in the QIIME2 package(Bokulich et al., 2018).

2.4 Data analysis

Files of the filtered feature table, phylogenetic tree, and taxonomic classifications were imported into RStudio 4.2.1. The vegan package was utilized to normalize the feature table to scale based on each sample’s library size that transformed the feature table into a relative feature table, aiming to remove technical bias caused by variations in sample collection, library preparation, or sequencing manifesting as uneven sampling depth and sparsity, which could not reflect the true difference in the underlying biology(Weiss et al., 2017). Then, the microeco package was used to calculate the alpha diversity indices including richness (Observed Species, Chao1, and ACE), diversity (Shannon, Simpson, Invsimpson, and Fisher), and phylogenetic diversity (PD) to evaluate the overall structure of the gut microbiota. Beta diversity was analyzed with the vegan package, based on principal coordinates (PCoA) and the non-metric multidimensional scaling (NMDS) analyses. Objects that are ordinated closer together have smaller dissimilarity values than those ordinated further farther apart. Next, the permutational multivariate analysis of variance (PERMANOVA) test and analysis of similarities (ANOSIM) based on the Bray-Curtis distance were performed to evaluate the similarities between the groups. Besides, based on the distribution of ASVs, the UpSet plot was used to show the ASV intersections among different groups through the UpSet package.

The SpiecEasi package was employed to identify the different genus and species through the linear discriminant analysis effect size (LEfSe) method (Linear Discriminant Analysis [LDA] score ≥2). The Kyoto Encyclopedia of Genes and Genomes (KEGG) functional pathways related to the microbial community were predicted with the Tax4Fun2 package based on 16S rRNA gene sequencing data classified with the SILVA database. Finally, the area under the curve (AUC) of the receiver operating characteristic (ROC) curve was calculated to evaluate the prediction effectiveness of the alpha diversity indices and differential microbial taxonomy.

The forest plots of the comparisons of the alpha diversity between the case and control groups were generated in Review Manager 5.3. A fixed-effects or random-effects model was selected according to the deviance information criterion (DIC). The standardized mean difference (SMD) and the corresponding 95% credible interval (CI) were used to evaluate the results. The beta diversity and LEfSe analyses were visualized by R 4.1.1. Graphs of the differential KEGG functional pathways were generated in GraphPad Prism 8.0.2. ROC curve analysis was conducted in GraphPad Prism 8.0.2. A model with an AUC score of > 0.7 would be considered acceptable. All tests were two-sided with a P-value of 0.05 set as the threshold for significance.

3 Results

3.1 Study characteristics

Of the 1,317 studies found in PubMed and the 73 records found in BioProject, 13 records related to HBV, HCV, and HEV were eventually included for subsequent analysis. Note that all the records concerning HAV and HDV were excluded (Figure 1). Among the 12 included studies from PubMed, 1 was about HEV infection(Wu et al., 2020), 4 were about HCV infection(Aly et al., 2016; Taylor et al., 2020; Sultan et al., 2021; Ali et al., 2023), and 7 were HBV-related liver disorders (Wang et al., 2017; Liu et al., 2018; Liu et al., 2019; Ni et al., 2019; Chen et al., 2020; Zheng et al., 2020; Li et al., 2022). The remaining unpublished data was from BioProject and concerned HBV (BioProject accession number PRJEB32568). These studies were from China, the United States, and Egypt, and contained a total of 950 individual samples (HBV, n = 546; HCV, n = 86; HEV, n = 24; healthy control [HC], n = 294). More comprehensive details of the included studies are presented in Table 1.

Figure 1

Table 1

Author, Year [Ref]PMIDBioProject accession numberCountryStudy settingsStudy periodSample sizeSample type16S rRNA Variable
Region, Sequencing
Platform, Sequencing technology
Wu et al., 202032500937ERP119119ChinaThe First Affiliated Hospital (College of Medicine) of Zhejiang University2018.05-2019.0512 HEV, 12 HEV-ALFStool samplesV3-V5, Ion Torrent S5 XL, Single read sequencing
Taylor et al., 202033051377ERP122366the United StatesThe University of California San DiegoNA111 HC, 13 HCVStool samplesV4, Illumina MiSeq, Single read sequencing
Liu et al., 201930675188PRJNA428932ChinaNanjing Medical
University Affiliated Cancer Hospital
2016.09-2017.0533 HC, 35 HBV-HCCStool samplesV4, Illumina HiSeq 2500, Single read sequencing
NANAPRJEB32568ChinaNANA5 HC, 12 CHB, 11 HBV-LC, 9 HBV-HCCStool samplesNA, Illumina HiSeq 2500, Single read sequencing
Aly et al., 201627625705PRJNA328966EgyptFaculty of Pharmacy, Cairo University2015.02-2016.098 HC, 6 HCVStool samplesV4, Illumina MiSeq, Paired-end sequencing
Wang et al., 201729180991PRJNA382861ChinaThe affiliated hospital of Shanghai University of Traditional Chinese Medicine, The Infectious Disease Hospital of Ningbo, and the Sixth of People’s Hospital of Shaoxing ZhejiangNA25 HC, 206 HBVStool samplesV3-V4, Illumina MiSeq, Paired-end sequencing
Chen et al., 202032265857PRJNA558158ChinaZhongshan Hospital, affiliated with Xiamen
University
2017.12-2018.0521 HC, 23 HBV, 28 CHB, 25 HBV-LCStool samplesV3-V4, Illumina HiSeq
2500, Paired-end sequencing
Sultan et al., 202133119247PRJNA634402EgyptDepartment of Endemic Hepatology and Gastroenterology, Mansoura University HospitalsNA38 HC, 38 HCVStool samplesV3-V4, Illumina MiSeq, Paired-end sequencing
Liu et al., 201829780327PRJNA445763Chinathe First Affiliated Hospital of Harbin Medical University2015.12-2016.1220 HC, 28 HBV-LCStool samplesV3-V4, Illumina HiSeq
2500, Paired-end sequencing
Ni et al., 201931293562PRJNA478823ChinaZhujiang Hospital, Southern Medical University2017.08-2017.1118 HC, 61 HBV-HCCStool samplesV4-V5, Illumina MiSeq, Paired-end sequencing
Zheng et al., 202032281295PRJNA540574Chinathe First Hospital of Jilin University2017.03-2018.0420 HC, 8 CHB, 35 HBV-HCCStool samplesV4, Illumina HiSeq
2500, Paired-end sequencing
Ali et al., 202336522461PRJNA727609the United Statesthe National Institutes of Health Clinical Center2015.06-2017.0213 HCV, 16 HCV-LCStool samplesV4, Illumina MiSeq, Paired-end sequencing
Li et al., 202235733959PRJNA838083ChinaGuangzhou Panyu Central Hospital2020.10-2021.715 HC, 23 CHB, 20 HBV-LC, 22 HBV-HCCStool samplesV4, Illumina Nova6000, Paired-end sequencing

Study characteristics of the included studies.

NA, Not Applicable.

3.2 Hepatitis virus infection and progression significantly reduce gut microbial diversity

As estimated by the Observed Species, Chao1, Shannon, Simpson, Fisher, and PD indexes, the alpha diversity was reduced significantly in HBV-infected patients (P < 0.05; Figure 2). As HBV infection progressed, there was a significant trend for downregulation in alpha diversity according to the Observed Species, Chao1, ACE, Shannon, Simpson, Invsimpson, and Fisher indexes (P < 0.05; Figure 3). The pooled estimate showed significant decreases in chronic hepatitis B (CHB) (SMD = -0.28; 95%CI, -0.45 to -0.11; P < 0.05; Figure 4) and HCV (SMD = -0.25; 95%CI, -0.39 to -0.10; P < 0.05; Figure 5), compared with the HC group. There was a numerical but nonsignificant trend for downregulation in patients with HBV-related LC (HBV-LC) and HBV-HCC compared with the HC group (Figures S1, S2). The comparisons of the alpha diversity between groups that could not be analyzed comprehensively were shown in the form of boxplots. Compared with HCV-infected patients, the alpha diversity showed significant decreases in the (HCV-LC) group (P < 0.05; Figure 6). Meanwhile, the results of boxplots showing the changes in alpha diversity between HEV-infected individuals and patients with HEV-related acute liver failure (HEV-ALF) displayed a nonsignificant upward trend among the patients with HEV-ALF (Figure S3).

Figure 2

Figure 3

Figure 4

Figure 5

Figure 6

The PCoA and NMDS analyses were conducted with PERMANOVA and ANOSIM, respectively, to evaluate the similarities in the gut microbiome composition among the groups in all the included studies. The PERMANOVA and ANOSIM results for the datasets of PRJEB32568 and PRJNA838083 (Li et al., 2022) reporting four separate groups (HC, CHB, HBV-LC, and HBV-HCC) indicated that as HBV infection progressed, there were significant differences in species composition between the groups (P < 0.01; Figures 7, 8). The PERMANOVA results for the PRJNA558158 (Chen et al., 2020) and PRJNA540574 (Zheng et al., 2020) datasets showed the same trend (P < 0.01; Figure 7). Based on the ASV level, there were differences between the groups— according to the UpSet plots showing the ASV intersections between groups—including the total abundance of ASVs, the number of shared ASVs, and the number of unique ASVs for each patient group (Figure S4).

Figure 7

Figure 8

3.3 Crucial microbiota and microbial functions associated with viral hepatitis

LEfSe analysis was used to identify the dominant microbiota in each group. Based on the LDA section, the abundance of 19 genera, including Alloprevotella, Butyricimonas, and Colidextribacter, was significantly up-regulated in HBV-infected patients, while Bacteroides, Parabacteroides, and Sutterella were down-regulated, compared with those in the HC group (P < 0.05; Figure 9). Among the HCV-infected patients, 10 taxa including Desulfovibrio, Eubacterium eligens, and Prevotalla were increased significantly, while 11 genera including Barnesiella, Colidextribacter, and Dorea were decreased significantly, compared with those in the HC group (P < 0.05; Figure 10). More details of LEfSe analysis were shown in Tables S2-S4.

Figure 9

Figure 10

Based on 16S rRNA gene sequences, the KEGG profile was constructed to predict the microbial community function. The dominant microbial functions related to viral hepatitis were summarized in Figure 11. The results showed that 88 microbial functions including tryptophan metabolism, fatty acid biosynthesis, and lipopolysaccharide (LPS) biosynthesis were remarkably increased in HBV-related liver disorders (Figure 11A), while 14 microbial functions including lipid metabolism and thiamine metabolism were significantly enriched with HCV infection and progression (Figure 11B).

Figure 11

3.4 Alpha diversity and Butyricimonas, Veillonella, Escherichia-Shigella, and Lactobacillus may serve as potential gut microbial markers to predict the risk for viral hepatitis

ROC curve analysis was conducted to evaluate the potential to use the gut microbiota as a non-invasive marker to predict the risk for viral hepatitis. Those with an AUC score > 0.7 are considered to have a high risk. In the model of HBV/HCV infection and progression, all alpha diversity metrics reached AUC scores of > 0.7 (P < 0.05; Figure 12). The Observed Species, Fisher, and InvSimpson indexes reached AUC values of 0.824, 0.841, and 0.943, respectively, in the model predicting the risk for the occurrence of HBV-LC, HBV-HCC, and HCV (Figures 12B, C, E), and the Observed Species and Shannon indexes reached AUC values of 0.752 and 0.843, respectively, for predicting the progression of HBV and HCV (Figures 12D, F).

Figure 12

Among the crucial microbiota screened out based on the LDA section, 10 genera including Butyricimonas, Veillonella, Escherichia-Shigella, and Lactobacillus had a high potential to predict the risk of HBV progression (AUC > 0.7; Figure 13C). Among them, the AUC score of the particular model of Butyricimonas was 0.775 to predict the risk of HBV infection (Figure 13A). The model of Butyricimonas, Veillonella, and Escherichia-Shigella had AUC scores of 0.917, 0.797, and 0.794, respectively, to predict HBV-LC (Figure 13B). Besides, Veillonella and Lactobacillus reached AUC scores of 0.857 and 0.745, respectively, between HCV-infected patients and HC (Figure 13D), while genera including Clostridia_UCG-014, Dorea, Monoglobus, and Ruminococcus decreased in the HCV-infected group, with AUC values of 0.900, 0.857, 0.886, and 0.857, respectively, compared with the HC group (Figure 14).

Figure 13

Figure 14

4 Discussion

In recent years, an increasing number of studies have explored the correlation between gut microbiota perturbations and the occurrence and progression of viral hepatitis. For example, a reduction in gut microbial diversity has been linked to the severity of the disease and poor health (Wilmanski et al., 2019). Changes in the abundance of some gut bacteria can lead to the secretion of anti-inflammatory factors, thus accelerating disease progression (Peng et al., 2022). Moreover, gut microbiome-derived metabolites such as LPS and bile acids could interact with liver immune cells leading to pathological effects (Milosevic et al., 2021). However, the existing research could not identify consistent microbial taxa that respond to the disease (Wang et al., 2017; Chen et al., 2020). Few studies have systematically assessed the association between the gut microbiome and viral hepatitis, and the implication of the gut microbiota on the progression of viral hepatitis remains unclear. Therefore, we comprehensively assessed gut microbiota perturbations across a spectrum of the occurrence and progression of viral hepatitis to evaluate the reproducibility and specificity of potential gut microbial biomarkers.

Alpha diversity summarizes the structure of an ecological community regarding its richness, evenness, or both (Willis, 2019). It is a validated marker of gastrointestinal health and metabolic disorders (Plassais et al., 2021). We identified showed significant decreases in gut microbial alpha diversity after HBV/HCV infection as well as lower diversity was detected as the disease progressed, as has been observed in previous studies (Aly et al., 2016; Chen et al., 2020). This result followes the general assumption that higher microbial diversity is more beneficial to host health (Nikolova et al., 2021). Meanwhile, the AUC values of the diversity indexes were high, indicating that the change in the gut microbial alpha diversity could be a potential indicator to predict the risk for HBV/HCV infection and progression. An nonsignificant trend for upregulation in the patients with HEV-ALF versus the HEV-infected individuals demonstrated that there might be compensatory modulation in the early stage of disease progression to reestablish gut homeostasis. However, it is unclear what occurs during HAV and HDV infection due to the absence of datasets. Regarding beta diversity, we noted consistent significant dissimilarities among communities at different HBV progression stages in two datasets according to both PCoA and NMDS analyses, indicating that the gut microbiome composition was significantly alteredas HBV infection progressed. However, there were insignificant differences in beta diversity between HBV/HCV-infected patients and HC. Thus the suitability of beta diversity as a biomarker needs to be further confirmed.

LEfSe analysis showed that compared with the HC, 19 genera were dominant in HBV-infected individuals and 10 genera were dominant in HCV-infected individuals, suggesting that the changes in the abundance of crucial taxa led to significant alterations in the gut microbiome composition in viral hepatitis. Among these dominant genera, Prevotella could lead to a reduction in short-chain acids (SCFAs) (Iljazovic et al., 2021), which play a crucial role in delaying the progression of HBV-related HCC (Mcbrearty et al., 2021). Lactobacillus, Escherichia-Shigella, and Veillonella could contribute to the production of pro-inflammatory factors such as LPS and tumor necrosis factor (TNF) (Liu et al., 2022, Mata Forsberg et al., 2019). These processes may lead to the development of viral hepatitis by affecting the immune responses in vivo (Ji et al., 2019; Radzikowska et al., 2019; Joo et al., 2021). However, we were unable to identify the dominant microbiota related to HAV, HDV, and HEV infection due to the lack of a dataset or insufficient information provided for comprehensive analysis. At the same time, among HCV-infected individuals, we identified four decreased genera Clostridia_UCG-014, Dorea, Monoglobus, and Ruminococcus as probiotics with great potential for the prevention and treatment of HCV. Besides, we found that 10 genera, including Butyricimonas, Escherichia-Shigella, Lactobacillus, and Veillonella, could have great potential to distinguish HBV- or HCV-infected individuals from HC and to predict the risk for the development of HBV infection. The upregulation of the crucial bacteria could influence the production of gut metabolites, including tryptophan, LPS, fatty acids, and lipids (Aujoulat et al., 2014; Huang et al., 2019; Mao et al., 2020). Changes in metabolites could be associated with some microbial functions that contribute to the pathogenesis and progression of the disease (Ren et al., 2020).

In this study, functions including tryptophan metabolism, LPS biosynthesis, fatty acid biosynthesis, and lipid metabolism were significantly enriched in the HBV- or HCV-infected group. Tryptophan is an essential amino acid that possesses diverse metabolic, neurological, and immunological roles, and it is involved in viral infections including HBV (Mehraj and Routy, 2015; Fiore and Murray, 2021). A high level of LPS can promote the secretion of inflammatory cytokines including TNF and interleukin-6 (IL-6), stimulating immune cells and finally may lead to the progression of viral hepatitis (Feng et al., 2020). Fatty acids are known to play diverse roles in immune cells. It can activate the inflammatory cell signaling pathways via cell surface or intracellular receptors (Calder, 2011). Abnormal fatty acid levels in the liver can result in synergistic induction of HBV-related proteins and liver inflammatory factors, which might affect HBV progression (Cho et al., 2014). Lipid metabolism is intimately connected to every step of the HCV life cycle, and HCV enhances its replication by modulating lipid metabolism in the host cell (Popescu et al., 2014). All of these findings suggest that changes in specific taxa can alter the production of metabolites, which may contribute to the development of viral hepatitis through diverse microbial functions. However, few studies have been able to clarify the implication of alteration in the gut microbiome and microbial functions in viral hepatitis, which is of great significance in understanding the occurrence and progression of viral hepatitis.

This study has several limitations. First, we failed to collect all of the data from gut microbiome—related viral hepatitis studies due to the non-availability of data or inadequate information about the detailed characteristics such as grouping or disease stages. Additionally, the existing evidence limits our analysis: For example, no studies have evaluated the relationship between the gut microbiota and HAV or HDV. Finally, the authors of the included studies applied inconsistent variable regions and instruments for 16S rRNA gene sequencing, which may have generated potential bias during analysis.

5 Conclusions

In conclusion, our study demonstrates that the occurrence and progression of viral hepatitis are accompanied by a significant decrease in gut microbial diversity. The decreased gut microbial alpha diversity as well as the increased abundance of genera including Butyricimonas, Escherichia-Shigella, Veillonella, and Lactobacillus have the greatest potential to serve as biomarkers to predict the risk for viral hepatitis. Meanwhile, Clostridia_UCG-014, Dorea, Monoglobus, and Ruminococcus have potential value as probiotics for the prevention and treatment of viral hepatitis. Crucial microbial functions including tryptophan metabolism, fatty acid biosynthesis, LPS biosynthesis, and lipid metabolism—related to the significantly upregulated microbial community as viral hepatitis progresses—play diverse roles in the activation of immune cells and inflammatory responses. These processes provide a valuable direction to confirm the association between the gut microbiota and the occurrence and progression of viral hepatitis.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author contributions

JH, HL, and PC designed the study. XinY and HM performed the study mainly including writing the manuscript, interpreting the result, and preparing the report for publication. LY supervised the process. JZ, ZL, QW, LL, FL, XipY, and BG participated in data acquisition. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by funds from the National Natural Science Foundation of China (NSFC, Grant No. 82060366, 82273694, and 82160385) and the Guangxi Natural Science Foundation (Grant No. 2018GXNSFAA050099).

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/fcimb.2023.1119875/full#supplementary-material

Supplementary Figure 1

The forest plots of alpha diversity comparing HBV-LC patients to HC individuals. Results showed a numerical but no significant decrease among HBV-LC patients (n=74) in each diversity index, compared with HC group (n=41).

Supplementary Figure 2

The forest plots of alpha diversity comparing HBV-HCC patients to HC individuals. Results showed a numerical but no significant decrease among HBV-HCC patients (n=126) in each diversity index, compared with HC group (n=65).

Supplementary Figure 3

Boxplots of comparisons of gut microbial alpha diversity between HEV-infected patients and HEV-ALF groups. Results based on the dataset of ERP119119 indicated that the diversity metrics displayed a nonsignificant upward trend among the HEV-ALF patients (n=12) compared with those in the HEV-infected patients (n=12).

Supplementary Figure 4

The UpSet plots show the ASV intersections among groups in the included studies. For example, in the dataset of PRJEB32568, the total number of ASVs was 939, and 360 ASVs were shared in both groups, while 26, 105, and 43 ASVs were unique for HBV-HCC, HBV-LC, and CHB groups, respectively.

References

  • 1

    AliR. O.QuinnG. M.UmarovaR.HaddadJ. A.ZhangG. Y.TownsendE. C.et al. (2023). Longitudinal multi-omics analyses of the gut-liver axis reveals metabolic dysregulation in hepatitis c infection and cirrhosis. Nat. Microbiol.8, 1227. doi: 10.1038/s41564-022-01273-y

  • 2

    AlyA. M.AdelA.El-GendyA. O.EssamT. M.AzizR. K. (2016). Gut microbiome alterations in patients with stage 4 hepatitis c. Gut Pathog.8, 42. doi: 10.1186/s13099-016-0124-2

  • 3

    AslanA. T.BalabanH. Y. (2020). Hepatitis e virus: epidemiology, diagnosis, clinical manifestations, and treatment. World J. Gastroenterol.26, 55435560. doi: 10.3748/wjg.v26.i37.5543

  • 4

    AujoulatF.BouvetP.Jumas-BilakE.Jean-PierreH.MarchandinH. (2014). Veillonella seminalis sp. nov., a novel anaerobic gram-stain-negative coccus from human clinical samples, and emended description of the genus veillonella. Int. J. Syst. Evol. Microbiol.64, 35263531. doi: 10.1099/ijs.0.064451-0

  • 5

    BarathanM.MohamedR.YongY. K.KannanM.VadiveluJ.SaeidiA.et al. (2018). Viral persistence and chronicity in hepatitis c virus infection: role of T-cell apoptosis, senescence and exhaustion. Cells7, 165. doi: 10.3390/cells7100165

  • 6

    BokulichN. A.KaehlerB. D.RideoutJ. R.DillonM.BolyenE.KnightR.et al. (2018). Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin. Microbiome6, 90. doi: 10.1186/s40168-018-0470-z

  • 7

    CalderP. C. (2011). Fatty acids and inflammation: the cutting edge between food and pharma. Eur. J. Pharmacol.668 Suppl 1, S50S58. doi: 10.1016/j.ejphar.2011.05.085

  • 8

    CallahanB. J.McmurdieP. J.RosenM. J.HanA. W.JohnsonA. J.HolmesS. P. (2016). DADA2: high-resolution sample inference from illumina amplicon data. Nat. Methods13, 581583. doi: 10.1038/nmeth.3869

  • 9

    ChenB.HuangH.PanC. Q. (2022). The role of gut microbiota in hepatitis b disease progression and treatment. J. Viral Hepat29, 94106. doi: 10.1111/jvh.13595

  • 10

    ChenZ.XieY.ZhouF.ZhangB.WuJ.YangL.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

  • 11

    ChiH.ZhaoS.YangJ.GaoX.PengG.ZhangJ.et al. (2023). T-Cell exhaustion signatures characterize the immune landscape and predict HCC prognosis via integrating single-cell RNA-seq and bulk RNA-sequencing. Front. Immunol.14, 1137025. doi: 10.3389/fimmu.2023.1137025

  • 12

    ChoH. K.KimS. Y.YooS. K.ChoiY. H.CheongJ. (2014). Fatty acids increase hepatitis b virus X protein stabilization and HBx-induced inflammatory gene expression. FEBS J.281, 22282239. doi: 10.1111/febs.12776

  • 13

    DatfarT.DoulberisM.PapaefthymiouA.HinesI. N.ManziniG. (2021). Viral hepatitis and hepatocellular carcinoma: state of the art. Pathogens10, 1366. doi: 10.3390/pathogens10111366

  • 14

    EstakiM.JiangL.BokulichN. A.McdonaldD.GonzalezA.KosciolekT.et al. (2020). QIIME 2 enables comprehensive end-to-End analysis of diverse microbiome data and comparative studies with publicly available data. Curr. Protoc. Bioinf.70, e100. doi: 10.1002/cpbi.100

  • 15

    FengS.AchouteL.MargolskeeR. F.JiangP.WangH. (2020). Lipopolysaccharide-induced inflammatory cytokine expression in taste organoids. Chem. Senses45, 187194. doi: 10.1093/chemse/bjaa002

  • 16

    FioreA.MurrayP. J. (2021). Tryptophan and indole metabolism in immune regulation. Curr. Opin. Immunol.70, 714. doi: 10.1016/j.coi.2020.12.001

  • 17

    ForsbergM. M.BjorkanderS.PangY.LundqvistL.NdiM.OttM.et al. (2019). Extracellular membrane vesicles from lactobacilli dampen IFN-gamma responses in a monocyte-dependent manner. Sci. Rep.9, 17109. doi: 10.1038/s41598-019-53576-6

  • 18

    HouX. J.YeF.LiX. Y.LiuW. T.JingY. Y.HanZ. P.et al. (2018). Immune response involved in liver damage and the activation of hepatic progenitor cells during liver tumorigenesis. Cell Immunol.326, 5259. doi: 10.1016/j.cellimm.2017.08.004

  • 19

    HuangN.HuaD.ZhanG.LiS.ZhuB.JiangR.et al. (2019). Role of actinobacteria and coriobacteriia in the antidepressant effects of ketamine in an inflammation model of depression. Pharmacol. Biochem. Behav.176, 93100. doi: 10.1016/j.pbb.2018.12.001

  • 20

    HuangH.RenZ.GaoX.HuX.ZhouY.JiangJ.et al. (2020). Integrated analysis of microbiome and host transcriptome reveals correlations between gut microbiota and clinical outcomes in HBV-related hepatocellular carcinoma. Genome Med.12, 102. doi: 10.1186/s13073-020-00796-5

  • 21

    IljazovicA.RoyU.GalvezE. J. C.LeskerT. R.ZhaoB.GronowA.et al. (2021). Perturbation of the gut microbiome by prevotella spp. enhances host susceptibility to mucosal inflammation. Mucosal Immunol.14, 113124. doi: 10.1038/s41385-020-0296-4

  • 22

    JefferiesM.RauffB.RashidH.LamT.RafiqS. (2018). Update on global epidemiology of viral hepatitis and preventive strategies. World J. Clin. cases6, 589599. doi: 10.12998/wjcc.v6.i13.589

  • 23

    JiY.YinY.LiZ.ZhangW. (2019). Gut microbiota-derived components and metabolites in the progression of non-alcoholic fatty liver disease (NAFLD). Nutrients11, 1712. doi: 10.3390/nu11081712

  • 24

    JooE. J.CheongH. S.KwonM. J.SohnW.KimH. N.ChoY. 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

  • 25

    LiR.YiX.YangJ.ZhuZ.WangY.LiuX.et al. (2022). Gut microbiome signatures in the progression of hepatitis b virus-induced liver disease. Front. Microbiol.13, 916061. doi: 10.3389/fmicb.2022.916061

  • 26

    LinK. Y.ChenG. J.LeeY. L.HuangY. C.ChengA.SunH. Y.et al. (2017). Hepatitis a virus infection and hepatitis a vaccination in human immunodeficiency virus-positive patients: a review. World J. Gastroenterol.23, 35893606. doi: 10.3748/wjg.v23.i20.3589

  • 27

    LiuY.JinY.LiJ.ZhaoL.LiZ.XuJ.et al. (2018). Small bowel transit and altered gut microbiota in patients with liver cirrhosis. Front. Physiol.9, 470. doi: 10.3389/fphys.2018.00470

  • 28

    LiuQ.LiF.ZhuangY.XuJ.WangJ.MaoX.et al. (2019). Alteration in gut microbiota associated with hepatitis b and non-hepatitis virus related hepatocellular carcinoma. Gut Pathog.11, 1. doi: 10.1186/s13099-018-0281-6

  • 29

    LiuL.LiuY.GuoX.JinX.YanW.LinB.et al. (2022). Activation of p38 mitogen-activated protein kinase pathway by lipopolysaccharide aggravates postoperative ileus in colorectal cancer patients. J. Gastroenterol. Hepatol.37, 518530. doi: 10.1111/jgh.15760

  • 30

    MaoJ.QiS.CuiY.DouX.LuoX. M.LiuJ.et al. (2020). Lactobacillus rhamnosus GG attenuates lipopolysaccharide-induced inflammation and barrier dysfunction by regulating MAPK/NF-kappaB signaling and modulating metabolome in the piglet intestine. J. Nutr.150, 13131323. doi: 10.1093/jn/nxaa009

  • 31

    McbreartyN.ArzumanyanA.BichenkovE.MeraliS.MeraliC.FeitelsonM. (2021). Short chain fatty acids delay the development of hepatocellular carcinoma in HBx transgenic mice. Neoplasia23, 529538. doi: 10.1016/j.neo.2021.04.004

  • 32

    MehrajV.RoutyJ. P. (2015). Tryptophan catabolism in chronic viral infections: handling uninvited guests. Int. J. Tryptophan Res.8, 4148. doi: 10.4137/IJTR.S26862

  • 33

    MilosevicI.RussoE.VujovicA.BaracA.StevanovicO.GittoS.et al. (2021). Microbiota and viral hepatitis: state of the art of a complex matter. World J. Gastroenterol.27, 54885501. doi: 10.3748/wjg.v27.i33.5488

  • 34

    NearingJ. T.DouglasG. M.ComeauA. M.LangilleM. G. I. (2018). Denoising the denoisers: an independent evaluation of microbiome sequence error-correction approaches. PeerJ6, e5364. doi: 10.7717/peerj.5364

  • 35

    NiJ.HuangR.ZhouH.XuX.LiY.CaoP.et al. (2019). Analysis of the relationship between the degree of dysbiosis in gut microbiota and prognosis at different stages of primary hepatocellular carcinoma. Front. Microbiol.101458. doi: 10.3389/fmicb.2019.01458

  • 36

    NikolovaV. L.SmithM. R. B.HallL. J.CleareA. J.StoneJ. M.YoungA. H. (2021). Perturbations in gut microbiota composition in psychiatric disorders: a review and meta-analysis. JAMA Psychiatry78, 13431354. doi: 10.1001/jamapsychiatry.2021.2573

  • 37

    PengY. C.XuJ. X.ZengC. F.ZhaoX. H.LiL. Q.QiL. N. (2022). Gut microbiome dysbiosis in patients with hepatitis b virus-related hepatocellular carcinoma after extended hepatectomy liver failure. Ann. Transl. Med.10, 549. doi: 10.21037/atm-22-1958

  • 38

    PlassaisJ.Gbikpi-BenissanG.FigarolM.ScheperjansF.GorochovG.DerkinderenP.et al. (2021). Gut microbiome alpha-diversity is not a marker of parkinson’s disease and multiple sclerosis. Brain Commun.3, fcab113. doi: 10.1093/braincomms/fcab113

  • 39

    Policy, C. F. I. D. R. A (2017) WHO: viral hepatitis deaths increasing. Available at: https://www.cidrap.umn.edu/news-perspective/2017/04/who-viral-hepatitis-deaths-increasing.

  • 40

    PopescuC. I.RivaL.VlaicuO.FarhatR.RouilleY.DubuissonJ. (2014). Hepatitis c virus life cycle and lipid metabolism. Biol. (Basel)3, 892921. doi: 10.3390/biology3040892

  • 41

    RadzikowskaU.RinaldiA. O.Celebi SozenerZ.KaraguzelD.WojcikM.CyprykK.et al. (2019). The influence of dietary fatty acids on immune responses. Nutrients11, 2990. doi: 10.3390/nu11122990

  • 42

    RenZ.FanY.LiA.ShenQ.WuJ.RenL.et al. (2020). Alterations of the human gut microbiome in chronic kidney disease. Adv. Sci. (Weinh)7, 2001936. doi: 10.1002/advs.202001936

  • 43

    RinninellaE.RaoulP.CintoniM.FranceschiF.MiggianoG. A. D.GasbarriniA.et al. (2019). What is the healthy gut microbiota composition? a changing ecosystem across age, environment, diet, and diseases. Microorganisms7, 14. doi: 10.3390/microorganisms7010014

  • 44

    SuK.GuoL.MaW.WangJ.XieY.RaoM.et al. (2022a). PD-1 inhibitors plus anti-angiogenic therapy with or without intensity-modulated radiotherapy for advanced hepatocellular carcinoma: a propensity score matching study. Front. Immunol.13, 972503. doi: 10.3389/fimmu.2022.972503

  • 45

    SuK.LiuY.WangP.HeK.WangF.ChiH.et al. (2022b). Heat-shock protein 90alpha is a potential prognostic and predictive biomarker in hepatocellular carcinoma: a large-scale and multicenter study. Hepatol. Int.16, 12081219. doi: 10.1007/s12072-022-10391-y

  • 46

    SuK.ShenQ.TongJ.GuT.XuK.LiH.et al. (2023). Construction and validation of a nomogram for HBV-related hepatocellular carcinoma: a large, multicenter study. Ann. Hepatol.28, 101109. doi: 10.1016/j.aohep.2023.101109

  • 47

    SultanS.El-MowafyM.ElgamlA.El-MeseryM.El ShabrawiA.ElegezyM.et al. (2021). Alterations of the treatment-naive gut microbiome in newly diagnosed hepatitis c virus infection. ACS Infect. Dis.7, 10591068. doi: 10.1021/acsinfecdis.0c00432

  • 48

    TaylorB. C.WeldonK. C.EllisR. J.FranklinD.McdonaldD.HumphreyG.et al. (2020). Reduced independence in daily living is associated with the gut microbiome in people with HIV and HCV. mSystems5, e00528-20. doi: 10.1128/mSystems.00528-20

  • 49

    WangJ.WangY.ZhangX.LiuJ.ZhangQ.ZhaoY.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

  • 50

    WeissS.XuZ. Z.PeddadaS.AmirA.BittingerK.GonzalezA.et al. (2017). Normalization and microbial differential abundance strategies depend upon data characteristics. Microbiome5, 27. doi: 10.1186/s40168-017-0237-y

  • 51

  • 52

    WillisA. D. (2019). Rarefaction, alpha diversity, and statistics. Front. Microbiol.10, 2407. doi: 10.3389/fmicb.2019.02407

  • 53

    WilmanskiT.RappaportN.EarlsJ. C.MagisA. T.ManorO.LovejoyJ.et al. (2019). Blood metabolome predicts gut microbiome alpha-diversity in humans. Nat. Biotechnol.37, 12171228. doi: 10.1038/s41587-019-0233-9

  • 54

    WinerB. Y.HuangT. S.PludwinskiE.HellerB.WojcikF.LipkowitzG. E.et al. (2017). Long-term hepatitis b infection in a scalable hepatic co-culture system. Nat. Commun.8, 125. doi: 10.1038/s41467-017-00200-8

  • 55

    WuJ.HuangF.LingZ.LiuS.LiuJ.FanJ.et al. (2020). Altered faecal microbiota on the expression of Th cells responses in the exacerbation of patients with hepatitis e infection. J. Viral Hepat27, 12431252. doi: 10.1111/jvh.13344

  • 56

    XueJ.AjuwonK. M.FangR. (2020). Mechanistic insight into the gut microbiome and its interaction with host immunity and inflammation. Anim. Nutr.6, 421428. doi: 10.1016/j.aninu.2020.05.007

  • 57

    YangW.CongY. (2021). Gut microbiota-derived metabolites in the regulation of host immune responses and immune-related inflammatory diseases. Cell Mol. Immunol.18, 866877. doi: 10.1038/s41423-021-00661-4

  • 58

    YangR.XuY.DaiZ.LinX.WangH. (2018). The immunologic role of gut microbiota in patients with chronic HBV infection. J. Immunol. Res.2018, 2361963. doi: 10.1155/2018/2361963

  • 59

    YaoX.YuH.FanG.XiangH.LongL.XuH.et al. (2021). Impact of the gut microbiome on the progression of hepatitis b virus related acute-on-Chronic liver failure. Front. Cell Infect. Microbiol.11, 573923. doi: 10.3389/fcimb.2021.573923

  • 60

    ZhengR.WangG.PangZ.RanN.GuY.GuanX.et al. (2020). Liver cirrhosis contributes to the disorder of gut microbiota in patients with hepatocellular carcinoma. Cancer Med.9, 42324250. doi: 10.1002/cam4.3045

Summary

Keywords

gut microbiota, 16s ribosomal RNA gene amplicon sequencing, viral hepatitis, liver disease, microbial markers

Citation

Yang X, Mai H, Zhou J, Li Z, Wang Q, Lan L, Lu F, Yang X, Guo B, Ye L, Cui P, Liang H and Huang J (2023) Alterations of the gut microbiota associated with the occurrence and progression of viral hepatitis. Front. Cell. Infect. Microbiol. 13:1119875. doi: 10.3389/fcimb.2023.1119875

Received

09 December 2022

Accepted

22 May 2023

Published

05 June 2023

Volume

13 - 2023

Edited by

Mahmoud Salami, Kashan University of Medical Sciences, Iran

Reviewed by

Sayyed Alireza Talaei, Kashan University of Medical Sciences, Iran; Gang Tian, Southwest Medical University, China; Jinyan Yang, Southwest Medical University, China

Updates

Copyright

*Correspondence: Jiegang Huang, ; Hao Liang, ; Ping Cui,

†These authors have contributed equally to this work and share first authorship

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics