Abstract

Background and Aim. Dysbiosis of gut microbiota is important in the development of hepatocellular carcinoma (HCC). However, little is known about whether and how dysbiosis impacts HCC progression. This cross-sectional study is aimed at evaluating microbiome dysbiosis, gut damage, and microbial translocation in different stages of HCC. Method. This study included 74 Chinese male patients with HCC. They were divided into early (), intermediate (), and terminal () groups, referred to as Barcelona Clinic Liver Cancer stage 0+A, B, and C+D, respectively. Paired fecal and plasma samples were collected. Microbial composition and profiles were analyzed by 16S rRNA gene sequencing. The levels of gut damage marker (regenerating islet-derived protein 3α (REG3α)) and microbial translocation markers (soluble CD14 (sCD14), lipopolysaccharide-binding protein (LBP), peptidoglycan recognition proteins (PGRPs)) were determined in plasma samples of patients by ELISA. Twenty plasma cytokine and chemokines were determined by Luminex. Results. In early, intermediate, and terminal groups, the abundance of the Bifidobacteriaceae family decreased significantly (3.52%, 1.55%, and 0.56%, respectively, ), while the abundance of the Enterococcaceae family increased significantly (1.6%, 2.9%, and 13.4%, respectively, ). Levels of REG3α and sCD14 were markedly elevated only in the terminal group compared with the early ( and ) and intermediate groups ( and ). The level of LBP significantly increased in the intermediate () and terminal () groups compared with the early group. The PGRP levels were elevated only in the terminal group compared with the early group (). The ratio of Enterococcaceae to Bifidobacteriaceae was significantly associated with the levels of REG3α, LBP, sCD14, and PGRPs. With HCC progression, increased levels of inflammatory cytokines accompanied by a T cell-immunosuppressive response and microbial translocation were observed. Conclusion. Gut microbiota compositional and functional shift, together with elevated gut damage and microbial translocation, may promote HCC development by stimulating inflammatory response and suppressing T cell response.

1. Introduction

Hepatocellular carcinoma (HCC) is the third leading cause of cancer mortality worldwide [1]. However, the underlying mechanism is poorly understood. Recently, increasing evidence has indicated the relationship between dysbiosis of gut microbiome and HCC development [24]. However, there has been only one report in which the abundance of proinflammatory bacteria Enterobacteriaceae tended to increase parallelly with the progression of HCC [5]. The findings demonstrated that microbiota dysbiosis could regulate the innate immune system and advance liver disease [68].

Besides microbiome dysbiosis, gut permeability was also found to correlate with HCC. A recent study reported that gut permeability was significantly higher in patients with HCC compared with healthy individuals [9]. The effect of microbial translocation on cirrhosis was well studied [10], but was not researched in HCC. Therefore, this cross-sectional study is aimed at evaluating microbiome dysbiosis, gut damage, and microbial translocation in 74 male patients in different stages of HCC. The levels of cytokines, which reflect inflammation and immune response, were also measured. It was speculated that dysbiosis together with increased microbial translocation might impact HCC progression by regulating systemic inflammation and immune response.

2. Materials and Methods

2.1. Patient Selection

Patients with HCC admitted to Beijing Youan Hospital from June 2018 to September 2019 were recruited in this study. The study was approved by the ethics committee of Beijing Youan Hospital, Capital Medical University (No. 2018-038). All patients provided written informed consents.

The diagnosis of HCC was made following the international guidelines [11]. The inclusion criteria were as follows: HCC; years; male sex; and hepatitis B virus (HBV) surface antigen-positive or anti-hepatitis C virus (anti-HCV) antibody positive. The exclusion criteria were as follows: patients suffering from other causes of liver diseases, such as nonalcoholic fatty liver disease, alcoholic liver disease, and cholestatic liver disease. Patients with irritable bowel syndrome or inflammatory bowel disease; other cancers besides HCC; and autoimmune disease or serious cardiac, kidney, and respiratory diseases were also excluded. No antibiotic, probiotic, and immunosuppressive drugs were taken within 2 weeks before sample collection.

HCC was staged according to the Barcelona Clinic Liver Cancer (BCLC) staging system [12]. Stage 0 referred to single and Child-Pugh A. Stage A was defined as single or three and Child-Pugh A or B. Stage B referred to patients with multinodules and Child-Pugh A or B. Stage C comprised patients with portal invasion or extrahepatic spread and Child-Pugh A or B. Stage D included patients with HCC and Child-Pugh C. In the present study, patients with HCC were divided into early, intermediate, and terminal groups, which referred to stage 0+A, B, and C+D, respectively.

2.2. Fecal Sample Collection, DNA Extraction, and Polymerase Chain Reaction Sequencing

Each fresh fecal sample was split into two tubes containing bacterial RNA LOCKER (Youkang, Nanjing, China) and stored at -80°C. Total bacterial DNA was isolated from stool samples using the QIAamp Fast DNA Stool Mini Kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions and diluted to 1 ng/μL using sterile water. The V4-V5 regions of the prokaryotic 16S rRNA gene were amplified using the universal primer pair 515F (5-GTGYCAGCMGCCGCGGTA-3) and 909R (5-CCCCGYCAATTCMTTTRAGT-3) with barcode and then sequenced and analyzed [13]. The reason to choose the V4-V5 regions was that we also wanted to detect low abundant bacteria in the samples [5, 14]. Paired stool samples from 20 patients were sequenced to validate the results. Sequencing libraries were generated using TruSeq DNA Polymerase Chain Reaction- (PCR-) Free Sample Preparation Kit (Illumina, USA), and index codes were added. The amplification products were sequenced on an Illumina HiSeq 2500 platform, and 250 bp paired-end reads were generated. Sequencing data were analyzed using Quantitative Insights into Microbial Ecology (QIIME) platform version 1.9 and R v3.3.1 [15]. The paired-end reads from the DNA fragments were merged using FLASH (v1.2.7) [16].

2.3. Operational Taxonomy Unit Clustering and Taxonomy Annotation

Sequences with more than 97% similarity were assigned to the same operational taxonomy units (OTUs). OTU classification was conducted by running a Basic Local Alignment Search Tool (BLAST) search against the Greengenes database using the representative sequence set as a query. A representative sequence for each OTU was screened for further annotation. The sequences were annotated using RDP classifier V.2.6 according to the developer’s documents (http://rdp.cme.msu.edu/classifier/class_help.jsp#conf) [17].

2.4. Bacterial Diversity Analysis

Bacterial diversity was determined by sampling-based OTU analysis. It is presented by Observed species index, Chao1 index, Shannon index, and Simpson index, all of which were calculated using the R vegan program package [18].

2.5. Structure of Microbial Communities and OTU Biomarker Identification

Principal component analysis (PCA) was conducted by R package (http://www.R-project.org/). The weighted and unweighted UniFrac distances were calculated with the phyloseq package [19]. The linear discriminant analysis (LDA) effect size (LEfSe) model was used to identify differences in microbiota composition for phylotypes [20]. Based on the normalized relative abundance matrix, taxa with significantly different abundances were determined by LEfSe using Kruskal-Wallis rank sum test.

2.6. Determination of Plasma Biomarkers of Gut Damage and Microbial Translocation

The plasma marker of gut damage, regenerating islet-derived protein (REG3α), was determined using the ELISA kit (R&D Systems, ON, Canada). Soluble CD14 (sCD14), lipopolysaccharide-binding protein (LBP), and peptidoglycan recognition proteins (PGRPs) were selected as microbial translocation markers. The ELISA kits were Human CD14 Quantikine ELISA kit (Bio-Techne Ltd., Abingdon, United Kingdom), Human LBP DuoSet ELISA kit (Bio-Techne Ltd.), and Human PGRPs ELISA kit (Thermo Fisher, Merelbeke, Belgium). All assays were performed in duplicate following the manufacturer’s protocol.

2.7. Determination of Plasma Cytokines and Chemokines

An inflammation 20-Plex Human ProcartaPlex Panel (EPX200-12185-901, Affymetrix eBioscience, Vienna, Austria) was used to detect the expression of 20 cytokines and chemokines. They were granulocyte colony-stimulating factor (GM-CSF), intercellular cell adhesion molecule-1 (ICAM-1), interferon-α (IFN-α), IFN-γ, interleukin-1α (IL-1α), IL-1β, IL-4, IL-6, IL-8, IL-10, IL-12p70, IL-13, IL-17A, interferon-induced protein-10 (IP-10), monocyte chemoattractant protein-1 (MCP-1), macrophage inflammatory protein 1α (MIP-1α), MIP-1β, sE-selectin, sP-selectin, and tumor necrosis factor-α (TNF-α).

2.8. Statistical Analysis

One-way analysis of variance was used to evaluate the differences among the three groups. Continuous variables were compared using the Wilcoxon rank sum test between the two groups. Fisher’s exact test compared categorical variables. Correlations were performed using a nonparametric Spearman test. A two-sided value < 0.05 indicated a significant difference. The threshold logarithmic LDA score for discriminative features was 2. Statistical analyses were conducted using GraphPad Prism 6.0 (La Jolla, CA, USA). Multivariate analysis was performed using SPSS 24.0 (IBM SPSS). Correlations analysis was performed by R i386 software (v.4.0.5).

In addition, an index was introduced to measure the degree of dysbiosis. This index was calculated based on the ratio of the relative abundance of Enterococcaceae family to that of Bifidobacteriaceae family. The ratio was calculated as follows:

3. Results

3.1. Patient Characteristics

Altogether, 74 male patients with HCC were included. The average age was 61.00 (53.00–65.00) years. Of them, the early group included 7 patients with stage 0 and 12 with stage A HCC; the intermediate group included 37 patients with stage B HCC; and the terminal group included 8 patients with stage C and 10 with stage D HCC. Among them, 37 patients had HBV and the others had HCV. The demographics and laboratory results are listed in Table 1 and Supplementary Table 1.

3.2. Gut Microbial Diversity Decreased with HCC Progression

The OTU number in early, intermediate, and terminal groups was , , and , respectively. Observed species and Chao1 indexes decreased significantly from early to intermediate and terminal groups (Observed species, and 0.038; Chao1, and 0.042; Figures 1(a) and 1(b)), while evenness reduced (Shannon, and 0.039; Simpson, and 0.028; Figures 1(c) and 1(d)). The results suggested that the alterations of OTUs were evenly contributed with HCC progression.

3.3. Gut Microbial Profiles and Compositions Shifted with HCC Progression

The PCA of the β-diversity index allowed a separation between any two groups, suggesting that the bacterial profiles of each group were distinguished (Figure 1(e)). The different taxa among the three groups identified by LDA analysis are summarized in Figure 1(f). And the ratio of Enterococcaceae to Bifidobacteriaceae family significantly increased from the early group to the terminal group (1.5 vs. 9.6, ; Figure 1(g)). The abundance of Actinobacteria phylum (4.46% vs. 0.88%, ; Figure 2(a)), Bifidobacteriaceae family (3.52% vs. 0.56%, ; Figure 2(b)), and Bifidobacterium genus (3.51% vs. 0.56%, ; Figure 2(c)) significantly decreased in the terminal group compared with the early group. Meanwhile, the abundance of Enterococcaceae (1.6% vs. 13.4%, ), Enterococcus genus (1.6% vs. 13.4%, ), and Enterobacteriaceae family (7.7% vs. 10.3%, ) markedly increased with disease progression. Lachnospiraceae, Peptostreptococcaceae, unidentified Clostridiales, Coriobacteriaceae, and Christensenellaceae families were enriched in the early group; however, they significantly decreased or were even undetectable in the terminal group (all values < 0.05). Correspondingly, at the genus level, five genera belonging to the Lachnospiraceae family (Blautia, Fusicatenibacter, Agathobacter, Anaerostipes, and Dorea), one genus belonging to the Clostridiales family (unidentified Clostridiales), and two genera belonging to the Peptostreptococcaceae family (Romboutsia and Intestinibacter) were also enriched in the early group and decreased or were even undetectable in the terminal group. All the alterations of the short-chain fatty acid- (SCFA-) producing bacteria mentioned earlier contributed to the decrease in richness along with HCC progression.

3.4. Gut Damage Increased in Patients in the Terminal Group

REG3α is a well-known marker of gut damage [21]. The plasma levels of REG3α significantly increased in patients in the terminal group () than in the early (, ) and intermediate groups (, ; Figure 3(a)).

3.5. Plasma Levels of Microbial Translocation Markers Elevated in Patients in the Terminal Group

Plasma sCD14 elevated in the terminal group () compared with the early (, ) and intermediate groups (, ; Figure 3(b)). The LBP level elevated in the intermediate () and terminal groups () compared with the early group (, and 0.025; Figure 3(c)). The PGRP level significantly increased in the terminal group () compared with the early group (, ; Figure 3(d)).

3.6. Cytokine and Chemokine Levels

The T cell immune response-related cytokines and chemokines, including IFN-γ, IL-4, IL-12p70, IL-13, and IL-17A, showed a decreased tendency from the early to the intermediate and terminal groups, but with no significant differences. The MCP-1 level decreased significantly in the terminal group (74.69 ng/mL) compared with the early (108.26 ng/mL, ) and intermediate groups (90.70 ng/mL, ; Figure 4(a)). The levels of chemokines, including GM-CSF, sE-selectin, and sP-selectin, tended to decrease in the terminal group without significant differences. The levels of other proinflammatory cytokines, including IL-1α, IL-6, IL-10, IP-10, MIP-1α, and MIP-1β, increased in the terminal group, but showed no significant difference. The IL-8 level significantly increased in the terminal group (16.31 ng/mL) compared with the early (7.42 ng/mL, ) and intermediate groups (6.80 ng/mL, ; Figure 4(b)). These results suggested that HCC progression was characterized by an elevated inflammatory response and suppressed T cell immune responses. The plasma levels of 20 cytokines and chemokines are summarized in Supplementary Table 2.

3.7. Relationship between HCC Progression, Gut Damage, Microbial Translocation, and Inflammatory Cytokines in Patients with HCC

The IL-6 and IL-8 levels positively correlated with markers of microbial translocation, but did not correlate with REG3α or bacterial changes. The IL-6 level positively correlated with the sCD14 level (, ) and the LBP level (, ; Figures 5(a) and 5(b)). They were markers of Gram-negative bacterial translocation. The IL-8 level correlated with the sCD14 (, ) and PGRP (, ; Figures 5(c) and 5(d)) levels, which was a marker of Gram-positive bacterial translocation. The MCP-1 level inversely correlated with the LBP (, ) and PGRP levels (, ; Figures 5(e) and 5(f)). No correlation was found between plasma TNF-α level and any of the aforementioned indexes.

3.8. The Ratio of Enterococcaceae to Bifidobacteriaceae Correlated with the Levels of REG3α and Markers of Bacterial Translocation

Both α-diversity and β-diversity could distinguish terminal HCC groups from the other two groups, so did gut damage and microbial translocation markers. Hence, whether gut-associated microbiota changes correlated with the loss of gut barrier integrity and circulating microbial translocation during HCC progression was evaluated. The LDA level revealed a decreased abundance of Bifidobacteriaceae family and increased abundance of Enterococcaceae family along with HCC progression. The ratio of Enterococcaceae to Bifidobacteriaceae was associated with the expression of gut damage marker REG3α (, ) and bacterial translocation markers, sCD14 (, ), LBP (, ), and PGRP (, ; Table 2). However, the ratio did not correlate with levels of IL-6 and IL-8.

4. Discussion

The findings of this study first showed that alterations in gut microbiota, elevated gut damage, and bacterial translocation were associated with HCC progression. Gut dysbiosis of HCC was characterized by increased abundance of Enterococcaceae and Enterobacteriaceae and decreased abundance of Bifidobacteriaceae and SCFA-producing bacteria. Elevation of inflammatory response and suppression of T cell immunity induced by microbial products might be one of the mechanisms of HCC progression.

The data in the present study showed that gut microbial richness decreased significantly in the terminal group compared with the early group. The decrease in richness was mainly due to the reduction in the abundance of SCFA-producing bacteria. The result was consistent with that of a previous study. The study showed that α-diversity significantly reduced in patients with BCLC stage C and D HCC compared with patients in stage A HCC [5]. Another similar study showed that the abundance of SCFA-producing family members declined significantly in patients with HCC compared with patients with hepatitis B [4]. SCFA-producing bacteria are known to participate in the process of fermenting diet fibers to SCFAs [22]. SCFAs are a major energy source of the intestinal enterocytes and are essential for maintaining the tight junction and intestinal barrier integrity [23]. However, in this study, the abundance of SCFA-producing bacteria did not correlate with the expression of the serum markers of gut damage and microbial translocation. It was speculated that multiple mechanisms might be involved in SCFA production and did not necessarily depend only on the change in the abundance of SCFA-producing bacteria.

Besides the changes in gut microbiota mentioned earlier, another predominant change was that the abundance of Bifidobacteriaceae decreased along with HCC progression accompanied by an increase in the abundance of Enterococcaceae. Bifidobacteria can reinforce gut barrier function, reduce mucosa inflammation, and protect the liver from injury [24]. Administration of Bifidobacteria could mitigate diethylnitrosamine- (DEN-) induced hepatocarcinogenesis in mouse models [25]. Enterococcus is regarded as pathogenic bacteria. The overgrowth of Enterococcus leads to the release of large amounts of polysaccharide A and LPS, which, in turn, increases gut permeability and facilitates microbial translocation. In this study, the ratio of Enterococcaceae to Bifidobacteriaceae correlated significantly with the expression of the serum markers of gut damage, microbial translocation, and HCC progression. Therefore, the ratio was considered an ideal quantitative marker of gut dysbiosis in patients with HCC. The ratio of abundance of Firmicutes to Bacteroidetes [26] or ratio of abundance of genus Bifidobacterium to the family Enterobacteriaceae was associated with the progression of chronic liver diseases [27]. However, both the ratios were similar among the three groups in this study. More studies should be performed in the future to confirm these results.

The changes in gut damage and microbial translocation biomarkers were subsequently detected due to the close relationship between gut dysbiosis and microbial translocation. REG3α is a C-type lectin antimicrobial peptide secreted into the gut lumen by Paneth cells [28]. REG3α invades into the bloodstream when gut barrier integrity is damaged [29]. Hence, the circulating level of REG3α is well accepted as a marker of gut damage. Previous studies reported that increasing REG3α level was always observed in patients with enteropathies, such as Crohn’s and celiac diseases, ulcerative colitis, and graft-versus-host disease [2931]. The results of this study suggested that gut damage worsened gradually along with HCC progression and facilitated the transfer of more microbial products into the bloodstream [32, 33]. This study was the first to show the relationship between gut damage and HCC progression.

This study also found that bacterial translocation increased with HCC progression. LBP is a well-known marker of bacterial translocation [34]. sCD14 is a marker of Gram-negative bacterial translocation, which is secreted by activated myeloid cells after LPS stimulation [35]. PGRP is a peptidoglycan recognition protein and a marker of Gram-positive bacterial translocation [10]. The present study results showed an elevation in the levels of the aforementioned markers with HCC progression. Moreover, increased LBP and sCD14 levels correlated with IL-6 levels, while increased sCD14 and PGRP levels correlated with IL-8 levels. In theory, elevated sCD14 and LBP levels indicated increased circulation of LPS. LPS resulted in the release of proinflammatory cytokines, such as IL-6 and IL-8, by T cells. Chronic inflammation contributed to liver damage and HCC progression [36, 37]. Therefore, it is reasonable that bacterial translocation correlated with HCC progression. However, no relationship was found between gut permeability and microbial translocation. This might be due to bacterial translocation via different mechanisms other than gut integrity. On the contrary, the damage of gut barrier alone might also not be sufficient to cause bacterial translocation, unless additional mechanisms in the complex interactions between host and microbes in the gut failed. A similar result was also observed in patients with alcohol use disorder. In this study, no correlation with intestinal permeability marker was observed, although the levels of all the three markers elevated significantly [10].

Another finding of this study was the elevation of inflammatory response profile during HCC progression, represented by IL-6 and IL-8. IL-8 is a key factor in neutrophil recruitment and activation, which also promoted angiogenesis and metastasis [38, 39]. One previous study found that high serum IL-8 was one of the meaningful predictive cytokines in patients with HBV- and HCV-related HCC [40]. Moreover, as a T cell-attracting chemokine, MCP-1 level decreased significantly in the terminal group, indicating a reduction in T cell immune response. The levels of other cytokines, such as IFN-γ, the signature cytokine of the Th1 cell; IL-4 and IL-13, the signature cytokines of the Th2 cell; and IL-17, the signature cytokine of the Th17 cell, reduced simultaneously but without significant differences. These results reflected comprehensive immunity suppression in terminal HCC. The results were in accordance with an ex vivo study. The study showed that bacterial extracts from patients with HCC microbiota elicited a T cell-immunosuppressive phenotype [41]. It is well known that SCFAs could promote and regulate the differentiation and apoptosis of Treg cells [42]. According to the results of a previous study and this study, a distinctive microbiome, characterized by decreased abundance of SCFA-producing bacteria, could modulate the peripheral immune response and result in HCC progression.

One limitation of this study was that no causation relationship could be deduced. Whether the alterations of gut microbial profiles, gut permeability, or bacterial translocation was the reason or the result of HCC progression needs further validation in animal models. The other limitation of this study was the lack of geographical representation of the study cohort. In a big hospital in Beijing, most patients came from North China. The third limitation was the lack of age-matched healthy controls. However, this study was designed to evaluate and compare microbiome dysbiosis in different stages of HCC and their relationship with gut damage and microbial translocation. In consideration of the effect of sex on gut microbiota, only male patients were included. The cause of HCC was only HBV and HCV, so the results need further validation in terms of other causes. Finally, the T cell function should be detected directly to confirm the results in this study.

5. Conclusions

This study indicated that the dysbiosis of gut microbiota, together with increased gut permeability and microbial translocation, was associated with elevated circulating inflammatory response and reduction in T cell response during HCC progression. The findings obtained in this study also suggested that gut microbiota might be a potential target for HCC treatment and intervention, especially in terminal HCC.

Abbreviations

HBV:Hepatitis B virus
HCC:Hepatocellular carcinoma
HCV:Hepatitis C virus
LBP:Lipopolysaccharide-binding protein
LPS:Lipopolysaccharide
MT:Microbial translocation
OTUs:Operational taxonomy units
PGRPs:Peptidoglycan recognition proteins
REG3α:Regenerating islet-derived protein 3α
sCD14:Soluble CD14
SCFAs:Short-chain fatty acids.

Data Availability

The data of clinical characteristics, relative abundance of gut microbiota, and plasma levels of REG3α, sCD14, LBP, PGRPs, cytokines, and chemokines of individuals used to support the findings of this study are included within the Supplementary Tables.

Additional Points

Lay Summary/Key Points. This study revealed that gut microbiota compositional and functional shifts correlated with hepatocellular carcinoma (HCC) progression, especially in the terminal stage. The mechanism might include elevated gut permeability and microbial translocation relevant to dysbiosis. The accompanying systemic inflammatory increase and immune suppression appeared to promote HCC progression. An increased ratio of Enterococcaceae to Bifidobacteriaceae could serve as a biomarker of dysbiosis in HCC progression.

Disclosure

The funding bodies had no role in the study design, data collection, data analysis or interpretation, or writing of the report.

Conflicts of Interest

All authors report no conflicts of interest relevant to this study.

Authors’ Contributions

Q.H., N.Z., and J.Z. conceptualized and designed the study. Y.G., N.Z., H.Y., S.L., N.C., Y.L., and J.Z. performed sample collection and data acquisition. N.Z., Y.G., J.Z., and Q.H. analyzed and interpreted the data. N.Z., Y.G., Q.H., and J.Z. prepared the original draft. Q.H. and J.Z. wrote and reviewed the manuscript. All authors read and agreed to the published version of the manuscript. Nan Zhang and Yusong Gou contributed equally to this study.

Acknowledgments

The authors thank all patients who agreed to participate in this study. This study was partly supported by the open funding of the Key Laboratory of Translational Research on Hepatitis B and Hepatocellular Carcinoma of Beijing Youan Hospital (BJYAHKF2018002); the fund of Clinical Research of Beijing Youan Hospital, CMU (YNKTTS20180108); and the fund of National Microbial Resource Center (No. NMRC-2021-2) and National Pathogen Resource Center (NPRC-32).

Supplementary Materials

Supplementary 1. Supplementary Table 1: detailed information of study subjects.

Supplementary 2. Supplementary Table 2: concentrations of 20 plasma cytokines and chemokines in the early, intermediated, and terminal groups of HCC patients.