Wide Distribution of Genes for Tetrahydromethanopterin/Methanofuran-Linked C1 Transfer Reactions Argues for Their Presence in the Common Ancestor of Bacteria and Archaea

Genomic profiling and characteristics of a C1 degrading heterotrophic fresh-water bacterium Paracoccus sp. strain DMF

Paracoccus species are metabolically versatile gram-negative, aerobic facultative methylotrophic bacteria showing enormous promise for environmental and bioremediation studies. Here we report, the complete genome analysis of Paracoccus sp. strain DMF (P. DMF) that was isolated from a domestic wastewater treatment plant in Kanpur, India (26.4287 °N, 80.3891 °E) based on its ability to degrade a recalcitrant organic solvent N, N-dimethylformamide (DMF). The results reveal a genome size of 4,202,269 base pairs (bp) with a G + C content of 67.9%. The assembled genome comprises 4141 coding sequences (CDS), 46 RNA sequences, and 2 CRISPRs. Interestingly, catabolic operons related to the conventional marine-based methylated amines (MAs) degradation pathway were functionally annotated within the genome of an obligated aerobic heterotroph that is P. DMF. The genomic data-based characterization presented here for the novel heterotroph P. DMF aims to improve the understanding of the phenotypic gene products, enzymes, and pathways involved with greater emphasis on facultative methylotrophic motility-based latent pathogenicity.

A diversified, widespread microbial gene cluster encodes homologs of methyltransferases involved in methanogenesis

Analyses of microbial genomes have revealed unexpectedly wide distributions of enzymes from specialized metabolism, including methanogenesis, providing exciting opportunities for discovery. Here, we identify a family of gene clusters (the type 1 mlp gene clusters (MGCs)) that encodes homologs of the soluble coenzyme M methyltransferases (SCMTs) involved in methylotrophic methanogenesis and is widespread in bacteria and archaea. Type 1 MGCs are expressed and regulated in medically, environmentally, and industrially important organisms, making them likely to be physiologically relevant. Enzyme annotation, analysis of genomic context, and biochemical experiments suggests these gene clusters play a role in methyl-sulfur and/or methyl-selenide metabolism in numerous anoxic environments, including the human gut microbiome, potentially impacting sulfur and selenium cycling in diverse, anoxic environments.

Structural diversity of the coenzyme methylofuran and identification of enzymes for the biosynthesis of its polyglutamate side chain

Methylofuran (MYFR) is a formyl-carrying coenzyme essential for the oxidation of formaldehyde in most methylotrophic bacteria. In Methylorubrum extorquens, MYFR contains a large and branched polyglutamate side chain of up to 24 glutamates. These glutamates play an essential role in interfacing the coenzyme with the formyltransferase/hydrolase complex, an enzyme that generates formate. To date, MYFR has not been identified in other methylotrophs, and it is unknown whether its structural features are conserved. Here, we examined nine bacterial strains for the presence and structure of MYFR using high-resolution liquid chromatography-mass spectrometry (LC-MS). Two of the strains produced MYFR as present in M. extorquens, while a modified MYFR containing tyramine instead of tyrosine in its core structure was detected in six strains. When M. extorquens was grown in the presence of tyramine, the compound was readily incorporated into MYFR, indicating that the biosynthetic enzymes are unable to discriminate tyrosine from tyramine. Using gene deletions in combination with LC-MS analyses, we identified three genes, orf5, orfY, and orf17 that are essential for MYFR biosynthesis. Notably, the orfY and orf5 mutants accumulated short MYFR intermediates with only one and two glutamates, respectively, suggesting that these enzymes catalyze glutamate addition. Upon homologous overexpression of orf5, a drastic increase in the number of glutamates in MYFR was observed (up to 40 glutamates), further corroborating the function of Orf5 as a glutamate ligase. We thus renamed OrfY and Orf5 to MyfA and MyfB to highlight that these enzymes are specifically involved in MYFR biosynthesis.

Metagenomic Profiling and Microbial Metabolic Potential of Perdido Fold Belt (NW) and Campeche Knolls (SE) in the Gulf of Mexico

The Gulf of Mexico (GoM) is a particular environment that is continuously exposed to hydrocarbon compounds that may influence the microbial community composition. We carried out a metagenomic assessment of the bacterial community to get an overall view of this geographical zone. We analyzed both taxonomic and metabolic markers profiles to explain how the indigenous GoM microorganims participate in the biogeochemical cycling. Two geographically distant regions in the GoM, one in the north-west (NW) and one in the south-east (SE) of the GoM were analyzed and showed differences in their microbial composition and metabolic potential. These differences provide evidence the delicate equilibrium that sustains microbial communities and biogeochemical cycles. Based on the taxonomy and gene groups, the NW are more oxic sediments than SE ones, which have anaerobic conditions. Both water and sediments show the expected sulfur, nitrogen, and hydrocarbon metabolism genes, with particularly high diversity of the hydrocarbon-degrading ones. Accordingly, many of the assigned genera were associated with hydrocarbon degradation processes, Nitrospira and Sva0081 were the most abundant in sediments, while Vibrio, Alteromonas, and Alcanivorax were mostly detected in water samples. This basal-state analysis presents the GoM as a potential source of aerobic and anaerobic hydrocarbon degradation genes important for the ecological dynamics of hydrocarbons and the potential use for water and sediment bioremediation processes.

Diversity, ecology and evolution of Archaea

Compared to bacteria, our knowledge of archaeal biology is limited. Historically, microbiologists have mostly relied on culturing and single-gene diversity surveys to understand Archaea in nature. However, only six of the 27 currently proposed archaeal phyla have cultured representatives. Advances in genomic sequencing and computational approaches are revolutionizing our understanding of Archaea. The recovery of genomes belonging to uncultured groups from the environment has resulted in the description of several new phyla, many of which are globally distributed and are among the predominant organisms on the planet. In this Review, we discuss how these genomes, together with long-term enrichment studies and elegant in situ measurements, are providing insights into the metabolic capabilities of the Archaea. We also debate how such studies reveal how important Archaea are in mediating an array of ecological processes, including global carbon and nutrient cycles, and how this increase in archaeal diversity has expanded our view of the tree of life and early archaeal evolution, and has provided new insights into the origin of eukaryotes.

Multidomain ribosomal protein trees and the planctobacterial origin of neomura (eukaryotes, archaebacteria)

Palaeontologically, eubacteria are > 3× older than neomura (eukaryotes, archaebacteria). Cell biology contrasts ancestral eubacterial murein peptidoglycan walls and derived neomuran N-linked glycoprotein coats/walls. Misinterpreting long stems connecting clade neomura to eubacteria on ribosomal sequence trees (plus misinterpreted protein paralogue trees) obscured this historical pattern. Universal multiprotein ribosomal protein (RP) trees, more accurate than rRNA trees, are taxonomically undersampled. To reduce contradictions with genically richer eukaryote trees and improve eubacterial phylogeny, we constructed site-heterogeneous and maximum-likelihood universal three-domain, two-domain, and single-domain trees for 143 eukaryotes (branching now congruent with 187-protein trees), 60 archaebacteria, and 151 taxonomically representative eubacteria, using 51 and 26 RPs. Site-heterogeneous trees greatly improve eubacterial phylogeny and higher classification, e.g. showing gracilicute monophyly, that many 'rDNA-phyla' belong in Proteobacteria, and reveal robust new phyla Synthermota and Aquithermota. Monoderm Posibacteria and Mollicutes (two separate wall losses) are both polyphyletic: multiple outer membrane losses in Endobacteria occurred separately from Actinobacteria; neither phylum is related to Chloroflexi, the most divergent prokaryotes, which originated photosynthesis (new model proposed). RP trees support an eozoan root for eukaryotes and are consistent with archaebacteria being their sisters and rooted between Filarchaeota (=Proteoarchaeota, including 'Asgardia') and Euryarchaeota sensu-lato (including ultrasimplified 'DPANN' whose long branches often distort trees). Two-domain trees group eukaryotes within Planctobacteria, and archaebacteria with Planctobacteria/Sphingobacteria. Integrated molecular/palaeontological evidence favours negibacterial ancestors for neomura and all life. Unique presence of key pre-neomuran characters favours Planctobacteria only as ancestral to neomura, which apparently arose by coevolutionary repercussions (explained here in detail, including RP replacement) of simultaneous outer membrane and murein loss. Planctobacterial C-1 methanotrophic enzymes are likely ancestral to archaebacterial methanogenesis and β-propeller-α-solenoid proteins to eukaryotic vesicle coats, nuclear-pore-complexes, and intraciliary transport. Planctobacterial chaperone-independent 4/5-protofilament microtubules and MamK actin-ancestors prepared for eukaryote intracellular motility, mitosis, cytokinesis, and phagocytosis. We refute numerous wrong ideas about the universal tree.

Multidomain ribosomal protein trees and the planctobacterial origin of neomura (eukaryotes, archaebacteria)

Palaeontologically, eubacteria are > 3× older than neomura (eukaryotes, archaebacteria). Cell biology contrasts ancestral eubacterial murein peptidoglycan walls and derived neomuran N-linked glycoprotein coats/walls. Misinterpreting long stems connecting clade neomura to eubacteria on ribosomal sequence trees (plus misinterpreted protein paralogue trees) obscured this historical pattern. Universal multiprotein ribosomal protein (RP) trees, more accurate than rRNA trees, are taxonomically undersampled. To reduce contradictions with genically richer eukaryote trees and improve eubacterial phylogeny, we constructed site-heterogeneous and maximum-likelihood universal three-domain, two-domain, and single-domain trees for 143 eukaryotes (branching now congruent with 187-protein trees), 60 archaebacteria, and 151 taxonomically representative eubacteria, using 51 and 26 RPs. Site-heterogeneous trees greatly improve eubacterial phylogeny and higher classification, e.g. showing gracilicute monophyly, that many 'rDNA-phyla' belong in Proteobacteria, and reveal robust new phyla Synthermota and Aquithermota. Monoderm Posibacteria and Mollicutes (two separate wall losses) are both polyphyletic: multiple outer membrane losses in Endobacteria occurred separately from Actinobacteria; neither phylum is related to Chloroflexi, the most divergent prokaryotes, which originated photosynthesis (new model proposed). RP trees support an eozoan root for eukaryotes and are consistent with archaebacteria being their sisters and rooted between Filarchaeota (=Proteoarchaeota, including 'Asgardia') and Euryarchaeota sensu-lato (including ultrasimplified 'DPANN' whose long branches often distort trees). Two-domain trees group eukaryotes within Planctobacteria, and archaebacteria with Planctobacteria/Sphingobacteria. Integrated molecular/palaeontological evidence favours negibacterial ancestors for neomura and all life. Unique presence of key pre-neomuran characters favours Planctobacteria only as ancestral to neomura, which apparently arose by coevolutionary repercussions (explained here in detail, including RP replacement) of simultaneous outer membrane and murein loss. Planctobacterial C-1 methanotrophic enzymes are likely ancestral to archaebacterial methanogenesis and β-propeller-α-solenoid proteins to eukaryotic vesicle coats, nuclear-pore-complexes, and intraciliary transport. Planctobacterial chaperone-independent 4/5-protofilament microtubules and MamK actin-ancestors prepared for eukaryote intracellular motility, mitosis, cytokinesis, and phagocytosis. We refute numerous wrong ideas about the universal tree.

Methylofuran is a prosthetic group of the formyltransferase/hydrolase complex and shuttles one-carbon units between two active sites

Methylotrophs play a crucial role in the global carbon cycle as they oxidize reduced one-carbon compounds such as methanol or methane to CO₂. The step-wise conversion of these substrates generally takes place while the one-carbon units are bound to diffusible carrier molecules. However, our crystal structure of the formyltransferase/hydrolase complex from Methylorubrum extorquens demonstrates that the one-carbon carrier methylofuran tightly binds to the enzyme via its extended and branched polyglutamate chain. A swinging motion of the coenzyme then shuttles one-carbon units between the two active sites of the enzyme complex over a distance of 50 Å. The structure further highlights how the bacterial formate generation system relates to the archaeal CO₂ fixation system.

Methylotrophy, the ability of microorganisms to grow on reduced one-carbon substrates such as methane or methanol, is a feature of various bacterial species. The prevailing oxidation pathway depends on tetrahydromethanopterin (H₄MPT) and methylofuran (MYFR), an analog of methanofuran from methanogenic archaea. Formyltransferase/hydrolase complex (Fhc) generates formate from formyl-H₄MPT in two consecutive reactions where MYFR acts as a carrier of one-carbon units. Recently, we chemically characterized MYFR from the model methylotroph Methylorubrum extorquens and identified an unusually long polyglutamate side chain of up to 24 glutamates. Here, we report on the crystal structure of Fhc to investigate the function of the polyglutamate side chain in MYFR and the relatedness of the enzyme complex with the orthologous enzymes in archaea. We identified MYFR as a prosthetic group that is tightly, but noncovalently, bound to Fhc. Surprisingly, the structure of Fhc together with MYFR revealed that the polyglutamate side chain of MYFR is branched and contains glutamates with amide bonds at both their α- and γ-carboxyl groups. This negatively charged and branched polyglutamate side chain interacts with a cluster of conserved positively charged residues of Fhc, allowing for strong interactions. The MYFR binding site is located equidistantly from the active site of the formyltransferase (FhcD) and metallo-hydrolase (FhcA). The polyglutamate serves therefore an additional function as a swinging linker to shuttle the one-carbon carrying amine between the two active sites, thereby likely increasing overall catalysis while decreasing the need for high intracellular MYFR concentrations.

Role of the malic enzyme in metabolism of the halotolerant methanotroph Methylotuvimicrobium alcaliphilum 20Z

The bacteria utilizing methane as a growth substrate (methanotrophs) are important constituents of the biosphere. Methanotrophs mitigate the emission of anthropogenic and natural greenhouse gas methane to the environment and are the promising agents for future biotechnologies. Many aspects of CH₄ bioconversion by methanotrophs require further clarification. This study was aimed at characterizing the biochemical properties of the malic enzyme (Mae) from the halotolerant obligate methanotroph Methylotuvimicrobium alcaliphilum 20Z. The His₆-tagged Mae was obtained by heterologous expression in Escherichia coli BL21 (DE3) and purified by affinity metal chelating chromatography. As determined by gel filtration and non-denaturating gradient gel electrophoresis, the molecular mass of the native enzyme is 260 kDa. The homotetrameric Mae (65x4 kDa) catalyzed an irreversible NAD⁺-dependent reaction of L-malate decarboxylation into pyruvate with a specific activity of 32 ± 2 units mg^-1 and K_m value of 5.5 ± 0.8 mM for malate and 57 ± 5 μM for NAD⁺. The disruption of the mae gene by insertion mutagenesis resulted in a 20-fold increase in intracellular malate level in the mutant compared to the wild type strain. Based on both enzyme and mutant properties, we conclude that the malic enzyme is involved in the control of intracellular L-malate level in Mtm. alcaliphilum 20Z. Genomic analysis has revealed that Maes present in methanotrophs fall into two different clades in the amino acid-based phylogenetic tree, but no correlation of the division with taxonomic affiliations of the host bacteria was observed.

Comparative genomics and transcriptomics insights into the C1 metabolic model of a formaldehyde-degrading strain Methylobacterium sp. XJLW

A formaldehyde-degrading strain Methylobacterium sp. XJLW was isolated and exhibited a special phenotype for formaldehyde utilization. The accumulation of formic acid in large quantities and lower cell growth was detected when XJLW utilized formaldehyde as the sole carbon source, suggesting XJLW has a potentially novel pathway to transfer formaldehyde to methanol and then enter the serine cycle for C1 metabolism. This mechanism requires exploration via molecular omics. Thus, the complete genome of XJLW was sequenced, and the transcriptome difference was also analyzed based on the RNA-seq data of strain XJLW cultivated with methanol and glucose, respectively. XJLW has a chromosome DNA and a mega-plasmid DNA. Ten percent of genes on chromosome DNA are strain-specific in genus Methylobacterium. Transcriptome analysis results showed that 623 genes were significantly up-regulated and that 207 genes were significantly down-regulated for growth in methanol. Among the up-regulated genes, 90 genes belong to strain-specific regions and are densely distributed in three areas. A specific gene (A3862_27225) annotated as methyltransferase was found ranking in the top 4 of up-regulated genes. This methyltransferase may play a role in the specific C1 metabolism of XJLW. Methylobacterium sp. XJLW should contain a potential methyl transport pathway via the novel methyltransferase, which is different from known pathways. These findings provide the basis for additional possibilities, which improve the formaldehyde-degrading ability of Methylobacterium sp. XJLW.

Methanogenesis on Early Stages of Life: Ancient but Not Primordial

Of the six known autotrophic pathways, the Wood-Ljungdahl pathway (WL) is the only one present in both the acetate producing Bacteria (homoacetogens) and the methane producing Archaea (hydrogenotrophic methanogens), and it has been suggested that WL is one of the oldest metabolic pathways. However, only the so-called carbonyl branch is shared by Archaea and Bacteria, while the methyl branch is different, both in the number of reactions and enzymes, which are not homologous among them. In this work we show that some parts of the methyl branch of archaeal Wood-Ljungdahl pathway (MBWL) are present in bacteria as well as in non-methanogen archaea, although the tangled evolutionary history of MBWL cannot be traced back to the Last Common Ancestor. We have also analyzed the different variants of methanogenesis (hydrogenotrophic, acetoclastic and methylotrophic pathways), and concluded that each of these pathways, and every different enzyme or subunit (in the case of multimeric enzymes), has their own intricate evolutionary history. Our study supports the scenario of hydrogenotrophic methanogenesis being older than the other variants, albeit not old enough to be present in the last archaeal common ancestor.

Evolutionary History of Copper Membrane Monooxygenases

Copper membrane monooxygenases (CuMMOs) oxidize ammonia, methane and some short-chain alkanes and alkenes. They are encoded by three genes, usually in an operon of xmoCAB. We aligned xmo operons from 66 microbial genomes, including members of the Alpha-, Beta-, and Gamma-proteobacteria, Verrucomicrobia, Actinobacteria, Thaumarchaeota and the candidate phylum NC10. Phylogenetic and compositional analyses were used to reconstruct the evolutionary history of the enzyme and detect potential lateral gene transfer (LGT) events. The phylogenetic analyses showed at least 10 clusters corresponding to a combination of substrate specificity and bacterial taxonomy, but with no overriding structure based on either function or taxonomy alone. Adaptation of the enzyme to preferentially oxidize either ammonia or methane has occurred more than once. Individual phylogenies of all three genes, xmoA, xmoB and xmoC, closely matched, indicating that this operon evolved or was consistently transferred as a unit, with the possible exception of the methane monooxygenase operons in Verrucomicrobia, where the pmoB gene has a distinct phylogeny from pmoA and pmoC. Compositional analyses indicated that some clusters of xmoCAB operons (for example, the pmoCAB in gammaproteobacterial methanotrophs and the amoCAB in betaproteobacterial nitrifiers) were compositionally very different from their genomes, possibly indicating recent lateral transfer of these operons. The combined phylogenetic and compositional analyses support the hypothesis that an ancestor of the nitrifying bacterium Nitrosococcus was the donor of methane monooxygenase (pMMO) to both the alphaproteobacterial and gammaproteobacterial methanotrophs, but that before this event the gammaproteobacterial methanotrophs originally possessed another CuMMO (Pxm), which has since been lost in many species.

Substrate Specificity Analysis of Dihydrofolate/Dihydromethanopterin Reductase Homologs in Methylotrophic α-Proteobacteria

Methane-producing archaea and methylotrophic bacteria use tetrahydromethanopterin (H₄MPT) and/or tetrahydrofolate (H₄F) as coenzymes in one-carbon (C1) transfer pathways. The α-proteobacterium Methylobacterium extorquens AM1 contains a dihydromethanopterin reductase (DmrA) and two annotated dihydrofolate reductases (DfrA and DfrB). DmrA has been shown to catalyze the final step of H₄MPT biosynthesis; however, the functions of DfrA and DfrB have not been examined biochemically. Moreover, sequence alignment (BLAST) searches have recognized scores of proteins that share up to 99% identity with DmrA but are annotated as diacylglycerol kinases (DAGK). In this work, we used bioinformatics and enzyme assays to provide insight into the phylogeny and substrate specificity of selected Dfr and DmrA homologs. In a phylogenetic tree, DmrA and homologs annotated as DAGKs grouped together in one clade. Purified histidine-tagged versions of the annotated DAGKs from Hyphomicrobium nitrativorans and M. nodulans (respectively, sharing 69 and 84% identity with DmrA) showed only low activity in phosphorylating 1,2-dihexanoyl-sn-glycerol when compared with a commercial DAGK from Escherichia coli. However, the annotated DAGKs successfully reduced a dihydromethanopterin analog (dihydrosarcinapterin, H₂SPT) with kinetic values similar to those determined for M. extorquens AM1 DmrA. DfrA and DfrB showed little or no ability to reduce H₂SPT under the conditions studied; however, both catalyzed the NADPH-dependent reduction of dihydrofolate. These results provide the first evidence that DfrA and DfrB function as authentic dihydrofolate reductases, while DAGKs with greater than 69% identity to DmrA may be misannotated and are likely to function in H₄MPT biosynthesis.

Current Trends in Methylotrophy

Methylotrophy is a field of study dealing with microorganisms capable of utilization of compounds devoid of carbon-carbon bonds (C1 compounds). In this review, we highlight several emerging trends in methylotrophy. First, we discuss the significance of the recent discovery of lanthanide-dependent alcohol dehydrogenases for understanding both the occurrence and the distribution of methylotrophy functions among bacteria, and then we discuss the newly appreciated role of lanthanides in biology. Next, we describe the detection of other methylotrophy pathways across novel bacterial taxa and insights into the evolution of methylotrophy. Further, data are presented on the occurrence and activity of aerobic methylotrophs in hypoxic and anoxic environments, questioning the prior assumptions on niche separation of aerobic and anaerobic methylotrophy. The concept of communal function in aerobic methane oxidation is also briefly discussed. Finally, we review recent research in engineering methylotrophs for biotechnological applications as well as recent progress in engineering synthetic methylotrophy.

Transposon Sequencing Uncovers an Essential Regulatory Function of Phosphoribulokinase for Methylotrophy

Methylotrophy is the ability of organisms to grow at the expense of reduced one-carbon compounds, such as methanol or methane. Here, we used transposon sequencing combining hyper-saturated transposon mutagenesis with high-throughput sequencing to define the essential methylotrophy genome of Methylobacterium extorquens PA1, a model methylotroph. To distinguish genomic regions required for growth only on methanol from general required genes, we contrasted growth on methanol with growth on succinate, a non-methylotrophic reference substrate. About 500,000 insertions were mapped for each condition, resulting in a median insertion distance of five base pairs. We identified 147 genes and 76 genes as specific for growth on methanol and succinate, respectively, and a set of 590 genes as required under both growth conditions. For the integration of metabolic functions, we reconstructed a genome-scale metabolic model and performed in silico essentiality analysis. In total, the approach uncovered 95 genes not previously described as crucial for methylotrophy, including genes involved in respiration, carbon metabolism, transport, and regulation. Strikingly, regardless of the absence of the Calvin cycle in the methylotroph, the screen led to the identification of the gene for phosphoribulokinase as essential during growth on methanol, but not during growth on succinate. Genetic experiments in addition to metabolomics and proteomics revealed that phosphoribulokinase serves a key regulatory function. Our data support a model according to which ribulose-1,5-bisphosphate is an essential metabolite that induces a transcriptional regulator driving one-carbon assimilation.