Characterization of metal(loid)s and antibiotic resistance in bacteria of human gut microbiota from chronic kidney disease subjects
Biological Research volume 55, Article number: 23 (2022)
Human Gut Microbiota (HGM) is composed of more than one thousand species, playing an important role in the health status of individuals. Dysbiosis (an HGM imbalance) is augmented as chronic kidney disease (CKD) progresses, as loss of kidney function accelerates. Increased antibiotic use in CKD subjects and consumption of nephrotoxic heavy metals and metalloids such as lead, cadmium, arsenic, and mercury in tap water increases the dysbiosis state. Studies in people with stage 3 CKD are complex to carry out, mainly because patients are self-reliant who rarely consult a specialist. The current work focused on this type of patient.
Lead and arsenic-resistant bacteria were obtained from self-reliant (that stands on its own) stage 3 CKD subjects. Pathogen-related Firmicutes and Proteobacteria genus bacteria were observed. Resistance and potentiation of antibiotic effects in the presence of metal(loid)s in vitro were found. Furthermore, the presence of the following genes markers for antibiotic and metal(loid) resistance were identified by qPCR: oxa10, qnrB1, mphB, ermB, mefE1, arr2, sulll, tetA, floR, strB, dhfr1, acrB, cadA2k, cadA3k, arsC, pbrA. We observed a decrease in the number of metal resistance markers.
The presence of cadA and arsC genetic markers of antibiotics and metal(loid)s resistance were detected in samples from stage 3 CKD subjects. Lower gene amplification in advanced stages of CKD were also observed, possibly associated with a decrease in resident HGM during kidney disease progression.
The human gut microbiota (HGM) is acquired rapidly after birth and evolves with age. A healthy intestinal tract and HGM are relatively stable throughout adulthood, but as a person ages, external elements such as antibiotic use, diet, and cellular stress, among others [1,2,3] trigger the so-called dysbiosis by altering HGM homeostasis. Specifically, antibiotic use influences HGM composition by eliminating sensitive bacteria and generating selection pressure, which influences the multi-resistance antibiotic ability of microorganisms .
In terms of bacterial species, between 1000 and 1150, mostly uncultivable, species have been described in healthy adults. The HGM is composed of 6 predominant phyla: Firmicutes, Bacteroidetes, Actinobacteria, Proteobacteria, Fusobacteria, and Verrucomicrobia; with Firmicutes and Bacteroidetes corresponding to 90% of the intestinal bacterial population .
Although carbohydrates are the preferred substrates for fermentation in the proximal colon, proteins may also be used as an energy source . Amino acid degradation by the HGM occurs by oxidation deamination, producing NH3, short-chain fatty acids (SCFA), hydrogen sulfide, and metabolites including phenol, p-cresol, indol, and scatol derivative compounds . These compounds, in turn, produce microbial-derived uremic toxins that affect the progression of chronic non-communicable diseases, such as chronic kidney disease (CKD) [8, 9] and cardiovascular disease . Firmicutes and Bacteroidetes participate in SCFA synthesis, reducing the inflammatory state and immune responses .
Chronic Kidney Disease (CKD) affects approximately 10% of the world population . The progression is defined as the loss of function maintained for ≥ 3 months, expressed by a glomerular filtration rate (GFR) < 60 mL/min/1.73 m2. . CKD is closely associated with cardiovascular disease; cardiovascular mortality can be up to 15 times higher in CKD patients compared to the general population . Stage 3 CKD (CKD3) is particularly complex to manage, mainly because patients are self-reliant and rarely consult a nephrology specialist, therefore the number of CKD3 patients is unknown.
Daily exposure to arsenic, cadmium, mercury, and lead in tap water raises the possibility of increased microorganism resistance or tolerance to these toxicants . The presence of metals also contributes to increased antibiotic resistance through co-selection. Co-selection can occur due to co-resistance phenomena, where 2 or more resistance genes are in the same genetic element and regulated by the same promoter, or cross-resistance, where a single resistance mechanism confers resistance to different components simultaneously, such as the cadmium efflux pump that can expel beta-lactam antibiotics . These issues may be particularly important for people with CKD compared to healthy subjects. In Chile, maximum concentrations for metals permitted in tap water is defined by Chilean standard NCH 409/1 (0.01 mg/L arsenic, 0.01 mg/L cadmium, 0.001 mg/L mercury and 0.05 mg/L lead) .
Recently, profiles of antibiotic resistome and the microbial community in groundwater have been analyzed in areas where CKD of unknown etiology (CKDu) appear. Researchers have shown a positive correlation between CKDu, resistome, and microbial communities found in drinking water . The presence of CKDu in the population has also been associated, among other factors, with consumption of cadmium, arsenic, lead, and aluminum .
In this work, we explored the interaction of metal(oids) and antibiotics in CKD patients, which would reflect both metal(loid)s and antibiotic resistance and influences the number of genetic determinants associated with the metabolism of antibiotics and metal(loid)s present in microbiota. Our main goal was to analyze changes on HGM over CKD course and to characterize the resistance to metal(loid)s and antibiotics of microorganisms obtained from fecal samples. We tested whether there were changes in resistance in the different stages of disease progression. To our knowledge, this is the first work exploring the relationship between daily exposure to metal(loid)s and the presence of bacterial genetic markers of resistance to metal(loid)s and/or antibiotics in patients. We hypothesized that there would be a higher number of resistant microorganisms and genetic markers associated with metal(loid)s and antibiotics resistance, and this number increase over the course of CKD, we infer that the change would be related to the augmented presence of these toxicant substances because of kidney failure and increased antibiotic due to recurrent infections in CKD patients.
Materials and methods
Sixteen patients participated in this study, with 4 individuals per group: stage 3, stage 4, stage 5 and healthy controls. To select study subjects, clinical records and previous laboratory tests were reviewed. We selected patients who had values of glomerular filtration rate (eGFR) that allowed for stage classification based on serum creatinine concentration. Classification of subjects was defined between 30 and 60 mL/min/1.73 m2 of eGFR (KDIGO). We followed guidelines for defining CKD including (a) structural or functional kidney abnormalities persisting over 3 months; (b) persistent proteinuria including microalbuminuria (MAU), A1-A2 categories < 300 mg/g; and (c) decreased filtration rate (less than 60 mL/min per 1.73 m2 . The inclusion criteria included healthy controls and CKD subjects between 30 and 65 years of age. Kidney replacement therapy (hemodialysis or peritoneal dialysis) only among stage 5 CKD subjects. The exclusion criteria included persons with a history of chronic digestive and intestinal pathologies and/or malabsorption syndrome, autoimmune or chronic inflammatory disease, or cancer. Drug consumption including tobacco, alcohol, antibiotics (< 6 months of non-consumption), anti-inflammatories, and laxatives, except those prescribed for treatment of CKD and pregnancy were considered exclusion criteria. The sample size was calculated using the G*power 3.1 software considering study power of 0.95 and 0.05 error, using the effect size of 75%.
After reviewing the protocol, participants signed a consent form. All documents and protocols were reviewed and approved by the scientific ethics committee of the Eastern Metropolitan Health Service and San Sebastián University.
Nutritional and biochemical parameters
Nutritional parameters were measured by a nutritionist. Height and weight were measured using a stadiometer (Seca 217) and scale (Seca 876), respectively, and BMI calculated. The waist and hip circumference were measured using the tape measure and waist to hip ratio calculated. Bioimpedanciometry was performed using the XCONTACT 350 device and percent of lean body mass was measured (D1000-3). The following biochemical parameters were analyzed by an external laboratory: albuminuria (mg/24 h), creatinine (mg/L), urea (g/L), glucose (mg/dL), albumin (g/dL), calcium (mg/dL), and phosphorous (mg/dL) (Additional file 5: Table S1).
Collection and preservation of samples
Stool samples were collected separately according to group (stage 3, 4, or 5 or healthy controls). Samples were processed in class II biosafety cabinets, aliquoted in 1.5 mL tubes, and stored at − 20 °C until processing.
E. coli isolation protocol and determination of the minimum inhibitory concentration (MIC) of metal(loid)s
In 200 mg of previously aliquoted fecal samples, 1 mL of YCFAm [21,22,23] was added to into a 1.5 mL tube. The mix was homogenized by pipetting and samples centrifuged for 30 s at 13,000×g. The resulting supernatant was used for bacterial isolation procedures, growth curves, resistance/tolerance to metal(loid)s and antibiotics resistance assay.
Two hundred μL of the mix was spread in Eosin methylene blue Levine agar to grow at 37 °C for 16 h. Bacteria that showed a metallic green coloration was isolated in YCFAm medium, the colonies were used in the minimum inhibitory concentration (MIC) assays. 20 mM of Mercury (Hg), arsenic (As) and lead (Pb), and 50 mM cadmium (Cd) stock solution was used for preparing tenfold serial dilutions. Later, 10 µL of a saturated culture of bacteria was inoculated into each dilution under aerobic conditions. Because E. coli is considered a sensitive bacterium for metal(loid)s toxicant effects, the MIC of E. coli was used as a cut-off point to determine resistant/tolerant microorganisms for metal(loid)s in stool samples.
Determination of resistant microorganisms
Plates with 2% agar YCFAm plus E. coli MIC of lead and arsenic were used to seed 100 µL of stool samples and select resistant/tolerant for further experiment. Colonies that grew in the presence of each of the metal(loid)s analyzed were stored in 30% glycerol at − 80 °C. The number of colony-forming units per gram (CFU/g) was determined after incubating at 37 °C for 24 h. The count of viable cells was performed by Image pixel analysis in ImageJ 1.8 software.
Determination of antibiotic resistance in bacteria
Microorganisms that showed resistance/tolerance to indicated metal(loid)s were analyzed using the Kirby Bauer method . Antibiotic paper disks of gentamicin (10 µg), ampicillin (10 µg), ciprofloxacin (5 µg), cefazolin (30 µg), ceftazidime (30 µg), and gentamicin (120 µg) were placed on bacteria previously seeded on agar plates in presence of indicated metal(loid)s and incubated at 37 °C for 24 h. Growth inhibition zones were measured and compared with the standardized diameter by The Standards for Antimicrobial Susceptibility Testing (CLSI) .
Bacteria grown for 16 h in culture were used to take inoculum and diluted by 1:1000 in fresh YCFAm broth and cultivated. When necessary, the medium was supplemented with ½ MIC of E. coli of each metal(oid) toxicant tested in 48-well plates (SPL life sciences) at 37 °C for 18 h with orbital shaking. The dynamics of grown bacteria were estimated by optical density at 600nm.
DNA isolation procedure
Stool samples were thawed at 4 °C and 200 mg of sample was aliquoted in a 1.5 mL tube containing 200 mg of beads for mechanical disruption. TE buffer was added for resuspension, a vortex was applied at 2500 rpm for 10 min, then the sample was centrifuged for 1 min at maximum speed and transferred to a clean tube. Proteinase K was added to a final concentration of 1 mg/mL and incubated at 50 °C for 10 min. Later stages were performed according to manufacturer instructions (Power Fecal kit protocol Qiagen) and, finally, the eluted DNA was stored at – 20 °C until use.
16S rRNA amplicon sequencing
16S rRNA V4 region sequencing was carried out using Miseq sequencing (Illumina). Primers 515F (5′-GTGCCAGCMGCCGCGGTAA-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′) were used to generate a product of 250 bp with a high ratio between the observed and predicted communities (R2 = 0.95) . Taxonomic characterization of bacterial populations present in YCFAm colonies able to grow in metal(oid)s plates from the microbiota of CKD3 patients was performed by phylogenetic analysis on Rstudio (v3.53) . The quality of the sequence obtained by forward and reverse primers was determined using the DADA2 (v1.12.1)  package and the SILVA database  for phylogenetic assignation. Only operational taxonomic units (OTUs) that obtained a relative abundance > 0.2% of total readings and > 1000 readings per taxa were considered.
Selection and synthesis of primers for genetic determinant analysis
Primers that amplify antibiotic resistance genes were obtained from the literature  and an alignment was carried out to find main hits of these oligos in the MEGARes , CARD  and GenBank  databases. Selected primers were subsequently evaluated by “in silico” PCR using UGENE software (Unipro)  (Fig. 1). For metal(loid) resistance primers design, principal resistance genes belonging to each of the 6 phyla of HGM were aligned to determine the consensus sequence. Sequence primers were designed with Primer3 software [35, 36]. Software conditions included: amplicon between 150 and 300 bp, GC percentage from 40 to 80%, Tm between 57 and 62 °C, and starter size between 20 and 27 bp. Selected primers were analyzed with “in silico” PCR (Fig. 1). The list of end primers used in this work is provided in the supplementary material (Additional file 6: Table S2).
Total DNA was analyzed by real-time PCR (qPCR), using a bright green 2 × qPCR Mastermix PCR without ROX (ABM). PCR was performed using Quantstudio 3 equipment and melting curves were analyzed for each amplicon to ensure the presence of a unique product of the expected size.
Statistical analysis was performed using GraphPad prism V9.2 software, except those corresponding to genomic analysis, which was performed using ggplot and analyzed on Rstudio V 3.51. We used the ANOVA test, considering the significance of P < 0.05 (*), P < 0.01 (**) or P < 0.001 (***).
Results and discussion
In this work we carried out exploratory research regarding the presence of microorganisms resistant to meta(oid)s and antibiotics in patients with CKD. We also performed an analysis of total DNA of stool samples to find resistance determinants in each patient, depending on CKD stage.
To analyze resistance/tolerance of metal(loid) and resistance to antibiotics in bacteria, samples from patients with CKD were analyzed in YCFAm culture medium using the MIC calculated from E. coli isolated from fecal samples (Additional file 7: Table S3). Our results exhibited colony growth in the presence of arsenic and lead under aerobic conditions (Additional file 1: Figure S1). No colonies were found on plates containing YCFAm medium in the presence of Hg and Cd (data not shown). Colonies harvested from each plate were used to perform an antibiotic resistance analysis.
Bacterial isolates were analyzed using the Kirby Bauer method in the presence of paper disks with antibiotics and characterized according to CLSI parameters for resistant, intermediate, or sensitive bacteria (Fig. 2A and B) in the presence or absence of ½ MIC of metal(oid)s, a sublethal concentration. For healthy subjects, a decreased appearance of gentamicin resistance in the presence of metal(oid)s, ceftazidime, and ciprofloxacin was observed compared with healthy subjects, which indicates a significant increase in resistance to antibiotics in the presence of metal(loid)s. For CKD3 patients, no significant difference in antibiotic resistance was observed compared with the healthy controls and ½ MIC of E. coli for metal(loid)s (arsenic and mercury both). However, higher numbers of microorganisms resistant to antibiotics ampicillin, cefazolin, and ciprofloxacin were shown in plates containing ½ MIC of E. coli metal(loid)s. No differences between groups (CKD3 versus healthy controls) and toxicant presence were observed when growth curves of stool samples in the presence of ½ MIC of metal(loid)s (Additional file 2: Figure S2).
To identify microorganisms in fecal samples that showed co-resistance to antibiotics and metal(loid)s, sequences of isolated colonies were analyzed by 16S amplicon sequencing. All isolates of healthy subjects and CKD subjects from each disease stage (3, 4, and 5) were compared after antibiotic treatments in the presence of ½ MIC of E. coli of arsenic or lead. For CKD3 patients and lead treatment, Firmicutes and Proteobacteria were mainly found (Fig. 3A). Genus analysis showed Pseudomonas spp., Janibacter spp., Escherichia/Shigella spp., and Bacillus spp.; resistant bacteria to gentamicin, cefazolin, ceftazidime, and ampicillin were also found. For CKD3 patients and arsenic treatment Pseudomonas spp., Escherichia spp./Shigella spp., Bacillus spp. were observed; and Enterococcus spp. showed resistance to ampicillin, ciprofloxacin, and gentamicin. Furthermore, some colonies of interest exhibited resistance to all tested antibiotics (general resistance). On the other hand, for healthy controls and lead treatment, we observed mainly Bacillus spp. and Pseudomonas spp. and resistance to ampicillin, cefazolin, gentamicin, and colonies with a general resistance. Healthy controls and arsenic treatment exhibited Bacillus spp., Escherichia spp./Shigella spp. and Pseudomonas spp. (Fig. 3A–C). It was impossible to detect microorganisms of Phylum Bacteroidetes, mainly because of other selection phenomena such as the culture medium used (YCFAm) and aerobic conditions. Parameters of microorganism abundance in samples grown in the synthetic and selective media were performed, Chao1, Shannon, and Simpson indices were calculated [S2A for metal(loid)s and S2B for antibiotics]. Alpha diversity indexes were associated with the selection pressures presented in healthy patients and CKD3, although, these parameters were not analyzed in depth.
We found that the Bacillus spp. in the presence of metal(loid)s. showed multiple resistance to ampicillin, gentamicin, cefazolin, ceftazidime and ciprofloxacin Pseudomonas spp. showed the same behavior in the presence of metal(loid)s, exhibiting the same antibiotic resistance, but less representation.
On the other hand, Enterococcus spp. and Escherichia spp./Shigella spp. exhibited resistance to ampicillin and gentamicin, respectively. We noted that both healthy controls and CKD3 patients showed some bacteria with general resistance to antibiotics, but resistance differed mainly in the type of metal (arsenic for CKD stage 3 and lead for a healthy group). Most studies that have found this co-resistance phenomena in stool samples have been conducted in wastewater systems . The phenomenon of co-selection allows conserving and promoting resistance to antibiotics and resistance/tolerance of metal(loid)s in bacteria, even in the absence of antibiotics through different co-resistance mechanisms, which opens the way to future investigations .
In taxonomic profiling, we identified the presence of Janibacter spp., which has been recently reported to be involved in the virulence process related to antibiotic resistance in HGM of celiac and immunocompromised subjects , expressing proteins such as β-lactamase and other 7 genes related to metal(loid) resistance/tolerance including arsenic . This led us to perform further analysis to characterize the isolated strain and identify the genes when this gut pathogen was exposed to metal(loid)s and antibiotics.
Beta diversity analysis between control and CKD stage 3 isolated bacteria showed no significant differences when stool samples were cultured under aerobic conditions and when exposed to metalloids and different antibiotics (Fig. 3C). This could be explained by a continuous process of the multi-resistance capacity of bacteria. Interestingly, this result suggests that the early stages of CKD do not influence the development of multi-resistant bacteria. To explore this idea in depth, we performed an analysis to detect resistance markers by qPCR.
It was initially suggested that CKD patients have augmented resistance genes in HGM because of a possible accumulation or transcendent increase of toxicant compounds such as metal(loid)s related to decreased renal filtration rate and the use of antibiotics for disease treatment. To determine the presence of resistance genes and genetic determinants related to metal(loid)s and antibiotics resistance in stool samples of healthy patients and to compare with subjects at different stages of CKD (3, 4, and 5), qPCR was performed (Additional file 6: Table S2). Primers were constructed to amplify regions of selected genes for metal(loid) resistance analysis corresponding to arsenic (arsC, arsA), lead (pbrA), mercury (merA), and cadmium (cadA). Specifically, for cadmium, 2 primers cadA2k and cadA3k were designed because in silico PCR results against MEGARes and CARD databases showed amplicons of different sizes and quantities. To increase coverage, both primers were synthesized.
To choose primers for antibiotic resistance genes, 140 described genes were analyzed  and 13 genes that represented each family were selected. For genes such as floR, a gene product resistance to chloramphenicol, a second product was initially detected in melting curves (not shown). To correct this bias, qPCR was performed with an annealing temperature gradient (60, 65, and 70 °C). It was observed that at 65 °C, only a single specific product was obtained. sulll primers of Sulfonamides dihydropteroate synthetase described by the group of Szczepanowski et al.  could not be used due to the detection of nonspecific products when analyzing melting curves in our samples (not shown). This could mean nonspecific interactions in other DNA regions, and therefore nonspecific amplification in our samples, producing false-positive results.
Pie charts (Fig. 4) describe the presence of antibiotic and metal(loid) resistance genes. Interestingly, the healthy control group (Fig. 4A) showed a higher number of resistance markers of metal(loid)s and antibiotics. In the same group, all analyzed genes for resistance to antibiotics and metals were found. The lack of some markers could not be evidenced in this group in comparison with the different stages of CKD. From CKD stage 3 samples (Fig. 4B), the presence of strB, dhfr1, floR, acrB, and arr2 for antibiotic resistance genes and cadA3k, arsC, and cadA2k for metal(loid)s resistance genes were observed. In the same group mefE1, catB4 and qnrB1 were not found. For stage 4 samples (Fig. 4C), the presence of acrB, arr2, qnfrB1, strB, dhfr1, floR, ermB, and tetA antibiotics were found and cadA3k, arsC, cadA2k, and pbrA for metal(loid) resistance genes, ermB was found again in stage 4. Concerning metal resistance genes, arsA and pbrA were not observed in CK3 subjects, however pbrA was identified in stage 4 patients. For stage 5 CKD samples (Fig. 4D), which correspond to subjects who no longer have renal function and require kidney replacement therapies, such as hemodialysis or peritoneodialysis and Firmicutes versus Bacteroidetes ratio tends to decrease because of the uremic state and increased proteolytic metabolism [7, 8], both groups of genes, for resistance to antibiotics and metal(loid)s were diminished compared with healthy subjects. The presence of qnrB1, dhfr1, and floR genes for resistance to antibiotics and cadA2k and merA for metal resistance was observed. We consider that the typical state of intestinal microbiota dysbiosis associated with CKD produces changes in microbiota that are associated with a smaller number of genetic determinants.
The study of the presence of metal(loid) resistance genes and antibiotics from the early stages of the disease is important because, at that stage, there are alterations in HGM composition -or dysbiosis- that could explain the decrease in the appearance of the analyzed genetic markers. Other alterations associated with the progression of disease include intestinal transit alterations, decreased protein absorption, decreased consumption of dietary fiber, and frequent use of antibiotics. These accumulative factors will contribute to the production of systemic inflammation and the bioaccumulation of uremic toxins and presumably metal(loid)s, which are absorbed by the intestine and eliminated by the kidney, these toxicants may play a central role in the physiopathology of CKD .
A high presence of resistance genes for cadmium and arsenic metabolism were observed in analyzed samples (Fig. 5A). Cadmium (Cd) is a heavy metal, usually present in soils . It is toxic to living organisms, is carcinogenic to humans [41,42,43,44], and lethal . The human body can absorb Cd in small portions through food intake, especially in grains, water, or air , where it accumulates and remains for a long time, causing health problems [43,44,45,46,47]. Cd accumulates in the liver and kidneys and has a long biological half-life, 17–30 years in humans. Toxicity involves 2 organ systems, the renal and skeletal, and is largely the consequence of the interactions between Cd and essential metals, particularly calcium [43, 45, 47]. The severity and damage of these metals depend on time, level of exposure, and susceptibility of the person, and which metal is absorbed .
Mining activities in the second region of Chile represents one of the main commercial and labor activities of the country. This industry emits waste rich in arsenic and other minerals, which affects both the environment and residents of the area [48, 49]. Arsenic is found food, water, and air. The main damage produced by arsenic exposure is cardiovascular, kidney, neurological, respiratory, skin cancer, and reproductive effects .
High levels of the dhfr1 gene (dihydrofolate reductase) (Fig. 5B) can be observed on the heat diagram for both the healthy and CKD group, which corresponds to a resistance gene present in some Escherichia coli strains (https://www.uniprot.org/uniprot/Q0MQM2), the main causative microorganism of urinary tract infection (UTI) in healthy and CKD population .
To determine whether the observed changes in the content of antibiotic and metal resistance genes were significant between groups, a correlation (Pearson R) of Contingency of Prospective data (chi-square test) was performed and plotted by heat map (Additional file 4: Figure S4A) and principal component analysis PCA (Additional file 4: Figure S4B). Our results showed that healthy subjects (control), CKD3, and stage 4 CKD subjects are considered similar groups with a positive correlation. On the contrary, the stage 5 CKD group did not show a correlation with the previous groups. The results show that the gene content for resistance to antibiotics and metals is substantially lower than in the previous groups, Thus, further studies are needed to determine how late stages of CKD relate to a decrease in resistance markers.
The incidence of UTI in CKD patients increases with disease progression and occurs mainly due to the presence of risk factors. For example, uremia causes alterations in the humoral response, lymphocyte function, macrophages, and polymorphonuclear cells [51,52,53].
The underlying symptom of CKD is sometimes a condition that compromises normal elimination of urine and integrity of the urinary tract or implies its manipulation due to medical procedures (e.g., vesicoureteral reflux, neurogenic bladder, urethral valves, prostatism, bladder catheterization, renal catheterization, complicated lithiasis, and polycystic disease). In other cases, diabetes is an underlying symptom of both CKD and greater susceptibility to the appearance of UTI and its evolution and occurs especially in elderly female patients .
This study represents a slightly closer look at the association between the microorganisms present in the intestine and the evolution of a chronic disease. Even though the number of patients included in this study is small, it was possible to show an important modification in the content of microorganisms and resistance markers present in the microbiota. We also demonstrated the possibility of using microorganisms as biomarkers of a person’s health status or biosensors of any toxicant exposure through the study of microbial genetic markers. These genetic markers are indicative of disease progression and/or can indicate the type of metal(loid) or antibiotic to which the patient has been exposed due to genetic determinants in the microorganism of HGM. This description is particularly important in patients suffering an important decrease in one of the most important detoxification systems of the body, the liver.
To our knowledge, this is the first work that shows a relationship between the presence or absence of antibiotic resistance markers, determined by qPCR, among subjects at different stages of CKD. Our results suggest modification of intestinal microbial resistome as a function of disease progression. The presence of antibiotic and metal(loid) resistance gene markers were found. cadA and arsC were detected in CKD3 and decreased gene amplification in advanced stages of CKD were observed, possibly associated with a decrease in resident HGM.
Although associated morbidities such as diabetes can result in CKD, there are other associated factors, such as metal(loid)s consumption . According to recent studies, the resistome could influence the development and progression of CKDu , metal(loid)s exposure and the presence of antibiotic resistance genes could change over the course of CKD, thus future research is needed to advance the understanding of this relationship.
Availability of data and materials
Chronic non-communicable diseases
Human gut microbiota
Chronic kidney disease
Chronic kidney disease of unknown etiology
Stage 3 chronic kidney disease
Glomerular filtration rate
Urinary tract infection
Renal replacement therapy
Short-chain fatty acids
Modified yeast chain fatty acids medium
Estimated glomerular filtration rate
Vemuri R, Gundamaraju R, Shastri MD, Shukla SD, Kalpurath K, Ball M, et al. Gut microbial changes, interactions, and their implications on human lifecycle: an ageing perspective. BioMed Res Int. 2018. https://doi.org/10.1155/2018/4178607.
Jandhyala SM, Talukdar R, Subramanyam C, Vuyyuru H, Sasikala M, Reddy DN. Role of the normal gut microbiota. World J Gastroenterol. 2015;21:8836–47.
Adak A, Khan MR. An insight into gut microbiota and its functionalities. Cell Mol Life Sci. 2019;76(3):473–93.
Jernberg C, Löfmark S, Edlund C, Jansson JK, Lofmark S, Edlund C, et al. Long-term impacts of antibiotic exposure on the human intestinal microbiota. Microbiology. 2010;156:3216–23.
Kim S, Covington A, Pamer EG. The intestinal microbiota: antibiotics, colonization resistance, and enteric pathogens. Immunol Rev. 2017;279(1):90–105.
Macfarlane GT, Macfarlane S. Factors affecting fermentation reactions in the large bowel. Proc Nutr Soc. 1993;52:367–73.
Smith EA, Macfarlane GT. Dissimilatory amino acid metabolism in human colonic bacteria. Anaerobe. 1997;3:327–37.
Meijers BK, Evenepoel P. The gut-kidney axis: indoxyl sulfate, p-cresyl sulfate and CKD progression. Nephrol Dial Transplant. 2011;26:759–61.
Evenepoel P, Poesen R, Meijers B. The gut-kidney axis. Pediatr Nephrol. 2017;32:2005–14.
Lekawanvijit S. Cardiotoxicity of uremic toxins: a driver of cardiorenal syndrome. Toxins. 2018;10(9):352.
Chang PV, Hao L, Offermanns S, Medzhitov R. The microbial metabolite butyrate regulates intestinal macrophage function via histone deacetylase inhibition. Proc Natl Acad Sci USA. 2014;111(6):2247–52.
Nasri H. World kidney day 2014; chronic kidney disease and aging: a global health alert. Iran J Pub Health. 2014;43(1):126–7.
Ministerio de Salud, Chile. Guía Clínica Prevención Enfermedad Renal Crónica; 2010.
Ministerio de Salud, Chile. Guías de prácticas clínicas GES prevención secundaria de la enfermedad renal crónica; 2017.
Breton JÔ, Daniel C, Dewulf J, Pothion S, Froux N, Sauty M, et al. Gut microbiota limits heavy metals burden caused by chronic oral exposure. Toxicol Lett. 2013;222(2):132–8.
Pal C, Asiani K, Arya S, Rensing C, Stekel DJ, Larsson DGJ, Robert KP, et al. Chapter seven—metal resistance and its association with antibiotic resistance. Adv Microb Physiol. 2017;70:261–313.
Salazar C, Etcheberrigaray M, Cáceres A, Echeverría E, Pizarro J, Bruna A, Lisboa E, Lozano V, Parra X, Santa-Ana S, Givovich A. Norma Chilena Oficial. NCh409/1.Of.2005—Parte 1. Requisitos para Agua Potable. Chile: Instituto Nacional de Normalización; 2006. p. 1–13.
Cooray T, Zhang J, Zhong H, Zheng L, Wei Y, Weragoda SK, et al. Profiles of antibiotic resistome and microbial community in groundwater of CKDu prevalence zones in Sri Lanka. J Hazard Mater. 2021. https://doi.org/10.1016/j.jhazmat.2020.123816.
Lunyera J, Mohottige D, von Isenburg M, Jeuland M, Patel UD, Stanifer JW. CKD of uncertain etiology: a systematic review. Clin J Am Soc Nephrol. 2016;11:379–85.
Willis K, Cheung M, Slifer S. KDIGO 2012 Clinical practice guideline for evaluation & management of CKD. Kidney Int Suppl. 2013. https://doi.org/10.1038/kisup.2012.77.
Duncan SH, Barcenilla A, Stewart CS, Pryde SE, Flint HJ. Acetate utilization and butyryl coenzyme A (CoA): acetate-CoA transferase in butyrate-producing bacteria from the human large intestine. Appl Environ Microbiol. 2002;68(10):5186–90.
Alou MT, Naud S, Khelaifia S, Bonnet M, Lagier JC, Raoult D. State of the art in the culture of the human microbiota: new interests and strategies. Clin Microbiol Rev. 2020;34(1):e00129-e219.
Naud S, Khelaifia S, Mbogning Fonkou MD, Dione N, Lagier JC, Raoult D. Proof of concept of culturomics use of time of care. Front Cell Infect Microbiol. 2020;10: 524769.
Dunkelberg WE. Kirby-Bauer disk diffusion method. Am J Clin Pathol. 1981;75:273.
Clinical and Laboratory Standards Institute. Clinical and Laboratory Standards Institute (CLSI). Performance standards for antimicrobial susceptibility testing. 30th ed. Wayne: CLSI Supplement M100; 2020.
Parada AE, Needham DM, Fuhrman JA. Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ Microbiol. 2016;18(5):1403–14.
Team Rs. RStudio: integrated development for R studio Inc. Boston: R for Beginners; 2005.
Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13(7):581–3.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.
Szczepanowski R, Linke B, Krahn I, Gartemann KH, Gützkow T, Eichler W, et al. Detection of 140 clinically relevant antibiotic-resistance genes in the plasmid metagenome of wastewater treatment plant bacteria showing reduced susceptibility to selected antibiotics. Microbiology. 2009;155(Pt 7):2306–19.
Lakin SM, Dean C, Noyes NR, Dettenwanger A, Ross AS, Doster E, et al. MEGARes: an antimicrobial resistance database for high throughput sequencing. Nucleic Acids Res. 2017;45:D574–80.
Alcock BP, Raphenya AR, Lau TTY, Tsang KK, Bouchard M, Edalatmand A, et al. CARD 2020: antibiotic resistome surveillance with the comprehensive antibiotic resistance database. Nucleic Acids Res. 2020;48(D1):D517–25.
Sayers EW, Cavanaugh M, Clark K, Pruitt KD, Schoch CL, Sherry ST, et al. GenBank. Nucleic Acids Res. 2021;49:D92–6.
Okonechnikov K, Golosova O, Fursov M, Varlamov A, Vaskin Y, Efremov I, et al. Unipro UGENE: a unified bioinformatics toolkit. Bioinformatics. 2012;28(8):1166–7.
Koressaar T, Remm M. Enhancements and modifications of primer design program Primer3. Bioinformatics. 2007;23(10):1289–91.
Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3-new capabilities and interfaces. Nucleic Acids Res. 2012;40(15): e115.
Yamina B, Tahar B, Laure FM. Isolation and screening of heavy metal resistant bacteria from wastewater: a study of heavy metal co-resistance and antibiotics resistance. Water Sci Technol. 2012;66(10):2041–8.
Fernández-Natal MI, Sáez-Nieto JA, Medina-Pascual MJ, Valdezate-Ramos S, Guerra-Laso JM, Rodríguez-Pollán RH, et al. First report of bacteremia by Janibacter terrae in humans. Infection. 2015;43(1):103–6.
Chander AM, Kochhar R, Dhawan DK, Bhadada SK, Mayilraj S. Genome sequence and comparative genomic analysis of a clinically important strain CD11-4 of Janibacter melonis isolated from celiac disease patient. Gut Pathogens. 2018;10:2.
Cigarran Guldris S, González Parra E, Cases AA. Microbiota intestinal en la enfermedad renal crónica. Nefrol Soc Esp Nefrol. 2017;37:9–19.
McDowell RW, Taylor MD, Stevenson BA. Natural background and anthropogenic contributions of cadmium to New Zealand soils. Agric Ecosyst Environ. 2013;165:80–87. https://doi.org/10.1016/j.agee.2012.12.011.
Fowler BA. Monitoring of human populations for early markers of cadmium toxicity: a review. Toxicol Appl Pharmacol. 2009;238(3):294–300.
Chang YF, Wen JF, Cai JF, Xiao-Ying W, Yang L, Guo YD. An investigation and pathological analysis of two fatal cases of cadmium poisoning. Forensic Sci Int. 2012;220(1–3):e5-8.
Liu Y, Xiao T, Baveye PC, Zhu J, Ning Z, Li H. Potential health risk in areas with high naturally-occurring cadmium background in southwestern China. Ecotoxicol Environ Saf. 2015. https://doi.org/10.1016/j.ecoenv.2014.10.022.
Goyer RA. Toxic and essential metal interactions. Annu Rev Nutr. 1997;17:37–50.
Grant CA, Buckley WT, Bailey LD, Selles F. Cadmium accumulation in crops. Can J Plant Sci. 1998. https://doi.org/10.4141/P96-100.
Wångstrand H, Eriksson J, Öborn I. Cadmium concentration in winter wheat as affected by nitrogen fertilization. Eur J Agron. 2007. https://doi.org/10.1016/j.eja.2006.09.010.
Ferreccio C, Sancha AM. Arsenic exposure and its impact on health in Chile. J Health Popul Nutr. 2006;24:164–75.
Ferrante M, Conti GO, Rasic-Milutinovic Z, Jovanovic D. Health effects of metals and related substances in drinking water. Water Intell Online. 2013. https://doi.org/10.2166/9781780405988.
Becar HC, Becar HC, Venegas G. Impacto y consecuencias del Arsénico en la salud y el medio ambiente en el Norte de Chile. Rev Interam Ambiente Tur RIAT. 2012. https://doi.org/10.4067/riatvol6iss1pp53-60%250718-235X.
García-Agudo R, Panizo N, Proy Vega B, García Martos P, Fernández RA. Infección del tracto urinario en la enfermedad renal crónica. Rev Colomb Nefrol. 2020. https://doi.org/10.22265/acnef.7.1.387.
Muacevic A, Adler J. Urinary tract infection in chronic kidney disease population: a clinical observational study. Cureus. 2021;13(1): e12486.
Kuo IC, Lee JJ, Hwang DY, Lim LM, Lin HY, Hwang SJ, Chen HC, Hung CC. Pyuria, urinary tract infection and renal outcome in patients with chronic kidney disease stage 3-5. Sci Rep. 2020;10(1):19460. https://doi.org/10.1038/s41598-020-76520-5
We want to thank the collaborators and institutions that participated in this project, particularly those who facilitated patient recruitment: Dialysis Center Ñuñoa, Providencia, Chile; Family Health Center Santa Julia, Macul, Chile; Family Health Center Dr. Felix de Amesti, Macul, Chile; Family Health Center Padre Alberto Hurtado, Macul, Chile.
This work was funded by FONDECYT 11170885, (Principal Investigator Dr. Waldo Díaz Vásquez).
Ethics approval and consent to participate
This project was approved by the scientific ethics committee of the Eastern Metropolitan Health Service and San Sebastian University, Patients signed an informed consent document.
Consent for publication
The authors declare their consent for this publication.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Colony count of viable bacteria (CFU/mL) of healthy and CKD3 stool samples seeded on YCFAm agar cultivated under aerobic conditions in the presence of E. coli MIC of arsenic or lead.
Growth curves of bacteria from the microbiota of stool samples from subjects with CKD3 and healthy controls. Growth curves from bacteria of fecal samples obtained from subjects with CKD3 and healthy subjects were analyzed in YCFAm medium with and without metal(loid) supplementation of ½ MIC of E. coli. No significant differences were found between treatments in each group, n = 12.
Measures of alpha-diversity. Observed Shannon and Simpson indexes were calculated for A Metal(loid)s, and B antibiotic resistance isolate colonies obtained from healthy and CKD3 stool samples. Indexes were calculated using the phyloseq package in RStudio, significance was calculated by ANOVA and Wilcoxon test (p = 0.3227) for the Shannon index.
Pearson correlations of contingency of prospective data (chi-square test) for the total appearance of genes of metal(loid)s and antibiotics resistance among healthy controls and, stage 3, 4, and 5 CKD patients. Results were expressed by a heat map (Additional file 4: Figure S4A) and principal component analysis PCA (Additional file 4: Figure S4B).
Anthropometric and biochemical parameters of study participants.
Primers used for antibiotic and metal(loid) resistance gene determination by qPCR.
Minimum inhibitory concentration of isolated E. coli for cadmium, mercury, and arsenic under aerobic and anaerobic conditions. Assays were carried out with 3 technical repetitions and 3 biological repetitions for each metal.
About this article
Cite this article
Miranda, M.V., González, F.C., Paredes-Godoy, O.S. et al. Characterization of metal(loid)s and antibiotic resistance in bacteria of human gut microbiota from chronic kidney disease subjects. Biol Res 55, 23 (2022). https://doi.org/10.1186/s40659-022-00389-z