Skip to main content

Relationship between gut microbiome diversity and hepatitis B viral load in patients with chronic hepatitis B



Hepatitis B virus (HBV) infection is associated with a reduced risk of developing dyslipidemia and non-alcoholic fatty liver diseases. Given that the gut microbiota plays a significant role in cholesterol metabolism, we compared the differences in gut microbial diversity and composition between HBV-infected and uninfected subjects.


A prospective case–control study was designed comprising healthy controls (group A) and HBV-infected individuals (group B) in a 1:1 ratio (57 participants each; total = 114). The patients in group B were divided into two subgroups according to their HBV DNA loads: B1 < 2000 IU/mL (N = 40) and B2 ≥ 2000 IU/mL (N = 17). In a pairwise comparison of HBV-infected individuals and controls, higher alpha diversity was noted in group B, and the difference was significant only in patients in group B1. Alloprevotella and Eubacterium coprostanoligenes were predominant in group B1 compared to the control, whereas the abundance of Bacteroides fragilis and Prevotella 2 was lower.


The gut microbiome in HBV-infected individuals with a low viral load is highly diverse and is dominated by specific taxa involved in fatty acid and lipid metabolism. To our knowledge, this is the first demonstration of a correlation between the presence of certain bacterial taxa and chronic HBV infection depending on the load of HBV DNA.


Chronic hepatitis B virus (HBV) infection is a notable health concern worldwide, especially in Asia and Africa, and is a major risk factor for fatality due to liver diseases such as cirrhosis and hepatocellular carcinoma (HCC). Although HBV is associated with several types of liver disease, previous cohort studies have provided significant evidence that HBV infections are linked to a reduction in serum cholesterol levels [1, 2], thereby reducing the risk of developing non-alcoholic fatty liver disease (NAFLD) [3]. HBV is considered a ‘metabolovirus,’ because the gene expression of HBV and the expression of key metabolic genes in hepatocytes were found to be similarly regulated [4, 5]. However, the underlying mechanisms of HBV-induced diseases currently unclear; therefore, no complete cure is available for chronic HBV infection [6]. A clearer understanding of the metabolic pathways involved in the pathological process of HBV infection may provide a new perspective on the potential novel targets of HBV treatment.

The gut microbiome, when considered a metabolic organ, has several roles in digestion, vitamin synthesis, immunomodulation, cardiovascular conditions, and the brain–gut axis. Consequently, it has an impact on the pathogenesis of gastrointestinal, hepatic, respiratory, cardiovascular, and endocrine diseases [7]. HBV infection has been found to alter dynamic changes in the gut microbiota [8, 9]. Given the significant role of the gut microbiota in cholesterol metabolism, it is possible that the effects of the altered gut microbiome in chronic HBV infection involve the regulation of microbe–host interactions at the gut interface. Here, we consider that the cholesterol-lowering properties of HBV infection may be attributable to the change of the gut microbiota by HBV infection. To evaluate the effect of the gut microbiome in chronic HBV infection on host metabolic health, we compared the differences in gut microbial diversity and composition between two groups of subjects: HBV-infected and uninfected subjects. Since gut microbiota is affected by different phases of disease in chronic hepatitis B and antiviral agents, we enrolled HBV-infected subjects in an inactive state who were not receiving antiviral treatment. We also investigated the differences in the gut microbiota with respect to viral loads among HBV-infected patients. A high viral load was defined as HBV DNA level ≥ 2000 IU/mL and a low viral load was defined as HBV DNA < 2000 IU/mL [10].


Demographic characteristics of the study subjects

A total of 134 individuals participated in the study. Faecal samples were collected from 114 subjects, of which 57 (50%) were healthy subjects (controls, group A) and 57 (50%) were subjects with chronic HBV (group B). We excluded 20 patients from the initial group of participants; 14 of these had not sent in their stool samples; three had incomplete information on HBV DNA, hepatitis B e antigen (HBeAg), or hepatitis B e antibody (HBeAb); and three showed seroclearance of hepatitis B surface antigen (HBsAg) (Fig. 1). Group B was divided into two subgroups according to their HBV DNA loads: 40 patients were had a low viral load of HBV DNA (< 2000 IU/mL, group B1) and 17 had a high viral load of HBV DNA (≥ 2000 IU/mL, group B2).

Fig. 1
figure 1

Enrolment of study participants

The baseline characteristics of the study subjects are summarised in Table 1. The mean age and body mass index of the study participants were 45.59 and 23.21, respectively. Males accounted for 52.6% of the study participants. Compared with group A, old age, male sex, and current smoking were positively associated with group B. Laboratory findings showed that white blood cell and platelet counts were negatively associated and elevated hepatic enzymes were positively associated with group B compared to that in group A. There was a decreasing tendency in the low-density lipoprotein-cholesterol (LDL-C) and high-density lipoprotein-cholesterol (HDL-C) levels of those with HBV infection compared to group A. When comparing groups A and B1, similar trends were observed as when comparing groups A and B, except for the presence of lower HDL-C in group B1. Further subgroup analysis in group B showed that subjects in group B1 were associated with old age and low HDL-C, compared to those in group B2.

Table 1 Demographic and clinical characteristics of study participants

Comparisons of gut microbial diversity among subjects with and without chronic HBV infections: alpha and beta diversity

We used 16S rRNA gene sequencing to determine the differences in gut microbiota composition according to HBV DNA load. The sequencing depth ranged from 21,459 to 102,066 reads per sample (mean = 47,752) and 1550 features in 114 subjects after contingency-based filtering of features.

Group B showed significantly higher alpha diversity than group A in both the non-phylogenetic and phylogenetic indices (Fig. 2A and Additional file 1: Table S1). To investigate the extent of similarity and differences among microbial communities, we used both non-phylogenetic and phylogenetic indices for beta diversity analysis. Considering the presence or absence of the features, there were significant differences in the gut microbial community between groups A and B (p = 0.004 in Jaccard dissimilarity and p = 0.012 in unweighted UniFrac distance). However, no significant differences between groups A and B were observed in terms of the abundance of the features (Table 2).

Fig. 2
figure 2

Alpha diversity. Comparison between subjects with and without HBV infection (A) and among three groups by hepatitis B viral load (B). *q < 0.05, **q < 0.01 (Mann–Whitney U test, Bonferroni correction). The boxes indicate the interquartile range. The median value is shown as a line within the box. Whiskers extend to the most extreme value within 1.5 × interquartile range

Table 2 Statistical significance among groups by HBV DNA load using distance matrices for beta-diversity

We then analysed the microbial diversity of each of the three groups: groups A, B1, and B2. The differences in alpha diversity among the three groups were mostly between groups A and B1, except for Pielou’s evenness (Fig. 2B and Additional file 1: Table S1). The pairwise test showed higher richness in group B1 than in group A (Bonferroni-corrected q = 0.039 in observed amplicon sequence variants (ASV) and q = 0.012 in Shannon’s index). Regarding phylogenetic diversity, Faith’s PD also revealed significantly higher diversity in group B1 than in group A (q = 0.039). When group B2 was compared with other groups, no statistically significant differences were observed in either non-phylogenetic or phylogenetic diversity. The beta-diversity analysis using permutational multivariate analysis of variance (PERMANOVA) revealed significant differences in microbial composition across the three groups in terms of Jaccard dissimilarity and unweighted UniFrac distance (p = 0.019 and p = 0.044, respectively; Table 2). The results of the pairwise PERMANOVA also demonstrated significant community differences between groups A and B1 (q = 0.030 in Jaccard dissimilarity and q = 0.024 in unweighted UniFrac distance), suggesting dissimilarities in the microbial composition of groups A and B1. However, owing to the high number of samples and inter-individual variation, the faecal microbiota in groups A and B1 could not be separated clearly by principal coordinate analysis (Additional file 1: Fig. S1); however, the differences in microbial community composition were significant between the two groups (Table 2). As in the results comparing groups A and B, the differences in bacterial communities among the three groups were also significant only in the indices of beta diversity in terms of the presence or absence of the features, whereas there was no significant difference in the indices regarding the abundance of features. Similar to the results of alpha diversity, the microbial community structure of group B2 was not significantly different from that of groups A or B1 (Table 2).

Composition of gut microbiota among subjects with and without HBV infection

To identify the microbial taxa associated with HBV load, we performed association analysis using Analysis of composition of microbiomes (ANCOM). In total, 1550 ASVs were mapped to 1265 distinct taxa, including 14 phyla, 23 classes, 45 orders, 94 families, 326 genera, and 751 species. Bacteriodetes and Firmicutes were the most dominant phyla, comprising 48.4 and 46.6% of the total sequences, respectively; however, there were no significant differences in the levels of these phyla among the groups.

We investigated the microbial composition of group A compared to that of group B. The ANCOM models showed several targets for HBV, including three genera (Alloprevotella, Paraprevotella, and Hungatella) and 17 species, after adjusting for age, sex, and smoking (Additional file 1: Table S2). We then compared the microbial compositions of groups A and B1 according to the results of the diversity analysis. There were no significant differences in the levels of phyla to order. However, significant differences between the two groups were found in two families, seven genera, and sixteen species (Additional file 1: Table S2). In Fig. 3, we present the relative abundances of the identified taxa. The genera Alloprevotella, Paraprevotella, Hungatella, Mitsuokella, and Family XIII AD3011 and an uncultured genus of the family Lachnospiraceae showed significant differences compared to other genera when groups A and B1 were compared (W ≥ 230; ANCOM). The W statistics indicate that a taxon has significantly different compositions between groups when compared to over 70% of the total taxa, such as 230 genera of the total 326 genera. DADA2 resolved fine-scale variation better than the operational taxonomic unit clustering method and accurately reconstructed amplicon-sequenced communities at the highest resolution [11]; however, the amplicon sequencing of the 16S rRNA gene was best at the genus-level resolution [12]. Although we identified the taxa associated with HBV at the species level, most species with high W values were unclassified bacteria belonging to several genera such as Alloprevotella, Paraprevotella, Family XIII AD3011, Ruminococcus, Hungatella, and Eubacterium coprostanoligenes (Additional file 1: Table S2), except for Bacteroides fragilis. Further, when group B2 was compared with groups A or B1, no statistically significant taxon was found in the analysis of the relative abundance of bacterial taxonomic groups.

Fig. 3
figure 3

Relative abundance of taxa associated with HBV infection. Bar plots showing the relative abundance of family (A), genus (B), and species (C) levels in health controls (group A) and HBV-infected individuals with a low viral load (group B1)

Based on the above results, we performed linear discriminant analysis effect size (LEfSe) to detect potential bacterial markers that most likely explain the differences between groups A and B1 by coupling standard tests for statistical significance with additional tests encoding biological consistency and effect relevance [13]. A comparison of groups A and B1 identified 16 taxa with linear discriminant analysis scores > 3 (Additional file 1: Fig. S2). The LEfSe results confirmed the significant enrichment of group Eubacterium coprostanoligenes in group B1 and found an enrichment of several species in Alloprevotella compared to the controls (group A). We also found that Prevotella 2 spp., Hungatella spp., and Bacteroides fragilis were more abundant in the control than in group B1.

Predicted functional microbiota in HBV-infected individuals with a low viral load (group B1)

To understand the gut microbial functions related to a low viral load (group B1), we used Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2 (PICRUSt2) to infer putative metagenomes from the 16S rRNA gene profiles. Statistical analysis of taxonomic and functional profiles (STAMP) was used to identify microbially relevant functions linked to HBV. However, after multiple comparison correction, we could not find any pathway that passed the significant threshold (Benjamini–Hochberg correction, FDR q < 0.05). Instead, Fig. 4 shows the suggested pathways with nominal significance (p < 0.05). The saturated fatty acid elongation pathway was lower in B1 than in A (p = 0.006), whereas peptidoglycan maturation, methylophosphonate degradation I, and the superpathway of thiamin diphosphate biosynthesis II pathways were enriched in B1 compared to A (p < 0.05).

Fig. 4
figure 4

Predicted metabolic pathways that changed in group B1 compared with group A. Bar plots on the left side display the mean proportion of each pathway. Dot plots on the right show the differences in mean proportions between the two indicated groups using p-value. The p-values were not corrected with any multiple comparison method


In this prospective case–control study, we found that chronic HBV infection might be associated with a decreasing tendency in cholesterol levels, compared to the control group. A similar trend was observed in the groups of healthy controls and HBV-infected patients with a low HBV DNA level of < 2000 IU/mL. In a pairwise comparison between HBV-infected individuals and healthy controls, higher alpha diversity was noted in HBV carriers with a low viral load (group B1). Significant enrichment of Alloprevotella spp. and the Eubacterium coprostanoligenes was found in the low HBV DNA group subjects (group B1) compared to the controls (group A), whereas Prevotella 2 spp., Hungatella spp., and Bacteroides fragilis were more abundant in the controls than in group B1. This study demonstrated a correlation between the presence of certain bacterial taxa and chronic HBV infection according to the load of HBV DNA.

Gut microbiota dysbiosis was noted in HBV-infected individuals, with differential microbial communities evident in the healthy controls compared with each clinical stage of HBV infection such as chronic hepatitis B, cirrhosis, and HCC [14]. Based on previous findings in cohort studies that beneficial effects on lipid metabolism were identified in healthy HBsAg-positive individuals compared to controls [1, 2], the study participants in this study were those with HBV infection who had not received antiviral treatment and excluding cases of active hepatitis, cirrhosis, or HCC. It is of interest that HBV-infected patients with a low viral load of < 2000 IU/mL had a higher diversity in the gut microbiome compared to the healthy controls and patients with a high viral load of ≥ 2000 IU/mL. Patients with inactive chronic HBV infection, as indicated by HBeAg-negative, HbeAb-positive, persistently normal ALT, and a low viral load of < 2000 IU/mL, have an excellent long-term prognosis with low rates of histological progression, high rates of HbsAg clearance, and a very low risk of cirrhosis and HCC [15]. Considering that an inactive HBV carrier state is related to a better prognosis against progression to cirrhosis or HCC, it is notable that the gut microbiome in patients with HBV infection with low viral load has different properties, such as high diversity and microbial composition associated with lipid metabolism, suggesting that the gut microbiome in an inactive state of HBV infection may interplay between viral replication and the host immune and metabolic responses. Further investigations are needed to evaluate the relationship between viral factors and metabolic profiles in patients with chronic HBV.

Alloprevotella was positively associated with HBV-infected individuals with a low viral DNA load in comparison with healthy controls. In a previous study, Alloprevotella was negatively associated with lifetime cardiovascular disease risk, but Prevotella 2 and Prevotella 7 were positively associated with lifetime cardiovascular disease risk [16]. In another study, exposure to a high dose of diethyl phosphate enriched butyrate-producing genera such as Alloprevotella, leading to an increase in estradiol and a resulting decrease in triglycerides and LDL-C [17]. These bacteria are short-chain fatty acid producers and are negatively correlated with NAFLD [18]. Short chain fatty acids can be digested by the intestine and indirectly regulate energy metabolism and insulin sensitivity through specific receptors. Considering that Alloprevotella is associated with a reduced risk of cardiovascular disease and NAFLD, its predominance in HBV-infected subjects with a low viral DNA level might be associated with a decreased incidence of NAFLD in this population.

Eubacterium coprostanoligenes are anaerobic, non-sporing, gram-positive rod-shaped bacteria that reduce approximately 90% of the cholesterol in growth media to coprostanol [19, 20]. Unlike cholesterol, coprostanol is poorly absorbed by the human intestine, suggesting that it may have an impact on cholesterol metabolism and modulation of serum cholesterol levels [21]. It has been reported that most strains assigned to the genus Eubacterium exhibit cholesterol-reducing properties [22]. In this study, Eubacterium coprostanoligenes were predominantly found in HBV-infected individuals with low viral DNA loads, compared to the controls, which might mean that it is the main contributor to the cholesterol-lowering effects in HBV-infected individuals.

It is important to understand the role of the gut microbiome in the interactions between HBV infection and metabolic profiles. We found that the inferential pathway of saturated fatty acid elongation was enriched in non-hepatitis B individuals compared to that in hepatitis B-infected subjects with a low DNA load. There is evidence that long-chain saturated fatty acids activate certain transcription factors that target lipogenic target genes, sometimes leading to negative health effects. During type 2 diabetes progression, long-chain saturated fatty acids trigger cell apoptosis and increase oxidative stress [23]. Chronic exposure to long-chain saturated fatty acids is a risk factor for beta-cell apoptosis. The specific toxic effects of long-chain fatty acids may be related to ceramide formation and persistent endoplasmic reticulum stress. It is well established that excessive levels are contributory factors in the development of obesity, type 2 diabetes, and cardiovascular disease [24, 25]. These findings suggest that gut microbiota in inactive HBV carriers, which are characterized by less saturated fatty acid production, might play a major role in preventing hyperlipidemia compared to those in the controls.

Our study has several limitations. First, we excluded chronic hepatitis B patients who were receiving antiviral treatments such as nucleos(t)ide analogues because of active hepatitis with increased alanine aminotransferase (ALT) levels and high levels of HBV DNA. Antiviral treatment generally decreases HBV DNA to very low or undetectable levels. Therefore, a reduced viral level in patients undergoing antiviral therapy makes it difficult to assess the direct relationship between HBV DNA levels and the diversity of gut microbiomes. Second, although TG contributes to one component of metabolic syndrome, triglyceride levels were not measured in this study. Instead, we compared other lipid profiles, including total, LDL, and HDL cholesterol levels, among the study subgroups. Third, the sample size of group B2 was small compared to those of the other groups; therefore, we cannot exclude the possibility that the null results for group B2 might have been affected by the statistical power. Therefore, further investigation is warranted on the role of the gut microbiota in individuals with high HBV DNA load. Finally, the current study used 16S methods which cannot provide the taxonomic resolution at strain and information on several bacterial genes. Whole metagenomic sequencing might allow a wider range of understanding related to definitive taxonomic classification and functional profiling.


The gut microbiome in HBV-infected individuals with a low viral load is highly diverse and is dominated by specific taxa involved in fatty acid and lipid metabolism. To the best of our knowledge, this is the first demonstration of a correlation between the presence of certain bacterial taxa and chronic HBV infection according to the load of HBV DNA.


Study participants and study design

The study included participants aged 18 years or older with at least one follow-up visit, who underwent serologic testing for HBsAg and lipid profiles, including total cholesterol, LDL-C, and HDL-C levels, during the study period. A prospective case–control study was designed with participants in a 1:1 ratio: healthy controls (N = 57) who underwent comprehensive health check-ups and patients (N = 57) with chronic HBV infection who visited the outpatient department of Gastroenterology and Hepatology at Kangbuk Samsung Hospital. Faecal samples were collected from the study participants for 7 months between 1 November 2018 and 31 May 2019.

We excluded subjects with previous or current use of antiviral agents for HBV infection, antibiotics, probiotics, proton-pump inhibitors, or cholesterol-lowering medication, and those with diabetes mellitus. In the case of patients with chronic HBV infection, any subject with liver cirrhosis or cancer related to HBV infection was also excluded.

Chronic HBV infection and group definition.

Chronic HBV infection was defined as HBsAg-positive for > 6 months. HBsAg-negative participants were defined as healthy controls. The healthy subjects (controls) were designated as group A, whereas the patients with HBV infection were designated as group B. In the subjects with chronic HBV infection (group B), additional laboratory tests were also performed to evaluate HBV DNA levels, HBeAg, and HBeAb. This study included HBV-infected subjects with ALT levels < 2 upper limit of normal (80 IU/mL). To further distinguish inactive carriers from reactivated patients, we used a HBV DNA level of 2000 IU/mL to categorise HBV-infected subjects into two subgroups according to the viral load: low viral load of HBV DNA < 2000 IU/mL (group B1) and high viral load of HBV DNA ≥ 2000 IU/mL (group B2) [26, 27].

This study was approved by the Institutional Review Board of Kangbuk Samsung Hospital (KBSMC 2018–07-027, registered in July 2018). All study participants provided written informed consent to participate in the study. The study was conducted in accordance with the ethical guidelines of the World Medical Association Declaration of Helsinki.


We collected data on the medical history, medication use, smoking habits, and sociodemographic characteristics of the study subjects. Body mass index was calculated as the weight in kilogrammes divided by the height in meters squared. A complete blood count was performed using a haematology analyser XE‐2100 (Sysmex, Kobe, Japan). Serum ALT, aspartate aminotransferase (AST), total cholesterol, LDL-C, and HDL-C levels were determined using an automated chemistry analyser (Cobas c702; Roche Diagnostics, Tokyo, Japan). Haemoglobin A1c (HbA1c) levels were measured using an immunoturbidimetry immunoassay (Cobas c513; Roche Diagnostics, Tokyo, Japan). HBeAg and HBeAb were determined using an electrochemiluminescence immunoassay (Cobas e801; Roche Diagnostics, Tokyo, Japan) based on the sandwich principle. HBV DNA was quantified via real-time PCR using the COBAS AmpliPrep-COBAS TaqMan HBV assay (Roche Molecular Systems, Basel, Switzerland).

DNA extraction from faecal samples and sequencing of the 16S rRNA gene

The OMNIgene-GUT collection kit (OMR-200, DNA Genotek, Ottawa, Canada) was used to collect stool samples. DNA extraction from faecal samples was performed within 1 month of storage using the MOBio PowerSoil DNA Isolation Kit (MO BIO Laboratories, Carlsbad, CA, USA) according to the manufacturer’s instructions. Amplification and sequencing were performed to analyse the bacterial communities as described previously [28]. Genomic DNA was amplified using fusion primers targeting the V3 and V4 regions of the 16S rRNA gene. Samples were pooled for sequencing on the Illumina MiSeq platform (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions [29]. The DADA2 [11] plugin of the QIIME2 package (version 2019.7, [30] was used to perform sequence quality control and construct a feature table of ASVs. For taxonomic structure analysis, taxonomy was assigned to ASVs using a pre-trained naïve Bayes classifier and the q2-feature-classifier plugin with the database Silva 132 release in the QIIME2 package.

16S rRNA gene compositional analysis and statistical analysis

Basic statistical analyses were performed using RStudio version 0.98.983 (Boston, MA, USA). For diversity analysis, the feature table was rarefied to 21,459 sequences per sample by random subsampling in QIIME2. To evaluate alpha diversity, we computed the number of ASVs observed in each sample, Shannon’s index, Pielou’s evenness, and Faith’s phylogenetic diversity (PD). The Mann–Whitney U test was used to test the differences between pairwise groups. The Kruskal–Wallis test was used as a non-parametric statistical test to compare the three groups. To measure beta diversity, we used the UniFrac distance [31] to estimate the dissimilarity among group members by incorporating the phylogenetic distances between ASVs. The unweighted and weighted UniFrac distances were calculated for the presence/absence and abundance of ASVs, respectively. Non-phylogenetic methods were also used with Bray–Curtis dissimilarities and Jaccard dissimilarity for the presence/absence and abundance of ASVs, respectively. Pairwise PERMANOVA with 999 random permutations was used to test the significance between groups.

Analysis of the composition of the microbiome via the ANCOM test [32] was performed to determine whether there were significant differences in the relative abundance of any taxa among the groups across multiple taxonomic levels. The final significance was expressed in the empirical distribution of W from the analyses of the two groups. We also conducted LEfSe analysis to detect potential HPV-specific bacterial markers [13].

For the functional inferences of the microbial community, we conducted the PICRUSt2 (v2.2.0-b) [33] with ASVs according to a recent manual PICRUSt2 predictions based on the following gene families are supported by enzyme classification numbers (EC numbers, as of 21 January 2016). We generated PICRUSt2 EC gene family and metabolic pathway database (Metacyc) abundance predictions [34]. The results were visualised using STAMP version 2.1.3 [35] and tested using Welch’s t-test for two groups. All predictions were adjusted by multiple testing correction (Benjamini–Hochberg correction, FDR q-value < 0.05).

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.



Hepatitis B virus


Hepatocellular carcinoma


Non-alcoholic fatty liver disease


Low-density lipoprotein-cholesterol


High-density lipoprotein-cholesterol


Hepatitis B e antigen


Hepatitis B e antibody


Hepatitis B surface antigen


Amplicon sequence variants


Permutational multivariate analysis of variance


Analysis of composition of microbiomes


Linear discriminant analysis effect size


Statistical analysis of taxonomic and functional profiles


False positive rate


Alanine aminotransferase


Aspartate aminotransferase


Haemoglobin A1c


  1. 1.

    Joo EJ, Chang Y, Yeom JS, Cho YK, Ryu S. Chronic hepatitis B virus infection and risk of dyslipidaemia: a cohort study. J Viral Hepat. 2019;26(1):162–9.

    CAS  Article  Google Scholar 

  2. 2.

    Chen JY, Wang JH, Lin CY, Chen PF, Tseng PL, Chen CH, et al. Lower prevalence of hypercholesterolemia and hyperglyceridemia found in subjects with seropositivity for both hepatitis B and C strains independently. J Gastroenterol Hepatol. 2010;25(11):1763–8.

    Article  Google Scholar 

  3. 3.

    Joo EJ, Chang Y, Yeom JS, Ryu S. Hepatitis B virus infection and decreased risk of nonalcoholic fatty liver disease: a cohort study. Hepatology. 2017;65(3):828–35.

    CAS  Article  Google Scholar 

  4. 4.

    Geier A. Hepatitis B virus: the “metabolovirus” highjacks cholesterol and bile acid metabolism. Hepatology. 2014;60(5):1458–60.

    CAS  Article  Google Scholar 

  5. 5.

    Jarcuska P, Drazilova S, Fedacko J, Pella D, Janicko M. Association between hepatitis B and metabolic syndrome: current state of the art. World J Gastroenterol. 2016;22(1):155–64.

    CAS  Article  Google Scholar 

  6. 6.

    Shi YX, Huang CJ, Yang ZG. Impact of hepatitis B virus infection on hepatic metabolic signaling pathway. World J Gastroenterol. 2016;22(36):8161–7.

    CAS  Article  Google Scholar 

  7. 7.

    O’Hara AM, Shanahan F. The gut flora as a forgotten organ. EMBO Rep. 2006;7(7):688–93.

    CAS  Article  Google Scholar 

  8. 8.

    Zhu Q, Xia P, Zhou X, Li X, Guo W, Zhu B, et al. Hepatitis B virus infection alters gut microbiota composition in mice. Front Cell Infect Microbiol. 2019;9:377.

    CAS  Article  Google Scholar 

  9. 9.

    Yun Y, Chang Y, Kim HN, Ryu S, Kwon MJ, Cho YK, et al. Alterations of the Gut microbiome in chronic hepatitis B virus infection associated with alanine aminotransferase level. J Clin Med. 2019;8(2):173.

    CAS  Article  Google Scholar 

  10. 10.

    Croagh CM, Lubel JS. Natural history of chronic hepatitis B: phases in a complex relationship. World J Gastroenterol. 2014;20(30):10395–404.

    Article  Google Scholar 

  11. 11.

    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.

    CAS  Article  Google Scholar 

  12. 12.

    Poretsky R, Rodriguez RL, Luo C, Tsementzi D, Konstantinidis KT. Strengths and limitations of 16S rRNA gene amplicon sequencing in revealing temporal microbial community dynamics. PLoS ONE. 2014;9(4):e93827.

    Article  Google Scholar 

  13. 13.

    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.

    Article  Google Scholar 

  14. 14.

    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.

    CAS  Article  Google Scholar 

  15. 15.

    Invernizzi F, Vigano M, Grossi G, Lampertico P. The prognosis and management of inactive HBV carriers. Liver Int. 2016;36:100–4.

    CAS  Article  Google Scholar 

  16. 16.

    Kelly TN, Bazzano LA, Ajami NJ, He H, Zhao J, Petrosino JF, et al. Gut Microbiome associates with lifetime cardiovascular disease risk profile among bogalusa heart study participants. Circ Res. 2016;119(8):956–64.

    CAS  Article  Google Scholar 

  17. 17.

    Yang F, Li J, Pang G, Ren F, Fang B. Effects of diethyl phosphate, a non-specific metabolite of organophosphorus pesticides, on serum lipid, hormones, inflammation, and gut microbiota. Molecules. 2019;24(10):2003.

    CAS  Article  Google Scholar 

  18. 18.

    Shang Q, Song G, Zhang M, Shi J, Xu C, Hao J, et al. Dietary fucoidan improves metabolic syndrome in association with increased Akkermansia population in the gut microbiota of high-fat diet-fed mice. J Funct Foods. 2017;28:138–46.

    CAS  Article  Google Scholar 

  19. 19.

    Kenny DJ, Plichta DR, Shungin D, Koppel N, Hall AB, Fu B, et al. Cholesterol metabolism by uncultured human gut bacteria influences host cholesterol level. Cell Host Microbe. 2020;28(2):245-57 e6.

    CAS  Article  Google Scholar 

  20. 20.

    Mukherjee A, Lordan C, Ross RP, Cotter PD. Gut microbes from the phylogenetically diverse genus Eubacterium and their various contributions to gut health. Gut Microbes. 2020;12(1):1802866.

    Article  Google Scholar 

  21. 21.

    Gerard P. Metabolism of cholesterol and bile acids by the gut microbiota. Pathogens. 2013;3(1):14–24.

    Article  Google Scholar 

  22. 22.

    Kriaa A, Bourgin M, Mkaouar H, Jablaoui A, Akermi N, Soussou S, et al. Microbial reduction of cholesterol to coprostanol: an old concept and new insights. Catalysts. 2019;9(2):167.

    Article  Google Scholar 

  23. 23.

    Nesterova AP, Klimov EA, Zharkova M, Sozin S, Sobolev V, Ivanikova NV, et al. Endocrine, nutritional, and metabolic diseases. In: Nesterova AP, Klimov EA, Zharkova M, Sozin S, Sobolev V, Ivanikova NV, et al., editors. Disease pathways: an atlas of human disease signaling pathways. 1st ed. Amsterdam: Elsevier; 2019. p. 121–218.

    Google Scholar 

  24. 24.

    van Dijk SJ, Feskens EJ, Bos MB, Hoelen DW, Heijligenberg R, Bromhaar MG, et al. A saturated fatty acid-rich diet induces an obesity-linked proinflammatory gene expression profile in adipose tissue of subjects at risk of metabolic syndrome. Am J Clin Nutr. 2009;90(6):1656–64.

    Article  Google Scholar 

  25. 25.

    French MA, Sundram K, Clandinin MT. Cholesterolaemic effect of palmitic acid in relation to other dietary fatty acids. Asia Pac J Clin Nutr. 2002;11(Suppl):S401–7.

    CAS  Article  Google Scholar 

  26. 26.

    Invernizzi F, Vigano M, Grossi G, Lampertico P. The prognosis and management of inactive HBV carriers. Liver Int. 2016;36(Suppl 1):100–4.

    CAS  Article  Google Scholar 

  27. 27.

    Liu J, Yang HI, Lee MH, Jen CL, Batrla-Utermann R, Lu SN, et al. Serum levels of hepatitis B surface antigen and DNA can predict inactive carriers with low risk of disease progression. Hepatology. 2016;64(2):381–9.

    CAS  Article  Google Scholar 

  28. 28.

    Kim HN, Yun Y, Ryu S, Chang Y, Kwon MJ, Cho J, et al. Correlation between gut microbiota and personality in adults: a cross-sectional study. Brain Behav Immun. 2018;69:374–85.

    Article  Google Scholar 

  29. 29.

    Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol. 2013;79(17):5112–20.

    CAS  Article  Google Scholar 

  30. 30.

    Bolyen E, Rideout J, Dillon M, Bokulich N, Abnet C, Al-Ghalith G, et al. Reproducible, interactive, scalable, and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):852–7.

    CAS  Article  Google Scholar 

  31. 31.

    Lozupone C, Lladser ME, Knights D, Stombaugh J, Knight R. UniFrac: an effective distance metric for microbial community comparison. ISME J. 2011;5(2):169–72.

    Article  Google Scholar 

  32. 32.

    Mandal S, Van Treuren W, White RA, Eggesbo M, Knight R, Peddada SD. Analysis of composition of microbiomes: a novel method for studying microbial composition. Microb Ecol Health Dis. 2015;26:27663.

    PubMed  Google Scholar 

  33. 33.

    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.

    CAS  Article  Google Scholar 

  34. 34.

    Caspi R, Billington R, Fulcher CA, Keseler IM, Kothari A, Krummenacker M, et al. The MetaCyc database of metabolic pathways and enzymes. Nucleic Acids Res. 2018;46(D1):D633–9.

    CAS  Article  Google Scholar 

  35. 35.

    Parks DH, Tyson GW, Hugenholtz P, Beiko RG. STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics. 2014;30(21):3123–4.

    CAS  Article  Google Scholar 

Download references


We thank all the patients who participated in this study.


This study was supported by a grant from the Basic Science Research Program through the National Research Foundation of Korea (NRF), funded by the Ministry of Education (NRF-2018R1C1B6009358).

Author information




EJJ, HNK, and YKC were contributed to the concept and design of the study. MJK carried out laboratory processing. HSC and WS interpretated the data. EJJ and HNK contributed data analysis and interpretation and drafted the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Han-Na Kim or Yong Kyun Cho.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Institutional Review Board of Kangbuk Samsung Hospital (KBSMC 2018–07-027, registered in July 2018). All study participants provided written informed consent to participate in the study. The study was conducted in accordance with the ethical guidelines of the World Medical Association Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no conflict of interest.

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

. Alpha diversity among groups of A1, B1 and B2. Table S2. Detection of differentially abundant taxa in all HBV carriers and HBV carriers with low viral DNA level compared to healthy controls. Fig. S1. Principal coordinates analysis of beta-diversity indices. Fig. S2. Low HBV DNA group-specific biomarkers.

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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Joo, EJ., Cheong, H.S., Kwon, MJ. et al. Relationship between gut microbiome diversity and hepatitis B viral load in patients with chronic hepatitis B. Gut Pathog 13, 65 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Hepatitis B
  • Chronic
  • Viral load
  • Gastrointestinal microbiome