Growth- and substrate-dependent transcription of the formate dehydrogenase (fdhCAB) operon in Methanobacterium thermoformicicum Z-245

Diverse electron carriers drive syntrophic interactions in an enriched anaerobic acetate-oxidizing consortium

In many anoxic environments, syntrophic acetate oxidation (SAO) is a key pathway mediating the conversion of acetate into methane through obligate cross-feeding interactions between SAO bacteria (SAOB) and methanogenic archaea. The SAO pathway is particularly important in engineered environments such as anaerobic digestion (AD) systems operating at thermophilic temperatures and/or with high ammonia. Despite the widespread importance of SAOB to the stability of the AD process, little is known about their in situ physiologies due to typically low biomass yields and resistance to isolation. Here, we performed a long-term (300-day) continuous enrichment of a thermophilic (55 °C) SAO community from a municipal AD system using acetate as the sole carbon source. Over 80% of the enriched bioreactor metagenome belonged to a three-member consortium, including an acetate-oxidizing bacterium affiliated with DTU068 encoding for carbon dioxide, hydrogen, and formate production, along with two methanogenic archaea affiliated with Methanothermobacter_A. Stable isotope probing was coupled with metaproteogenomics to quantify carbon flux into each community member during acetate conversion and inform metabolic reconstruction and genome-scale modeling. This effort revealed that the two Methanothermobacter_A species differed in their preferred electron donors, with one possessing the ability to grow on formate and the other only consuming hydrogen. A thermodynamic analysis suggested that the presence of the formate-consuming methanogen broadened the environmental conditions where ATP production from SAO was favorable. Collectively, these results highlight how flexibility in electron partitioning during SAO likely governs community structure and fitness through thermodynamic-driven mutualism, shedding valuable insights into the metabolic underpinnings of this key functional group within methanogenic ecosystems.

Evolving understanding of rumen methanogen ecophysiology

Production of methane by methanogenic archaea, or methanogens, in the rumen of ruminants is a thermodynamic necessity for microbial conversion of feed to volatile fatty acids, which are essential nutrients for the animals. On the other hand, methane is a greenhouse gas and its production causes energy loss for the animal. Accordingly, there are ongoing efforts toward developing effective strategies for mitigating methane emissions from ruminant livestock that require a detailed understanding of the diversity and ecophysiology of rumen methanogens. Rumen methanogens evolved from free-living autotrophic ancestors through genome streamlining involving gene loss and acquisition. The process yielded an oligotrophic lifestyle, and metabolically efficient and ecologically adapted descendants. This specialization poses serious challenges to the efforts of obtaining axenic cultures of rumen methanogens, and consequently, the information on their physiological properties remains in most part inferred from those of their non-rumen representatives. This review presents the current knowledge of rumen methanogens and their metabolic contributions to enteric methane production. It also identifies the respective critical gaps that need to be filled for aiding the efforts to mitigate methane emission from livestock operations and at the same time increasing the productivity in this critical agriculture sector.

An integrated systems biology approach reveals differences in formate metabolism in the genus Methanothermobacter

Methanogenesis allows methanogenic archaea to generate cellular energy for their growth while producing methane. Thermophilic hydrogenotrophic species of the genus Methanothermobacter have been recognized as robust biocatalysts for a circular carbon economy and are already applied in power-to-gas technology with biomethanation, which is a platform to store renewable energy and utilize captured carbon dioxide. Here, we generated curated genome-scale metabolic reconstructions for three Methanothermobacter strains and investigated differences in the growth performance of these same strains in chemostat bioreactor experiments with hydrogen and carbon dioxide or formate as substrates. Using an integrated systems biology approach, we identified differences in formate anabolism between the strains and revealed that formate anabolism influences the diversion of carbon between biomass and methane. This finding, together with the omics datasets and the metabolic models we generated, can be implemented for biotechnological applications of Methanothermobacter in power-to-gas technology, and as a perspective, for value-added chemical production.

Formate and hydrogen in hydrothermal vents and their use by extremely thermophilic methanogens and heterotrophs

Extremely thermophilic methanogens in the Methanococci and heterotrophs in the Thermococci are common in deep-sea hydrothermal vents. All Methanococci use H₂ as an electron donor, and a few species can also use formate. Most Methanococci have a coenzyme F₄₂₀-reducing formate dehydrogenase. All Thermococci reduce S⁰ but have hydrogenases and produce H₂ in the absence of S⁰. Some Thermococci have formate hydrogenlyase (Fhl) that reversibly converts H₂ and CO₂ to formate or an NAD(P)⁺-reducing formate dehydrogenase (Nfd). Questions remain if Methanococci or Thermococci use or produce formate in nature, why only certain species can grow on or produce formate, and what the physiological role of formate is? Formate forms abiotically in hydrothermal fluids through chemical equilibrium with primarily H₂, CO₂, and CO and is strongly dependent upon H₂ concentration, pH, and temperature. Formate concentrations are highest in hydrothermal fluids where H₂ concentrations are also high, such as in ultramafic systems where serpentinization reactions occur. In nature, Methanococci are likely to use formate as an electron donor when H₂ is limiting. Thermococci with Fhl likely convert H₂ and CO₂ to formate when H₂ concentrations become inhibitory for growth. They are unlikely to grow on formate in nature unless formate is more abundant than H₂ in the environment. Nearly all Methanococci and Thermococci have a gene for at least one formate dehydrogenase catalytic subunit, which may be used to provide free formate for de novo purine biosynthesis. However, only species with a membrane-bound formate transporter can grow on or secrete formate. Interspecies H₂ transfer occurs between Thermococci and Methanococci. This and putative interspecies formate transfer may support Methanococci in low H₂ environments, which in turn may prevent growth inhibition of Thermococci by its own H₂. Future research directions include understanding when, where, and how formate is used and produced by these organisms in nature, and how transcription of Thermococci genes encoding formate-related enzymes are regulated.

A Shuttle-Vector System Allows Heterologous Gene Expression in the Thermophilic Methanogen Methanothermobacter thermautotrophicus ΔH

Thermophilic Methanothermobacter spp. are used as model microbes to study the physiology and biochemistry of the conversion of molecular hydrogen and carbon dioxide into methane (i.e., hydrogenotrophic methanogenesis). Yet, a genetic system for these model microbes was missing despite intensive work for four decades. Here, we report the successful implementation of genetic tools for Methanothermobacter thermautotrophicus ΔH. We developed shuttle vectors that replicated in Escherichia coli and M. thermautotrophicus ΔH. For M. thermautotrophicus ΔH, a thermostable neomycin resistance cassette served as the selectable marker for positive selection with neomycin, and the cryptic plasmid pME2001 from Methanothermobacter marburgensis served as the replicon. The shuttle-vector DNA was transferred from E. coli into M. thermautotrophicus ΔH via interdomain conjugation. After the successful validation of DNA transfer and positive selection in M. thermautotrophicus ΔH, we demonstrated heterologous gene expression of a thermostable β-galactosidase-encoding gene (bgaB) from Geobacillus stearothermophilus under the expression control of four distinct synthetic and native promoters. In quantitative in-vitro enzyme activity assay, we found significantly different β-galactosidase activity with these distinct promoters. With a formate dehydrogenase operon-encoding shuttle vector, we allowed growth of M. thermautotrophicus ΔH on formate as the sole growth substrate, while this was not possible for the empty-vector control.

On the diversity of F420 -dependent oxidoreductases: A sequence- and structure-based classification

The F₄₂₀ deazaflavin cofactor is an intriguing molecule as it structurally resembles the canonical flavin cofactor, although behaves as a nicotinamide cofactor due to its obligate hydride-transfer reactivity and similar low redox potential. Since its discovery, numerous enzymes relying on it have been described. The known deazaflavoproteins are taxonomically restricted to Archaea and Bacteria. The biochemistry of the deazaflavoenzymes is diverse and they exhibit great structural variability. In this study a thorough sequence and structural homology evolutionary analysis was performed in order to generate an overarching classification of the F₄₂₀ -dependent oxidoreductases. Five different deazaflavoenzyme Classes (I-V) are described according to their structural folds as follows: Class I encompassing the TIM-barrel F₄₂₀ -dependent enzymes; Class II including the Rossmann fold F₄₂₀ -dependent enzymes; Class III comprising the β-roll F₄₂₀ -dependent enzymes; Class IV which exclusively gathers the SH3 barrel F₄₂₀ -dependent enzymes and Class V including the three layer ββα sandwich F₄₂₀ -dependent enzymes. This classification provides a framework for the identification and biochemical characterization of novel deazaflavoenzymes.

Comparative genomic analysis reveals starvation survival systems in Methanothermobacter thermautotrophicus ΔH

Methanothermobacter thermautotrophicus ΔH (MTH) is a thermophilic hydrogenotrophic methanogenic archaeon capable of reducing CO₂ with H₂ to produce methane gas. It is the potential candidate in the biomethanation of CO₂ and CO in anaerobic reactors and biogas upgrading process. However, systematic studies addressing its genome conservation and function remain scant in this genome. In this study, we have evaluated its evolutionary resemblance and metabolic discrepancy, particularly in starvation survival systems by comparing the genomic contexts with Methanothermobacter marburgensis str. Marburg (MMG) and Methanobacterium formicicum DSM 1535 (MFO). The phylogenomic analysis of this study indicated that there was a strong phylogenomic signal among MTH, MMG, and MFO in the whole-genome tree. DNA replication machinery was conserved in the MTH genome and might have evolved at different evolution rates. Genome synteny analysis observed collinearity of either gene orders or gene families has to be maintained with syntenic blocks located in the syntenic out-paralogs. A genome-wide metabolic analysis identified some unique putative metabolic subsystems in MTH, which are proposed to determine its growth characteristics in diverse environments. MTH genome comprised of 93 unique genes-coding for starvation survival and stress-response proteins. These proteins confer its adaptation to nutritional deprivation and other abiotic stresses. MTH has a typical system to withstand its growth and cell viability during stable operation and recovery after prolonged starvation. Thus, the present work will provide an insight to improve the genome refinement and metabolic reconstruction in parallel to other closely related species.

Microbial vanadate and nitrate reductions coupled with anaerobic methane oxidation in groundwater

Vanadate contaminant in groundwater receives increasing attentions, but little is known on its biogeochemical transformation with gaseous electron donors. This study investigated bio-reduction of vanadate coupled with anaerobic methane oxidation and its relationship with nitrate reduction. Results showed 95.8 ± 3.1% of 1 mM vanadate was removed within 7 days using methane as the sole electron donor. Tetravalent vanadium compounds were the main reduction products, which precipitated naturally in groundwater environment. The introduction of nitrate inhibited vanadate reduction, though both were reduced in parallel. Accumulations of volatile fatty acids (VFAs) were observed from methane oxidation. Preliminary microbial community structure and metabolite analyses indicated that vanadate was likely reduced via Methylomonas coupled with methane oxidation or through synergistic relationships between methane oxidizing bacteria and heterotrophic vanadate reducers with VFAs served as the intermediates.

Concerted Metabolic Shifts Give New Insights Into the Syntrophic Mechanism Between Propionate-Fermenting Pelotomaculum thermopropionicum and Hydrogenotrophic Methanocella conradii

Microbial syntrophy is a thermodynamically-based cooperation between microbial partners that share the small amounts of free energy for anaerobic growth. To gain insights into the mechanism by which syntrophic microorganisms coordinate their metabolism, we constructed cocultures of propionate-oxidizing Pelotomaculum thermopropionicum and hydrogenotrophic Methanocella conradii and compared them to monocultures. Transcriptome analysis was performed on these cultures using strand-specific mRNA sequencing (RNA-Seq). The results showed that in coculture both P. thermopropionicum and M. conradii significantly upregulated the expression of genes involved in catabolism but downregulated those for anabolic biosynthesis. Specifically, genes coding for the methylmalonyl-CoA pathway in P. thermopropionicum and key genes for methanogenesis in M. conradii were substantially upregulated in coculture compared to monoculture. The putative flavin-based electron bifurcation/confurcation systems in both organisms were also upregulated in coculture. Formate dehydrogenase encoding genes in both organisms were markedly upregulated, indicating that formate was produced and utilized by P. thermopropionicum and M. conradii, respectively. The inhibition of syntrophic activity by formate and 2-bromoethanesulphonate (2-BES) but not H₂/CO₂ also suggested that formate production was used by P. thermopropionicum for the recycling of intracellular redox mediators. Finally, flagellum-induced signal transduction and amino acids exchange was upregulated for syntrophic interactions. Together, our study suggests that syntrophic organisms employ multiple strategies including global metabolic shift, utilization of electron bifurcation/confurcation and employing formate as an alternate electron carrier to optimize their metabolisms for syntrophic growth.

Discovery of a new metal and NAD+-dependent formate dehydrogenase from Clostridium ljungdahlii

Over the next decades, with the growing concern of rising atmospheric carbon dioxide (CO₂) levels, the importance of investigating new approaches for its reduction becomes crucial. Reclamation of CO₂ for conversion into biofuels represents an alternative and attractive production method that has been studied in recent years, now with enzymatic methods gaining more attention. Formate dehydrogenases (FDHs) are NAD(P)H-dependent oxidoreductases that catalyze the conversion of formate into CO₂ and have been extensively used for cofactor recycling in chemoenzymatic processes. A new FDH from Clostridium ljungdahlii (ClFDH) has been recently shown to possess activity in the reverse reaction: the mineralization of CO₂ into formate. In this study, we show the successful homologous expression of ClFDH in Escherichia coli. Biochemical and kinetic characterization of the enzyme revealed that this homologue also demonstrates activity toward CO₂ reduction. Structural analysis of the enzyme through homology modeling is also presented.

Anion-selective Formate/nitrite transporters: taxonomic distribution, phylogenetic analysis and subfamily-specific conservation pattern in prokaryotes

Background

The monovalent anions formate, nitrite and hydrosulphide are main metabolites of bacterial respiration during anaerobic mixed-acid fermentation. When accumulated in the cytoplasm, these anions become cytotoxic. Membrane proteins that selectively transport these monovalent anions across the membrane have been identified and they belong to the family of Formate/Nitrite Transporters (FNTs). Individual members that selectively transport formate, nitrite and hydrosulphide have been investigated. Experimentally determined structures of FNTs indicate that they share the same hourglass helical fold with aquaporins and aquaglyceroporins and have two constriction regions, namely, cytoplasmic slit and central constriction. Members of FNTs are found in bacteria, archaea, fungi and protists. However, no FNT homolog has been identified in mammals. With FNTs as potential drug targets for many bacterial diseases, it is important to understand the mechanism of selectivity and transport across these transporters.

Results

We have systematically searched the sequence databases and identified 2206 FNT sequences from bacteria, archaea and eukaryotes. Although FNT sequences are very diverse, homology modeling followed by structure-based sequence alignment revealed that nearly one third of all the positions within the transmembrane region exhibit high conservation either as a group or at the level of individual residues across all three kingdoms. Phylogenetic analysis of prokaryotic FNT sequences revealed eight different subgroups. Formate, nitrite and hydrosulphide transporters respectively are clustered into two (FocA and FdhC), three (NirC-α, NirC-β and NirC-γ) and one (HSC) subfamilies. We have also recognized two FNT subgroups (YfdC-α and YfdC-β) with unassigned function. Analysis of taxonomic distribution indicates that each subfamily prefers specific taxonomic groups. Structure-based sequence alignment of individual subfamily members revealed that certain positions in the two constriction regions and some residues facing the interior show subfamily-specific conservation. We have also identified examples of FNTs with the two constriction regions formed by residues that are less frequently observed. We have developed dbFNT, a database of FNT models and associated details. dbFNT is freely available to scientific community.

Conclusions

Taxonomic distribution and sequence conservation of FNTs exhibit subfamily-specific features. The conservation pattern in the central constriction and cytoplasmic slit in the open and closed states are distinct for YfdC and NirC subfamilies. The same is true for some residues facing the interior of the transporters. The specific residues in these positions can exert influence on the type of solutes that are transported by these proteins. With FNTs found in many disease-causing bacteria, the knowledge gained in this study can be used in the development and design of anti-bacterial drugs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3947-4) contains supplementary material, which is available to authorized users.

Intergenomic evolution and metabolic cross-talk between rumen and thermophilic autotrophic methanogenic archaea

Methanobrevibacter ruminantium M1 (MRU) is a rumen methanogenic archaean that can be able to utilize formate and CO₂/H₂ as growth substrates. Extensive analysis on the evolutionary genomic contexts considered herein to unravel its intergenomic relationship and metabolic adjustment acquired from the genomic content of Methanothermobacter thermautotrophicus ΔH. We demonstrated its intergenomic distance, genome function, synteny homologs and gene families, origin of replication, and methanogenesis to reveal the evolutionary relationships between Methanobrevibacter and Methanothermobacter. Comparison of the phylogenetic and metabolic markers was suggested for its archaeal metabolic core lineage that might have evolved from Methanothermobacter. Orthologous genes involved in its hydrogenotrophic methanogenesis might be acquired from intergenomic ancestry of Methanothermobacter via Methanobacterium formicicum. Formate dehydrogenase (fdhAB) coding gene cluster and carbon monoxide dehydrogenase (cooF) coding gene might have evolved from duplication events within Methanobrevibacter-Methanothermobacter lineage, and fdhCD gene cluster acquired from bacterial origins. Genome-wide metabolic survey found the existence of four novel pathways viz. l-tyrosine catabolism, mevalonate pathway II, acyl-carrier protein metabolism II and glutathione redox reactions II in MRU. Finding of these pathways suggested that MRU has shown a metabolic potential to tolerate molecular oxygen, antimicrobial metabolite biosynthesis and atypical lipid composition in cell wall, which was acquainted by metabolic cross-talk with mammalian bacterial origins. We conclude that coevolution of genomic contents between Methanobrevibacter and Methanothermobacter provides a clue to understand the metabolic adaptation of MRU in the rumen at different environmental niches.

Response of a rice paddy soil methanogen to syntrophic growth as revealed by transcriptional analyses

Members of Methanocellales are widespread in paddy field soils and play the key role in methane production. These methanogens feature largely in these organisms’ adaptation to low H2 and syntrophic growth with anaerobic fatty acid oxidizers. The adaptive mechanisms, however, remain unknown. In the present study, we determined the transcripts of 21 genes involved in the key steps of methanogenesis and acetate assimilation of Methanocella conradii HZ254, a strain recently isolated from paddy field soil. M. conradii was grown in monoculture and syntrophically with Pelotomaculum thermopropionicum (a propionate syntroph) or Syntrophothermus lipocalidus (a butyrate syntroph). Comparison of the relative transcript abundances showed that three hydrogenase-encoding genes and all methanogenesis-related genes tested were upregulated in cocultures relative to monoculture. The genes encoding formylmethanofuran dehydrogenase (Fwd), heterodisulfide reductase (Hdr), and the membrane-bound energy-converting hydrogenase (Ech) were the most upregulated among the evaluated genes. The expression of the formate dehydrogenase (Fdh)-encoding gene also was significantly upregulated. In contrast, an acetate assimilation gene was downregulated in cocultures. The genes coding for Fwd, Hdr, and the D subunit of F420-nonreducing hydrogenase (Mvh) form a large predicted transcription unit; therefore, the Mvh/Hdr/Fwd complex, capable of mediating the electron bifurcation and connecting the first and last steps of methanogenesis, was predicted to be formed in M. conradii. We propose that Methanocella methanogens cope with low H2 and syntrophic growth by (i) stabilizing the Mvh/Hdr/Fwd complex and (ii) activating formatedependent methanogenesis.

Assembly and catalysis of molybdenum or tungsten-containing formate dehydrogenases from bacteria

The global carbon cycle depends on the biological transformations of C1 compounds, which include the reductive incorporation of CO₂into organic molecules (e.g. in photosynthesis and other autotrophic pathways), in addition to the production of CO₂from formate, a reaction that is catalyzed by formate dehydrogenases (FDHs). FDHs catalyze, in general, the oxidation of formate to CO₂and H⁺. However, selected enzymes were identified to act as CO₂reductases, which are able to reduce CO₂to formate under physiological conditions. This reaction is of interest for the generation of formate as a convenient storage form of H₂for future applications. Cofactor-containing FDHs are found in anaerobic bacteria and archaea, in addition to facultative anaerobic or aerobic bacteria. These enzymes are highly diverse and employ different cofactors such as the molybdenum cofactor (Moco), FeS clusters and flavins, or cytochromes. Some enzymes include tungsten (W) in place of molybdenum (Mo) at the active site. For catalytic activity, a selenocysteine (SeCys) or cysteine (Cys) ligand at the Mo atom in the active site is essential for the reaction. This review will focus on the characterization of Mo- and W-containing FDHs from bacteria, their active site structure, subunit compositions and its proposed catalytic mechanism. We will give an overview on the different mechanisms of substrate conversion available so far, in addition to providing an outlook on bio-applications of FDHs. This article is part of a Special Issue entitled: Cofactor-dependent proteins: evolution, chemical diversity and bio-applications.

Physiological and genetic basis for self-aggregation of a thermophilic hydrogenotrophic methanogen, Methanothermobacter strain CaT2

Several thermophilic hydrogenotrophic methanogens naturally aggregate in their habitats in association with hydrogen-producing bacteria for efficient transfer of the methane fermentation intermediates to produce methane. However, physiology of aggregation and the identity of aggregation-specific genes remain to be elucidated. Here, we isolated and characterized a hydrogen and formate-utilizing Methanothermobacter sp. CaT2 that is capable of self-aggregation and utilizing formate. CaT2 produced methane from propionate oxidation in association with a syntrophic propionate-oxidizing bacterium faster than other methanogens, including Methanothermobacter thermautotrophicus ΔH and Methanothermobacter thermautotrophicus Z-245. CaT2 also aggregated throughout the culture period and was coated with polysaccharides, which was not found on the ΔH and Z-245 cells. Sugar content (particularly of rhamnose and mannose) was also higher in the CaT2 cells than the ΔH and Z-245 cells. Comparative genomic analysis of CaT2 indicated that four candidate genes, all of which encode glycosyltransferase, were involved in aggregation of CaT2. Transcriptional analysis showed that one glycosyltransferase gene was expressed at relatively high levels under normal growth conditions. The polysaccharide layer on the CaT2 cell surface, which is probably assembled by these glycosyltransferases, may be involved in cell aggregation.

Concerted movement in pH-dependent gating of FocA from molecular dynamics simulations

FocA, a member of the formate-nitrite transporter (FNT) family, transports formate and nitrite across biological membranes in cellular organisms. The export and uptake of formate in bacteria are both mediated by FocA, which undergoes a pH-dependent functional switch. Recently, the crystal structures of Escherichia coli FocA (EcFocA), Vibrio cholerae FocA (VcFocA), and Salmonella typhimurium FocA (StFocA) were reported. We performed molecular dynamics (MD) on StFocA and EcFocA with different states of His209 (protonated and unprotonated), representing different pH conditions of FocA. The N-terminal helix in each protomer of StFocA covers and blocks the formate channel. At neutral or high pH (MD simulations with unprotonated His209), the concerted movement of the N-terminal helices of pairs of protomers of StFocA opens its formate channel. At low pH (MD simulations with protonated His209), protonated His209 interacts tightly with its neighboring residue Asn262, and the channel becomes narrower, so that the formate can hardly pass through the channel. We obtained similar results for EcFocA. Our study shows that pairs of protomers of FocA move in a concerted way to achieve its pH-dependent gating function, which provides information on the dynamics of the gating mechanism of FNT proteins and aquaporins.

ColorPhylo: A Color Code to Accurately Display Taxonomic Classifications

Color may be very useful to visualise complex data. As far as taxonomy is concerned, color may help observing various species' characteristics in correlation with classification. However, choosing the number of subclasses to display is often a complex task: on the one hand, assigning a limited number of colors to taxa of interest hides the structure imbedded in the subtrees of the taxonomy; on the other hand, differentiating a high number of taxa by giving them specific colors, without considering the underlying taxonomy, may lead to unreadable results since relationships between displayed taxa would not be supported by the color code. In the present paper, an automatic color coding scheme is proposed to visualise the levels of taxonomic relationships displayed as overlay on any kind of data plot. To achieve this goal, a dimensionality reduction method allows displaying taxonomic "distances" onto a Euclidean two-dimensional space. The resulting map is projected onto a 2D color space (the Hue, Saturation, Brightness colorimetric space with brightness set to 1). Proximity in the taxonomic classification corresponds to proximity on the map and is therefore materialised by color proximity. As a result, each species is related to a color code showing its position in the taxonomic tree. The so called ColorPhylo displays taxonomic relationships intuitively and can be combined with any biological result. A Matlab version of ColorPhylo is available at http://sy.lespi.free.fr/ColorPhylo-homepage.html. Meanwhile, an ad-hoc distance in case of taxonomy with unknown edge lengths is proposed.

Structure of the formate transporter FocA reveals a pentameric aquaporin-like channel

FocA is a representative member of the formate-nitrite transporter family, which transports short-chain acids in bacteria, archaea, fungi, algae and parasites. The structure and transport mechanism of the formate-nitrite transporter family remain unknown. Here we report the crystal structure of Escherichia coli FocA at 2.25 A resolution. FocA forms a symmetric pentamer, with each protomer consisting of six transmembrane segments. Despite a lack of sequence homology, the overall structure of the FocA protomer closely resembles that of aquaporin and strongly argues that FocA is a channel, rather than a transporter. Structural analysis identifies potentially important channel residues, defines the channel path and reveals two constriction sites. Unlike aquaporin, FocA is impermeable to water but allows the passage of formate. A structural and biochemical investigation provides mechanistic insights into the channel activity of FocA.

Functionally distinct genes regulated by hydrogen limitation and growth rate in methanogenic Archaea

The use of molecular hydrogen as electron donor for energy generation is a defining characteristic of the hydrogenotrophic methanogens, an ancient group that dominates the phylum Euryarchaeota. We present here a global study of changes in mRNA abundance in response to hydrogen availability for a hydrogenotrophic methanogen. Cells of Methanococcus maripaludis were grown by using continuous culture to deconvolute the effects of hydrogen limitation and growth rate, and microarray analyses were conducted. Hydrogen limitation markedly increased mRNA levels for genes encoding enzymes of the methanogenic pathway that reduce or oxidize the electron-carrying deazaflavin, coenzyme F(420). F(420)-dependent redox functions in energy-generating metabolism are characteristic of the methanogenic Archaea, and the results show that their regulation is distinct from other redox processes in the cell. Rapid growth increased mRNA levels of the gene for an unusual hydrogenase, the hydrogen-dependent methylenetetrahydromethanopterin dehydrogenase.

Methanol-dependent gene expression demonstrates that methyl-coenzyme M reductase is essential in Methanosarcina acetivorans C2A and allows isolation of mutants with defects in regulation of the methanol utilization pathway

Methanosarcina acetivorans C2A is able to convert several substrates to methane via at least four distinct methanogenic pathways. A common step in each of these pathways is the reduction of methyl-coenzyme M (CoM) to methane catalyzed by methyl-CoM reductase (MCR). Because this enzyme is used in each of the known pathways, the mcrBDCGA operon, which encodes MCR, is expected to be essential. To validate this prediction, a system for conditional gene inactivation was developed. A heterologous copy of the mcrBDCGA operon was placed under the control of the highly regulated mtaC1 promoter, which directs the expression of genes involved in methanol utilization, and recombined onto the M. acetivorans chromosome. This allowed for disruption of the endogenous mcr operon in the presence of methanol. Because the PmtaC1 promoter is transcribed only during growth on methanol, mcrBDCGA was rendered methanol dependent and the strain was unable to grow in trimethylamine media, strongly suggesting that mcrBDCGA is essential. Upon prolonged incubation, suppressed mutants which expressed mcrBDCGA constitutively could be selected. Expression analysis of PmtaC1::uidA gene fusions in several isolated suppressed mutants suggests that they carry trans-active mutations leading to deregulation of all genes under control of this promoter. Subsequently, proteome analysis of one such suppressed mutant revealed that all known proteins derived from mtaC1 promoter-dependent expression were constitutively expressed in this mutant. This genetic system can therefore be employed for the testing of essential genes and for the identification of genes under a common regulatory mechanism by making regulatory mutations phenotypically selectable.

Expression profiling-based identification of CO2-responsive genes regulated by CCM1 controlling a carbon-concentrating mechanism in Chlamydomonas reinhardtii

Photosynthetic acclimation to CO2-limiting stress is associated with control of genetic and physiological responses through a signal transduction pathway, followed by integrated monitoring of the environmental changes. Although several CO2-responsive genes have been previously isolated, genome-wide analysis has not been applied to the isolation of CO2-responsive genes that may function as part of a carbon-concentrating mechanism (CCM) in photosynthetic eukaryotes. By comparing expression profiles of cells grown under CO2-rich conditions with those of cells grown under CO2-limiting conditions using a cDNA membrane array containing 10,368 expressed sequence tags, 51 low-CO2 inducible genes and 32 genes repressed by low CO2 whose mRNA levels were changed more than 2.5-fold in Chlamydomonas reinhardtii Dangeard were detected. The fact that the induction of almost all low-CO2 inducible genes was impaired in the ccm1 mutant suggests that CCM1 is a master regulator of CCM through putative low-CO2 signal transduction pathways. Among low-CO2 inducible genes, two novel genes, LciA and LciB, were identified, which may be involved in inorganic carbon transport. Possible functions of low-CO2 inducible and/or CCM1-regulated genes are discussed in relation to the CCM.

Function and regulation of the formate dehydrogenase genes of the methanogenic archaeon Methanococcus maripaludis

Methanococcus maripaludis is a mesophilic species of Archaea capable of producing methane from two substrates: hydrogen plus carbon dioxide and formate. To study the latter, we identified the formate dehydrogenase genes of M. maripaludis and found that the genome contains two gene clusters important for formate utilization. Phylogenetic analysis suggested that the two formate dehydrogenase gene sets arose from duplication events within the methanococcal lineage. The first gene cluster encodes homologs of formate dehydrogenase alpha (FdhA) and beta (FdhB) subunits and a putative formate transporter (FdhC) as well as a carbonic anhydrase analog. The second gene cluster encodes only FdhA and FdhB homologs. Mutants lacking either fdhA gene exhibited a partial growth defect on formate, whereas a double mutant was completely unable to grow on formate as a sole methanogenic substrate. Investigation of fdh gene expression revealed that transcription of both gene clusters is controlled by the presence of H(2) and not by the presence of formate.

Molecular and biochemical characterization of two tungsten- and selenium-containing formate dehydrogenases from Eubacterium acidaminophilum that are associated with components of an iron-only hydrogenase

Two gene clusters encoding similar formate dehydrogenases (FDH) were identified in Eubacterium acidaminophilum. Each cluster is composed of one gene coding for a catalytic subunit ( fdhA-I, fdhA-II) and one for an electron-transferring subunit ( fdhB-I, fdhB-II). Both fdhA genes contain a TGA codon for selenocysteine incorporation and the encoded proteins harbor five putative iron-sulfur clusters in their N-terminal region. Both FdhB subunits resemble the N-terminal region of FdhA on the amino acid level and contain five putative iron-sulfur clusters. Four genes thought to encode the subunits of an iron-only hydrogenase are located upstream of the FDH gene cluster I. By sequence comparison, HymA and HymB are predicted to contain one and four iron-sulfur clusters, respectively, the latter protein also binding sites for FMN and NAD(P). Thus, HymA and HymB seem to represent electron-transferring subunits, and HymC the putative catalytic subunit containing motifs for four iron-sulfur clusters and one H-cluster specific for Fe-only hydrogenases. HymD has six predicted transmembrane helices and might be an integral membrane protein. Viologen-dependent FDH activity was purified from serine-grown cells of E. acidaminophilum and the purified protein complex contained four subunits, FdhA and FdhB, encoded by FDH gene cluster II, and HymA and HymB, identified after determination of their N-terminal sequences. Thus, this complex might represent the most simple type of a formate hydrogen lyase. The purified formate dehydrogenase fraction contained iron, tungsten, a pterin cofactor, and zinc, but no molybdenum. FDH-II had a two-fold higher K(m) for formate (0.37 mM) than FDH-I and also catalyzed CO(2) reduction to formate. Reverse transcription (RT)-PCR pointed to increased expression of FDH-II in serine-grown cells, supporting the isolation of this FDH isoform. The fdhA-I gene was expressed as inactive protein in Escherichia coli. The in-frame UGA codon for selenocysteine incorporation was read in the heterologous system only as stop codon, although its potential SECIS element exhibited a quite high similarity to that of E. coli FDH.

Nitrite transport to the chloroplast in Chlamydomonas reinhardtii: molecular evidence for a regulated process

Nitrite transport to the chloroplast is not a well documented process in spite of being a central step in the nitrate assimilation pathway. The lack of molecular evidence, as well as the easy diffusion of nitrite through biological membranes, have made this physiological process difficult to understand in plant nutrition. The aim of this review is to illustrate that nitrite transport to the chloroplast is a regulated step, intimately related to the efficiency of nitrate utilization. In Chlamydomonas reinhardtii, the Nar1;1 gene has been shown to have this role in nitrate assimilation. NAR1;1 corresponds to a plastidic membrane transporter protein related to the bacterial formate/nitrite transporters. At least four Nar1 genes might exist in Chlamydomonas. The existence of orthologous Nar1 genes in plants is discussed.

An Lrp-like transcriptional regulator from the archaeon Pyrococcus furiosus is negatively autoregulated

The archaeal transcriptional initiation machinery closely resembles core elements of the eukaryal polymerase II system. However, apart from the established basal archaeal transcription system, little is known about the modulation of gene expression in archaea. At present, no obvious eukaryal-like transcriptional regulators have been identified in archaea. Instead, we have previously isolated an archaeal gene, the Pyrococcus furiosus lrpA, that potentially encodes a bacterial-like transcriptional regulator. In the present study, we have for the first time addressed the actual involvement of an archaeal Lrp homologue in transcription modulation. For that purpose, we have produced LrpA in Escherichia coli. In a cell-free P. furiosus transcription system we used wild-type and mutated lrpA promoter fragments to demonstrate that the purified LrpA negatively regulates its own transcription. In addition, gel retardation analyses revealed a single protein-DNA complex, in which LrpA appeared to be present in (at least) a tetrameric conformation. The location of the LrpA binding site was further identified by DNaseI and hydroxyl radical footprinting, indicating that LrpA binds to a 46-base pair sequence that overlaps the transcriptional start site of its own promoter. The molecular basis of the transcription inhibition by LrpA is discussed.

Prokaryotic carbonic anhydrases

Carbonic anhydrases catalyze the reversible hydration of CO(2) [CO(2)+H(2)Oright harpoon over left harpoon HCO(3)(-)+H(+)]. Since the discovery of this zinc (Zn) metalloenzyme in erythrocytes over 65 years ago, carbonic anhydrase has not only been found in virtually all mammalian tissues but is also abundant in plants and green unicellular algae. The enzyme is important to many eukaryotic physiological processes such as respiration, CO(2) transport and photosynthesis. Although ubiquitous in highly evolved organisms from the Eukarya domain, the enzyme has received scant attention in prokaryotes from the Bacteria and Archaea domains and has been purified from only five species since it was first identified in Neisseria sicca in 1963. Recent work has shown that carbonic anhydrase is widespread in metabolically diverse species from both the Archaea and Bacteria domains indicating that the enzyme has a more extensive and fundamental role in prokaryotic biology than previously recognized. A remarkable feature of carbonic anhydrase is the existence of three distinct classes (designated alpha, beta and gamma) that have no significant sequence identity and were invented independently. Thus, the carbonic anhydrase classes are excellent examples of convergent evolution of catalytic function. Genes encoding enzymes from all three classes have been identified in the prokaryotes with the beta and gamma classes predominating. All of the mammalian isozymes (including the 10 human isozymes) belong to the alpha class; however, only nine alpha class carbonic anhydrase genes have thus far been found in the Bacteria domain and none in the Archaea domain. The beta class is comprised of enzymes from the chloroplasts of both monocotyledonous and dicotyledonous plants as well as enzymes from phylogenetically diverse species from the Archaea and Bacteria domains. The only gamma class carbonic anhydrase that has thus far been isolated and characterized is from the methanoarchaeon Methanosarcina thermophila. Interestingly, many prokaryotes contain carbonic anhydrase genes from more than one class; some even contain genes from all three known classes. In addition, some prokaryotes contain multiple genes encoding carbonic anhydrases from the same class. The presence of multiple carbonic anhydrase genes within a species underscores the importance of this enzyme in prokaryotic physiology; however, the role(s) of this enzyme is still largely unknown. Even though most of the information known about the function(s) of carbonic anhydrase primarily relates to its role in cyanobacterial CO(2) fixation, the prokaryotic enzyme has also been shown to function in cyanate degradation and the survival of intracellular pathogens within their host. Investigations into prokaryotic carbonic anhydrase have already led to the identification of a new class (gamma) and future research will undoubtedly reveal novel functions for carbonic anhydrase in prokaryotes.

Biochemical and molecular characterization of the [NiFe] hydrogenase from the hyperthermophilic archaeon, Thermococcus litoralis

Thermococcus litoralis is a hyperthermophilic archaeon that grows at temperatures up to 98 degrees C by fermentative metabolism and reduces elemental sulfur (S0) to H2S. A [NiFe] hydrogenase, responsible for H2S or H2 production, has been purified and characterized. The enzyme is composed of four subunits with molecular mass 46, 42, 34 and 32 kDa. Elemental analyses gave approximate values of 22 Fe, 22 S and 1 Ni per hydrogenase. EPR spectra at 70 and 5 K indicated the presence of four or five [4Fe-4S] and one [2Fe-2S] type clusters. The optimal temperature for both H2 evolution and oxidation, using artificial electron carriers, was around 80 degrees C. The operon encoding the T. litoralis enzyme is composed of four genes forming one transcriptional unit, and transcription is not regulated by S0. An unusual transcription-initiation site is located 139 bp upstream from the translational start point. Sequence analyses indicated the presence of new putative nucleotide-binding domains. Upstream from the hydrogenase operon, ORFs probably encoding a molybdopterin oxidoreductase enzyme have been identified. Based on sequence, biochemical and biophysical analyses, a model of the enzyme and the pathway of electron flow during catalysis is proposed.

Transcriptional regulation in the hyperthermophilic archaeon Pyrococcus furiosus: coordinated expression of divergently oriented genes in response to beta-linked glucose polymers

The genetic organization, expression, and regulation of the celB locus of the hyperthermophilic archaeon Pyrococcus furiosus were analyzed. This locus includes the celB gene, which codes for an intracellular beta-glucosidase, and a divergently orientated gene cluster, adhA-adhB-lamA, which codes for two alcohol dehydrogenases and an extracellular beta-1,3-endoglucanase that is transcribed as a polycistronic messenger (the lamA operon). During growth of P. furiosus on either the beta-1,4-linked glucose dimer cellobiose or the beta-1,3-linked glucose polymer laminarin, the activities of both beta-glucosidase and endoglucanase were increased at least fivefold compared with levels during growth on maltose or pyruvate. Northern blot analysis revealed an enhanced transcription of both the celB gene and the lamA operon in the presence of these glucose-containing substrates. The in vivo and in vitro transcription initiation sites of both the celB gene and the lamA operon were identified 25 nucleotides downstream of conserved TATA box motifs. A number of repeating sequences have been recognized in the celB-adhA intergenic region, some of which might be part of a transcriptional regulator-binding site.

Phylogenetic characterization of novel transport protein families revealed by genome analyses

As a result of recent genome sequencing projects as well as detailed biochemical, molecular genetic and physiological experimentation on representative transport proteins, we have come to realize that all organisms possess an extensive but limited array of transport protein types that allow the uptake of nutrients and excretion of toxic substances. These proteins fall into phylogenetic families that presumably reflect their evolutionary histories. Some of these families are restricted to a single phylogenetic group of organisms and may have arisen recently in evolutionary time while others are found ubiquitously and may be ancient. In this study we conduct systematic phylogenetic analyses of 26 families of transport systems that either had not been characterized previously or were in need of updating. Among the families analyzed are some that are bacterial-specific, others that are eukaryotic-specific, and others that are ubiquitous. They can function by either a channel-type or a carrier-type mechanism, and in the latter case, they are frequently energized by coupling solute transport to the flux of an ion down its electrochemical gradient. We tabulate the currently sequenced members of the 26 families analyzed, describe the properties of these families, and present partial multiple alignments, signature sequences and phylogenetic trees for them all.

Biochemistry of methanogenesis: a tribute to Marjory Stephenson. 1998 Marjory Stephenson Prize Lecture

Max-Planck-Institut für terrestrische Mikrobiologie, Karl-von-Frisch-Straße, D-35043 Marburg, and Laboratorium für Mikrobiologie, Fachbereich Biologie, Philipps-Universität, Karl-von-Frisch-Straße, D-35032 Marburg, Germany

In 1933, Stephenson & Stickland (1933a) published that they had isolated from river mud, by the single cell technique, a methanogenic organism capable of growth in an inorganic medium with formate as the sole carbon source.

Clustered genes encoding the methyltransferases of methanogenesis from monomethylamine

Coenzyme M (CoM) is methylated during methanogenesis from monomethyamine in a reaction catalyzed by three proteins. Using monomethylamine, a 52-kDa polypeptide termed monomethylamine methyltransferase (MMAMT) methylates the corrinoid cofactor bound to a second polypeptide, monomethylamine corrinoid protein (MMCP). Methylated MMCP then serves as a substrate for MT2-A, which methylates CoM. The genes for these proteins are clustered on 6.8 kb of DNA in Methanosarcina barkeri MS. The gene encoding MMCP (mtmC) is located directly upstream of the gene encoding MMAMT (mtmB). The gene encoding MT2-A (mtbA) was found 1.1 kb upstream of mtmC, but no obvious open reading frame was found in the intergenic region between mtbA and mtmC. A single monocistronic transcript was found for mtbA that initiated 76 bp from the translational start. Separate transcripts of 2.4 and 4.7 kb were detected, both of which carried mtmCB. The larger transcript also encoded mtmP, which is homologous to the APC family of cationic amine permeases and may therefore encode a methylamine permease. A single transcriptional start site was found 447 bp upstream of the translational start of mtmC. MtmC possesses the corrinoid binding motif found in corrinoid proteins involved in dimethylsulfide- and methanol-dependent methanogenesis, as well as in methionine synthase. The open reading frame of mtmB was interrupted by a single in-frame, midframe, UAG codon which was also found in mtmB from M. barkeri NIH. A mechanism that circumvents UAG-directed termination of translation must operate during expression of mtmB in this methanogen.

Complete genome sequence of Methanobacterium thermoautotrophicum deltaH: functional analysis and comparative genomics

The complete 1,751,377-bp sequence of the genome of the thermophilic archaeon Methanobacterium thermoautotrophicum deltaH has been determined by a whole-genome shotgun sequencing approach. A total of 1,855 open reading frames (ORFs) have been identified that appear to encode polypeptides, 844 (46%) of which have been assigned putative functions based on their similarities to database sequences with assigned functions. A total of 514 (28%) of the ORF-encoded polypeptides are related to sequences with unknown functions, and 496 (27%) have little or no homology to sequences in public databases. Comparisons with Eucarya-, Bacteria-, and Archaea-specific databases reveal that 1,013 of the putative gene products (54%) are most similar to polypeptide sequences described previously for other organisms in the domain Archaea. Comparisons with the Methanococcus jannaschii genome data underline the extensive divergence that has occurred between these two methanogens; only 352 (19%) of M. thermoautotrophicum ORFs encode sequences that are >50% identical to M. jannaschii polypeptides, and there is little conservation in the relative locations of orthologous genes. When the M. thermoautotrophicum ORFs are compared to sequences from only the eucaryal and bacterial domains, 786 (42%) are more similar to bacterial sequences and 241 (13%) are more similar to eucaryal sequences. The bacterial domain-like gene products include the majority of those predicted to be involved in cofactor and small molecule biosyntheses, intermediary metabolism, transport, nitrogen fixation, regulatory functions, and interactions with the environment. Most proteins predicted to be involved in DNA metabolism, transcription, and translation are more similar to eucaryal sequences. Gene structure and organization have features that are typical of the Bacteria, including genes that encode polypeptides closely related to eucaryal proteins. There are 24 polypeptides that could form two-component sensor kinase-response regulator systems and homologs of the bacterial Hsp70-response proteins DnaK and DnaJ, which are notably absent in M. jannaschii. DNA replication initiation and chromosome packaging in M. thermoautotrophicum are predicted to have eucaryal features, based on the presence of two Cdc6 homologs and three histones; however, the presence of an ftsZ gene indicates a bacterial type of cell division initiation. The DNA polymerases include an X-family repair type and an unusual archaeal B type formed by two separate polypeptides. The DNA-dependent RNA polymerase (RNAP) subunits A', A", B', B" and H are encoded in a typical archaeal RNAP operon, although a second A' subunit-encoding gene is present at a remote location. There are two rRNA operons, and 39 tRNA genes are dispersed around the genome, although most of these occur in clusters. Three of the tRNA genes have introns, including the tRNAPro (GGG) gene, which contains a second intron at an unprecedented location. There is no selenocysteinyl-tRNA gene nor evidence for classically organized IS elements, prophages, or plasmids. The genome contains one intein and two extended repeats (3.6 and 8.6 kb) that are members of a family with 18 representatives in the M. jannaschii genome.

Hydrogen regulation of growth, growth yields, and methane gene transcription in Methanobacterium thermoautotrophicum deltaH

Changes in growth rate, methanogenesis, growth yield (Y(CH4)), and methane gene transcription have been correlated with changes in the supply of H2 to Methanobacterium thermoautotrophicum deltaH cells growing on H2 plus CO2 in fed-batch cultures. Under conditions of excess H2, biomass and methanogenesis increased exponentially and in parallel, resulting in cultures with a constant Y(CH4) and transcription of the mth and mrt genes that encode the H2-dependent N5,N10-methenyltetrahydromethanopterin (methenyl-H4MPT) reductase (MTH) and methyl coenzyme M reductase II (MRII), respectively. Reducing the H2 supply, by decreasing the percentage of H2 in the input gas mixture or by reducing the mixing speed of the fermentor impeller, decreased the growth rate and resulted in lower and constant rates of methanogenesis. Under such H2-limited growth conditions, cultures grew with a continuously increasing Y(CH4) and the mtd and mcr genes that encode the reduced coenzyme F420-dependent N5,N10-methenyl-H4MPT reductase (MTD) and methyl coenzyme M reductase I (MRI), respectively, were transcribed. Changes in the kinetics of growth, methanogenesis, and methane gene transcription directed by reducing the H2 supply could be reversed by restoring a high H2 supply. Methane production continued, but at a low and constant rate, and only mcr transcripts could be detected when the H2 supply was reduced to a level insufficient for growth. ftsA transcripts, which encode coenzyme F390 synthetase, were most abundant in cells growing with high H2 availability, consistent with coenzyme F390 synthesis signaling a high exogenous supply of reductant.