Microbial genetic and transcriptional contributions to oxalate degradation by the gut microbiota in health and disease
Abstract
Over-accumulation of oxalate in humans may lead to nephrolithiasis and nephrocalcinosis. Humans lack endogenous oxalate degradation pathways (ODP), but intestinal microbes can degrade oxalate using multiple ODPs and protect against its absorption. The exact oxalate-degrading taxa in the human microbiota and their ODP have not been described. We leverage multi-omics data (>3000 samples from >1000 subjects) to show that the human microbiota primarily uses the type II ODP, rather than type I. Furthermore, among the diverse ODP-encoding microbes, an oxalate autotroph, Oxalobacter formigenes, dominates this function transcriptionally. Patients with inflammatory bowel disease (IBD) frequently suffer from disrupted oxalate homeostasis and calcium oxalate nephrolithiasis. We show that the enteric oxalate level is elevated in IBD patients, with highest levels in Crohn’s disease (CD) patients with both ileal and colonic involvement consistent with known nephrolithiasis risk. We show that the microbiota ODP expression is reduced in IBD patients, which may contribute to the disrupted oxalate homeostasis. The specific changes in ODP expression by several important taxa suggest that they play distinct roles in IBD-induced nephrolithiasis risk. Lastly, we colonize mice that are maintained in the gnotobiotic facility with O. formigenes, using either a laboratory isolate or an isolate we cultured from human stools, and observed a significant reduction in host fecal and urine oxalate levels, supporting our in silico prediction of the importance of the microbiome, particularly O. formigenes in host oxalate homeostasis.
Introduction
Over-accumulation of oxalate in humans leads to toxicity (Asplin et al., 1998; Beck et al., 2013). The most common oxalate toxicity is calcium oxalate nephrolithiasis, which accounts for more than 70% of overall nephrolithiasis, affecting 9% of the US population with a 20% 5-year recurrence rate (Lieske et al., 2014; Rule et al., 2014; Saran et al., 2018; Scales et al., 2012; Stamatelou et al., 2003). Oxalate toxicity can also induce chronic kidney disease, an illness affecting more than 30 million Americans, via multiple mechanisms including the activation of the NALP3 inflammasome pathway (Knauf et al., 2013; Mulay et al., 2017; Mulay et al., 2013; Saran et al., 2018; Waikar et al., 2019), RIPK3-MLKL-mediated necroptosis (Mulay et al., 2016a), and oxidative stress-induced cell injury (Khan et al., 2006). In extreme cases, life-threatening systemic oxalosis occurs.
Humans lack endogenous oxalate-degrading enzymes. By contrast, the mammalian intestinal microbes can degrade oxalate, partially protecting their hosts against toxicity (Allison and Cook, 1981; Allison et al., 1986; Allison et al., 1977; Barber and Gallimore, 1940; Miller et al., 2014; Azcarate-Peril et al., 2006) Oxalobacter formigenes,Lactobacillus sp., Bifidobacterium sp., Enterobacteriaceae, and others can degrade oxalate in vitro (Mogna et al., 2014), and colonization with these taxa in rodent hyperoxaluria models showed reduction in urinary oxalate indicating oxalate degradation in vivo (Hatch et al., 2006; Klimesova et al., 2015; Kwak et al., 2006). Of these oxalate degraders, only O. formigenes is a specialist that uses oxalate as its sole energy source (Cornick and Allison, 1996). O. formigenes also is unique because it induces host oxalate secretion into the colonic lumen (Arvans et al., 2017).
The oxalate-degrading microbes in the human microbiota in vivo have not been characterized. Previous studies on oxalate-degrading microbes have been chiefly done in vitro or in animal models, and the relevance of those microbes to human health remains undefined (Fontenot et al., 2013; Klimesova et al., 2015; Kullin et al., 2014; Turroni et al., 2010; Turroni et al., 2007). This gap limits our understanding of the role of microbiota in diseases induced by oxalate toxicity. Here, we leveraged multi-omics data of the healthy human microbiome to characterize the oxalate-degrading microbes in vivo.
Inflammatory bowel disease (IBD) patients are at increased risk for oxalate toxicity, due to a condition called enteric hyperoxaluria (EH). In EH, enhanced bioavailability and hyperabsorption of intestinal oxalate result in oxalate nephrolithiasis (Corica and Romano, 2016; Liu and Nazzal, 2019). In the USA, >50,000 IBD patients suffer from EH and recurrent calcium oxalate kidney stones (Corica and Romano, 2016; McConnell et al., 2002). EH in IBD patients may reflect lipid malabsorption and increased gut permeability. However, an alternate hypothesis is that microbiota-based oxalate degradation is impaired in dysbiotic IBD patients, leading to increased oxalate absorption (Allison et al., 1986). We interrogated the microbiota of IBD patients to understand shifts in microbiota-based oxalate degradation functions and their metabolic consequences.
Results
Type I and type II microbial oxalate degradation pathways
To determine the oxalate degradation pathways (ODPs) used by human gut bacteria, we curated all experimentally validated microbial ODP from literature review and database searches (see Materials and methods) (Allison et al., 1985; Anand et al., 2002; Blackmore and Quayle, 1970; Daniel et al., 2004; Dumas et al., 1993; Foster et al., 2012; Pierce et al., 2010). We classified those ODPs into two types based on their enzymatic mechanisms and co-factor requirements. Type I ODPs cleave the oxalate carbon-carbon (C-C) bond in a single step (Figure 1A). The two recognized type I enzymes, oxalate oxidase and oxalate decarboxylase, are indistinguishable at the amino acid level; therefore, we refer to them jointly as oxalate oxidase/decarboxylase (OXDD) (Figure 1A). Type II ODPs consists of two enzymatic reactions requiring coenzyme A as co-factor (Figure 1A). First, a coenzyme A molecule is added to oxalate to form oxalyl-CoA via enzymes including formyl-CoA transferase (FRC) (Figure 1A). In the second step, oxalyl-CoA decarboxylase (OXC) metabolizes oxalyl-CoA into CO2 and formyl-CoA (Figure 1A).
Knowing the relevant oxalate degradation enzymes (ODEs), we then acquired all available protein homologs of the three ODE OXDD (n = 2836), FRC (n = 1947), and OXC (n = 1284), which enable homology search from UniProt Interpro (Mitchell et al., 2019; Mulder et al., 2005). By tracing the taxonomic origin of the genes encoding those homologs, we found that OXDD-coding taxa can be fungal or bacterial, whereas FRC- and OXC-coding taxa are strictly bacterial (Figure 1B). The frequent co-occurrence of FRC and OXC in individual genomes indicates encoding complete type II ODP (Figure 1B). As expected, OXDD, FRC, and OXC each are conserved within the same microbial class, but exhibit substantial divergence across classes (Figure 1—figure supplement 1). These data provide both a comprehensive inventory of ODPs and a reference set of ODP-encoding microbes, which enable analyses to elucidate those relevant to humans.
ODPs utilized by the gut microbiota of healthy humans
Next, we asked whether those ODPs are encoded or expressed by microbes within the intestinal tract of healthy humans. To do so, we computationally examined the presence of the ODE within the gut metagenome and metatranscriptome from publicly available samples of healthy humans. From five studies, we analyzed a total of 2359 metagenome and 1053 transcriptome samples from 660 and 165 healthy individuals, respectively (Figure 2—figure supplement 1, Supplementary file 1a). After quality-filtering (see Materials and methods), the sequencing reads were aligned to the unique OXDD, FRC, and OXC homologous proteins we had identified, using DIAMOND Blastx (Buchfink et al., 2015). Alignment pairs with >90% identity were retained for downstream analyses. The alignment cutoff was based on the protein identity of the inter- and intra-species ODEs and determined to be robust for distinguishing ODEs originating from differing microbial species (Figure 2—figure supplement 2).
We found that the majority of the healthy gut microbiomes include ODEs with at least one ODE detected in the metagenome of 607 (92%) of 660 subjects and the metatranscriptome of 132 (80%) of 165 subjects. In the metagenomes, the type II frc and oxc genes were more common (Figure 2A) and more abundant (Figure 2B) than the type I oxdd gene. Similarly, in metatranscriptomes, expressions of type II genes were more common and abundant (Figure 2C,D). Expression of oxdd was only detected in 10 (6%) of the 165 subjects, and the median RPKM was 2-log10 lower than those of frc or oxc (Figure 2C,D). Furthermore, frc and oxc – the coding genes of the two enzymes that catalyze the separate steps in type II PDP – were frequently co-expressed within the same microbiota, indicating expression of the complete type II ODPs (Figure 2—figure supplement 3).
These data indicate that microbes utilizing type II rather than type I ODPs predominate in the human intestine. Such finding is consistent in all studies, despite the differences in source populations and sample preparation methods (Supplementary file 1a). For the remaining analyses, we focused on the type II ODPs.
Microbial species that transcribe the type II ODPs in vivo
Although multiple human commensal microbes are known to encode frc and oxc, whether they transcribe those genes in vivo has not been studied (Abratt and Reid, 2010; Cho et al., 2015; Fontenot et al., 2013; Mogna et al., 2014). We next characterized the microbial species transcribing these type II ODP genes in the microbiota of healthy humans. In the metagenomes of 660 individuals, oxc gene of multiple species, including Escherichia coli, O. formigenes, and several Muribaculaceae, Bifidobacterium, and Lactobacillus sp., was detected (Figure 3A, left); E. coli oxc was the most common (56% of subjects), followed by O. formigenes oxc (39% of subjects).
In the metatranscriptomes of 167 individuals, oxc gene expression did not directly correlate with corresponding gene abundance or prevalence (Figure 3A, right). O. formigenes oxc expression was both most abundant and most prevalent (in 61% of subjects) in the metatranscriptomes (Figure 3A, right). Despite the detection of E. coli oxc in 56% of subjects, its transcript was present in only 12% of the subjects (Figure 3A). For Bifidobacterium and Lactobacillus species, for which oxalate degradation activity was reported in vitro and in animal models (Federici et al., 2004; Klimesova et al., 2015; Turroni et al., 2010), oxc expression was minimal (<5%) or not prevalent (Figure 3A right). The dichotomy between metagenomic and metatranscriptomic oxc was consistent across different studies (Figure 3—figure supplement 1, Figure 3—figure supplement 2), and also present for the other type II ODP gene frc (Figure 3—figure supplement 3).
To more rigorously examine ODP expression by individual taxa, we further co-analyzed the presence of ODP genes and transcripts matched by subject (Figure 3B). O. formigenes frc and oxc were transcribed in nearly all subjects in whom the genes were detected, as well as in others in whom the gene was not detected, indicating that their expression is common in vivo (Figure 3B). The under-detection of O. formigenes ODP genes in the metagenomes may reflect the highly variable abundance of the organism, often below the lower detection limit using gene-based methods (Barnett et al., 2016; Guo et al., 2017; Kelly et al., 2011; PeBenito et al., 2019). In contrast, E. coli frc and oxc were expressed in only a few subjects even when the corresponding genes were detected metagenomically (Figure 3B). These data indicate that in vivo E. coli rarely transcribe ODP, a pathway used for defense against oxalate-induced stress (Fontenot et al., 2013). In total, these findings demonstrate that ODP transcription varies widely in individual hosts, and by species.
The contributions of individual species to the global microbiota ODP
We then assessed the impact of individual species on global ODP by quantifying their population-level contributions (see Materials and methods). The contribution of O. formigenes to ODP increased from 17% to 63% from the metagenomic to the metatranscriptomic level, greater than the transcriptomic contributions of all other species combined (Figure 3C, Supplementary file 1b). Conversely, the E. coli contribution to ODP was markedly reduced from the metagenomic (36%) to the metatranscriptomic (7%) level (Figure 3C). Other species had low but varied contributions (Figure 3C). A parallel pattern was observed for frc (Figure 3—figure supplement 4, Supplementary file 1c). With the low activity of non-O. formigenes species, network analysis did not yield significant species-species interactions related to oxc transcription (data not shown). In summary, we found that the type II ODP genes, frc and oxc, are encoded by multiple gut microbes, but O. formigenes dominated this pathway at the transcriptional level. These data provide a baseline for examining disease-associated changes.
Validation of ODP detection using ShortBRED
ShortBRED (Kaminski et al., 2015) is a tool for microbiome functional profiling, which clusters protein homologs into clusters, and identifies marker peptide for each cluster, thus potentially achieving high specificity. We reanalyzed all samples using ShortBRED to validate our bioinformatics findings (see Materials and methods). Based on ShortBRED, FRC and OXC were significantly more abundant and more prevalent than OXDD; their abundances were significantly correlated with our previous results (Figure 3—figure supplement 5A). Consistently, O. formigenes was the species with the highest transcriptional activity for FRC (Figure 3—figure supplement 5B). Specifically, each of the three O. formigenes FRC homologs (C3 × 9Y2, C3 × 762, and C3 × 2D4) are distinct from other homologs and from each other; thus, each formed a singleton family (Supplementary file 2) with unique peptide markers (Supplementary file 3). The three O. formigenes FRCs are the most commonly transcribed among FRCs encoded by any taxon. They are present in the metatranscriptome of 50, 52, and 41% of the tested subjects, respectively (Figure 3—figure supplement 5B, right), and also are among the most abundant (ranks 1, 2, and 5) among all FRC proteins (Figure 3—figure supplement 5B, right).
The taxonomic contributions based on BURRITO (McNally et al., 2018) (see Materials and methods) (Supplementary file 1d) are highly concordant with findings using our prior method (Supplementary file 1b, metagenome, Healthy). E. coli and O. formigenes are consistently the two largest contributors, and several Bifidobacterium spp and Lactobacillus spp contributed to a lesser extent.
ShortBRED did not detect O. formigenes OXC in any sample. We reasoned that this under-detection is due to poorly-selected marker peptides. Indeed, for the protein cluster of O. formigenes OXC (consists of two O. formigenes OXC homologs: C3XBB9 and C3 × 545, with C3 × 545 as the centroid), the makers selected are as follows: (1)short (median length 16.2 amino acids-in comparison, the marker for O. formigenes FRC is 300 amino acids) and (2) improperly broken up from long continuous regions (Figure 3—figure supplement 6A, Marker #1–20). Therefore, we manually generated new markers (Markers #21–24) by combining markers that were one amino acid apart. Using these longer markers, ShortBRED successfully detected OXC using three of the four markers (except for the shortest one) with the expected patterns (Figure 3—figure supplement 6B,C).
In summary, all of our major conclusions examined were confirmed with ShortBRED. The diamond-mapping method utilizes reference full length and relies on best alignment scores to assign reads with higher sensitivity, while ShortBRED uses unique regions of reference with higher specificity. Because of their different advantages, the two methods are complementary and both were used for downstream analyses.
Increased enteric oxalate levels and reduced microbial ODP expression in IBD patients
IBD patients, particularly patients suffering from CD, frequently have EH, with oxalate hyperabsorption and calcium oxalate nephrolithiasis (Corica and Romano, 2016; Liu and Nazzal, 2019; McConnell et al., 2002). Ulcerative colitis (UC), regardless of severity and location, is associated with stone formation (Cury et al., 2013), but ileocolonic CD is associated with greater nephrolithiasis risk than either ileal or colonic involvement alone (Cury et al., 2013). We hypothesized that oxalate degradation by the intestinal microbiome may be impaired in IBD patients, leading to more luminal oxalate in the host available for passive absorption.
We tested this hypothesis, using the multi-omics data of IBD patients and healthy subjects from the iHMP-IBD study (Franzosa et al., 2019; Lloyd-Price et al., 2019). The patients were stratified by illness: UC (N = 30 subjects) and CD (N = 54), and the CD group was further divided into CD-L3 (N = 25), with ileocolonic phenotype at baseline and CD-nonL3 (N = 29) without, as defined by the L3 IBD Montreal classification (Satsangi et al., 2006). Consistent with the clinical nephrolithiasis risk, fecal oxalate relative abundances were elevated in both the UC (p=0.005) and CD (p=0.06) patients compared to healthy controls (Figure 4A). All of the CD risk was in the CD-L3 patients (p<0.001), and not in the non-CD-L3 patients, indicating that IBD location, particularly ileocolonic involvement, is key for EH risk (Figure 4A, Figure 4—figure supplement 1). The fecal oxalate levels were not clearly different in relation to disease activity (Damms and Bischoff, 2008; Jowett et al., 2001; Manz et al., 2012; Pathirana et al., 2018; Figure 4—figure supplement 2). We observed increased fecal oxalate in association with higher inflammation levels (fecal calprotectin > 50 µg/g) (Figure 4—figure supplement 2C); since our number of study subjects was small, this relationship needs to be examined with larger patient cohorts in future studies.
The global transcripts of frc and oxc were reduced in the IBD patients compared to the controls in all four studies analyzed (Figure 4B, Figure 4—figure supplement 3). Expression of oxc was detected in 57% UC, 54% CD and 48% CD-L3 subjects, lower than the 71% in healthy individuals (Figure 4B). ODP expression was least impacted in the CD-nonL3 group (Figure 4B, Figure 4—figure supplement 3), in which the fecal oxalate increase was not observed (Figure 4A). The total transcripts of frc and oxc were also significantly lower in IBD patients (Figure 4C, Figure 4—figure supplement 3B) (p<0.001 for frc for all groups compared to healthy, p<0.01 for oxc in UC patients compared to healthy). We also observed that frc expression was significantly inversely correlated with oxalate relative abundance in the UC group (Figure 4—figure supplement 4A). The same inverse trend was marginally observed for oxc but was not statistically significant (Figure 4—figure supplement 4B); thus, this question will need to be examined in larger future cohorts. In contrast to reduced frc and oxc transcripts in IBD patients, the frc and oxc genes were significantly more abundant (Figure 4—figure supplement 5), indicating that the frc- and oxc-encoding taxa (e.g. Enterobascteriaceae) are enriched in the IBD gut but do not actively express ODP.
Collectively, these data showing the reduction of ODP-associated transcripts in the IBD patient microbiota suggest their role in the elevated intestinal oxalate levels, and possibly increased susceptibility to nephrolithiasis.
Loss of O. formigenes and its ODP-associated transcripts in IBD patients
We next sought to identify the microbial species accounting for the reduced ODP transcripts in the IBD patients. E. coli and O. formigenes with the largest ODP contributions at the genomic and transcriptional level, respectively, were notable. Using gene and transcript jointly as markers, ODP expression by O. formigenes was detected in ~25% of UC and CD patients (Figure 5A), significantly less than in healthy persons (~70%), either when studies were combined (Figure 5A) or separate (Figure 3B). In contrast, E. coli ODP was detected in nearly all IBD subjects and was transcribed more frequently compared to healthy subjects (Figure 5A). Consistent with the low overall prevalence, transcripts for O. formigenes ODP expression were less abundant (lower RPKM values) in UC, CD, and in CD-L3 patients compared with controls (p<0.01 for all groups) (Figure 5B). In the IBD subjects, the observed O. formigenes ODP genes were always actively expressed (Figure 5A). Significantly elevated fecal calprotectin levels were observed when O. formigenes was absent in healthy individuals and in CD-L3 patients (Figure 4—figure supplement 2D).
In total, these data indicate that the absence of O. formigenes colonization or colonization below the level of detection is responsible for the reduction in global ODP transcripts. E. coli, and Lactobacillus and Bifidobacterium spp. use ODPs to defend against oxalate-induced acid stress; their upregulation in IBD appears secondary to the elevated oxalate levels present.
Effect on O. formigenes colonization on the host urinary and fecal oxalate levels
To validate our bioinformatics prediction that O. formigenes is an important oxalate-degrading organism that can influence oxalate homeostasis in vivo, we examined whether the colonization of O. formigenes results in significant reduction in urinary oxalate in mice (Figure 6A). Mice from our gnotobiotic facility which had a total microbiota 2–3 log10 lower than conventional mice were fed with diet supplemented with 1% sodium oxalate and 0.5% calcium. The 1% dietary oxalate is lower than previous studies (Hatch et al., 2006; Hatch et al., 2011), in order to reduce chronic kidney damage described previously (Mulay et al., 2016b). This diet resulted in significant increase of oxalate (23.3% and 174.8% increase in fecal and urinary oxalate, respectively) compared to normal chow (data not shown). We then colonized these mice with either of two O. formigenes strains, including a widely studied lab strain (OXCC13) and a human O. formigenes isolate from the stool sample of a patient with primary hyperoxaluria type 1 (PH1) isolated as described (Allison et al., 1985). Continued colonization with each was established by qPCR (Figure 6—figure supplement 1). Using targeted oxalate assays, the O. formigenes-colonized mice had significantly lower urinary and fecal oxalate levels, compared to the non-colonized mice (Figure 6B–C); the values from the two measurements tended to correlate in individual mice (Figure 6D).
In our model, fecal O. formigenes averaged ~5×106 by qPCR (Figure 6—figure supplement 1). In contrast, O. formigenes was below the detection limit (<102) in the jejunum and ileum of most mice (Figure 6—figure supplement 2). These data suggest that O. formigenes primarily colonizes and degrades oxalate in the host colon, which is known to be an important site for oxalate absorption.
Taken together, these findings suggest that O. formigenes is sufficient to reduce host urinary and fecal oxalate levels, and support our in silico prediction of the importance of O. formigenes to host oxalate homeostasis, with protection against oxalate-induced toxicity.
Discussion
Oxalate degradation by the human microbiota has been known since the 1940s (Allison and Cook, 1981; Allison et al., 1986; Allison et al., 1985; Allison et al., 1977; Barber and Gallimore, 1940), but the taxa involved in vivo has not been systematically described. We present the first comprehensive study of human oxalate-degrading microbes and define their individual contributions. We successfully distinguished the taxa that are actively transcribing ODP from those that encode the pathway, but with low expression, by co-analyzing metagenome and metatranscriptome data.
Our finding that multiple human gut microbes encode ODPs is consistent with prior studies (Abratt and Reid, 2010; Klimesova et al., 2015; Mogna et al., 2014; Stern et al., 2016; Ticinesi et al., 2018). But surprisingly, at the transcriptional level, O. formigenes dominates the global ODP, which is consistent with O. formigenes being an oxalate autotroph, as well as FRC and OXC being the most abundant proteins during both exponential and lag stages (Ellis et al., 2016). These findings were consistent using our method and ShortBRED, which have their unique analytic strengths. The diamond-mapping method utilizes protein full-length information and is completely agnostic. By setting a stringent identity cutoff of 90%, we retain only the high-confidence alignment pairs. The finding generated using this 90% cutoff is consistent with the observations generated through the second method ShortBRED, which uses ODE-specific marker peptides, which provide evidence that the cutoff is highly effective at preventing spurious alignments. By contrast, ShortBRED generated highly specific markers but is limited by several parameter choices (CD-Hit, centroid protein selection, length and identity of short identical regions and final marker). Because of their specific advantages, the two methods are complementary.
The contrasting genetic and transcriptional ODP differences we observed highlight the importance of analysis beyond the gene level for microbiome studies. There are two critical advantages of an approach that co-analyzes metagenome and metatranscriptome: (1) Detection of genes in the metagenome does not ensure that they are being actively utilized by the indicated taxon. Thus, analyzing metatranscriptomic data allows distinguishing the taxa that are actively contributing to a biological process by generating the relevant proteins. In contrast, for those taxa that are not, the analysis provides a new view of their metabolism, in that although they have the potential for using a certain pathway or metabolite, it is not active under the particular condition tested. (2) Having metatranscriptomic data or other functional readout is particularly useful for cross-comparisons of the microbiome of different cohorts (i.e., diseased versus healthy subjects, in humans or experimental animals). For example, in our study, ODP genes were increased in IBD patients, whereas transcripts were significantly reduced. Specifically, the increased ODP gene abundances were due to the over-representation of E. coli strains in IBD patients, which uncommonly transcribe ODP in vivo. In contrast, the decrease of oxalate-degrading gene expression is caused by the loss of O. formigenes, which is the dominant microbe that transcribes this pathway. Having observed the differential abundance for a gene does not necessarily indicate functional shifts, as gene abundance is driven by the most abundant taxa, which could be transcriptionally silent. As such, the transcriptional evidence is a better indicator for evaluation of microbiome functional differences. The integrative multi-omics analysis framework built for this study, (now deposited on Github), can be extended to a broad range of microbiome functions.
EH is frequent in IBD patients, particularly in CD, and in those who underwent Roux-en-Y gastric bypass; the latter population is optimal for study since pre- and post-treatment samples can be easily collected. In the present study, we examine IBD patients, a risk population for EH and nephrolithiasis, and show that the impaired metabolic activity of the microbiota is correlated with enteric oxalate levels, thus nephrolithiasis risk. That fecal oxalate relative abundances were elevated in the iHMP-IBD UC patients, and in the CD patients with ileocolonic involvement, is consistent with their high nephrolithiasis risk (Cury et al., 2013; Hylander et al., 1979; Hylander et al., 1978). In this population at risk for EH, the ODP gene abundances are increased, but the extent of gene expression is decreased. This contrast between gene abundance and expression is due to differential ODP transcription in particular microbes.
IBD patients are known to have low levels of O. formigenes colonization (Kumar et al., 2004). The oxalate-degrading specialist O. formigenes is a gram-negative anaerobe susceptible to multiple antibiotics including macrolides, tetracyclines, metronidazole, and clindamycin (Kharlamb et al., 2011; Mittal et al., 2005). Thus, the reduced colonization of O. formigenes in IBD patients might reflect the frequent antibiotic treatments they receive. Human data showed that a single course of antibiotics to eradicate Helicobacter pylori results in the long-term suppression of O. formigenes colonization (Kharlamb et al., 2011). This supports the notion that antibiotics exacerbate microbiome dysbiosis and may lead to secondary conditions such as EH. These findings raise a further potential cost of antibiotic treatment in IBD patients, which should be weighed against potential benefits, especially when infection is not clearly demonstrated. Other changes in the intestinal milieu of IBD patients affecting pH, oxygen, and bile acid levels have been shown to impact O. formigenes in vitro (Duncan et al., 2002; Allison et al., 1985).
Our findings that O. formigenes is the main contributor to oxalate degradation in the healthy state, but is diminished in the IBD population, provides a strong rationale for O. formigenes-based restoration therapy. Restoring O. formigenes has been tested in primary hyperoxaluria (PH) patients, but yielded mixed results (Hoppe et al., 2011; Hoppe et al., 2017). However, PH patients might not be the best subjects for restoration because their oxalate toxicity is caused by hepatic overproduction of oxalate, and microbial degradation is restricted to the oxalate secreted into the gut lumen. In EH, the oxalate overload is in the intestinal lumen, which provides a nutrient-rich environment for microbes such as O. formigenes to (re)colonize and degrade oxalate prior to absorption into the host circulation. A recent study (Canales and Hatch, 2017) of a surgery-induced EH rat model via Roux-en-Y gastric bypass (RYGB), indicated the potential of O. formigenes for treating EH in the post-RYGB rat, O. formigenes strain OXWR achieved 100% colonization and decreased urinary oxalate by 74% compared to 39% in the sham-operated group. However, levels of O. formigenes colonization in EH patients need to be established in future studies. Although PH is a rare disease, EH in IBD patients is common and growing more so, with currently limited treatment options.
Although our study has significant findings with the potential for translational and mechanistic studies, we acknowledge study limitations. The fact that we only used sequences of proteins instead of the whole genomes during read mapping could lead to false taxa assignments due to possible horizontal gene transfer events. Our analysis was limited to the currently characterized ODPs; therefore, we cannot rule out the existence of other enzymes in the human microbiota to degrade oxalate. Our observations largely rely on metatranscriptomics data, and therefore could be limited by technical biases. It is more difficult to acquire high-quality metatranscriptome, as RNA is less stable and subject to degradation during sample preparation. However, the high consistency across the separate studies (Figure 2, Figure 2—figure supplement 3, Figure 3—figure supplement 1) done using different sample collection methods, library preparation, and sequencing methods (Supplementary file 1a), suggests that our findings are robust with reference to technical variation.
Also, we did not have access to urinary oxalate or ascertain kidney stone history in our IBD cohort, both critical variables to extend our conclusions. However, previous studies demonstrated intestinal oxalate is predominantly absorbed paracellularly (Binder, 1974; Knauf et al., 2011; Saunders et al., 1975), and oxalate absorption is determined by concentration gradient, gut permeability, and oxalate bioavailability. Thus, colonic oxalate levels should correlate with both renal oxalate and with risk for CaOx stones. In addition to the untargeted metabolomics data used in this study, targeted measurement of fecal and urinary oxalate is desired in future controlled human studies. Oxalate transport in the human intestine has not been completely elucidated, nor has the relative importance of the small and large intestine to oxalate absorption been affirmatively determined. Therefore, colonization site and biogeography may be important factors for oxalate-degrading microbes to reduce host oxalate absorption. The ODP-transcribing microbes identified in the present study are based on fecal samples, which may bias towards organisms colonizing the colon, but not the small intestine. Furthermore, we showed in mice that O. formigenes colonization significantly affects oxalate homeostasis in vivo but further colonization with synthetic communities of different oxalate degraders is needed to identify which oxalate-degrading species have the largest impact on oxalate degradation and overall oxalate homeostasis.
In prior human studies, O. formigenes alone can effectively decrease host fecal or urinary oxalate levels, when host are exposed to high oxalate levels (Canales and Hatch, 2017; Jiang et al., 2011; Li et al., 2016; Li et al., 2015). In our study, O. formigenes colonization was below the level of detection in the jejunum and ileum; therefore, it is likely that O. formigenes predominately colonizes and performs oxalate degradation in the colon, reducing dietary oxalate being absorbed into the host. A study with germ-free mice (Li et al., 2016) observed that colonization of O. formigenes strain OXCC13 decreased the mouse fecal oxalate, which is consistent with our findings, but not urinary oxalate. However, several differences were noted between their study and ours, most importantly the use of different mouse strain and their use of a lower oxalate and higher calcium diet (0.1% oxalate and 1% calcium). In our study, both O. formigenes strains we tested reduced host urinary and fecal oxalate, although PH1 that was isolated from human showed a greater effect. This difference suggests that establishing the metabolic activity and host adaptability is critical for evaluating the therapeutic potential of individual O. formigenes strains.
Materials and methods
Meta-omics data of the human microbiome
Request a detailed protocolMetagenomic and metatranscriptomic data of healthy human subjects were collected from five and four studies, respectively (Abu-Ali et al., 2018; Ehrlich, 2011; Le Chatelier et al., 2013 Franzosa et al., 2014; Lloyd-Price et al., 2019; Petersen et al., 2017; Schirmer et al., 2018). Metagenomic and metatranscriptomic data of healthy humans and IBD subjects were collected from the iHMP-IBD study (Lloyd-Price et al., 2019). Each sample was cleaned by KneaData to remove low-quality reads and host-associated reads. The metabolic profiles of each sample were surveyed using HUMAnN2 v0.11.1 (Franzosa et al., 2018) under parameters --prescreen-threshold 0.01, --pathways-database metacyc_reactions_level4, metacyc_pathways_structured, and --protein-database uniref50, for the comparison in Figure 2—figure supplement 1. Fecal oxalate relative abundance was determined from untargeted metabolomics data from iHMP-IBD; measurements related to oxalate were selected and normalized against total metabolite (percent abundance of all observed metabolites) for analysis.
Homologous proteins of ODE
Request a detailed protocolThe protein homolog families of OXDD, FRC, and OXC were characterized by UniProt Interpro (Mitchell et al., 2019; Mulder et al., 2005) (V70) in protein families IPR017774, IPR017659, and IPR017660, respectively. We acquired the taxonomic origin and amino acid sequences of 2699 OXDD, 1947 FRC, and 1284 OXC homologs. Protein homologs that are 100% identical were then removed, leaving 2519 OXDD, 1556 FRC, and 1037 OXC unique homologs, which were used as a reference database of ODEs for subsequent query against the meta’omics data. Oxalate oxidoreductase (Figure 1A), a recently discovered enzyme for which there only is limited information (Anand et al., 2002; Dumas et al., 1993; Grąz et al., 2016; Kumar et al., 2011; Svedruzić et al., 2007; Tanner et al., 2001), was not considered in this present study.
Pairwise identity between ODE protein homologs
Request a detailed protocolMultiple sequence alignments were performed among the unique protein homologs separately, by muscle (Edgar, 2004) in seaview v4.7 (Gouy et al., 2010), and alignments were trimmed and imported into R. The pairwise alignment distance was calculated using function dist.alignment in the seqinR package (Charif et al., 2005) based on identity or Fitch matrix (Fitch, 1966). The alignment distance was subsequently converted to percent protein identity , following the documentation of dist.alignment.
Detection of ODE in the meta-omics data
Request a detailed protocolWe used two different approaches with complementary methodologies to ensure the conclusions are robust.
The first method we used is an agnostic approach that uses full protein sequences with high sensitivity. The quality-filtered meta’omics data were aligned against the reference protein databases consisting of the unique ODE homologs, by diamond blastx (Buchfink et al., 2015), with best hit returned (--max-target-seqs 1). Alignments with identity <90% were arbitrarily filtered out to prevent non-ODE reads from aligning to the ODE-specific reference proteins due to local similarity (Figure 2—figure supplement 2) By setting a stringent identity cutoff of 90%, we retain only the high-confidence alignment pairs. The abundance of each ODE protein homolog was calculated as reads per kilobase per million (RPKM) in each sample. When multiple timepoints were available, each subject was represented by the mean measurements across all samples provided.
We also used ShortBRED, which compares ODEs with all other known proteins to identify highly specific marker peptides, and thus can achieve high specificity. Using ShortBRED-identify, FRC, OXC, and OXDD were clustered into 202, 190, and 846 families respectively, with a centroid/representative protein selected for each family (Supplementary file 2). Then those centroids were compared against the uniref90 reference protein database (the ODE homologs were excluded), to remove the common region for identification of short peptide markers (Supplementary file 3). Lastly, filtered meta-omics reads were mapped against those peptide markers using ShortBRED-quantify with parameters --pctlength 0.5 and --id 0.9, to calculate the RPKM for each protein family. Default parameters were used unless noted otherwise.
Population-level contribution to ODE
Request a detailed protocolThe population-level contribution of a species to ODE was designed as a measurement to take both prevalence and abundance information into consideration. It is calculated for each ODE separately, based on their abundances (RPKM values). Using oxc as the example, suppose there are oxc-coding species and samples. In any given sample , the contribution of species to OXC, , is represented by its relative oxc abundance, calculated from
where the denotes the RPKMoxc of species in sample . In this way, we normalize across samples with the total contribution in any OXC-positive samples fixed to 1, and to 0 in any OXC-negative samples.
The population-level contribution of species : , can be subsequently calculated from summating contribution of species in samples, as follows
Note that population-level contribution of species monotonically increases with sample size . Therefore, it has been transformed to relative scale when being compared across different populations or different sample types (metagenome vs. metatranscriptome), such as in Figure 3C, Figure 3—figure supplement 4 and 12.
BURRITO for linking function to taxa
Request a detailed protocolBURRITO (McNally et al., 2018) was used to deconvolve ODE genes into taxa at the species level. Taxonomic profiles of healthy metagenome (n = 2539 samples) were generated by Metaphlan2 (Truong et al., 2015) under default parameters. In each sample, RPKM of all FRC or OXC homologs were summed as total abundance of frc or oxc genes and then supplied to BURRITO. The genomic content file was derived from taxonomic annotation of protein homologs from UniProt Interpro. We assumed one copy of frc and oxc for each genome to minimize bias, as such information is not available for all species.
Network analysis
Request a detailed protocolA network analysis of oxc or frc expression from microbial species used SpiecEasi (Kurtz et al., 2015). The raw RPKM values were used, and networks were constructed under default parameters method='mb', sel.criterion='bstars', lambda.min.ratio = 2e-2, nlambda = 100, and pulsar.params=list(rep.num = 20, ncores = 2).
Code availability
Request a detailed protocolSource code of the pipeline can be found on Github via https://github.com/ml3958/FindTaxaCtrbt (Liu, 2021; copy archived at swh:1:rev:13bbc4662f458bff327e348162bf51d875ed34d3). Downstream analysis scripts are available per request.
O. formigenes isolation and culturing
Request a detailed protocolO. formigenes strain PH1 was isolated from the stool sample of a primary hyperoxaluria patient using methods described in Allison et al., 1985. Both strains were cultured in defined oxalate broth as described (Allison et al., 1985), except that the concentration of oxalate is 50 mM. O. formigenes was cultured at 37°C in an anaerobic chamber before each mouse gavage.
O. formigenes colonization of mice
Request a detailed protocolA total of 17 C57Bl/6 male mice maintained in our gnotobiotic facility were used. The median of mouse baseline microbiome 16S qPCR measurements was 2.6 × 105, which is significantly lower than that of SPF mice (usually 108 or 109 in our experiments, representing <0.1% of the absolute abundance). The taxa detected in the baseline samples (Supplementary file 1e) did not lead to colonization resistance to the introduced O. formigenes strains.
At days 0, 2, and 7, mice were gavaged with 100 µl from a 24 hr growing culture of O. formigenes strain OXCC13 (n = 5), O. formigenes freshly isolated from a PH1 subject (n = 6), and O. formigenes culture medium alone (n = 6). Mice were fed normal chow supplemented with 1% sodium oxalate from study commencement (day −5) until mouse sacrifice at day 35. We did not perform an explicit power analysis since there is no previous data on our animal model and experimental conditions. We thought that five to six mice per group is a reasonable number of mice per group to show a significant change in urinary oxalate. Biological replicates are parallel measurements of biologically distinct samples that capture random biological variation, which may itself be a subject of study or a noise source. In each mouse group, we performed five to six biological replicates by gavaging material (culture media, growing culture of strain OXCC13 or PH1) into five to six mice.
DNA extraction and O. formigenes quantitative PCR
Request a detailed protocolMurine fecal pellets or intestinal contents were collected and stored at −80°C until DNA was extracted, using the MoBio 96-well extraction kit, following the manufacturer’s instructions. We confirmed that mice were colonized with O. formigenes at sacrifice using qPCR of fecal samples (Figure S16). qPCR was used to quantitate the number of copies of the oxc mRNA using the LightCycler 480 SYBR Green I Master Mix and run using the LightCycler 480 system. Paired primers (forward 5′-TGT-TTG-TCG-GCA-TTC-CTA-TC-3′, reverse 5′-TTG-GGA-AGC-AGT-TGG-TGG-3′) were used under the conditions: 95°C for 10 min, followed by 40 cycles of: 95°C for 23 s, 63°C for 20 s, 70°C for 40 s, and a final 30 s at 40°C as described (Pebenito et al., 2019).
Fecal and urinary collections and oxalate measurements
Request a detailed protocolMice were housed (n = 1 or 2) in metabolic cages, and 48 hr collections (1–2 mice per urine pool) were made under mineral oil into vessels containing crystal thymol as a preservative. Urinary oxalate (mg/dl) and creatinine (mg/dl) concentrations were determined in acidified (HCl) samples collected from all mice over a 48 hr period by Litholink Corp (Chicago, IL). Fecal pellets were collected at the end of the urine collections. Fecal pellets were acidified using 2M HCl, vortexed for 20 min, and then centrifuged at 21,000 g at room temperature, and supernatant fecal water collected using described methods (Jiang et al., 2011). Fecal water oxalate was measured using an oxalate calorimetric assay (Abcam, ab196990, Cambridge, UK) per the manufacturer’s instructions. Technical replicates are repeated measurements of the same sample that represent independent measures of the random noise associated with protocols or equipment. Fecal oxalate and qPCR were measured in duplicates.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files. Source data files have been provided for Figures 2–5.
-
NCBI BioProjectID PRJNA305507. Athlete Microbiome Project (AMP).
-
EBI European Nucleotide ArchiveID PRJEB4336. MetaHIT.
References
-
Oxalate-degrading Bacteria of the human gut as probiotics in the management of kidney stone diseaseAdvances in Applied Microbiology 72:63–87.https://doi.org/10.1016/S0065-2164(10)72003-7
-
Changes in ruminal oxalate degradation rates associated with adaptation to oxalate ingestionJournal of Animal Science 45:1173–1179.https://doi.org/10.2527/jas1977.4551173x
-
Oxalate degradation by gastrointestinal Bacteria from humansThe Journal of Nutrition 116:455–460.https://doi.org/10.1093/jn/116.3.455
-
Oxalobacter formigenes-Derived Bioactive Factors Stimulate Oxalate Transport by Intestinal Epithelial CellsJournal of the American Society of Nephrology 28:876–887.https://doi.org/10.1681/ASN.2016020132
-
Supersaturation and stone composition in a network of dispersed treatment sitesJournal of Urology 159:1821–1825.https://doi.org/10.1016/S0022-5347(01)63164-1
-
Transcriptional and functional analysis of oxalyl-coenzyme A (CoA) decarboxylase and formyl-CoA transferase genes from Lactobacillus acidophilusApplied and Environmental Microbiology 72:1891–1899.https://doi.org/10.1128/AEM.72.3.1891-1899.2006
-
The metabolism of oxalic acid in the animal bodyBiochemical Journal 34:144–148.https://doi.org/10.1042/bj0340144
-
The presence of Oxalobacter formigenes in the microbiome of healthy young adultsThe Journal of Urology 195:499–506.https://doi.org/10.1016/j.juro.2015.08.070
-
Hyperoxaluria and systemic oxalosis: an update on current therapy and future directionsExpert Opinion on Investigational Drugs 22:117–129.https://doi.org/10.1517/13543784.2013.741587
-
Intestinal oxalate absorptionGastroenterology 67:441–446.https://doi.org/10.1016/S0016-5085(19)32846-X
-
Microbial growth on oxalate by a route not involving glyoxylate carboligaseBiochemical Journal 118:53–59.https://doi.org/10.1042/bj1180053
-
Fast and sensitive protein alignment using DIAMONDNature Methods 12:59–60.https://doi.org/10.1038/nmeth.3176
-
Oxalobacter formigenes colonization normalizes oxalate excretion in a gastric bypass model of hyperoxaluriaSurgery for Obesity and Related Diseases 13:1152–1157.https://doi.org/10.1016/j.soard.2017.03.014
-
Assessment of in vitro oxalate degradation by Lactobacillus species cultured from veterinary probioticsAmerican Journal of Veterinary Research 76:801–806.https://doi.org/10.2460/ajvr.76.9.801
-
Renal involvement in inflammatory bowel diseasesJournal of Crohn's and Colitis 10:226–235.https://doi.org/10.1093/ecco-jcc/jjv138
-
Assimilation of oxalate, acetate, and CO2 by Oxalobacter formigenesCanadian Journal of Microbiology 42:1081–1086.https://doi.org/10.1139/m96-138
-
Nephrolithiasis in patients with inflammatory bowel disease in the communityInternational Journal of Nephrology and Renovascular Disease 6:139–142.https://doi.org/10.2147/IJNRD.S45466
-
Validation and clinical significance of a new calprotectin rapid test for the diagnosis of gastrointestinal diseasesInternational Journal of Colorectal Disease 23:985–992.https://doi.org/10.1007/s00384-008-0506-0
-
Oxalate metabolism by the acetogenic bacterium Moorella thermoaceticaFEMS Microbiology Letters 231:39–43.https://doi.org/10.1016/S0378-1097(03)00924-8
-
Identification of barley oxalate oxidase as a germin-like proteinComptes Rendus De l'Academie Des Sciences. Serie III, Sciences De La Vie 316:793–798.
-
Oxalobacter formigenes and its potential role in human healthApplied and Environmental Microbiology 68:3841–3847.https://doi.org/10.1128/AEM.68.8.3841-3847.2002
-
BookMetaHIT: The European Union Project on Metagenomics of the Human Intestinal TractIn: Nelson K, editors. Metagenomics of the Human Body. Springer. pp. 307–316.https://doi.org/10.1007/978-1-4419-7089-3_15
-
Proteome dynamics of the specialist oxalate degrader Oxalobacter formigenesJournal of Proteomics & Bioinformatics 9:19–24.https://doi.org/10.4172/jpb.1000384
-
Characterization and heterologous expression of the oxalyl coenzyme A decarboxylase gene from Bifidobacterium lactisApplied and Environmental Microbiology 70:5066–5073.https://doi.org/10.1128/AEM.70.9.5066-5073.2004
-
An improved method of testing for evolutionary homologyJournal of Molecular Biology 16:9–16.https://doi.org/10.1016/S0022-2836(66)80258-9
-
YfdW and YfdU are required for oxalate-induced acid tolerance in Escherichia coli K-12Journal of Bacteriology 195:1446–1455.https://doi.org/10.1128/JB.01936-12
-
SeaView version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree buildingMolecular Biology and Evolution 27:221–224.https://doi.org/10.1093/molbev/msp259
-
Oxalic acid degradation by a novel fungal oxalate oxidase from Abortiporus biennisActa Biochimica Polonica 63:1282.https://doi.org/10.18388/abp.2016_1282
-
Enteric oxalate elimination is induced and oxalate is normalized in a mouse model of primary hyperoxaluria following intestinal colonization with OxalobacterAmerican Journal of Physiology-Gastrointestinal and Liver Physiology 300:G461–G469.https://doi.org/10.1152/ajpgi.00434.2010
-
Efficacy and safety of Oxalobacter formigenes to reduce urinary oxalate in primary hyperoxaluriaNephrology Dialysis Transplantation 26:3609–3615.https://doi.org/10.1093/ndt/gfr107
-
Enteric hyperoxaluria: dependence on small intestinal resection, Colectomy, and steatorrhoea in chronic inflammatory bowel diseaseScandinavian Journal of Gastroenterology 13:577–588.https://doi.org/10.3109/00365527809181767
-
Urolithiasis and hyperoxaluria in chronic inflammatory bowel diseaseScandinavian Journal of Gastroenterology 14:475–479.
-
High-Specificity targeted functional profiling in microbial communities with ShortBREDPLOS Computational Biology 11:e1004557.https://doi.org/10.1371/journal.pcbi.1004557
-
Factors related to colonization with Oxalobacter formigenes in U.S. adultsJournal of Endourology 25:673–679.https://doi.org/10.1089/end.2010.0462
-
Net intestinal transport of oxalate reflects passive absorption and SLC26A6-mediated secretionJournal of the American Society of Nephrology 22:2247–2255.https://doi.org/10.1681/ASN.2011040433
-
A functional analysis of the formyl-coenzyme A (frc) gene from Lactobacillus reuteri 100-23CJournal of Applied Microbiology 116:1657–1667.https://doi.org/10.1111/jam.12500
-
Infrequency of colonization with Oxalobacter formigenes in inflammatory bowel disease: possible role in renal stone formationJournal of Gastroenterology and Hepatology 19:1403–1409.https://doi.org/10.1111/j.1440-1746.2004.03510.x
-
Sparse and compositionally robust inference of microbial ecological networksPLOS Computational Biology 11:e1004226.https://doi.org/10.1371/journal.pcbi.1004226
-
Oxalobacter formigenes colonization and oxalate dynamics in a mouse modelApplied and Environmental Microbiology 81:5048–5054.https://doi.org/10.1128/AEM.01313-15
-
Response of germ-free mice to colonization with O. formigenes and altered schaedler floraApplied and Environmental Microbiology 82:6952–6960.https://doi.org/10.1128/AEM.02381-16
-
Stone composition as a function of age and sexClinical Journal of the American Society of Nephrology 9:2141–2146.https://doi.org/10.2215/CJN.05660614
-
SoftwareFindTaxaCtrbt, version swh:1:rev:13bbc4662f458bff327e348162bf51d875ed34d3Software Heritage.
-
Enteric hyperoxaluria: role of Microbiota and antibioticsCurrent Opinion in Nephrology and Hypertension 28:352–359.https://doi.org/10.1097/MNH.0000000000000518
-
The gastrointestinal tract of the white-throated woodrat (Neotoma albigula) harbors distinct consortia of oxalate-degrading BacteriaApplied and Environmental Microbiology 80:1595–1601.https://doi.org/10.1128/AEM.03742-13
-
InterPro in 2019: improving coverage, classification and access to protein sequence annotationsNucleic Acids Research 47:D351–D360.https://doi.org/10.1093/nar/gky1100
-
Effect of antibiotics on Oxalobacter formigenes colonization of human gastrointestinal tractJournal of Endourology 19:102–106.https://doi.org/10.1089/end.2005.19.102
-
Screening of different probiotic strains for their in vitro ability to metabolise oxalates: any prospective use in humans?Journal of Clinical Gastroenterology 48 Suppl 1:S91–S95.https://doi.org/10.1097/MCG.0000000000000228
-
Calcium oxalate crystals induce renal inflammation by NLRP3-mediated IL-1β secretionJournal of Clinical Investigation 123:236–246.https://doi.org/10.1172/JCI63679
-
Cytotoxicity of crystals involves RIPK3-MLKL-mediated necroptosisNature Communications 7:10274.https://doi.org/10.1038/ncomms10274
-
Oxalate-induced chronic kidney disease with its uremic and cardiovascular complications in C57BL/6 miceAmerican Journal of Physiology-Renal Physiology 310:F785–F795.https://doi.org/10.1152/ajprenal.00488.2015
-
Hyperoxaluria requires TNF receptors to initiate crystal adhesion and kidney stone diseaseJournal of the American Society of Nephrology 28:761–768.https://doi.org/10.1681/ASN.2016040486
-
InterPro, progress and status in 2005Nucleic Acids Research 33:gki106.https://doi.org/10.1093/nar/gki106
-
Development of a humanized murine model for the study of Oxalobacter formigenes intestinal colonizationThe Journal of Infectious Diseases 220:1848–1858.https://doi.org/10.1093/infdis/jiz370
-
The ROKS nomogram for predicting a second symptomatic stone episodeJournal of the American Society of Nephrology 25:2878–2886.https://doi.org/10.1681/ASN.2013091011
-
US renal data system 2017 annual data report: epidemiology of kidney disease in the united statesAmerican Journal of Kidney Diseases 71:A7.https://doi.org/10.1053/j.ajkd.2018.01.002
-
Prevalence of kidney stones in the united statesEuropean Urology 62:160–165.https://doi.org/10.1016/j.eururo.2012.03.052
-
Time trends in reported prevalence of kidney stones in the United States: 1976-1994Kidney International 63:1817–1823.https://doi.org/10.1046/j.1523-1755.2003.00917.x
-
Investigating the roles of putative active site residues in the oxalate decarboxylase from Bacillus subtilisArchives of Biochemistry and Biophysics 464:36–47.https://doi.org/10.1016/j.abb.2007.03.016
-
Oxalate decarboxylase requires manganese and dioxygen for activity. Overexpression and characterization of Bacillus subtilis YvrK and YoaNThe Journal of Biological Chemistry 276:43627–43634.https://doi.org/10.1074/jbc.M107202200
-
MetaPhlAn2 for enhanced metagenomic taxonomic profilingNature Methods 12:902–903.https://doi.org/10.1038/nmeth.3589
Article and author information
Author details
Funding
National Institute of Allergy and Infectious Diseases (NIAID) (U01AI22285)
- Martin J Blaser
National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK) (R01DK110014)
- Huilin Li
Rare Kidney Stone Consortium (U54 DK083908)
- Lama Nazzal
The C & D and Zlinkoff Funds
- Martin J Blaser
Oxalosis and Hyperoxaluria Foundation (career development grant)
- Lama Nazzal
TransAtlantic Partnership of the Fondation LeDucq
- Martin J Blaser
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Dr. David Goldfarb and Xuhui Zheng for their helpful comments. We thank Dr. Tim Borbet for the helpful discussion and input. Funding: This study was supported in part by U01AI22285, R01DK110014, and the Rare Kidney Stone Consortium (U54 DK083908) from the National Institutes of Health, by the C and D and Zlinkoff Funds, Oxalosis and Hyperoxaluria Foundation-American Society of Nephrology career development grant, and the TransAtlantic Partnership of the Fondation LeDucq.
Ethics
Animal experimentation: This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All of the animals were handled according to approved institutional animal care and use committee (IACUC) protocols (#IA16-00822) of the New York University Langone Medical Center.
Copyright
© 2021, Liu et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 3,224
- views
-
- 393
- downloads
-
- 40
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Neuroscience
The basolateral amygdala (BLA) is a key site where fear learning takes place through synaptic plasticity. Rodent research shows prominent low theta (~3–6 Hz), high theta (~6–12 Hz), and gamma (>30 Hz) rhythms in the BLA local field potential recordings. However, it is not understood what role these rhythms play in supporting the plasticity. Here, we create a biophysically detailed model of the BLA circuit to show that several classes of interneurons (PV, SOM, and VIP) in the BLA can be critically involved in producing the rhythms; these rhythms promote the formation of a dedicated fear circuit shaped through spike-timing-dependent plasticity. Each class of interneurons is necessary for the plasticity. We find that the low theta rhythm is a biomarker of successful fear conditioning. The model makes use of interneurons commonly found in the cortex and, hence, may apply to a wide variety of associative learning situations.
-
- Cancer Biology
- Computational and Systems Biology
Effects from aging in single cells are heterogenous, whereas at the organ- and tissue-levels aging phenotypes tend to appear as stereotypical changes. The mammary epithelium is a bilayer of two major phenotypically and functionally distinct cell lineages: luminal epithelial and myoepithelial cells. Mammary luminal epithelia exhibit substantial stereotypical changes with age that merit attention because these cells are the putative cells-of-origin for breast cancers. We hypothesize that effects from aging that impinge upon maintenance of lineage fidelity increase susceptibility to cancer initiation. We generated and analyzed transcriptomes from primary luminal epithelial and myoepithelial cells from younger <30 (y)ears old and older >55y women. In addition to age-dependent directional changes in gene expression, we observed increased transcriptional variance with age that contributed to genome-wide loss of lineage fidelity. Age-dependent variant responses were common to both lineages, whereas directional changes were almost exclusively detected in luminal epithelia and involved altered regulation of chromatin and genome organizers such as SATB1. Epithelial expression of gap junction protein GJB6 increased with age, and modulation of GJB6 expression in heterochronous co-cultures revealed that it provided a communication conduit from myoepithelial cells that drove directional change in luminal cells. Age-dependent luminal transcriptomes comprised a prominent signal that could be detected in bulk tissue during aging and transition into cancers. A machine learning classifier based on luminal-specific aging distinguished normal from cancer tissue and was highly predictive of breast cancer subtype. We speculate that luminal epithelia are the ultimate site of integration of the variant responses to aging in their surrounding tissue, and that their emergent phenotype both endows cells with the ability to become cancer-cells-of-origin and represents a biosensor that presages cancer susceptibility.