- Research
- Open access
- Published:
Diversity and composition of gut microbiota in healthy individuals and patients at different stages of hepatitis B virus-related liver disease
Gut Pathogens volume 15, Article number: 24 (2023)
Abstract
Background
Hepatitis B virus (HBV) causes chronic hepatitis B (CHB), liver cirrhosis, and hepatocellular carcinoma. The evolution of human gut microbiota during the progression of HBV-related liver diseases remains unclear. Therefore, we prospectively enrolled patients with HBV-related liver diseases and healthy individuals. Through 16S ribosomal RNA amplicon sequencing, we characterized the gut microbiota of the participants and predicted the functions of microbial communities.
Results
We analyzed the gut microbiota of 56 healthy controls and 106 patients with HBV-related liver disease [14 with resolved HBV infection, 58 with CHB, and 34 with advanced liver disease (15 with liver cirrhosis and 19 with hepatocellular carcinoma)]. Patients with HBV-related liver disease exhibited a higher degree of bacterial richness (all P < 0.05) than did healthy controls. Beta diversity analyses revealed a distinct clustering pattern between healthy controls and patients with HBV-related liver disease (all P < 0.05). The composition of bacteria (from the phylum level to the genus level) varied across the stages of liver disease. Linear discriminant analysis effect size revealed multiple taxa that differ significantly in abundance between healthy controls and patients with HBV-related liver disease; however, fewer differences were observed among patients with resolved HBV infection, those with CHB, and those with advanced liver disease. The ratio of Firmicutes to Bacteroidetes was increased in all three patient groups compared with the ratio in healthy controls (all P < 0.001). The analysis of the sequencing data by using PICRUSt2 revealed the changes in microbial functions with disease progression.
Conclusions
The diversity and composition of gut microbiota appear to vary significantly between healthy controls and patients at different stages of HBV-related liver disease. The understanding of gut microbiota may provide novel therapeutic options in these patients.
Background
An estimated 296 million people are chronically infected with hepatitis B virus (HBV) worldwide, with 1.5 million new infections each year [1]. Thus, HBV infection remains a major threat to public health. Patients infected with HBV can develop acute or chronic hepatitis, liver cirrhosis (LC), and hepatocellular carcinoma (HCC) [1]. The term gut microbiota is referred to the microorganisms that colonize in the intestinal tract. The number of genes in the gut microbiota genome is 150-fold higher than in the human genome [2]. Accumulating evidence suggests a strong association between the liver and the intestines [3, 4]. The liver secretes primary bile acids and antimicrobial molecules into the intestinal lumen to facilitate digestion and control intestinal bacterial overgrowth. The portal system carries intestinal microbial products such as microbe-associated molecular patterns and endotoxins to the liver, facilitating the generation of a proinflammatory state [5].
Multiple studies have explored the alterations in the composition of microbiota in various liver diseases, including LC, alcoholic liver disease, and nonalcoholic fatty liver disease [6]. Gut microbiota may play a role in the pathogenesis of various HBV-related liver diseases [6]. Alterations have been demonstrated in the composition of gut microbiota in patients with HBV-related liver disease; these alterations are associated with a reduction in the abundance of beneficial bacteria and an increase in the growth rate of pathogenic species [7,8,9,10,11]. Nonetheless, the composition of intestinal microbiota may be affected by various host and environmental factors such as physical and mental health, medication, diet, and exposome [12]; therefore, the literature remains inconclusive.
A functional cure for chronic hepatitis B (CHB) is achieved through the seroclearance of hepatitis B surface antigen (HBsAg) [1]. The composition of microbiota in patients with resolved HBV infection (resolved HBV) remains unknown. Therefore, we compared the gut microbiota characteristics at different stages of HBV-related liver diseases (resolved HBV, CHB, and advanced stages of liver diseases) to identify potential taxonomic biomarkers and microbial functional profiles and thus improve our understanding of the gut microbiota in these patients.
Results
Clinical characteristics of the study population
In this study, we analyzed the gut microbiota of 56 healthy controls and 106 patients with HBV-related liver diseases. Among these patients, 14 had resolved HBV, 58 had CHB, and 34 had advanced liver disease (15 had LC, whereas 19 had HCC). Of the patients with HCC, 12 (63%) had a history of LC. The baseline characteristics of the participants are summarized in Table 1. The median age was 28 years for healthy controls, 57 years for resolved HBV patients, 51 years for CHB patients, and 62 years for advanced liver disease patients (P < 0.0001); the corresponding proportions of men were 32%, 50%, 66%, and 82% (P < 0.0001), respectively. In general, the levels of alanine aminotransferase (ALT), aspartate aminotransferase (AST), blood urea nitrogen, creatinine, and fasting blood glucose were higher in patients with resolved HBV, those with CHB, and those with advanced liver disease than in healthy controls (all P < 0.0001). Patients with advanced liver disease had lower platelet counts than did those with resolved HBV or CHB (P = 0.0019).
Overview of the sequencing data
The 16S ribosomal RNA (rRNA) gene V3–V4 amplicon sequencing yielded 15,801,571 reads, and the average number of sequence reads per sample was 97,540 (minimum: 46,986; maximum: 183,523). To account for sample size–related variations, rarefaction was performed to normalize the reads to 46,986 reads per sample. A total of 7,504 amplicon sequence variants (ASVs) were retained after rarefaction.
Differences in gut microbiome structures
The common and unique ASVs among the different stages of liver diseases and the number of observed ASVs are presented in Fig. 1A, B. Alpha diversity matrices were calculated on the basis of rarefied ASVs. The estimated values of Simpson and Shannon diversities revealed higher degrees of bacterial richness and evenness in patients with resolved HBV than in healthy controls (Fig. 1C, D). Chao1 and Chao2 indices indicated a higher degree of species richness in patients with HBV-related liver disease than in healthy controls (all P < 0.05; Fig. 1E, F).
Principal coordinate analysis (PCoA) performed using weighted and unweighted UniFrac distances revealed distinct clusters between healthy controls and the aforementioned three patient groups (Fig. 2A, B). The adonis2 function, which was used to evaluate between-group dissimilarities and adjust for the possible confounding effect of age in the multivariate model, revealed a nonsignificant effect of age (P = 0.671 for weighted UniFrac distance and P = 0.146 for unweighted UniFrac distance) on microbiota composition compared with the significant effect of disease status (P = 0.001 for both distances). Pairwise comparisons performed using the permutational multivariate analysis of variance (PERMANOVA) indicated significant differences in bacterial composition between patients with resolved HBV and healthy controls, between patients with CHB and healthy controls, and between patients with advanced liver disease and healthy controls (all q < 0.05 for weighted and unweighted UniFrac distances; Additional file 1). The composition of fecal microbiota communities was relatively similar among the three patient groups.
Bacterial abundance and taxonomic distribution
The most abundant phyla in the study cohort were Bacteroidetes (44.7%), Firmicutes (41.4%), Proteobacteria (6.0%), Actinobacteriota (5.5%), and Verrucomicrobiota (1.0%). The most abundant genera were Bacteroides (30.3%), Prevotella (6.5%), Bifidobacterium (3.4%), Faecalibacterium (3.3%), Blautia (3.1%), Parabacteroides (3.1%), Escherichia-Shigella (2.6%), and Megamonas (2.0%). The taxonomic distribution of the predominant bacteria among the study groups is shown in Fig. 3A–E. The ratio of Firmicutes to Bacteroidetes (F/B ratio) was higher in all three patient groups than in healthy controls (all P < 0.001; Fig. 3F).
Differences in gut microbiota compositions
Linear discriminant analysis (LDA) effect size (LEfSe) was used to identify specific taxonomic biomarkers for the study groups. The cutoff value of the logarithmic linear discriminant analysis (LDA) score was set at 3.0 to determine the major taxonomic differences between the following groups: (1) healthy controls and patients with resolved HBV, (2) healthy controls and patients with CHB, (3) healthy controls and patients with advanced liver disease, (4) patients with CHB and those with resolved HBV, and (5) patients with CHB and those with advanced liver disease. Significant differences were observed between healthy controls and patients with HBV-related liver disease (Fig. 4A–E).
Regarding phylum-level comparisons, Firmicutes, Verrucomicrobiota, and Fusobacteria were enriched and Bacteroidetes was depleted in the three patient groups compared with the findings in healthy controls. Furthermore, Proteobacteria and Actinobacteriota were enriched in patients with advanced liver disease. The abundance of Proteobacteria was higher in patients with advanced liver disease than in those with CHB (Additional file 2).
Genus-level comparisons between patients with resolved HBV and healthy controls (Fig. 4A) revealed that the most enriched genera in patients with resolved HBV were Veillonella (LDA score = 4.08; P < 0.001), Blautia (LDA score = 4.05; P = 0.003), Akkermansia (LDA score = 4.01; P = 0.004), Streptococcus (LDA score = 3.74; P < 0.001), and Clostridia UCG-014 (LDA score = 3.73; P = 0.017), whereas the most enriched genera in healthy controls were Bacteroides (LDA score = 4.94; P = 0.004), Parabacteroides (LDA score = 4.04; P = 0.008), and Sutterella (LDA score = 3.84; P < 0.001).
Genus-level comparisons between patients with CHB and healthy controls (Fig. 4B) revealed that the most enriched genera in patients with CHB were Streptococcus (LDA score = 4.15; P < 0.001), Blautia (LDA score = 4.01; P = 0.001), Veillonella (LDA score = 3.87; P < 0.001), Fusobacterium (LDA score = 3.87; P = 0.007), and Akkermansia (LDA score = 3.66; P = 0.012), whereas the most enriched genera in healthy controls were Bacteroides (LDA score = 4.89; P < 0.001), Megamonas (LDA score = 4.17; P = 0.005), Parabacteroides (LDA score = 3.95; P < 0.001), Sutterella (LDA score = 3.65; P = 0.001), and Lachnoclostridium (LDA score = 3.44; P = 0.005).
Genus-level comparisons between patients with advanced liver disease and healthy controls (Fig. 4C) revealed that the most enriched genera in patients with advanced liver disease were Escherichia-Shigella (LDA score = 4.26; P = 0.018), Veillonella (LDA score = 4.13; P < 0.001), Streptococcus (LDA score = 4.10; P < 0.001), Blautia (LDA score = 4.04; P < 0.001), and Collinsella (LDA score = 3.96; P = 0.002), whereas the most enriched genera in healthy controls were Bacteroides (LDA score = 4.87; P = 0.001), Prevotella (LDA score = 4.56; P < 0.001), Megamonas (LDA score = 4.09; P = 0.003), Parabacteroides (LDA score = 3.98; P < 0.001), and Faecalibacterium (LDA score = 3.92; P = 0.002).
Genus-level comparisons between patients with resolved HBV and those with CHB (Fig. 4D) revealed that Clostridium sensu stricto 1 (LDA score = 3.39; P = 0.012) and Romboutsia (LDA score = 3.20; P = 0.041) were enriched in patients with resolved HBV, whereas no taxon was enriched in patients with CHB.
Finally, genus-level comparisons between patients with CHB and those with advanced liver disease (Fig. 4E) revealed that Prevotella (LDA score = 4.36; P = 0.013), Faecalibacterium (LDA score = 3.81; P = 0.026), Dialister (LDA score = 3.44; P = 0.003), and Eubacterium ruminantium group (LDA score = 3.15; P = 0.039) were enriched in patients with CHB, whereas Lactobacillus (LDA score = 3.60; P = 0.028) and CAG-352 (LDA score = 3.03; P = 0.027) were enriched in patients with advanced liver disease.
Fig. 5 shows the mean relative abundance of the following 11 selected genera: Bacteroides, Megamonas, Parabacteroides, and Sutterella, which were enriched in healthy controls; Akkermansia, Blautia, Lactobacillus, Streptococcus, and Veillonella, which were enriched in the patient groups (all P for trend < 0.01); and Clostridium sensu stricto 1 and Romboutsia, which had higher abundances in patients with resolved HBV than in those with CHB.
Functional analysis based on 16S rRNA sequencing data
We identified 45 metabolic functions that were classified as Level 2 of the Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthologs (Additional file 3). A metagenome function prediction analysis performed using the Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2 (PICRUSt2) tool revealed alterations in bacterial functions with disease progression. A total of 28 differentially expressed pathways were identified (Additional file 4). A heatmap was generated to visualize the differential expression levels of Level 2 KEGG pathways (Fig. 6). The pathways enriched with disease progression were involved in membrane transport, cancer, bacterial infection, viral infection, transcription, endocrine and metabolic disease, and cell motility. By contrast, pathways enriched in healthy controls were involved in cell growth and death, endocrine system, carbohydrate, glycan, lipid, terpenoids, and polyketides metabolism, transport and catabolism, energy metabolism, circulatory system, nervous system, and immune disease. Notably, pathways involved in the metabolism of cofactors and vitamins and the development of antineoplastic drug resistance exhibited the highest levels of enrichment in patients with resolved HBV.
Subgroup analysis in patients with CHB
Patients with CHB were further divided into the following four groups (representing the four phases of HBV natural history) according to the presence of the hepatitis B e antigen (HBeAg) and the level of ALT: HBeAg (+) chronic HBV infection, HBeAg (+) chronic hepatitis B, HBeAg (−) chronic HBV infection, and HBeAg (−) chronic hepatitis B. The alpha and beta diversities were compared between the four groups and patients with resolved HBV. The values of the alpha diversity indices, such as Shannon diversity, Simpson diversity, Chao1, and Chao2, and the number of observed ASVs were lower in patients with HBeAg (−) chronic hepatitis B (Additional file 5). PCoA plots generated on the basis of weighted and unweighted UniFrac distances revealed no distinct separation between the aforementioned four groups and patients with resolved HBV (all P > 0.05; Additional file 6 and Additional file 7).
Discussion
HBV infection remains difficult to cure despite the development of antiviral therapies [13]. Therefore, targeting gut microbiota has emerged as a novel therapy for HBV infection and its complications [14]. In this study, we identified the alterations in the diversity and composition of gut microbiota between healthy controls and patients with resolved HBV, CHB, or advanced liver disease to elucidate the relationship between the composition of gut microbiota and the progression of HBV-related liver disease.
In this study, healthy controls were younger than patients with HBV-related liver disease. A large cohort-based study revealed no association between microbial alpha diversity and age in a cohort from China, in contrast to the findings in cohorts from the United States, the United Kingdom, and Colombia [15]. Our findings also indicated minimal confounding effect of age on the microbiota composition. Thus, the differences in microbiota composition between healthy controls and the three patient groups might have been due to the different disease statuses, rather than the baseline characteristics of the patients.
In our study, the species richness estimated using Chao1 and Chao2 indices and the observed ASVs were increased in the three patient groups compared with those in healthy controls. Shannon and Simpson diversities were similar between the healthy controls, patients with CHB, and patients with advanced liver disease but higher in patients with resolved HBV than in healthy controls. By contrast, a study reported lower Chao1 index values and smaller operational taxonomic unit (OTU) numbers in patients with LC than in healthy controls and asymptomatic HBV carriers [8]. Furthermore, another study revealed lower Shannon diversity, Simpson diversity, and Chao1 index values in patients with LC than in those with early-stage HBV-related disease [7]. However, other studies have reported no significant between-group differences in alpha diversity [9, 16]. These discrepancies might be attributed to the heterogeneity in patient population.
The gut microbiome structure, assessed in terms of beta diversity, varied qualitatively and quantitatively between healthy controls and patients with HBV-related liver disease. Our finding is consistent with that reported by Zeng et al.[9] Beta diversity analyses performed using Bray − Curtis-based nonmetric multidimensional scaling also revealed alterations in microbiota composition in patients with HBV infection, those with CHB, and those with LC in another study [7].
The phyla consistently enriched in the three patient groups were Firmicutes, Verrucomicrobiota, and Fusobacteria; by contrast, Bacteroidetes was depleted in the patient groups. Proteobacteria and Actinobacteriota were enriched in patients with advanced liver disease. Firmicutes and Bacteroidetes are two predominant phyla found in the gut environment. Bacteroidetes are gram-negative bacteria; the members of this phylum perform metabolic functions such as the degradation of proteins and complex sugar polymers. By contrast, Firmicutes are gram-positive bacteria. This phylum includes probiotic species such as Lactobacillus sp. and opportunistic pathogens such as Clostridium difficile.[17] A marked reduction in the abundance of Bacteroidetes was noted in patients with HBV-related LC [18, 19]. This reduction might have been because of the overgrowth of the pathogenic strains of Proteobacteria and Fusobacteria [19]. An increase in the Firmicutes to Bacteroidetes ratio, which is a marker of dysbiosis, was indicated to be associated with HBV-related LC [20]. In our study, the F/B ratio observed in patients with resolved HBV was more similar to the ratios observed in patients with CHB and those with advanced liver disease than to the ratio observed in healthy controls. This finding suggests the presence of dysbiosis in the patient groups.
The genera Akkermansia, Veillonella, Blautia, and Streptococcus were consistently enriched, whereas the genera Sutterella, Parabacteroides, and Bacteroides were depleted in patients with HBV-related liver disease. Bacteroides participates in carbohydrate and lipid metabolism and contributes to short-chain fatty acid (SCFA) synthesis [21]. In addition, Bacteroides plays a role in the intestinal immune network by activating intestinal dendritic cells, which help control the overgrowth of pathogenic strains and maintain the integrity of the intestinal epithelial barrier [22, 23]. Changes in the microenvironment of the intestinal tract with the progression of liver diseases may reduce the abundance of Bacteroides and facilitate the colonization of opportunistic pathogens. Veillonella has been reported to be enriched in patients with advanced-stage HBV infection [7,8,9,10, 24]. This bacterial genus is associated with inflammation and can thus serve as a pathogenic biomarker of chronic liver inflammation [25]. The abundance of another opportunistic pathogen, Streptococcus, was reported to be higher in patients with LC than in healthy controls [8, 10].
Akkermansia was reported to be most enriched in patients with LC, with relatively low abundance in the healthy controls, patients with CHB, and patients with HCC [9]. However, we found that Akkermansia was generally enriched in patients with HBV-related liver disease. A key species of Akkermansia is A. muciniphila, which is a mucin-degrading bacterium that is associated with a lower risk of obesity. This species also performs several metabolic functions and maintains gastrointestinal mucosal integrity [26, 27]. Future studies should aim to identify the factors contributing to the increased abundance of Akkermansia in patients with HBV-related CHB, LC, or HCC and the role of this species in the progression of the aforementioned disease.
Prevotella, Faecalibacterium, Dialister, and Eubacterium ruminantium group were enriched in patients with CHB, whereas Lactobacillus and CAG-352 were enriched in those with advanced liver disease. Lactobacillus spp. are commonly used probiotics and are effective in ameliorating gastrointestinal problems such as gastrointestinal tract infections and diarrhea, inflammatory bowel disease, nonalcoholic fatty liver disease, and nonalcoholic steatohepatitis [28]. However, consistent with our findings, studies have reported an association between the increased abundance of Lactobacillus and the progression of HBV-related liver diseases, such as decompensated LC and HCC [8, 16, 24]. Further studies on the metabolic functions of Lactobacillus in HBV-related liver disease are needed.
In this study, the abundances of Clostridium sensu stricto 1 and Romboutsia were higher in patients with resolved HBV than in those with CHB. Lu et al. demonstrated an increase of the abundance of Clostridium sensu stricto 1 in patients with CHB receiving entecavir for 8Â weeks [29]. Furthermore, an animal study revealed that the abundance of Clostridium sensu stricto was negatively correlated with the levels of HBsAg and HBeAg after entecavir treatment [30]. The family Clostridiaceae includes bacteria that produce SCFAs, such as acetate, butyrate, and propionate, through carbohydrate fermentation [31]. SCFAs participate in multiple metabolic functions and function as anti-inflammatory agents [32]. The genus Romboutsia has recently been separated from the genus Clostridium and has been described as a nonpathogenic commensal found in the gastrointestinal tract [33]. Thus, the increased abundances of Clostridium and Romboutsia species may indicate reduced liver inflammation and HBsAg seroclearance.
Notably, although patients with resolved HBV achieved the seroclearance of HBsAg and normalization of hepatic inflammation either spontaneously or through antiviral therapies, the alterations in their gut microbiota persisted and were more similar to the microbial profiles noted in patients with chronic liver diseases than to the profiles noted in healthy controls. This might be because the long-term chronic infection had already shaped the composition of patients’ gut microbiota and the alterations in gut microbiota composition after HBsAg seroclearance might not have been sufficiently significant to reach the preset LDA cutoff score of 3.0.
PICRUSt2 analysis revealed that the metabolic function of viral infectious disease was increasingly activated with disease severity in a dose-dependent manner. Unlike two relevant studies [8, 9], our study indicated that the lipid metabolism pathway was enriched in healthy controls, possibly because of the higher abundance of Bacteroidetes in these individuals than in the patient groups. Furthermore, we identified 28 differentially expressed metabolic pathways; most of these pathways varied between healthy controls and the three patient groups. Therefore, the interactions between microbiota composition and metagenomic functional profiles might have contributed to the progression of HBV-related disease.
Our study has some limitations. First, because of the cross-sectional design of this study, we could obtain limited data regarding the role of microbiota in the progression of HBV-related liver diseases. Second, the smaller sample size of patients with resolved HBV allowed us to draw only preliminary conclusions. Finally, the results obtained through the random sampling of fecal samples and the subsequent rarefaction of reads might not have represented the actual microbiota composition in the gastrointestinal tract. To offer in-depth species-level information and clarify metabolic functions, shotgun metagenomic sequencing should be performed in future studies aimed at elucidating the mechanisms underlying the pathogenesis of HBV-related liver diseases.
Conclusions
The diversity and composition of gut microbiota appear to vary significantly between healthy individuals and patients at different stages of HBV-related liver disease. Gut dysbiosis persists even after the seroclearance of HBsAg (resolved HBV). In the future, studies should be conducted on the shift and evolution of microbiota in CHB, LC, HCC, and resolved HBV patients to elucidate the mechanisms underlying the pathogenesis of HBV-related liver disease and develop effective microbiome-based treatment strategies for these diseases.
Methods
Patients
Patients with HBV-related liver disease (resolved HBV, CHB, or advanced liver disease [LC or HCC]) were prospectively enrolled at National Taiwan University Hospital between 2019 and 2021. Healthy individuals with no known history of liver or gastrointestinal diseases, such as inflammatory bowel disorders, irritable bowel syndrome, or colitis, were enrolled in 2016.
Chronic HBV infection was confirmed on the basis of the presence of HBsAg for at least 6 months. Resolved HBV was defined as seronegativity for HBsAg but seropositivity for hepatitis B core antibody (anti-HBc) and normal ALT levels in patients with a known history of CHB. LC was diagnosed on the basis of histological findings, a fibroscan result of ≥ 12 kPa, or the observation of a nodular liver surface, a coarse liver parenchymal texture, and narrowed vessels with irregular intrahepatic vessel contours on ultrasonography [34]. HCC was diagnosed on the basis of pathological findings or one typical dynamic imaging study according to the American Association for the Study of Liver Diseases guideline [35]. Patients were excluded if (1) they had a history of hepatitis C, hepatitis D, or human immunodeficiency virus coinfection; (2) they had autoimmune disorders or used any immunomodulatory drug; (3) they had a history of antibiotics use in the previous 30 days; or (4) they had hepatobiliary malignancy other than HCC.
This study was conducted in accordance with the ethical principles for medical research involving human subjects of the Declaration of Helsinki (updated version, 2013), and recommendations outlined in the International Conference on Harmonization guideline on Good Clinical Practice Guideline. Our study was approved by the Institutional Review Board of National Taiwan University Hospital (approval numbers: 201809039RIND and 201912040RINA). Written informed consent was obtained from all participants.
Laboratory measurements
Baseline biochemical values including alanine aminotransferase (ALT), aspartate aminotransferase (AST), blood urea nitrogen (BUN), creatinine (CRE), fasting blood glucose (GLU-AC), platelet count, total cholesterol (T-CHO), high-density lipoprotein (HDL), low-density lipoprotein (LDL), and triglyceride (TG) were collected at enrollment. Laboratory assays were performed using commercially standardized automated techniques. The serum levels of HBeAg, anti-HBe, anti-HBs, anti-HCV, and anti-HDV were tested using commercial kits (Abbott Laboratories, Abbott Park, IL, USA).
Fecal samples collection and DNA extraction
All fecal samples were frozen at − 20 °C immediately after collection by the participants at enrollment and were transported to the laboratory under frozen status. Bacterial DNA was extracted using the QIAamp PowerFecal Pro DNA Kit (Qiagen, Valencia, CA, USA), according to the manufacturer’s instructions. After extraction, aliquoted DNA samples were stored at – 80 ℃ until retrieval for this study. DNA concentration and integrity were evaluated using NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Hudson, NH, United States) and through agarose gel electrophoresis, respectively.
16S rRNA gene sequencing
Two-step polymerase chain reaction was performed for library preparation according to the 16S library preparation guide published by Illumina. Primer pair sequences (forward: 5′- TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGCCTACGGGNGGCWGCAG-3′; reverse: 5′- GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGGACTACHVGGGTATCTAATCC-3′) were used to sequence the V3 and V4 regions of the bacterial 16S rRNA gene. Dual indices and Illumina sequencing adapters were attached through polymerase chain reaction using Nextera XT DNA Library Preparation Kit according to the manufacturer’s protocol. Library quantification was performed for quality control before sequencing by Agilent Technologies 2100 Bioanalyzer. Then, the pooled libraries were sequenced on the Illumina MiSeq platform with v3 reagents (paired-end sequencing; 2 × 300 bps) [36].
Bioinformatics and statistical analysis
The DADA2 and QIIME 2 (version 2019.10) pipelines were used for sequence denoising and further processing [37, 38]. The taxonomy was assigned using a naïve Bayes classifier trained on the SILVA 138 99% full-length 16S rRNA gene sequence database. Samples were rarefied to 46,986 reads through random sampling, and alpha and beta diversities were measured using the R (version 4.1.3; the R Project for Statistical Computing) package phyloseq [39].
Alpha diversity parameters, such as observed ASVs, Shannon diversity, Simpson diversity, Chao1 index, and Chao2 index, were assessed to estimate species evenness and richness. The R package ggvenn was used to create the corresponding Venn diagram. The Kruskal–Wallis test was performed to compare the study groups in terms of alpha diversity. Beta diversity matrices (weighted and unweighted UniFrac distances) were generated, and principal coordinate analysis (PCoA) was performed to estimate the between-group dissimilarities in microbiota composition. The permutational multivariate analysis of variance (PERMANOVA) with 999 permutations was performed using the Adonis2 function and pairwiseAdonis package in R was used to analyze the statistical significance of beta diversity. The results of alpha and beta diversity analyses were visualized using the R package ggplot2.
Linear discriminant analysis (LDA) effect size (LEfSe) was performed using the Galaxy browser (http://huttenhower.sph.harvard.edu/galaxy/) to identify the taxa that were most likely to be associated with between-group differences [40]. Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2 (PICRUSt2) was used to predict microbial community functions on the basis of the Kyoto Encyclopedia of Genes and Genomes (KEGG) database [41, 42]. Kruskal–Wallis test was used for the between-group comparisons of the differential expression of Level 2 KEGG pathways.
Continuous data are presented in terms of median (interquartile range) values, and categorical data are presented in terms of number (percentage). Between-group differences were assessed using the Kruskal–Wallis, analysis of variance, or chi-square tests, as appropriate. The Benjamini–Hochberg procedure was used to control the false discovery rate. The Jonckheere–Terpstra test was used to evaluate trends of microbial abundance across different disease statuses. Statistical analysis was performed using STATA (version 16.0; Stata Corporation, College Station, TX, USA) and R. All tests were two-sided and a P value < 0.05 indicated statistical significance.
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the National Center for Biotechnology Information (NCBI) database repository, [https://www.ncbi.nlm.nih.gov/bioproject/PRJNA872871].
Change history
17 June 2023
A Correction to this paper has been published: https://doi.org/10.1186/s13099-023-00555-y
Abbreviations
- ALT:
-
Alanine aminotransferase
- AST:
-
Aspartate aminotransferase
- ASV:
-
Amplicon sequence variant
- BMI:
-
Body mass index
- BUN:
-
Blood urea nitrogen
- CHB:
-
Chronic hepatitis B
- CRE:
-
Creatinine
- F/B ratio:
-
Firmicutes/Bacteroidetes ratio
- GLU-AC:
-
Fasting blood glucose
- HBeAg:
-
Hepatitis B e antigen
- HBsAg:
-
Hepatitis B surface antigen
- HBV:
-
Hepatitis B virus
- HCC:
-
Hepatocellular carcinoma
- HCV:
-
Hepatitis C virus
- HDL:
-
High-density lipoprotein
- HDV:
-
Hepatitis D virus
- LC:
-
Liver cirrhosis
- LDL:
-
Low-density lipoprotein
- OTU:
-
Operational taxonomic unit
- SCFA:
-
Short-chain fatty acid
- T-CHO:
-
Total cholesterol
- TG:
-
Triglyceride
References
Jeng WJ, Papatheodoridis GV, Lok ASF. Hepatitis B. Lancet. 2023;401(10381):1039–52.
Chassaing B, Etienne-Mesmin L, Gewirtz AT. Microbiota-liver axis in hepatic disease. Hepatology. 2014;59(1):328–39.
Albillos A, de Gottardi A, Rescigno M. The gut-liver axis in liver disease: pathophysiological basis for therapy. J Hepatol. 2020;72(3):558–77.
Tilg H, Adolph TE, Trauner M. Gut-liver axis: pathophysiological concepts and clinical implications. Cell Metab. 2022;34(11):1700–18.
Tripathi A, Debelius J, Brenner DA, Karin M, Loomba R, Schnabl B, et al. The gut-liver axis and the intersection with the microbiome. Nat Rev Gastroenterol Hepatol. 2018;15(7):397–411.
Albhaisi SAM, Bajaj JS, Sanyal AJ. Role of gut microbiota in liver disease. Am J Physiol Gastrointest Liver Physiol. 2020;318(1):G84-g98.
Chen Z, Xie Y, Zhou F, Zhang B, Wu J, Yang L, et al. Featured gut microbiomes associated with the progression of chronic hepatitis B disease. Front Microbiol. 2020;11:383.
Yang XA, Lv F, Wang R, Chang Y, Zhao Y, Cui X, et al. Potential role of intestinal microflora in disease progression among patients with different stages of Hepatitis B. Gut Pathog. 2020;12:50.
Zeng Y, Chen S, Fu Y, Wu W, Chen T, Chen J, et al. Gut microbiota dysbiosis in patients with hepatitis B virus-induced chronic liver disease covering chronic hepatitis, liver cirrhosis and hepatocellular carcinoma. J Viral Hepat. 2020;27(2):143–55.
Tang Y, Zhou H, Xiang Y, Cui F. The diagnostic potential of gut microbiome for early hepatitis B virus-related hepatocellular carcinoma. Eur J Gastroenterol Hepatol. 2021;33(1S Suppl 1):e167–75.
Lu H, Wu Z, Xu W, Yang J, Chen Y, Li L. Intestinal microbiota was assessed in cirrhotic patients with hepatitis B virus infection. Intestinal microbiota of HBV cirrhotic patients. Microb Ecol. 2011;61(3):693–703.
Gacesa R, Kurilshikov A, Vich Vila A, Sinha T, Klaassen MAY, Bolte LA, et al. Environmental factors shaping the gut microbiome in a Dutch population. Nature. 2022;604(7907):732–9.
Zeisel MB, Lucifora J, Mason WS, Sureau C, Beck J, Levrero M, et al. Towards an HBV cure: state-of-the-art and unresolved questions–report of the ANRS workshop on HBV cure. Gut. 2015;64(8):1314–26.
Paratore M, Santopaolo F, Cammarota G, Pompili M, Gasbarrini A, Ponziani FR. Fecal microbiota transplantation in patients with HBV infection or other chronic liver diseases: update on current knowledge and future perspectives. J Clin Med. 2021;10:12.
de la Cuesta-Zuluaga J, Kelley ST, Chen Y, Escobar JS, Mueller NT, Ley RE, et al. Age- and sex-dependent patterns of gut microbial diversity in human adults. mSystems. 2019;4:4.
Shu W, Shanjian C, Jinpiao L, Qishui O. Gut microbiota dysbiosis in patients with hepatitis B virus-related cirrhosis. Ann Hepatol. 2022;27(2):100676.
Rajilić-Stojanović M, de Vos WM. The first 1000 cultured species of the human gastrointestinal microbiota. FEMS Microbiol Rev. 2014;38(5):996–1047.
Wei X, Yan X, Zou D, Yang Z, Wang X, Liu W, et al. Abnormal fecal microbiota community and functions in patients with hepatitis B liver cirrhosis as revealed by a metagenomic approach. BMC Gastroenterol. 2013;13:175.
Chen Y, Yang F, Lu H, Wang B, Chen Y, Lei D, et al. Characterization of fecal microbial communities in patients with liver cirrhosis. Hepatology. 2011;54(2):562–72.
Deng YD, Peng XB, Zhao RR, Ma CQ, Li JN, Yao LQ. The intestinal microbial community dissimilarity in hepatitis B virus-related liver cirrhosis patients with and without at alcohol consumption. Gut Pathog. 2019;11:58.
Jandhyala SM, Talukdar R, Subramanyam C, Vuyyuru H, Sasikala M, Nageshwar RD. Role of the normal gut microbiota. World J Gastroenterol. 2015;21(29):8787–803.
Lutgendorff F, Akkermans LM, Söderholm JD. The role of microbiota and probiotics in stress-induced gastro-intestinal damage. Curr Mol Med. 2008;8(4):282–98.
He B, Xu W, Santini PA, Polydorides AD, Chiu A, Estrella J, et al. Intestinal bacteria trigger T cell-independent immunoglobulin A(2) class switching by inducing epithelial-cell secretion of the cytokine APRIL. Immunity. 2007;26(6):812–26.
Liu Q, Li F, Zhuang Y, Xu J, Wang J, Mao X, et al. Alteration in gut microbiota associated with hepatitis B and non-hepatitis virus related hepatocellular carcinoma. Gut Pathog. 2019;11:1.
Kummen M, Holm K, Anmarkrud JA, Nygård S, Vesterhus M, Høivik ML, et al. The gut microbial profile in patients with primary sclerosing cholangitis is distinct from patients with ulcerative colitis without biliary disease and healthy controls. Gut. 2017;66(4):611–9.
Zhang T, Li Q, Cheng L, Buch H, Zhang F. Akkermansia muciniphila is a promising probiotic. Microb Biotechnol. 2019;12(6):1109–25.
Zhou Q, Zhang Y, Wang X, Yang R, Zhu X, Zhang Y, et al. Gut bacteria Akkermansia is associated with reduced risk of obesity: evidence from the American Gut Project. Nutr Metab (Lond). 2020;17:90.
Goldstein EJ, Tyrrell KL, Citron DM. Lactobacillus species: taxonomic complexity and controversial susceptibilities. Clin Infect Dis. 2015;60(Suppl 2):S98-107.
Lu YX, He CZ, Wang YX, Ai ZS, Liang P, Yang CQ. Effect of entecavir on the intestinal microflora in patients with chronic hepatitis B: a controlled cross-sectional and longitudinal real-world study. Infect Dis Ther. 2021;10(1):241–52.
Li X, Wu S, Du Y, Yang L, Li Y, Hong B. Entecavir therapy reverses gut microbiota dysbiosis induced by hepatitis B virus infection in a mouse model. Int J Antimicrob Agents. 2020;56(1):106000.
Louis P, Flint HJ. Diversity, metabolism and microbial ecology of butyrate-producing bacteria from the human large intestine. FEMS Microbiol Lett. 2009;294(1):1–8.
Blaak EE, Canfora EE, Theis S, Frost G, Groen AK, Mithieux G, et al. Short chain fatty acids in human gut and metabolic health. Benef Microbes. 2020;11(5):411–55.
Gerritsen J, Fuentes S, Grievink W, van Niftrik L, Tindall BJ, Timmerman HM, et al. Characterization of romboutsia ilealis gen. nov., sp. nov., isolated from the gastro-intestinal tract of a rat, and proposal for the reclassification of five closely related members of the genus clostridium into the genera romboutsia gen. nov., intestinibacter gen. nov., terrisporobacter gen. nov. and asaccharospora gen. nov. Int J Syst Evol Microbiol. 2014;64(5):1600–16.
Su TH, Liao CH, Liu CH, Huang KW, Tseng TC, Yang HC, et al. Acoustic radiation force impulse US imaging: liver stiffness in patients with chronic hepatitis B with and without antiviral therapy. Radiology. 2018;288(1):293–9.
Marrero JA, Kulik LM, Sirlin CB, Zhu AX, Finn RS, Abecassis MM, et al. Diagnosis, staging, and management of hepatocellular carcinoma: 2018 practice guidance by the American association for the study of liver diseases. Hepatology. 2018;68(2):723–50.
Lin CH, Chen CC, Chiang HL, Liou JM, Chang CM, Lu TP, et al. Altered gut microbiota and inflammatory cytokine responses in patients with Parkinson’s disease. J Neuroinflammation. 2019;16(1):129.
Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJ, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13(7):581–3.
Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):852–7.
McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8(4):e61217.
Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12(6):R60.
Douglas GM, Maffei VJ, Zaneveld JR, Yurgel SN, Brown JR, Taylor CM, et al. PICRUSt2 for prediction of metagenome functions. Nat Biotechnol. 2020;38(6):685–8.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Acknowledgements
We thank Mr. Bin-Chen Hsieh, Miss Pei-Ying Chiang, Yun-Yi Lin, and Pei-Ying Yang for their laboratory and administrative assistance, and the staff of the Department of Medical Research, National Taiwan University Hospital for the Integrated Medical Database (NTUH-iMD). We thank the 7th and the 8th Core Laboratory of the Department of Medical Research, National Taiwan University Hospital for technical assistance.
Funding
This work was supported by grants from the National Science and Technology Council, Taiwan (grant numbers MOST 109-2326-B-002-012-MY3, MOST 110-2326-B-400-004, MOST 110-2628-B-002-041), National Taiwan University Hospital (grant numbers VN-112-13, 112-TMU151, 112-S0245, 111-L1005), and the Liver Disease Prevention & Treatment Research Foundation, Taiwan.
Author information
Authors and Affiliations
Contributions
MJL contributed to study concept and design, analysis and interpretation of data, and drafting the manuscript. THS contributed to study concept and design, acquisition, analysis and interpretation of data, drafting the manuscript, and obtaining funding. CCC contributed to study concept and design, acquisition, analysis and interpretation of data. WKW contributed to acquisition, analysis and interpretation of data. SJH contributed to acquisition, analysis and interpretation of data. TCT contributed to acquisition, analysis and interpretation of data. SHL contributed to acquisition, analysis and interpretation of data. CMH contributed to acquisition, analysis and interpretation of data. HCY contributed to acquisition, analysis and interpretation of data. CJL contributed to acquisition, analysis and interpretation of data. MSW contributed to study concept and design, acquisition, analysis and interpretation of data. JHK contributed to study concept and design, analysis and interpretation of data, obtaining funding, critical revision of the manuscript for important intellectual content, and study supervision. All authors read and approved the final manuscript.Â
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
This study was conducted in accordance with the ethical principles for medical research involving human subjects of the Declaration of Helsinki (updated version, 2013), and recommendations outlined in the International Conference on Harmonization guideline on Good Clinical Practice Guideline. Our study was approved by the Institutional Review Board of National Taiwan University Hospital (approval numbers: 201809039RIND and 201912040RINA). Written informed consent was obtained from all participants.
Consent for publication
Not applicable.
Competing interests
T.-H. S. received research grant from Gilead Sciences, and was on speaker’s bureaus for Abbvie, Bayer, Bristol-Myers Squibb, Gilead Sciences, Merck Sharp and Dohme, and Takeda. J.-H. K. has served as a consultant for Abbvie, Gilead Sciences, Merck Sharp and Dohme, and Roche and on speaker’s bureaus for Abbvie, Bristol-Myers Squibb, Gilead Sciences, Merck Sharp and Dohme. Others declare no conflict of interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1:
Table S1. Statistical significance of beta diversity distance matrices.
Additional file 2:
Table S2. Differentially abundant phyla between groups identified by LEfSe.
Additional file 3:
Table S3. The average abundance of predicted Level 2 KEGG pathways in each group.
Additional file 4:
Table S4. Comparison of Level 2 KEGG pathways.
Additional file 5: Figure S1.
Comparisons of bacterial diversity and richness of patients with HBeAg (+) chronic HBV infection, HBeAg (+) chronic hepatitis B, HBeAg (−) chronic HBV infection, HBeAg (−) chronic hepatitis B, and resolved HBV. A A Venn diagram displays the unique and shared ASVs among the four groups. B The HBeAg (−) chronic hepatitis B group had the least observed ASVs among the five groups. The HBeAg (−) chronic hepatitis B group had the lowest alpha diversity indices, including C Shannon diversity, D Simpson diversity, E Chao1 index, and F Chao2 index. * means P < 0.05 and ** means P < 0.01.
Additional file 6: Figure S2.
Beta diversity indices of patients with HBeAg (+) chronic HBV infection, HBeAg (+) chronic hepatitis B, HBeAg (−) chronic HBV infection, HBeAg (−) chronic hepatitis B, and resolved HBV. A PCoA plot of bacterial beta diversity based on the weighted UniFrac distance. B PCoA plot of bacterial beta diversity based on the unweighted UniFrac distance. No separate cluster was found between HBeAg (+) chronic HBV infection, HBeAg (+) chronic hepatitis B, HBeAg (−) chronic HBV infection, HBeAg (−) chronic hepatitis B, and resolved HBV.
Additional file 7:
Table S5. Statistical significance of beta diversity distance matrices among patients with HBeAg (+) chronic HBV infection, HBeAg (+) chronic hepatitis B, HBeAg (−) chronic HBV infection, HBeAg (−) chronic hepatitis B, and resolved HBV.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Lin, MJ., Su, TH., Chen, CC. et al. Diversity and composition of gut microbiota in healthy individuals and patients at different stages of hepatitis B virus-related liver disease. Gut Pathog 15, 24 (2023). https://doi.org/10.1186/s13099-023-00549-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13099-023-00549-w