Microbiota analysis with next-generation 16S rDNA gene sequencing in recurrent common bile duct stones.

Background: Endoscopic retrograde cholangiopancreatography (ERCP) is the main remedy for gallstones, but the postoperative recurrence rate is high. Recent research has indicated that the biliary microbiome takes part in the pathogenesis of cholelithiasis. However, it is not yet known whether biliary microbiome dysbiosis is relevant to recurrent cholelithiasis.
Methods: Thus, we investigated the bacterial communities of the biliary microbiomes of patients with recurrent common bile duct (CBD) stones and analyzed the relationship between recurrent CBD stones and biliary microbiota. The bile specimens of 5 patients with recurrent CBD stones (FF) and 45 patients with primary CBD stones (YF) were collected during the ERCP process. The microbiota was analyzed using 16S ribosomal DNA (rDNA) high-throughput sequencing. We also identified the link between recurrent CBD stones and biliary microbiota.
Results: Our results showed that at the phylum level, proteobacteria and firmicutes were the main two genera groups, and proteobacteria was high in FF patients. Additionally, synergistetes were high, but Bacteroidetes and actinobacteria were low in FF patients. The microbiomes in the bile of the YF patients were more evenly distributed than those in the bile of the FF patients. We also discovered that FF patients had decreased microbial bile diversity. At the genus level, klebsiella dominated in the FF patients, while Escherichia-shigella dominated in the YF patients. Additionally, klebsiella was higher in the FF patients than the YF patients. Conclusions: The observed differences in the genera between the recurrent CBD stone FF patients and the YF patients provide novel insights into the link between biliary microbiota changes and recurrent CBD stones. 2022 Annals of Translational Medicine. All rights reserved.

Introduction

Endoscopic retrograde cholangiopancreatography (ERCP) is widely used to remove bile duct stones. Many authors consider choledocholithiasis a late complication of ERCP (1-3). The recurrence rates of common bile duct (CBD) stones vary across different studies, but range from 4–24% (4). Research has confirmed that the reasons for gallstone formation are aberrant metabolism and the secretion of cholesterol and bile acids (5). Currently, a sphincter of Oddi laxity (SOL) is considered an important factor in cholangiolithiasis occurrence. Most patients with recurrent cholangiolithiasis have a history of biliary surgery, such as ERCP and the endoscopic stone extraction technique. After surgery, enterobiliary reflux caused by the loss of the sphincter function of Oddi leads to intestinal flora ectopia (6-8). Consequently, the bile duct microenvironment is thought to be another important factor in cholangiolithiasis occurrence.

The enhanced reflux of intestinal contents changes the microenvironment and promotes the process of gallstone formation. The number of molecular studies focusing on the association between the biliary microbiome and gallstones has continued to increase in the past few years. For example, a study by Liang et al. revealed that patients with SOL had a more severe bacterially infected bile duct microenvironment and stronger lithogenicity (9). Research continues to confirm the role of the biliary microbiome in the formation of gallstones. For example, urease positive helicobacter species and enterohepatic helicobacters have been shown to promote the process of gallstone formation (10,11). There is a close relationship between the formation of gallstones and microorganisms. For example, β-glucuronidase secreted by microorganisms contributes to the formation of calcium bilirubinate stones (12). Phospholipase and mucin also play an important role in the formation of gallstones (13). However, it is not yet known whether biliary microbiome dysbiosis is connected to recurrent cholelithiasis.

Next-generation sequencing (NGS) is a high-throughput method that allows rapid sequencing of base pairs in DNA or RNA samples. Supporting a variety of applications, including gene expression profiling, chromosome counting, detection of epigenetic changes, and molecular analysis, NGS facilitates discovery and enables the future of personalized medicine (14). In our study, by using NGS, we conducted 16S ribosomal DNA (rDNA) sequencing to study the structure and composition of bile microbial communities in primary CBD stone (YF) and recurrent CBD stone (FF) patients and then explored differences in the bacterial communities in bile between the YF and FF groups. We also sought to study the characteristics of the biliary microbiome in FF patients and its potential connection to recurrent CBD. According to the results, we hypothesized that there was an underlying association between the biliary microbiome and recurrent CBD stones. This study revealed the differences of microbial structure in bile between primary and recurrent bile duct stones for the first time, and this study provided clues for the early diagnosis of recurrent bile duct stones. We present the following article in accordance with the MDAR reporting checklist (available at https://atm.amegroups.com/article/view/10.21037/atm-22-2247/rc).

Methods

Sample collection

Data were collected from 5 recurrent CBD stone patients and 45 primary CBD stone patients who had been admitted to the Guangdong Second Provincial General Hospital from October 2019 to March 2020. None of the patients had received antibiotic treatment at least 3 months before the ERCP procedure. Written informed consent was obtained from all the patients, and the study conformed to the Declaration of Helsinki (as revised in 2013). The study was also approved by the Ethics Committee of the Guangdong Second Provincial General Hospital (Approval No. 2021-ET-15). All the samples were collected by ERCP.

Experimental procedures

The collected samples were sequenced for 16S rDNA by a basic process that comprised DNA extraction, DNA detection, polymerase chain reaction (PCR) amplification, purification, library preparation and quality testing, and sequencing (see Figure S1). These processes are described below in detail.

DNA isolation

The total DNA in all the bile samples was extracted using the E. Z. N. A.® Stool DNA Kit (D4015, Omega, Inc., USA) and stored at −80 ℃. Agarose gel electrophoresis was applied to detect the DNA quality, and an ultraviolet (UV) spectrophotometer was used to research the DNA quantity (15).

PCR amplification

The hypervariable V3–V4 regions of the bacteria 16S rDNA in the genomic DNA were amplified with the following primers: 341F (5'-CCTACGGGNGGCWGCAG-3') and 805R (5'-GACTACHVGGGTATCTAATCC-3'). PCR amplification was performed to amplify the bacteria 16S fragments.

16S rDNA sequencing

In the PCR process, ultrapure water was chosen as a negative control (NC). AMPure XT beads (Beckman Coulter Genomics, Danvers, MA, USA) and Qubit (Invitrogen, USA) were applied to purify and quantify the PCR products. Next, the amplicon pools were applied for sequencing. The size and quantity of the amplicon libraries were determined using an Agilent 2100 Bioanalyzer (Agilent, USA) and Illumina library quantitative kits (Kapa Biosciences, Woburn, MA, USA). Finally, the libraries were also sequenced on the 250 PE MiSeq platform.

Bioinformatics

The 16S rDNA amplicons were sequenced on the Illumina MiSeq platform (LC-Bio). After quality filtering, high-quality clean tags were obtained using FQTRIM (v0.94). Next, the operational taxonomic units (OTUs) were clustered using the Ribosomal Database Project classifier. Next, the alpha diversity analysis was calculated with QIIME (version 18.0), the beta diversity analysis was calculated by a principal coordinate analysis (PCoA), and a cluster analysis was used to analyze the complexity of the diversity of the species. Additionally, a redundancy analysis (RDA) was also conducted.

Statistical analysis

The independent duplicate data are represented as the mean ± standard deviation (SD), and were examined using SPSS 20.0 (SPSS, Chicago, IL, USA) with the Student’s t-test or Mann-Whitney U-test. A P value <0.05 was considered statistically significant.

Results

Diversity analysis of the microbiota in the bile of the YF and FF patients

Through 16S rDNA sequencing, we further analyzed the microbial composition difference of 45 samples from the YF patients and 5 samples from the FF patients. First, a phylogenetic analysis was performed by tags obtained from the YF and FF patients. The tags were clustered into OTUs with a 97% threshold after filtering out the chimeras and singletons. A Venn diagram of the OTUs was generated at a level of 97% similarity to compare the microbiota compositions of the YF and FF groups. From the Venn diagram of the two groups, we observed that 480 of the 4,141 OTUs were shared by the two groups, while 3,432 were specific to the YF group, and 229 were specific to the FF group (see ).

Figure 1

Diversity analysis of microbiota in the bile of the YF and FF patients. (A) Venn diagram representing the unique and shared OTUs in the YF and FF patient groups. (B) Shannon-Wiener curves of the microbiota in the TF and FF groups. (C) The abundance and diversity were determined by a UPGMA hierarchical clustering analysis (unweighted unifrac and weighted unifrac). (D) Diversity index differences between the YF and FF groups were counted and exhibited in a representative box plot (unweighted and weighted). (E) PCoA. Each dot represents 1 sample. The green colors represent the YF group, and the red colors represent the FF group. FF: recurrent CBD stones patients; YF: primary CBD stones patients. PCA, principal component analysis; PCoA, principal coordinate analysis; OTUs, operational taxonomic units; CBD, common bile duct.

The Shannon index was adopted to signify the change in biliary microbiota diversity and richness. The saturation plateaus of the rarefaction curves were obtained and the Shannon diversity of all the sequences was stable. Consequently, we suggested that most of the diversity had already been discovered (see ).

We also used hierarchical clustering to visualize relationships among the TF and FF samples according to the average linkage. The relationships among the samples based on the similarity of the microbial composition, which was calculated by the UPGMA, are presented in a dendrogram (see ). Additionally, a box plot shows the differences between the two groups using a METAGE Nassist pipeline, and the data revealed that the distance was significantly lower in the FF group than the YF group (see ). Conversely, the PCoA of the weighted UniFrac distances was used to compare the bacterial communities between samples from the FF group (red circles; n=5) and YF group (blue circles; n=45). As shows, while the bile of the YF and FF communities were similar for a large portion of intersection, which indicates that the YF and FF samples shared bacterial communities, there was a clear separation between the two groups (see ). Based on the above analyses, we confirmed that microbial signatures were available for the patient groups and could be used to distinguish among the samples.

Species analysis of the different microbiota in the bile of the YF and FF patients

Subsequently, the different microorganisms were analyzed to examine the species. Based on the feature annotation results and feature table, a species classification abundance table was obtained. Additionally, based on the species abundance and annotation table, we selected the 30 species with the highest abundance and presented the relative abundance of each sample in different forms. Our research focused on classification at the phylum and genus levels. As shows, at the phylum level, 1,276 OTUs belonged to 30 phyla, while at the phylum level, proteobacteria and firmicutes were the main two genera in both the YF and FF groups, which is consistent with previous research results (16,17). Compared to the YF patients, the FF patients had notably higher synergistetes and dramatically lower Bacteroidetes and actinobacteria. Additionally, the microbiomes in the bile of the YF patients were more evenly distributed than those in the bile of the FF patients. The heat-map analysis led us to the same conclusion; that is, that certain bacteria, such as Bacteroidetes and actinobacteria, were more abundant in the YF than the FF patients; however, synergistetes were more abundant in the FF patients than the YF patients (see ).

Figure 2

Species analysis of different microbiota in the bile of the YF and FF patients. (A) The relative abundance and distribution of different bacteria are displayed using stacked bar charts at the phylum and genus levels in each sample. (B) The relative abundance of different bacteria is exhibited in the heat-map analysis at the phylum and genus levels. (C) The cladogram of the microbiota is displayed. (D) The bar-plot difference analysis shows the relative abundance of genes through 16S rDNA sequencing in the YF and FF groups. FF: recurrent CBD stones patients; YF: primary CBD stones patients. rDNA, ribosomal DNA; CBD, common bile duct.

Additionally, to study the phylogenetic relationships between the microbiota, we used a phylogenetic tree to display the evolutionary relationship. Different branches represent different levels of genus classification. Different genera of the same color indicate that branches belong to the same phylum. The genetic similarity of branches is negatively correlated with their distance in the clade (see ). Additionally, a bar plot displays the differences of the microorganisms at each level, and indicates whether there was a significant difference between the groups of the species (P<0.05). As Figure 2D shows, the bacteria that exhibited a relatively high abundance included klebsiella, rhodococcus, prevotella, alloprevotella, acinetobacter, stenotrophomonas, curvibacter, novosphingobium, and raoultella.

Microbiota alterations in the bile of the YF and FF patients

We also performed an advanced analysis of microbiota alterations in the bile of the YF and FF patients. First, the bubble chart represents the relative abundance of the 16S rDNA gene (%) in the YF and FF patients at the phylum and genus levels. The size of the bubble relates to the relative abundance of the 16S rDNA gene (%) at the genus level. The color of each bubble represents different phylum. According to the results, we found that the biliary microbiome is mainly facultative anaerobe proteobacteria, of which the genus klebsiella and Escherichia-shigella are the most two abundant. Additionally, klebsiella was distinctly higher in the FF patients than the YF patients (see ).

Figure 3

Microbiota alterations in the bile of the YF and FF patients. (A) A bubble plot shows the species annotation and abundance of two-dimensional analysis. (B) We analyzed the distribution proportion of each dominant species in different groups using Circos diagrams at the phylum and genus levels. (C) The potential biomarkers were defined by LEfSe. (D) Sankey plots of different microbiota. The relative abundance of different microbiota is shown at the phylum level (middle) and genus level (right) for the different samples (left). FF: recurrent CBD stones patients; YF, primary CBD stones patients. LEfSe, linear discriminant analysis effect size; CBD, common bile duct.

Second, the distribution of the top 5 most abundant bacteria is shown in a Circos diagram. The left (large) semicircle denotes the species composition (the inner bands signify the different taxa, and the outer bands signify the species), and the right (small) semicircle denotes the taxonomic distribution ratio of the species in different samples (the inner bands signify the species, and the outer bands signify the different taxa), and the length denotes the distribution of the samples. Compared to the YF patients, the FF patients had observably higher proteobacteria, and prominently lower Bacteroidetes and actinobacteria (see ).

Third, the linear discriminant analysis effect size (LEfSe) analysis revealed that several differentially abundant taxa were significantly mor present in specific fractions. The differentially abundant taxa were highlighted in a phylogenetic tree. The linear discriminant analysis scores were used to test the differentiation between the YF and FF patients. The colored nodes from the inner to the outer circles denote the taxa from the phylum to the genus levels. The evidently different taxa are represented by different colors. The taxa were defined by LEfSe. The cladogram shows the potential biomarkers of the YF and FF groups (see ). Next, Sankey plots show the complex pairing interactions between the phylum level (middle) and the genus level (right) in FF and YF. According to the results, the biliary microbiome in the YF and FF groups are mainly facultative anaerobe proteobacteria, of which the genus klebsiella and Escherichia-shigella are the top 2 most abundant. Additionally, klebsiella was significantly higher in the FF patients than the YF patients (see ).

Distribution of a single species in the bile of the YF and FF patients

Additionally, we further predicted the differential flora at the phylum level and the microbial-associated phenotypes. The Manhattan plots show the abundance of the microbiota enriched in the FF bile relative to the YF bile. The most different bacterial community was found in the proteobacteria (see ). Next, scatter plots were used to predict the different phenotypes in the YF and FF groups, and the results showed that compared to the YF group, the relative abundance of the aerobic was reduced, while the relative abundances of the stress tolerance, content of mobile elements, formation of biofilms, potentially pathogenic, and facultatively anaerobic were elevated in FF the group; however, the relative abundance of anaerobic, gram negative and gram positive did not change between the two groups (see ).

Figure 4

Distribution of a single species in the bile of the YF and FF patients. (A) Manhattan plots showing the abundance of microbiota in the bile of the YF and FF patients. (B) Bacterial phenotypes as predicted by BugBase. FF: recurrent CBD stones patients; YF: primary CBD stones patients. CBD, common bile duct.

Discussion

The bile duct microenvironment and SOL are two important factors in cholangiolithiasis occurrence, as the enhanced reflux of intestinal contents change the microenvironment and promote the process of gallstone formation. After ERCP, enterobiliary regurgitation due to a constriction disorder of the Oddi sphincter lead to intestinal flora ectopia, which provides the conditions for stone recurrence. Recently, culture-independent 454 pyrosequencing of the bacterial 16S rDNA gene has been used to study components of the microbiome, and this approach has greatly enhanced understandings of how microbiome are associated with human diseases (18). Many studies have focused on the gut microbiota. Research has confirmed that gut microbiota are linked with multiple diseases, such as type 2 diabetes, non-alcoholic fatty liver disease, Crohn’s disease, obesity, and cirrhosis with minimal hepatic encephalopathy (19-23). Additionally, gut microbiota contribute to the formation of gallstones (16,24). For example, gut microbiota in the lithogenic bile indicate an increase in intestinal permeability during biliary obstruction (25,26). Gut microbiota have also been found to be associated with elevated inflammation and stone formation (27). Further, previous research has identified the microbiota of the CBD (28-30). A few numbers of studies have addressed alterations in the microbial makeup of recurrent gallstones. Ye et al. showed that the invasion of E. coli during ERCP may be one of the potential factors for gallstone recurrence (4). Escherichia coli was also detected in the samples involved in this study, but whether it is related to the recurrence of gallstones needs further exploration.

Keizman et al. found that the recurrence rate of CBD stones in elderly patients is about 20%, and noted that the incidence of periampullary diverticulum is increased in elderly patients (31). Periampullary diverticulum, which increases biliary pressure by compressing the peripheral bile duct, is thought to be a risk factor for the recurrence of CBD stones (32). Based on the results, we hypothesized that there is a close relationship between the damaged homeostasis of biliary flora environment after biliary tract surgery and the recurrence of CBD stones.

We aimed to study the bacterial communities of the biliary microbiome in FF patients and investigate the relationship between recurrent CBD stones and biliary microbiota. A better understanding of differences in the biliary microbiome of recurrent and primary patients is key to extending understandings of the relationship between biliary microbiome and cholangiolithiasis recurrence. A 16S rDNA gene-based sequencing analysis was conducted to compare the composition and activity of the bile microbiota of an FF group to a YF group. We also conducted a heat-map analysis, a UniFrac-based PCoA, and an analysis of variance to identify the differences in microbial composition and diversity between the two groups. We discovered significant variations in the biliary microbial components between the YF and FF patients. The 4,141 OTUs belonged to 30 phyla. Additionally, alpha and beta diversity analyses were conducted to compare microbial diversity and species complexity in the bile samples of the two groups. The saturation plateaus of the rarefaction curves were obtained, and the Shannon diversity of all the sequences was stable, indicating that most of the diversity in the bile of the patients had already been discovered. The results of the beta diversity analysis confirmed that microbial signatures can be used to group patients and to distinguish among samples. At the phylum level, we found that proteobacteria and firmicutes were the main two genera groups. Proteobacteria was significantly higher in the FF group than the YF group. Eckburg et al. found that proteobacteria is a minor component in normal intestinal microbial communities (33). However, facultative anaerobic proteobacteria was shown to be significantly increased because of the destruction of intestinal microbiota after antibiotic treatment (34). Shin et al. found that the continuous increase of proteobacteria abundance is a signal of ecological imbalance and disease (35). Consequently, we speculated that patients with recurrent choledocholithiasis have a significant imbalance of biliary flora.

Additionally, relative to the YF patients, the FF patients had observably higher synergistetes, and significantly lower Bacteroidetes and actinobacteria. However, the microbiomes in the bile of the YF patients were more evenly distributed than those in the bile of the FF patients. We also observed a decrease in the microbial diversity of the bile of the FF patients compared to that of patients with cholecystolithiasis for the first time.

A reduction in biodiversity destroys the resilience of a natural ecosystem and may lead to the serious degradation of that ecosystem (36). Thus, we hypothesized that the reduction of biodiversity caused by the destruction of ecological stability provides the necessary conditions for the recurrence of choledocholithiasis. According to the results, we found that the biliary microbiomes were mainly facultative anaerobe proteobacteria, of which the genus klebsiella and Escherichia-shigella were the top 2 most abundant in the two groups. Additionally, klebsiella was significantly higher in the FF patients than the YF patients (see ). At the genus level, klebsiella dominated in the FF group, while Escherichia-shigella dominated in the YF group. Escherichia-shigella and klebsiella are both opportunistic pathogens, and their growth increases the risk of host disease. A previous study has found that Escherichia coli is the main bacteria in the bile samples of patients with gallstones (25). Liu et al. suggested that Escherichia coli is prominently correlated with bile endotoxin and plays a key role in the progress of acute cholecystitis (37). Razaghi et al. found that choledocholithiasis is associated with Escherichia coli colonization in the biliary tract (38). Liang et al. found that aerobe and facultative anaerobe are the most common bacteria in the bile ducts of choledocholithiasis patients. Among them, proteobacteria and firmicutes were the most common phylum, and Enterobacteriaceae was the most abundant family (9). These findings are basically consistent with those of our study. However, few studies on the biliary microbiota of recurrent gallstones appear to have been conducted. Klebsiella together with Escherichia coli, the yersinia species, salmonella species, and shigella species are the main members of the Enterobacteriaceae family. For a long time, klebsiella was considered an opportunistic pathogen, as it was usually found in immunocompromised individuals and hospitalized patients (39).

Gorrie et al. and Martin et al. reported that the rate of gastrointestinal colonization in hospitalized patients was around 20% (40,41). Research has shown that the intestinal colonizing klebsiella invades the intestinal mucous membrane and then enters the portal vein, and eventually leads to liver abscesses (42). Keizman et al. found that the recurrence rate of CBD stones in elderly patients is about 20% (31). We know that elderly patients have low levels of immunity. According to the present study, klebsiella was significantly higher in the FF patients than the YF patients. We hypothesized that the intestinal colonizing of klebsiella increases after biliary tract surgery and plays a key role in cholangiolithiasis recurrence.

In conclusion, we demonstrated that patients with recurrent cholelithiasis have a potential imbalance of biliary microbiota, and that biliary microbiota may contribute to the formation of biliary calculi in patients with recurrent cholelithiasis. Thus, our findings have extended understandings of the biliary microbiota of FF patients and may be helpful in the prevention of cholangiolithiasis recurrence. In this study, klebsiella was more abundant in the FF than in the YF. In view of the reported association between klebsiella and gallstones (25), we suggest that the detection of klebsiella may be a potentially effective method in monitoring the recurrence of gallstones.

The article’s supplementary files as