Characterizing the metabolism of Dehalococcoides with a constraint-based model

Taxonomic Re-Classification and Expansion of the Phylum Chloroflexota Based on over 5000 Genomes and Metagenome-Assembled Genomes

The phylum Chloroflexota (formerly Chloroflexi) encompasses metabolically diverse bacteria that often have high prevalence in terrestrial and aquatic habitats, some even with biotechnological application. However, there is substantial disagreement in public databases which lineage should be considered a member of the phylum and at what taxonomic level. Here, we addressed these issues through extensive phylogenomic analyses. The analyses were based on a collection of >5000 Chloroflexota genomes and metagenome-assembled genomes (MAGs) from public databases, novel environmental sites, as well as newly generated MAGs from publicly available sequence reads via an improved binning approach incorporating covariance information. Based on calculated relative evolutionary divergence, we propose that Candidatus Dormibacterota should be listed as a class (i.e., Ca. Dormibacteria) within Chloroflexota together with the classes Anaerolineae, Chloroflexia, Dehalococcoidia, Ktedonobacteria, Ca. Limnocylindria, Thermomicrobia, and two other classes containing only uncultured members. All other Chloroflexota lineages previously listed at the class rank appear to be rather orders or families in the Anaerolineae and Dehalococcoidia, which contain the vast majority of genomes and exhibited the strongest phylogenetic radiation within the phylum. Furthermore, the study suggests that a common ecophysiological capability of members of the phylum is to successfully cope with low energy fluxes.

Metagenomic Analysis Reveals Microbial Interactions at the Biocathode of a Bioelectrochemical System Capable of Simultaneous Trichloroethylene and Cr(VI) Reduction

Bioelectrochemical systems (BES) are attractive and versatile options for the bioremediation of organic or inorganic pollutants, including trichloroethylene (TCE) and Cr(VI), often found as co-contaminants in the environment. The elucidation of the microbial players' role in the bioelectroremediation processes for treating multicontaminated groundwater is still a research need that attracts scientific interest. In this study, 16S rRNA gene amplicon sequencing and whole shotgun metagenomics revealed the leading microbial players and the primary metabolic interactions occurring in the biofilm growing at the biocathode where TCE reductive dechlorination (RD), hydrogenotrophic methanogenesis, and Cr(VI) reduction occurred. The presence of Cr(VI) did not negatively affect the TCE degradation, as evidenced by the RD rates estimated during the reactor operation with TCE (111±2 μeq/Ld) and TCE/Cr(VI) (146±2 μeq/Ld). Accordingly, Dehalococcoides mccartyi, the primary biomarker of the RD process, was found on the biocathode treating both TCE (7.82E+04±2.9E+04 16S rRNA gene copies g-1 graphite) and TCE/Cr(VI) (3.2E+07±2.37E+0716S rRNA gene copies g-1 graphite) contamination. The metagenomic analysis revealed a selected microbial consortium on the TCE/Cr(VI) biocathode. D. mccartyi was the sole dechlorinating microbe with H2 uptake as the only electron supply mechanism, suggesting that electroactivity is not a property of this microorganism. Methanobrevibacter arboriphilus and Methanobacterium formicicum also colonized the biocathode as H2 consumers for the CH4 production and cofactor suppliers for D. mccartyi cobalamin biosynthesis. Interestingly, M. formicicum also harbors gene complexes involved in the Cr(VI) reduction through extracellular and intracellular mechanisms.

Impact of iron- and/or sulfate-reduction on a cis-1,2-dichloroethene and vinyl chloride respiring bacterial consortium: experiments and model-based interpretation

Process understanding of microbial communities containing organohalide-respiring bacteria (OHRB) is important for effective bioremediation of chlorinated ethenes. The impact of iron and sulfate reduction on cis-1,2-dichloroethene (cDCE) and vinyl chloride (VC) dechlorination by a consortium containing the OHRB Dehalococcoides spp. was investigated using multiphase batch experiments. The OHRB consortium was found to contain endogenous iron- and sulfate-reducing bacteria (FeRB and SRB). A biogeochemical model was developed and used to quantify the mass transfer, aquatic geochemical, and microbial processes that occurred in the multiphase batch system. It was determined that the added SRB had the most significant impact on contaminant degradation. Addition of the SRB increased maximum specific substrate utilization rates, k_max, of cDCE and VC by 129% and 294%, respectively. The added FeRB had a slight stimulating effect on VC dechlorination when exogenous SRB were absent, but when cultured with the added SRB, FeRB moderated the SRB's stimulating effect. This study demonstrates that subsurface microbial community interactions are more complex than categorical, guild-based competition for resources such as electron donor.

Effects of Arsenic on Trichloroethene-Dechlorination Activities of Dehalococcoides mccartyi 195

Arsenic and trichloroethene (TCE) are among the most prevalent groundwater contaminants in the United States. Co-contamination of these two compounds has been detected at 63% of current TCE-contaminated National Priorities List sites. When in situ TCE reductive dechlorination is stimulated by the addition of fermentable substrates to generate a reducing environment, the presence of arsenic can be problematic because of the potential for increased mobilization and toxicity caused by the reduction of arsenate [As(V)] to arsenite [As(III)]. This study assesses the effects of arsenic exposure on the TCE-dechlorinating activities of Dehalococcoides mccartyi strain 195. Our results indicate that 9.1 μM As(III) caused a 50% decrease in D. mccartyi cell growth. While As(V) concentrations up to 200 μM did not initially impact TCE dechlorination, inhibition was observed in cultures amended with 200 μM As(V) and 100 μM As(V) in 12 and 17 days, respectively, corresponding with the accumulation of As(III). Transcriptomic and metabolomic analyses were performed to evaluate cellular responses to both As(V) and As(III) stress. Amendment of amino acids enhanced arsenic tolerance of D. mccartyi. Results from this study improve our understanding of potential inhibitions of D. mccartyi metabolism caused by arsenic and can inform the design of bioremediation strategies at co-contaminated sites.

Heteroexpression and biochemical characterization of thermostable citrate synthase from the cyanobacteria Anabaena sp. PCC7120

The present study recombinantly expressed a citrate synthase from cyanobacteria Anabaena sp. PCC7120 (AnCS) in Escherichia coli and characterized its enzymatic activity. The molecular mass of native AnCS was 88,533.1 Da containing two 44,162.7 Da subunits. Recombinant AnCS revealed the highest activity at pH 9.0 and 25 °C. AnCS displayed high thermal stability with a half-life time (t_1/2) of approximately 6.5 h at 60 °C, which was more thermostable than most CS from general organisms, but less than those from hyperthermophilic bacteria. The K_m values of oxaloacetate and acetyl-CoA were 138.50 and 18.15 μM respectively, suggesting a higher affinity to acetyl-CoA than oxaloacetate. Our inhibition assays showed that AnCS activity was not severely affected by most metal ions, but was strongly inhibited by Cu²⁺ and Zn²⁺. Treatments with ATP, ADP, AMP, NADH, and DTT depressed the AnCS activity. Overall, our results provide information on the enzymatic properties of AnCS, which contributes to the basic knowledge on CS selection for industrial utilizations.

Microbial synergistic interactions for reductive dechlorination of polychlorinated biphenyls

Dehalococcoides usually work closely with other beneficial microorganisms for removal of halogenated organic compounds at contaminated sites. Traditional microbial cultivation is necessary but not enough to gain insights into key microbial populations and their interactions in complex communities. In this study, we cultivated and characterized two D. mccartyi strains (CG3 and SG1), and further revealed interspecies synergistic interactions in PCB-dechlorinating microbial communities via metagenomic analysis. Strain CG3 and SG1 originated from distinct geographic sites employ reductive dehalogenase CG3-RD11 (PcbA1-like) and SG1-RD28 (PcbA4/5-like), respectively, to catalyze chlorine-removal from PCBs. In their parent mixed cultures CG-3 and SG-1, as well as in previously enriched PCB-dechlorinating cultures CG-1, CG-4 and CG-5, Methanosarcina and Desulfovibrio were found as major non-dechlorinating populations which may play roles in mediating acetate- and H₂-sources for D. mccartyi. They together form a stable microbial community for interspecies carbon- and electron-transfers to facilitate organohalide respiration of D. mccartyi, being confirmed in a synthetic microbial community consisting of the Dehalococcoides, Methanosarcina and Desulfovibrio. The results provide insights into which and how other microorganisms support D. mccartyi to dechlorinate PCBs, and suggest that Methanosarcina may play a larger role in PCB-dechlorinating communities than currently appreciated.

An integrative overview of genomic, transcriptomic and proteomic analyses in organohalide respiration research

Organohalide respiration (OHR) is a crucial process in the global halogen cycle and of interest for bioremediation. However, investigations on OHR are hampered by the restricted genetic accessibility and the poor growth yields of many organohalide-respiring bacteria (OHRB). Therefore, genomics, transcriptomics and proteomics are often used to investigate OHRB. In general, these gene expression studies are more useful when the data of the different 'omics' approaches are integrated and compared among a wide range of cultivation conditions and ideally involve several closely related OHRB. Despite the availability of a couple of proteomic and transcriptomic datasets dealing with OHRB, such approaches are currently not covered in reviews. Therefore, we here present an integrative and comparative overview of omics studies performed with the OHRB Sulfurospirillum multivorans, Dehalococcoides mccartyi, Desulfitobacterium spp. and Dehalobacter restrictus. Genes, transcripts, proteins and the regulatory and biochemical processes involved in OHR are discussed, and a comprehensive view on the unusual metabolism of D. mccartyi, which is one of the few bacteria possibly using a quinone-independent respiratory chain, is provided. Several 'omics'-derived theories on OHRB, e.g. the organohalide-respiratory chain, hydrogen metabolism, corrinoid biosynthesis or one-carbon metabolism are critically discussed on the basis of this integrative approach.

Metabolic models and gene essentiality data reveal essential and conserved metabolism in prokaryotes

Essential metabolic reactions are shaping constituents of metabolic networks, enabling viable and distinct phenotypes across diverse life forms. Here we analyse and compare modelling predictions of essential metabolic functions with experimental data and thereby identify core metabolic pathways in prokaryotes. Simulations of 15 manually curated genome-scale metabolic models were integrated with 36 large-scale gene essentiality datasets encompassing a wide variety of species of bacteria and archaea. Conservation of metabolic genes was estimated by analysing 79 representative genomes from all the branches of the prokaryotic tree of life. We find that essentiality patterns reflect phylogenetic relations both for modelling and experimental data, which correlate highly at the pathway level. Genes that are essential for several species tend to be highly conserved as opposed to non-essential genes which may be conserved or not. The tRNA-charging module is highlighted as ancestral and with high centrality in the networks, followed closely by cofactor metabolism, pointing to an early information processing system supplied by organic cofactors. The results, which point to model improvements and also indicate faults in the experimental data, should be relevant to the study of centrality in metabolic networks and ancient metabolism but also to metabolic engineering with prokaryotes.

If we tried to list every known chemical reaction within an organism–human, plant or even bacteria–we would get quite a long and confusing read. But when this information is represented in so-called genome-scale metabolic networks, we have the means to access computationally each of those reactions and their interconnections. Some parts of the network have alternatives, while others are unique and therefore can be essential for growth. Here, we simulate growth and compare essential reactions and genes for the simplest type of unicellular species–prokaryotes–to understand which parts of their metabolism are universally essential and potentially ancestral. We show that similar patterns of essential reactions echo phylogenetic relationships (this makes sense, as the genome provides the building plan for the enzymes that perform those reactions). Our computational predictions correlate strongly with experimental essentiality data. Finally, we show that a crucial step of protein synthesis (tRNA charging) and the synthesis and transformation of small molecules that enzymes require (cofactors) are the most essential and conserved parts of metabolism in prokaryotes. Our results are a step further in understanding the biology and evolution of prokaryotes but can also be relevant in applied studies including metabolic engineering and antibiotic design.

Analysis of Piscirickettsia salmonis Metabolism Using Genome-Scale Reconstruction, Modeling, and Testing

Piscirickettsia salmonis is an intracellular bacterial fish pathogen that causes piscirickettsiosis, a disease with highly adverse impact in the Chilean salmon farming industry. The development of effective treatment and control methods for piscireckttsiosis is still a challenge. To meet it the number of studies on P. salmonis has grown in the last couple of years but many aspects of the pathogen's biology are still poorly understood. Studies on its metabolism are scarce and only recently a metabolic model for reference strain LF-89 was developed. We present a new genome-scale model for P. salmonis LF-89 with more than twice as many genes as in the previous model and incorporating specific elements of the fish pathogen metabolism. Comparative analysis with models of different bacterial pathogens revealed a lower flexibility in P. salmonis metabolic network. Through constraint-based analysis, we determined essential metabolites required for its growth and showed that it can benefit from different carbon sources tested experimentally in new defined media. We also built an additional model for strain A1-15972, and together with an analysis of P. salmonis pangenome, we identified metabolic features that differentiate two main species clades. Both models constitute a knowledge-base for P. salmonis metabolism and can be used to guide the efficient culture of the pathogen and the identification of specific drug targets.

Systems Biology Approach to Bioremediation of Nitroaromatics: Constraint-Based Analysis of 2,4,6-Trinitrotoluene Biotransformation by Escherichia coli

Microbial remediation of nitroaromatic compounds (NACs) is a promising environmentally friendly and cost-effective approach to the removal of these life-threating agents. Escherichia coli (E. coli) has shown remarkable capability for the biotransformation of 2,4,6-trinitro-toluene (TNT). Efforts to develop E. coli as an efficient TNT degrading biocatalyst will benefit from holistic flux-level description of interactions between multiple TNT transforming pathways operating in the strain. To gain such an insight, we extended the genome-scale constraint-based model of E. coli to account for a curated version of major TNT transformation pathways known or evidently hypothesized to be active in E. coli in present of TNT. Using constraint-based analysis (CBA) methods, we then performed several series of in silico experiments to elucidate the contribution of these pathways individually or in combination to the E. coli TNT transformation capacity. Results of our analyses were validated by replicating several experimentally observed TNT degradation phenotypes in E. coli cultures. We further used the extended model to explore the influence of process parameters, including aeration regime, TNT concentration, cell density, and carbon source on TNT degradation efficiency. We also conducted an in silico metabolic engineering study to design a series of E. coli mutants capable of degrading TNT at higher yield compared with the wild-type strain. Our study, therefore, extends the application of CBA to bioremediation of nitroaromatics and demonstrates the usefulness of this approach to inform bioremediation research.

Refined experimental annotation reveals conserved corrinoid autotrophy in chloroform-respiring Dehalobacter isolates

Two novel chlorinated alkane-respiring Dehalobacter restrictus strains CF and DCA were isolated from the same enrichment culture, ACT-3, and characterized. The closed genomes of these highly similar sister strains were previously assembled from metagenomic sequence data and annotated. The isolation of the strains enabled experimental verification of predicted annotations, particularly focusing on irregularities or predicted gaps in central metabolic pathways and cofactor biosynthesis. Similar to D. restrictus strain PER-K23, strains CF and DCA require arginine, histidine and threonine for growth, although the corresponding biosynthesis pathways are predicted to be functional. Using strain CF to experimentally verify annotations, we determined that the predicted defective serine biosynthesis pathway can be rescued with a promiscuous serine hydroxymethyltransferase. Strain CF grew without added thiamine although the thiamine biosynthesis pathway is predicted to be absent; intracellular thiamine diphosphate, the cofactor of carboxylases in central metabolism, was not detected in cell extracts. Thus, strain CF may use amino acids to replenish central metabolites, portending entangled metabolite exchanges in ACT-3. Consistent with annotation, strain CF possesses a functional corrinoid biosynthesis pathway, demonstrated by increasing corrinoid content during growth and guided cobalamin biosynthesis in corrinoid-free medium. Chloroform toxicity to corrinoid-producing methanogens and acetogens may drive the conservation of corrinoid autotrophy in Dehalobacter strains. Heme detection in strain CF cell extracts suggests the 'archaeal' heme biosynthesis pathway also functions in anaerobic Firmicutes. This study reinforces the importance of incorporating enzyme promiscuity and cofactor availability in genome-scale functional predictions and identifies essential nutrient interdependencies in anaerobic dechlorinating microbial communities.

Sister Dehalobacter Genomes Reveal Specialization in Organohalide Respiration and Recent Strain Differentiation Likely Driven by Chlorinated Substrates

The genomes of two closely related Dehalobacter strains (strain CF and strain DCA) were assembled from the metagenome of an anaerobic enrichment culture that reductively dechlorinates chloroform (CF), 1,1,1-trichloroethane (1,1,1-TCA) and 1,1-dichloroethane (1,1-DCA). The 3.1 Mbp genomes of strain CF (that dechlorinates CF and 1,1,1-TCA) and strain DCA (that dechlorinates 1,1-DCA) each contain 17 putative reductive dehalogenase homologous (rdh) genes. These two genomes were systematically compared to three other available organohalide-respiring Dehalobacter genomes (Dehalobacter restrictus strain PER-K23, Dehalobacter sp. strain E1 and Dehalobacter sp. strain UNSWDHB), and to the genomes of Dehalococcoides mccartyi strain 195 and Desulfitobacterium hafniense strain Y51. This analysis compared 42 different metabolic and physiological categories. The genomes of strains CF and DCA share 90% overall average nucleotide identity and >99.8% identity over a 2.9 Mbp alignment that excludes large insertions, indicating that these genomes differentiated from a close common ancestor. This differentiation was likely driven by selection pressures around two orthologous reductive dehalogenase genes, cfrA and dcrA, that code for the enzymes that reduce CF or 1,1,1-TCA and 1,1-DCA. The many reductive dehalogenase genes found in the five Dehalobacter genomes cluster into two small conserved regions and were often associated with Crp/Fnr transcriptional regulators. Specialization is on-going on a strain-specific basis, as some strains but not others have lost essential genes in the Wood-Ljungdahl (strain E1) and corrinoid biosynthesis pathways (strains E1 and PER-K23). The gene encoding phosphoserine phosphatase, which catalyzes the last step of serine biosynthesis, is missing from all five Dehalobacter genomes, yet D. restrictus can grow without serine, suggesting an alternative or unrecognized biosynthesis route exists. In contrast to D. mccartyi, a complete heme biosynthesis pathway is present in the five Dehalobacter genomes. This pathway corresponds to a newly described alternative heme biosynthesis route first identified in Archaea. This analysis of organohalide-respiring Firmicutes and Chloroflexi reveals profound evolutionary differences despite very similar niche-specific metabolism and function.

PSAMM: A Portable System for the Analysis of Metabolic Models

The genome-scale models of metabolic networks have been broadly applied in phenotype prediction, evolutionary reconstruction, community functional analysis, and metabolic engineering. Despite the development of tools that support individual steps along the modeling procedure, it is still difficult to associate mathematical simulation results with the annotation and biological interpretation of metabolic models. In order to solve this problem, here we developed a Portable System for the Analysis of Metabolic Models (PSAMM), a new open-source software package that supports the integration of heterogeneous metadata in model annotations and provides a user-friendly interface for the analysis of metabolic models. PSAMM is independent of paid software environments like MATLAB, and all its dependencies are freely available for academic users. Compared to existing tools, PSAMM significantly reduced the running time of constraint-based analysis and enabled flexible settings of simulation parameters using simple one-line commands. The integration of heterogeneous, model-specific annotation information in PSAMM is achieved with a novel format of YAML-based model representation, which has several advantages, such as providing a modular organization of model components and simulation settings, enabling model version tracking, and permitting the integration of multiple simulation problems. PSAMM also includes a number of quality checking procedures to examine stoichiometric balance and to identify blocked reactions. Applying PSAMM to 57 models collected from current literature, we demonstrated how the software can be used for managing and simulating metabolic models. We identified a number of common inconsistencies in existing models and constructed an updated model repository to document the resolution of these inconsistencies.

The broad application of genome-scale metabolic modeling has made it a useful technique for tackling fundamental questions in biological research and engineering. Today over 100 models have been constructed for organisms that carry out a diverse array of metabolic activities spanning all three kingdoms of life. These models, however, have been curated independently following different conventions. The maintenance of model consistency has been challenging due to the lack of consensus in model representation and the absence of integrated modeling software for associating mathematical simulations with the annotation and biological interpretation of metabolic models. To solve this problem, we developed a new software package, PSAMM, and a new model format that incorporates heterogeneous, model-specific annotation information into modular representations of model definitions and simulation settings. PSAMM provides significant advances in standardizing the workflow of model annotation and consistency checking. Compared to existing tools, PSAMM supports more flexible configurations and is more efficient in running constraint-based simulations. All functions of PSAMM are freely available for academic users and can be downloaded from a public Git repository (https://zhanglab.github.io/psamm/) under the GNU General Public License.

Experimental validation of in silico model-predicted isocitrate dehydrogenase and phosphomannose isomerase from Dehalococcoides mccartyi

Gene sequences annotated as proteins of unknown or non‐specific function and hypothetical proteins account for a large fraction of most genomes. In the strictly anaerobic and organohalide respiring Dehalococcoides mccartyi, this lack of annotation plagues almost half the genome. Using a combination of bioinformatics analyses and genome‐wide metabolic modelling, new or more specific annotations were proposed for about 80 of these poorly annotated genes in previous investigations of D. mccartyi metabolism. Herein, we report the experimental validation of the proposed reannotations for two such genes (KB1_0495 and KB1_0553) from D. mccartyi strains in the KB‐1 community. KB1_0495 or DmIDH was originally annotated as an NAD⁺‐dependent isocitrate dehydrogenase, but biochemical assays revealed its activity primarily with NADP⁺ as a cofactor. KB1_0553, also denoted as DmPMI, was originally annotated as a hypothetical protein/sugar isomerase domain protein. We previously proposed that it was a bifunctional phosphoglucose isomerase/phosphomannose isomerase, but only phosphomannose isomerase activity was identified and confirmed experimentally. Further bioinformatics analyses of these two protein sequences suggest their affiliation to potentially novel enzyme families within their respective larger enzyme super families.

The bacterial pangenome as a new tool for analysing pathogenic bacteria

The bacterial pangenome was introduced in 2005 and, in recent years, has been the subject of many studies. Thanks to progress in next-generation sequencing methods, the pangenome can be divided into two parts, the core (common to the studied strains) and the accessory genome, offering a large panel of uses. In this review, we have presented the analysis methods, the pangenome composition and its application as a study of lifestyle. We have also shown that the pangenome may be used as a new tool for redefining the pathogenic species. We applied this to the Escherichia coli and Shigella species, which have been a subject of controversy regarding their taxonomic and pathogenic position.

•

Pangenome is a new way of studying pathogenic bacteria.

•

Pangenome can be used as a taxonomic tool.

•

This review describes pangenome in the world of pathogenic bacteria.

New insights into Dehalococcoides mccartyi metabolism from a reconstructed metabolic network-based systems-level analysis of D. mccartyi transcriptomes

Organohalide respiration, mediated by Dehalococcoides mccartyi, is a useful bioremediation process that transforms ground water pollutants and known human carcinogens such as trichloroethene and vinyl chloride into benign ethenes. Successful application of this process depends on the fundamental understanding of the respiration and metabolism of D. mccartyi. Reductive dehalogenases, encoded by rdhA genes of these anaerobic bacteria, exclusively catalyze organohalide respiration and drive metabolism. To better elucidate D. mccartyi metabolism and physiology, we analyzed available transcriptomic data for a pure isolate (Dehalococcoides mccartyi strain 195) and a mixed microbial consortium (KB-1) using the previously developed pan-genome-scale reconstructed metabolic network of D. mccartyi. The transcriptomic data, together with available proteomic data helped confirm transcription and expression of the majority genes in D. mccartyi genomes. A composite genome of two highly similar D. mccartyi strains (KB-1 Dhc) from the KB-1 metagenome sequence was constructed, and operon prediction was conducted for this composite genome and other single genomes. This operon analysis, together with the quality threshold clustering analysis of transcriptomic data helped generate experimentally testable hypotheses regarding the function of a number of hypothetical proteins and the poorly understood mechanism of energy conservation in D. mccartyi. We also identified functionally enriched important clusters (13 for strain 195 and 11 for KB-1 Dhc) of co-expressed metabolic genes using information from the reconstructed metabolic network. This analysis highlighted some metabolic genes and processes, including lipid metabolism, energy metabolism, and transport that potentially play important roles in organohalide respiration. Overall, this study shows the importance of an organism's metabolic reconstruction in analyzing various "omics" data to obtain improved understanding of the metabolism and physiology of the organism.

Incomplete Wood-Ljungdahl pathway facilitates one-carbon metabolism in organohalide-respiring Dehalococcoides mccartyi

The acetyl-CoA "Wood-Ljungdahl" pathway couples the folate-mediated one-carbon (C1) metabolism to either CO2 reduction or acetate oxidation via acetyl-CoA. This pathway is distributed in diverse anaerobes and is used for both energy conservation and assimilation of C1 compounds. Genome annotations for all sequenced strains of Dehalococcoides mccartyi, an important bacterium involved in the bioremediation of chlorinated solvents, reveal homologous genes encoding an incomplete Wood-Ljungdahl pathway. Because this pathway lacks key enzymes for both C1 metabolism and CO2 reduction, its cellular functions remain elusive. Here we used D. mccartyi strain 195 as a model organism to investigate the metabolic function of this pathway and its impacts on the growth of strain 195. Surprisingly, this pathway cleaves acetyl-CoA to donate a methyl group for production of methyl-tetrahydrofolate (CH3-THF) for methionine biosynthesis, representing an unconventional strategy for generating CH3-THF in organisms without methylene-tetrahydrofolate reductase. Carbon monoxide (CO) was found to accumulate as an obligate by-product from the acetyl-CoA cleavage because of the lack of a CO dehydrogenase in strain 195. CO accumulation inhibits the sustainable growth and dechlorination of strain 195 maintained in pure cultures, but can be prevented by CO-metabolizing anaerobes that coexist with D. mccartyi, resulting in an unusual syntrophic association. We also found that this pathway incorporates exogenous formate to support serine biosynthesis. This study of the incomplete Wood-Ljungdahl pathway in D. mccartyi indicates a unique bacterial C1 metabolism that is critical for D. mccartyi growth and interactions in dechlorinating communities and may play a role in other anaerobic communities.

Kinetics of dechlorination by Dehalococcoides mccartyi using different carbon sources

Stimulated anaerobic dechlorination is generally considered a valuable step for the remediation of aquifers polluted with chlorinated ethenes (CEs). Correct simulation and prediction of this process in situ, however, require good knowledge of the associated biological reactions. The aim of this study was to evaluate the dechlorination reaction in an aquifer contaminated with trichloroethene (TCE) and its daughter products, discharging into the Zenne River. Different carbon sources were used in batch cultures and these were related to the dechlorination reaction, together with the monitored biomarkers. Appropriate kinetic formulations were assessed. Reductive dechlorination of TCE took place only when external carbon sources were added to microcosms, and occurred concomitant with a pronounced increase in the Dehalococcoides mccartyi cell count as determined by 16S rRNA gene-targeted qPCR. This indicates that native dechlorinating bacteria are present in the aquifer of the Zenne site and that the oligotrophic nature of the aquifer prevents a complete degradation to ethene. The type of carbon source, the cell number of D. mccartyi or the reductive dehalogenase genes, however, did not unequivocally explain the observed differences in degradation rates or the extent of dechlorination. Neither first-order, Michaelis-Menten nor Monod kinetics could perfectly simulate the dechlorination reactions in TCE spiked microcosms. A sensitivity analysis indicated that the inclusion of donor limitation would not significantly enhance the simulations without a clear process understanding. Results point to the role of the supporting microbial community but it remains to be verified how the complexity of the microbial (inter)actions should be represented in a model framework.

Simple topological properties predict functional misannotations in a metabolic network

Motivation: Misannotation in sequence databases is an important obstacle for automated tools for gene function annotation, which rely extensively on comparison with sequences with known function. To improve current annotations and prevent future propagation of errors, sequence-independent tools are, therefore, needed to assist in the identification of misannotated gene products. In the case of enzymatic functions, each functional assignment implies the existence of a reaction within the organism’s metabolic network; a first approximation to a genome-scale metabolic model can be obtained directly from an automated genome annotation. Any obvious problems in the network, such as dead end or disconnected reactions, can, therefore, be strong indications of misannotation.

Results: We demonstrate that a machine-learning approach using only network topological features can successfully predict the validity of enzyme annotations. The predictions are tested at three different levels. A random forest using topological features of the metabolic network and trained on curated sets of correct and incorrect enzyme assignments was found to have an accuracy of up to 86% in 5-fold cross-validation experiments. Further cross-validation against unseen enzyme superfamilies indicates that this classifier can successfully extrapolate beyond the classes of enzyme present in the training data. The random forest model was applied to several automated genome annotations, achieving an accuracy of in most cases when validated against recent genome-scale metabolic models. We also observe that when applied to draft metabolic networks for multiple species, a clear negative correlation is observed between predicted annotation quality and phylogenetic distance to the major model organism for biochemistry (Escherichia coli for prokaryotes and Homo sapiens for eukaryotes).

Contact:j.pinney@imperial.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

Community genomic analyses constrain the distribution of metabolic traits across the Chloroflexi phylum and indicate roles in sediment carbon cycling

BACKGROUND

Sediments are massive reservoirs of carbon compounds and host a large fraction of microbial life. Microorganisms within terrestrial aquifer sediments control buried organic carbon turnover, degrade organic contaminants, and impact drinking water quality. Recent 16S rRNA gene profiling indicates that members of the bacterial phylum Chloroflexi are common in sediment. Only the role of the class Dehalococcoidia, which degrade halogenated solvents, is well understood. Genomic sampling is available for only six of the approximate 30 Chloroflexi classes, so little is known about the phylogenetic distribution of reductive dehalogenation or about the broader metabolic characteristics of Chloroflexi in sediment.

RESULTS

We used metagenomics to directly evaluate the metabolic potential and diversity of Chloroflexi in aquifer sediments. We sampled genomic sequence from 86 Chloroflexi representing 15 distinct lineages, including members of eight classes previously characterized only by 16S rRNA sequences. Unlike in the Dehalococcoidia, genes for organohalide respiration are rare within the Chloroflexi genomes sampled here. Near-complete genomes were reconstructed for three Chloroflexi. One, a member of an unsequenced lineage in the Anaerolinea, is an aerobe with the potential for respiring diverse carbon compounds. The others represent two genomically unsampled classes sibling to the Dehalococcoidia, and are anaerobes likely involved in sugar and plant-derived-compound degradation to acetate. Both fix CO2 via the Wood-Ljungdahl pathway, a pathway not previously documented in Chloroflexi. The genomes each encode unique traits apparently acquired from Archaea, including mechanisms of motility and ATP synthesis.

CONCLUSIONS

Chloroflexi in the aquifer sediments are abundant and highly diverse. Genomic analyses provide new evolutionary boundaries for obligate organohalide respiration. We expand the potential roles of Chloroflexi in sediment carbon cycling beyond organohalide respiration to include respiration of sugars, fermentation, CO2 fixation, and acetogenesis with ATP formation by substrate-level phosphorylation.

Complete genome sequence of Dehalobacter restrictus PER-K23(T.)

Dehalobacter restrictus strain PER-K23 (DSM 9455) is the type strain of the species Dehalobacter restrictus. D. restrictus strain PER-K23 grows by organohalide respiration, coupling the oxidation of H₂ to the reductive dechlorination of tetra- or trichloroethene. Growth has not been observed with any other electron donor or acceptor, nor has fermentative growth been shown. Here we introduce the first full genome of a pure culture within the genus Dehalobacter. The 2,943,336 bp long genome contains 2,826 protein coding and 82 RNA genes, including 5 16S rRNA genes. Interestingly, the genome contains 25 predicted reductive dehalogenase genes, the majority of which appear to be full length. The reductive dehalogenase genes are mainly located in two clusters, suggesting a much larger potential for organohalide respiration than previously anticipated.

The restricted metabolism of the obligate organohalide respiring bacterium Dehalobacter restrictus: lessons from tiered functional genomics

Dehalobacter restrictus strain PER-K23 is an obligate organohalide respiring bacterium, which displays extremely narrow metabolic capabilities. It grows only via coupling energy conservation to anaerobic respiration of tetra- and trichloroethene with hydrogen as sole electron donor. Dehalobacter restrictus represents the paradigmatic member of the genus Dehalobacter, which in recent years has turned out to be a major player in the bioremediation of an increasing number of organohalides, both in situ and in laboratory studies. The recent elucidation of the D. restrictus genome revealed a rather elaborate genome with predicted pathways that were not suspected from its restricted metabolism, such as a complete corrinoid biosynthetic pathway, the Wood-Ljungdahl (WL) pathway for CO2 fixation, abundant transcriptional regulators and several types of hydrogenases. However, one important feature of the genome is the presence of 25 reductive dehalogenase genes, from which so far only one, pceA, has been characterized on genetic and biochemical levels. This study describes a multi-level functional genomics approach on D. restrictus across three different growth phases. A global proteomic analysis allowed consideration of general metabolic pathways relevant to organohalide respiration, whereas the dedicated genomic and transcriptomic analysis focused on the diversity, composition and expression of genes associated with reductive dehalogenases.

Network motifs provide signatures that characterize metabolism

Motifs are repeating patterns that determine the local properties of networks. In this work, we characterized all 3-node motifs using enzyme commission numbers of the International Union of Biochemistry and Molecular Biology to show that motif abundance is related to biochemical function. Further, we present a comparative analysis of motif distributions in the metabolic networks of 21 species across six kingdoms of life. We found the distribution of motif abundances to be similar between species, but unique across cellular organelles. Finally, we show that motifs are able to capture inter-species differences in metabolic networks and that molecular differences between some biological species are reflected by the distribution of motif abundances in metabolic networks.

Metagenome analysis reveals yet unexplored reductive dechlorinating potential of Dehalobacter sp. E1 growing in co-culture with Sedimentibacter sp

The importance of Dehalobacter species in bioremediation as dedicated degraders of chlorinated organics has been well recognized. However, still little is known about Dehalobacter's full genomic repertoires, including the genes involved in dehalogenation. Here we report the first insights into the genome sequence of Dehalobacter sp. E1 that grows in strict co-culture with Sedimentibacter sp. B4. Based on the co-culture metagenome and the genome of strain B4 (4.2 Mbp) we estimate the genome sequence of strain E1 to be 2.6 Mbp. Ten putative reductive dehalogenase homologue (Rdh)-encoding gene clusters were identified. One cluster has a putative tetrachloroethene Rdh-encoding gene cluster, similar to the pceABCT operon previously identified in Dehalobacter restrictus. Metagenome analysis indicated that the inability of strain E1 to synthesize cobalamin, an essential cofactor of reductive dehalogenases, is complemented by Sedimentibacter. The metagenomic exploration described here maps the extensive dechlorinating potential of Dehalobacter, and paves way for elucidation of the interactions with its co-cultured Sedimentibacter.

Ecogenomics of microbial communities in bioremediation of chlorinated contaminated sites

Organohalide compounds such as chloroethenes, chloroethanes, and polychlorinated benzenes are among the most significant pollutants in the world. These compounds are often found in contamination plumes with other pollutants such as solvents, pesticides, and petroleum derivatives. Microbial bioremediation of contaminated sites, has become commonplace whereby key processes involved in bioremediation include anaerobic degradation and transformation of these organohalides by organohalide respiring bacteria and also via hydrolytic, oxygenic, and reductive mechanisms by aerobic bacteria. Microbial ecogenomics has enabled us to not only study the microbiology involved in these complex processes but also develop tools to better monitor and assess these sites during bioremediation. Microbial ecogenomics have capitalized on recent advances in high-throughput and -output genomics technologies in combination with microbial physiology studies to address these complex bioremediation problems at a system level. Advances in environmental metagenomics, transcriptomics, and proteomics have provided insights into key genes and their regulation in the environment. They have also given us clues into microbial community structures, dynamics, and functions at contaminated sites. These techniques have not only aided us in understanding the lifestyles of common organohalide respirers, for example Dehalococcoides, Dehalobacter, and Desulfitobacterium, but also provided insights into novel and yet uncultured microorganisms found in organohalide respiring consortia. In this paper, we look at how ecogenomic studies have aided us to understand the microbial structures and functions in response to environmental stimuli such as the presence of chlorinated pollutants.

Comparative metagenomics of three Dehalococcoides-containing enrichment cultures: the role of the non-dechlorinating community

Background

The Dehalococcoides are strictly anaerobic bacteria that gain metabolic energy via the oxidation of H₂ coupled to the reduction of halogenated organic compounds. Dehalococcoides spp. grow best in mixed microbial consortia, relying on non-dechlorinating members to provide essential nutrients and maintain anaerobic conditions.

A metagenome sequence was generated for the dechlorinating mixed microbial consortium KB-1. A comparative metagenomic study utilizing two additional metagenome sequences for Dehalococcoides-containing dechlorinating microbial consortia was undertaken to identify common features that are provided by the non-dechlorinating community and are potentially essential to Dehalococcoides growth.

Results

The KB-1 metagenome contained eighteen novel homologs to reductive dehalogenase genes. The metagenomes obtained from the three consortia were automatically annotated using the MG-RAST server, from which statistically significant differences in community composition and metabolic profiles were determined. Examination of specific metabolic pathways, including corrinoid synthesis, methionine synthesis, oxygen scavenging, and electron-donor metabolism identified the Firmicutes, methanogenic Archaea, and the ∂-Proteobacteria as key organisms encoding these pathways, and thus potentially producing metabolites required for Dehalococcoides growth.

Conclusions

Comparative metagenomics of the three Dehalococcoides-containing consortia identified that similarities across the three consortia are more apparent at the functional level than at the taxonomic level, indicating the non-dechlorinating organisms’ identities can vary provided they fill the same niche within a consortium. Functional redundancy was identified in each metabolic pathway of interest, with key processes encoded by multiple taxonomic groups. This redundancy likely contributes to the robust growth and dechlorination rates in dechlorinating enrichment cultures.

Stable isotope peptide mass spectrometry to decipher amino acid metabolism in Dehalococcoides strain CBDB1

Dehalococcoides species are key players in the anaerobic transformation of halogenated solvents at contaminated sites. Here, we analyze isotopologue distributions in amino acid pools from peptides of Dehalococcoides strain CBDB1 after incubation with (13)C-labeled acetate or bicarbonate as a carbon source. The resulting data were interpreted with regard to genome annotations to identify amino acid biosynthesis pathways. In addition to using gas chromatography-mass spectrometry (GC-MS) for analyzing derivatized amino acids after protein hydrolysis, we introduce a second, much milder method, in which we directly analyze peptide masses after tryptic digest and peptide fragments by nano-liquid chromatography-electrospray ionization-tandem mass spectrometry (nano-LC-ESI-MS/MS). With this method, we identify isotope incorporation patterns for 17 proteinaceous amino acids, including proline, cysteine, lysine, and arginine, which escaped previous analyses in Dehalococcoides. Our results confirmed lysine biosynthesis via the α-aminoadipate pathway, precluding lysine formation from aspartate. Similarly, the isotopologue pattern obtained for arginine provided biochemical evidence of its synthesis from glutamate. Direct peptide MS/MS analysis of the labeling patterns of glutamine and asparagine, which were converted to glutamate and aspartate during protein hydrolysis, gave biochemical evidence of their precursors and confirmed glutamate biosynthesis via a Re-specific citrate synthase. By addition of unlabeled free amino acids to labeled cells, we show that in strain CBDB1 none of the 17 tested amino acids was incorporated into cell mass, indicating that they are all synthesized de novo. Our approach is widely applicable and provides a means to analyze amino acid metabolism by studying specific proteins even in mixed consortia.

The design of long-term effective uranium bioremediation strategy using a community metabolic model

Acetate amendment at uranium contaminated sites in Rifle, CO. leads to an initial bloom of Geobacter accompanied by the removal of U(VI) from the groundwater, followed by an increase of sulfate-reducing bacteria (SRBs) which are poor reducers of U(VI). One of the challenges associated with bioremediation is the decay in Geobacter abundance, which has been attributed to the depletion of bio-accessible Fe(III), motivating the investigation of simultaneous amendments of acetate and Fe(III) as an alternative bioremediation strategy. In order to understand the community metabolism of Geobacter and SRBs during artificial substrate amendment, we have created a genome-scale dynamic community model of Geobacter and SRBs using the previously described Dynamic Multi-species Metabolic Modeling framework. Optimization techniques are used to determine the optimal acetate and Fe(III) addition profile. Field-scale simulation of acetate addition accurately predicted the in situ data. The simulations suggest that batch amendment of Fe(III) along with continuous acetate addition is insufficient to promote long-term bioremediation, while continuous amendment of Fe(III) along with continuous acetate addition is sufficient to promote long-term bioremediation. By computationally minimizing the acetate and Fe(III) addition rates as well as the difference between the predicted and target uranium concentration, we showed that it is possible to maintain the uranium concentration below the environmental safety standard while minimizing the cost of chemical additions. These simulations show that simultaneous addition of acetate and Fe(III) has the potential to be an effective uranium bioremediation strategy. They also show that computational modeling of microbial community is an important tool to design effective strategies for practical applications in environmental biotechnology.

Identification and characterization of a re-citrate synthase in Dehalococcoides strain CBDB1

The genome annotations of all sequenced Dehalococcoides strains lack a citrate synthase, although physiological experiments have indicated that such an activity should be encoded. We here report that a Re face-specific citrate synthase is synthesized by Dehalococcoides strain CBDB1 and that this function is encoded by the gene cbdbA1708 (NCBI accession number CAI83711), previously annotated as encoding homocitrate synthase. Gene cbdbA1708 was heterologously expressed in Escherichia coli, and the recombinant enzyme was purified. The enzyme catalyzed the condensation of oxaloacetate and acetyl coenzyme A (acetyl-CoA) to citrate. The protein did not have homocitrate synthase activity and was inhibited by citrate, and Mn2+ was needed for full activity. The stereospecificity of the heterologously expressed citrate synthase was determined by electrospray ionization liquid chromatography-mass spectrometry (ESI LC/MS). Citrate was synthesized from [2-(13)C]acetyl-CoA and oxaloacetate by the Dehalococcoides recombinant citrate synthase and then converted to acetate and malate by commercial citrate lyase plus malate dehydrogenase. The formation of unlabeled acetate and 13C-labeled malate proved the Re face-specific activity of the enzyme. Shotgun proteome analyses of cell extracts of strain CBDB1 demonstrated that cbdbA1708 is expressed in strain CBDB1.

Selective utilization of exogenous amino acids by Dehalococcoides ethenogenes strain 195 and its effects on growth and dechlorination activity

Bacteria of the genus Dehalococcoides are important members of bioremediation communities because of their ability to detoxify chloroethenes to the benign end product ethene. Genome-enabled studies conducted with Dehalococcoides ethenogenes 195 have revealed that two ATP-binding cassette (ABC)-type amino acid transporters are expressed during its exponential growth stages. In light of previous findings that Casamino Acids enhanced its dechlorination activity, we hypothesized that strain 195 is capable of importing amino acids from its environment to facilitate dechlorination and growth. To test this hypothesis, we applied isotopomer-based dilution analysis with (13)C-labeled acetate to differentiate the amino acids that were taken up by strain 195 from those synthesized de novo and to determine the physiological changes caused by the significantly incorporated amino acids. Our results showed that glutamate/glutamine and aspartate/asparagine were almost exclusively synthesized by strain 195, even when provided in excess in the medium. In contrast, phenylalanine, isoleucine, leucine, and methionine were identified as the four most highly incorporated amino acids, at levels >30% of respective proteinogenic amino acids. When either phenylalanine or all four highly incorporated amino acids were added to the defined mineral medium, the growth rates, dechlorination activities, and yields of strain 195 were enhanced to levels similar to those observed with supplementation with 20 amino acids. However, genes for the putative ABC-type amino acids transporters and phenylalanine biosynthesis exhibited insignificant regulation in response to the imported amino acids. This study also demonstrates that using isotopomer-based metabolite analysis can be an efficient strategy for optimizing nutritional conditions for slow-growing microorganisms.

Design and verification of a pangenome microarray oligonucleotide probe set for Dehalococcoides spp

Dehalococcoides spp. are an industrially relevant group of Chloroflexi bacteria capable of reductively dechlorinating contaminants in groundwater environments. Existing Dehalococcoides genomes revealed a high level of sequence identity within this group, including 98 to 100% 16S rRNA sequence identity between strains with diverse substrate specificities. Common molecular techniques for identification of microbial populations are often not applicable for distinguishing Dehalococcoides strains. Here we describe an oligonucleotide microarray probe set designed based on clustered Dehalococcoides genes from five different sources (strain DET195, CBDB1, BAV1, and VS genomes and the KB-1 metagenome). This "pangenome" probe set provides coverage of core Dehalococcoides genes as well as strain-specific genes while optimizing the potential for hybridization to closely related, previously unknown Dehalococcoides strains. The pangenome probe set was compared to probe sets designed independently for each of the five Dehalococcoides strains. The pangenome probe set demonstrated better predictability and higher detection of Dehalococcoides genes than strain-specific probe sets on nontarget strains with <99% average nucleotide identity. An in silico analysis of the expected probe hybridization against the recently released Dehalococcoides strain GT genome and additional KB-1 metagenome sequence data indicated that the pangenome probe set performs more robustly than the combined strain-specific probe sets in the detection of genes not included in the original design. The pangenome probe set represents a highly specific, universal tool for the detection and characterization of Dehalococcoides from contaminated sites. It has the potential to become a common platform for Dehalococcoides-focused research, allowing meaningful comparisons between microarray experiments regardless of the strain examined.

In situ to in silico and back: elucidating the physiology and ecology of Geobacter spp. using genome-scale modelling

There is a wide diversity of unexplored metabolism encoded in the genomes of microorganisms that have an important environmental role. Genome-scale metabolic modelling enables the individual reactions that are encoded in annotated genomes to be organized into a coherent whole, which can then be used to predict metabolic fluxes that will optimize cell function under a range of conditions. In this Review, we summarize a series of studies in which genome-scale metabolic modelling of Geobacter spp. has resulted in an in-depth understanding of their central metabolism and ecology. A similar iterative modelling and experimental approach could accelerate elucidation of the physiology and ecology of other microorganisms inhabiting a diversity of environments, and could guide optimization of the practical applications of these species.