- Department of Radiology, Shengjing Hospital of China Medical University, Shenyang, China
Introduction: Pyogenic liver abscess (PLA) patients combined with diabetes mellitus (DM) tend to have more severe clinical manifestations than without DM. The mechanism responsible for this phenomenon is not entirely clear. The current study therefore aimed to comprehensively analyze the microbiome composition and metabolome in pus from PLA patients with and without DM, to determine the potential reasons for these differences.
Methods: Clinical data from 290 PLA patients were collected retrospectively. We analyzed the pus microbiota using 16S rDNA sequencing in 62 PLA patients. In addition, the pus metabolomes of 38 pus samples were characterized by untargeted metabolomics analysis. Correlation analyses of microbiota, metabolites and laboratory findings were performed to identify significant associations.
Results: PLA patients with DM had more severe clinical manifestations than PLA patients without DM. There were 17 discriminating genera between the two groups at the genus level, among which Klebsiella was the most discriminating taxa. The ABC transporters was the most significant differential metabolic pathway predicted by PICRUSt2. Untargeted metabolomics analysis showed that concentrations of various metabolites were significantly different between the two groups and seven metabolites were enriched in the ABC transporters pathway. Phosphoric acid, taurine, and orthophosphate in the ABC transporters pathway were negatively correlated with the relative abundance of Klebsiella and the blood glucose level.
Discussion: The results showed that the relative abundance of Klebsiella in the pus cavity of PLA patients with DM was higher than those without DM, accompanied by changes of various metabolites and metabolic pathways, which may be associated with more severe clinical manifestations.
Introduction
Pyogenic liver abscess (PLA) is a common intraperitoneal infection with several risk factors, such as diabetes mellitus (DM) and underlying hepatobiliary or pancreatic disease (Chadwick et al., 2018; Mukthinuthalapati et al., 2020). Although the mortality rate has improved, a significant percentage of PLA patients develop complications including extrahepatic migratory infection (EMI), abscess rupture and so on (Mavilia et al., 2016). DM is a predisposing factor for PLA (Thomsen et al., 2007; Mavilia et al., 2016); PLA patients with DM tend to have a higher prevalence of metastatic infection, bacteremia, multi-organ dysfunction syndrome and intensive care unit (ICU) admission (Foo et al., 2010; Tian et al., 2012). The mechanism underlying this phenomenon is not fully understood, so it is important to better understand the relevant mechanism, which can be used to develop preventive and treatment strategies.
Pathogenic microorganisms and their virulence levels are important factors affecting the prognosis of PLA (Liu et al., 2018). The gut microbiota play an important role in maintaining intestinal homeostasis and inducing the development of the intestinal immune system (Bouskra et al., 2008; Sonnenberg and Artis, 2012). Intestinal microorganisms entering the liver through the portal vein is one of the main causes of PLA (Lardière-Deguelte et al., 2015). Studies have reported a significant association between changes in the composition profile of gut microbiota and the development of DM (Gurung et al., 2020; Iatcu et al., 2021). In addition, studies have focused on the associations between metabolites and diabetes (Tang et al., 2019; Qi et al., 2022). Therefore, we speculated that the microbiome composition and metabolites in the pus cavity of PLA patients with DM were different from those of PLA patients without DM. A combination of 16SrDNA sequencing used for pathogen identification and untargeted metabolomics analyses for metabolite identification provided technical support to verify our hypothesis (Moon et al., 2021; Yu et al., 2021).
Materials and methods
Study population
The Ethics Committee of Shengjing Hospital of China Medical University approved this study (2022PS1067K) and written informed consent was obtained from all patients involved in the study for the publication of any potentially identifiable data included in this article. We retrieved the records of patients with a diagnosis of PLA from January 2017 to December 2021 from the electronic medical database at our institution (n = 290). According to the presence or absence of DM, these patients were divided into DM group (n = 189) and Non-DM group (n = 101).
Demographic information, laboratory findings, underlying diseases, hospital length of stay (LOS), the incidence of intensive care unit (ICU) admission and extrahepatic migratory infection (EMI) were obtained from medical records. Diagnostic criteria for PLA, DM, EMI, biliary tract disease and digestive system cancer were as previously described (Chang et al., 2015; Ren et al., 2020; Wang et al., 2022).
Pus sample collection
From January 2020 to December 2021, a total of 62 pus samples were prospective and continuously collected for conventional culture and 16S rDNA sequencing. According to the presence or absence of DM, these samples were divided into DM group (n = 38) and Non-DM group (n = 24). Of these, 38 pus samples (DM group, n = 22; Non-DM group, n = 16) were collected for untargeted metabolomics analysis. All collected pus samples were analyzed and no samples were artificially removed. The workflow of this study is shown in Figure 1.
The inclusion criteria were: (a) presence of a focal lesion or lesions in the liver on contrast enhanced CT images; (b) frank pus aspirated from the abscess cavity through diagnostic, and/or surgical drainage procedures; and (c) positive microbiological culture results from liver abscess and/or blood cultures.
The methods used for abscess drainage and pus collection were as described in a previous study (Song et al., 2014). All percutaneous procedures which strictly follows the principle of asepsis were performed under ultrasound guidance. Eighteen gage P.T.C Chiba needles of varying lengths were used to puncture the abscesses. The aspirated isolates including cultures were sent for microbiological analysis, 16S rDNA gene sequencing and untargeted metabolomics analysis. Under fluoroscopic guidance, 2–4 mL of undiluted contrast media was instilled into the abscess cavity through the 18G needle and then a 0.035-inch wire was inserted into the abscess cavity. After serial dilatation, an 8.5F pigtail catheter was inserted into the abscess. After collection, pus samples were stored at −80°C in a biological sample bank until used for 16S rDNA sequencing and untargeted metabolomics analysis.
16S rDNA sequencing
The full details regarding microbiome methodology are provided in the Supporting Information. The 16S rDNA sequencing was performed to characterize the pus microbiome composition. Length-heterogeneity PCR fingerprinting was routinely used to rapidly survey our samples and standardize the community amplification. We then characterized the microbial taxa associated with the pus microbiome using Multitag Sequencing of samples. Statistical data analysis was performed based on the feature table and feature sequence obtained as described in the Supporting Information. The full analysis script is available in the Supporting Information. Alpha and beta diversity analyses were calculated using QIIME2 and all the graphs were constructed using the R package (R Foundation for Statistical Computing, Vienna, Austria). PICRUSt was used for inferred metagenomic analyses as previously described (Puri et al., 2018). Significant difference analysis was performed to compare the differential species at the genus level between the two groups (p < 0.05). Enrichment analysis of abundant taxa or functions between the two groups was performed using the linear discriminant analysis effect size (LEfSe) with a threshold logarithmic linear discriminant analysis (LDA) score set at 3.0. BugBase was used to predict the bacterial phenotype (Meyer et al., 2008; Liao et al., 2022).
Untargeted metabolomics analysis
The full details regarding untargeted metabolomics analysis are provided in the Supporting Information. Here, we briefly describe the sample analysis. Partial Least-Squares Discriminant Analysis (PLS-DA) was conducted using metaX to discriminate the different variables between groups and the variable importance for the projection (VIP) value was calculated. A VIP cut-off value of 1.0 was used to select important features. A volcano plot was used to compare differential metabolites between the two groups. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed on significantly different metabolites (satisfying ratio ≥ 2 or ratio ≤ 0.5; a value of p ≤ 0.05; VIP > 1). All graphs were constructed using the R package (R Foundation for Statistical Computing).
Statistical analyses
Continuous variables were expressed as the mean ± standard deviation or median (interquartile range). Categorical variables were expressed as absolute numbers with frequencies (%). Differences between continuous data were analyzed using Student’s t-test or the Mann–Whitney U-test. Differences between categorical data were analyzed using the chi-squared test or Fisher’s exact test, as appropriate. All statistical analyses were performed using SPSS statistical software for Windows, version 25.0 (IBM, Armonk, NY, United States). Two-tailed p-values < 0.05 were considered statistically significant. Correlation heatmaps with signs were performed using the OmicStudio tools. Scatter plots were constructed using Prism 8 software (GraphPad, San Diego, CA, United States).
Results
Clinical characteristics of participants
A total of 290 PLA patients (DM, n = 189; non-DM, n = 101) were studied (Table 1). As expected, there are a number of laboratory findings of PLA patients in both DM and non-DM groups falling outside the reference range. For example, leukocyte, neutrophils (NEUT),NEUT% and so on. Compared with the non-DM group, PLA patients with DM had a higher C-reactive protein (CRP, an inflammatory marker)level, higher EMI incidence and longer LOS, which showed that PLA patients with DM have more severe clinical manifestations.
Pus microbiota analysis
To determine the effects of DM on the pus microbiome composition in PLA patients, pus samples of 62 consecutive PLA patients (DM, n = 38; non-DM, n = 24) were performed using 16S rDNA sequencing. The clinical data of these 62 patients are shown in Table 2. As expected, compared with non-DM patients, PLA patients with DM had a higher CRP level. Although there was a statistical difference in platelets (PLT) in Table 2, the PLT in both groups were all within the reference range and had no clinical significance. Besides, there was no statistical difference in PLT in Table 1. Therefore, we speculate that the statistical difference in PLT in Table 2 may be due to the small sample size.
Supplementary Table S1 shows that 16SrDNA sequencing was able to detect the microbiome composition in conventionally cultured negative samples and was able to detect more bacteria in samples that were conventionally cultured as a single strain.
Alpha and beta diversities of the pus microbiota
Alpha diversity represents the species evenness and richness within the microbiota, while beta diversity represents the shared diversity within the microbiota at different ecological distances (Yang et al., 2022). The Chao 1 index describing species richness revealed no significant difference (DM vs. non-DM; p = 0.12) (Figure 2A). In addition, the Shannon diversity index did not show significant difference between the DM vs. non-DM patients; p = 0.89 (Figure 2B). We then evaluated the beta diversity of the two groups. Using Bray Curtis PCoA, we found some differences in species classifications between DM and non-DM patients (Figure 2C), but there was no statistical difference. We then conducted a correlation analysis between Chao1, Simpson, Shannon and Pielou-e indices of α diversity and clinical data, which suggested (Supplementary Table S2) that the Simpson and Pielou-e indices were negatively correlated with glycosylated hemoglobin HbA1c levels (Figures 2D,E).
Figure 2. Alpha diversity, beta diversity, and taxonomy community of 16SrDNA sequencing and linear correlation analysis. (A) Chao 1 index of the diabetes mellitus (DM) and non-DM groups. (B) Shannon diversity index of the DM and non-DM groups. (C) Beta diversity analyzed by Bray Curtis PCoA. (D) Scatter plot of the correlation between Simpson diversity index and HbA1c. (E) Scatter plot of the correlation between the Pielou-e diversity index and HbA1c. (F) Taxonomic percentages according to compositions at the genus level.
Alterations in pus microbiome composition due to DM
To determine the effects of DM, we first compared the microbial taxa in DM vs non-DM patients. Compared with non-DM patients, at the genus level, DM patients showed significantly greater relative abundance of Klebsiella (Figure 2F). Additionally, DM patients showed lesser levels of the genera, Acidovorax, Enterococcus, Bacteroides, Citrobacter and Streptococcus, when compared with the non-DM group.
Microbiome signatures for PLA patients with DM were distinct from non-DM patients
The statistical results of LEfSe included two parts, namely, a histogram of the LDA value distribution and a cladogram. The histogram of LDA values (Figure 3A) showed that there were 39 differentially abundant taxa at different taxonomic levels, of which 24 were from the DM group and 15 were from the non-DM group. The genus Klebsiella had the largest LDA score. The cladogram (Figure 3B) showed that the relative abundance of genus Klebsiella was greater in DM patients, when compared with non-DM patients.
Figure 3. LEfSe analysis, microbial difference analysis, and heatmap of correlation analysis. (A) Histogram of the logarithmic linear discriminant analysis (LD) value distribution (LDA score > 3.0; p < 0.05). A longer bar size corresponds to a larger enrichment of the group. (B) Cladogram. (C) Analysis results of microbial differences at the genus level in the DM and non-DM groups. (D) Heatmap of correlation analysis between the relative abundance of the top 10 microbiota at the genus level and the laboratory findings.
Taxa relative abundance in the DM and non-DM groups for the discriminating genera is shown in Figure 3C. There were 17 discriminating genera between the two groups at the genus level, among which Klebsiella was the most discriminating taxa.
Correlation analysis between pus microbiome composition and laboratory findings
To determine the functional relationships of microbiome composition and laboratory findings, we performed correlation analysis between the relative abundance of the top 10 microbiota at the genus level and using laboratory findings (Figure 3D). The results showed that the relative abundance of Klebsiella was positively correlated with blood glucose and CRP levels, which suggested that the relative abundance of Klebsiella was related to a high glucose environment and severe inflammatory reaction. The relative abundance of Acidovorax and Faecalibacterium was negatively correlated with CRP levels, but there was no correlation with blood glucose levels.
Functional prediction using BugBase and PICRUSt2
As shown by the potential predictions of phenotypic functions in the bacterial community in PLA, nine potential microbial phenotypes were detected. We observed a significantly higher representation of those containing mobile elements bacteria, facultatively anaerobic bacteria, potential pathogenic bacteria, bacteria related to biofilm formation and stress-tolerant bacteria in the DM group, when compared with the non-DM group (Table 3).
PICRUSt2 was used for functional profiling of the microbial community by using the KEGG database. Finally, we identified 38 KEGG pathways at level 3, which were significantly changed in percentage means between the two groups (p < 0.01). The top 30 pathway names sorted by p-values are included in Figure 4, with the most significant differential metabolic pathway involving ABC transporters.
Figure 4. Prediction of microbial function in diabetic mellitus (DM) and non-DM groups based on the Kyoto Encyclopedia of Genes and Genomes database. The ABC transporters pathway was the most significant differential metabolic pathway in the top 30 pathways shown in the figure.
Pus metabolic differences between PLA patients with and without DM
To further identify metabolomic changes in PLA patients with DM, we performed a non-targeted metabolomics study of 38 pus samples (DM, n = 22; non-DM, n = 16). PLS-DA, which has good discrimination ability, was used to analyze the metabolic profiles based on class information (Yu et al., 2021). The PLS-DA model showed that there was separation between the two groups (Figure 5A).
Figure 5. Partial least squares–discriminate analysis PLS-DA model, volcano plot, and heatmap of correlation analysis. (A) Score plot of the PLS-DA model in the diabetic mellitus (DM) and non-DM groups. (B) Volcano plot of differential metabolites between the two groups. The horizontal coordinate represents the fold change (log2FoldChange) of metabolites in different groups, and the vertical coordinate represents the significance level of the difference (−log10 p value). Each point in the volcano plot represents a metabolite, in which the metabolites with significantly up-regulated expressions are represented by red dots, and those with significantly down-regulated expressions are represented by green dots. The dot size variable importance for the projection values are represented by dot size. (C) Volcano plot of differential metabolites in the ABC transporters pathway between the two groups. (D) Heatmap of correlation analysis between the top 10 relative abundance microbiota at the genus level and the seven differentially accumulated metabolites in the ABC transporters pathway. (E) Heatmap of correlation analysis between the laboratory findings and the seven differentially accumulated metabolites in the ABC transporters pathway.
A volcano plot (Figure 5B) was constructed to show trends in differentially abundant metabolites between the two groups. Using KEGG enrichment analysis (Supplementary Figure S1), seven main pathways (Cellular Processes, Drug Development, Environmental Information Processing, Genetic Information Processing, Human Diseases Metabolism and Organismal Systems) were identified, with the most abundant being Metabolism. In this pathway, we found that the differential metabolites enriched in the DM group were mainly involved in amino acid biosynthesis and metabolism, lipid metabolism, glucose metabolism and polycyclic aromatic hydrocarbon degradation.
In the environmental information processing pathway, we found the ABC transporters pathway showed differences between the two groups, which was consistent with the 16SrDNA prediction. Based on variable importance in the projection values >1 in the loading plot, FC ≥ 2 or FC ≤ 0.5, and p < 0.05 with KEGG annotations, seven differentially accumulated metabolites were identified in the ABC transporters pathway, of which four metabolites (2′-deoxyuridine, 2′-deoxyinosine, 3-dehydroxycarnitine, and inosine) were enriched and three metabolites (phosphoric acid, taurine, and orthophosphate) were depleted (Supplementary Table S3; Figure 5C).
Correlation analysis between differential metabolites in the ABC transporters pathway and pus microbiota composition
To identify functional relationships of pus microbiota and differentially accumulated metabolites in the ABC transporters pathway, we performed correlation analysis based on Pearson’s correlation coefficients. The top 10 relative abundance microbiota at the genus level and the above seven differentially-accumulated metabolites in the ABC transporters pathway were included for analysis, which showed that the metabolites was correlated with pus microbiota. The heatmap of the correlation is shown in Figure 5D. The results showed that phosphoric acid, taurine, and orthophosphate were positively correlated with Acidovorax, Caulobacter, Citrobacter, and Pseudomonas, while they were negatively correlated with Klebsiella. The 2′-deoxyuridine, 2′-deoxyinosine, and 3-dehydroxycarnitine were negatively correlated with Acidovorax and Caulobacter.
Correlation analysis between the differential metabolites in the ABC transporters pathway and laboratory findings
Figure 5E shows the correlation heatmap for the seven differential metabolites in the ABC transporters pathway and the laboratory findings. The results suggested that all these metabolites were closely related to fasting blood glucose (FBG). The 2′-deoxyuridine, 2′-deoxyinosine, 3-dehydroxycarnitine and inosine were positively related to FBG, while phosphoric acid, taurine and orthophosphate were negatively correlated with FBG. These metabolites were also related to albumin (ALB), alanine aminotransferase (ALT) and albumin-bilirubin (ALBI), which are closely correlated with liver function.
Discussion
The current study showed that PLA patients with DM had a higher CRP level, a higher EMI incidence and longer LOS, which is consistent with previous studies (Foo et al., 2010; Tian et al., 2012). Diabetes-induced functional defects in neutrophil responses to pathogenic microorganism, including impairments in phagocytosis, bacterial killing, neutrophil migration, cytokine production, apoptosis and neutrophil extracellular trap (NET) formation may contribute to the incidence of EMI in PLA patients, when combined with DM (Hodgson et al., 2015; Tay et al., 2018; Njeim et al., 2020). However, Lee et al. (2018) reported that type 2 DM did not overtly affect neutrophil intra- and extra-cellular killing of hypervirulent Klebsiella pneumoniae. Thus, the mechanism of severity of PLA patients with diabetes still needs to be fully determined.
Host immune defenses and pathogenic microorganisms affect the prognosis and the development of PLA. Identifying the differential microbiome composition differences between PLA patients with and without DM is therefore helpful in identifying potential causes. To the best of our knowledge, the current study is the first to use 16S rDNA sequencing to identify differences of microbiota composition between PLA patients with and without DM. Although the alpha diversity did not distinguish between the two groups, the results showed that the Simpson index linked with species richness and the Pielou-e index linked with species evenness were negatively correlated with the HbA1c levels, suggesting that a high glucose level may inhibit microbial diversity in the pus cavity. The influence of DM on intestinal flora abundance is still controversial. Some studies reported lower alpha diversity indices (Shannon and Chao) of gut microbiota in diabetic patients, when compared with healthy controls (Nuli et al., 2019; Zhao et al., 2019), while other studies reported increased microbial richness of gut microbiota in diabetic patients (Bai et al., 2022). The effect of DM on intestinal flora and microorganisms in pus may also differ.
Klebsiella, whose relative abundance was higher in the DM group, was the most discriminating genus between the two groups and the relative abundance of Klebsiella was positively correlated with blood glucose and CRP levels. Therefore, we speculate that the severe inflammation in the DM group was partially due to the higher relative abundance of Klebsiella. Previous studies reported that Klebsiella pneumonia liver abscess tended to have a higher incidence of EMI, especially in DM patients (Yoon et al., 2014; Xu et al., 2020). Lee et al. (2018) postulated that the mechanism may be that Klebsiella pneumonia might activate platelet activity especially in DM patients with poor glycemic control and form a septic thrombosis, which can lead to seeding of microorganisms with metastatic infections.
Based on the prediction of BugBase, we observed a significantly higher representation of those containing mobile elements bacteria, facultatively anaerobic bacteria, potential pathogenic bacteria, bacteria related to biofilm formation and stress-tolerant bacteria in the DM group, when compared with the non-DM group. Biofilm formation plays an important role in antibiotic resistance and escape of microbes from the body’s immune surveillance system (Rather et al., 2021). Stress tolerance is a characteristic that improves resistance toward adverse environmental abiotic and biotic stress factors (Vassilev et al., 2012). Differences in these bacterial phenotypes may be one of the reasons that PLA patients with DM have more serious clinical manifestations.
Another noteworthy finding was that the ABC transporters pathway was the most significant differential metabolic pathway in KEGG pathways predicted by PICRUSt2. The ABC transporters superfamily is among the largest and most broadly expressed protein superfamily, is found in all living organisms and is involved in almost every cellular, biological and physiological system (Mercado-Lubo and McCormick, 2010). ABC transporters in the accessory genomes of bacterial pathogens significantly influence both virulence and antimicrobial resistance (Theodoulou and Kerr, 2015; Farzand et al., 2021). There was a large number of differential metabolites between the two groups, of which seven belonged to the ABC transporters pathway. All of these seven metabolites were closely related to FBG. Some of them were correlated with liver function and the relative abundance of microbiota, suggesting that these metabolites may contribute to the more severe clinical manifestations in the DM group. Unfortunately, there are few reports on the functions of these metabolites, which need further study in the future.
There were several limitations to this study. First, the current study is an exploratory study, so we did not perform a power analysis to estimate the sample size. Second, our study only detected the metabolic situation in pus samples without untargeted metabolomics analysis on serum samples, which can better reflect the metabolic situation of the body. In addition, it was difficult to determine whether the differential metabolites between the two groups detected in pus samples came from the host or microorganisms in pus. It is therefore necessary to collect both pus samples and blood samples in future research.
In conclusion, PLA patients with DM had more severe clinical manifestations than PLA patients without DM. The relative abundance of Klebsiella in the pus cavity of PLA patients with DM was higher than in PLA patients without DM and was accompanied by changes of various metabolites and metabolic pathways, which may be associated with the more severe clinical manifestations.
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.
Ethics statement
The studies involving human participants were reviewed and approved by The Ethics Committee of Shengjing Hospital of China Medical University. The patients/participants provided their written informed consent to participate in this study.
Author contributions
ZC and YG designed the study and analyzed and interpreted the data. YG and HW collected data. YG, ZL, and ZC wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
The present study was supported by the National Natural Science Foundation of China (81901856 and 82272097) and the 345 Talent Project in Shengjing Hospital of China Medical 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.1211835/full#supplementary-material
References
Bai, J., Wan, Z., Zhang, Y., Wang, T., Xue, Y., and Peng, Q. (2022). Composition and diversity of gut microbiota in diabetic retinopathy. Front. Microbiol. 13:926926. doi: 10.3389/fmicb.2022.926926
Bouskra, D., Brézillon, C., Bérard, M., Werts, C., Varona, R., Boneca, I. G., et al. (2008). Lymphoid tissue genesis induced by commensals through NOD1 regulates intestinal homeostasis. Nature 456, 507–510. doi: 10.1038/nature07450
Chadwick, M., Shamban, L., and Neumann, M. (2018). Pyogenic liver abscess with no predisposing risk factors. Case Rep. Gastrointest. Med. 2018, 9509356–9509354. doi: 10.1155/2018/9509356
Chang, Z., Zheng, J., Ma, Y., and Liu, Z. (2015). Analysis of clinical and CT characteristics of patients with Klebsiella pneumoniae liver abscesses: an insight into risk factors of metastatic infection. Int. J. Infect. Dis. 33, 50–54. doi: 10.1016/j.ijid.2014.12.041
Farzand, R., Rajakumar, K., Barer, M. R., Freestone, P. P. E., Mukamolova, G. V., Oggioni, M. R., et al. (2021). A virulence associated Siderophore importer reduces antimicrobial susceptibility of Klebsiella pneumoniae. Front. Microbiol. 12:607512. doi: 10.3389/fmicb.2021.607512
Foo, N. P., Chen, K. T., Lin, H. J., and Guo, H. R. (2010). Characteristics of pyogenic liver abscess patients with and without diabetes mellitus. Am. J. Gastroenterol. 105, 328–335. doi: 10.1038/ajg.2009.586
Gurung, M., Li, Z., You, H., Rodrigues, R., Jump, D. B., Morgun, A., et al. (2020). Role of gut microbiota in type 2 diabetes pathophysiology. EBioMedicine 51:102590. doi: 10.1016/j.ebiom.2019.11.051
Hodgson, K., Morris, J., Bridson, T., Govan, B., Rush, C., and Ketheesan, N. (2015). Immunological mechanisms contributing to the double burden of diabetes and intracellular bacterial infections. Immunology 144, 171–185. doi: 10.1111/imm.12394
Iatcu, C. O., Steen, A., and Covasa, M. (2021). Gut microbiota and complications of Type-2 diabetes. Nutrients 14:166. doi: 10.3390/nu14010166
Lardière-Deguelte, S., Ragot, E., Amroun, K., Piardi, T., Dokmak, S., Bruno, O., et al. (2015). Hepatic abscess: diagnosis and management. J. Visc. Surg. 152, 231–243. doi: 10.1016/j.jviscsurg.2015.01.013
Lee, C. H., Chuah, S. K., Tai, W. C., and Chen, I. L. (2018). Platelet reactivity in diabetic patients with invasive Klebsiella pneumoniae liver abscess syndrome. Infect. Drug Resist. 11, 1669–1676. doi: 10.2147/idr.S174913
Liao, Y., Zeng, X., Xie, X., Liang, D., Qiao, H., Wang, W., et al. (2022). Bacterial signatures of cerebral thrombi in large vessel occlusion stroke. MBio 13:e0108522. doi: 10.1128/mbio.01085-22
Liu, C., Shi, J., and Guo, J. (2018). High prevalence of hypervirulent Klebsiella pneumoniae infection in the genetic background of elderly patients in two teaching hospitals in China. Infect. Drug Resist. 11, 1031–1041. doi: 10.2147/idr.S161075
Mavilia, M. G., Molina, M., and Wu, G. Y. (2016). The evolving nature of hepatic abscess: a review. J. Clin. Transl. Hepatol. 4, 158–168. doi: 10.14218/JCTH.2016.00004
Mercado-Lubo, R., and McCormick, B. A. (2010). The interaction of gut microbes with host ABC transporters. Gut Microbes 1, 301–306. doi: 10.4161/gmic.1.5.12925
Meyer, F., Paarmann, D., D'Souza, M., Olson, R., Glass, E. M., Kubal, M., et al. (2008). The metagenomics RAST server - a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformat. 9:386. doi: 10.1186/1471-2105-9-386
Moon, J., Kim, N., Lee, H. S., Lee, S. T., Jung, K. H., Park, K. I., et al. (2021). Nanopore 16S amplicon sequencing enhances the understanding of pathogens in medically intractable diabetic foot infections. Diabetes 70, 1357–1371. doi: 10.2337/db20-0907
Mukthinuthalapati, V., Attar, B. M., Parra-Rodriguez, L., Cabrera, N. L., Araujo, T., and Gandhi, S. (2020). Risk factors, management, and outcomes of pyogenic liver abscess in a US safety net hospital. Dig. Dis. Sci. 65, 1529–1538. doi: 10.1007/s10620-019-05851-9
Njeim, R., Azar, W. S., Fares, A. H., Azar, S. T., Kfoury Kassouf, H., and Eid, A. A. (2020). NETosis contributes to the pathogenesis of diabetes and its complications. J. Mol. Endocrinol. 65, R65–r76. doi: 10.1530/jme-20-0128
Nuli, R., Azhati, J., Cai, J., Kadeer, A., Zhang, B., and Mohemaiti, P. (2019). Metagenomics and Faecal metabolomics integrative analysis towards the impaired glucose regulation and type 2 diabetes in Uyghur-related Omics. J. Diabetes Res. 2019, 2893041–2893015. doi: 10.1155/2019/2893041
Puri, P., Liangpunsakul, S., Christensen, J. E., Shah, V. H., Kamath, P. S., Gores, G. J., et al. (2018). The circulating microbiome signature and inferred functional metagenomics in alcoholic hepatitis. Hepatology 67, 1284–1302. doi: 10.1002/hep.29623
Qi, Q., Li, J., Yu, B., Moon, J. Y., Chai, J. C., Merino, J., et al. (2022). Host and gut microbial tryptophan metabolism and type 2 diabetes: an integrative analysis of host genetics, diet, gut microbiome and circulating metabolites in cohort studies. Gut 71, 1095–1105. doi: 10.1136/gutjnl-2021-324053
Rather, M. A., Gupta, K., and Mandal, M. (2021). Microbial biofilm: formation, architecture, antibiotic resistance, and control strategies. Braz. J. Microbiol. 52, 1701–1718. doi: 10.1007/s42770-021-00624-x
Ren, Y., Wang, H., Chang, Z., and Liu, Z. (2020). Clinical and computed tomography features of extended-spectrum β-lactamase-producing Klebsiella pneumoniae liver abscess. BMC Infect. Dis. 20:416. doi: 10.1186/s12879-020-05142-z
Song, Y. G., Shim, S. G., Kim, K. M., Lee, D. H., Kim, D. S., Choi, S. H., et al. (2014). Profiling of the bacteria responsible for pyogenic liver abscess by 16S rRNA gene pyrosequencing. J. Microbiol. 52, 504–509. doi: 10.1007/s12275-014-4241-7
Sonnenberg, G. F., and Artis, D. (2012). Innate lymphoid cell interactions with microbiota: implications for intestinal health and disease. Immunity 37, 601–610. doi: 10.1016/j.immuni.2012.10.003
Tang, X., You, J., Liu, D., Xia, M., He, L., and Liu, H. (2019). 5-Hydroxyhexanoic acid predicts early renal functional decline in type 2 diabetes patients with microalbuminuria. Kidney Blood Press. Res. 44, 245–263. doi: 10.1159/000498962
Tay, H. M., Dalan, R., Li, K. H. H., Boehm, B. O., and Hou, H. W. (2018). A novel microdevice for rapid neutrophil purification and Phenotyping in type 2 diabetes mellitus. Small 14:1702832. doi: 10.1002/smll.201702832
Theodoulou, F. L., and Kerr, I. D. (2015). ABC transporter research: going strong 40 years on. Biochem. Soc. Trans. 43, 1033–1040. doi: 10.1042/bst20150139
Thomsen, R. W., Jepsen, P., and Sørensen, H. T. (2007). Diabetes mellitus and pyogenic liver abscess: risk and prognosis. Clin. Infect. Dis. 44, 1194–1201. doi: 10.1086/513201
Tian, L. T., Yao, K., Zhang, X. Y., Zhang, Z. D., Liang, Y. J., Yin, D. L., et al. (2012). Liver abscesses in adult patients with and without diabetes mellitus: an analysis of the clinical characteristics, features of the causative pathogens, outcomes and predictors of fatality: a report based on a large population, retrospective study in China. Clin. Microbiol. Infect. 18, E314–E330. doi: 10.1111/j.1469-0691.2012.03912.x
Vassilev, N., Eichler-Löbermann, B., and Vassileva, M. (2012). Stress-tolerant P-solubilizing microorganisms. Appl. Microbiol. Biotechnol. 95, 851–859. doi: 10.1007/s00253-012-4224-8
Wang, H., Guo, Y., Yan, B., Zhang, Q., Pan, T., Liu, Z., et al. (2022). Development and validation of a prediction model based on clinical and CT features for invasiveness of K. pneumoniae liver abscess. Eur. Radiol. 32, 6397–6406. doi: 10.1007/s00330-022-08740-4
Xu, S. B., Yang, S. Y., Weng, S. S., Chen, C., Chen, S., Zhang, W. H., et al. (2020). Clinical characteristics of patients with Klebsiella pneumoniae pyogenic liver abscess. Zhonghua Nei Ke Za Zhi 59, 439–444. doi: 10.3760/cma.j.cn112138-20190610-00403
Yang, L., Xiang, Z., Zou, J., Zhang, Y., Ni, Y., and Yang, J. (2022). Comprehensive analysis of the relationships between the gut microbiota and fecal Metabolome in individuals with primary Sjogren's syndrome by 16S rRNA sequencing and LC-MS-based metabolomics. Front. Immunol. 13:874021. doi: 10.3389/fimmu.2022.874021
Yoon, J. H., Kim, Y. J., Jun, Y. H., Kim, S. I., Kang, J. Y., Suk, K. T., et al. (2014). Liver abscess due to Klebsiella pneumoniae: risk factors for metastatic infection. Scand. J. Infect. Dis. 46, 21–26. doi: 10.3109/00365548.2013.851414
Yu, L., Lai, Q., Feng, Q., Li, Y., Feng, J., and Xu, B. (2021). Serum metabolic profiling analysis of chronic gastritis and gastric Cancer by untargeted metabolomics. Front. Oncol. 11:636917. doi: 10.3389/fonc.2021.636917
Keywords: bacterial composition, pyogenic liver abscess, 16S rDNA sequencing, untargeted metabolomics analysis, metabolome
Citation: Guo Y, Wang H, Liu Z and Chang Z (2023) Comprehensive analysis of the microbiome and metabolome in pus from pyogenic liver abscess patients with and without diabetes mellitus. Front. Microbiol. 14:1211835. doi: 10.3389/fmicb.2023.1211835
Edited by:
George Grant, University of Aberdeen, United KingdomReviewed by:
Scarlet Louis-Jean, Anne Arundel Medical Center, United StatesVishal G. Shelat, Tan Tock Seng Hospital, Singapore
Whitney Crossland, Texas Tech University, United States
Copyright © 2023 Guo, Wang, Liu and Chang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhihui Chang, changzh@sj-hospital.org