Evolutionary origin and ecological implication of a unique nif island in free-living Bradyrhizobium lineages

Biodegradable microplastics affect tomato (Solanum lycopersicum L.) growth by interfering rhizosphere key phylotypes

Biodegradable microplastics (BMPs), which form as biodegradable plastics degrade in agricultural settings, may influence plant growth and soil health. This study investigates the effects of BMPs on tomato growth and the microbial mechanisms involved. A greenhouse experiment applied BMPs-polyhydroxyalkanoate (PHA), polylactic acid (PLA), poly(butylene succinate-co-butylene adipate) (PBSA), and poly(butylene-adipate-co-terephthalate) (PBAT)-to tomato plants. The study analyzed their effects on plant growth, soil properties, and rhizosphere microbial communities. BMP treatments significantly reduced tomato biomass, height, and chlorophyll content compared to the control. PLA0.1 decreased the chlorophyll a/b ratio, while PLA1 increased it. Elemental analysis showed PLA1 increased phosphorus, calcium, and potassium in leaves, whereas all BMPs reduced nitrogen levels. BMPs also altered soil nitrogen and DOC levels, significantly shifting rhizosphere microbial communities, with a notable increase in Betaproteobacteria abundance. Ecological network analysis revealed that BMPs disrupted key microbial modules linked to plant growth. Beneficial modules positively associated with biomass and nutrient uptake were reduced under BMP treatments, whereas harmful microbial taxa in module 3, associated to poor plant health, were promoted. These shifts suggest that BMPs disrupt microbial ecological relationships critical for optimal plant growth. The findings highlight the potential negative impacts of BMPs on tomato growth through changes in microbial dynamics and soil properties.

Metagenomic analysis revealed community-level metabolic differences between full-scale EBPR and S2EBPR systems

Side-Stream Enhanced Biological Phosphorus Removal (S2EBPR) has emerged as a promising technology addressing certain challenges of conventional Enhanced Biological Phosphorus Removal (EBPR), notably stability in phosphorus removal, yet the underlying mechanisms are not fully understood. Metagenomic analysis presents a powerful approach to elucidate community-level metabolic differences between EBPR and S2EBPR configurations. In this study, we compared three EBPR and three S2EBPR activated sludge communities using metagenomic analysis at taxonomy, key functional pathways/genes, and polyphosphate-metabolism marker genes. Our analysis revealed larger genus-level diversity variance in S2EBPR communities, indicating distinct microbial community compositions influenced by different operational configurations. A higher diversity index in the S2EBPR than the EBPR was observed, and a higher Ca. Accumulibacter abundance was detected in EBPRs, whereas the fermentative candidate PAOs genera, including Ca. Phosphoribacter and Ca. Promineifilum, were more abundant in S2EBPR systems. EBPR and S2EBPR groups displayed similar gene and pathway abundance patterns related to core metabolisms essential for carbon and nitrogen metabolism. PolyP-metabolism marker gene phylogeny analysis suggested that exopolyphosphatase gene (ppx) showed better distinctions between EBPR and S2EBPR communities than polyphosphate kinase gene (ppk). This also highlighted the needs in fine-cale microdiversity analysis and finding novel Ca. Accumulibacter clades and species as resolved using the ppk gene. These findings provide valuable insights into AS community dynamics and metabolic functionalities, paving the way for further research into optimizing phosphorus removal processes in wastewater treatment systems.

High-quality PacBio draft genome sequences of 17 free-living Bradyrhizobium and four related Nitrobacteraceae strains isolated from arid soils in the Santa Catalina Mountains of Southern Arizona

Non-symbiotic Bradyrhizobium are among the most abundant and ubiquitous microbes in bulk soils globally. Despite this, most available genomic resources for Bradyrhizobium are derived from plant-associated strains. We present high-quality draft genomes for 17 Bradyrhizobium and four Nitrobacteraceae cultures isolated from bulk semiarid soils in Arizona, USA. The genome sizes range from 5.99 to 10.4 Mbp. Phylogenomic analysis of the 21 genomes indicates they fall into four clades. Two of the clades are nested within the Bradyrhizobium genus. The other two clades were associated with Nitrobacteraceae outgroups basal to Bradyrhizobium. All genomes lack genes coding for molybdenum or vanadium nitrogenases, and nod genes that code for proteins involved in nodulation, suggesting these isolates are free-living, non-symbiotic and do not fix dinitrogen gas. These genomes offer new resources for investigating free-living Bradyrhizobium lineages.

Correlating phylogenetic and functional diversity of the nod-free but nodulating Bradyrhizobium phylogroup

Bradyrhizobium is a main rhizobial lineage of which most members nodulate legume plants using Nod factors synthetized by the nod genes. However, members of the Photosynthetic supergroup (phylogroup) within Bradyrhizobium are nod-free, but still capable of establishing nitrogen-fixing nodules with some tropical legumes of the Aeschynomene genus. These unusual findings are based on the genomic sequences of only 13 Photosynthetic Bradyrhizobium strains, and almost all were isolated from Aeschynomene nodules. Here, we report that Photosynthetic Bradyrhizobium supergroup members are more abundantly associated with rice root (endosphere and rhizosphere) compared to grassland, forest, and maize samples based on rpoB amplicon sequence analyses. We sequenced 263 new isolates of this supergroup mostly from two main subspecies of cultivated rice (Oryza sativa L. spp. indica and japonica). The extended supergroup comprises three major clades with their diversity broadly covering the natural community of this supergroup: a basal clade with significant expansion of its diversity, a clade composed by two phylogenetically diverse strains including one newly isolated, and a new clade exclusively represented by our new strains. Although this supergroup members universally lack the canonical nod genes, all 28 assayed strains covering the broad diversity induced nodules on Aeschynomene indica. The three clades displayed important differences in the efficiency of symbiosis, aligning well with their phylogenetic divergence. With this expanded ecological, phylogenetic, and functional diversity, we conclude that the nod factor-independent nodulation of Aeschynomene is a common trait of this supergroup, in contrast to the photosynthetic trait originally thought of as its unifying feature.

Great diverse rhizobial community nodulating Astragalus mongholicus in the northeastern region of China

Introduction

Astragalus mongholicus Bunge is an important medicinal legume species widely cultivated in northeastern China (NEC) and northwestern China (NWC) and can establish a symbiotic relationship with nitrogen-fixing rhizobial strains. However, there are limited reports comparing the genetic diversity, differentiation, and gene flow of rhizobial strains associated with this plant in different geographic regions.

Methods

We used multilocus sequence analysis (MLSA) to investigate the phylogeny and genetic diversity of rhizobia and to estimate their intra- and inter-regional gene flow and genetic differentiation based on the analysis of concatenated core genes (recA, atpD, and glnII) and the critical symbiotic gene nodC.

Results

We isolated eight known and three novel genospecies representing four genera, among which Rhizobium yanglingense was the most predominant microsymbiont. Phylogenetic analysis revealed a highly diverse rhizobial community nodulating Astragalus mongholicus in NEC, consisting of the four genera Rhizobium, Bradyrhizobium, Sinorhizobium, and Mesorhizobium. This community differed markedly from the rhizobial community found in NWC. Various rhizobial genospecies with different symbiotic gene nodC sequences were capable of nodulating A. mongholicus in NEC. Therefore, A. mongholicus exhibits promiscuity in its association with symbionts in the natural environment, showing no strong preference for either the species-defining core genes or the symbiotic genes of rhizobia. We also found that the Glyco_tranf_GTA_type superfamily (Glycosyltransferase family A) is the most highly conserved and essential domain in the NodC protein, which is encoded by the symbiotic nodC gene, across nodulating rhizobia. In addition, we found independent genetic differentiation among rhizobial communities geographically, and the frequency of gene flow among microsymbionts between NEC and NWC was low. We speculate that the formation of the highly diverse rhizobial community in NEC resulted from the independent evolution of each ancestral lineage. This diversity likely arose from intraregional genetic differentiation driven by mutations rather than recombination.

Conclusion

Ecogeographical isolation between NEC and NWC restricted inter-regional genetic drift and gene flow. Therefore, intraregional genetic differentiation is the major evolutionary force underlying the genetic diversity of rhizobia.

Dry season irrigation promotes nutrient cycling by reorganizing Eucalyptus rhizosphere microbiome

In southern China, seasonal droughts and low soil phosphorus content constrain the productivity of Eucalyptus trees. To understand the rhizosphere microbiome response to the dry season, metagenomic sequencing analysis was used to investigate the 6-year-old Eucalyptus rhizosphere microbiome under four different irrigation and fertilization treatments. The results showed that irrigation and fertilization during the dry season significantly altered the composition of microbiome in the rhizosphere soil of Eucalyptus plantations. The soil physicochemical properties and enzyme activity explained 30.73 % and 29.75 % of the changes in bacterial and fungal community structure in Eucalyptus rhizosphere soil, respectively. Irrigation and fertilization during the dry season significantly altered the physicochemical properties of rhizosphere soil. Compared with the seasonal drought without fertilizer treatment (CK), the dry season irrigation with fertilizer treatment (WF) significantly increased the content of total nitrogen (46.34 %), available nitrogen (37.72 %), available phosphorus (440.9 %), and organic matter (35.34 %). Soil organic matter (OM), pH, and available phosphorus (AP) were key environmental factors influencing the microbial community composition. Moreover, irrigation and fertilization promoted carbon fixation and nitrogen and phosphorus mineralization, increasing soil OM content and the availability of inorganic nitrogen and phosphorus. Meanwhile, compared to the CK, the increase of acid phosphatase (16.81 %), invertase (146.89 %)and urease (59.45 %) in rhizosphere soil under irrigation (W) treatment further proves that dry season irrigation promote the soil carbon, nitrogen and phosphorus cycles. Irrigation and fertilization treatment alleviated the constraints of low phosphorus in southern China's soil, which promoted Eucalyptus productivity. In conclusion, we suggest implementing reasonable irrigation and fertilization strategies in the production practice of Eucalyptus and utilizing microbial resources to improve soil fertility and Eucalyptus productivity.

Microbial community composition and co-occurrence network analysis of the rhizosphere soil of the main constructive tree species in Helan Mountain of Northwest China

To understand the microbial diversity and community composition within the main constructive tree species, Picea crassifolia, Betula platyphylla, and Pinus tabuliformis, in Helan Mountain and their response to changes in soil physicochemical factors, a high throughput sequencing technology was used to analyze the bacterial and fungal diversity and community structure. RDA (Redundancy Analysis) and Pearson correlation analysis were used to explore the influence of soil physicochemical factors on microbial community construction, and co-occurrence network analysis was conducted on the microbial communities. The results showed that the fungal and bacterial diversity was highest in B. platyphylla, and lowest in P. crassifolia. Additionally, the fungal/bacterial richness was greatest in the rhizosphere soils of P. tabuliformis and B. platyphylla. RDA and Pearson correlation analysis revealed that NN (nitrate nitrogen) and AP (available phosphorus) were the main determining factors of the bacterial community, while NN and SOC (soil water content) were the main determining factors of the fungal community. Pearson correlation analysis between soil physicochemical factors and the alpha diversity of the microbial communities revealed a significant positive correlation between pH and the bacterial and fungal diversity, while SOC, TN (total nitrogen), AP, and AN (available nitrogen) were significantly negatively correlated with the bacterial and fungal diversity. Co-occurrence network analysis revealed that the soil bacterial communities exhibit richer network nodes, edges, greater diversity, and greater network connectivity. Indicating that bacterial communities exhibit more complex and stable interaction patterns in soil. This study reveals the complex interactive relationship between microbial communities and soil physicochemical factors in forest ecosystems. By analyzing the response of rhizosphere microbial communities of major tree species in Helan Mountain to nutrient dynamics and pH changes, we can deepen our understanding of the role of microorganisms in regulating ecosystem functions and provide theoretical basis for soil improvement and ecological restoration strategies.

Rhizosphere microbial community structure in high-producing, low-input switchgrass families

Switchgrass (Panicum virgatum L.) is a native, low-input North American perennial crop primarily grown for bioenergy, livestock forage, and industrial fiber. To achieve no-input switchgrass production that meets biomass needs, several switchgrass genotypes have been identified that have a low or negative response to nitrogen fertilizer, i.e., the biomass accumulation with added nitrogen is less than or equal to that when grown without nitrogen. In order to improve the viability of low-input switchgrass production, a more detailed understanding of the biogeochemical mechanisms active in these select genotypes is needed. 16S and ITS amplicon sequencing and qPCR of key functional genes were applied to switchgrass rhizospheres to elucidate microbial community structure in high-producing, no-input switchgrass families. Rhizosphere microbial community structure differed strongly between sites, and nitrogen responsiveness.

Species pool and soil properties in mangrove habitats influence the species-immigration process of diazotrophic communities across southern China

Microbial immigration is an ecological process in natural environments; however, the ecological trade-off mechanisms that govern the balance between species extinction and migration are still lacking. In this study, we investigated the mechanisms underlying the migration of diazotrophic communities from soil to leaves across six natural mangrove habitats in southern China. The results showed that the diazotrophic alpha and beta diversity exhibited significant regional and locational variations. The diazotrophic species pool gradually increased from the leaves to nonrhizosphere soil at each site, exhibiting a vertical distribution pattern. Mantel test analyses suggested that climate factors, particularly mean annual temperature, significantly influenced the structure of the diazotrophic community. The diazotrophic community assembly was mainly governed by dispersal limitation in soil and root samples, whereas dispersal limitation and ecological drift were dominant in leaves. Partial least squares path modeling revealed that the species pool and soil properties, particularly the oxidation-reduction potential and pH, were closely linked to the species-immigration ratio of diazotrophic communities. Our study provides novel insights for understanding the ecological trait diversity patterns and spread pathways of functional microbial communities between below- and aboveground habitats in natural ecosystems.IMPORTANCEEnvironmental selection plays key roles in microbial transmission. In this study, we have provided a comprehensive framework to elucidate the driving patterns of the ecological trade-offs in diazotrophic communities across large-scale mangrove habitats. Our research revealed that Bradyrhizobium japonicum, Marinobacterium lutimaris, and Agrobacterium tumefaciens were more abundant in root-associated soil than in leaves by internal and external pathways. The nonrhizospheric and rhizospheric soil samples harbored the most core amplicon sequence variants, indicating that these dominant diazotrophs could adapt to broader ecological niches. Correlation analysis indicated that the diversities of the diazotrophic community were regulated by biotic and abiotic factors. Furthermore, this study found a lower species immigration ratio in the soil than in the leaves. Both species pool and soil properties regulate the species-immigration mechanisms of the diazotrophic community. These results suggest that substantial species immigration is a widespread ecological process, leading to alterations in local community diversity across diverse host environments.

Effects of humic acid fertilizer on the growth and microbial network stability of Panax notoginseng from the forest understorey

Humic acid (HA) can substantially enhance plant growth and improve soil health. Currently, the impacts of HA concentrations variation on the development and soil quality of Panax notoginseng (Sanqi) from the forest understorey are still unclear. In this study, exogenous HA was administered to the roots of Sanqi at varying concentrations (2, 4, and 6 ml/L). Subsequently, the diversity and community structure of bacteria and fungi were assessed through high-throughput sequencing technology. The investigation further involved analyzing the interplay among the growth of sanqi, soil edaphic factors, and the microbial network stability. Our finding revealed that moderate concentrations (4 ml/L) of HA improved the fresh/dry weight of Sanqi and NO3--N levels. Compared with control, the moderate concentrations of HA had a notable impact on the bacterial and fungal communities compositions. However, there was no significant difference in the α and β diversity of bacteria and fungi. Moreover, the abundance of beneficial bacteria (Bradyrhizobium) and harmful bacteria (Xanthobacteraceae) increased and decreased at 4 ml/L HA, respectively, while the bacterial and fungal network stability were enhanced. Structural equation model (SEM) revealed that the fresh weight of Sanqi and bacterial and fungal communities were the factors that directly affected the microbial network stability at moderate concentrations of HA. In conclusion, 4 ml/L of HA is beneficial for promoting Sanqi growth and soil quality. Our study provides a reference for increasing the yield of Sanqi and sustainable development of the Sanqi-pine agroforestry system.

Optimizing intermittent micro-aeration as a strategy for enhancing aniline anaerobic biodegradation: kinetic, ecotoxicity, and microbial community dynamics analyses

Groundwater and soil contamination by aromatic amines (AAs), used in the production of polymers, plastics, and pesticides, often results from improper waste disposal and accidental leaks. These compounds are resistant to anaerobic degradation; however, micro-aeration can enhance this process by promoting microbial interactions. In batch assays, anaerobic degradation of aniline (0.14 mM), a model AA, was tested under three micro-aeration conditions: T30, T15, and T10 (30, 15, and 10 min of micro-aeration every 2 h, respectively). Aniline degradation occurred in all conditions, producing both aerobic (catechol) and anaerobic (benzoic acid) byproducts. The main genera involved in T30 and T15 were Comamonas, Clostridium, Longilinea, Petrimonas, Phenylobacterium, Pseudoxanthomonas, and Thiobacillus. In contrast, in T10 were Pseudomonas, Delftia, Leucobacter, and Thermomonas. While T30 and T15 promoted microbial cooperation for anaerobic degradation and facultative respiration, T10 resulted in a competitive environment due to dominance and oxygen scarcity. Despite aniline degradation in 9.4 h under T10, this condition was toxic to Allium cepa seeds and exhibited cytogenotoxic effects. Therefore, T15 emerged as the optimal condition, effectively promoting anaerobic degradation without accumulating toxic byproducts. Intermittent micro-aeration emerges as a promising strategy for enhancing the anaerobic degradation of AA-contaminated effluents.

Fertilization regime changes rhizosphere microbial community assembly and interaction in Phoebe bournei plantations

Fertilizer input is one of the effective forest management practices, which improves soil nutrients and microbial community compositions and promotes forest productivity. However, few studies have explored the response of rhizosphere soil microbial communities to various fertilization regimes across seasonal dynamics. Here, we collected the rhizosphere soil samples from Phoebe bournei plantations to investigate the response of community assemblages and microbial interactions of the soil microbiome to the short-term application of four typical fertilizer practices (including chemical fertilizer (CF), organic fertilizer (OF), compound microbial fertilizer (CMF), and no fertilizer control (CK)). The amendments of organic fertilizer and compound microbial fertilizer altered the composition of rhizosphere soil bacterial and fungal communities, respectively. The fertilization regime significantly affected bacterial diversity rather than fungal diversity, and rhizosphere fungi responded more sensitively than bacteria to season. Fertilization-induced fungal networks were more complex than bacterial networks. Stochastic processes governed both rhizosphere soil bacterial and fungal communities, and drift and dispersal limitation dominated soil fungal and bacterial communities, respectively. Collectively, these findings demonstrate contrasting responses to community assemblages and interactions of rhizosphere bacteria and fungi to fertilizer practices. The application of organic fertilization strengthens microbial interactions and changes the succession of key taxa in the rhizosphere habitat. KEY POINTS: • Fertilization altered the key taxa and microbial interaction • Organic fertilizer facilitated the turnover of rhizosphere microbial communities • Stochasticity governed soil fungal and bacterial community assembly.

Arbuscular mycorrhizal fungal interactions bridge the support of root-associated microbiota for slope multifunctionality in an erosion-prone ecosystem

The role of diverse soil microbiota in restoring erosion-induced degraded lands is well recognized. Yet, the facilitative interactions among symbiotic arbuscular mycorrhizal (AM) fungi, rhizobia, and heterotrophic bacteria, which underpin multiple functions in eroded ecosystems, remain unclear. Here, we utilized quantitative microbiota profiling and ecological network analyses to explore the interplay between the diversity and biotic associations of root-associated microbiota and multifunctionality across an eroded slope of a Robinia pseudoacacia plantation on the Loess Plateau. We found explicit variations in slope multifunctionality across different slope positions, associated with shifts in limiting resources, including soil phosphorus (P) and moisture. To cope with P limitation, AM fungi were recruited by R. pseudoacacia, assuming pivotal roles as keystones and connectors within cross-kingdom networks. Furthermore, AM fungi facilitated the assembly and composition of bacterial and rhizobial communities, collectively driving slope multifunctionality. The symbiotic association among R. pseudoacacia, AM fungi, and rhizobia promoted slope multifunctionality through enhanced decomposition of recalcitrant compounds, improved P mineralization potential, and optimized microbial metabolism. Overall, our findings highlight the crucial role of AM fungal-centered microbiota associated with R. pseudoacacia in functional delivery within eroded landscapes, providing valuable insights for the sustainable restoration of degraded ecosystems in erosion-prone regions.

Precipitation affects soil nitrogen fixation by regulating active diazotrophs and nitrate nitrogen in an alpine grassland of Qinghai-Tibetan Plateau

Soil asymbiotic nitrogen (N) fixation provides a critical N source to support plant growth in alpine grasslands, and precipitation change is expected to lead to shifts in soil asymbiotic N fixation. However, large gaps remain in understanding the response of soil asymbiotic N fixation to precipitation gradients. Here we simulated five precipitation gradients (10 % (0.1P), 50 % (0.5P), 70 % (0.7P), 100 % (1.0P) and 150 % (1.5P) of the natural precipitation) in an alpine grassland of Qinghai-Tibetan Plateau and examined the soil nitrogenase activity and N fixation rate for each gradient. Quantitative PCR and high-throughput sequencing were used to measure the abundance and community composition of the soil nifH DNA (total diazotrophs) and nifH RNA reverse transcription (active diazotrophs) gene. Our results showed that the soil diazotrophic abundance, diversity and nifH gene expression rate peaked under the 0.5P. Soil nitrogenase activity and N fixation rate varied in the range 0.032-0.073 nmol·C2H4·g-1·h-1 and 0.008-0.022 nmol·N2·g-1·h-1 respectively, being highest under the 0.5P. The 50 % precipitation reduction enhanced the gene expression rates of Azospirillum and Halorhodospira which were likely responsible for the high N fixation potential. The 0.5P treatment also possessed a larger and more complex active diazotrophic network than the other treatments, which facilitated the resistance of diazotrophic community to environmental stress and thus maintained a high N fixation potential. The active diazotrophic abundance had the largest positive effect on soil N fixation, while nitrate nitrogen had the largest negative effect. Together, our study suggested that appropriate precipitation reduction can enhance soil N fixation through promoting the abundance of the soil active diazotrophs and decreasing soil nitrate nitrogen, and soil active diazotrophs and nitrate nitrogen should be considered in predicting soil N inputs in the alpine grassland of Qinghai-Tibetan Plateau under precipitation change.

Comparative genomics analysis reveals genetic characteristics and nitrogen fixation profile of Bradyrhizobium

Bradyrhizobium is a genus of nitrogen-fixing bacteria, with some species producing nodules in leguminous plants. Investigations into Bradyrhizobium have recently revealed its substantial genetic resources and agricultural benefits, but a comprehensive survey of its genetic diversity and functional properties is lacking. Using a panel of various strains (N = 278), this study performed a comparative genomics analysis to anticipate genes linked with symbiotic nitrogen fixation. Bradyrhizobium's pan-genome consisted of 84,078 gene families, containing 824 core genes and 42,409 accessory genes. Core genes were mainly involved in crucial cell processes, while accessory genes served diverse functions, including nitrogen fixation and nodulation. Three distinct genetic profiles were identified based on the presence/absence of gene clusters related to nodulation, nitrogen fixation, and secretion systems. Most Bradyrhizobium strains from soil and non-leguminous plants lacked major nif/nod genes and were evolutionarily more closely related. These findings shed light on Bradyrhizobium's genetic features for symbiotic nitrogen fixation.

Genetic diversity into a novel free-living species of Bradyrhizobium from contaminated freshwater sediment

A free-living Bradyrhizobium strain isolated from a contaminated sediment sample collected at a water depth of 4 m from the Hongze Lake in China was characterized. Phylogenetic investigation of the 16S rRNA gene, concatenated housekeeping gene sequences, and phylogenomic analysis placed this strain in a lineage distinct from all previously described Bradyrhizobium species. The sequence similarities of the concatenated housekeeping genes support its distinctiveness with the type strains of the named species. The complete genome of strain S12-14-2 consists of a single chromosome of size 7.3M. The strain lacks both a symbiosis island and important nodulation genes. Based on the data presented here, the strain represents a new species, for which the name Bradyrhizobium roseus sp. nov. is proposed for the type strain S12-14-2T. Several functional differences between the isolate and other published genomes indicate that the genus Bradyrhizobium is extremely heterogeneous and has functions within the community, such as non-symbiotic nitrogen fixation. Functional denitrification and nitrogen fixation genes were identified on the genomes of strain S12-14-2T. Genes encoding proteins for sulfur oxidation, sulfonate transport, phosphonate degradation, and phosphonate production were also identified. Lastly, the B. roseus genome contained genes encoding ribulose 1,5-bisphosphate carboxylase/oxygenase, a trait that presumably enables autotrophic flexibility under varying environmental conditions. This study provides insights into the dynamics of a genome that could enhance our understanding of the metabolism and evolutionary characteristics of the genus Bradyrhizobium and a new genetic framework for future research.

Marine nitrogen-fixers in the Canadian Arctic Gateway are dominated by biogeographically distinct noncyanobacterial communities

We describe diazotrophs present during a 2015 GEOTRACES expedition through the Canadian Arctic Gateway (CAG) using nifH metabarcoding. In the less studied Labrador Sea, Bradyrhizobium sp. and Vitreoscilla sp. nifH variants were dominant, while in Baffin Bay, a Stutzerimonas stutzeri variant was dominant. In comparison, the Canadian Arctic Archipelago (CAA) was characterized by a broader set of dominant variants belonging to Desulfobulbaceae, Desulfuromonadales, Arcobacter sp., Vibrio spp., and Sulfuriferula sp. Although dominant diazotrophs fell within known nifH clusters I and III, only a few of these variants were frequently recovered in a 5-year weekly nifH times series in the coastal NW Atlantic presented herein, notably S. stutzeri and variants belonging to Desulfobacterales and Desulfuromonadales. In addition, the majority of dominant Arctic nifH variants shared low similarity (< 92% nucleotide identities) to sequences in a global noncyanobacterial diazotroph catalog recently compiled by others. We further detected UCYN-A throughout the CAG at low-levels using quantitative-PCR assays. Temperature, depth, salinity, oxygen, and nitrate were most strongly correlated to the Arctic diazotroph diversity observed, and we found a stark division between diazotroph communities of the Labrador Sea versus Baffin Bay and the CAA, hence establishing that a previously unknown biogeographic community division can occur for diazotrophs in the CAG.

Complex evolutionary history of photosynthesis in Bradyrhizobium

Bradyrhizobium comprises a diverse group of bacteria with various lifestyles. Although best known for their nodule-based nitrogen-fixation in symbiosis with legumes, a select group of bradyrhizobia are also capable of photosynthesis. This ability seems to be rare among rhizobia, and its origin and evolution in these bacteria remain a subject of substantial debate. Therefore, our aim here was to investigate the distribution and evolution of photosynthesis in Bradyrhizobium using comparative genomics and representative genomes from closely related taxa in the families Nitrobacteraceae, Methylobacteriaceae, Boseaceae and Paracoccaceae . We identified photosynthesis gene clusters (PGCs) in 25 genomes belonging to three different Bradyrhizobium lineages, notably the so-called Photosynthetic, B. japonicum and B. elkanii supergroups. Also, two different PGC architectures were observed. One of these, PGC1, was present in genomes from the Photosynthetic supergroup and in three genomes from a species in the B. japonicum supergroup. The second cluster, PGC2, was also present in some strains from the B. japonicum supergroup, as well as in those from the B. elkanii supergroup. PGC2 was largely syntenic to the cluster found in Rhodopseudomonas palustris and Tardiphaga . Bayesian ancestral state reconstruction unambiguously showed that the ancestor of Bradyrhizobium lacked a PGC and that it was acquired horizontally by various lineages. Maximum-likelihood phylogenetic analyses of individual photosynthesis genes also suggested multiple acquisitions through horizontal gene transfer, followed by vertical inheritance and gene losses within the different lineages. Overall, our findings add to the existing body of knowledge on Bradyrhizobium ’s evolution and provide a meaningful basis from which to explore how these PGCs and the photosynthesis itself impact the physiology and ecology of these bacteria.

Natural variation of GmRj2/Rfg1 determines symbiont differentiation in soybean

Symbiotic nitrogen fixation (SNF) provides much of the N utilized by leguminous plants throughout growth and development. Legumes may simultaneously establish symbiosis with different taxa of microbial symbionts. Yet, the mechanisms used to steer associations toward symbionts that are most propitious across variations in soil types remain mysterious. Here, we demonstrate that GmRj2/Rfg1 is responsible for regulating symbiosis with multiple taxa of soybean symbionts. In our experiments, the GmRj2/Rfg1SC haplotype favored association with Bradyrhizobia, which is mostly distributed in acid soils, whereas the GmRj2/Rfg1HH haplotype and knockout mutants of GmRj2/Rfg1SC associated equally with Bradyrhizobia and Sinorhizobium. Association between GmRj2/Rfg1 and NopP, furthermore, appeared to be involved in symbiont selection. Furthermore, geographic distribution analysis of 1,821 soybean accessions showed that GmRj2/Rfg1SC haplotypes were enriched in acidic soils where Bradyrhizobia were the dominant symbionts, whereas GmRj2/Rfg1HH haplotypes were most prevalent in alkaline soils dominated by Sinorhizobium, and neutral soils harbored no apparent predilections toward either haplotype. Taken together, our results suggest that GmRj2/Rfg1 regulates symbiosis with different symbionts and is a strong determinant of soybean adaptability across soil regions. As a consequence, the manipulation of the GmRj2/Rfg1 genotype or application of suitable symbionts according to the haplotype at the GmRj2/Rfg1 locus might be suitable strategies to explore for increasing soybean yield through the management of SNF.

Metagenomics analysis reveals differences in rumen microbiota in cows with low and high milk protein percentage

Variation exists in milk protein concentration of dairy cows of the same breed that are fed and managed in the same environment, and little information was available on this variation which might be attributed to differences in rumen microbial composition as well as their fermentation metabolites. This study is aimed at investigating the difference in the composition and functions of rumen microbiota as well as fermentation metabolites in Holstein cows with high and low milk protein concentrations. In this study, 20 lactating Holstein cows on the same diet were divided into two groups (10 cows each), high degree of milk protein group (HD), and low degree of milk protein (LD) concentrations based on previous milk composition history. Rumen content samples were obtained to explore the rumen fermentation parameters and rumen microbial composition. Shotgun metagenomics sequencing was employed to investigate the rumen microbial composition and sequences were assembled via the metagenomics binning technique. Metagenomics revealed that 6 Archaea genera, 5 Bacteria genera, 7 Eukaryota genera, and 7 virus genera differed significantly between the HD and LD group. The analysis of metagenome-assembled genomes (MAGs) showed that 2 genera (g__Eubacterium_H and g__Dialister) were significantly enriched (P < 0.05, linear discriminant analysis (LDA) > 2) in the HD group. However, the LD group recorded an increased abundance (P < 0.05, LDA > 2) of 8 genera (g__CAG-603, g__UBA2922, g__Ga6A1, g__RUG13091, g__Bradyrhizobium, g__Sediminibacterium, g__UBA6382, and g__Succinivibrio) when compared to the HD group. Furthermore, investigation of the KEGG genes revealed an upregulation in a higher number of genes associated with nitrogen metabolism and lysine biosynthesis pathways in the HD group as compared to the LD group. Therefore, the high milk protein concentration in the HD group could be explained by an increased ammonia synthesis by ruminal microbes which were converted to microbial amino acids and microbial protein (MCP) in presence of an increased energy source made possible by higher activities of carbohydrate-active enzymes (CAZymes). This MCP gets absorbed in the small intestine as amino acids and might be utilized for the synthesis of milk protein. KEY POINTS: • Rumen microbiota and their functions differed between cows with high milk protein % and those with low milk protein %. • The rumen microbiome of cows with high milk protein recorded a higher number of enriched genes linked to the nitrogen metabolism pathway and lysine biosynthesis pathway. • The activities of carbohydrate-active enzymes were found to be higher in the rumen of cows with high milk protein %.

Warmer winters result in reshaping of the European beech forest soil microbiome (bacteria, archaea and fungi)-With potential implications for ecosystem functioning

In temperate regions, climate warming alters temperature and precipitation regimes. During winter, a decline in insulating snow cover changes the soil environment, where especially frost exposure can have severe implications for soil microorganisms and subsequently for soil nutrient dynamics. Here, we investigated winter climate change responses in European beech forests soil microbiome. Nine study sites with each three treatments (snow exclusion, insolation, and ambient) were investigated. Long-term adaptation to average climate was explored by comparing across sites. Triplicated treatment plots were used to evaluate short-term (one single winter) responses. Community profiles of bacteria, archaea and fungi were created using amplicon sequencing. Correlations between the microbiome, vegetation and soil physicochemical properties were found. We identify core members of the forest-microbiome and link them to key processes, for example, mycorrhizal symbiont and specialized beech wood degraders (fungi) and nitrogen cycling (bacteria, archaea). For bacteria, the shift of the microbiome composition due to short-term soil temperature manipulations in winter was similar to the community differences observed between long-term relatively cold to warm conditions. The results suggest a strong link between the changes in the microbiomes and changes in environmental processes, for example, nitrogen dynamics, driven by variations in winter climate.

Adaptive Evolution of Rhizobial Symbiosis beyond Horizontal Gene Transfer: From Genome Innovation to Regulation Reconstruction

There are ubiquitous variations in symbiotic performance of different rhizobial strains associated with the same legume host in agricultural practices. This is due to polymorphisms of symbiosis genes and/or largely unexplored variations in integration efficiency of symbiotic function. Here, we reviewed cumulative evidence on integration mechanisms of symbiosis genes. Experimental evolution, in concert with reverse genetic studies based on pangenomics, suggests that gain of the same circuit of key symbiosis genes through horizontal gene transfer is necessary but sometimes insufficient for bacteria to establish an effective symbiosis with legumes. An intact genomic background of the recipient may not support the proper expression or functioning of newly acquired key symbiosis genes. Further adaptive evolution, through genome innovation and reconstruction of regulation networks, may confer the recipient of nascent nodulation and nitrogen fixation ability. Other accessory genes, either co-transferred with key symbiosis genes or stochastically transferred, may provide the recipient with additional adaptability in ever-fluctuating host and soil niches. Successful integrations of these accessory genes with the rewired core network, regarding both symbiotic and edaphic fitness, can optimize symbiotic efficiency in various natural and agricultural ecosystems. This progress also sheds light on the development of elite rhizobial inoculants using synthetic biology procedures.

Soil bacterial community changes along elevation gradients in karst graben basin of Yunnan-Kweichow Plateau

Elevation gradients could provide natural experiments to examine geomorphological influences on biota ecology and evolution, however little is known about microbial community structures with soil depths along altitudinal gradients in karst graben basin of Yunnan-Kweichow Plateau. Here, bulk soil in A layer (0 ~ 10 cm) and B layer (10 ~ 20 cm) from two transect Mounts were analyzed by using high-throughput sequencing coupled with physicochemical analysis. It was found that the top five phyla in A layer were Proteobacteria, Acidobacteria, Actinobacteria, Bacteroidetes, and Verrucomicrobia, and the top five phyla in B layer were Proteobacteria, Acidobacteria, Actinobacteria, Verrucomicrobia, and Chloroflexi in a near-neutral environment. Edaphic parameters were different in two layers along altitudinal gradients. Besides that, soil microbial community compositions varied along altitudinal gradient, and soil organic carbon (SOC) and total nitrogen (TN) increased monotonically with increasing elevation. It was further observed that Shannon indexes with increasing altitudes in two transect Mounts decreased monotonically with significant difference (p = 0.001), however beta diversity followed U-trend with significant difference (p = 0.001). The low proportions of unique operational taxonomic units (OTUs) appeared at high altitude areas which impact the widely accepted elevation Rapoport's rules. The dominant Bradyrhizobium (alphaproteobacterial OTU 1) identified at high altitudes in two layers constitutes the important group of free-living diazotrophs and could bring fixed N into soils, which simultaneously enhances SOC and TN accumulation at high altitudes (p < 0.01). Due to different responses of bacterial community to environmental changes varying with soil depths, altitudinal gradients exerted negative effects on soil bacterial communities via soil physical properties and positive effects on soil bacterial diversities via soil chemical properties in A layer, however the results in B layer were opposite. Overall, our study is the first attempt to bring a deeper understanding of soil microbial structure patterns along altitudinal gradients at karst graben basin areas.

Response of cbbL-harboring microorganisms to precipitation changes in a naturally-restored grassland

The impact of the long-term uneven precipitation distribution model on the diversity and community composition of soil C-fixing microorganisms in arid and semiarid grasslands remains unclear. In 2015, we randomly set up five experimental plots with precipitation gradients on the natural restoration grassland of the Loess Plateau (natural precipitation, NP; ± 40% natural precipitation: decreased precipitation (DP), DP40; increased precipitation (IP), IP40; ± 80% natural precipitation: DP80; IP80). In the third and fifth years after the experimental layout (spanned two years), we explored the cbbL-genes, which are functional genes in the Calvin cycle, harboring microbial diversity and community composition under different precipitation treatments. The results showed that the increase in mean annual precipitation significantly changed the cbbL-harboring microbial alpha diversity, especially when controlling for 40% natural precipitation. The response of the dominant microbial communities to interannual increased precipitation variation shifted from Gammaproteobacteria (Bradyrhizobium) to Betaproteobacteria (Variovorax). The structural equation model showed that precipitation directly affected the cbbL-harboring microbial diversity and community composition and indirectly by affecting soil NO₃^- (mg N kg ^-1), soil organic matter, dissolved organic N content, and above- and underground biomass. In conclusion, studying how cbbL-harboring microbial diversity and community composition respond to uneven precipitation variability provides new insights into the ecological processes of C-fixing microbes in semi-arid naturally-restored grasslands dominated by the Calvin cycle.

Assembly of root-associated N2O-reducing communities of annual crops is governed by selection for nosZ clade I over clade II

The rhizosphere is a hotspot for denitrification. The nitrous oxide (N2O) reductase among denitrifiers and nondenitrifying N2O reducers is the only known N2O sink in the biosphere. We hypothesized that the composition of root-associated N2O-reducing communities when establishing on annual crops depend on soil type and plant species, but that assembly processes are independent of these factors and differ between nosZ clades I and II. Using a pot experiment with barley and sunflower and two soils, we analyzed the abundance, composition, and diversity of soil and root-associated N2O reducing communities by qPCR and amplicon sequencing of nosZ. Clade I was more abundant on roots compared to soil, while clade II showed the opposite. In barley, this pattern coincided with N2O availability, determined as potential N2O production rates, but for sunflower no N2O production was detected in the root compartment. Root and soil nosZ communities differed in composition and phylogeny-based community analyses indicated that assembly of root-associated N2O reducers was driven by the interaction between plant and soil type, with inferred competition being more influential than habitat selection. Selection between clades I and II in the root/soil interface is suggested, which may have functional consequences since most clade I microorganisms can produce N2O.

Comparative Genomics Reveal the High Conservation and Scarce Distribution of Nitrogen Fixation nif Genes in the Plant-Associated Genus Herbaspirillum

The genus Herbaspirillum gained the spotlight due to the several reports of diazotrophic strains and promising results in plant-growth field assays. However, as diversity exploration of Herbaspirillum species gained momentum, it became clearer that the plant beneficial lifestyle was not the only form of ecological interaction in this genus, due to reports of phytopathogenesis and nosocomial infections. Here we performed a deep search across all publicly available Herbaspirillum genomes. Using a robust core genome phylogeny, we have found that all described species are well delineated, being the only exception H. aquaticum and H. huttiense clade. We also uncovered that the nif genes are only highly prevalent in H. rubrisubalbicans; however, irrespective to the species, all nif genes share the same gene arrangement with high protein identity, and are present in only two main types, in inverted strands. By means of a NifHDKENB phylogenetic tree, we have further revealed that the Herbaspirillum nif sequences may have been acquired from the same last common ancestor belonging to the Nitrosomonadales order.

Genotypic and symbiotic diversity studies of rhizobia nodulating Acacia saligna in Tunisia reveal two novel symbiovars within the Rhizobium leguminosarum complex and Bradyrhizobium

Acacia saligna is an invasive alien species that has the ability to establish symbiotic relationships with rhizobia. In the present study, genotypic and symbiotic diversity of native rhizobia associated with A. saligna in Tunisia were studied. A total of 100 bacterial strains were selected and three different ribotypes were identified based on rrs PCR-RFLP analysis. Sequence analyses of rrs and four housekeeping genes (recA, atpD, gyrB and glnII) assigned 30 isolates to four putative new lineages and a single strain to Sinorhizobium meliloti. Thirteen slow-growing isolates representing the most dominant IGS (intergenic spacer) profile clustered distinctly from known rhizobia species within Bradyrhizobium with the closest related species being Bradyrhizobium shewense and Bradyrhizobium niftali, which had 95.17% and 95.1% sequence identity, respectively. Two slow-growing isolates, 1AS28L and 5AS6L, had B. frederekii as their closest species with a sequence identity of 95.2%, an indication that these strains could constitute a new lineage. Strains 1AS14I, 1AS12I and 6AS6 clustered distinctly from known rhizobia species but within the Rhizobium leguminosarum complex (Rlc) with the most closely related species being Rhizobium indicum with 96.3% sequence identity. Similarly, the remaining 11 strains showed 96.9 % and 97.2% similarity values with R. changzhiense and R. indicum, respectively. Based on nodC and nodA phylogenies and cross inoculation tests, these 14 strains of Rlc species clearly diverged from strains of Sinorhizobium and Rlc symbiovars, and formed a new symbiovar for which the name sv. "salignae" is proposed. Bacterial strains isolated in this study that were taxonomically assigned to Bradyrhizobium harbored different symbiotic genes and the data suggested a new symbiovar, for which sv. "cyanophyllae" is proposed. Isolates formed effective nodules on A. saligna.

Long-term fertilization lowers the alkaline phosphatase activity by impacting the phoD-harboring bacterial community in rice-winter wheat rotation system

PhoD-harboring bacteria and the secreted alkaline phosphatases (ALP) are crucial in the regulation of soil phosphorus (P) cycling. However, the influential factors of these crucial indicators and their internal interactions remain controversial. Here, a long-term field experiment containing different fertilization regimes was conducted (chemical, organic, and no fertilizer applied). The results indicated that the richness and diversity of phoD-harboring bacterial community were significantly decreased after long-term fertilization. The applied fertilizer promoted the growth of competitive species, while phoD-harboring bacteria lost the advantage to outcompete other microorganisms after long-term fertilization. The decreased ALP activity was caused by the declined phoD gene abundance, which is attributed to the comprehensive effects of soil organic C (SOC), total nitrogen (TN), and various forms of P. The random forest models identified SOC, TN, and available P (AP) to be the dominant environmental factors in shaping the phoD-harboring bacterial community. In addition, some other forms of P such as organic P (Po), inorganic P (Pi) or total P (TP) also exerted significant effects. Different fertilization regimes changed the keystone genera that contributed significantly to soil ALP activities, while Pseudolabrys and Pseudomonas were predicted to be the most important genera regardless of different fertilization regimes. This study extends the understanding of the main process and mechanisms of P mobilization in response to different fertilization regimes.

The Space-Exposed Kombucha Microbial Community Member Komagataeibacter oboediens Showed Only Minor Changes in Its Genome After Reactivation on Earth

Komagataeibacter is the dominant taxon and cellulose-producing bacteria in the Kombucha Microbial Community (KMC). This is the first study to isolate the K. oboediens genome from a reactivated space-exposed KMC sample and comprehensively characterize it. The space-exposed genome was compared with the Earth-based reference genome to understand the genome stability of K. oboediens under extraterrestrial conditions during a long time. Our results suggest that the genomes of K. oboediens IMBG180 (ground sample) and K. oboediens IMBG185 (space-exposed) are remarkably similar in topology, genomic islands, transposases, prion-like proteins, and number of plasmids and CRISPR-Cas cassettes. Nonetheless, there was a difference in the length of plasmids and the location of cas genes. A small difference was observed in the number of protein coding genes. Despite these differences, they do not affect any genetic metabolic profile of the cellulose synthesis, nitrogen-fixation, hopanoid lipids biosynthesis, and stress-related pathways. Minor changes are only observed in central carbohydrate and energy metabolism pathways gene numbers or sequence completeness. Altogether, these findings suggest that K. oboediens maintains its genome stability and functionality in KMC exposed to the space environment most probably due to the protective role of the KMC biofilm. Furthermore, due to its unaffected metabolic pathways, this bacterial species may also retain some promising potential for space applications.

High Abundance of Thaumarchaeota Found in Deep Metamorphic Subsurface in Eastern China

Members of the Thaumarchaeota phylum play a key role in nitrogen cycling and are prevalent in a variety of environments including soil, sediment, and seawater. However, few studies have shown the presence of Thaumarchaeota in the terrestrial deep subsurface. Using high-throughput 16S rRNA gene sequencing, this study presents evidence for the high relative abundance of Thaumarchaeota in a biofilm sample collected from the well of Chinese Continental Scientific Drilling at a depth of 2000 m. Phylogenetic analysis showed a close relationship of these thaumarchaeotal sequences with known ammonia-oxidizing archaea (AOA) isolates, suggesting the presence of AOA in the deep metamorphic environment of eastern China which is believed to be oxic. Based on fluid geochemistry and FAProTax functional prediction, a pathway of nitrogen cycling is proposed. Firstly, heterotrophic nitrogen fixation is executed by diazotrophic bacteria coupled with methane oxidation. Then, ammonia is oxidized to nitrite by AOA, and nitrite is further oxidized to nitrate by bacteria within the phylum Nitrospirae. Denitrification and anaerobic ammonia oxidation occur slowly, leading to nitrate accumulation in the subsurface. With respect to biogeochemistry, the reaction between downward diffusing O2 and upward diffusing CH4 potentially fuels the ecosystem with a high relative abundance of Thaumarchaeota.