- 1Department of Urology, The Second Affiliated Hospital of Soochow University, Suzhou, China
- 2Department of Urology, Changshu Hospital Affiliated to Soochow University, Changshu, China
- 3Department of Rheumatology, Changshu Hospital Affiliated to Soochow University, Changshu, China
- 4Department of Nephrology, Taizhou People’s Hospital, Taizhou, China
Previous studies suggest that patients with nephrolithiasis exhibit dysbiosis in their gut microbiota, but those studies were conducted in calcium oxalate stone patients. We aimed to explore the association of gut microbiota and biochemical features of renal uric acid stone (UAS) patients in a Chinese population and identify the related bacteria that may affect the pathopoiesis of UAS. A case-control study of 117 patients with UAS, 123 patients with gout, and 135 healthy controls were included from January 2014 to October 2020. For each subject, data on demographics, biochemical parameters of blood and urine were analyzed. Fifteen patients with gout, 16 patients with UAS, 17 UAS patients with gout, and 17 healthy subjects were enrolled and provided fecal samples. The characteristics of gut microbiota were explored by using 16S ribosomal RNA (rRNA) gene sequencing and analyzed by using a combination of software mother and R. Hyperuricemia was the main risk factor for the development of gout and UAS. Obesity, dyslipidemia, and aciduria were unique risk factors for UAS patients. The richness, diversity, and relative abundance of dominant bacteria at the phylum and genus levels of gut microbiota in UAS patients were significantly distinct from other subjects. Abundance of Bacteroides and Fusobacterium was significantly positively correlated with the serum uric acid (UA) level of UAS patients. Fusobacteria was involved in the metabolism and degradation of certain short-chain fatty acids, amino acids, and sugars in pathopoiesis of UAS, and inhibited their synthesis pathways. Fusobacteria may be related to the pathogenesis of UAS, and this finding contributes to the personalized treatment of UAS from the perspective of maintaining micro-ecological equilibrium in gut.
Introduction
Nephrolithiasis is a global disease across regions and ethnicities, which is considered an important public health problem. It has a severe impact on human health and causes a huge social and economic burden (Sakhaee, 2008). To make matters worse, the recurrence rate of nephrolithiasis remains high, with a recurrence rate of 50% within 10 years. Renal uric acid stone (UAS) is indicated in the presence of a radiolucent stone (Ganesan et al., 2018), accounting for 10–15% of all urinary calculi, and its prevalence rate varies globally (Abou-Elela, 2017; Trinchieri and Montanari, 2017). The recurrence rate of patients with UAS is close to 100% (Wang et al., 2017). Although no large-scale epidemiological studies have been conducted, the prevalence of UAS in China has increased significantly in the past 30 years. In the economically developed areas of southern China, the prevalence of UAS is as high as 12–18% (Ma et al., 2018).
The microbiota of the human intestinal tract is a complex community composed of more than 100 trillion microbial cells, including more than 1,000 different kinds of species (De Sordi et al., 2017). High-throughput sequencing has facilitated great advances in our understanding of gut microbiome. Gut microbiota is increasingly linked to the development of various metabolic diseases such as obesity, diabetes, dyslipidemia, kidney disease, and nephrolithiasis (Langille et al., 2013; Stern et al., 2016; Stanford et al., 2020). The discovery of Oxalobacter formigenes, which is a kind of oxalate degrading bacteria, makes it possible that the gut microbiota affects absorption and secretion of solutes relevant to kidney stone formation (Siva et al., 2009). To date, more studies on the relationship between gut microbiota and nephrolithiasis are reported (Tang et al., 2018). They found that the cause of calcium oxalate stone was related to a group of bacteria involved in oxalate degradation and transportation, rather than a single kind of bacteria (Lee and Stern, 2019). The relationship of gut microbiota and calcium oxalate stone has been investigated in a limited amount (Stern et al., 2016; Lee and Stern, 2019), especially with no study of gut microbiota and UAS, to the best of our knowledge.
Although hyperuricemia is the common physiological and pathological bases of UAS and gout, only part of the population will develop UAS or gout, and the mechanism is still being studied. Patients with a history of gout have greater risk of forming UAS, and patients with obesity, diabetes, dyslipidemia, or other metabolic syndrome (Joosten et al., 2020). However, many patients with UAS are not accompanied by gout. The loss of bicarbonate leads to the formation of acidic urine, which is the most direct risk factor for UAS. Gastrointestinal diseases such as inflammatory bowel disease have been shown to cause bicarbonate loss clearly, and we have the reason to believe that gut microbiota and its metabolites regulate urine pH. We attempted to explore the differences in biochemical and gut microbiota features between UAS and gout patients, and find out bacteria genera, which involved in the pathogenesis of UAS. Through the analysis and detection of target bacteria by computational intelligence, it facilitates a new thinking for the personalized diagnosis and treatment of UAS patients.
Materials and Methods
Human Study Designs, Subjects, and Sampling
This study was approved by the Ethics Review Committee of Changshu Hospital affiliated to Soochow University. All the participants were local residents of unrelated southern Han Chinese. A case-control study of 123 gout patients (Gout group), 87 UAS patients without gout (UAS group), 30 UAS patients complicated with gout (Gout + UAS group), and 135 stone-free healthy people (Control Group) who received physical examination in Changshu Hospital affiliated to Soochow University from January 2014 to October 2020 were conducted. All the stone patients were diagnosed by ultrasound of the urinary system or abdominal-computed tomography, and received ureteroscopy or percutaneous nephrolithotomy. The stones obtained after operation were identified as pure or mixed uric acid calculus (uric acid content is greater than 50%) by Fourier-transformed infrared spectrophotometry LIIR-20 (Lanmode scientific instrument Co., Ltd., Tianjin, China). Patients who had a history of statins, chronic liver insufficiency, malignant tumors, and thyroid or parathyroid diseases were excluded from the study. Controls with a history or any evidence of nephrolithiasis, self-reported history of dyslipidemia or use of statins were excluded from the study.
General body measurements (height, weight, and blood pressure), past medical history, and test of blood and urine were carried out in all subjects. About 5 ml fasting venous blood sample and mid-stream urine of the first urine in the morning were drawn from each subject. Under the matching of age, sex, and BMI, 15 patients in the Gout group, 16 patients in the UAS group, 17 patients in the Gout + UAS group, and 17 patients in the controls were selected (Supplementary Table S1). All subjects were long-term residents with similar dietary habits, and underwent food frequency questionnaire before enrollment. After 3 days of eating uniform diet provided by the hospital (Supplementary Table S2), fresh fecal samples from the aforementioned people were collected, froze immediately, and stored under −80°C until analysis. Subjects who support fecal samples were excluded if they used antibiotics within 3 months, had a history of chronic diarrhea or constipation, chronic enteritis, irritable bowel syndrome, gastrointestinal tumor, or intestinal surgery.
16s DNA Extraction, PCR Amplification, and Target Gene Sequencing
The DNA was extracted from 200 mg samples using the QIAamp DNA fecal mini kit (QIAGEN, Hilden, Germany) following the manufacturer’s instructions. DNA was checked by running the samples on 1.2% agarose gels.
The V3–V4 hyper-variable regions of the bacteria 16SrRNA gene were amplified with primers 357F (5′-ACTCCTACGGRAGGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′) by the polymerase chain reaction (PCR) system. Prior to library pooling, the barcoded PCR products were purified using a DNA gel extraction kit (Axygen, United States) and quantified using the FTC-3000 TM real-time PCR (Funglyn Shanghai). The PCR products from different samples were indexed and mixed at equal ratios for sequencing on the illumina platform at TinyGen Bio-Tech (Shanghai) Co., Ltd.
Statistical and Bioinformatics Analysis
The general characteristics and biochemical parameters of the subjects were analyzed by SPSS22.0 and 2-sided p < 0.05 was defined as statistically significant. Continuous variables were summarized with mean and standard deviation if they satisfied the homogeneity of normality and variance. Otherwise, they were reported with median and inter quartile range (IQR). The quantitative variables were tested by the analysis of variance (ANOVA), Kruskal–Wallis test, and Mann–Whitney test. The chi-square test was used for categorical variables. Multiple logistic regression analysis was used to analyze the risk factors of disease in each group, and OR values were adjusted according to age and sex.
16S sequences were analyzed using a combination of software mother and R. The demultiplexed reads were clustered at 97% sequence identity into operational taxonomic units (OTUs) and the singleton OTUs were deleted using the UPARSE pipeline. OTU taxonomies were determined based on the NCBI. Based on taxonomic information, the community structure was statistically analyzed from the classification level of phylum, class, order, family, genus, and species. For the alpha-diversity analysis, Shannon, Simpson, Chao, ACE index, and rarefaction curves were calculated using mothur and plotted by R. The Kruskal–Wallis test was used to detect the significant changes of Shannon, Simpson, Chao, and ACE index between each group. For the beta-diversity metrics, the weighted UniFrac distance matrix were calculated using mothur and visualized with principal coordinate analysis (PCoA) by ape package in R. Linear discriminant analysis effect size (LEfSe) analysis was performed. Linear discriminant analysis was performed on samples with different grouping conditions according to the taxonomic composition. LDA was used to screen the communities or species which had significant influence on the differences in sample division with a cutoff of 2.0. The Spearman correlation analysis was used to evaluate association of dominant bacteria genera and biochemical parameters. The metabolic pathways of Kyoto Encyclopedia of Genes and Genomes (KEGG) were predicted by PICRUST software package.
Results
General Characteristics of the Study Population
The proportion of hypertensive patients and body mass index (BMI) in observation groups were significantly higher than the controls. Age (p = 0.571), gender (p = 0.310), and prevalence of diabetes mellitus (p = 0.413) were not significantly different among each group. However, BMI of the two stone groups, especially the Gout + UAS group was significantly higher than the other groups (p < 0.001). More hypertension individuals were found among the three observation groups and the Gout + UAS group is still with the highest prevalence rate (p < 0.001) (Table1).
Biochemical Features and Risk Factors Analysis Among the Groups
Hyperuricemia was the main risk factor for the development of gout as well as UAS. Obesity, hypertriglyceridemia, and aciduria were unique risk factors for UAS patients. Main biochemical features are summarized in Table 2 or Supplementary Figure S1. Compared to the controls, the triglycerides (TG) level in the three observation groups was significantly higher, and the level of TG in the Gout + UAS group was the highest (p < 0.05). The level of high-density lipoprotein cholesterol (HDL-C) in Gout as well as Gout + UAS group was significantly lower than that in the control group. The concentration of potassium in the Gout and UAS groups was significantly lower than that in the Gout + UAS group (p < 0.05), and the serum calcium concentration in the Gout group was also significantly lower than the Gout + UAS group and the controls (p < 0.05). All the three observation groups have significantly higher serum uric acid (UA) value than the controls (p < 0.001), and among them the Gout + UAS group’s UA value was the highest (p < 0.001). All the patients in three observation groups demonstrated obvious aciduria, especially in the Gout + UAS group (Table2, Supplementary Figure S1).
Univariate logistic analysis suggested that hypertension, low HDL-cholesterolemiamia, hyperuricemia, and aciduria were independent risk factors for gout. Obesity, hypertension, hypertriglyceridemia, low HDL-cholesterolemiamia, hyperuricemia, and aciduria were common independent risk factors for all UAS groups (Supplementary Table S3). Adjusted for age and sex, multivariate logistic analysis indicated that only hyperuricemia was the common risk factor for the development of gout as well as UAS. For patients with UAS alone, obesity, hypertriglyceridemia, and aciduria were the other three risk factors. While obesity was another risk factor for patients with UAS complicated with gout besides hyperuricemia (Table 3).
16S Sequencing Depth and Analysis of the Sample Size
The study was carried out with reasonable sample collection and high species rich, and the amount of data was reasonable and could reflect the majority of gut microbiota information objectively among groups. After eliminating the sequences of repetitive and fuzzy bases that affect the quality of analysis, a total of 2,329,290 high-quality sequences with an average length of 459 were received. The number of total OTUs in this study was 700, including one kingdom, 12 phylum, 21 class, 29 order, 46 family, 132 genus, and 181 species. The Venn diagram shows the common and unique OUTs among each group. Among the 648 known OUTs, the number of OUTs in the Gout + UAS group was lower than other groups (Figure 1A). With the curve reached the saturation plateau, more data contribute little to the discovery of new OUTs (Figure 1B), which means the depth of sequencing was reasonable. ANOSIM analysis calculated the relationship ranking among samples through variables, and performed the substitution test to determine whether the difference between groups was significantly different from the difference within groups. The results showed that the difference of the gut microbiota structure among the four groups was statistically significant (R = 0.155, p = 0.001) (Figure 1C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/888883/fphar-13-888883-HTML-r1/image_m/fphar-13-888883-g001.jpg)
FIGURE 1. Analysis of gut microbiota among each group by using 16s rRNA. (A) Venn diagram for indicating the common and unique OTUs among four groups. (B) Multi-sample rarefaction curves for comparing the abundance of diverse species in each sample. (C) ANOSIM analysis for identifying the existence of differences between each group. (D) Comparison of alpha diversity of gut microbiota between each group. The ACE, Chao, Shannon, and Simpson indices at the operational taxonomic units (OTUs) level were compared between each group by Kruskal–Wallis test (*p < 0.05**p < 0.01).
Richness and Diversity of Gut Microbiota Among Each Group
The richness and diversity of gut microbiota in UAS patients, especially in stone patients complicated with gout, were significantly lower than gout patients and controls. The richness and bacterial diversity of gut microbiota in the Gout + UAS group were significantly lower than that in the Gout and Control groups (Figure 1D) from the ACE, Chao, Shannon, and Simpson indices by the Kruskal–Wallis test. PCoA analysis indicated that there are significant differences in microbiota communities in specific evolutionary lineages among the four groups. (Supplementary Figure S2).
Analysis Based on Species Information of Microbiota Between Each Group
The species diversity and relative abundance of dominant bacteria at the phylum and genus levels among the four groups were different, and certain bacteria were dominant in UAS patients. At the level of phylum, Bacteroidetes, Firmicutes, Proteobacteria, Fusobacteria, Actinobacteria were the most five abundant bacterial in each group. In the Gout + UAS group, the relative abundance of Fusobacteria was significantly higher than other groups (p = 0.02), while Tenericutes was significantly lower than other groups (p < 0.001) (Figures 2A,B). At the level of genus, the most five abundant bacterial were Bacteroides, Prevotella, Megamonas, Fusobacterium, and Faecalibacterium. The relative abundance of Bacteroides and Fusobacterium in the Gout + UAS group was significantly higher than other groups (Figures 2C,D). The linear discriminant analysis effect size diagram indicated that there were the following specific genera in each group on the phylum and genus levels Fusobacteria in the Gout + UAS group on the phylum level, Tenericutes in the Gout Group on the phylum level, Bacteroides and Fusobacterium in the Gout + UAS group on the genus level, Streptococcus, Lactobacillus, Weissella, Gemella, and Campylobacter in the Gout group on the genus level, Dialister in the UAS group on the genus level, Subdoligranulum in controls on the genus level (Figure 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/888883/fphar-13-888883-HTML-r1/image_m/fphar-13-888883-g002.jpg)
FIGURE 2. (A,C) Composition of gut microbiota between each group at the phylum and genus level. (B,D) Histograms of differences in the abundance of species analyzed by Kruskal–Wallis test among four groups at the levels of phylum and genus (*p < 0.05, **p < 0.01, and ***p < 0.001).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/888883/fphar-13-888883-HTML-r1/image_m/fphar-13-888883-g003.jpg)
FIGURE 3. Histogram (A) and cladogram (B) of Linear discriminant analysis effect size (LEfSe) analysis based on OTUs characterizes microbiota among the controls, Gout patients, UAS patients and Gout + UAS patients.
Association of Dominant Bacteria Genera and Biochemical Parameters
Bacteroides and Fusobacterium were significantly positively correlated with the serum UA level of patients with UAS. We looked for the dominant species whose LDA value was more than two and the relative abundance was in the top 20 in each group at the genus level, and screened out the following six species: Bacteroides, Fusobacterium, Subdoligranulum, Streptococcus, Dialister, and Lactobacillus. The correlation of the aforementioned genera abundance and the significant different biochemical parameters were analyzed. Bacteroides was positively correlated with the UA level in the Gout + UAS group (r = 0.520, p = 0.033), and positively correlated with the serum potassium level in the controls (r = 0.829, p < 0.001). Fusobacterium was positively correlated with the UA level in the UAS group (r = 0.560, p = 0.024). Subdoligranulum was positively correlated with the HDL-C level in the controls (r = 0.493, p = 0.044). Streptococcus was positively correlated with urinary pH in the Gout group (r = 0.536, p = 0.039), and inversely correlated with the TG level in the Gout + UAS group (r = -0.517, p = 0.034). Dialister was positively correlated with the HDL-C level (r = 0.726, p = 0.001), but was inversely correlated with the UA level in the UAS group (r = -0.629, p = 0.009) (Table 4, Supplementary Figure S3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/888883/fphar-13-888883-HTML-r1/image_m/fphar-13-888883-t004.jpg)
TABLE 4. | Spearman correlation coefficient of significant different biochemical parameters and major bacterial genera in each group.
Prediction of Metabolic Function of Gut Microbiota
Fusobacteria was involved in the metabolism and degradation of certain short-chain fatty acids, amino acids, and sugars in pathogenetic of UAS, and inhibited their synthesis pathways. The PICRUSt tool was used to predict the differences in metabolic pathways between all the stone patients (Group S) and the controls (Group N) (Figure 4). Fusobacteria played a dominant role in microbiota metabolism in UAS patients and had a significant positive correlation in propanoate and butanoate metabolism, beta-alanine and tryptophan metabolism, amino sugar and nucleotide sugar metabolism, fructose and mannose metabolism, lysine, valine, leucine and isoleucine degradation and other metabolism or degradation pathway. The biosynthesis of short-chain fatty acids, the biosynthesis of phenylalanine, tyrosine, tryptophan and arginine, the biosynthesis of peptidoglycan and glucosinolate, the metabolism of cysteine, methionine, histidine, alanine, aspartate, and glutamate were inversely correlated with Fusobacteria while they were positively correlated with probiotic bacteria genera in controls.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/888883/fphar-13-888883-HTML-r1/image_m/fphar-13-888883-g004.jpg)
FIGURE 4. Correlation of different bacterial abundance and Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathways in the controls (Group N) and the renal uric acid stone patients (Group S). The red cell indicates a positive correlation and the blue cell indicates a negative correlation. Stars indicate the degree of significant correlations (*,p < 0.05,**,p < 0.01).
Discussion
In the present study, there were significant differences in biochemical features between UAS formers and gout patients. Whether for gout or UAS patients, hyperuricemia was always a risk factor. Obesity, hypertriglyceridemia, and aciduria were the risk factors of UAS, while obesity and hyperuricemia were more associated with UAS complicated with gout. The richness and diversity of gut microbiota in UAS patients were significantly lower than those in controls. Bacteroides and Fusobacterium were the dominant species that distinguish the intestinal flora of UAS patients from normal population, and they had a significant positive correlation with the level of serum UA. Fusobacteria was mainly involved in the metabolism and degradation of certain short-chain fatty acids, amino acids, and sugars in patients with UAS, and played an important role in inhibiting their synthesis pathways.
It is well known that uric acid is a weak organic acid with an ionization constant (PKA) of 5.5, the solubility of uric acid crystals decreases sharply in acidic urine with pH less than 5.5 (Tran and Maalouf, 2020). The disturbance of Na+/H+ exchange and decrease of ammonia secretion in renal proximal tubule will reduce the urinal pH value, promoting the formation of UAS and depositing of monosodium urate monohydrate, as known as the gout crystals (Friedlander et al., 2014; Rajkumar and Pluznick, 2018).
Our study indicated that the urine pH of UAS patients, especially stone formers with gout, was significantly lower than other subjects. The diseases that cause acidic urine included obesity, cardiovascular disease, insulin resistance, and chronic diarrhea (Bobulescu et al., 2019). The baseline characteristics of UAS patients in our study were partly consistent with the aforementioned clinic features. Epidemiological investigation has demonstrated that increased metabolic syndrome related factors significantly promoted the risk of UAS (Boyd et al., 2018). Several researchers have reported and confirmed that there was an inverse correlation between BMI and urine pH (Nakajima et al., 2016). Multivariate analysis was carried out in 459 24-h urine samples of 183 adult patients with nephrolithiasis and found that hyperglycemia was associated with lower pH and higher urinary saturation with respect to uric acid, suggesting that glycemic control may be considered a target for the treatment of UAS (Maciolek et al., 2018). In addition, a cross-sectional study analyzed 24-h urine parameters and abdominal CT scan results of 98 patients with nephrolithiasis and reported that the greater abdominal visceral fat content was associated with lower urinary pH and higher risk of UAS (Patel et al., 2017). Our results suggested that the incidence of diabetes in UAS patients with gout was higher than that in other groups, but there was no statistical significance, possibly due to an insufficient sample size.
Our study revealed significant differences in the levels of TG and HDL-C between the Gout + UAS group and the other three groups, which indicated a positive correlation between the severity of dyslipidemia and the risk of metabolic disorders in patients with UAS. This conclusion was roughly similar to our previous study’s results (Ding et al., 2019).
Nephrolithiasis is a known risk factor for chronic kidney disease (CKD), while UAS formation is associated with greater CKD risk compared with other stone types. Studies of patients with nephrolithiasis in Taiwan and Saudi Arabia found that the glomerular filtration rate (GFR) in UAS patients were significantly lower than those with other types of stones (Li et al., 2018; Nassir et al., 2018). While the serum potassium concentration of the Gout + UAS group was higher than other groups, suggesting decreased GFR and impaired renal tubule potassium excretion. On the flip side, the serum calcium concentration of the gout patients was the lowest, which was related to the calcium and phosphorus metabolism disorder caused by gouty nephropathy. Interestingly, there was no significant difference in the serum potassium level between gout patients and the controls. Glucocorticoid, as one of the treatment drugs for acute gout, can promote the sodium retention and potassium excretion of distal convoluted tubules and collecting tubules.
Under normal circumstances, the gut microbiota and host remains in a dynamic balance, which participates in various physiological processes of host nutrition absorption and metabolism. The gut microbiota can also affect the urine composition of the host, so the destruction of gut microbiota may lead to the occurrence of nephrolithiasis. Ticinesi et al. studied 52 patients with calcium oxalate nephrolithiasis and 48 healthy subjects, and found that the fecal microbial diversity in the stone group was significantly lower than that of the controls (Ticinesi et al., 2018). Stern et al. collected feces and urine samples of stone and non-stone patients, concluded that Bacteroides genus was the most abundant in stone formers and Prevotella genus was the most abundant in the control group. They also indicated that the composition of 24-h urine seemed to be related to the abundance of gut microbiota (Stern et al., 2016). A larger sample size and open-source project on the relationship between Oxalobacter formigenes and nephrolithiasis is being carried out in the United States (Liu et al., 2017).
Different from previous studies, we focused on the characteristics of gut microbiota in UAS patients for the first time. The healthy people rely on rich micro-ecosystem so as to maintain the necessary nutrition and metabolism, but with the aggravation of metabolic disorder, the intestinal microecology is destroyed. Bacteroides and Firmicutes accounted for more than 85% in the total amount of gut microbiota in our subjects, consistent with other experimental results at home and abroad (Ruan et al., 2020). However, there was a significant difference in the specific gravity of Fusobacteria at the phylum level. The relative abundance of Fusobacteria increased with the aggravation of metabolic disorders, and the specific gravity in Gout + UAS patients was the highest. At the genus level, Fusobacterium and Bacteroides were the dominant bacteria in patients with UAS, especially in the Gout + UAS group. Furthermore, Spearman correlation analysis results indicated that the abundance of Bacteroides and Fusobacterium was positively correlated with the serum UA level in UAS patients.
The involvement of the gut microbiota in multiple metabolic pathways in the host is widely recognized. Shotgun metagenomics sequencing technology is of course the most powerful approach to identify specific bacteria involved in metabolic pathways and identify gene function using the KEGG database. Liu et al. found a highly expressed bacterial gene in patients with recurrent calcium oxalate nephrolithiasis, which was involved in oxalate degradation and oxalate synthesis, and was related to high levels of urinary oxalate and acetic acid excretion (Liu et al., 2020). In this study, our PICRUSt analysis showed that Fusobacteria altered the microbial community functions, especially participated in the metabolism and degradation of certain short-chain fatty acids, amino acids, and sugars in UAS patients, but significantly inhibited their synthetic pathway. Future work including shotgun metagenomics analysis would help to confirm the specific gene functions.
Fusobacteria is a kind of anaerobic Gram-negative bacteria, which is widely colonized in human intestinal and oral mucosa and is related to the invasion of tumor cells (Kelly et al., 2018; Wang et al., 2020). There was a strong correlation between the abundance of Fusobacteria and the expression of pro-inflammatory markers such as COX-2, suggesting that Fusobacteria could create a pro-inflammatory micro-environment conducive to colorectal cancer by recruiting tumor-infiltrating immune cells (Kostic et al., 2013; Sears, 2018; King et al., 2020). Fusobacteria might potentially enhance the invasiveness of cancer cells when it existed in the micro-environment of oral tumor (Harrandah et al., 2020). In addition, studies have shown that the enrichment of Fusobacteria is related to acute appendicitis in children (Zhong et al., 2014; Rogers et al., 2016).
Our study indicated that the enrichment of Fusobacteria in the intestinal tract of UAS formers was higher than that of normal population for the first time. On one hand, hyperuricemia can improve the excretion level of uric acid in urine; on the other hand, it can cause dysfunction of the function of renal tubules to secrete acid and then form stones. The unique risk factors for UAS include persistent acidic urine and hyperuricuria, which can also be verified in biochemical parameters (Table 2) and multivariate logistic analysis (Table 3) in this study. Fusobacterium may be involved in the degradation and metabolism of certain short-chain fatty acids (Figure 4), and these metabolic pathways are related to the acid-secreting function of renal tubules. The specific expression of Fusobacterium may induce a pro-inflammatory micro-environment that gout patients do not have.
Whether Fusobacteria promotes the expression of inflammatory factors in renal tubular epithelial cells, whether these inflammatory factors promote the formation of stones caused by urine acidification, and whether the metabolic pathway involved in Fusobacteria can be targeted for pharmacological intervention? With the rapid development of computational intelligence, the application scenarios of artificial intelligence (AI) are becoming richer and richer. In the future, computational intelligence can be applied to pharmacological intervention of gut microbiota to treat metabolic diseases such as UAS. The computer obtained the basic data group of gut microbiota of target patients through meta-analysis of a large number of literatures. Through algorithm improvement, a suitable retrieval model was established to complete the application research of the AI-based gut microbiota intervention model for UAS. The computer intelligently learned the basic data processing in data mining, data processing, data cleaning and so on. Through the analysis, it can be concluded whether pharmacological intervention for gut microbiota of UAS is feasible based on data mining, so that computational intelligence can learn more pharmacological intervention plans on the basis of disease association rules. Research and development of drugs through AI can greatly shorten the time of drug development, improve the efficiency of research, control the cost of research, and improve the efficiency of doctors’ personalized treatment.
Until now, the number of subjects included in this study was not large enough, and the results need to be further verified by expanding the sample size. Second, there was no 24-h urine composition analysis for stone formers in this study. The last but not least, due to financial constraints, collected fecal samples were not tested for short-chain fatty acids and lacked metabonomic analysis.
Conclusion
Renal uric acid stone formers, especially complicated with gout, often have various types of dyslipidemia, persistent hyperuricemia, and aciduria. The richness and diversity of their gut microbiota were different from the gout patients as well as normal population. Bacteroides and Fusobacterium were positively correlated with the serum UA level of patients with UAS. The pro-inflammatory bacteria, Fusobacteria, may be related to the pathogenesis of UAS and has the potential to become the biomarker. These findings may provide a novel and non-invasive target for the prevention and treatment of UAS, which requires further large-scale investigation.
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 below: https://www.ncbi.nlm.nih.gov/Traces/study; SRP332626.
Author Contributions
CC, JZ, and D-RY contributed to conception and design of the study. CC analyzed the data and wrote the manuscript. Data collection was performed by BF, J-YC, and ZN. D-RY and JZ revised the manuscript. All authors read and approved the final manuscript.
Funding
This study was funded by the Suzhou Youth Science and Technology Project (KJXW2019060), the Changshu Medical Science and Technology Program (cswsq202002), the Suzhou Science and Technology Project (SKJY2021008) and the Suzhou Medical Research Key Talents Program (GSWS2021016).
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/fphar.2022.888883/full#supplementary-material
References
Abou-Elela, A. (2017). Epidemiology, Pathophysiology, and Management of Uric Acid Urolithiasis: A Narrative Review. J. Adv. Res. 8 (5), 513–527. doi:10.1016/j.jare.2017.04.005
Bobulescu, I. A., Park, S. K., Xu, L. H. R., Blanco, F., Poindexter, J., Adams-Huet, B., et al. (2019). Net Acid Excretion and Urinary Organic Anions in Idiopathic Uric Acid Nephrolithiasis. Clin. J. Am. Soc. Nephrol. 14 (3), 411–420. doi:10.2215/cjn.10420818
Boyd, C., Wood, K., Whitaker, D., and Assimos, D. G. (2018). The Influence of Metabolic Syndrome and its Components on the Development of Nephrolithiasis. Asian J. Urol. 5 (4), 215–222. doi:10.1016/j.ajur.2018.06.002
De Sordi, L., Khanna, V., and Debarbieux, L. (2017). The Gut Microbiota Facilitates Drifts in the Genetic Diversity and Infectivity of Bacterial Viruses. Cell Host Microbe 22 (6), 801–e3. doi:10.1016/j.chom.2017.10.010
Ding, Q., Ouyang, J., Fan, B., Cao, C., Fan, Z., Ding, L., et al. (2019). Association between Dyslipidemia and Nephrolithiasis Risk in a Chinese Population. Urol. Int. 103 (2), 156–165. doi:10.1159/000496208
Friedlander, J. I., Moreira, D. M., Hartman, C., Elsamra, S. E., Smith, A. D., and Okeke, Z. (2014). Comparison of the Metabolic Profile of Mixed Calcium Oxalate/uric Acid Stone Formers to that of Pure Calcium Oxalate and Pure Uric Acid Stone Formers. Urology 84 (2), 289–294. doi:10.1016/j.urology.2014.04.019
Ganesan, V., De, S., Shkumat, N., Marchini, G., and Monga, M. (2018). Accurately Diagnosing Uric Acid Stones from Conventional Computerized Tomography Imaging: Development and Preliminary Assessment of a Pixel Mapping Software. J. Urol. 199 (2), 487–494. doi:10.1016/j.juro.2017.09.069
Harrandah, A. M., Chukkapalli, S. S., Bhattacharyya, I., Progulske-Fox, A., and Chan, E. K. L. (2020). Fusobacteria Modulate Oral Carcinogenesis and Promote Cancer Progression. J. Oral Microbiol. 13 (1), 1849493. doi:10.1080/20002297.2020.1849493
Joosten, L. A. B., Crişan, T. O., Bjornstad, P., and Johnson, R. J. (2020). Asymptomatic Hyperuricaemia: a Silent Activator of the Innate Immune System. Nat. Rev. Rheumatol. 16 (2), 75–86. doi:10.1038/s41584-019-0334-3
Kelly, D., Yang, L., and Pei, Z. (2018). Gut Microbiota, Fusobacteria, and Colorectal Cancer. Diseases 6 (4). doi:10.3390/diseases6040109
King, M., Hurley, H., Davidson, K. R., Dempsey, E. C., Barron, M. A., Chan, E. D., et al. (2020). The Link between Fusobacteria and Colon Cancer: a Fulminant Example and Review of the Evidence. Immune Netw. 20 (4), e30. doi:10.4110/in.2020.20.e30
Kostic, A. D., Chun, E., Robertson, L., Glickman, J. N., Gallini, C. A., Michaud, M., et al. (2013). Fusobacterium Nucleatum Potentiates Intestinal Tumorigenesis and Modulates the Tumor-Immune Microenvironment. Cell Host Microbe 14 (2), 207–215. doi:10.1016/j.chom.2013.07.007
Langille, M. G., Zaneveld, J., Caporaso, J. G., McDonald, D., Knights, D., Reyes, J. A., et al. (2013). Predictive Functional Profiling of Microbial Communities Using 16S rRNA Marker Gene Sequences. Nat. Biotechnol. 31 (9), 814–821. doi:10.1038/nbt.2676
Lee, J. A., and Stern, J. M. (2019). Understanding the Link between Gut Microbiome and Urinary Stone Disease. Curr. Urol. Rep. 20 (5), 19. doi:10.1007/s11934-019-0882-8
Li, C. C., Chien, T. M., Wu, W. J., Huang, C. N., and Chou, Y. H. (2018). Uric Acid Stones Increase the Risk of Chronic Kidney Disease. Urolithiasis 46 (6), 543–547. doi:10.1007/s00240-018-1050-1
Liu, M., Koh, H., Kurtz, Z. D., Battaglia, T., PeBenito, A., Li, H., et al. (2017). Oxalobacter Formigenes-Associated Host Features and Microbial Community Structures Examined Using the American Gut Project. Microbiome 5 (1), 108. doi:10.1186/s40168-017-0316-0
Liu, Y., Jin, X., Hong, H. G., Xiang, L., Jiang, Q., Ma, Y., et al. (2020). The Relationship between Gut Microbiota and Short Chain Fatty Acids in the Renal Calcium Oxalate Stones Disease. FASEB J. 34 (8), 11200–11214. doi:10.1096/fj.202000786R
Ma, Q., Fang, L., Su, R., Ma, L., Xie, G., and Cheng, Y. (2018). Uric Acid Stones, Clinical Manifestations and Therapeutic Considerations. Postgrad. Med. J. 94 (1114), 458–462. doi:10.1136/postgradmedj-2017-135332
Maciolek, K. A., Penniston, K. L., Jhagroo, R. A., and Best, S. L. (2018). Successful Diabetic Control as Measured by Hemoglobin A1c Is Associated with Lower Urine Risk Factors for Uric Acid Calculi. J. Endourol 32 (8), 771–776. doi:10.1089/end.2018.0376
Nakajima, K., Oda, E., and Kanda, E. (2016). Latent Association between Low Urine pH and Low Body Weight in an Apparently Healthy Population. Scand. J. Clin. Lab. Invest. 76 (1), 58–63. doi:10.3109/00365513.2015.1092049
Nassir, A., Saada, H., Alnajjar, T., Nasser, J., Jameel, W., Elmorsy, S., et al. (2018). The Impact of Stone Composition on Renal Function. Urol. Ann. 10 (2), 215–218. doi:10.4103/ua.Ua_85_17
Patel, N. D., Ward, R. D., Calle, J., Remer, E. M., and Monga, M. (2017). Computerized Tomography Based Diagnosis of Visceral Obesity and Hepatic Steatosis Is Associated with Low Urine pH. J. Urol. 198 (5), 1085–1090. doi:10.1016/j.juro.2017.05.079
Rajkumar, P., and Pluznick, J. L. (2018). Acid-base Regulation in the Renal Proximal Tubules: Using Novel pH Sensors to Maintain Homeostasis. Am. J. Physiol. Ren. PhysiolRenal Physiol. 315 (5), F1187–F1190. doi:10.1152/ajprenal.00185.2018
Rogers, M. B., Brower-Sinning, R., Firek, B., Zhong, D., and Morowitz, M. J. (2016). Acute Appendicitis in Children Is Associated with a Local Expansion of Fusobacteria. Clin. Infect. Dis. 63 (1), 71–78. doi:10.1093/cid/ciw208
Ruan, W., Engevik, M. A., Spinler, J. K., and Versalovic, J. (2020). Healthy Human Gastrointestinal Microbiome: Composition and Function after a Decade of Exploration. Dig. Dis. Sci. 65 (3), 695–705. doi:10.1007/s10620-020-06118-4
Sakhaee, K. (2008). Nephrolithiasis as a Systemic Disorder. Curr. Opin. Nephrol. Hypertens. 17 (3), 304–309. doi:10.1097/MNH.0b013e3282f8b34d
Sears, C. L. (2018). The Who, where and How of Fusobacteria and colon Cancer. eLife 7. doi:10.7554/eLife.28434
Siva, S., Barrack, E. R., Reddy, G. P., Thamilselvan, V., Thamilselvan, S., Menon, M., et al. (2009). A Critical Analysis of the Role of Gut Oxalobacter Formigenes in Oxalate Stone Disease. BJU Int. 103 (1), 18–21. doi:10.1111/j.1464-410X.2008.08122.x
Stanford, J., Charlton, K., Stefoska-Needham, A., Ibrahim, R., and Lambert, K. (2020). The Gut Microbiota Profile of Adults with Kidney Disease and Kidney Stones: a Systematic Review of the Literature. BMC Nephrol. 21 (1), 215. doi:10.1186/s12882-020-01805-w
Stern, J. M., Moazami, S., Qiu, Y., Kurland, I., Chen, Z., Agalliu, I., et al. (2016). Evidence for a Distinct Gut Microbiome in Kidney Stone Formers Compared to Non-stone Formers. Urolithiasis 44 (5), 399–407. doi:10.1007/s00240-016-0882-9
Tang, R., Jiang, Y., Tan, A., Ye, J., Xian, X., Xie, Y., et al. (2018). 16S rRNA Gene Sequencing Reveals Altered Composition of Gut Microbiota in Individuals with Kidney Stones. Urolithiasis 46 (6), 503–514. doi:10.1007/s00240-018-1037-y
Ticinesi, A., Milani, C., Guerra, A., Allegri, F., Lauretani, F., Nouvenne, A., et al. (2018). Understanding the Gut-Kidney axis in Nephrolithiasis: an Analysis of the Gut Microbiota Composition and Functionality of Stone Formers. Gut 67 (12), 2097–2106. doi:10.1136/gutjnl-2017-315734
Tran, T. V. M., and Maalouf, N. M. (2020). Uric Acid Stone Disease: Lessons from Recent Human Physiologic Studies. Curr. Opin. Nephrol. Hypertens. 29 (4), 407–413. doi:10.1097/mnh.0000000000000610
Trinchieri, A., and Montanari, E. (2017). Prevalence of Renal Uric Acid Stones in the Adult. Urolithiasis 45 (6), 553–562. doi:10.1007/s00240-017-0962-5
Wang, P., Ding, S., Sun, L., Feng, Y., Guo, K., Zhu, Y., et al. (2020). Characteristics and Differences of Gut Microbiota in Patients with Different Traditional Chinese Medicine Syndromes of Colorectal Cancer and normal Population. J. Cancer 11 (24), 7357–7367. doi:10.7150/jca.50318
Wang, W., Fan, J., Huang, G., Li, J., Zhu, X., Tian, Y., et al. (2017). Prevalence of Kidney Stones in mainland China: A Systematic Review. Sci. Rep. 7, 41630. doi:10.1038/srep41630
Keywords: gut microbiota, uric acid stone, nephrolithiasis, 16srRNA, biomarker, computational intelligence
Citation: Cao C, Fan B, Zhu J, Zhu N, Cao J-Y and Yang D-R (2022) Association of Gut Microbiota and Biochemical Features in a Chinese Population With Renal Uric Acid Stone. Front. Pharmacol. 13:888883. doi: 10.3389/fphar.2022.888883
Received: 03 March 2022; Accepted: 01 April 2022;
Published: 19 May 2022.
Edited by:
Khin Wee Lai, University of Malaya, MalaysiaReviewed by:
Saneera Hemantha Kulathilake, Rajarata University of Sri Lanka, Sri LankaJing Xue, Wuxi People’s Hospital Affiliated to Nanjing Medical University, China
Copyright © 2022 Cao, Fan, Zhu, Zhu, Cao and Yang. 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: Dong-Rong Yang, ydrsuda@163.com; Jin Zhu, oceanzhu79@qq.com