Phylogeny of Nitrogenase Structural and Assembly Components Reveals New Insights into the Origin and Distribution of Nitrogen Fixation across Bacteria and Archaea

The phylogeny of nitrogenase has only been analyzed using the structural proteins NifHDK. As nifHDKENB has been established as the minimum number of genes necessary for in silico prediction of diazotrophy, we present an updated phylogeny of diazotrophs using both structural (NifHDK) and cofactor assembly proteins (NifENB). Annotated Nif sequences were obtained from InterPro from 963 culture-derived genomes. Nif sequences were aligned individually and concatenated to form one NifHDKENB sequence. Phylogenies obtained using PhyML, FastTree, RapidNJ, and ASTRAL from individuals and concatenated protein sequences were compared and analyzed. All six genes were found across the Actinobacteria, Aquificae, Bacteroidetes, Chlorobi, Chloroflexi, Cyanobacteria, Deferribacteres, Firmicutes, Fusobacteria, Nitrospira, Proteobacteria, PVC group, and Spirochaetes, as well as the Euryarchaeota. The phylogenies of individual Nif proteins were very similar to the overall NifHDKENB phylogeny, indicating the assembly proteins have evolved together. Our higher resolution database upheld the three cluster phylogeny, but revealed undocumented horizontal gene transfers across phyla. Only 48% of the 325 genera containing all six nif genes are currently supported by biochemical evidence of diazotrophy. In addition, this work provides reference for any inter-phyla comparison of Nif sequences and a quality database of Nif proteins that can be used for identifying new Nif sequences.

The Evolution of Molybdenum Dependent Nitrogenase in Cyanobacteria

Nitrogen fixation plays a crucial role in the nitrogen cycle by helping to convert nitrogen into a form usable by other organisms. Bacteria capable of fixing nitrogen are found in six phyla including Cyanobacteria. Molybdenum dependent nitrogenase (nif) genes are thought to share a single origin as they have homologs in various phyla. However, diazotrophic bacteria have a mosaic distribution within the cyanobacterial lineage. Therefore, the aim of this study was to determine the cause of this mosaic distribution. We identified nif gene operon structures in the genomes of 85 of the 179 cyanobacterial strains for which whole genome sequences were available. Four nif operons were conserved in each diazotroph Cyanobacterium, although there were some gene translocations and insertions. Phylogenetic inference of these genes did not reveal horizontal gene transfer from outside the phylum Cyanobacteria. These results support the hypothesis that the mosaic distribution of diazotrophic bacteria in the cyanobacterial lineage is the result of the independent loss of nif genes inherited from common cyanobacterial ancestors in each lineage.

Comparative Genomics Provides Insights into the Taxonomy of Azoarcus and Reveals Separate Origins of Nif Genes in the Proposed Azoarcus and Aromatoleum Genera

Among other attributes, the Betaproteobacterial genus Azoarcus has biotechnological importance for plant growth-promotion and remediation of petroleum waste-polluted water and soils. It comprises at least two phylogenetically distinct groups. The “plant-associated” group includes strains that are isolated from the rhizosphere or root interior of the C4 plant Kallar Grass, but also strains from soil and/or water; all are considered to be obligate aerobes and all are diazotrophic. The other group (now partly incorporated into the new genus Aromatoleum) comprises a diverse range of species and strains that live in water or soil that is contaminated with petroleum and/or aromatic compounds; all are facultative or obligate anaerobes. Some are diazotrophs. A comparative genome analysis of 32 genomes from 30 Azoarcus-Aromatoleum strains was performed in order to delineate generic boundaries more precisely than the single gene, 16S rRNA, that has been commonly used in bacterial taxonomy. The origin of diazotrophy in Azoarcus-Aromatoleum was also investigated by comparing full-length sequences of nif genes, and by physiological measurements of nitrogenase activity using the acetylene reduction assay. Based on average nucleotide identity (ANI) and whole genome analyses, three major groups could be discerned: (i) Azoarcus comprising Az. communis, Az. indigens and Az. olearius, and two unnamed species complexes, (ii) Aromatoleum Group 1 comprising Ar. anaerobium, Ar. aromaticum, Ar. bremense, and Ar. buckelii, and (iii) Aromatoleum Group 2 comprising Ar. diolicum, Ar. evansii, Ar. petrolei, Ar. toluclasticum, Ar. tolulyticum, Ar. toluolicum, and Ar. toluvorans. Single strain lineages such as Azoarcus sp. KH32C, Az. pumilus, and Az. taiwanensis were also revealed. Full length sequences of nif-cluster genes revealed two groups of diazotrophs in Azoarcus-Aromatoleum with nif being derived from Dechloromonas in Azoarcus sensu stricto (and two Thauera strains) and from Azospira in Aromatoleum Group 2. Diazotrophy was confirmed in several strains, and for the first time in Az. communis LMG5514, Azoarcus sp. TTM-91 and Ar. toluolicum TT. In terms of ecology, with the exception of a few plant-associated strains in Azoarcus (s.s.), across the group, most strains/species are found in soil and water (often contaminated with petroleum or related aromatic compounds), sewage sludge, and seawater. The possession of nar, nap, nir, nor, and nos genes by most Azoarcus-Aromatoleum strains suggests that they have the potential to derive energy through anaerobic nitrate respiration, so this ability cannot be usefully used as a phenotypic marker to distinguish genera. However, the possession of bzd genes indicating the ability to degrade benzoate anaerobically plus the type of diazotrophy (aerobic vs. anaerobic) could, after confirmation of their functionality, be considered as distinguishing phenotypes in any new generic delineations. The taxonomy of the Azoarcus-Aromatoleum group should be revisited; retaining the generic name Azoarcus for its entirety, or creating additional genera are both possible outcomes.