Differential substrate specificity and kinetic behavior of Escherichia coli YfdW and Oxalobacter formigenes formyl coenzyme A transferase

Engineered microorganisms: A new direction in kidney stone prevention and treatment

Numerous studies have shown that intestinal and urinary tract flora are closely related to the formation of kidney stones. The removal of probiotics represented by lactic acid bacteria and the colonization of pathogenic bacteria can directly or indirectly promote the occurrence of kidney stones. However, currently existing natural probiotics have limitations. Synthetic biology is an emerging discipline in which cells or living organisms are genetically designed and modified to have biological functions that meet human needs, or even create new biological systems, and has now become a research hotspot in various fields. Using synthetic biology approaches of microbial engineering and biological redesign to enable probiotic bacteria to acquire new phenotypes or heterologous protein expression capabilities is an important part of synthetic biology research. Synthetic biology modification of microorganisms in the gut and urinary tract can effectively inhibit the development of kidney stones by a range of means, including direct degradation of metabolites that promote stone production or indirect regulation of flora homeostasis. This article reviews the research status of engineered microorganisms in the prevention and treatment of kidney stones, to provide a new and effective idea for the prevention and treatment of kidney stones.

Acyltransferase families that act on thioesters: Sequences, structures, and mechanisms

Acyltransferases (AT) are enzymes that catalyze the transfer of acyl group to a receptor molecule. This review focuses on ATs that act on thioester-containing substrates. Although many ATs can recognize a wide variety of substrates, sequence similarity analysis allowed us to classify the ATs into fifteen distinct families. Each AT family is originated from enzymes experimentally characterized to have AT activity, classified according to sequence similarity, and confirmed with tertiary structure similarity for families that have crystallized structures available. All the sequences and structures of the AT families described here are present in the thioester-active enzyme (ThYme) database. The AT sequences and structures classified into families and available in the ThYme database could contribute to enlightening the understanding acyl transfer to thioester-containing substrates, most commonly coenzyme A, which occur in multiple metabolic pathways, mostly with fatty acids.

Metabolic reprogramming and biosensor-assisted mutagenesis screening for high-level production of L-arginine in Escherichia coli

L-arginine is a value-added amino acid with promising applications in the pharmaceutical and nutraceutical industries. Further unleashing the potential of microbial cell factories to make L-arginine production more competitive remains challenging due to the sophisticated intracellular interaction networks and the insufficient knowledge of global metabolic regulation. Here, we combined multilevel rational metabolic engineering with biosensor-assisted mutagenesis screening to exploit the L-arginine production potential of Escherichia coli. First, multiple metabolic pathways were systematically reprogrammed to redirect the metabolic flux into L-arginine synthesis, including the L-arginine biosynthesis, TCA cycle, and L-arginine export. Specifically, a toggle switch responding to special cellular physiological conditions was designed to dynamically control the expression of sucA and pull more carbon flux from the TCA cycle toward L-arginine biosynthesis. Subsequently, a biosensor-assisted high-throughput screening platform was designed and applied to further exploit the L-arginine production potential. The best-engineered ARG28 strain produced 132 g/L L-arginine in a 5-L bioreactor with a yield of 0.51 g/g glucose and productivity of 2.75 g/(L ⋅ h), which were the highest values reported so far. Through whole genome sequencing and reverse engineering, Frc frameshift mutant, PqiB A78P mutant, and RpoB P564T mutant were revealed for enhancing the L-arginine biosynthesis. Our study exhibited the power of coupling rational metabolic reprogramming and biosensor-assisted mutagenesis screening to unleash the cellular potential for value-added metabolite production.

Expanding the taxonomic and environmental extent of an underexplored carbon metabolism-oxalotrophy

Oxalate serves various functions in the biological processes of plants, fungi, bacteria, and animals. It occurs naturally in the minerals weddellite and whewellite (calcium oxalates) or as oxalic acid. The environmental accumulation of oxalate is disproportionately low compared to the prevalence of highly productive oxalogens, namely plants. It is hypothesized that oxalotrophic microbes limit oxalate accumulation by degrading oxalate minerals to carbonates via an under-explored biogeochemical cycle known as the oxalate-carbonate pathway (OCP). Neither the diversity nor the ecology of oxalotrophic bacteria is fully understood. This research investigated the phylogenetic relationships of the bacterial genes oxc, frc, oxdC, and oxlT, which encode key enzymes for oxalotrophy, using bioinformatic approaches and publicly available omics datasets. Phylogenetic trees of oxc and oxdC genes demonstrated grouping by both source environment and taxonomy. All four trees included genes from metagenome-assembled genomes (MAGs) that contained novel lineages and environments for oxalotrophs. In particular, sequences of each gene were recovered from marine environments. These results were supported with marine transcriptome sequences and description of key amino acid residue conservation. Additionally, we investigated the theoretical energy yield from oxalotrophy across marine-relevant pressure and temperature conditions and found similar standard state Gibbs free energy to "low energy" marine sediment metabolisms, such as anaerobic oxidation of methane coupled to sulfate reduction. These findings suggest further need to understand the role of bacterial oxalotrophy in the OCP, particularly in marine environments, and its contribution to global carbon cycling.

Crystal Structure of an Intramolecular Mesaconyl-Coenzyme A Transferase From the 3-Hydroxypropionic Acid Cycle of Roseiflexus castenholzii

Coenzyme A (CoA) transferases catalyze reversible transfer of CoA groups from CoA-thioesters to free acids, playing important roles in the metabolism of carboxylic acids in all organisms. An intramolecular CoA transferase, Mesaconyl-CoA C1-C4 CoA transferase (MCT) was identified in the autotrophic CO2 fixation pathway, 3-hydroxypropionic acid cycle of filamentous anoxygenic phototrophs (FAPs). Different from the well-known CoA transferases that catalyze CoA transfer between two distinct substrates, MCT specifically catalyzes the reversible transformation of mesaconyl-C1-CoA to mesaconyl-C4-CoA, a key reaction intermediate for carbon fixation. However, the molecular mechanism of MCT in employing one substrate is enigmatic. Here we determined the crystal structure of MCT from a chlorosome-less FAP Roseiflexus castenholzii at 2.5 Å resolution, and characterized the catalytic mechanisms through structural analyses and molecular dynamic simulations. The structure of R. castenholzii MCT consists of a Rossmann fold larger domain and a small domain that are connected by two linkers. Two MCT subunits are cross interlocked at the linker regions to form a functional dimer in solution, in which the substrate binding pockets are located at the interface of the Rossmann fold larger domain from one subunit and the small domain from the other subunit. In the simulated binding structures, both the substrate mesaconyl-C1-CoA and product mesaconyl-C4-CoA form extensive electrostatic and hydrogen bonding interactions with MCT. But some differences exist in the binding mode of these two CoA analogs, Arg314’ from the second subunit of the dimer presenting dramatic conformational changes in binding with mesaconyl-C4-CoA. Together with Arg47 and one water molecule, a strictly conserved residue Asp165 are essential for catalyzing the reversible intramolecular CoA transfer reaction, through the electrostatic and hydrogen bonding interactions with the mesaconic tail of both the substrate and product. This study revealed a previously unrecognized mechanism for the uncommon intramolecular CoA transfer reaction, which will not only broaden the knowledge on the catalytic mechanisms of CoA transferases, but also contribute to enzyme engineering or biosynthetic applications of the 3-HP cycle for synthesis of fine chemicals and important metabolites.

Redefining the coenzyme A transferase superfamily with a large set of manually-annotated proteins

The coenzyme A (CoA) transferases are a superfamily of proteins central to the metabolism of acetyl-CoA and other CoA thioesters. They are diverse group, catalyzing over a hundred biochemical reactions and spanning all three domains of life. A deeply rooted idea, proposed two decades ago, is these enzymes fall into three families (I, II, III). Here we find they fall into different families, which we achieve by analyzing all CoA transferases characterized to date. We manually annotated 94 CoA transferases with functional information (including rates of catalysis for 208 reactions) from 97 publications. This represents all enzymes we could find in the primary literature, and it is double the number annotated in four protein databases (BRENDA, KEGG, MetaCyc, UniProt). We found family I transferases are not closely related to each other in terms of sequence, structure, and reactions catalyzed. This family is not even monophyletic. These problems are solved by regrouping the three families into six, including one family with many non-CoA transferases. The problem (and solution) became apparent only by analyzing our large set of manually-annotated proteins. It would have been missed if we had used the small number of proteins annotated in UniProt and other databases. Our work is important to understanding the biology of CoA transferases. It also warns investigators doing phylogenetic analyses of proteins to go beyond information in databases. This article is protected by copyright. All rights reserved.

Probiotics in the Prevention of the Calcium Oxalate Urolithiasis

Nephrolithiasis ranks third among urological diseases in terms of prevalence, making up about 15% of cases. The continued increase in the incidence of nephrolithiasis is most probably due to changes in eating habits (high protein, sodium, and sugar diets) and lifestyle (reduced physical activity) in all developed countries. Some 80% of all kidney stones cases are oxalate urolithiasis, which is also characterized by the highest risk of recurrence. Frequent relapses of nephrolithiasis contribute to severe complications and high treatment costs. Unfortunately, there is no known effective way to prevent urolithiasis at present. In cases of diet-related urolithiasis, dietary changes may prevent recurrence. However, in some patients, the condition is unrelated to diet; in such cases, there is evidence to support the use of stone-related medications. Interestingly, a growing body of evidence indicates the potential of the microbiome to reduce the risk of developing renal colic. Previous studies have primarily focused on the use of Oxalobacterformigenes in patients with urolithiasis. Unfortunately, this bacterium is not an ideal probiotic due to its antibiotic sensitivity and low pH. Therefore, subsequent studies sought to find bacteria which are capable of oxalate degradation, focusing on well-known probiotics including Lactobacillus and Bifidobacterium strains, Eubacterium lentum, Enterococcus faecalis, and Escherichia coli.

Rv3272 encodes a novel Family III CoA transferase that alters the cell wall lipid profile and protects mycobacteria from acidic and oxidative stress

The availability of complete genome sequence of Mycobacterium tuberculosis has provided an important tool to understand the mycobacterial biology with respect to host-pathogen interaction, which is an unmet need of the hour owing to continuous increasing drug resistance. Hypothetical proteins are often an overlooked pool though half the genome encodes for such proteins of unknown function that could potentially play vital roles in mycobacterial biology. In this context, we report the structural and functional characterization of the hypothetical protein Rv3272. Sequence analysis classifies Rv3272 as a Family III CoA transferase with the classical two domain structure and conserved Aspartate residue (D175). The crystal structure of the wild type protein (2.2 Å) demonstrated the associated inter-locked dimer while that of the D175A mutant co-crystallized with octanoyl-CoA demonstrated relative movement between the two domains. Isothermal titration calorimetry studies indicate that Rv3272 binds to fatty acyl-CoAs of varying carbon chain lengths, with palmitoyl-CoA (C16:0) exhibiting maximum affinity. To determine the functional relevance of Rv3272 in mycobacterial biology, we ectopically expressed Rv3272 in M. smegmatis and assessed that its expression encodes significant alteration in cell surface with marked differences in triacylglycerol accumulation. Additionally, Rv3272 expression protects mycobacteria from acidic, oxidative and antibiotic stress under in vitro conditions. Taken together, these studies indicate a significant role for Rv3272 in host-pathogen interaction.

Characterization of mutations in the PAS domain of the EvgS sensor kinase selected by laboratory evolution for acid resistance in Escherichia coli

Laboratory-based evolution and whole-genome sequencing can link genotype and phenotype. We used evolution of acid resistance in exponential phase Escherichia coli to study resistance to a lethal stress. Iterative selection at pH 2.5 generated five populations that were resistant to low pH in early exponential phase. Genome sequencing revealed multiple mutations, but the only gene mutated in all strains was evgS, part of a two-component system that has already been implicated in acid resistance. All these mutations were in the cytoplasmic PAS domain of EvgS, and were shown to be solely responsible for the resistant phenotype, causing strong upregulation at neutral pH of genes normally induced by low pH. Resistance to pH 2.5 in these strains did not require the transporter GadC, or the sigma factor RpoS. We found that EvgS-dependent constitutive acid resistance to pH 2.5 was retained in the absence of the regulators GadE or YdeO, but was lost if the oxidoreductase YdeP was also absent. A deletion in the periplasmic domain of EvgS abolished the response to low pH, but not the activity of the constitutive mutants. On the basis of these results we propose a model for how EvgS may become activated by low pH.

A functional analysis of the formyl-coenzyme A (frc) gene from Lactobacillus reuteri 100-23C

AIM

To examine the role of the Lactobacillus reuteri 100-23C frc gene product in oxalate metabolism, host colonization and the acid stress response.

METHODS AND RESULTS

Genes encoding putative formyl-CoA transferase (frc) and oxalyl-CoA decarboxylase (oxc) enzymes are present in the genome sequences of Lact. reuteri strains. Two strains isolated from humans harboured an IS200 insertion sequence in the frc ORF and a group 2 intron-associated transposase downstream of the frc gene, both of which were lacking in two strains of animal origin, which contained intact frc and oxc genes. An frc(-) insertional mutant of Lact. reuteri 100-23C was compared with the parent strain with respect to oxalate degradation, colonization of an RLF-mouse host model and growth in the presence of acids. Neither parent nor mutant degraded oxalate in vitro or in vivo. However, the parent outcompeted the frc(-) mutant in the mouse intestine during co-colonization and the frc(-) mutant showed a reduced growth rate in the presence of hydrochloric acid.

CONCLUSIONS

Intact oxc and frc genes do not ensure oxalate degradation under the conditions tested. The frc gene product is important during host colonization and survival of acid stress by Lact. reuteri 100-23C.

SIGNIFICANCE AND IMPACT OF THE STUDY

Oxalate metabolism by oxalate-degrading intestinal bacterial strains may be important in preventing urolithiasis and might lead to the derivation of probiotic products. To produce safe and efficacious probiotics, however, an understanding of the genetic characteristics of potential oxalate degraders must be obtained, together with knowledge of their functional ramifications.

Function and X-ray crystal structure of Escherichia coli YfdE

Many food plants accumulate oxalate, which humans absorb but do not metabolize, leading to the formation of urinary stones. The commensal bacterium Oxalobacter formigenes consumes oxalate by converting it to oxalyl-CoA, which is decarboxylated by oxalyl-CoA decarboxylase (OXC). OXC and the class III CoA-transferase formyl-CoA:oxalate CoA-transferase (FCOCT) are widespread among bacteria, including many that have no apparent ability to degrade or to resist external oxalate. The EvgA acid response regulator activates transcription of the Escherichia coli yfdXWUVE operon encoding YfdW (FCOCT), YfdU (OXC), and YfdE, a class III CoA-transferase that is ~30% identical to YfdW. YfdW and YfdU are necessary and sufficient for oxalate-induced protection against a subsequent acid challenge; neither of the other genes has a known function. We report the purification, in vitro characterization, 2.1-Å crystal structure, and functional assignment of YfdE. YfdE and UctC, an orthologue from the obligate aerobe Acetobacter aceti, perform the reversible conversion of acetyl-CoA and oxalate to oxalyl-CoA and acetate. The annotation of YfdE as acetyl-CoA:oxalate CoA-transferase (ACOCT) expands the scope of metabolic pathways linked to oxalate catabolism and the oxalate-induced acid tolerance response. FCOCT and ACOCT active sites contain distinctive, conserved active site loops (the glycine-rich loop and the GNxH loop, respectively) that appear to encode substrate specificity.

YfdW and YfdU are required for oxalate-induced acid tolerance in Escherichia coli K-12

Escherichia coli has several mechanisms for surviving low-pH stress. We report that oxalic acid, a small-chain organic acid (SCOA), induces a moderate acid tolerance response (ATR) in two ways. Adaptation of E. coli K-12 at pH 5.5 with 50 mM oxalate and inclusion of 25 mM oxalate in pH 3.0 minimal challenge medium separately conferred protection, with 67% ± 7% and 87% ± 17% survival after 2 h, respectively. The combination of oxalate adaptation and oxalate supplementation in the challenge medium resulted in increased survival over adaptation or oxalate in the challenge medium alone. The enzymes YfdW, a formyl coenzyme A (CoA) transferase, and YfdU, an oxalyl-CoA decarboxylase, are required for the adaptation effect but not during challenge. Unlike other SCOAs, this oxalate ATR is not a part of the RpoS regulon but appears to be linked to the signal protein GadE. We theorize that this oxalate ATR could enhance the pathogenesis of virulent E. coli consumed with oxalate-containing foods like spinach.

Formyl-coenzyme A (CoA):oxalate CoA-transferase from the acidophile Acetobacter aceti has a distinctive electrostatic surface and inherent acid stability

Bacterial formyl-CoA:oxalate CoA-transferase (FCOCT) and oxalyl-CoA decarboxylase work in tandem to perform a proton-consuming decarboxylation that has been suggested to have a role in generalized acid resistance. FCOCT is the product of uctB in the acidophilic acetic acid bacterium Acetobacter aceti. As expected for an acid-resistance factor, UctB remains folded at the low pH values encountered in the A. aceti cytoplasm. A comparison of crystal structures of FCOCTs and related proteins revealed few features in UctB that would distinguish it from nonacidophilic proteins and thereby account for its acid stability properties, other than a strikingly featureless electrostatic surface. The apparently neutral surface is a result of a "speckled" charge decoration, in which charged surface residues are surrounded by compensating charges but do not form salt bridges. A quantitative comparison among orthologs identified a pattern of residue substitution in UctB that may be a consequence of selection for protein stability by constant exposure to acetic acid. We suggest that this surface charge pattern, which is a distinctive feature of A. aceti proteins, creates a stabilizing electrostatic network without stiffening the protein or compromising protein-solvent interactions.

The evolution of metabolic networks of E. coli

Background

Despite the availability of numerous complete genome sequences from E. coli strains, published genome-scale metabolic models exist only for two commensal E. coli strains. These models have proven useful for many applications, such as engineering strains for desired product formation, and we sought to explore how constructing and evaluating additional metabolic models for E. coli strains could enhance these efforts.

Results

We used the genomic information from 16 E. coli strains to generate an E. coli pangenome metabolic network by evaluating their collective 76,990 ORFs. Each of these ORFs was assigned to one of 17,647 ortholog groups including ORFs associated with reactions in the most recent metabolic model for E. coli K-12. For orthologous groups that contain an ORF already represented in the MG1655 model, the gene to protein to reaction associations represented in this model could then be easily propagated to other E. coli strain models. All remaining orthologous groups were evaluated to see if new metabolic reactions could be added to generate a pangenome-scale metabolic model (iEco1712_pan). The pangenome model included reactions from a metabolic model update for E. coli K-12 MG1655 (iEco1339_MG1655) and enabled development of five additional strain-specific genome-scale metabolic models. These additional models include a second K-12 strain (iEco1335_W3110) and four pathogenic strains (two enterohemorrhagic E. coli O157:H7 and two uropathogens). When compared to the E. coli K-12 models, the metabolic models for the enterohemorrhagic (iEco1344_EDL933 and iEco1345_Sakai) and uropathogenic strains (iEco1288_CFT073 and iEco1301_UTI89) contained numerous lineage-specific gene and reaction differences. All six E. coli models were evaluated by comparing model predictions to carbon source utilization measurements under aerobic and anaerobic conditions, and to batch growth profiles in minimal media with 0.2% (w/v) glucose. An ancestral genome-scale metabolic model based on conserved ortholog groups in all 16 E. coli genomes was also constructed, reflecting the conserved ancestral core of E. coli metabolism (iEco1053_core). Comparative analysis of all six strain-specific E. coli models revealed that some of the pathogenic E. coli strains possess reactions in their metabolic networks enabling higher biomass yields on glucose. Finally the lineage-specific metabolic traits were compared to the ancestral core model predictions to derive new insight into the evolution of metabolism within this species.

Conclusion

Our findings demonstrate that a pangenome-scale metabolic model can be used to rapidly construct additional E. coli strain-specific models, and that quantitative models of different strains of E. coli can accurately predict strain-specific phenotypes. Such pangenome and strain-specific models can be further used to engineer metabolic phenotypes of interest, such as designing new industrial E. coli strains.

Identification and characterization of two bile acid coenzyme A transferases from Clostridium scindens, a bile acid 7α-dehydroxylating intestinal bacterium

The human bile acid pool composition is composed of both primary bile acids (cholic acid and chenodeoxycholic acid) and secondary bile acids (deoxycholic acid and lithocholic acid). Secondary bile acids are formed by the 7α-dehydroxylation of primary bile acids carried out by intestinal anaerobic bacteria. We have previously described a multistep biochemical pathway in Clostridium scindens that is responsible for bile acid 7α-dehydroxylation. We have identified a large (12 kb) bile acid inducible (bai) operon in this bacterium that encodes eight genes involved in bile acid 7α-dehydroxylation. However, the function of the baiF gene product in this operon has not been elucidated. In the current study, we cloned and expressed the baiF gene in E. coli and discovered it has bile acid CoA transferase activity. In addition, we discovered a second bai operon encoding three genes. The baiK gene in this operon was expressed in E. coli and found to encode a second bile acid CoA transferase. Both bile acid CoA transferases were determined to be members of the type III family by amino acid sequence comparisons. Both bile acid CoA transferases had broad substrate specificity, except the baiK gene product, which failed to use lithocholyl-CoA as a CoA donor. Primary bile acids are ligated to CoA via an ATP-dependent mechanism during the initial steps of 7α-dehydroxylation. The bile acid CoA transferases conserve the thioester bond energy, saving the cell ATP molecules during bile acid 7α-dehydroxylation. ATP-dependent CoA ligation is likely quickly supplanted by ATP-independent CoA transfer.

Use of an isothermal microcalorimetry assay to characterize microbial oxalotrophic activity

Isothermal microcalorimetry (IMC) has been used in the past to monitor metabolic activities in living systems. A few studies have used it on ecological research. In this study, IMC was used to monitor oxalotrophic activity, a widespread bacterial metabolism found in the environment, and particularly in soils. Six model strains were inoculated in solid angle media with K-oxalate as the sole carbon source. Cupriavidus oxalaticus, Cupriavidus necator, and Streptomyces violaceoruber presented the highest activity (91, 40, and 55 μW, respectively) and a maximum growth rate (μmax h(-1) ) of 0.264, 0.185, and 0.199, respectively, among the strains tested. These three strains were selected to test the incidence of different oxalate sources (Ca, Cu, and Fe-oxalate salts) in the metabolic activity. The highest activity was obtained in Ca-oxalate for C. oxalaticus. Similar experiments were carried out with a model soil to test whether this approach can be used to measure oxalotrophic activity in field samples. Although measuring oxalotrophic activity in a soil was challenging, there was a clear effect of the amendment with oxalate on the metabolic activity measured in soil. The correlation between heat flow and growth suggests that IMC analysis is a powerful method to monitor bacterial oxalotrophic activity.

Sampling long time scale protein motions: OSRW simulation of active site loop conformational free energies in formyl-CoA:oxalate CoA transferase

X-ray crystallographic snapshots have shown that conformational changes of a tetraglycine loop in the active site of formyl-CoA:oxalate CoA transferase (FRC) play an important role in the catalytic cycle of the enzyme. Orthogonal space random walk (OSRW) simulations have been applied to obtain quantitative computational estimates of the relative free energy of the "open" and "closed" conformations of this loop together with the energetic barrier for interconversion of these states in wild type FRC. These OSRW calculations not only show that the two conformations have similar free energies but also predict a barrier that is consistent with the observed turnover number of the enzyme. In an effort to quantitate the importance of specific residues in the tetraglycine loop, OSRW simulations have also been performed on the G258A, G259A, G260A, and G261A FRC variants both to examine the energetic effects of replacing each glycine residue and to correlate the computed energies with kinetic and structural observations. In enzymes with substantially reduced catalytic efficiency (k(cat)/K(M)), the OSRW simulations reveal the adoption of additional low energy loop conformations. In the case of the G260A FRC variant, the new conformation identified by simulation is similar to that observed in the X-ray crystal structure of the protein. These results provide further evidence for the power of the OSRW method in sampling conformational space and, hence, in providing quantitative free energy estimates for the conformations adopted by functionally important active site loops. In addition, these simulations model the motions of side chains that are correlated with changes in loop conformation thereby permitting access to long time-scale motions through the use of nanosecond simulations.

Use of the frc gene as a molecular marker to characterize oxalate-oxidizing bacterial abundance and diversity structure in soil

Oxalate catabolism, which can have both medical and environmental implications, is performed by phylogenetically diverse bacteria. The formyl-CoA-transferase gene was chosen as a molecular marker of the oxalotrophic function. Degenerated primers were deduced from an alignment of frc gene sequences available in databases. The specificity of primers was tested on a variety of frc-containing and frc-lacking bacteria. The frc-primers were then used to develop PCR-DGGE and real-time SybrGreen PCR assays in soils containing various amounts of oxalate. Some PCR products from pure cultures and from soil samples were cloned and sequenced. Data were used to generate a phylogenetic tree showing that environmental PCR products belonged to the target physiological group. The extent of diversity visualised on DGGE pattern was higher for soil samples containing carbonate resulting from oxalate catabolism. Moreover, the amount of frc gene copies in the investigated soils was detected in the range of 1.64x10(7) to 1.75x10(8)/g of dry soil under oxalogenic tree (representing 0.5 to 1.2% of total 16S rRNA gene copies), whereas the number of frc gene copies in the reference soil was 6.4x10(6) (or 0.2% of 16S rRNA gene copies). This indicates that oxalotrophic bacteria are numerous and widespread in soils and that a relationship exists between the presence of the oxalogenic trees Milicia excelsa and Afzelia africana and the relative abundance of oxalotrophic guilds in the total bacterial communities. This is obviously related to the accomplishment of the oxalate-carbonate pathway, which explains the alkalinization and calcium carbonate accumulation occurring below these trees in an otherwise acidic soil. The molecular tools developed in this study will allow in-depth understanding of the functional implication of these bacteria on carbonate accumulation as a way of atmospheric CO(2) sequestration.