The anaerobic ammonium oxidation (anammox) bacteria can transform ammonium and nitrite to dinitrogen gas, and this obligate anaerobic process accounts for up to half of the global nitrogen loss in surface environments. Yet its origin and evolution, which may give important insights into the biogeochemistry of early Earth, remains enigmatic. Here, we performed comprehensive phylogenomic and molecular clock analysis of anammox bacteria within the phylum Planctomycetes. After accommodating the uncertainties and factors influencing time estimates, which includes implementing both a traditional cyanobacteria-based and a recently developed mitochondria-based molecular dating approach, we estimated a consistent origin of anammox bacteria at early Proterozoic and most likely around the so-called Great Oxidation Event (GOE; 2.32 to 2.5 billion years ago [Ga]) which fundamentally changed global biogeochemical cycles. We further showed that during the origin of anammox bacteria, genes involved in oxidative stress adaptation, bioenergetics and anammox granules formation were recruited, which might have contributed to their survival on an increasingly oxic Earth. Our findings suggest the rising levels of atmospheric oxygen, which made nitrite increasingly available, was a potential driving force for the emergence of anammox bacteria. This is one of the first studies that link the GOE to the evolution of obligate anaerobic bacteria.

Untangling the complex variations of microbiome associated with large-scale host phenotypes or environment types challenges the currently available analytic methods. Here, we present tmap, an integrative framework based on topological data analysis for population-scale microbiome stratification and association studies. The performance of tmap in detecting nonlinear patterns is validated by different scenarios of simulation, which clearly demonstrate its superiority over the most commonly used methods. Application of tmap to several population-scale microbiomes extensively demonstrates its strength in revealing microbiome-associated host or environmental features and in understanding the systematic interrelations among their association patterns. tmap is available at https://github.com/GPZ-Bioinfo/tmap.

Nitrate is one of the major inorganic nitrogen sources for microbes. Many bacterial and archaeal lineages have the capacity to express assimilatory nitrate reductase (NAS), which catalyzes the rate-limiting reduction of nitrate to nitrite. Although a nitrate assimilatory pathway in mycobacteria has been proposed and validated physiologically and genetically, the putative NAS enzyme has yet to be identified. Here, we report the characterization of a novel NAS encoded by Mycolicibacterium smegmatis Msmeg_4206, designated NasN, which differs from the canonical NASs in its structure, electron transfer mechanism, enzymatic properties, and phylogenetic distribution. Using sequence analysis and biochemical characterization, we found that NasN is an NADPH-dependent, diflavin-containing monomeric enzyme composed of a canonical molybdopterin cofactor-binding catalytic domain and an FMN–FAD/NAD-binding, electron-receiving/transferring domain, making it unique among all previously reported hetero-oligomeric NASs. Genetic studies revealed that NasN is essential for aerobic M. smegmatis growth on nitrate as the sole nitrogen source and that the global transcriptional regulator GlnR regulates nasN expression. Moreover, unlike the NADH-dependent heterodimeric NAS enzyme, NasN efficiently supports bacterial growth under nitrate-limiting conditions, likely due to its significantly greater catalytic activity and oxygen tolerance. Results from a phylogenetic analysis suggested that the nasN gene is more recently evolved than those encoding other NASs and that its distribution is limited mainly to Actinobacteria and Proteobacteria. We observed that among mycobacterial species, most fast-growing environmental mycobacteria carry nasN, but that it is largely lacking in slow-growing pathogenic mycobacteria because of multiple independent genomic deletion events along their evolution.

The alphaproteobacterial genus Bradyrhizobium has been best known as N2-fixing members that nodulate legumes, supported by the nif and nod gene clusters. Recent environmental surveys show that Bradyrhizobium represents one of the most abundant free-living bacterial lineages in the world’s soils. However, our understanding of Bradyrhizobium comes largely from symbiotic members, biasing the current knowledge of their ecology and evolution. Here, we report the genomes of 88 Bradyrhizobium strains derived from diverse soil samples, including both nif-carrying and non-nif-carrying free-living (nod free) members. Phylogenomic analyses of these and 252 publicly available Bradyrhizobium genomes indicate that nif-carrying free-living members independently evolved from symbiotic ancestors (carrying both nif and nod) multiple times. Intriguingly, the nif phylogeny shows that the vast majority of nif-carrying free-living members comprise an independent cluster, indicating that horizontal gene transfer promotes nif expansion among the free-living Bradyrhizobium. Comparative genomics analysis identifies that the nif genes found in free-living Bradyrhizobium are located on a unique genomic island of ~50 kb equipped with genes potentially involved in coping with oxygen tension. We further analyze amplicon sequencing data to show that Bradyrhizobium members presumably carrying this nif island are widespread in a variety of environments. Given the dominance of Bradyrhizobium in world’s soils, our findings have implications for global nitrogen cycles and agricultural research.

The anaerobic ammonium oxidation (anammox) bacteria transform ammonium and nitrite to dinitrogen gas, and this obligate anaerobic process accounts for nearly half of global nitrogen loss. Yet its origin and evolution, which may give important insights into the biogeochemistry of early Earth, remains enigmatic. Here, we compile a comprehensive sequence data set of anammox bacteria, and confirm their single origin within the phylum Planctomycetes. After accommodating the uncertainties and factors influencing time estimates with different statistical methods, we estimate that anammox bacteria originated at around the so-called Great Oxidation Event (GOE; 2.32 to 2.5 billion years ago [Gya]) which is thought to have fundamentally changed global biogeochemical cycles. We further show that during the origin of anammox bacteria, genes involved in oxidative stress, bioenergetics and anammox granules formation were recruited, which may have contributed to their survival in an increasingly oxic Earth. Our findings suggest the rising level of atmospheric oxygen, which made nitrite increasingly available, as a potential driving force for the emergence of anammox bacteria. This is one of the first studies that link GOE to the evolution of obligate anaerobic bacteria.

Population-scale microbiome study poses specific challenges in data analysis, from enterotype analysis, identification of driver species, to microbiome-wide association of host covariates. Application of advanced data mining techniques to high-dimensional complex dataset is expected to meet the rapid advancement in large scale and integrative microbiome research. Here, we present tmap, a topological data analysis framework for population-scale microbiome study. This framework can capture complex shape of large scale microbiome data into a compressive network representation. We also develop network-based statistical analysis for driver species identification and microbiome-wide association analysis. tmap can be used for exploring variations in a population-scale microbiome landscape to study host-microbiome association.Availability and implementationtmap is available at GitHub (https://github.com/GPZ-Bioinfo/tmap), accompanied with online documentation and tutorial (http://tmap.readthedocs.io).Contacthttp://hk.zhou@siat.ac.cn

The alphaproteobacterial genus Bradyrhizobium has been best known as N2-fixing members that nodulate legumes, supported by the nif and nod gene clusters. Recent environmental surveys show that Bradyrhizobium represents one of the most abundant free-living bacterial lineages in the world's soils. However, our understanding of Bradyrhizobium comes largely from symbiotic members, biasing the current knowledge of their ecology and evolution. Here, we report the genomes of 88 Bradyrhizobium strains derived from diverse soil samples, including both nif-carrying and non-nif-carrying free-living (nod free) members. Phylogenomic analyses of these and 252 publicly available Bradyrhizobium genomes indicate that nif-carrying free-living members independently evolved from symbiotic ancestors (carrying both nif and nod) multiple times. Intriguingly, the nif phylogeny shows that all nif-carrying free-living members comprise a cluster which branches off earlier than most symbiotic lineages. These results indicate that horizontal gene transfer (HGT) promotes nif expansion among the free-living Bradyrhizobium and that the free-living nif cluster represents a more ancestral version compared to that in symbiotic lineages. Further evidence for this rampant HGT is that the nif in free-living members consistently co-locate with several important genes involved in coping with oxygen tension which are missing from symbiotic members, and that while in free-living Bradyrhizobium nif and the co-locating genes show a highly conserved gene order, they each have distinct genomic context. Given the dominance of Bradyrhizobium in world's soils, our findings have implications for global nitrogen cycles and agricultural research.

No abstract