Insight into the mechanism of gallstone disease by proteomic and metaproteomic characterization of human bile

Introduction: Cholesterol gallstone disease is a prevalent condition that has a significant economic impact. However, the role of the bile microbiome in its development and the host’s responses to it remain poorly understood.

Methods: In this study, we conducted a comprehensive analysis of microbial and human bile proteins in 40 individuals with either gallstone disease or gallbladder polyps. We employed a combined proteomic and metaproteomic approach, as well as meta-taxonomic analysis, functional pathway enrichment, and Western blot analyses.

Results: Our metaproteomic analysis, utilizing the lowest common ancestor algorithm, identified 158 microbial taxa in the bile samples. We discovered microbial taxa that may contribute to gallstone formation, including β-glucuronidase-producing bacteria such as Streptococcus, Staphylococcus, and Clostridium, as well as those involved in biofilm formation like Helicobacter, Cyanobacteria, Pseudomonas, Escherichia coli, and Clostridium. Furthermore, we identified 2,749 human proteins and 87 microbial proteins with a protein false discovery rate (FDR) of 1% and at least 2 distinct peptides. Among these proteins, we found microbial proteins crucial to biofilm formation, such as QDR3, ompA, ndk, pstS, nanA, pfIB, and dnaK. Notably, QDR3 showed a gradual upregulation from chronic to acute cholesterol gallstone disease when compared to polyp samples. Additionally, we discovered other microbial proteins that enhance bacterial virulence and gallstone formation by counteracting host oxidative stress, including sodB, katG, rbr, htrA, and ahpC. We also identified microbial proteins like lepA, rtxA, pckA, tuf, and tpiA that are linked to bacterial virulence and potential gallstone formation, with lepA being upregulated in gallstone bile compared to polyp bile. Furthermore, our analysis of the host proteome in gallstone bile revealed enhanced inflammatory molecular profiles, including innate immune molecules against microbial infections. Gallstone bile exhibited overrepresented pathways related to blood coagulation, folate metabolism, and the IL-17 pathway. However, we observed suppressed metabolic activities, particularly catabolic metabolism and transport activities, in gallstone bile compared to polyp bile. Notably, acute cholelithiasis bile demonstrated significantly impaired metabolic activities compared to chronic cholelithiasis bile.

Conclusion: Our study provides a comprehensive metaproteomic analysis of bile samples related to gallstone disease, offering new insights into the microbiome-host interaction and gallstone formation mechanism.

1 Introduction

Cholesterol gallstone disease is a prevalent condition and a significant economic burden worldwide. The prevalence rates of gallstone disease vary across regions, with approximately 30–40% in South America, around 20% in Europe, and 5–8% in Asia (Chen et al., 2006; Lammert et al., 2016; Li and Gao, 2019). Gallstones can occur in various parts of the human biliary system, including the gallbladder, extrahepatic bile duct, or intrahepatic duct. Gallstones are commonly classified as either cholesterol or bilirubin gallstones, with cholesterol gallstones being the most predominant type. While the exact cause of gallstone formation remains unclear, various factors contribute to the development of gallstone disease, including genetic predisposition, cholesterol crystallization, and bile supersaturation. Additionaly, several risk factors have also been identified, such as fatty liver, diabetes, obesity, and dysbiosis of the microbiome (Di Ciaula et al., 2018). Notably, microbial infections and the interaction between the host and the microbiome play a significant role in the development of cholecystitis (CHL) and gallstone disease (Grigor'eva and Romanova, 2020).

Dysbiosis, characterized by altered bacterial composition and abundance, has been implicated in the pathogenesis of various diseases, including inflammatory conditions, cancers, cardiovascular diseases, and metabolic disorders like diabetes and obesity (Thursby and Juge, 2017). To investigate the changes in the microbiota and its association with gallstone formation, researchers have utilized 16S rRNA gene sequencing to analyze microbial components in bile, gut, and gallstones of patients with cholelithiasis (Shen et al., 2015; Kose et al., 2018; Molinero et al., 2019). The microbiota in the gut, biliary tract, and gallbladder may contribute to gallstone development. For instance, studies have demonstrated that gut microbes, such as Desulfovibrionales, can influence the formation of cholesterol-type gallstones by regulating bile acid production and cholesterol secretion (Hu et al., 2022). Furthermore, dysbiosis in the biliary microbiota has been associated with the recurrence of bile duct stones (Tan et al., 2022). Specific components of the bile microbiome, such as Enterobacteriaceae, have also been linked to the development of CHL and gallbladder cancer (Choe et al., 2021). Therefore, 16S rRNA gene sequencing serves as a valuable tool for identifying and assessing the relative abundance of microbiota in bile samples, including at the species level, in both diseased individuals and healthy controls. It is important to acknowledge that DNA sequencing technologies alone cannot directly elucidate the protein function of microbial communities within these bile samples.

Mass spectrometry (MS)-based proteomics has emerged as a powerful tool for investigating protein expression patterns in human samples. Several proteogenomic studies have been performed to analyze the composition of human protein in bile samples obtained various clinical conditions, such as sclerosing cholangitis (Reinhard et al., 2012; Rupp et al., 2018), biliary stenosis (Farina et al., 2009), cholangiocarcinoma (Kristiansen et al., 2004; Farid et al., 2011; Lankisch et al., 2011; Shen et al., 2012; Navaneethan et al., 2015; Laohaviroj et al., 2017; Voigtlander et al., 2017; Ren et al., 2019; Son et al., 2020), gangrenous CHL (without gallstones) (Wang L. et al., 2018) and gallstone disease (Zhou et al., 2005; Barbhuiya et al., 2011; Zhang et al., 2013). However, these studies were constrained by the utilization of pooled-sample strategies, low throughput methods, or outdated platforms and technologies, which limited the scope of protein identification. Recently, a study analyzed the bile metaproteome of gallbladder cancer, gallstones, and healthy controls (Sharma et al., 2022). Although the study successfully explored microbial taxa based on microbial peptides, further investigation into microbial proteins remains incomplete.

In the present population-scale study, we conducted a comprehensive proteomic and metaproteomic analysis of gallbladder bile samples obtained from 40 patients diagnosed with gallstone disease and gallbladder polyp. We employed trapped ion mobility spectrometry coupled with time-of-flight mass spectrometry (TIMS-TOF-MS) to explore the microbial composition and microbiome characteristics linked to gallstone disease, as well as the host immune responses that could contribute to gallstone formation and cholelithiasis.

2 Materials and methods

2.1 Human bile sample collection

A total of 40 human bile samples were analyzed, consisting of 31 gallstone bile samples and 9 control bile samples (gallbladder polyp or normal bile). Among the 31 gallstone patients, 14 had acute CHL, while 17 had chronic CHL (Supplementary Table S1). The 9 gallstone-free controls consisted of individuals with polyps but no visible gallstone or CHL symptoms. Detailed information regarding the samples can be found in Supplementary Table S1. The bile samples from patients diagnosed with cholesterol gallstone disease were collected at the Department of General Surgery of Shanghai Tongren Hospital (Shanghai, China). Control bile samples were obtained from individuals with gallbladder polyps who underwent thorough examinations to confirm the absence of gallstones and gallbladder inflammation. Bile extraction from the gallbladder was performed during cholecystectomy using a sterile syringe, and subsequently, the samples were stored at −80°C until further analysis. Written consent was obtained from all patients or their relatives, and the study was approved by the Institutional Research Ethics Committee of the Institutes of Biomedical Sciences, Fudan University.

2.2 Experimental design and statistical rationale

The proteins from each sample were separated using SDS-PAGE, followed by in-gel tryptic digestion. The resulting digestion products were then subjected to analysis using reverse-phase liquid chromatography-mass spectrometry (RPLC–MS). A total of 40 LC–MS analyses were conducted.

2.3 Extraction of bile proteins

To extract proteins from the bile samples, 200 μL of raw bile was mixed with 400 μL of chloroform and 10 μL of acetic acid. The mixture was vortexed for 1 min and then chilled on ice. This cycle was repeated 10 times. Subsequently, the tube was centrifuged at 12,000 rpm for 10 min at 4°C, and the enriched lower organic phase containing bilirubin was discarded. The upper and middle layers containing proteins were collected and mixed with chloroform and acetic acid. The extraction and centrifugation steps were repeated as described above. The resulting supernatant was then added to 1.2 mL of ice-cold acetonitrile (ACN) and placed in a −20°C freezer to promote precipitation. After centrifugation at 12,000 rpm for 10 min at 4°C, the obtained pellet was suspended with 1 mL of ice-chilled acetone and sonicated to ensure complete suspension. Subsequently, the tube was centrifuged, and the pellet was suspended and washed with 1 mL of isopropanol. Another round of centrifugation was performed, followed by resuspension and washing with 1 mL of dehydrated alcohol. Finally, the resulting pellet was dissolved using 100 μL 0.5 M triethylammonium bicarbonate buffer (TEAB, cat no. T7408, Sigma Life Science, Shanghai, China) with 0.1% SDS, aided by sonication. Protein concentrations were determined using the BCA protein assay kit (Beyotime, LTD Inc., Shanghai, China) according to the manufactory’s instructions.

2.4 SDS-PAGE and in-gel tryptic digestion

The proteins of the mixed samples (30 μg) were separated using mini-size 10% SDS-PAGE. Each lane was cut into one slice (Figure 1) and subjected to in-gel tryptic digestion, following a previously described protocol with slight modifications (Cui et al., 2013). Gel slices were cut into approximately 1 mm³ granules using a surgical knife. A destaining solution containing 50% ACN was added to the tube, which was then agitated at 1,250 rpm for 5–10 min. Subsequently, the tube was centrifuged and the supernatant was discarded. This destaining step was repeated three times to ensure the removal of residual dye. The gel granules were shrunk by adding 1 mL of 100% ACN and agitated. The tube was then centrifuged, and the supernatant was aspirated. The granules were vacuum dried for 20 min. For the reduction reaction, 300 μL of freshly prepared 64 mM DTT/25 mM NH₄HCO₃ was added, followed by incubation at 56°C for 1 h with gentle agitation. The supernatant was aspirated, and the free sulfhydryl groups were blocked with 300 μL iodoacetamide (IAA) solutioin (130 mM in 25 mM NH₄HCO₃) for 45 min at room temperature in a dark place. The tubes were centrifuged to remove the supernatant, and the pellets were incubated in 500 μL of 25 mM NH₄HCO₃ for 3 min with agitation. The supernatant was removed by centrifugation, and the pellets were shrunk using 100% ACN. The granules were then swollen with the enzyme mixture containing 16 ng/μL trypsin (cat no. V5111, Promega, Beijing, China) and 3 ng/μL rLysC (HLS LYS001C, HUA LISHI SCIENTIFIC, Beijing, China). Excess enzyme solution was removed, and digestion was carried out overnight at 37°C with gentle agitation. After digestion, 300 μL of 50% ACN/5% formic acid (FA) solution was added to the tubes, which were incubated in a water bath with sonication for 10–20 min. Following centrifugation, the supernatant was collected. The extraction step was repeated four times, and the supernatants were combined and lyophilized.

FIGURE 1

Figure 1. Proteomic and metaproteomic analysis of human bile with or without gallstone. (A) The cohort enrolled in this study contains 14 acute CHLs, 17 chronic CHLs and 9 gallstone-free patients with gallbladder polyps. The bile fluids were collected during gallbladder resection of gallstone patients or polyp remove of patients. The proteins were extracted and separated using SDS-PAGE, followed by in-gel tryptic digestion. LC–MS analysis was performed using TIMS-TOF pro mass spectrometer. (B) The peptides were identified by two-step database searching using Metalab against a comprehensive protein database (UniProtDB) and the peptides were further subjected to taxonomic lowest common ancestor (LCA) analysis. Metataxonomic analysis using public 16S rRNA datasets of human bile by Kraken2 program was performed independently to verify the metaproteomic results. (C) The correlation plot of the bile samples based on protein identification. (D) PLS-DA of bile samples based on identified microbial taxa. The ellipses illustrated the confidence interval.

2.5 Desalting

The peptides were dissolved in 100 μL of 0.1% FA and desalted using MonoSpin C18 (Cat. No. 5010–21701, GL Sciences Inc.) according to the manufacturer’s instruction with minor modifications. In brief, the blank spin column was conditioned with 100 μL of 100% ACN/0.1% FA solution, followed by a conditioning step with 200 μL of water. Next, the 100 μL peptide sample was loaded onto the spin column for adsorption. After centrifugation to remove the flow-through, the spin column was washed with 100 μL of water. Subsequently, the clean peptides were eluted with 150 μL of 50% ACN/0.1% FA solution and lyophilized.

2.6 Peptide quantitation

Peptide quantification was conducted employing the Pierce Quantitative Colorimetric Peptide Assay kit (Pierce# 23275, Thermo Fisher Scientific, Shanghai, China), following the manufacturer’s guidelines. We used standard peptides to create a dilution series for generating a standard curve. Peptide samples underwent a 10-fold dilution with ultrapure water. The reaction process involved the addition of working reagents, followed by an incubation period of 15 min at 37°C. We measured the absorbance at a wavelength of 480 nm using a microplate reader. The concentration of peptides in each sample was calculated referencing the standard curve. Lastly, we adjusted the concentration of each peptide sample was adjusted to 200 ng/μL before proceeding with mass spectrometry analysis.

2.7 Label-free LC–MS analysis of bile proteome

Bile protein digests from each sample were subjected to a single-shot analysis using a TIMS-TOF Pro mass spectrometer (Bruker Daltonics GmbH, Bremen, Germany) equipped with a Nano electrospray ionization (ESI) source. We separated 200 ng of peptides using an Ionopticks Aurora Series emitter column (25 mm x 75 μm ID, 1.6 μm C18). The mobile phase A (H₂O/0.1% formic acid) and mobile phase B (100% acetonitrile/0.1% formic acid) were delivered through a liquid chromatograph nanoElute (Bruker Daltonics GmbH, Bremen, Germany) at a flow rate of 200 nL/min, maintaining the column temperature at 50°C. The 30-min gradient elution program was set as follows: an increment from 5% B to 22% B over 20 min, a rise to 37% B for 4 min, a subsequent increase to 80% B for 2 min, followed by an isocratic hold at 80% B for 4 min. The survey MS scans ranging from m/z 100 to 1,700 were succeeded by 10 Parallel Accumulation-Serial Fragmentation (PASEF) tandem MS scans with a cycle time of 1.16 s. The intensity threshold for online PASEF analysis was set at 5000, with the accumulation and ramp times fixed at 100 ms. The ion source voltage and temperature were 1,500 v and 180°C, respectively. The drying gas flow rate was set to 3 L/min, and the ionic mobility range was configured to 0.7–1.3 Vs/cm².

2.8 Metaproteomic database searching

The raw data generated by TIMS-TOF Pro were converted into mxXML files using ProteoWizard Version: 3.0.21352-04cc5bb¹ (Chambers et al., 2012). These mxXML files served as input for the downstream metaproteomic database searching, which was automated using Metalab version 2.2.1 (Cheng et al., 2017). Metalab incorporates X!Tandem (version 20151215) and MaxQuant (version 1.6.5.0) for the first and second round of database searching, respectively (Craig and Beavis, 2004; Tyanova et al., 2016). In the first round of X!Tandem searching, the database search was conducted against a custom microbial database (407,878 entries) comprising protein sequences from Bacteria, Virus, Archaea, and fungus. These sequences were extracted from the UniProtKB/Swiss-Prot database of Release 2022_01 (566,996 entries). To construct the custom database, we obtained all taxonomy IDs associated with Bacteria (taxon id 2), Virus (taxon id 10239), Archaea (taxon id 2157) and Fungus (taxon id 4751) using NCBI Entrez Edirect. The protein sequences were then extracted from the UniProtKB/Swiss-Prot database based on the organism annotation (“OX”) provided in each sequence title. The database search was performed in parallel against the Human (taxon id 9606) protein database, which comprised 20,377 sequences extracted from UniProtKB/Swiss-Prot database of Release 2022_01. The fixed post-translation modification (PTM) was carbamidomethylation of cysteine (C), while the variable modifications included acetylation of protein N-terminus, deamidation of asparagine (N) and glutamine (Q), and oxidation of methionine (M). The enzyme used was trypsin, allowing for two missed cleavage sites. No spectra clustering was performed. A sample-specific database was generated at the end of X!Tandem searching.

In the second round searching, MaxQuant was utilized to search against the reduced sample-specific database. The false discovery rates (FDR) for peptide-spectrum matches (PSMs), protein, and modification site was set at 0.01. The search was operated in parallel using revert decoy mode and included contaminant sequences. The match-between-run function was applied with the following parameters: Match ion mobility window, 0.05; Alignment time window, 20 min; Alignment ion mobility, 1; Match unidentified features, true. Protein expression was quantified using the label-free quantification (LFQ) module, with the following parameters: LFQ min. Ratio count, 1; Label min. Ration count, 1. Proteins that shared the same set of peptides were automatically grouped into a protein group. Default settings were accepted for all other parameters not specified here. The resulted proteinGroups.txt file was filtered using R scripts, and proteins labeled as “Potential contaminant,” “Reverse” or “Only identified by site” were discarded. Additionaly, proteins with fewer than 2 unique peptides were discarded.

Using the peptides identified by MaxQuant, Metalab performed taxonomic lowest common ancestor (LCA) analysis against the built-in taxonomic database. The analysis revealed the top nodes, including Eukaryota, Bacteria, Archaea, and Viruses. Taxa were determined based on a minimum of three unique peptide counts. The quantification of each taxonomic node was performed using the LFQ value.

2.9 Functional enrichment analysis

We obtained the Gene ontology (GO) annotations for all proteins from the Uniprot database.² To performed function enrichment analysis of differential human proteins, we utilized g:Profiler³ (Raudvere et al., 2019). In the g:Profiler analysis, the Benjamini-Hochberg FDR method was employed for multiple testing correction, with a threshold set at 0.0001. No electronic GO annotations were allowed. We analyzed the Gene ontology of molecular function (GO_MF), cellular component (GO_CC), and biological process (GO_BP), as well as the biological pathways of KEGG, Reactome (REAC), and WikiPathways (WP).

2.10 Linear discriminant analysis effect size analysis

To identify the key phyla that contribute significantly to the differences between gallstone disease and polyps, we conducted linear discriminant analysis effect size (LEfSe) analysis using the online Galaxy pipeline.⁴ The factorial Kruskal-Wallis test was employed to calculate LDA Effect Size, with an alpha value of 0.05 for the test among classes. Discriminative features were determined based on a logarithmic LDA score threshold of 2.0. Prior to the analysis, the LFQ-intensity values were converted to percentages. This conversion involved subtracting the LFQ-intensity value of a taxon from the total sum of LFQ-intensity values of all taxa in a specific sample. In cases where values were missing, they were substituted with zero.

2.11 Other bioinformatic analyses

Pearson correlation between bile samples were analyzed using the R package Corrplot v0.84. For hierarchical clustering and heatmap generation, the R package pheatmap v1.0.12 was employed with Euclidean distance and complete linkage. To conduct partial Least-Squares Discriminant Analysis (PLS-DA) using taxonomic data, we utilized R package mixOmics v6.22.0.

2.12 Metagenomic analysis of 16S rRNA sequencing datasets

We downloaded the raw sequences (filtered and clipped, fastq format) of six datasets of 16S rRNA sequencing from the Sequence Read Archive (SRA) of NCBI database. The accessions for these datasets were PRJNA439241, PRJEB15501, PRJEB12755, PRJNA580086, PRJNA793871, and PRJNA543184. These datasets consisted of 27, 126, 59, 45, 122, and 100 human bile samples, respectively, which were subjected to 16S rRNA gene sequencing (Supplementary Table S3). To performed taxonomic classification, we utilized Kraken2 program⁵ and matched the sequences against the Kraken2-microbial database (September 2018, 30GB) (downloaded from https://lomanlab.github.io/mockcommunity/mc_databases.html). This database contained genomic sequences from archaea, bacteria, fungi, protozoa, viral and Univec_Core. Sequences from taxa other than bacteria were served as a performance control, ensuring that the 16S rRNA sequencing reads would not identify viral hits. The Kraken2-microbial database consisted of a total of 21,207 taxonomy nodes and had a table size of 5,631,915,627. Kraken2 employs an exact k-mer method to accurately and efficiently match and classify reads to the lowest common ancestor (LCA) of microbial genomes (Lu et al., 2022). The Kraken2 reports from different samples were further analyzed and integrated using Pavian⁶ (Breitwieser and Salzberg, 2020).

2.13 Dot blotting analysis

Proteins from 5 μL of bile samples were transferred to PVDF membrane (Merck Millipore, Shanghai, China) using the Bio-Dot microfiltration apparatus according to the manufacturer’s instructions (Bio-Rad Laboratories, Inc., United States). The immunoblot was conducted according to previously described methods (Wang et al., 2021). The CFI rabbit polyclonal antibody (Cat. XC5623) and C5 rabbit polyclonal antibody (Cat. XC8104) were purchased from Shanghai Outdo Biotech Co., Ltd.

2.14 Availability of data and materials

The raw MS data, which underpins the findings of this article, can be accessed at the ProteomeXchange Consortium⁷ via the iProX partner repository.⁸ The dataset identifier for this study is IPX0004827000 or PXD035915. Additionally, the dataset that supports the conclusions of this article is included within the article itself.

2.15 Statistics

Statistical significance was determined using a two-sided Welch’s t test. A p-value < 0.05 was considered statistically significant. For the g:Profiler analysis, the Benjamini-Hochberg FDR method was employed for multiple testing correction, with a threshold set at 0.0001. All measured values were reported as mean ± standard deviation (SD).

3 Results

3.1 Proteomics and metaproteomic analyses of human bile from gallstone and gallstone-free patients

We performed a comprehensive analysis of 40 human bile samples obtained from gallstone patients and gallstone-free controls (Figure 1A). To facilitate the analysis, we employed a two-stage database search strategy using the MaxQuant workflow integrated into the Metalab program (Cheng et al., 2017), as depicted in Figure 1B. In the first round of searching, we searched over 2.6 million MSMS spectra against the primary database. This first-stage searching identified 30,575 protein sequences, which were subsequently used as a secondary database for the second round searching by MaxQuant (Tyanova et al., 2016). MaxQuant identified 23,033 peptides, which were then subjected to taxonomy LCA analysis by Metalab. Ultimately, we identified 2,836 protein groups, including 87 microbial and 2,749 human proteins. The correlation coefficients for the bile components between polyp samples were found to be 0.75 ± 0.16. Similarly, the correlation coefficients between chronic CHL samples were 0.73 ± 0.13, while the correlation coefficients between acute CHL samples were 0.72 ± 0.15. These results indicate that there is a greater degree of variation in individual bile components within samples from patients with CHL compared to the samples from individuals with polyps (Figure 1C).

3.2 Microbial taxa in human gallstone bile identified by microbial peptides

We successfully identified 711 distinct peptides that were specific to microorganisms. Utilizing these peptides, the LCA algorithm enabled the identification of 142 microbial taxa, each supported by a minimum of 3 distinct peptides (Supplementary Table S2). PLS-DA analysis using taxonomic data shared features between chronic and acute CHLs, while the polyp samples exhibited distinct characteristics from CHL (Figure 1D). We re-analyzed six 16S-rRNA gene sequencing datasets comprising a total of 479 bile samples and validated the presence of 86 of these bacterial phyla in bile samples, as shown in Supplementary Tables S2, S3. Additionally, through data mining of public references, we discovered that the majority of the taxa (139 out of 142) were found in various human body sites, except for bile (Supplementary Table S2), such as the gut/feces, scalp/skin, milk, oral tissues, vagina, etc. (Figure 2D). Among all the distinct peptides identified in our samples (a total of 23,033), microbial-specific peptides accounted for approximately 3% of the total (Figure 2A). This suggests that the microbial content in bile is relatively low in abundance.

FIGURE 2

Figure 2. Identification of microbes from gallstone and polyp bile samples. (A) The proportion of unique peptides identified in human, fungi, bacteria, archaea, and virus. (B) The structure and abundance of the top 10 phyla in both polyp and gallstone bile samples were analyzed using the intensity of identified microbial peptides. CHL, cholecystitis. (C) The structure and abundance of the top 10 microbial species in polyp and gallstone bile samples. (D) The taxa we identified were found in other Human body sites through data mining of public references. The details were deposited in Supplementary Table S2. (E) The tandem MS spectrum of peptide LYCEEVGWVICGIK identified from lepA protein (LEPA_BUCCC) of Buchnera aphidicola. (F) The Sankey plot shows the metataxonomic analysis of human bile samples. The dataset was retrieved from GenBank with the BioProject accession number PRJNA580086 (https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA580086). In this project, 45 biliary bile samples were analyzed using 16S rRNA sequencing. We re-analyzed the dataset against kraken2-microbial database (https://lomanlab.github.io/mockcommunity/mc_databases.html) using kraken2 program (Lu et al., 2022). The red arrowhead represents the clade of Buchnera aphidicola, which has been detected in 44 out of 45 bile samples, comprising a total of 7,313 reads.

Our analysis at the Phylum level revealed the presence of several bacterial phyla in bile samples, including Proteobacteria, Firmicutes, Actinobacteria, Bacteroidetes, Tenericutes, Chlorobi, Cyanobacteria, Chloroflexi, and Spirochaetes (Figure 2B; Supplementary Table S2). The identification of some phyla in the bile samples, like Chlorobi and Spirochaetes, is consistent with previous studies (Saltykova et al., 2016; Molinero et al., 2019). Among the fungal phyla, Ascomycota was the most abundant in bile samples. Euryarchaeota was the dominant Phylum among Archaea, and it has also been found in the human oral cavity (Dewhirst et al., 2010). In the human gut, Ascomycota is the most common fungal phylum, while Euryarchaeota is the most common archaeal Phylum (Xiong et al., 2022). We also identified the presence of the Genus Candida, belonging to the Phylum Ascomycota (Supplementary Table S2). Candida spp. is an infectious agent associated with CHL and has serious prognostic implications (Diebel et al., 1996). In terms of Viruses, Pisuviricota exhibited the highest abundance (Figure 2B; Supplementary Table S2). Pisuviricota is a Phylum of RNA viruses, including the Genus Cosavirus, which commonly causes acute gastroenteritis infections in pediatric patients (Razizadeh et al., 2021).

At the Class level, the most abundant Classes for Bacteria, Eukaryota, Archaea and viruses were Gammaproteobacteria, Saccharomycetes, Methanomicrobia, and Pisoniviricetes, respectively (Supplementary Figure S1A). Supplementary Figures S1B–D show the most abundant taxa at the Order, Family, and Genus levels. Among these taxa, Streptococcus and Clostridium are commonly found in gallstone samples and have been associated with acute CHL (Ploszaj et al., 2021; Shigemori et al., 2022). Notably, we identified that Staphylococcus aureus as the third most abundant Species in the bile samples (Figure 2C). Both Streptococcus and S. aureus are known to be β-glucuronidase-producing bacteria, and infections caused by these bacteria have been associated with the development of gallstone formation (Hancke and Marklein, 1983).

At the Species level, Vibrio vulnificus was identified as the most prevalent species, and its identification was supported by the metataxonomic analysis of bile datasets (Figure 2C; Supplementary Table S3). V. vulnificus is a pathogenic agent typically associated with exposure to seawater or consumption of contaminated seafood. Interestingly, its presence in human gallstone bile fluid suggests a potential link between V. vulnificus infection and gallstone formation (Zhao et al., 2015). While there is no specific evidence linking Schizosaccharomyces pombe to CHL, this yeast Species has been identified in the gut mycobiota of both healthy individuals and colorectal cancer patients (Chin et al., 2018). Buchnera aphidicola is typically recognized as an obligate endosymbiont of aphids, establishing a mutualistic relationship with their hosts. However, it is noteworthy that peptides of B. aphidicola were detected in all the bile samples in our study. This unexpected finding suggests the presence of B. aphidicola or its remnants in the bile samples, potentially indicating a broader ecological role or a possible association with gallstone disease beyond its traditional symbiotic relationship with aphids. Figure 2E displayed a representative tandem MS spectrum of the peptide LYCEEVGWVICGIK, which was identified from the lepA protein (LEPA_BUCCC) of B. aphidicola. Metataxonomic analysis of 16S rRNA gene sequencing datasets further validates its presence in human bile samples (Figure 2F). In fact, a study utilizing metagenomic shotgun sequencing has confirmed the presence of B. aphidicola in the human skin microbiome (Oh et al., 2016). The fifth Species we have identified is Clostridium perfringens, a bacterium that produces β-glucuronidase, and this particular bacterium has been linked to the formation of pigment gallstones (Skar et al., 1986). Eremothecium gossypii is a filamentous fungus known to be an infectious agent capable of producing riboflavin and glucosylceramide, and the accumulation of glucosylceramide it produces could contribute to cholelithiasis (Taddei et al., 2010). It is known that riboflavin plays a role in lipid metabolism and can influence cholesterol accumulation. Neosartorya fumigata (Aspergillus fumigatus), an opportunistic fungal pathogen and a major allergen, is the most common cause of invasive aspergillosis, particularly for immunocompromised individuals (O'Gorman et al., 2009).

3.3 Microbial markers of gallstone disease

In the CHL group, we observed significant higher abundance of 34 phyla, including 3 fungi, 5 archaea, 12 bacteria, and 14 viruses, compared to the polyp group (Figure 3A; Supplementary Table S2). We identified 44 differential phyla between acute CHL and polyp bile samples, with 29 phyla common to the CHL and polyp comparison (Supplementary Figure S2A; Supplementary Table S2). Additionally, we found 27 differential phyla between chronic CHL and polyp bile samples, all of which were common to the CHL vs. polyp comparison (Supplementary Figure S2B; Supplementary Table S2). However, only four differential phyla were identified between acute and chronic CHL (Supplementary Figure S2C; Supplementary Table S2). To identify potential microbial biomarkers for CHL, we performed a LEfSe analysis. Interestingly, four clades were identified as microbial biomarkers in CHL (LDA score > 2), including the Thermococcaceae clade, Helicobacteraceae clade, Methanosarcina clade, and Reoviridae clade (Figures 3B,C). Among the Archaea we identified, Methanosarcina and the Thermococcaceae clade have been found in the human skin microbiome (Oh et al., 2016). In the bacterial phyla upregulated in CHL, Helicobacteraceae is a Family that belongs to the Order Campylobacterales and the Class Epsilonproteobacteria. Interestingly, a human cohort study reveals that Helicobacter infection is associated with an increased risk of gallstone disease (Cen et al., 2023). Our analysis detected Helicobacteraceae peptides in 39 out of 40 bile samples (Figure 3D). Furthermore, upon re-analyzing the dataset PRJNA439241, we detected 16S rRNA reads of Helicobacteraceae in all 27 bile samples (Figure 3E; Supplementary Tables S2, S3). Within the Chlorobium clade, Chlorobia, Chlorobiaceae, Chlorobiales, and Chlorobium have been identified in human skin microbiome (Oh et al., 2016), while Chlorobi (Phylum) was identified in bile (Molinero et al., 2019). The presence of Beijerinckiaceae in bile was confirmed by reanalysis of 16S rRNA datasets of bile (4 out of 6) (Supplementary Tables S2, S3), and it was also identified in human scalp (Futatsuya et al., 2020). Cyanobacteria was widely distributed in the human body, not only in gallstone bile (Wu et al., 2013), but also in colon (Stearns et al., 2011), scalp (Futatsuya et al., 2020), and stool (Morales et al., 2022). Notably, Cyanobacteria, well-known for its biofilm production capability, may contribute to gallstone formation within the bile environment (Maeda et al., 2021). Rhodobacterales was found to be a bile biomarker in a study analyzing bile and gallstone samples (Hu et al., 2023). Mycoplasma was the only bacteria downregulated in CHL. As a human pathogen, Mycoplasma has been detected in various human locations, including bile, subgingival and buccal mucosa (Griffen et al., 2012; Kraal et al., 2014; Chen et al., 2019). Its presence in bile was further validated by reanalysis of bile 16S rRNA datasets (Supplementary Tables S2, S3). We also revealed several differential virus phyla in CHL samples. Bamfordvirae, Nucleocytoviricota, Megaviricetes, Imitervirales, and Mimiviridae were identified based on the same set of peptides (Supplementary Table S2). Nucleocytoviricota, although not identified in our analysis of 16S rRNA datasets, has been found in the human ocular surface (Shao et al., 2023), while Mimiviridae or Minivirus has been identified in skin microbiome (Oh et al., 2016) or pulmonary samples (Saadi et al., 2013). The Reoviridae clade, as well as Pisuviricota (Pisoniviricetes), belongs to the Kindom Orthornavirae. Reoviridae (Duplornaviricota, Reovirales, Resentoviricetes, Rotavirus) has been commonly implicated in biliary atresia, a childhood disease (Averbukh and Wu, 2018). Among the fungal taxa, Debaryomyces has been previously identified in the human digestive tract and associated with Crohn’s disease (Olaisen et al., 2022), while Sordariomycetes (Hypocreales) was found in skin microbiome (Oh et al., 2016).

FIGURE 3

Figure 3. Taxon changes associated with gallstone formation determined by the microbial peptides. (A) Hierarchical cluster analysis of the microbial taxa significantly changed in gallstone bile samples (n = 31) comparing to polyp bile samples (n = 9). (B,C) Linear discriminant analysis Effect Size (LEfSe) analysis of microbial biomarkers of CHL (n = 31) and polyp (n = 9) using an LDA threshold of 2. The numbers indicate the major marker clades in CHL. (D) The abundance of Helicobacteraceae between CHL (n = 31) and polyp (n = 9) patients based on the proteomic data. The p values were calculated using Welch’s t test. MS, mass spectrometry. (E) The abundance of Helicobacteraceae between CHL (n = 14) and healthy (n = 13) patients based on the 16S rRNA data (dataset PRJNA439241). 16S, 16S rRNA dataset analysis. Chole, cholecystitis. Healt, healthy.

3.4 Identification of microbial protein groups

We identified a total of 87 microbial proteins (Figure 4A; Supplementary Table S4), with 47 (54%) of them associated with the Phylum Proteobacteria. It is noteworthy that the ability of bacteria to cause biliary infection plays a crucial role in gallstone formation by promoting the production of biofilms, also known as slime or glycocalyx, which facilitates bacterial colonization and solidification of gallstones. Among the identified microbial proteins, several are known to be involved in biofilm formation. We observed a gradual upregulation of quinidine resistance protein 3 (QDR3), a quinidine drug resistance protein from Saccharomyces, from chronic CHL to acute CHL when compared to polyp samples (Figures 4B–E). This protein is instrumental in biofilm formation and enhances the virulence of fungal pathogens (Shah et al., 2014). Similarly, ompA, an outer membrane protein found in Escherichia. coli, is crucial for the complete formation of certain bacterial biofilms in a bile environment (Tsai et al., 2020). The nucleoside diphosphate kinase (ndk) protein in Pseudomonas not only influences cytotoxicity and host pathogenicity (Yu et al., 2016), but also plays a role in biofilm production (Sikarwar et al., 2022). Pseudomonas ndk is a critical host-responsive gene that regulates bacterial virulence during infection, as its expression was downregulated during acute pneumonia (Yu et al., 2016). The phosphate-binding protein pstS, identified in Pasteurella (Supplementary Table S4), has been shown to play an important role in biofilm formation in its homolog in Pseudomonas (Neznansky et al., 2014). N-acetylneuraminate pyruvate-lyase (nanA) in Clostridium is a sialic acid aldolase, and sialic acid is associated with bacterial pathogenicity. The ortholog of this enzyme in other bacteria has been demonstrated to participate in biofilm formation (Di Pasquale et al., 2016). Formate acetyltransferase 1 (pfIB) in E. coli also contributes to biofilm formation by regulating acytyl-CoA concentrations (Mugabi et al., 2012). Similarly, the chaperone protein dnaK in Clostridium enhances biofilm formation (Jain et al., 2017).

FIGURE 4

Figure 4. The microbial proteins identified from gallstone and polyp bile samples. (A) Abundance rank of identified proteins. The intensity of each protein was determined by summing up the LFQ values across all samples. Microbial proteins were highlighted. (B) The abundance of differential microbial proteins between CHL (n = 31) and polyp (n = 9) patients. (C) The differential microbial proteins between acute CHL (n = 14) and polyp (n = 9) patients. (D) The differential microbial proteins between chronic CHL (n = 17) and polyp (n = 9) patients. (E) Three microbial proteins upregulated in CHL compared to polyp. p values were calculated using Welch’s t test. (F) The differential microbial proteins between acute CHL (n = 14) and chronic CHL (n = 17) patients.

We have identified several microbial proteins that enhance bacterial virulence and gallstone formation by counteracting host oxidative stress. One of these proteins is superoxide dismutase (SOD), which acts as the primary scavenger of host oxidative superoxide anion. We found that Pseudomonas iron-cofactored SOD [Fe] (sodB/SODF_PSEPK) was upregulated in both acute and chronic CHL compared to polyp samples (Figures 4B–E). Pseudomonas aeruginosa has been revealed as the dominant agent associated with gallstone disease (Peng et al., 2015) and acute CHL (Sakurai et al., 1992). sodB plays a vital role in enabling Pseudomonas to survive against host oxidative stress (Cavinato et al., 2020). KatG, a catalase-peroxidase, plays a crucial role in protecting the biofilm of certain bacteria by neutralizing external oxidative stress, thereby contributing to biofilm growth on cholesterol gallstones (von Beek et al., 2019). Rubrerythrin (rbr) in Clostridium functions as a scavenger of oxygen radicals, aiding in the protection of bacteria against oxidative stress (Jean et al., 2004). The serine protease htrA, implicated in the virulence of many pathogens including Staphylococcus (Rigoulay et al., 2005), contributes to pathogenicity by regulating the expression of essential virulence factors and stress proteins. Salmonella has been associated with gallstone development due to its ability to form biofilms on the cholesterol-coated surface of gallstone (Gonzalez et al., 2019). We identified a protein in Salmonella, the alkyl hydroperoxide reductase subunit C (ahpC). This peroxiredoxin reduces organic peroxide, protecting Salmonella from host oxidative stress and enhancing its virulence (Hebrard et al., 2009).

We also identified microbial proteins that contribute to bacteria virulence and potential play a role in gallstone formation. Similar to QDR3, Buchnera lepA exhibited a gradual upregulation from chronic to acute CHL compared to polyp samples (Figures 4B–E). LepA also known as elongation factor 4 (EF-4), is a GTP-binding protein that acts a virulence factor in several bacteria (Kida et al., 2011). The RtxA gene, which encodes heat-stable enterotoxin, has been identified as an essential virulence factor of Vibrio vulnificus (Lee et al., 2007). The pckA gene, which encodes phosphoenolpyruvate carboxykinase (PEPCK), acts as a pathogenic factor for virulent bacteria (Liu et al., 2005). Tuf, the elongation factor Tu (EFTU) in Pseudomonas, was identified by 13 unique or razor peptides in bile (Supplementary Table S4). Tuf can serve as a virulence factor for Pseudomonas, enabling it to evade human complement attack by binding to the human complement regulator Factor H and plasminogen (Kunert et al., 2007). Pseudomonas, E. coli and Staphylococcus are biofilm-forming bacteria that have been identified in bile stone samples (Tajeddin et al., 2016). The capability of Pseudomonas to form biofilms on the surface of implanted gallstones has been confirmed in animal models (Sung et al., 1991). TpiA, a key metabolic triosephosphate isomerase of pathogenic bacteria, affects virulence, antibiotic resistance, and fitness within the host (Paterson et al., 2009; Xia et al., 2020), and its upregulation in CHL may indicate the adaptation of invading bacteria to the altered host environment.

The identification of microbial proteins in our study encompasses transporters that counteract bile salt stress or drug attack, thereby facilitating microbial survival in bile environment. One of these proteins is the outer membrane protein TolC, which functions as a crucial pump in E. coli, enabling the efflux of various substances like bile salts (Rosner and Martin, 2009). Additionally, QDR3, a multidrug transporter in yeast, has been identified in bile samples.

In comparing acute and chronic CHL, we observed differential expression of three microbial proteins: cysI (CYSI_NITWN), gltX1 (SYE1_RICTY), and MJ0811 (Y811_METJA) (Figure 4F). These findings suggest that there is only a minimal difference in microbial expression between the two sample types.

3.5 Enhanced host inflammatory in cholelithiasis

We identified a total of 2,749 host proteins across all bile samples (Supplementary Table S5), including 67 upregulated and 188 downregulated proteins in the cholelithiasis group (n = 31) compared to the polyp group (n = 9) (CHL_vs_polyp), with an absolute log2-fold change ≥1 and a p value < 0.05 (Supplementary Table S5). We then performed a functional enrichment analysis using g:Profiler (Raudvere et al., 2019). The analysis revealed that the upregulated proteins in CHL were significantly enriched in serine-type endopeptidase activity (p_adjusted = 7.0E-11), heparin binding (p_adjusted = 1.3E-07), and complement binding (p_adjusted = 0.0009) (Supplementary Table S6; Figure 5A). Moreover, the GO Biological Process (GO_BP) terms in CHL were primarily associated with various immune response processes, including humoral immune response (p_adjusted = 4.0E-19) and defense response (p_adjusted = 2.2E-19) (Supplementary Table S6). These proteins comprised immunoglobulins (Igs), C4BPA, RNASE3, CTSG, AZU1, ELANE, PRTN3, S100A9, DEFA3, FGA, and others. Additionally, the presence of complement-activation-associated proteins like C4BPA, C1QA, C5, CFP, CD5L, IGHG2, A2M, and CFH, was increased in CHL bile (Figure 5B). KEGG analysis indicated that the upregulated proteins in CHL were predominantly involved in pathways such as complement and coagulation cascades (p_adjusted = 1.4E-08) and neutrophil extracellular trap formation (p_adjusted = 6.1E-05). Consistently, Reactome (REAC) analysis revealed that CHL-upregulated proteins were significantly enriched in the innate immune system (p_adjusted = 3.1E-12), complement cascade (p_adjusted = 6.6E-07) and antimicrobial peptides (p_adjusted = 3.5E-08). Furthermore, the WikiPathways (WP) analysis indicated an overrepresentation of the complement system (p_adjusted = 5.2E-11) among the upregulated proteins in CHL (Supplementary Table S6; Figure 5A). In the GO cellular component (GO_CC) analysis, the upregulated proteins were primarily enriched in the extracellular space (p_adjusted = 3.0E-23), while the downregulated proteins showed enrichment in both the extracellular space (p_adjusted = 1.8E-44) and the cytoplasm (p_adjusted = 3.7E-14).

FIGURE 5

Figure 5. The host proteins identified from gallstone and polyp bile samples. (A) Gene ontology (GO) enrichment was performed on the upregulated human proteins in CHL compared to polyp using the g:Profiler tool (https://biit.cs.ut.ee/gprofiler/gost) (Raudvere et al., 2019). (B) Differential expression of complement proteins between CHL and polyp samples. p values were calculated using unpaired Welch’s t test. (C) Two mucin proteins were upregulated in CHL compared to polyp. (D) GO enrichment of the downregulated human proteins in CHL compared to polyp. GO_MF, Gene ontology of molecular function; GO_CC, cellular component; GO_BP biological process; REAC, Reactome; WP, WikiPathways. (E) Proteins associated with bile secretion were downregulated in CHL compared to polyp.

We identified 74 increased and 362 decreased host proteins in acute CHL (n = 14) compared to polyp (n = 9) (acute_vs_polyp) (Supplementary Figure S3; Supplementary Table S5). Similarly, we identified 51 upregulated and 95 downregulated host proteins in chronic CHL (n = 17) compared to polyp (n = 9) (chronic_vs_polyp) (Supplementary Table S5; Supplementary Figure S4). The enrichment analyses of these differential proteins yielded results to the CHL vs. polyp comparison, suggesting a potential role of complement activation in the cholelithiasis process. To validate the expression of selected complement proteins, we conducted a dot plot analysis. As shown in Supplementary Figure S5, consistent with the proteomic quantitation results, the expression of CFI remained constant between polyp and gallstone bile samples, while C5 was found to be increased in gallstone bile samples.

These findings indicate the presence of an inflammatory feature in CHL and the hyperactivation of both adaptive and innate immune pathways in the host bile environment.

3.6 Increased host defense against bacterial infection

Five GO:BP terms exhibited specific enrichment in proteins upregulated in CHL, including defense response to bacterium (p_adjusted = 6.8E-14) response to bacterium (p_adjusted = 2.4E-11), antimicrobial humoral response (p_adjusted = 2.1E-13), antibacterial humoral response (p_adjusted = 9.5E-10), and defense response to fungus (p_adjusted = 5.1E-10) (Supplementary Table S6). Notably, 17 human proteins were upregulated in the defense response to bacterium, including C4BPA, RNASE3, CTSG, MMP9, RETN, AZU1, EPX, MPO, ELANE, SERPINB4, PRTN3, LCN2, CRP, S100A9, S100A8, HNRNPA0, DEFA3, APCS, FGA, NT5C3A, HRG, F12, LTF, FGB, FGG, RTN4R, C1QA, C5, CFP, TTN, IGHG2, A2M, TF, PPP1CB, and CFH. Furthermore, the upregulated proteins in CHL demonstrated enrichment in Staphylococcus aureus infection (p_adjusted = 0.001) based on the KEGG analysis, and in antimicrobial peptides (p_adjusted = 3.5E-8) according to the REAC analysis. These results strongly suggest a crucial role for the antimicrobial immune response in CHL.

3.7 Coagulation and other pathways activated in CHL

We observed a significant upregulation of proteins in CHL compared to polyp, specifically those involved in blood coagulation and fibrin clot formation (p_adjusted = 7.1E-08, GO_BP), formation of fibrin clot (clotting cascade) (p_adjusted = 4.5E-07, REAC), common pathway of fibrin clot formation (p_adjusted = 2.9E-05, REAC), and blood clotting cascade (p_adjusted = 0.0001, WP) (Supplementary Table S6; Figure 5A). Similar enrichments were observed in the differential proteins of acute_vs_polyp and chronic_vs_polyp. Notably, the proteins involved in these processes were FGA, F12, FGB, FGG, FN1, PRTN3, and A2M. Furthermore, we found an overrepresentation of the folate metabolism pathway in CHL (p_adjusted = 4.1E-04), which includes components MPO, CRP, FGA, FGB, and FGG. Additionally, the interleukin-17 (IL-17) signaling pathway (p_adjusted = 8.2E-06) was enriched in the upregulated protein in CHL, including MMP9, LCN2, S100A9, S100A8, MUC5AC, MUC5B, and MMP1. It is worth mentioning that among the seven mucin (MUC) proteins, MUC5A and MUC5B exhibited a significant increase in gallstone bile compared to polyp bile (Figure 5C). These findings are consistent with previous findings, which suggest that high-molecular-weight glycoproteins, such as intraluminal mucin deposition, may contribute to gallstone formation (Yoo et al., 2016).

3.8 Host catabolic metabolism activity was suppressed in gallstone bile

In gallstone bile (n = 31) compared to polyp (n = 9), a total of 188 human proteins exhibited downregulation. These proteins showed enrichment in GO_MF terms related to catalytic activity (299 proteins, 53%) (p_adjusted = 2.7E-14) and hydrolase activity (p_adjusted = 4.3E-14) (Supplementary Table S6; Figure 5D). GO_BP analysis revealed that many of the downregulated host proteins were associated with “metabolic” or catabolic,” such as catabolic process (p_adjusted = 1.6E-16, 71 proteins, 37.8%), small molecule metabolic process (p_adjusted = 1.1E-11), and organic substance catabolic process (p_adjusted = 2.0E-14). This finding indicates that gallstone disease significantly suppresses metabolic activity, particularly catabolic metabolism. Moreover, this conclusion was further supported by KEGG, REAC, and WP analyses (Figure 5D). Additionally, GO_CC analysis revealed that the downregulated proteins were not only enriched in extracellular/secreted proteins but also in cytoplasmic proteins, which aligns with the fact that many metabolism-associated proteins are primarily located in the cytoplasm.

3.9 Bile secretion and transportation were impaired in gallstone formation

The proteins downregulated in CHL, in comparison to polyps, demonstrated enrichment in the KEGG pathway of bile secretion (p_adjusted = 0.039), which included ATP1B1, ABCC2, ABCB11, ABCB4, ABCB1, and ABCG8 (Figures 5D,E; Supplementary Table S6). Furthermore, we observed a significant downregulation of host proteins related to disorders of bile acid synthesis and biliary transport (WP:WP5176) in gallstone bile (p_adjusted = 0.0001). Notable proteins in this category were AKR1D1, ABCC2, HSD17B4, ABCB11, ABCB4, and ATP8B1 (Supplementary Table S6). Additionally, the downregulated proteins in gallstone bile showed enrichment in bile acid and bile salt transport (p_adjusted = 0.004), which included AKR1C1, ABCC2, ABCB11, and ATP8B1. These findings indicate a disruption in bile secretion in CHL compared to polyp, potentially leading to bile supersaturation and the formation of cholesterol-type gallstones.

3.10 Catabolic processes were further suppressed in acute CHL compared to chronic CHL

We identified 8 upregulated and 95 downregulated proteins in acute CHL (n = 14) compared to chronic CHL (n = 17) (Supplementary Figure S6A; Supplementary Table S6). No enrichment was observed among the upregulated proteins (Supplementary Figure S6B). However, the downregulated proteins exhibited enrichment in metabolic pathways, particularly the catabolic process (p_adjusted = 0.0011), cellular catabolic process (p_adjusted = 0.0006), and cytoplasmic translation (p_adjusted = 3.2E-07) (Supplementary Figure S6C; Supplementary Table S6). These findings indicate a further suppression of catabolic activities in acute CHL compared to chronic CHL.

4 Discussion

In this study, we performed an integrated proteomic and metaproteomic analysis of human bile samples obtained from individuals with gallstone and polyp disease. Utilizing LCA algorithm, we performed a taxonomic analysis of the microbiome based on 711 identified microbial peptides, and identified 142 microbial taxa. The lower abundance of microbial proteins in bile contributed to the relatively small identification rate of microbial peptides compared to human peptides. While this may affect the composition analysis of microbial taxa, the lack of appropriate enrichment methods for microbial peptides limits the improvement in identification. To validate the presence of these taxa in bile samples, we re-analyzed six 16S rRNA gene sequencing datasets, confirming the presence of 60% of the identified taxa. Furthermore, publicly available references demonstrated that the majority of these taxa were also found in different body sites of humans. These findings provide robust evidence for the reliability of our metaproteome results and underscore the potential of metaproteomics in detecting previously undiscovered microbial entities that may have been missed by DNA sequencing techniques. Notably, our metaproteomic analysis uncovered the presence of microbes, such as B. aphidicola, an endosymbiont commonly found in aphids, which were previously not considered to have a relationship with humans. The predominant bacterial phyla identified in gallstone bile through metaproteomic analysis were Proteobacteria, Firmicutes, Actinobacteria, and Bacteroidetes, displaying similarities to the composition revealed by 16S rRNA sequencing, albeit with differences in taxonomic order (Molinero et al., 2019). Several factors, including the distinct nature of nucleotides and proteins, sample variation, and potential limitations in data acquisition, may contribute to this discrepancy.

Gallstones can be categorized into three major types: cholesterol stones, brown pigment stones, and black pigment stones (Wang Y. et al., 2018). In our study, we identified taxa known to produce β-glucuronidase, such as Streptococcus, S. aureus, and Clostridium. Bacteria-derived β-glucuronidase hydrolyzes bilirubin conjugates in gallbladder bile, leading to the release of free bilirubin and glucuronic acid, which preciptates as insoluble salts with ionized calcium. This process is considered the primary mechanism underlying the formation of brown pigment gallstones (Wang Y. et al., 2018). Through LEfSe analysis, we identified marker taxa in gallstone disease bile, including the Helicobacteraceae clade, which was present in almost all bile samples. Helicobacter infection is commonly considered as a potential risk factor for gallstone disease (Cen et al., 2023). However, experiments using animal models of H. pylori infection to study cholesterol gallstone formation have not provided conclusive evidence (Maurer et al., 2006). Bacterial biofilm plays a significant role in bacteria colonization and gallstone solidification, serving as an additional pathway for gallstone formation (Wang Y. et al., 2018). Bacteria with slime activity, such as Salmonella, Enterococcus, and Neisseria, have been observed to promote gallstone solidification (Sharma et al., 2020). Our identification of bacteria responsible for biofilm formation includes Cyanobacteria (Maeda et al., 2021), Pseudomonas (Sikarwar et al., 2022), E. coli (Tsai et al., 2020), and Clostridium (Jain et al., 2017).

In our analysis, we successfully identified 87 microbial proteins that shed light on the molecular mechanisms underlying gallstone formation. Firstly, we discovered several microbial proteins involved in biofilm production, such as ndk, QDR3, ompA, pstS, dnaK, nanA, and pflB. These proteins originated from microbes such as Pseudomonas, Saccharomyces, E. coli, Pasteurella, or Clostridium. Secondly, we identified microbial proteins that enhance bacterial virulence and contribute to gallstone formation by counteracting host oxidative stress. These proteins include sodB, KatG, rbr, HtrA, and ahpC. Gallstone patients exhibit high levels of oxidative stress in the gallbladder mucosa (Geetha, 2002). Oxidative stress facilitates the rapid oxidation of bilirubin and free radical polymerization, which play a crucial role in gallstone nucleation and formation (Sanikidze and Chikvaidze, 2016). Thirdly, we discovered several microbial proteins that contribute to bacterial virulence and potentially induce gallstone formation, such as Buchnera lepA, RtxA gene, pckA gene, Tuf, and TpiA. These proteins have diverse roles in bacterial virulence, toxins production, immune evasion, biofilm formation, and adaptation to host environment changes, which have been addressed through gene expression studies, bile stone samples analysis, and animal models. Lastly, we identified microbial proteins, including transporters like TolC and QDR3, that enable microbial survival in the bile environment by counteracting bile salt stress or drug attacks.

We identified a total of 2,749 human proteins, among which 67 were significantly upregulated and 188 were downregulated in gallstone bile compared to gallstone-free bile. The host protein profile revealed an enhanced immune or inflammatory response associated with gallstone and CHL diseases. The upregulated human bile proteins were predominantly enriched in immune or inflammation-related biological processes or pathways, indicating the inherent inflammatory state of CHL. Notably, the complement system and antimicrobial peptides were overrepresented in CHL bile, suggesting excessive activation of the innate immune system. This finding is consistent with gene RNA-seq analysis conducted on both gallstone and polyp patients (Li et al., 2014). The increased levels of these enzymes may indicate a pancreatic leak in the context of CHL inflammation and pancreatitis. Additionally, we observed that the increased host proteins in gallstone bile are associated with resistance against microbial infection. Gallbladder S100A9 and S100A8 have been implicated in proinflammatory functions in CHL, particularly during the chronic stage (Szmyt et al., 2013).

In addition to the previously mentioned pathways, several other crucial pathways were found to be overrepresented in CHL, including blood coagulation, folate metabolism, and the IL-17 pathway. In cases of inflammatory gallstone disease, bleeding and the formation of fibrin clots may occur in the biliary tract (Chang et al., 2020). WP analysis indicated a significant alteration in folate metabolism, potentially implicating its involvement in gallstone disease (Worthington et al., 2004). Notably, among the proteins associated with folate metabolism, MPO exhibited a significant increase during gallstone development and mucus layer thickening in a mouse model (39). Furthermore, C-reactive protein (CRP), a marker of chronic inflammation, showed a significant association with a higher risk of gallstone disease in population studies (Barahona Ponce et al., 2021). The activation of the IL-17 pathway, which leads to inflammation and plays a significant role in cholestasis potentially induced by gallstone disease, was also observed (Tedesco et al., 2018). IL-17 is a proinflammatory cytokine, and CD4-positive Th17 cells that secrete IL-17 have been implicated in the chronic inflammation of bile ducts (Harada and Nakanuma, 2010).

In our analysis of downregulated proteins, we observed a significant suppression of metabolic activities, particularly catabolic metabolism, and transport activities in gallstone bile. Furthermore, the downregulated proteins were enriched in pathways related to bile secretion, bile acid synthesis, and biliary transportation. These changes contribute to the gradual accumulation of cholesterol and subsequent deposition of gallstones. We also identified host proteins involved in the formation of cholesterol-type gallstones. Notably, proteins associated with cholesterol transportation exhibited significant downregulation in gallstone bile. ABCG5/ABCG8, which function as hepatic and intestinal cholesterol transporters, form dimers to export sterols out of enterocytes and hepatocytes (Rebholz et al., 2018). Their downregulation leads to decreased phosphatidylcholine levels and increased cholesterol concentration, thereby promoting gallstone formation. Similarly, downregulation of ABCB11 and ABCB4, which act as bile salt exporters and phosphatidylcholine translocators, respectively, reduce cholesterol solubility (Chan and Vandeberg, 2012; Wang Y. et al., 2018). The downregulation of AKR1C4 and AKR1D1, which are involved in the synthesis of bile acids and bile salts that enhance cholesterol solubility, may promote gallstone formation (Wang Y. et al., 2018). Additionally, downregulation of NPC1L1, ABCG5, ABCG8, and LDLR, which are involved in intestinal and hepatobiliary cholesterol absorption, reduces cholesterol reabsorption from bile. This leads to biliary cholesterol supersaturation and subsequent gallstone deposition (Mokhtar et al., 2022).

5 Conclusion

In this study, we conducted a proteomic and metaproteomic analysis to explore the microbial features of gallstone and gallstone-free bile samples. We identified microbes and microbial factors involved in β-glucuronidase activity, biofilm formation, virulence, and anti-oxidative stress, which may contribute to gallstone development. Gallstone bile, or cholelithiasis bile, exhibited an enhanced inflammatory molecular profile, including molecules associated with the innate immune system’s response to microbial infections. Overrepresentation of pathways related to blood coagulation, folate metabolism, and IL-17 was observed in gallstone bile. Conversely, host metabolic activities, particularly catabolic metabolism and transport activities, were significantly reduced in gallstone bile. Additionally, we found that acute cholelithiasis bile exhibited more pronounced impairment in metabolic activities compared to chronic cholelithiasis bile. In conclusion, our analysis provide insights into the dysbiosis of resident microbes and the molecular interplay between the microbiome and host, contributing to the understanding of gallstone formation.

Data availability statement

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

Author contributions

X-TY: Data curation, Formal analysis, Investigation, Visualization, Writing – review & editing. JW: Formal analysis, Methodology, Validation, Investigation, Writing – review & editing. Y-HJ: Formal analysis, Methodology, Investigation, Writing – review & editing. LZ: Formal analysis, Methodology, Investigation, Writing – review & editing. LD: Data curation, Resources, Investigation, Writing – review & editing. JL: Conceptualization, Resources, Data curation, Writing – review & editing. FL: Conceptualization, Formal analysis, Funding acquisition, Project administration, Resources, Supervision, Visualization, Writing – original draft, Writing – review & editing, Data curation.

Ethics statement

The studies involving humans were approved by Institutional Research Ethics Committee of the Institutes of Biomedical Sciences, Fudan University. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The work was supported by the National Natural Science Foundation of China (NSFC) (No. 81572833 and 32000505).

Acknowledgments

We kindly thank Ming-Qi Liu for his suggestion regarding data submission to ProteomeXchange. The bioinformatic analysis was supported by the Medical Research Data Center of Fudan University.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1276951/full#supplementary-material

Abbreviations

FDR, False discovery rate; GO, Gene ontology; GO_BP, GO biological process; GO_CC, GO cellular component; HP, Human phenotype ontology; LCA, Lowest common ancestor; LEfSe, LDA effect size; LFQ, Label-free quantitation; MUCs, Mucins; PLS-DA, Partial least-squares discriminant analysis; REAC, Reactome; RPLC-MS, Reverse-phase liquid chromatography-mass spectrometry; WP, WikiPathways.

Footnotes

1. ^https://proteowizard.sourceforge.io/

2. ^https://www.uniprot.org/

3. ^https://biit.cs.ut.ee/gprofiler/gost

4. ^http://huttenhower.sph.harvard.edu/galaxy/

5. ^https://ccb.jhu.edu/software/kraken2/

6. ^https://fbreitwieser.shinyapps.io/pavian/

7. ^http://proteomecentral.proteomexchange.org

8. ^https://iprox.org/

References

Averbukh, L. D., and Wu, G. Y. (2018). Evidence for viral induction of biliary atresia: a review. J. Clin. Transl. Hepatol. 6, 1–10. doi: 10.14218/JCTH.2018.00046

PubMed Abstract | CrossRef Full Text | Google Scholar

Barahona Ponce, C., Scherer, D., Brinster, R., Boekstegers, F., Marcelain, K., Garate-Calderon, V., et al. (2021). Gallstones, body mass index, C-reactive protein, and gallbladder Cancer: Mendelian randomization analysis of Chilean and European genotype data. Hepatology 73, 1783–1796. doi: 10.1002/hep.31537

PubMed Abstract | CrossRef Full Text | Google Scholar

Barbhuiya, M. A., Sahasrabuddhe, N. A., Pinto, S. M., Muthusamy, B., Singh, T. D., Nanjappa, V., et al. (2011). Comprehensive proteomic analysis of human bile. Proteomics 11, 4443–4453. doi: 10.1002/pmic.201100197

CrossRef Full Text | Google Scholar

Breitwieser, F. P., and Salzberg, S. L. (2020). Pavian: interactive analysis of metagenomics data for microbiome studies and pathogen identification. Bioinformatics 36, 1303–1304. doi: 10.1093/bioinformatics/btz715

PubMed Abstract | CrossRef Full Text | Google Scholar

Cavinato, L., Genise, E., Luly, F. R., Di Domenico, E. G., Del Porto, P., and Ascenzioni, F. (2020). Escaping the phagocytic oxidative burst: the role of SODB in the survival of Pseudomonas aeruginosa within macrophages. Front. Microbiol. 11:326. doi: 10.3389/fmicb.2020.00326

PubMed Abstract | CrossRef Full Text | Google Scholar

Cen, L., Wu, J., Zhu, S., Pan, J., Zhou, T., Yan, T., et al. (2023). The potential bidirectional association between Helicobacter pylori infection and gallstone disease in adults: a two-cohort study. Eur. J. Clin. Investig. 53:e13879. doi: 10.1111/eci.13879

PubMed Abstract | CrossRef Full Text | Google Scholar

Chambers, M. C., Maclean, B., Burke, R., Amodei, D., Ruderman, D. L., Neumann, S., et al. (2012). A cross-platform toolkit for mass spectrometry and proteomics. Nat. Biotechnol. 30, 918–920. doi: 10.1038/nbt.2377

PubMed Abstract | CrossRef Full Text | Google Scholar

Chan, J., and Vandeberg, J. L. (2012). Hepatobiliary transport in health and disease. Clin. Lipidol. 7, 189–202. doi: 10.2217/clp.12.12

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, C. W., Wu, Y. H., and Liu, H. L. (2020). Upper gastrointestinal bleeding with Hemobilia caused by gallstones. J. Emerg. Med. 58, e227–e229. doi: 10.1016/j.jemermed.2020.03.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, B., Fu, S. W., Lu, L., and Zhao, H. (2019). A preliminary study of biliary microbiota in patients with bile duct stones or distal cholangiocarcinoma. Biomed. Res. Int. 2019, 1–12. doi: 10.1155/2019/1092563

CrossRef Full Text | Google Scholar

Chen, C. H., Huang, M. H., Yang, J. C., Nien, C. K., Etheredge, G. D., Yang, C. C., et al. (2006). Prevalence and risk factors of gallstone disease in an adult population of Taiwan: an epidemiological survey. J. Gastroenterol. Hepatol. 21, 1737–1743. doi: 10.1111/j.1440-1746.2006.04381.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, K., Ning, Z., Zhang, X., Li, L., Liao, B., Mayne, J., et al. (2017). MetaLab: an automated pipeline for metaproteomic data analysis. Microbiome 5:157. doi: 10.1186/s40168-017-0375-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Chin, S. F., Megat Mohd Azlan, P. I. H., Mazlan, L., and Neoh, H. M. (2018). Identification of Schizosaccharomyces pombe in the guts of healthy individuals and patients with colorectal cancer: preliminary evidence from a gut microbiome secretome study. Gut Pathog. 10:29. doi: 10.1186/s13099-018-0258-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Choe, J. W., Lee, J. M., Hyun, J. J., and Lee, H. S. (2021). Analysis on Microbial Profiles & Components of bile in patients with recurrent CBD stones after endoscopic CBD stone removal: a preliminary study. J. Clin. Med. 10:3303. doi: 10.3390/jcm10153303

PubMed Abstract | CrossRef Full Text | Google Scholar

Craig, R., and Beavis, R. C. (2004). TANDEM: matching proteins with tandem mass spectra. Bioinformatics 20, 1466–1467. doi: 10.1093/bioinformatics/bth092

CrossRef Full Text | Google Scholar

Cui, S. J., Xu, L. L., Zhang, T., Xu, M., Yao, J., Fang, C. Y., et al. (2013). Proteomic characterization of larval and adult developmental stages in Echinococcus granulosus reveals novel insight into host-parasite interactions. J. Proteome 84, 158–175. doi: 10.1016/j.jprot.2013.04.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Dewhirst, F. E., Chen, T., Izard, J., Paster, B. J., Tanner, A. C., Yu, W. H., et al. (2010). The human oral microbiome. J. Bacteriol. 192, 5002–5017. doi: 10.1128/JB.00542-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Di Ciaula, A., Wang, D. Q., and Portincasa, P. (2018). An update on the pathogenesis of cholesterol gallstone disease. Curr. Opin. Gastroenterol. 34, 71–80. doi: 10.1097/MOG.0000000000000423

PubMed Abstract | CrossRef Full Text | Google Scholar

Di Pasquale, P., Caterino, M., Di Somma, A., Squillace, M., Rossi, E., Landini, P., et al. (2016). Exposure of E. coli to DNA-Methylating agents impairs biofilm formation and invasion of eukaryotic cells via down regulation of the N-Acetylneuraminate Lyase NanA. Front. Microbiol. 7:147. doi: 10.3389/fmicb.2016.00147

PubMed Abstract | CrossRef Full Text | Google Scholar

Diebel, L. N., Raafat, A. M., Dulchavsky, S. A., and Brown, W. J. (1996). Gallbladder and biliary tract candidiasis. Surgery 120, 760–765. doi: 10.1016/s0039-6060(96)80028-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Farid, S. G., Craven, R. A., Peng, J., Bonney, G. K., Perkins, D. N., Selby, P. J., et al. (2011). Shotgun proteomics of human bile in hilar cholangiocarcinoma. Proteomics 11, 2134–2138. doi: 10.1002/pmic.201000653

PubMed Abstract | CrossRef Full Text | Google Scholar

Farina, A., Dumonceau, J. M., Frossard, J. L., Hadengue, A., Hochstrasser, D. F., and Lescuyer, P. (2009). Proteomic analysis of human bile from malignant biliary stenosis induced by pancreatic cancer. J. Proteome Res. 8, 159–169. doi: 10.1021/pr8004925

PubMed Abstract | CrossRef Full Text | Google Scholar

Futatsuya, T., Ushigami, T., Nomura, F., Anzawa, K., Mochizuki, T., Cho, O., et al. (2020). Scalp microbiota in members of a Japanese high school judo team including Trichophyton tonsurans carriers. J. Dermatol. 47, 1020–1026. doi: 10.1111/1346-8138.15462

PubMed Abstract | CrossRef Full Text | Google Scholar

Geetha, A. (2002). Evidence for oxidative stress in the gall bladder mucosa of gall stone patients. J. Biochem. Mol. Biol. Biophys. 6, 427–432. doi: 10.1080/1025814021000036179

CrossRef Full Text | Google Scholar

Gonzalez, J. F., Tucker, L., Fitch, J., Wetzel, A., White, P., and Gunn, J. S. (2019). Human bile-mediated regulation of Salmonella Curli fimbriae. J. Bacteriol. 201:e00055-19. doi: 10.1128/JB.00055-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Griffen, A. L., Beall, C. J., Campbell, J. H., Firestone, N. D., Kumar, P. S., Yang, Z. K., et al. (2012). Distinct and complex bacterial profiles in human periodontitis and health revealed by 16S pyrosequencing. ISME J. 6, 1176–1185. doi: 10.1038/ismej.2011.191

PubMed Abstract | CrossRef Full Text | Google Scholar

Grigor'eva, I. N., and Romanova, T. I. (2020). Gallstone disease and microbiome. Microorganisms 8:835. doi: 10.3390/microorganisms8060835

PubMed Abstract | CrossRef Full Text | Google Scholar

Hancke, E., and Marklein, G. (1983). Experimental gallstone formation. Etiological significance of beta-glucuronidase producing bacteria and biliary obstruction. Langenbecks Arch. Chir. 359, 257–264. doi: 10.1007/BF01257312

CrossRef Full Text | Google Scholar

Harada, K., and Nakanuma, Y. (2010). Biliary innate immunity: function and modulation. Mediat. Inflamm. 2010, 1–9. doi: 10.1155/2010/373878

PubMed Abstract | CrossRef Full Text | Google Scholar

Hebrard, M., Viala, J. P., Meresse, S., Barras, F., and Aussel, L. (2009). Redundant hydrogen peroxide scavengers contribute to Salmonella virulence and oxidative stress resistance. J. Bacteriol. 191, 4605–4614. doi: 10.1128/JB.00144-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, H., Shao, W., Liu, Q., Liu, N., Wang, Q., Xu, J., et al. (2022). Gut microbiota promotes cholesterol gallstone formation by modulating bile acid composition and biliary cholesterol secretion. Nat. Commun. 13:252. doi: 10.1038/s41467-021-27758-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, J., Tang, J., Zhang, X., Yang, K., Zhong, A., Yang, Q., et al. (2023). Landscape in the gallbladder mycobiome and bacteriome of patients undergoing cholelithiasis with chronic cholecystitis. Front. Microbiol. 14:1131694. doi: 10.3389/fmicb.2023.1131694

PubMed Abstract | CrossRef Full Text | Google Scholar

Jain, S., Smyth, D., O'Hagan, B. M. G., Heap, J. T., McMullan, G., Minton, N. P., et al. (2017). Inactivation of the dnaK gene in Clostridium difficile 630 Deltaerm yields a temperature-sensitive phenotype and increases biofilm-forming ability. Sci. Rep. 7:17522. doi: 10.1038/s41598-017-17583-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Jean, D., Briolat, V., and Reysset, G. (2004). Oxidative stress response in Clostridium perfringens. Microbiology (Reading) 150, 1649–1659. doi: 10.1099/mic.0.27017-0

CrossRef Full Text | Google Scholar

Kida, Y., Shimizu, T., and Kuwano, K. (2011). Cooperation between LepA and PlcH contributes to the in vivo virulence and growth of Pseudomonas aeruginosa in mice. Infect. Immun. 79, 211–219. doi: 10.1128/IAI.01053-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Kose, S. H., Grice, K., Orsi, W. D., Ballal, M., and Coolen, M. J. L. (2018). Metagenomics of pigmented and cholesterol gallstones: the putative role of bacteria. Sci. Rep. 8:11218. doi: 10.1038/s41598-018-29571-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Kraal, L., Abubucker, S., Kota, K., Fischbach, M. A., and Mitreva, M. (2014). The prevalence of species and strains in the human microbiome: a resource for experimental efforts. PLoS One 9:e97279. doi: 10.1371/journal.pone.0097279

PubMed Abstract | CrossRef Full Text | Google Scholar

Kristiansen, T. Z., Bunkenborg, J., Gronborg, M., Molina, H., Thuluvath, P. J., Argani, P., et al. (2004). A proteomic analysis of human bile. Mol. Cell. Proteomics 3, 715–728. doi: 10.1074/mcp.M400015-MCP200

CrossRef Full Text | Google Scholar

Kunert, A., Losse, J., Gruszin, C., Huhn, M., Kaendler, K., Mikkat, S., et al. (2007). Immune evasion of the human pathogen Pseudomonas aeruginosa: elongation factor Tuf is a factor H and plasminogen binding protein. J. Immunol. 179, 2979–2988. doi: 10.4049/jimmunol.179.5.2979

PubMed Abstract | CrossRef Full Text | Google Scholar

Lammert, F., Gurusamy, K., Ko, C. W., Miquel, J. F., Mendez-Sanchez, N., Portincasa, P., et al. (2016). Gallstones. Nat. Rev. Dis. Primers. 2:16024. doi: 10.1038/nrdp.2016.24

CrossRef Full Text | Google Scholar

Lankisch, T. O., Metzger, J., Negm, A. A., Vosskuhl, K., Schiffer, E., Siwy, J., et al. (2011). Bile proteomic profiles differentiate cholangiocarcinoma from primary sclerosing cholangitis and choledocholithiasis. Hepatology 53, 875–884. doi: 10.1002/hep.24103

CrossRef Full Text | Google Scholar

Laohaviroj, M., Potriquet, J., Jia, X., Suttiprapa, S., Chamgramol, Y., Pairojkul, C., et al. (2017). A comparative proteomic analysis of bile for biomarkers of cholangiocarcinoma. Tumour Biol. 39:101042831770576. doi: 10.1177/1010428317705764

CrossRef Full Text | Google Scholar

Lee, J. H., Kim, M. W., Kim, B. S., Kim, S. M., Lee, B. C., Kim, T. S., et al. (2007). Identification and characterization of the Vibrio vulnificus rtxA essential for cytotoxicity in vitro and virulence in mice. J. Microbiol. 45, 146–152.

Google Scholar

Li, X., and Gao, P. (2019). Fatty liver increases gallstone disease risk in younger Chinese patients. Medicine (Baltimore) 98:e15940. doi: 10.1097/MD.0000000000015940

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Q., Ge, X., Xu, X., Zhong, Y., and Qie, Z. (2014). Comparison of the gene expression profiles between gallstones and gallbladder polyps. Int. J. Clin. Exp. Pathol. 7, 8016–8023.

PubMed Abstract | Google Scholar

Liu, P., Wood, D., and Nester, E. W. (2005). Phosphoenolpyruvate carboxykinase is an acid-induced, chromosomally encoded virulence factor in Agrobacterium tumefaciens. J. Bacteriol. 187, 6039–6045. doi: 10.1128/JB.187.17.6039-6045.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, J., Rincon, N., Wood, D. E., Breitwieser, F. P., Pockrandt, C., Langmead, B., et al. (2022). Metagenome analysis using the kraken software suite. Nat. Protoc. 17, 2815–2839. doi: 10.1038/s41596-022-00738-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Maeda, K., Okuda, Y., Enomoto, G., Watanabe, S., and Ikeuchi, M. (2021). Biosynthesis of a sulfated exopolysaccharide, synechan, and bloom formation in the model cyanobacterium Synechocystis sp. strain PCC 6803. elife 10:e66538. doi: 10.7554/eLife.66538

CrossRef Full Text | Google Scholar

Maurer, K. J., Rogers, A. B., Ge, Z., Wiese, A. J., Carey, M. C., and Fox, J. G. (2006). Helicobacter pylori and cholesterol gallstone formation in C57L/J mice: a prospective study. Am. J. Physiol. Gastrointest. Liver Physiol. 290, G175–G182. doi: 10.1152/ajpgi.00272.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Mokhtar, F. B. A., Plat, J., and Mensink, R. P. (2022). Genetic variation and intestinal cholesterol absorption in humans: a systematic review and a gene network analysis. Prog. Lipid Res. 86:101164. doi: 10.1016/j.plipres.2022.101164

CrossRef Full Text | Google Scholar

Molinero, N., Ruiz, L., Milani, C., Gutierrez-Diaz, I., Sanchez, B., Mangifesta, M., et al. (2019). The human gallbladder microbiome is related to the physiological state and the biliary metabolic profile. Microbiome 7:100. doi: 10.1186/s40168-019-0712-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Morales, C., Rojas, G., Rebolledo, C., Rojas-Herrera, M., Arias-Carrasco, R., Cuadros-Orellana, S., et al. (2022). Characterization of microbial communities from gut microbiota of hypercholesterolemic and control subjects. Front. Cell. Infect. Microbiol. 12:943609. doi: 10.3389/fcimb.2022.943609

PubMed Abstract | CrossRef Full Text | Google Scholar

Mugabi, R., Sandgren, D., Born, M., Leith, I., Horne, S. M., and Prubeta, B. M. (2012). The role of activated acetate intermediates in the control of Escherichia coli biofilm amounts. Webmedcentral 3:3577

Google Scholar

Navaneethan, U., Lourdusamy, V., Gk Venkatesh, P., Willard, B., Sanaka, M. R., and Parsi, M. A. (2015). Bile proteomics for differentiation of malignant from benign biliary strictures: a pilot study. Gastroenterol. Rep. (Oxf.) 3, 136–143. doi: 10.1093/gastro/gou066

PubMed Abstract | CrossRef Full Text | Google Scholar

Neznansky, A., Blus-Kadosh, I., Yerushalmi, G., Banin, E., and Opatowsky, Y. (2014). The Pseudomonas aeruginosa phosphate transport protein PstS plays a phosphate-independent role in biofilm formation. FASEB J. 28, 5223–5233. doi: 10.1096/fj.14-258293

CrossRef Full Text | Google Scholar

O'Gorman, C. M., Fuller, H., and Dyer, P. S. (2009). Discovery of a sexual cycle in the opportunistic fungal pathogen aspergillus fumigatus. Nature 457, 471–474. doi: 10.1038/nature07528

PubMed Abstract | CrossRef Full Text | Google Scholar

Oh, J., Byrd, A. L., Park, M., Program, N. C. S., Kong, H. H., and Segre, J. A. (2016). Temporal stability of the human skin microbiome. Cells 165, 854–866. doi: 10.1016/j.cell.2016.04.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Olaisen, M., Richard, M. L., Beisvag, V., Granlund, A. V. B., Royset, E. S., Rue, O., et al. (2022). The ileal fungal microbiota is altered in Crohn's disease and is associated with the disease course. Front. Med. (Lausanne) 9:868812. doi: 10.3389/fmed.2022.868812

PubMed Abstract | CrossRef Full Text | Google Scholar

Paterson, G. K., Cone, D. B., Northen, H., Peters, S. E., and Maskell, D. J. (2009). Deletion of the gene encoding the glycolytic enzyme triosephosphate isomerase (tpi) alters morphology of Salmonella enterica serovar typhimurium and decreases fitness in mice. FEMS Microbiol. Lett. 294, 45–51. doi: 10.1111/j.1574-6968.2009.01553.x

CrossRef Full Text | Google Scholar

Peng, Y., Yang, Y., Liu, Y., Nie, Y., Xu, P., Xia, B., et al. (2015). Cholesterol gallstones and bile host diverse bacterial communities with potential to promote the formation of gallstones. Microb. Pathog. 83-84, 57–63. doi: 10.1016/j.micpath.2015.05.002

CrossRef Full Text | Google Scholar

Ploszaj, T., Brauncajs, M., Traczyk-Borszynska, M., Matyjas, T., Pomorski, L., Wasiak, T., et al. (2021). The value of bacterial metagenomic analysis in post-surgical examination of gallstones. Arch. Microbiol. 203, 6323–6328. doi: 10.1007/s00203-021-02580-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Raudvere, U., Kolberg, L., Kuzmin, I., Arak, T., Adler, P., Peterson, H., et al. (2019). G:profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update). Nucleic Acids Res. 47, W191–W198. doi: 10.1093/nar/gkz369

PubMed Abstract | CrossRef Full Text | Google Scholar

Razizadeh, M. H., Khatami, A., and Zarei, M. (2021). Global status of Bufavirus, Cosavirus, and Saffold virus in gastroenteritis: a systematic review and Meta-analysis. Front. Med. (Lausanne) 8:775698. doi: 10.3389/fmed.2021.775698

PubMed Abstract | CrossRef Full Text | Google Scholar

Rebholz, C., Krawczyk, M., and Lammert, F. (2018). Genetics of gallstone disease. Eur. J. Clin. Investig. 48:e12935. doi: 10.1111/eci.12935

CrossRef Full Text | Google Scholar

Reinhard, L., Rupp, C., Riedel, H. D., Ruppert, T., Giese, T., Flechtenmacher, C., et al. (2012). S100A9 is a biliary protein marker of disease activity in primary sclerosing cholangitis. PLoS One 7:e29821. doi: 10.1371/journal.pone.0029821

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, H., Luo, M., Chen, J., Zhou, Y., Li, X., Zhan, Y., et al. (2019). Identification of TPD52 and DNAJB1 as two novel bile biomarkers for cholangiocarcinoma by iTRAQbased quantitative proteomics analysis. Oncol. Rep. 42, 2622–2634. doi: 10.3892/or.2019.7387

PubMed Abstract | CrossRef Full Text | Google Scholar

Rigoulay, C., Entenza, J. M., Halpern, D., Widmer, E., Moreillon, P., Poquet, I., et al. (2005). Comparative analysis of the roles of HtrA-like surface proteases in two virulent Staphylococcus aureus strains. Infect. Immun. 73, 563–572. doi: 10.1128/IAI.73.1.563-572.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosner, J. L., and Martin, R. G. (2009). An excretory function for the Escherichia coli outer membrane pore TolC: upregulation of marA and soxS transcription and rob activity due to metabolites accumulated in tolC mutants. J. Bacteriol. 191, 5283–5292. doi: 10.1128/JB.00507-09

CrossRef Full Text | Google Scholar

Rupp, C., Bode, K. A., Leopold, Y., Sauer, P., and Gotthardt, D. N. (2018). Pathological features of primary sclerosing cholangitis identified by bile proteomic analysis. Biochim. Biophys. Acta Mol. basis Dis. 1864, 1380–1389. doi: 10.1016/j.bbadis.2017.09.012

CrossRef Full Text | Google Scholar

Saadi, H., Pagnier, I., Colson, P., Cherif, J. K., Beji, M., Boughalmi, M., et al. (2013). First isolation of Mimivirus in a patient with pneumonia. Clin. Infect. Dis. 57, e127–e134. doi: 10.1093/cid/cit354

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakurai, S., Shinagawa, N., Fukui, T., and Yura, J. (1992). Bacterial adherence to human gallbladder epithelium. Surg. Today 22, 504–507. doi: 10.1007/BF00308894

PubMed Abstract | CrossRef Full Text | Google Scholar

Saltykova, I. V., Petrov, V. A., Logacheva, M. D., Ivanova, P. G., Merzlikin, N. V., Sazonov, A. E., et al. (2016). Biliary microbiota, gallstone disease and infection with Opisthorchis felineus. PLoS Negl. Trop. Dis. 10:e0004809. doi: 10.1371/journal.pntd.0004809

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanikidze, T., and Chikvaidze, E. (2016). Role of the free radicals in mechanisms of gallstone formation: an EPR study. Radiat. Prot. Dosim. 172, 317–324. doi: 10.1093/rpd/ncw237

PubMed Abstract | CrossRef Full Text | Google Scholar

Shah, A. H., Singh, A., Dhamgaye, S., Chauhan, N., Vandeputte, P., Suneetha, K. J., et al. (2014). Novel role of a family of major facilitator transporters in biofilm development and virulence of Candida albicans. Biochem. J. 460, 223–235. doi: 10.1042/BJ20140010

CrossRef Full Text | Google Scholar

Shao, Z., Shan, X., Jing, L., Wang, W., Li, W., Ren, Z., et al. (2023). Metagenome investigation of ocular microbiota of cataract patients with and without type 2 diabetes. Transl. Vis. Sci. Technol. 12:1. doi: 10.1167/tvst.12.6.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, R., Sachan, S. G., and Sharma, S. R. (2020). In vitro analysis of gallstone formation in the presence of bacteria. Indian J. Gastroenterol. 39, 473–480. doi: 10.1007/s12664-020-01055-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, N., Yadav, M., Tripathi, G., Mathew, B., Bindal, V., Falari, S., et al. (2022). Bile multi-omics analysis classifies lipid species and microbial peptides predictive of carcinoma of gallbladder. Hepatology 76, 920–935. doi: 10.1002/hep.32496

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, J., Wang, W., Wu, J., Feng, B., Chen, W., Wang, M., et al. (2012). Comparative proteomic profiling of human bile reveals SSP411 as a novel biomarker of cholangiocarcinoma. PLoS One 7:e47476. doi: 10.1371/journal.pone.0047476

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, H., Ye, F., Xie, L., Yang, J., Li, Z., Xu, P., et al. (2015). Metagenomic sequencing of bile from gallstone patients to identify different microbial community patterns and novel biliary bacteria. Sci. Rep. 5:17450. doi: 10.1038/srep17450

PubMed Abstract | CrossRef Full Text | Google Scholar

Shigemori, T., Hiasa, A., Inoue, Y., Oka, S., Yasuma, T., Nishiwaki, R., et al. (2022). Acute calculous cholecystitis caused by Streptococcus gallolyticus subspecies pasteurianus: a case report. Microorganisms 10:1929. doi: 10.3390/microorganisms10101929

PubMed Abstract | CrossRef Full Text | Google Scholar

Sikarwar, J., Singh, J., Singh, T. P., Sharma, P., and Sharma, S. (2022). The mechanism of action of Lactoferrin - nucleoside diphosphate kinase complex in combating biofilm formation. Protein Pept. Lett. 29, 839–850. doi: 10.2174/0929866529666220816160517

CrossRef Full Text | Google Scholar

Skar, V., Skar, A. G., Midtvedt, T., Lotveit, T., and Osnes, M. (1986). Beta-glucuronidase-producing bacteria in bile from the common bile duct in patients treated with endoscopic papillotomy for gallstone disease. Scand. J. Gastroenterol. 21, 253–256. doi: 10.3109/00365528609034656

PubMed Abstract | CrossRef Full Text | Google Scholar

Son, K. H., Ahn, C. B., Kim, H. J., and Kim, J. S. (2020). Quantitative proteomic analysis of bile in extrahepatic cholangiocarcinoma patients. J. Cancer 11, 4073–4080. doi: 10.7150/jca.40964

PubMed Abstract | CrossRef Full Text | Google Scholar

Stearns, J. C., Lynch, M. D., Senadheera, D. B., Tenenbaum, H. C., Goldberg, M. B., Cvitkovitch, D. G., et al. (2011). Bacterial biogeography of the human digestive tract. Sci. Rep. 1:170. doi: 10.1038/srep00170

PubMed Abstract | CrossRef Full Text | Google Scholar

Sung, J. Y., Leung, J. W., Olson, M. E., Lundberg, M. S., and Costerton, J. W. (1991). Demonstration of transient bacterobilia by foreign body implantation in feline biliary tract. Dig. Dis. Sci. 36, 943–948. doi: 10.1007/BF01297145

PubMed Abstract | CrossRef Full Text | Google Scholar

Szmyt, M., Kasprzak, A., Malkowski, W., Surdyk-Zasada, J., Przybyszewska, W., Siodla, E., et al. (2013). Tissue expression of S100 proteins in gallbladder mucosa of the patients with calculous cholecystitis. Folia Histochem. Cytobiol. 51, 141–148. doi: 10.5603/FHC.2013.0022

PubMed Abstract | CrossRef Full Text | Google Scholar

Taddei, T. H., Dziura, J., Chen, S., Yang, R., Hyogo, H., Sullards, C., et al. (2010). High incidence of cholesterol gallstone disease in type 1 Gaucher disease: characterizing the biliary phenotype of type 1 Gaucher disease. J. Inherit. Metab. Dis. 33, 291–300. doi: 10.1007/s10545-010-9070-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Tajeddin, E., Sherafat, S. J., Majidi, M. R., Alebouyeh, M., Alizadeh, A. H., and Zali, M. R. (2016). Association of diverse bacterial communities in human bile samples with biliary tract disorders: a survey using culture and polymerase chain reaction-denaturing gradient gel electrophoresis methods. Eur. J. Clin. Microbiol. Infect. Dis. 35, 1331–1339. doi: 10.1007/s10096-016-2669-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, W., Chen, R., Song, J., He, D., Wu, J., Chen, X., et al. (2022). Microbiota analysis with next-generation 16S rDNA gene sequencing in recurrent common bile duct stones. Ann. Transl. Med. 10:576. doi: 10.21037/atm-22-2247

PubMed Abstract | CrossRef Full Text | Google Scholar

Tedesco, D., Thapa, M., Chin, C. Y., Ge, Y., Gong, M., Li, J., et al. (2018). Alterations in intestinal microbiota Lead to production of interleukin 17 by intrahepatic gammadelta T-cell receptor-positive cells and pathogenesis of Cholestatic liver disease. Gastroenterology 154, 2178–2193. doi: 10.1053/j.gastro.2018.02.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Thursby, E., and Juge, N. (2017). Introduction to the human gut microbiota. Biochem. J. 474, 1823–1836. doi: 10.1042/BCJ20160510

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsai, M. H., Liang, Y. H., Chen, C. L., and Chiu, C. H. (2020). Characterization of Salmonella resistance to bile during biofilm formation. J. Microbiol. Immunol. Infect. 53, 518–524. doi: 10.1016/j.jmii.2019.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Tyanova, S., Temu, T., and Cox, J. (2016). The MaxQuant computational platform for mass spectrometry-based shotgun proteomics. Nat. Protoc. 11, 2301–2319. doi: 10.1038/nprot.2016.136

PubMed Abstract | CrossRef Full Text | Google Scholar

Voigtlander, T., Metzger, J., Schonemeier, B., Jager, M., Mischak, H., Manns, M. P., et al. (2017). A combined bile and urine proteomic test for cholangiocarcinoma diagnosis in patients with biliary strictures of unknown origin. United European Gastroenterol J 5, 668–676. doi: 10.1177/2050640616687836

PubMed Abstract | CrossRef Full Text | Google Scholar

von Beek, C., Waern, I., Eriksson, J., Melo, F. R., Robinson, C., Waller, A. S., et al. (2019). Streptococcal sagA activates a proinflammatory response in mast cells by a sublytic mechanism. Cell. Microbiol. 21:e13064. doi: 10.1111/cmi.13064

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Qi, M., Qin, C., and Hong, J. (2018). Role of the biliary microbiome in gallstone disease. Expert Rev. Gastroenterol. Hepatol. 12, 1193–1205. doi: 10.1080/17474124.2018.1533812

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Sun, W., Chang, Y., and Yi, Z. (2018). Differential proteomics analysis of bile between gangrenous cholecystitis and chronic cholecystitis. Med. Hypotheses 121, 131–136. doi: 10.1016/j.mehy.2018.07.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Yan, H. B., Zhang, Q., Liu, W. Y., Jiang, Y. H., Peng, G., et al. (2021). Enhancement of E-cadherin expression and processing and driving of cancer cell metastasis by ARID1A deficiency. Oncogene 40, 5468–5481. doi: 10.1038/s41388-021-01930-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Worthington, H. V., Hunt, L. P., McCloy, R. F., Ubbink, J. B., and Braganza, J. M. (2004). Dietary antioxidant lack, impaired hepatic glutathione reserve, and cholesterol gallstones. Clin. Chim. Acta 349, 157–165. doi: 10.1016/j.cccn.2004.06.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, T., Zhang, Z., Liu, B., Hou, D., Liang, Y., Zhang, J., et al. (2013). Gut microbiota dysbiosis and bacterial community assembly associated with cholesterol gallstones in large-scale study. BMC Genomics 14:669. doi: 10.1186/1471-2164-14-669

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, Y., Wang, D., Pan, X., Xia, B., Weng, Y., Long, Y., et al. (2020). TpiA is a key metabolic enzyme that affects virulence and resistance to aminoglycoside antibiotics through CrcZ in Pseudomonas aeruginosa. MBio 11:e02079-19. doi: 10.1128/mBio.02079-19

CrossRef Full Text | Google Scholar

Xiong, X., Rao, Y., Tu, X., Wang, Z., Gong, J., Yang, Y., et al. (2022). Gut archaea associated with bacteria colonization and succession during piglet weaning transitions. BMC Vet. Res. 18:243. doi: 10.1186/s12917-022-03330-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoo, K. S., Choi, H. S., Jun, D. W., Lee, H. L., Lee, O. Y., Yoon, B. C., et al. (2016). MUC expression in gallbladder epithelial tissues in cholesterol-associated gallbladder disease. Gut Liver 10, 851–858. doi: 10.5009/gnl15600

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, H., Xiong, J., Zhang, R., Hu, X., Qiu, J., Zhang, D., et al. (2016). Ndk, a novel host-responsive regulator, negatively regulates bacterial virulence through quorum sensing in Pseudomonas aeruginosa. Sci. Rep. 6:28684. doi: 10.1038/srep28684

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, D., Xiang, J., Wang, L., Xu, Z., Sun, L., Zhou, F., et al. (2013). Comparative proteomic analysis of gallbladder bile proteins related to cholesterol gallstones. PLoS One 8:e54489. doi: 10.1371/journal.pone.0054489

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, H., Xu, L., Dong, H., Hu, J., Gao, H., Yang, M., et al. (2015). Correlations between clinical features and mortality in patients with Vibrio vulnificus infection. PLoS One 10:e0136019. doi: 10.1371/journal.pone.0136019

CrossRef Full Text | Google Scholar

Zhou, H., Chen, B., Li, R. X., Sheng, Q. H., Li, S. J., Zhang, L., et al. (2005). Large-scale identification of human biliary proteins from a cholesterol stone patient using a proteomic approach. Rapid Commun. Mass Spectrom. 19, 3569–3578. doi: 10.1002/rcm.2207

CrossRef Full Text | Google Scholar