A late methanogen origin for molybdenum-dependent nitrogenase

Exploring the influence of atmospheric CO2 and O2 levels on the utility of nitrogen isotopes as proxy for biological N2 fixation

Biological N2 fixation (BNF) is traced to the Archean. The nitrogen isotopic fractionation composition (δ15N) of sedimentary rocks is commonly used to reconstruct the presence of ancient diazotrophic ecosystems. While δ15N has been validated mostly using organisms grown under present-day conditions; it has not under the pre-Cambrian conditions, when atmospheric pO2 was lower and pCO2 was higher. Here, we explore δ15N signatures under three atmospheres with (i) elevated CO2 and no O2 (Archean), (ii) present-day CO2, and O2 and (iii) future elevated CO2, in marine and freshwater, heterocytous cyanobacteria. Additionally, we augment our data set from literature for more generalized dependencies of δ15N and the associated fractionation factor epsilon (ε = δ15Nbiomass - δ15NN2) during BNF in Archaea and Bacteria, including cyanobacteria, and habitats. The ε ranges between 3.70‰ and -4.96‰ with a mean ε value of -1.38 ± 0.95‰, for all bacteria, including cyanobacteria, across all tested conditions. The expanded data set revealed correlations of isotopic fractionation of BNF with CO2 concentrations, toxin production, and light, although within 1‰. Moreover, correlation showed significant dependency of ε to species type, C/N ratios and toxin production in cyanobacteria, albeit it within a small range (-1.44 ± 0.89‰). We therefore conclude that δ15N is likely robust when applied to the pre-Cambrian-like atmosphere, stressing the strong cyanobacterial bias. Interestingly, the increased fractionation (lower ε) observed in the toxin-producing Nodularia and Nostoc spp. suggests a heretofore unknown role of toxins in modulating nitrogen isotopic signals that warrants further investigation.IMPORTANCENitrogen is an essential element of life on Earth; however, despite its abundance, it is not biologically accessible. Biological nitrogen fixation is an essential process whereby microbes fix N2 into biologically usable NH3. During this process, the enzyme nitrogenase preferentially uses light 14N, resulting in 15N depleted biomass. This signature can be traced back in time in sediments on Earth, and possibly other planets. In this paper, we explore the influence of pO2 and pCO2 on this fractionation signal. We find the signal is stable, especially for the primary producers, cyanobacteria, with correlations to CO2, light, and toxin-producing status, within a small range. Unexpectedly, we identified higher fractionation signals in toxin-producing Nodularia and Nostoc species that offer insight into why some organisms produce these N-rich toxic secondary metabolites.

A nitrogenase-like enzyme is involved in the novel anaerobic assimilation pathway of a sulfonate, isethionate, in the photosynthetic bacterium Rhodobacter capsulatus

Prokaryotes contribute to the global sulfur cycle by using diverse sulfur compounds as sulfur sources or electron acceptors. In this study, we report that a nitrogenase-like enzyme (NFL) and a radical SAM enzyme (RSE) are involved in the novel anaerobic assimilation pathway of a sulfonate, isethionate, in the photosynthetic bacterium Rhodobacter capsulatus. The nflHDK genes for NFL are localized at a locus containing genes for known sulfonate metabolism in the genome. A gene nflB encoding an RSE is present just upstream of nflH, forming a small gene cluster nflBHDK. Mutants lacking any nflBHDK genes are incapable of growing with isethionate as the sole sulfur source under anaerobic photosynthetic conditions, indicating that all four NflBHDK proteins are essential for the isethionate assimilation pathway. Heterologous expression of the islAB genes encoding a known isethionate lyase that degrades isethionate to sulfite and acetaldehyde restored the isethionate-dependent growth of a mutant lacking nflDK, indicating that the enzyme encoding nflBHDK is involved in an isethionate assimilation reaction to release sulfite. Furthermore, the heterologous expression of nflBHDK and ssuCAB encoding an isethionate transporter in the closely related species R. sphaeroides, which does not have nflBHDK and cannot grow with isethionate as the sole sulfur source, conferred isethionate-dependent growth ability to this species. We propose to rename nflBHDK as isrBHDK (isethionate reductase). The isrBHDK genes are widely distributed among various prokaryote phyla. Discovery of the isethionate assimilation pathway by IsrBHDK provides a missing piece for the anaerobic sulfur cycle and for understanding the evolution of ancient sulfur metabolism.IMPORTANCENitrogenase is an important enzyme found in prokaryotes that reduces atmospheric nitrogen to ammonia and plays a fundamental role in the global nitrogen cycle. It has been noted that nitrogenase-like enzymes (NFLs), which share an evolutionary origin with nitrogenase, have evolved to catalyze diverse reactions such as chlorophyll biosynthesis (photosynthesis), coenzyme F430 biosynthesis (methanogenesis), and methionine biosynthesis. In this study, we discovered that an NFL with unknown function in the photosynthetic bacterium Rhodobacter capsulatus is a novel isethionate reductase (Isr), which catalyzes the assimilatory degradation of isethionate, a sulfonate, releasing sulfite used as the sulfur source under anaerobic conditions. Isr is widely distributed among various bacterial phyla, including intestinal bacteria, and is presumed to play an important role in sulfur metabolism in anaerobic environments such as animal guts and microbial mats. This finding provides a clue for understanding ancient metabolism that evolved under anaerobic environments at the dawn of life.

Co-evolution of early Earth environments and microbial life

Two records of Earth history capture the evolution of life and its co-evolving ecosystems with interpretable fidelity: the geobiological and geochemical traces preserved in rocks and the evolutionary histories captured within genomes. The earliest vestiges of life are recognized mostly in isotopic fingerprints of specific microbial metabolisms, whereas fossils and organic biomarkers become important later. Molecular biology provides lineages that can be overlayed on geologic and geochemical records of evolving life. All these data lie within a framework of biospheric evolution that is primarily characterized by the transition from an oxygen-poor to an oxygen-rich world. In this Review, we explore the history of microbial life on Earth and the degree to which it shaped, and was shaped by, fundamental transitions in the chemical properties of the oceans, continents and atmosphere. We examine the diversity and evolution of early metabolic processes, their couplings with biogeochemical cycles and their links to the oxygenation of the early biosphere. We discuss the distinction between the beginnings of metabolisms and their subsequent proliferation and their capacity to shape surface environments on a planetary scale. The evolution of microbial life and its ecological impacts directly mirror the Earth's chemical and physical evolution through cause-and-effect relationships.

Structural comparison of (hyper-)thermophilic nitrogenase reductases from three marine Methanococcales

The nitrogenase reductase NifH catalyses ATP-dependent electron delivery to the Mo-nitrogenase, a reaction central to biological dinitrogen (N2) fixation. While NifHs have been extensively studied in bacteria, structural information about their archaeal counterparts is limited. Archaeal NifHs are considered more ancient, particularly those from Methanococcales, a group of marine hydrogenotrophic methanogens, which includes diazotrophs growing at temperatures near 92 °C. Here, we structurally and biochemically analyse NifHs from three Methanococcales, offering the X-ray crystal structures from meso-, thermo-, and hyperthermophilic methanogens. While NifH from Methanococcus maripaludis (37 °C) was obtained through heterologous recombinant expression, the proteins from Methanothermococcus thermolithotrophicus (65 °C) and Methanocaldococcus infernus (85 °C) were natively purified from the diazotrophic archaea. The structures from M. thermolithotrophicus crystallised as isolated exhibit high flexibility. In contrast, the complexes of NifH with MgADP obtained from the three methanogens are superposable, more rigid, and present remarkable structural conservation with their homologues. They retain key structural features of P-loop NTPases and share similar electrostatic profiles with the counterpart from the bacterial model organism Azotobacter vinelandii. In comparison to the NifH from the phylogenetically distant Methanosarcina acetivorans, these reductases do not cross-react significantly with Mo-nitrogenase from A. vinelandii. However, they associate with bacterial nitrogenase when ADP·AlF 4 - is added to mimic a transient reactive state. Accordingly, detailed surface analyses suggest that subtle substitutions would affect optimal binding during the catalytic cycle between the NifH from Methanococcales and the bacterial nitrogenase, implying differences in the N2-machinery from these ancient archaea.

Deciphering the Nitrogen Fixation Gene Cluster in Vibrio natriegens: A Study on Optimized Expression and Application of Nitrogenase

Microbial nitrogen fixation presents a viable alternative to chemical fertilizers, yet the limited colonization and specificity of naturally occurring nitrogen-fixing microorganisms present significant challenges to their widespread application. In this study, we identified a nitrogen fixation gene cluster (VNnif) in Vibrio natriegens (VN) and tested its nitrogenase activity through the acetylene reduction assay. We investigated the potential utilization of nitrogenase by incorporating the nitrogenase gene cluster from VN into plant growth-promoting rhizosphere bacteria Pseudomonas protegens CHA0 and enhancing its activity to 48.16 nmol C2H2/mg/h through promoter replacement and cluster rearrangement. The engineered strain CHA0-PVNnif was found to positively impact the growth of Arabidopsis thaliana col-0 and Triticum aestivum L. (wheat). This study expanded the role of plant growth-promoting rhizobacteria (PGPR) and provided a research foundation for enhancing nitrogenase activity.

Enigmatic evolution of microbial nitrogen fixation: insights from Earth's past

The evolution of nitrogen fixation undoubtedly altered nearly all corners of the biosphere, given the essential role of nitrogen in the synthesis of biomass. To date, there is no unified view on what planetary conditions gave rise to nitrogen fixation or how these conditions have sustained it evolutionarily. Intriguingly, the concentrations of metals that nitrogenases require to function have changed throughout Earth's history. In this review, we describe the interconnection of the metal and nitrogen cycles with nitrogenase evolution and the importance of ancient ecology in the formation of the modern nitrogen cycle. We argue that exploration of the nitrogen cycle's deep past will provide insights into humanity's immediate environmental challenges centered on nitrogen availability.

NFixDB (Nitrogen Fixation DataBase)-a comprehensive integrated database for robust 'omics analysis of diazotrophs

Biological nitrogen fixation is a fundamental biogeochemical process that transforms molecular nitrogen into biologically available nitrogen via diazotrophic microbes. Diazotrophs anaerobically fix nitrogen using the nitrogenase enzyme which is arranged in three different gene clusters: (i) molybdenum nitrogenase (nifHDK) is the most abundant, followed by it's alternatives, (ii) vanadium nitrogenase (vnfHDK) and (iii) iron nitrogenase (anfHDK). Multiple databases have been constructed as resources for diazotrophic 'omics analysis; however, an integrated database based on whole genome references does not exist. Here, we present NFixDB (Nitrogen Fixation DataBase), a comprehensive integrated whole genome based database for diazotrophs, which includes all nitrogenases (nifHDK, vnfHDK, anfHDK) and nitrogenase-like enzymes (e.g. nflHD) linked to ribosomal RNA operons (16S-5S-23S). NFixDB was computed using Hidden Markov Models (HMMs) against the entire whole genome based Genome Taxonomy Database (GTDB R214), providing searchable reference HMMs for all nitrogenase and nitrogenase-like genes, complete ribosomal RNA operons, both GTDB and NCBI/RefSeq taxonomy, and an SQL database for querying matches. We compared NFixDB to nifH databases from Buckley, Zehr, Mise and FunGene finding extensive evidence of nifH, in addition to vnfH and nflH. NFixDB contains >4000 verified nifHDK sequences contained on 50 unique phyla of bacteria and archaea. NFixDB provides the first comprehensive nitrogenase database available to researchers unlocking diazotrophic microbial potential.

Emergence of an Orphan Nitrogenase Protein Following Atmospheric Oxygenation

Molecular innovations within key metabolisms can have profound impacts on element cycling and ecological distribution. Yet, much of the molecular foundations of early evolved enzymes and metabolisms are unknown. Here, we bring one such mystery to relief by probing the birth and evolution of the G-subunit protein, an integral component of certain members of the nitrogenase family, the only enzymes capable of biological nitrogen fixation. The G-subunit is a Paleoproterozoic-age orphan protein that appears more than 1 billion years after the origin of nitrogenases. We show that the G-subunit arose with novel nitrogenase metal dependence and the ecological expansion of nitrogen-fixing microbes following the transition in environmental metal availabilities and atmospheric oxygenation that began ∼2.5 billion years ago. We identify molecular features that suggest early G-subunit proteins mediated cofactor or protein interactions required for novel metal dependency, priming ancient nitrogenases and their hosts to exploit these newly diversified geochemical environments. We further examined the degree of functional specialization in G-subunit evolution with extant and ancestral homologs using laboratory reconstruction experiments. Our results indicate that permanent recruitment of the orphan protein depended on the prior establishment of conserved molecular features and showcase how contingent evolutionary novelties might shape ecologically important microbial innovations.

Mapping Geological Events and Nitrogen Fixation Evolution Onto the Timetree of the Evolution of Nitrogen-Fixation Genes

Nitrogen is essential for all organisms, but biological nitrogen fixation (BNF) occurs only in a small fraction of prokaryotes. Previous studies divided nitrogenase-gene-carrying prokaryotes into Groups I to IV and provided evidence that BNF first evolved in bacteria. This study constructed a timetree of the evolution of nitrogen-fixation genes and estimated that archaea evolved BNF much later than bacteria and that nitrogen-fixing cyanobacteria evolved later than 1,900 MYA, considerably younger than the previous estimate of 2,200 MYA. Moreover, Groups III and II/I diverged ∼2,280 MYA, after the Kenorland supercontinent breakup (∼2,500-2,100 MYA) and the Great Oxidation Event (∼2,400-2,100 MYA); Groups III and Vnf/Anf diverged ∼2,086 MYA, after the Yarrabubba impact (∼2,229 MYA); and Groups II and I diverged ∼1,920 MYA, after the Vredefort impact (∼2,023 MYA). In summary, this study provided a timescale of BNF events and discussed the possible effects of geological events on BNF evolution.

Influence of sulfide on diazotrophic growth of the methanogen Methanococcus maripaludis and its implications for the origin of nitrogenase

Methanogens inhabit euxinic (sulfide-rich) or ferruginous (iron-rich) environments that promote the precipitation of transition metals as metal sulfides, such as pyrite, reducing metal or sulfur availability. Such environments have been common throughout Earth's history raising the question as to how anaerobes obtain(ed) these elements for the synthesis of enzyme cofactors. Here, we show a methanogen can synthesize molybdenum nitrogenase metallocofactors from pyrite as the source of iron and sulfur, enabling nitrogen fixation. Pyrite-grown, nitrogen-fixing cells grow faster and require 25-fold less molybdenum than cells grown under euxinic conditions. Growth yields are 3 to 8 times higher in cultures grown under ferruginous relative to euxinic conditions. Physiological, transcriptomic, and geochemical data indicate these observations are due to sulfide-promoted metal limitation, in particular molybdenum. These findings suggest that molybdenum nitrogenase may have originated in a ferruginous environment that titrated sulfide to form pyrite, facilitating the availability of sufficient iron, sulfur, and molybdenum for cofactor biosynthesis.

Mineral-Bound Trace Metals as Cofactors for Anaerobic Biological Nitrogen Fixation

Nitrogenase is the only known biological enzyme capable of reducing N₂ to bioavailable NH₃. Most nitrogenases use Mo as a metallocofactor, while alternative cofactors V and Fe are also viable. Both geological and bioinformatic evidence suggest an ancient origin of Mo-based nitrogenase in the Archean, despite the low concentration of dissolved Mo in the Archean oceans. This apparent paradox would be resolvable if mineral-bound Mo were bioavailable for nitrogen fixation by ancient diazotrophs. In this study, the bioavailability of mineral-bound Mo, V, and Fe was determined by incubating an obligately anaerobic diazotroph Clostridium kluyveri with Mo-, V-, and Fe-bearing minerals (molybdenite, cavansite, and ferrihydrite, respectively) and basalt under diazotrophic conditions. The results showed that C. kluyveri utilized mineral-associated metals to express nitrogenase genes and fix nitrogen, as measured by the reverse transcription quantitative polymerase chain reaction and acetylene reduction assay, respectively. C. kluyveri secreted chelating molecules to extract metals from the minerals. As a result of microbial weathering, mineral surface chemistry significantly changed, likely due to surface coating by microbial exudates for metal extraction. These results provide important support for the ancient origin of Mo-based nitrogenase, with profound implications for coevolution of the biosphere and geosphere.

Overview of physiological, biochemical, and regulatory aspects of nitrogen fixation in Azotobacter vinelandii

Understanding how Nature accomplishes the reduction of inert nitrogen gas to form metabolically tractable ammonia at ambient temperature and pressure has challenged scientists for more than a century. Such an understanding is a key aspect toward accomplishing the transfer of the genetic determinants of biological nitrogen fixation to crop plants as well as for the development of improved synthetic catalysts based on the biological mechanism. Over the past 30 years, the free-living nitrogen-fixing bacterium Azotobacter vinelandii emerged as a preferred model organism for mechanistic, structural, genetic, and physiological studies aimed at understanding biological nitrogen fixation. This review provides a contemporary overview of these studies and places them within the context of their historical development.

Nitrogenase resurrection and the evolution of a singular enzymatic mechanism

The planetary biosphere is powered by a suite of key metabolic innovations that emerged early in the history of life. However, it is unknown whether life has always followed the same set of strategies for performing these critical tasks. Today, microbes access atmospheric sources of bioessential nitrogen through the activities of just one family of enzymes, nitrogenases. Here, we show that the only dinitrogen reduction mechanism known to date is an ancient feature conserved from nitrogenase ancestors. We designed a paleomolecular engineering approach wherein ancestral nitrogenase genes were phylogenetically reconstructed and inserted into the genome of the diazotrophic bacterial model, Azotobacter vinelandii, enabling an integrated assessment of both in vivo functionality and purified nitrogenase biochemistry. Nitrogenase ancestors are active and robust to variable incorporation of one or more ancestral protein subunits. Further, we find that all ancestors exhibit the reversible enzymatic mechanism for dinitrogen reduction, specifically evidenced by hydrogen inhibition, which is also exhibited by extant A. vinelandii nitrogenase isozymes. Our results suggest that life may have been constrained in its sampling of protein sequence space to catalyze one of the most energetically challenging biochemical reactions in nature. The experimental framework established here is essential for probing how nitrogenase functionality has been shaped within a dynamic, cellular context to sustain a globally consequential metabolism.

Engineering Nitrogenases for Synthetic Nitrogen Fixation: From Pathway Engineering to Directed Evolution

Globally, agriculture depends on industrial nitrogen fertilizer to improve crop growth. Fertilizer production consumes fossil fuels and contributes to environmental nitrogen pollution. A potential solution would be to harness nitrogenases-enzymes capable of converting atmospheric nitrogen N2 to NH3 in ambient conditions. It is therefore a major goal of synthetic biology to engineer functional nitrogenases into crop plants, or bacteria that form symbiotic relationships with crops, to support growth and reduce dependence on industrially produced fertilizer. This review paper highlights recent work toward understanding the functional requirements for nitrogenase expression and manipulating nitrogenase gene expression in heterologous hosts to improve activity and oxygen tolerance and potentially to engineer synthetic symbiotic relationships with plants.

A naturalist perspective of microbiology: Examples from methanogenic archaea

Storytelling has been the primary means of knowledge transfer over human history. The effectiveness and reach of stories are improved when the message is appropriate for the target audience. Oftentimes, the stories that are most well received and recounted are those that have a clear purpose and that are told from a variety of perspectives that touch on the varied interests of the target audience. Whether scientists realize or not, they are accustomed to telling stories of their own scientific discoveries through the preparation of manuscripts, presentations, and lectures. Perhaps less frequently, scientists prepare review articles or book chapters that summarize a body of knowledge on a given subject matter, meant to be more holistic recounts of a body of literature. Yet, by necessity, such summaries are often still narrow in their scope and are told from the perspective of a particular discipline. In other words, interdisciplinary reviews or book chapters tend to be the rarity rather than the norm. Here, we advocate for and highlight the benefits of interdisciplinary perspectives on microbiological subjects.

Comparative Transcriptomics Sheds Light on Remodeling of Gene Expression during Diazotrophy in the Thermophilic Methanogen Methanothermococcus thermolithotrophicus

Some marine thermophilic methanogens are able to perform energy-consuming nitrogen fixation despite deriving only little energy from hydrogenotrophic methanogenesis. We studied this process in Methanothermococcus thermolithotrophicus DSM 2095, a methanogenic archaeon of the order Methanococcales that contributes to the nitrogen pool in some marine environments. We successfully grew this archaeon under diazotrophic conditions in both batch and fermenter cultures, reaching the highest cell density reported so far. Diazotrophic growth depended strictly on molybdenum and, in contrast to other diazotrophs, was not inhibited by tungstate or vanadium. This suggests an elaborate control of metal uptake and a specific metal recognition system for the insertion into the nitrogenase cofactor. Differential transcriptomics of M. thermolithotrophicus grown under diazotrophic conditions with ammonium-fed cultures as controls revealed upregulation of the nitrogenase machinery, including chaperones, regulators, and molybdate importers, as well as simultaneous upregulation of an ammonium transporter and a putative pathway for nitrate and nitrite utilization. The organism thus employs multiple synergistic strategies for uptake of nitrogen nutrients during the early exponential growth phase without altering transcription levels for genes involved in methanogenesis. As a counterpart, genes coding for transcription and translation processes were downregulated, highlighting the maintenance of an intricate metabolic balance to deal with energy constraints and nutrient limitations imposed by diazotrophy. This switch in the metabolic balance included unexpected processes, such as upregulation of the CRISPR-Cas system, probably caused by drastic changes in transcription levels of putative mobile and virus-like elements. IMPORTANCE The thermophilic anaerobic archaeon M. thermolithotrophicus is a particularly suitable model organism to study the coupling of methanogenesis to diazotrophy. Likewise, its capability of simultaneously reducing N₂ and CO₂ into NH₃ and CH₄ with H₂ makes it a viable target for biofuel production. We optimized M. thermolithotrophicus cultivation, resulting in considerably higher cell yields and enabling the successful establishment of N₂-fixing bioreactors. Improved understanding of the N₂ fixation process would provide novel insights into metabolic adaptations that allow this energy-limited extremophile to thrive under diazotrophy, for instance, by investigating its physiology and uncharacterized nitrogenase. We demonstrated that diazotrophic growth of M. thermolithotrophicus is exclusively dependent on molybdenum, and complementary transcriptomics corroborated the expression of the molybdenum nitrogenase system. Further analyses of differentially expressed genes during diazotrophy across three cultivation time points revealed insights into the response to nitrogen limitation and the coordination of core metabolic processes.

Gluconacetobacter diazotrophicus Gene Fitness during Diazotrophic Growth

Plant growth-promoting (PGP) bacteria are important to the development of sustainable agricultural systems. PGP microbes that fix atmospheric nitrogen (diazotrophs) could minimize the application of industrially derived fertilizers and function as a biofertilizer. The bacterium Gluconacetobacter diazotrophicus is a nitrogen-fixing PGP microbe originally discovered in association with sugarcane plants, where it functions as an endophyte. It also forms endophyte associations with a range of other agriculturally relevant crop plants. G. diazotrophicus requires microaerobic conditions for diazotrophic growth. We generated a transposon library for G. diazotrophicus and cultured the library under various growth conditions and culture medium compositions to measure fitness defects associated with individual transposon inserts (transposon insertion sequencing [Tn-seq]). Using this library, we probed more than 3,200 genes and ascertained the importance of various genes for diazotrophic growth of this microaerobic endophyte. We also identified a set of essential genes. IMPORTANCE Our results demonstrate a succinct set of genes involved in diazotrophic growth for G. diazotrophicus, with a lower degree of redundancy than what is found in other model diazotrophs. The results will serve as a valuable resource for those interested in biological nitrogen fixation and will establish a baseline data set for plant free growth, which could complement future studies related to the endophyte relationship.

Reductive biomining of pyrite by methanogens

Pyrite (FeS₂) is the most abundant iron sulfide mineral in Earth's crust. Until recently, FeS₂ has been considered a sink for iron (Fe) and sulfur (S) at low temperature in the absence of oxygen or oxidative weathering, making these elements unavailable to biology. However, anaerobic methanogens can transfer electrons extracellularly to reduce FeS₂ via direct contact with the mineral. Reduction of FeS₂ occurs through a multistep process that generates aqueous sulfide (HS^-) and FeS₂-associated pyrrhotite (Fe_1-xS). Subsequent dissolution of Fe_1-xS provides Fe(II)_(aq), but not HS^-, that rapidly complexes with HS^-_(aq) generated from FeS₂ reduction to form soluble iron sulfur clusters [nFeS_(aq)]. Cells assimilate nFeS_(aq) to meet Fe/S nutritional demands by mobilizing and hyperaccumulating Fe and S from FeS₂. As such, reductive dissolution of FeS₂ by methanogens has important implications for element cycling in anoxic habitats, both today and in the geologic past.

Nitrogenase Cofactor Maturase NifB Isolated from Transgenic Rice is Active in FeMo-co Synthesis

The engineering of nitrogen fixation in plants requires assembly of an active prokaryotic nitrogenase complex, which is yet to be achieved. Nitrogenase biogenesis relies on NifB, which catalyzes the formation of the [8Fe-9S-C] metal cluster NifB-co. This is the first committed step in the biosynthesis of the iron-molybdenum cofactor (FeMo-co) found at the nitrogenase active site. The production of NifB in plants is challenging because this protein is often insoluble in eukaryotic cells, and its [Fe-S] clusters are extremely unstable and sensitive to O₂. As a first step to address this challenge, we generated transgenic rice plants expressing NifB from the Archaea Methanocaldococcus infernus and Methanothermobacter thermautotrophicus. The recombinant proteins were targeted to the mitochondria to limit exposure to O₂ and to have access to essential [4Fe-4S] clusters required for NifB-co biosynthesis. M. infernus and M. thermautotrophicus NifB accumulated as soluble proteins in planta, and the purified proteins were functional in the in vitro FeMo-co synthesis assay. We thus report NifB protein expression and purification from an engineered staple crop, representing a first step in the biosynthesis of a functional NifDK complex, as required for independent biological nitrogen fixation in cereals.

Origin and evolution of nitrogen fixation in prokaryotes

The origin of nitrogen fixation is an important issue in evolutionary biology. While nitrogen is required by all living organisms, only a small fraction of bacteria and archaea can fix nitrogen. The prevailing view is that nitrogen fixation first evolved in archaea and was later transferred to bacteria. However, nitrogen-fixing (Nif) bacteria are far larger in number and far more diverse in ecological niches than Nif archaea. We, therefore, propose the bacteria-first hypothesis, which postulates that nitrogen fixation first evolved in bacteria and was later transferred to archaea. As >30,000 prokaryotic genomes have been sequenced, we conduct an in-depth comparison of the two hypotheses. We first identify the six genes involved in nitrogen fixation in all sequenced prokaryotic genomes and then reconstruct phylogenetic trees using the six Nif proteins individually or in combination. In each of these trees, the earliest lineages are bacterial Nif protein sequences and in the oldest clade (group) the archaeal sequences are all nested inside bacterial sequences, suggesting that the Nif proteins first evolved in bacteria. The bacteria-first hypothesis is further supported by the observation that the majority of Nif archaea carry the major bacterial Mo (molybdenum) transporter (ModABC) rather than the archaeal Mo transporter (WtpABC). Moreover, in our phylogeny of all available ModA and WtpA protein sequences, the earliest lineages are bacterial sequences while archaeal sequences are nested inside bacterial sequences. Furthermore, the bacteria-first hypothesis is supported by available isotopic data. In conclusion, our study strongly supports the bacteria-first hypothesis.

Designing Artificial Metalloenzymes by Tuning of the Environment beyond the Primary Coordination Sphere

Metalloenzymes catalyze a variety of reactions using a limited number of natural amino acids and metallocofactors. Therefore, the environment beyond the primary coordination sphere must play an important role in both conferring and tuning their phenomenal catalytic properties, enabling active sites with otherwise similar primary coordination environments to perform a diverse array of biological functions. However, since the interactions beyond the primary coordination sphere are numerous and weak, it has been difficult to pinpoint structural features responsible for the tuning of activities of native enzymes. Designing artificial metalloenzymes (ArMs) offers an excellent basis to elucidate the roles of these interactions and to further develop practical biological catalysts. In this review, we highlight how the secondary coordination spheres of ArMs influence metal binding and catalysis, with particular focus on the use of native protein scaffolds as templates for the design of ArMs by either rational design aided by computational modeling, directed evolution, or a combination of both approaches. In describing successes in designing heme, nonheme Fe, and Cu metalloenzymes, heteronuclear metalloenzymes containing heme, and those ArMs containing other metal centers (including those with non-native metal ions and metallocofactors), we have summarized insights gained on how careful controls of the interactions in the secondary coordination sphere, including hydrophobic and hydrogen bonding interactions, allow the generation and tuning of these respective systems to approach, rival, and, in a few cases, exceed those of native enzymes. We have also provided an outlook on the remaining challenges in the field and future directions that will allow for a deeper understanding of the secondary coordination sphere a deeper understanding of the secondary coordintion sphere to be gained, and in turn to guide the design of a broader and more efficient variety of ArMs.

Temperature and Geographic Location Impact the Distribution and Diversity of Photoautotrophic Gene Variants in Alkaline Yellowstone Hot Springs

Alkaline hot springs in Yellowstone National Park (YNP) provide a framework to study the relationship between photoautotrophs and temperature. Previous work has focused on studying how cyanobacteria (oxygenic phototrophs) vary with temperature, sulfide, and pH, but many questions remain regarding the ecophysiology of anoxygenic photosynthesis due to the taxonomic and metabolic diversity of these taxa. To this end, we examined the distribution of genes involved in phototrophy, carbon fixation, and nitrogen fixation in eight alkaline (pH 7.3-9.4) hot spring sites near the upper temperature limit of photosynthesis (71ºC) in YNP using metagenome sequencing. Based on genes encoding key reaction center proteins, geographic isolation plays a larger role than temperature in selecting for distinct phototrophic Chloroflexi, while genes typically associated with autotrophy in anoxygenic phototrophs, did not have distinct distributions with temperature. Additionally, we recovered Calvin cycle gene variants associated with Chloroflexi, an alternative carbon fixation pathway in anoxygenic photoautotrophs. Lastly, we recovered several abundant nitrogen fixation gene sequences associated with Roseiflexus, providing further evidence that genes involved in nitrogen fixation in Chloroflexi are more common than previously assumed. Together, our results add to the body of work on the distribution and functional potential of phototrophic bacteria in Yellowstone National Park hot springs and support the hypothesis that a combination of abiotic and biotic factors impact the distribution of phototrophic bacteria in hot springs. Future studies of isolates and metagenome assembled genomes (MAGs) from these data and others will further our understanding of the ecology and evolution of hot spring anoxygenic phototrophs. IMPORTANCE Photosynthetic bacteria in hot springs are of great importance to both microbial evolution and ecology. While a large body of work has focused on oxygenic photosynthesis in cyanobacteria in Mushroom and Octopus Springs in Yellowstone National Park, many questions remain regarding the metabolic potential and ecology of hot spring anoxygenic phototrophs. Anoxygenic phototrophs are metabolically and taxonomically diverse, and further investigations into their physiology will lead to a deeper understanding of microbial evolution and ecology of these taxa. Here, we have quantified the distribution of key genes involved in carbon and nitrogen metabolism in both oxygenic and anoxygenic phototrophs. Our results suggest that temperature >68ºC selects for distinct groups of cyanobacteria and that carbon fixation pathways associated with these taxa are likely subject to the same selective pressure. Additionally, our data suggest that phototrophic Chloroflexi genes and carbon fixation genes are largely influenced by local conditions as evidenced by our gene variant analysis. Lastly, we recovered several genes associated with potentially novel phototrophic Chloroflexi. Together, our results add to the body of work on hot springs in Yellowstone National Park and set the stage for future work on metagenome assembled genomes.

Functional Nitrogenase Cofactor Maturase NifB in Mitochondria and Chloroplasts of Nicotiana benthamiana

Engineering plants to synthesize nitrogenase and assimilate atmospheric N2 will reduce crop dependency on industrial N fertilizers. This technology can be achieved by expressing prokaryotic nitrogen fixation gene products for the assembly of a functional nitrogenase in plants. NifB is a critical nitrogenase component since it catalyzes the first committed step in the biosynthesis of all types of nitrogenase active-site cofactors. Here, we used a library of 30 distinct nifB sequences originating from different phyla and ecological niches to restore diazotrophic growth of an Azotobacter vinelandii nifB mutant. Twenty of these variants rescued the nifB mutant phenotype despite their phylogenetic distance to A. vinelandii. Because multiple protein interactions are required in the iron-molybdenum cofactor (FeMo-co) biosynthetic pathway, the maturation of nitrogenase in a heterologous host can be divided in independent modules containing interacting proteins that function together to produce a specific intermediate. Therefore, nifB functional modules composed of a nifB variant, together with the A. vinelandii NifS and NifU proteins (for biosynthesis of NifB [Fe4S4] clusters) and the FdxN ferredoxin (for NifB function), were expressed in Nicotiana benthamiana chloroplasts and mitochondria. Three archaeal NifB proteins accumulated at high levels in soluble fractions of chloroplasts (Methanosarcina acetivorans and Methanocaldococcus infernus) or mitochondria (M. infernus and Methanothermobacter thermautotrophicus). These NifB proteins were shown to accept [Fe4S4] clusters from NifU and were functional in FeMo-co synthesis in vitro. The accumulation of significant levels of soluble and functional NifB proteins in chloroplasts and mitochondria is critical to engineering biological nitrogen fixation in plants. IMPORTANCE Biological nitrogen fixation is the conversion of inert atmospheric dinitrogen gas into nitrogen-reactive ammonia, a reaction catalyzed by the nitrogenase enzyme of diazotrophic bacteria and archaea. Because plants cannot fix their own nitrogen, introducing functional nitrogenase in cereals and other crop plants would reduce our strong dependency on N fertilizers. NifB is required for the biosynthesis of the active site cofactors of all nitrogenases, which arguably makes it the most important protein in global nitrogen fixation. NifB functionality is therefore a requisite to engineer a plant nitrogenase. The expression of nifB genes from a wide range of prokaryotes into the model diazotroph Azotobacter vinelandii shows a surprising level of genetic complementation suggestive of plasticity in the nitrogenase biosynthetic pathway. In addition, we obtained NifB proteins from both mitochondria and chloroplasts of tobacco that are functional in vitro after reconstitution by providing [Fe4S4] clusters from NifU, paving the way to nitrogenase cofactor biosynthesis in plants.

Genomic Features of the Bundle-Forming Heliobacterium Heliophilum fasciatum

Eight species of heliobacteria have had their genomes sequenced. However, only two of these genomes have been analyzed in detail, those from the thermophilic Heliomicrobium (Hmi.) modesticaldum and the alkaliphilic Heliorestis (Hrs.) convoluta. Here we present analyses of the draft genome sequence of a species of heliobacterium that grows optimally at a moderate temperature and neutral pH. The organism, Heliophilum (Hph.) fasciatum, is phylogenetically unique among cultured heliobacteria and was isolated from rice soil, a common habitat for heliobacteria. The Hph. fasciatum genome contains 3.14 Mbp—similar to that of other reported heliobacteria—but has a G+C base ratio that lies between that of Hmi. modesticaldum and Hrs. convoluta. Many of the genomic features of Hmi. modesticaldum and Hrs. convoluta, such as the absence of genes encoding autotrophic pathways, the presence of a superoperonal cluster of photosynthesis-related genes, and genes encoding endospore-specific proteins, are also characteristic of the Hph. fasciatum genome. However, despite the fact that Hph. fasciatum is diazotrophic, classical nif genes encoding the alpha and beta subunits of dinitrogenase (nifDK) present in other heliobacteria could not be identified. Instead, genes encoding several highly divergent NifDK homologs were present, at least one of which likely encodes a functional dinitrogenase and another a methylthio-alkane reductase (MarDK) for sulfur assimilation. A classical NifH (dinitrogenase reductase) homolog was also absent in Hph. fasciatum, but a related protein was identified that likely carries out this function as well as electron delivery to MarDK. The N₂-fixing system of Hph. fasciatum is therefore distinct from that of other heliobacteria and may have unusual properties.

Genomic and Experimental Analysis of the Biostimulant and Antagonistic Properties of Phytopathogens of Bacillus safensis and Bacillus siamensis

The B. safensis RGM 2450 and B. siamensis RGM 2529 strains were isolated from the rhizosphere of plants presenting resilience to abiotic and biotic stress conditions. To understand the implications of bacteria in resilience, a genomic and experimental analysis was carried out on their biostimulant and phytopathogenic antagonist properties. Genome analyses of both strains indicated that they have the potential to synthesize bioactive compounds such as the battery of non-ribosomal peptides, polyketides, extracellular enzymes and phytohormones. These results were consistent with the antagonistic activities of both strains against the phytopathogens Botrytis cinerea, Colletotrichum acutatum, Fusarium oxysporum and Phytophtora cinnamomi. They also showed the capacity to solubilize phosphorus, fix nitrogen and produce indole acetic acid. This was observed in tomato seedlings grown from seeds inoculated with the mixture of strains which presented significantly greater length as well as wet and dry weight in comparison with the treatments individually inoculated with each strain and the control. Accordingly, the combination of B. safensis RGM 2450 and B. siamensis RGM 2529 showed synergistic biostimulant activity. These findings contribute new knowledge of the genomic and metabolomic properties taking part in the symbiotic interactions between these strains and the plants and uphold the combined use of both strains as a biostimulant.

Reconstruction of nitrogenase predecessors suggests origin from maturase-like proteins

The evolution of biological nitrogen fixation, uniquely catalyzed by nitrogenase enzymes, has been one of the most consequential biogeochemical innovations over life’s history. Though understanding the early evolution of nitrogen fixation has been a longstanding goal from molecular, biogeochemical, and planetary perspectives, its origins remain enigmatic. In this study, we reconstructed the evolutionary histories of nitrogenases, as well as homologous maturase proteins that participate in the assembly of the nitrogenase active-site cofactor but are not able to fix nitrogen. We combined phylogenetic and ancestral sequence inference with an analysis of predicted functionally divergent sites between nitrogenases and maturases to infer the nitrogen-fixing capabilities of their shared ancestors. Our results provide phylogenetic constraints to the emergence of nitrogen fixation and are consistent with a model wherein nitrogenases emerged from maturase-like predecessors. Though the precise functional role of such a predecessor protein remains speculative, our results highlight evolutionary contingency as a significant factor shaping the evolution of a biogeochemically essential enzyme.

Pathways of Iron and Sulfur Acquisition, Cofactor Assembly, Destination, and Storage in Diverse Archaeal Methanogens and Alkanotrophs

Archaeal methanogens, methanotrophs, and alkanotrophs have a high demand for iron (Fe) and sulfur (S); however, little is known of how they acquire, traffic, deploy, and store these elements. Here, we examined the distribution of homologs of proteins mediating key steps in Fe/S metabolism in model microorganisms, including iron(II) sensing/uptake (FeoAB), sulfide extraction from cysteine (SufS), and the biosynthesis of iron-sulfur [Fe-S] clusters (SufBCDE), siroheme (Pch2 dehydrogenase), protoheme (AhbABCD), cytochrome c (Cyt c) (CcmCF), and iron storage/detoxification (Bfr, FtrA, and IssA), among 326 publicly available, complete or metagenome-assembled genomes of archaeal methanogens/methanotrophs/alkanotrophs. The results indicate several prevalent but nonuniversal features, including FeoB, SufBC, and the biosynthetic apparatus for the basic tetrapyrrole scaffold, as well as its siroheme (and F₄₃₀) derivatives. However, several early-diverging genomes lacked SufS and pathways to synthesize and deploy heme. Genomes encoding complete versus incomplete heme biosynthetic pathways exhibited equivalent prevalences of [Fe-S] cluster binding proteins, suggesting an expansion of catalytic capabilities rather than substitution of heme for [Fe-S] in the former group. Several strains with heme binding proteins lacked heme biosynthesis capabilities, while other strains with siroheme biosynthesis capability lacked homologs of known siroheme binding proteins, indicating heme auxotrophy and unknown siroheme biochemistry, respectively. While ferritin proteins involved in ferric oxide storage were widespread, those involved in storing Fe as thioferrate were unevenly distributed. Collectively, the results suggest that differences in the mechanisms of Fe and S acquisition, deployment, and storage have accompanied the diversification of methanogens/methanotrophs/alkanotrophs, possibly in response to differential availability of these elements as these organisms evolved.

IMPORTANCE Archaeal methanogens, methanotrophs, and alkanotrophs, argued to be among the most ancient forms of life, have a high demand for iron (Fe) and sulfur (S) for cofactor biosynthesis, among other uses. Here, using comparative bioinformatic approaches applied to 326 genomes, we show that major differences in Fe/S acquisition, trafficking, deployment, and storage exist in this group. Variation in these characters was generally congruent with the phylogenetic placement of these genomes, indicating that variation in Fe/S usage and deployment has contributed to the diversification and ecology of these organisms. However, incongruency was observed among the distribution of cofactor biosynthesis pathways and known protein destinations for those cofactors, suggesting auxotrophy or yet-to-be-discovered pathways for cofactor biosynthesis.

Purple sulfur bacteria fix N2 via molybdenum-nitrogenase in a low molybdenum Proterozoic ocean analogue

Biological N₂ fixation was key to the expansion of life on early Earth. The N₂-fixing microorganisms and the nitrogenase type used in the Proterozoic are unknown, although it has been proposed that the canonical molybdenum-nitrogenase was not used due to low molybdenum availability. We investigate N₂ fixation in Lake Cadagno, an analogue system to the sulfidic Proterozoic continental margins, using a combination of biogeochemical, molecular and single cell techniques. In Lake Cadagno, purple sulfur bacteria (PSB) are responsible for high N₂ fixation rates, to our knowledge providing the first direct evidence for PSB in situ N₂ fixation. Surprisingly, no alternative nitrogenases are detectable, and N₂ fixation is exclusively catalyzed by molybdenum-nitrogenase. Our results show that molybdenum-nitrogenase is functional at low molybdenum conditions in situ and that in contrast to previous beliefs, PSB may have driven N₂ fixation in the Proterozoic ocean.

Functional analysis of multiple nifB genes of Paenibacillus strains in synthesis of Mo-, Fe- and V-nitrogenases

Background

Biological nitrogen fixation is catalyzed by Mo-, V- and Fe-nitrogenases that are encoded by nif, vnf and anf genes, respectively. NifB is the key protein in synthesis of the cofactors of all nitrogenases. Most diazotrophic Paenibacillus strains have only one nifB gene located in a compact nif gene cluster (nifBHDKENX(orf1)hesAnifV). But some Paenibacillus strains have multiple nifB genes and their functions are not known.

Results

A total of 138 nifB genes are found in the 116 diazotrophic Paenibacillus strains. Phylogeny analysis shows that these nifB genes fall into 4 classes: nifBI class including the genes (named as nifB1 genes) that are the first gene within the compact nif gene cluster, nifBII class including the genes (named as nifB2 genes) that are adjacent to anf or vnf genes, nifBIII class whose members are designated as nifB3 genes and nifBIV class whose members are named as nifB4 genes are scattered on genomes. Functional analysis by complementation of the ∆nifB mutant of P. polymyxa which has only one nifB gene has shown that both nifB1 and nifB2 are active in synthesis of Mo-nitrogenase, while nifB3 and nifB4 genes are not. Deletion analysis also has revealed that nifB1 of Paenibacillus sabinae T27 is involved in synthesis of Mo-nitrogenase, while nifB3 and nifB4 genes are not. Complementation of the P. polymyxa ∆nifBHDK mutant with the four reconstituted operons: nifB1anfHDGK, nifB2anfHDGK, nifB1vnfHDGK and nifB2vnfHDGK, has shown both that nifB1 and nifB2 were able to support synthesis of Fe- or V-nitrogenases. Transcriptional results obtained in the original Paenibacillus strains are consistent with the complementation results.

Conclusions

The multiple nifB genes of the diazotrophic Paenibacillus strains are divided into 4 classes. The nifB1 located in a compact nif gene cluster (nifBHDKENX(orf1)hesAnifV) and the nifB2 genes being adjacent to nif or anf or vnf genes are active in synthesis of Mo-, Fe and V-nitrogenases, but nifB3 and nifB4 are not. The reconstituted anf system comprising 8 genes (nifBanfHDGK and nifXhesAnifV) and vnf system comprising 10 genes (nifBvnfHDGKEN and nifXhesAnifV) support synthesis of Fe-nitrogenase and V-nitrogenase in Paenibacillus background, respectively.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-021-01629-9.

A New Analysis of Archaea-Bacteria Domain Separation: Variable Phylogenetic Distance and the Tempo of Early Evolution

Comparative genomics and molecular phylogenetics are foundational for understanding biological evolution. Although many studies have been made with the aim of understanding the genomic contents of early life, uncertainty remains. A study by Weiss et al. (Weiss MC, Sousa FL, Mrnjavac N, Neukirchen S, Roettger M, Nelson-Sathi S, Martin WF. 2016. The physiology and habitat of the last universal common ancestor. Nat Microbiol. 1(9):16116.) identified a number of protein families in the last universal common ancestor of archaea and bacteria (LUCA) which were not found in previous works. Here, we report new research that suggests the clustering approaches used in this previous study undersampled protein families, resulting in incomplete phylogenetic trees which do not reflect protein family evolution. Phylogenetic analysis of protein families which include more sequence homologs rejects a simple LUCA hypothesis based on phylogenetic separation of the bacterial and archaeal domains for a majority of the previously identified LUCA proteins (∼82%). To supplement limitations of phylogenetic inference derived from incompletely populated orthologous groups and to test the hypothesis of a period of rapid evolution preceding the separation of the domains, we compared phylogenetic distances both within and between domains, for thousands of orthologous groups. We find a substantial diversity of interdomain versus intradomain branch lengths, even among protein families which exhibit a single domain separating branch and are thought to be associated with the LUCA. Additionally, phylogenetic trees with long interdomain branches relative to intradomain branches are enriched in information categories of protein families in comparison to those associated with metabolic functions. These results provide a new view of protein family evolution and temper claims about the phenotype and habitat of the LUCA.

PPIT: an R package for inferring microbial taxonomy from nifH sequences

MOTIVATION

Linking microbial community members to their ecological functions is a central goal of environmental microbiology. When assigned taxonomy, amplicon sequences of metabolic marker genes can suggest such links, thereby offering an overview of the phylogenetic structure underpinning particular ecosystem functions. However, inferring microbial taxonomy from metabolic marker gene sequences remains a challenge, particularly for the frequently sequenced nitrogen fixation marker gene, nitrogenase reductase (nifH). Horizontal gene transfer in recent nifH evolutionary history can confound taxonomic inferences drawn from the pairwise identity methods used in existing software. Other methods for inferring taxonomy are not standardized and require manual inspection that is difficult to scale.

RESULTS

We present Phylogenetic Placement for Inferring Taxonomy (PPIT), an R package that infers microbial taxonomy from nifH amplicons using both phylogenetic and sequence identity approaches. After users place query sequences on a reference nifH gene tree provided by PPIT (n = 6317 full-length nifH sequences), PPIT searches the phylogenetic neighborhood of each query sequence and attempts to infer microbial taxonomy. An inference is drawn only if references in the phylogenetic neighborhood are: (1) taxonomically consistent and (2) share sufficient pairwise identity with the query, thereby avoiding erroneous inferences due to known horizontal gene transfer events. We find that PPIT returns a higher proportion of correct taxonomic inferences than BLAST-based approaches at the cost of fewer total inferences. We demonstrate PPIT on deep-sea sediment and find that Deltaproteobacteria are the most abundant potential diazotrophs. Using this dataset we show that emending PPIT inferences based on visual inspection of query sequence placement can achieve taxonomic inferences for nearly all sequences in a query set. We additionally discuss how users can apply PPIT to the analysis of other marker genes.

AVAILABILITY

PPIT is freely available to non-commercial users at https://github.com/bkapili/ppit. Installation includes a vignette that demonstrates package use and reproduces the nifH amplicon analysis discussed here. The raw nifH amplicon sequence data have been deposited in the GenBank, EMBL, and DDBJ databases under BioProject number PRJEB37167.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Expression, Isolation, and Characterization of Vanadium Nitrogenase from Azotobacter vinelandii

Nitrogenases are the sole enzymes known to mediate biological nitrogen fixation, an essential process for sustaining life on earth. Among the three known variants, molybdenum nitrogenase is the best-studied to date. Recent work on the alternative vanadium nitrogenase provided important insights into the mechanism of nitrogen fixation since this enzyme differs from its molybdenum counterpart in some important aspects. Here, we present a protocol to obtain unmodified vanadium nitrogenase in high yield and purity from the paradigmatic diazotroph Azotobacter vinelandii, including procedures for cell cultivation, purification, and protein characterization.

Radiation of nitrogen-metabolizing enzymes across the tree of life tracks environmental transitions in Earth history

Nitrogen is an essential element to life and exerts a strong control on global biological productivity. The rise and spread of nitrogen-utilizing microbial metabolisms profoundly shaped the biosphere on the early Earth. Here, we reconciled gene and species trees to identify birth and horizontal gene transfer events for key nitrogen-cycling genes, dated with a time-calibrated tree of life, in order to examine the timing of the proliferation of these metabolisms across the tree of life. Our results provide new insights into the evolution of the early nitrogen cycle that expand on geochemical reconstructions. We observed widespread horizontal gene transfer of molybdenum-based nitrogenase back to the Archean, minor horizontal transfer of genes for nitrate reduction in the Archean, and an increase in the proliferation of genes metabolizing nitrite around the time of the Mesoproterozoic (~1.5 Ga). The latter coincides with recent geochemical evidence for a mid-Proterozoic rise in oxygen levels. Geochemical evidence of biological nitrate utilization in the Archean and early Proterozoic may reflect at least some contribution of dissimilatory nitrate reduction to ammonium (DNRA) rather than pure denitrification to N₂ . Our results thus help unravel the relative dominance of two metabolic pathways that are not distinguishable with current geochemical tools. Overall, our findings thus provide novel constraints for understanding the evolution of the nitrogen cycle over time and provide insights into the bioavailability of various nitrogen sources in the early Earth with possible implications for the emergence of eukaryotic life.

Nitrogen-fixing Ability and Nitrogen Fixation-related Genes of Thermophilic Fermentative Bacteria in the Genus Caldicellulosiruptor

Fermentative nitrogen-fixing bacteria have not yet been examined in detail in thermal environments. In the present study, we isolated the thermophilic fermentative bacterium, strain YA01 from a hot spring. This strain grew at temperatures up to 78°C. A phylogenetic analysis based on its 16S rRNA gene sequence indicated that strain YA01 belonged to the genus Caldicellulosiruptor, which are fermentative bacteria in the phylum Firmicutes, with 97.7–98.0% sequence identity to its closest relatives. Strain YA01 clearly exhibited N₂-dependent growth at 70°C. We also confirmed N₂-dependent growth in the relatives of strain YA01, Caldicellulosiruptor hydrothermalis 108 and Caldicellulosiruptor kronotskyensis 2002. The nitrogenase activities of these three strains were examined using the acetylene reduction assay. Similar activities were detected for all tested strains, and were slightly suppressed by the addition of ammonium. A genome analysis revealed that strain YA01, as well as other Caldicellulosiruptor, possessed a gene set for nitrogen fixation, but lacked the nifN gene, which encodes a nitrogenase iron-molybdenum cofactor biosynthesis protein that is commonly detected in nitrogen-fixing bacteria. The amino acid sequences of nitrogenase encoded by nifH, nifD, and nifK shared 92–98% similarity in Caldicellulosiruptor. A phylogenetic tree of concatenated NifHDK sequences showed that NifHDK of Caldicellulosiruptor was in the deepest clade. To the best of our knowledge, this is the first study to demonstrate the nitrogen-fixing ability of fermentative bacteria at 70°C. Caldicellulosiruptor may have retained an ancient nitrogen-fixing enzyme system.

Evolutionary History of Bioessential Elements Can Guide the Search for Life in the Universe

Our understanding of life in the universe comes from one sample, life on Earth. Current and next-generation space missions will target exoplanets as well as planets and moons in our own solar system with the primary goal of detecting, interpreting and characterizing indications of possible biological activity. Thus, understanding life's fundamental characteristics is increasingly critical for detecting and interpreting potential biological signatures elsewhere in the universe. Astrobiologists have outlined the essential roles of carbon and water for life, but we have yet to decipher the rules governing the evolution of how living organisms use bioessential elements. Does the suite of life's essential chemical elements on Earth constitute only one possible evolutionary outcome? Are some elements so essential for biological functions that evolution will select for them despite low availability? How would this play out on other worlds that have different relative element abundances? When we look for life in the universe, or the conditions that could give rise to life, we must learn how to recognize it in extremely different chemical and environmental conditions from those on Earth. We argue that by exposing self-organizing biotic chemistries to different combinations of abiotic materials, and by mapping the evolutionary history of metalloenzyme biochemistry onto geological availabilities of metals, alternative element choices that are very different from life's present-day molecular structure might result. A greater understanding of the paleomolecular evolutionary history of life on Earth will create a predictive capacity for detecting and assessing life's existence on worlds where alternate evolutionary paths might have been taken.

Reactivity, Mechanism, and Assembly of the Alternative Nitrogenases

Biological nitrogen fixation is catalyzed by the enzyme nitrogenase, which facilitates the cleavage of the relatively inert triple bond of N2. Nitrogenase is most commonly associated with the molybdenum-iron cofactor called FeMoco or the M-cluster, and it has been the subject of extensive structural and spectroscopic characterization over the past 60 years. In the late 1980s and early 1990s, two "alternative nitrogenase" systems were discovered, isolated, and found to incorporate V or Fe in place of Mo. These systems are regulated by separate gene clusters; however, there is a high degree of structural and functional similarity between each nitrogenase. Limited studies with the V- and Fe-nitrogenases initially demonstrated that these enzymes were analogously active as the Mo-nitrogenase, but more recent investigations have found capabilities that are unique to the alternative systems. In this review, we will discuss the reactivity, biosynthetic, and mechanistic proposals for the alternative nitrogenases as well as their electronic and structural properties in comparison to the well-characterized Mo-dependent system. Studies over the past 10 years have been particularly fruitful, though key aspects about V- and Fe-nitrogenases remain unexplored.

The Spectroscopy of Nitrogenases

Nitrogenases are responsible for biological nitrogen fixation, a crucial step in the biogeochemical nitrogen cycle. These enzymes utilize a two-component protein system and a series of iron-sulfur clusters to perform this reaction, culminating at the FeMco active site (M = Mo, V, Fe), which is capable of binding and reducing N2 to 2NH3. In this review, we summarize how different spectroscopic approaches have shed light on various aspects of these enzymes, including their structure, mechanism, alternative reactivity, and maturation. Synthetic model chemistry and theory have also played significant roles in developing our present understanding of these systems and are discussed in the context of their contributions to interpreting the nature of nitrogenases. Despite years of significant progress, there is still much to be learned from these enzymes through spectroscopic means, and we highlight where further spectroscopic investigations are needed.

Transfer of Nitrogen Fixation (nif) Genes to Non-diazotrophic Hosts

Nitrogen is one of the most important nutrients for plant growth. To enhance crop productivity, chemical nitrogen fertilizer is commonly applied in agriculture. Biological nitrogen fixation, the conversion of atmospheric N2 to NH3 , is an important source of nitrogen input in agriculture and represents a promising substitute for chemical nitrogen fertilizers. However, nitrogen fixation is only sporadically distributed within bacteria and archaea (diazotrophs). Thus, many biologists hope to reconstitute a nitrogenase biosynthetic pathway in a eukaryotic host, with the final aim of developing N2 -fixing cereal crops. With the advent of synthetic biology and a deep understanding of the fundamental genetic determinants necessary to sustain nitrogen fixation in bacteria, much progress has been made toward this goal. Transfer of native and refactored nif (nitrogen fixation) genes to non-diazotrophs has been attempted in model bacteria, yeast, and plants. Specifically, nif genes from Klebsiella oxytoca, Azotobacter vinelandii, and Paenibacillus polymyxa have been successfully transferred and expressed in Escherichia coli, Saccharomyces cerevisiae, and even in the tobacco plant. These advances have laid the groundwork to enable cereal crops to "fix" nitrogen themselves to sustain their growth and yield.

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.

Reconstructing the evolutionary history of nitrogenases: Evidence for ancestral molybdenum-cofactor utilization

The nitrogenase metalloenzyme family, essential for supplying fixed nitrogen to the biosphere, is one of life's key biogeochemical innovations. The three forms of nitrogenase differ in their metal dependence, each binding either a FeMo-, FeV-, or FeFe-cofactor where the reduction of dinitrogen takes place. The history of nitrogenase metal dependence has been of particular interest due to the possible implication that ancient marine metal availabilities have significantly constrained nitrogenase evolution over geologic time. Here, we reconstructed the evolutionary history of nitrogenases, and combined phylogenetic reconstruction, ancestral sequence inference, and structural homology modeling to evaluate the potential metal dependence of ancient nitrogenases. We find that active-site sequence features can reliably distinguish extant Mo-nitrogenases from V- and Fe-nitrogenases and that inferred ancestral sequences at the deepest nodes of the phylogeny suggest these ancient proteins most resemble modern Mo-nitrogenases. Taxa representing early-branching nitrogenase lineages lack one or more biosynthetic nifE and nifN genes that both contribute to the assembly of the FeMo-cofactor in studied organisms, suggesting that early Mo-nitrogenases may have utilized an alternate and/or simplified pathway for cofactor biosynthesis. Our results underscore the profound impacts that protein-level innovations likely had on shaping global biogeochemical cycles throughout the Precambrian, in contrast to organism-level innovations that characterize the Phanerozoic Eon.

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.

Spectroscopic Description of the E1 State of Mo Nitrogenase Based on Mo and Fe X-ray Absorption and Mössbauer Studies

Mo nitrogenase (N2ase) utilizes a two-component protein system, the catalytic MoFe and its electron-transfer partner FeP, to reduce atmospheric dinitrogen (N2) to ammonia (NH3). The FeMo cofactor contained in the MoFe protein serves as the catalytic center for this reaction and has long inspired model chemistry oriented toward activating N2. This field of chemistry has relied heavily on the detailed characterization of how Mo N2ase accomplishes this feat. Understanding the reaction mechanism of Mo N2ase itself has presented one of the most challenging problems in bioinorganic chemistry because of the ephemeral nature of its catalytic intermediates, which are difficult, if not impossible, to singly isolate. This is further exacerbated by the near necessity of FeP to reduce native MoFe, rendering most traditional means of selective reduction inept. We have now investigated the first fundamental intermediate of the MoFe catalytic cycle, E1, as prepared both by low-flux turnover and radiolytic cryoreduction, using a combination of Mo Kα high-energy-resolution fluorescence detection and Fe K-edge partial-fluorescence-yield X-ray absorption spectroscopy techniques. The results demonstrate that the formation of this state is the result of an Fe-centered reduction and that Mo remains redox-innocent. Furthermore, using Fe X-ray absorption and 57Fe Mössbauer spectroscopies, we correlate a previously reported unique species formed under cryoreducing conditions to the natively formed E1 state through annealing, demonstrating the viability of cryoreduction in studying the catalytic intermediates of MoFe.

Geobiological feedbacks, oxygen, and the evolution of nitrogenase

Biological nitrogen fixation via the activity of nitrogenase is one of the most important biological innovations, allowing for an increase in global productivity that eventually permitted the emergence of higher forms of life. The complex metalloenzyme termed nitrogenase contains complex iron-sulfur cofactors. Three versions of nitrogenase exist that differ mainly by the presence or absence of a heterometal at the active site metal cluster (either Mo or V). Mo-dependent nitrogenase is the most common while V-dependent or heterometal independent (Fe-only) versions are often termed alternative nitrogenases since they have apparent lower activities for N₂ reduction and are expressed in the absence of Mo. Phylogenetic data indicates that biological nitrogen fixation emerged in an anaerobic, thermophilic ancestor of hydrogenotrophic methanogens and later diversified via lateral gene transfer into anaerobic bacteria, and eventually aerobic bacteria including Cyanobacteria. Isotopic evidence suggests that nitrogenase activity existed at 3.2 Ga, prior to the advent of oxygenic photosynthesis and rise of oxygen in the atmosphere, implying the presence of favorable environmental conditions for oxygen-sensitive nitrogenase to evolve. Following the proliferation of oxygenic phototrophs, diazotrophic organisms had to develop strategies to protect nitrogenase from oxygen inactivation and generate the right balance of low potential reducing equivalents and cellular energy for growth and nitrogen fixation activity. Here we review the fundamental advances in our understanding of biological nitrogen fixation in the context of the emergence, evolution, and taxonomic distribution of nitrogenase, with an emphasis placed on key events associated with its emergence and diversification from anoxic to oxic environments.

How to resurrect ancestral proteins as proxies for ancient biogeochemistry

Throughout the history of life, enzymes have served as the primary molecular mediators of biogeochemical cycles by catalyzing the metabolic pathways that interact with geochemical substrates. The byproducts of enzymatic activities have been preserved as chemical and isotopic signatures in the geologic record. However, interpretations of these signatures are limited by the assumption that such enzymes have remained functionally conserved over billions of years of molecular evolution. By reconstructing ancient genetic sequences in conjunction with laboratory enzyme resurrection, preserved biogeochemical signatures can instead be related to experimentally constrained, ancestral enzymatic properties. We may thereby investigate instances within molecular evolutionary trajectories potentially tied to significant biogeochemical transitions evidenced in the geologic record. Here, we survey recent enzyme resurrection studies to provide a reasoned assessment of areas of success and common pitfalls relevant to ancient biogeochemical applications. We conclude by considering the Great Oxidation Event, which provides a constructive example of a significant biogeochemical transition that warrants investigation with ancestral enzyme resurrection. This event also serves to highlight the pitfalls of facile interpretation of paleophenotype models and data, as applied to two examples of enzymes that likely both influenced and were influenced by the rise of atmospheric oxygen - RuBisCO and nitrogenase.

Combined Assembly and Targeted Integration of Multigene for Nitrogenase Biosynthetic Pathway in Saccharomyces cerevisiae

Biological nitrogen fixation, a process unique to diazotrophic prokaryote, is catalyzed by the nitrogenase complex. There has been a long-standing interest in reconstituting a nitrogenase biosynthetic pathway in a eukaryotic host with the final aim of developing N₂-fixing cereal crops. In this study, we report that a nitrogenase biosynthetic pathway (∼38 kb containing 15 genes) was assembled in two individual one-step methods via in vivo assembly and integrated at δ and HO sites in Saccharomyces cerevisiae chromosome. Of the 15 genes, 11 genes (nifB, nifH, nifD, nifK, nifE, nifN, nifX, hesA, nifV, groES, groEL) were from Paenibacillus polymyxa WLY78 and 4 genes (nifS, nifU, nifF, nifJ) were from Klebsiella oxytoca. The 15-gene nitrogenase biosynthetic pathway was correctly assembled and transcribed in the recombinant S. cerevisiae. The NifDK tetramer with an identical molecular weight as that of P. polymyxa was formed in yeast and the expressed NifH exhibited the activity of Fe protein. This study demonstrates that it will be possible to produce active nitrogenase in eukaryotic hosts.

Early Archean origin of Photosystem II

Photosystem II is a photochemical reaction center that catalyzes the light-driven oxidation of water to molecular oxygen. Water oxidation is the distinctive photochemical reaction that permitted the evolution of oxygenic photosynthesis and the eventual rise of eukaryotes. At what point during the history of life an ancestral photosystem evolved the capacity to oxidize water still remains unknown. Here, we study the evolution of the core reaction center proteins of Photosystem II using sequence and structural comparisons in combination with Bayesian relaxed molecular clocks. Our results indicate that a homodimeric photosystem with sufficient oxidizing power to split water had already appeared in the early Archean about a billion years before the most recent common ancestor of all described Cyanobacteria capable of oxygenic photosynthesis, and well before the diversification of some of the known groups of anoxygenic photosynthetic bacteria. Based on a structural and functional rationale, we hypothesize that this early Archean photosystem was capable of water oxidation to oxygen and had already evolved protection mechanisms against the formation of reactive oxygen species. This would place primordial forms of oxygenic photosynthesis at a very early stage in the evolutionary history of life.

Phylogenetic Diversity of Nitrogenase Reductase Genes and Possible Nitrogen-Fixing Bacteria in Thermophilic Chemosynthetic Microbial Communities in Nakabusa Hot Springs

Chemosynthetic microbial communities develop and form dense cell aggregates in slightly alkaline sulfidic hot springs in the temperature range of 70–86°C at Nakabusa, Japan. Nitrogenase activity has recently been detected in the microbial communities collected. To identify possible members capable of nitrogen fixation, we examined the diversities of 16S rRNA and nitrogenase reductase (NifH) gene sequences in four types of chemosynthetic communities with visually different colors and thicknesses. The results of a 16S rRNA gene analysis indicated that all four microbial communities had similar bacterial constituents; the phylum Aquificae was the dominant member, followed in abundance by Thermodesulfobacteria, Firmicutes, and Thermotogae. Most of the NifH sequences were related to sequences reported in hydrothermal vents and terrestrial hot springs. The results of a phylogenetic analysis of NifH sequences revealed diversity in this gene among the communities collected, distributed within 7 phylogenetic groups. NifH sequences affiliated with Aquificae (Hydrogenobacter/Thermocrinis) and Firmicutes (Caldicellulosiruptor) were abundant. At least two different energy metabolic pathways appeared to be related to nitrogen fixation in the communities analyzed; aerobic sulfur/hydrogen-oxidizing bacteria in Aquificae and fermentative bacteria in Firmicutes. The metabolic characteristics of these two dominant phyla differed from those previously inferred from nitrogenase activity assays on chemosynthetic communities, which were associated with hydrogen-dependent autotrophic sulfate reduction. These assays may correspond to the observed NifH sequences that are distantly related to the known species of Thermodesulfovibrio sp. (Nitrospirae) detected in the present study. The activities of nitrogen-fixing organisms in communities may depend on redox states as well as the availability of electron donors, acceptors, and carbon sources.

Nitrous oxide from chemodenitrification: A possible missing link in the Proterozoic greenhouse and the evolution of aerobic respiration

The potent greenhouse gas nitrous oxide (N₂ O) may have been an important constituent of Earth's atmosphere during Proterozoic (~2.5-0.5 Ga). Here, we tested the hypothesis that chemodenitrification, the rapid reduction of nitric oxide by ferrous iron, would have enhanced the flux of N₂ O from ferruginous Proterozoic seas. We empirically derived a rate law, d N 2 O d t = 7.2 × 10 - 5 [ Fe 2 + ] 0.3 [ NO ] 1 , and measured an isotopic site preference of +16‰ for the reaction. Using this empirical rate law, and integrating across an oceanwide oxycline, we found that low nM NO and μM-low mM Fe²⁺ concentrations could have sustained a sea-air flux of 100-200 Tg N₂ O-N year^-1 , if N₂ fixation rates were near-modern and all fixed N₂ was emitted as N₂ O. A 1D photochemical model was used to obtain steady-state atmospheric N₂ O concentrations as a function of sea-air N₂ O flux across the wide range of possible pO₂ values (0.001-1 PAL). At 100-200 Tg N₂ O-N year^-1 and >0.1 PAL O₂ , this model yielded low-ppmv N₂ O, which would produce several degrees of greenhouse warming at 1.6 ppmv CH₄ and 320 ppmv CO₂ . These results suggest that enhanced N₂ O production in ferruginous seawater via a previously unconsidered chemodenitrification pathway may have helped to fill a Proterozoic "greenhouse gap," reconciling an ice-free Mesoproterozoic Earth with a less luminous early Sun. A particularly notable result was that high N₂ O fluxes at intermediate O₂ concentrations (0.01-0.1 PAL) would have enhanced ozone screening of solar UV radiation. Due to rapid photolysis in the absence of an ozone shield, N₂ O is unlikely to have been an important greenhouse gas if Mesoproterozoic O₂ was 0.001 PAL. At low O₂ , N₂ O might have played a more important role as life's primary terminal electron acceptor during the transition from an anoxic to oxic surface Earth, and correspondingly, from anaerobic to aerobic metabolisms.

Study on the bacterial and archaeal community structure and diversity of activated sludge from three wastewater treatment plants

In this study, the bacterial and archaeal communities along with their functions of activated sludge from three wastewater treatment plants were investigated by Illumina MiSeq Platform. The treatment processes were modified A/A/O, DE oxidation ditch and pre-anaerobic carrousel oxidation ditch, respectively. The taxonomic analyses showed that Proteobacteria was the predominant bacterial phylum, and Nitrosospira was the dominant nitrification genus. Candidatus Accumulibacter was abundant in DE oxidation ditch process, and the main archaea communities were methanosaeta-like species which had the capability to anaerobic ammonia oxidation. The results illustrated that anaerobic ammonium oxidation played an important role in the nitrogen metabolism and there might be other unknown phosphate-accumulating organisms (PAOs) performing phosphorus removal in activated sludge. The predicted function analyses indicated that both bacteria and archaea were involved in nitrification, denitrification, ammonification and phosphorus removal processes, and their relative abundance varied metabolic modules differed from each other.

Origin and Evolution of Flavin-Based Electron Bifurcating Enzymes

Twelve evolutionarily unrelated oxidoreductases form enzyme complexes that catalyze the simultaneous coupling of exergonic and endergonic oxidation–reduction reactions to circumvent thermodynamic barriers and minimize free energy loss in a process known as flavin-based electron bifurcation. Common to these 12 bifurcating (Bf) enzymes are protein-bound flavin, the proposed site of bifurcation, and the electron carrier ferredoxin. Despite the documented role of Bf enzymes in balancing the redox state of intracellular electron carriers and in improving the efficiency of cellular metabolism, a comprehensive description of the diversity and evolutionary history of Bf enzymes is lacking. Here, we report the taxonomic distribution, functional diversity, and evolutionary history of Bf enzyme homologs in 4,588 archaeal, bacterial, and eukaryal genomes and 3,136 community metagenomes. Bf homologs were primarily detected in the genomes of anaerobes, including those of sulfate-reducers, acetogens, fermenters, and methanogens. Phylogenetic analyses of Bf enzyme catalytic subunits (oxidoreductases) suggest they were not a property of the Last Universal Common Ancestor of Archaea and Bacteria, which is consistent with the limited and unique taxonomic distributions of enzyme homologs among genomes. Further, phylogenetic analyses of oxidoreductase subunits reveal that non-Bf homologs predate Bf homologs. These observations indicate that multiple independent recruitments of flavoproteins to existing oxidoreductases enabled coupling of numerous new electron Bf reactions. Consistent with the role of these enzymes in the energy metabolism of anaerobes, homologs of Bf enzymes were enriched in metagenomes from subsurface environments relative to those from surface environments. Phylogenetic analyses of homologs from metagenomes reveal that the earliest evolving homologs of most Bf enzymes are from subsurface environments, including fluids from subsurface rock fractures and hydrothermal systems. Collectively, these data suggest strong selective pressures drove the emergence of Bf enzyme complexes via recruitment of flavoproteins that allowed for an increase in the efficiency of cellular metabolism and improvement in energy capture in anaerobes inhabiting a variety of subsurface anoxic habitats where the energy yield of oxidation-reduction reactions is generally low.