Geospatial variation in co-occurrence networks of nitrifying microbial guilds

Variation of N cycle guilds of the rye rhizosphere microbiome is driven by crop productivity along a tillage erosion catena

Tillage erosion poses threats to crop yields. A transition towards more sustainable agricultural practices may be advanced by harnessing ecosystem services provided by plant microbiomes. However, targeting microbiomes at the agroecosystem scale necessitates bridging the gap to microscale structures of microbiomes. We hypothesized that differences of microbial nitrogen (N) cycle guilds in the rhizosphere of rye align with a soil catena that has been formed by tillage erosion. The rhizosphere was sampled at four sites, which captured a complete tillage erosion gradient from extremely eroded to depositional soils. The gene abundances characteristic of microbial N cycle guilds were assessed via metagenomics. The eroded sites showed the lowest plant productivity and soil mineral N availability, which was associated with an enrichment of glnA in the rhizosphere. Genes associated with dissimilatory nitrate-to-ammonium reducers and diazotrophy prevailed in the eroded soil profiles. The strongest correlations of the biomasses of rye plants along the catena with N cycle functions were observed for norBC. Thus, tillage erosion as a legacy of agricultural management aligns with substantial differences in rhizosphere microbiome functionality in N cycling. These microbiome differences were linked to plant shoot properties. Thus, the dynamics of the microbiome can be indirectly assessed by remote sensing.

Biochar and peat amendments affect nitrogen retention, microbial capacity and nitrogen cycling microbial communities in a metal and polycyclic aromatic hydrocarbon contaminated urban soil

Soil contaminants may restrict soil functions. A promising soil remediation method is amendment with biochar, which has the potential to both adsorb contaminants and improve soil health. However, effects of biochar amendment on soil-plant nitrogen (N) dynamics and N cycling microbial guilds in contaminated soils are still poorly understood. Here, a metal- and polycyclic aromatic hydrocarbon (PAH) contaminated soil was amended with either biochar (0, 3, 6 % w/w) and/or peat (0, 1.5, 3 % w/w) in a full-factorial design and sown with perennial ryegrass in an outdoor field trial. After three months, N and the stable isotopic ratio δ15N was measured in soil, roots and leaves, along with microbial responses. Aboveground grass biomass decreased by 30 % and leaf N content by 20 % with biochar, while peat alone had no effect. Peat in particular, but also biochar, stimulated the abundance of microorganisms (measured as 16S rRNA gene copy number) and basal respiration. Microbial substrate utilization (MicroResp™) was altered differentially, as peat increased respiration of all carbon sources, while for biochar, respiration of carboxylic acids increased, sugars decreased, and was unaffected for amino acids. Biochar increased the abundance of ammonia oxidizing archaea, while peat stimulated ammonia oxidizing bacteria, Nitrobacter-type nitrite oxidizers and comB-type complete ammonia oxidizers. Biochar and peat also increased nitrous oxide reducing communities (nosZI and nosZII), while peat alone or combined with biochar also increased abundance of nirK-type denitrifiers. However, biochar and peat lowered leaf δ15N by 2-4 ‰, indicating that processes causing gaseous N losses, like denitrification and ammonia volatilization, were reduced compared to the untreated contaminated soil, probably an effect of biotic N immobilization. Overall, this study shows that in addition to contaminant stabilization, amendment with biochar and peat can increase N retention while improving microbial capacity to perform important soil functions.

Ammonia-oxidizing archaea adapted better to the dark, alkaline oligotrophic karst cave than their bacterial counterparts

Subsurface karst caves provide unique opportunities to study the deep biosphere, shedding light on microbial contribution to elemental cycling. Although ammonia oxidation driven by both ammonia-oxidizing bacteria (AOB) and ammonia-oxidizing archaea (AOA) is well explored in soil and marine environments, our understanding in the subsurface biosphere still remained limited to date. To address this gap, weathered rock and sediment samples were collected from the Xincuntun Cave in Guilin City, an alkaline karst cave, and subjected to high-throughput sequencing and quantification of bacterial and archaeal amoA, along with determination of the potential nitrification rates (PNR). Results revealed that AOA dominated in ammonia oxidation, contributing 48-100% to the PNR, and AOA amoA gene copies outnumbered AOB by 2 to 6 orders. Nitrososphaera dominated in AOA communities, while Nitrosopira dominated AOB communities. AOA demonstrated significantly larger niche breadth than AOB. The development of AOA communities was influenced by deterministic processes (50.71%), while AOB communities were predominantly influenced by stochastic processes. TOC, NH4+, and Cl- played crucial roles in shaping the compositions of ammonia oxidizers at the OTU level. Cross-domain co-occurrence networks highlighted the dominance of AOA nodes in the networks and positive associations between AOA and AOB, especially in the inner zone, suggesting collaborative effort to thrive in extreme environments. Their high gene copies, dominance in the interaction with ammonia oxidizing bacteria, expansive niche breadth and substantial contribution to PNR collectively confirmed that AOA better adapted to alkaline, oligotrophic karst caves environments, and thus play a fundamental role in nitrogen cycling in subsurface biosphere.

Assessing the activity of different plant-derived molecules and potential biological nitrification inhibitors on a range of soil ammonia- and nitrite-oxidizing strains

IMPORTANCE

Synthetic nitrification inhibitors are routinely used with nitrogen fertilizers to reduce nitrogen losses from agroecosystems, despite having drawbacks like poor efficiency, cost, and entry into the food chain. Plant-derived BNIs constitute a more environmentally conducive alternative. Knowledge on the activity of BNIs to soil nitrifiers is largely based on bioassays with a single Nitrosomonas europaea strain which does not constitute a dominant member of the community of ammonia-oxidizing microorganisms (AOM) in soil. We determined the activity of several plant-derived molecules reported as having activity, including the recently discovered maize-isolated BNI, zeanone, and its natural analog, 2-methoxy-1,4-naphthoquinone, on a range of ecologically relevant AOM and one nitrite-oxidizing bacterial culture, expanding our knowledge on the intrinsic inhibition potential of BNIs toward AOM and highlighting the necessity for a deeper understanding of the effect of BNIs on the overall soil microbiome integrity before their further use in agricultural settings.

Metagenomic and network analyses reveal key players in nitrification in upland agricultural soils

Nitrification, a key step in soil nitrogen cycling, is a biologically mediated process crucial to the ecological environment. However, how nitrifiers drive nitrification under different soil properties and climatic factors at large spatial scales is poorly understood. Here, using metagenomic sequencing and network-based approaches, we identified key nitrifying species of upland agricultural soils in northern China, which spans a wide range of climates and geographic distances. We found that potential nitrification rates (PNRs) varied in different soils and were positively correlated with soil pH (5.42-8.46) and mean annual temperature (MAT) and negatively correlated with the C/N ratio. Network analysis revealed that one module (module 3) was significantly correlated with PNR. In this module, 16 dominant nodes were associated with AOB Nitrosomonas and most nodes were significantly correlated with environmental factors, suggesting that abiotic conditions are important for determining the assembly of these key nitrifiers. Our study advanced the understanding of the key nitrifying populations and their environmental drivers in upland agricultural soil across different soil and climate types.

Fertilization altered co-occurrence patterns and microbial assembly process of ammonia-oxidizing microorganisms

Ammonia-oxidizing archaea and bacteria (AOA and AOB, respectively) are important intermediate links in the nitrogen cycle. Apart from the AOA and AOB communities in soil, we further investigated co-occurrence patterns and microbial assembly processes subjected to inorganic and organic fertilizer treatments for over 35 years. The amoA copy numbers and AOA and AOB communities were found to be similar for the CK and organic fertilizer treatments. Inorganic fertilizers decreased the AOA gene copy numbers by 0.75-0.93-fold and increased the AOB gene copy numbers by 1.89-3.32-fold compared to those of the CK treatment. The inorganic fertilizer increased Nitrososphaera and Nitrosospira. The predominant bacteria in organic fertilizer was Nitrosomonadales. Furthermore, the inorganic fertilizer increased the complexity of the co-occurrence pattern of AOA and decreased the complexity pattern of AOB comparing with organic fertilizer. Different fertilizer had an insignificant effect on the microbial assembly process of AOA. However, great difference exists in the AOB community assembly process: deterministic process dominated in organic fertilizer treatment and stochastic processes dominated in inorganic fertilizer treatment, respectively. Redundancy analysis indicated that the soil pH, NO₃^-N, and available phosphorus contents were the main factors affecting the changes in the AOA and AOB communities. Overall, this findings expanded our knowledge concerning AOA and AOB, and ammonia-oxidizing microorganisms were more disturbed by inorganic fertilizers than organic fertilizers.

Niche Differentiation Among Canonical Nitrifiers and N2O Reducers Is Linked to Varying Effects of Nitrification Inhibitors DCD and DMPP in Two Arable Soils

The efficacy of nitrification inhibitors (NIs) dicyandiamide (DCD) and 3,4-dimethylpyrazole phosphate (DMPP) varies with soil types. Understanding the microbial mechanisms for this variation may lead to better modelling of NI efficacy and therefore on-farm adoption. This study addressed the response patterns of mineral nitrogen, nitrous oxide (N₂O) emission, abundances of N-cycling functional guilds and soil microbiota characteristics, in relation to urea application with or without DCD or DMPP in two arable soils (an alkaline and an acid soil). The inhibition of nitrification rate and N₂O emission by NI application occurred by suppressing ammonia-oxidizing bacteria (AOB) abundances and increasing the abundances of nosZI-N₂O reducers; however, abundances of ammonia-oxidizing archaea (AOA) were also stimulated with NIs-added in these two arable soils. DMPP generally had stronger inhibition efficiency than DCD, and both NIs' addition decreased Nitrobacter, while increased Nitrospira abundance only in alkaline soil. N₂O emissions were positively correlated with AOB and negatively correlated with nosZI in both soils and AOA only in acid soil. Moreover, N₂O emissions were also positively correlated with nirK-type denitrifiers in alkaline soil, and clade A comammox in acid soil. Amendment with DCD or DMPP altered soil microbiota community structure, but had minor effect on community composition. These results highlight a crucial role of the niche differentiation among canonical ammonia oxidizers (AOA/AOB), Nitrobacter and Nitrospira, as well as nosZI- and nosZII-N₂O reducers in determining the varying efficacies of DCD and DMPP in different arable soils.

Diverse Bathyarchaeotal Lineages Dominate Archaeal Communities in the Acidic Dajiuhu Peatland, Central China

Bathyarchaeota are believed to have roles in the carbon cycle in marine systems. However, the ecological knowledge of Bathyarchaeota is limited in peatland ecosystems. Here, we investigated the vertical distribution of Bathyarchaeota community structure using quantitative PCR and high-throughput sequencing technology of ribosomal 16S rRNA gene integrated with detailed chemical profiling in the Dajiuhu Peatland, central China. Eight archaeal phyla were observed in peat samples, which mainly composed of Bathyarchaeota with a mean relative abundance about 88%, followed by Thaumarchaeota (9%). Bathyarchaeota were further split into 17 subgroups, and some subgroups showed habitat specificity to peat horizons with distinct lithological and physicochemical properties, for example, Bathy-6 and Bathy-15 had preference for the acrotelm, Bathy-5b, Bathy-16, and Bathy-19 were enriched in the catotelm, Bathy-5a, Bathy-8, and Bathy-11 were specific for the clay horizon. This spatial distribution pattern of archaeal communities along peat profile was mainly influenced by water content as indicated by RDA ordination and permutational MANOVA, whereas organic matter content exclusively affected Bathyarchaeota distribution along the peat profile significantly. The abundance of archaeal 16S rRNA genes ranged from 10⁵ to 10⁷ copies per gram dry sediment, and the highest archaeal biomass was observed in the periodically oxic mesotelm horizon with more dynamic archaeal interaction relationship as indicated by the network analysis. Bathyarchaeota dominated the archaeal interaction network with 82% nodes, 96% edges, and 71% keystone species. Our results provide an overview of the archaeal population, community structure, and relationship with environmental factors that affect the vertical distribution of archaeal communities and emphasize the ecology of bathyarchaeotal lineages in terrestrial peatland ecosystems.

The Effects of Quinone Imine, a New Potent Nitrification Inhibitor, Dicyandiamide, and Nitrapyrin on Target and Off-Target Soil Microbiota

Dicyandiamide (DCD) and nitrapyrin (NP) are nitrification inhibitors (NIs) used in agriculture for over 40 years. Recently, ethoxyquin (EQ) was proposed as a novel potential NI, acting through its derivative quinone imine (QI). Still, the specific activity of these NIs on the different groups of ammonia-oxidizing microorganisms (AOM), and mostly their effects on other soil microbiota remain unknown. We determined the impact of QI, and comparatively of DCD and NP, applied at two doses (regular versus high), on the function, diversity, and dynamics of target (AOM), functionally associated (nitrite-oxidizing bacteria-NOB), and off-target prokaryotic and fungal communities in two soils mainly differing in pH (5.4 versus 7.9). QI was equally effective to DCD but more effective than NP in inhibiting nitrification in the acidic soil, while in the alkaline soil QI was less efficient than DCD and NP. This was attributed to the higher activity of QI toward AOA prevailing in the acidic soil. All NIs induced significant effects on the composition of the AOB community in both soils, unlike AOA, which were less responsive. Beyond on-target effects, we noted an inhibitory effect of all NIs on the abundance of NOB in the alkaline soil, with Nitrobacter being more sensitive than Nitrospira. QI, unlike the other NIs, induced significant changes in the composition of the bacterial and fungal communities in both soils. Our findings have serious implications for the efficiency and future use of NIs on agriculture and provide unprecedented evidence for the potential off-target effects of NIs on soil microbiota. IMPORTANCE NIs could improve N use efficiency and decelerate N cycling. Still, we know little about their activity on the distinct AOM groups and about their effects on off-target soil microorganisms. Here, we studied the behavior of a new potent NI, QI, compared to established NIs. We show that (i) the variable efficacy of NIs across soils with different pH reflects differences in the inherent specific activity of the NIs to AOA and AOB; (ii) beyond AOM, NIs exhibit negative effects on other nitrifiers, like NOB; (iii) QI was the sole NI that significantly affected prokaryotic and fungal diversity. Our findings (i) highlight the need for novel NI strategies that consider the variable sensitivity of AOM groups to the different NIs (ii) identify QI as a potent AOA inhibitor, and (iii) stress the need for monitoring NIs' impact on off-target soil microorganisms to ensure sustainable N fertilizers use and soil ecosystem functioning.

Effects of agricultural land use on the differentiation of nitrifier communities and functional patterns from natural terrestrial ecosystems

Human activities severely affect the global nitrogen (N) cycle. Croplands receive intensive N fertilization; consequently, cropland and natural ecosystem differentiation often results in community and functional variation in N-transforming microbes, including nitrifiers, which perform nitrification central to N cycle. However, evidence of such variation is mostly limited to ammonia oxidizers (AO) in local fields, excluding soil heterogeneity and nitrite-oxidizing bacteria (NOB); the variation under diverse climatic and soil conditions is not comprehensively understood. We conducted a large-scale survey of 131 cropland and natural sites in China. The community patterns of ammonia-oxidizing bacteria (AOB) and NOB differed significantly between croplands and some natural ecosystems, whereas ammonia-oxidizing archaea (AOA) were not affected by ecosystem type. The AOB population and nitrification potential (NP) were significantly higher in agroecosystems than in natural systems except wetlands. Fewer co-occurrence interactions involving nitrifiers were observed in croplands than in natural ecosystems except forests, systematically indicating the ecological diversification of nitrifiers in potential microbial associations among these habitats. Ecosystem type, pH, organic matter (OM), total phosphorus (TP), mean annual temperature (MAT) and mean annual precipitation (MAP) were primary drivers of nitrifier community and functional shifts. This study provides the first large-scale evidence of overall nitrifier community (i.e., AOA, AOB and NOB) and potential functional shifts between agroecosystems and natural environments, enabling predictions of terrestrial N cycle under foreseeable natural land use conversions and global climate change.

From diversity to complexity: Microbial networks in soils

Network analysis has been used for many years in ecological research to analyze organismal associations, for example in food webs, plant-plant or plant-animal interactions. Although network analysis is widely applied in microbial ecology, only recently has it entered the realms of soil microbial ecology, shown by a rapid rise in studies applying co-occurrence analysis to soil microbial communities. While this application offers great potential for deeper insights into the ecological structure of soil microbial ecosystems, it also brings new challenges related to the specific characteristics of soil datasets and the type of ecological questions that can be addressed. In this Perspectives Paper we assess the challenges of applying network analysis to soil microbial ecology due to the small-scale heterogeneity of the soil environment and the nature of soil microbial datasets. We review the different approaches of network construction that are commonly applied to soil microbial datasets and discuss their features and limitations. Using a test dataset of microbial communities from two depths of a forest soil, we demonstrate how different experimental designs and network constructing algorithms affect the structure of the resulting networks, and how this in turn may influence ecological conclusions. We will also reveal how assumptions of the construction method, methods of preparing the dataset, and definitions of thresholds affect the network structure. Finally, we discuss the particular questions in soil microbial ecology that can be approached by analyzing and interpreting specific network properties. Targeting these network properties in a meaningful way will allow applying this technique not in merely descriptive, but in hypothesis-driven research. Analysing microbial networks in soils opens a window to a better understanding of the complexity of microbial communities. However, this approach is unfortunately often used to draw conclusions which are far beyond the scientific evidence it can provide, which has damaged its reputation for soil microbial analysis. In this Perspectives Paper, we would like to sharpen the view for the real potential of microbial co-occurrence analysis in soils, and at the same time raise awareness regarding its limitations and the many ways how it can be misused or misinterpreted.

Reactive nitrogen restructures and weakens microbial controls of soil N2O emissions

The global surplus of reactive nitrogen (Nr) in agricultural soils is accelerating nitrous oxide (N2O) emission rates, and may also strongly influence the microbial controls of this greenhouse gas resulting in positive feedbacks that further exacerbate N2O emissions. Yet, the link between legacy effects of Nr on microbial communities and altered regulation of N2O emissions is unclear. By examining soils with legacies of Nr-addition from 14 field experiments with different edaphic backgrounds, we show that increased potential N2O production is associated with specific phylogenetic shifts in communities of frequently occurring soil microbes. Inputs of Nr increased the complexity of microbial co-association networks, and altered the relative importance of biotic and abiotic predictors of potential N2O emissions. Our results provide a link between the microbial legacy of Nr addition and increased N2O emissions by demonstrating that biological controls of N2O emissions were more important in unfertilized soils and that these controls are weakened by increasing resource levels in soil.

Response of a methane-driven interaction network to stressor intensification

Microorganisms may reciprocally select for specific interacting partners, forming a network with interdependent relationships. The methanotrophic interaction network, comprising methanotrophs and non-methanotrophs, is thought to modulate methane oxidation and give rise to emergent properties beneficial for the methanotrophs. Therefore, microbial interaction may become relevant for community functioning under stress. However, empirical validation of the role and stressor-induced response of the interaction network remains scarce. Here, we determined the response of a complex methane-driven interaction network to a stepwise increase in NH4Cl-induced stress (0.5-4.75 g L-1, in 0.25-0.5 g L-1 increments) using enrichment of a naturally occurring complex community derived from a paddy soil in laboratory-scale incubations. Although ammonium and intermediates of ammonium oxidation are known to inhibit methane oxidation, methanotrophic activity was unexpectedly detected even in incubations with high ammonium levels, albeit rates were significantly reduced. Sequencing analysis of the 16S rRNA and pmoA genes consistently revealed divergent communities in the reference and stressed incubations. The 16S rRNA-based co-occurrence network analysis revealed that NH4Cl-induced stress intensification resulted in a less complex and modular network, likely driven by less stable interaction. Interestingly, the non-methanotrophs formed the key nodes, and appear to be relevant members of the community. Overall, stressor intensification unravels the interaction network, with adverse consequences for community functioning.

Depth Profile of Nitrifying Archaeal and Bacterial Communities in the Remote Oligotrophic Waters of the North Pacific

Nitrification is a vital ecosystem function in the open ocean that regenerates inorganic nitrogen and promotes primary production. Recent studies have shown that the ecology and physiology of nitrifying organisms is more complex than previously postulated. The distribution of these organisms in the remote oligotrophic ocean and their interactions with the physicochemical environment are relatively understudied. In this work, we aimed to evaluate the depth profile of nitrifying archaea and bacteria in the Eastern North Pacific Subtropical Front, an area with limited biological surveys but with intense trophic transferences and physicochemical gradients. Furthermore, we investigated the dominant physicochemical and biological relationships within and between ammonia-oxidizing archaea (AOA), ammonia-oxidizing bacteria (AOB), and nitrite-oxidizing bacteria (NOB) as well as with the overall prokaryotic community. We used a 16S rRNA gene sequencing approach to identify and characterize the nitrifying groups within the first 500 m of the water column and to analyze their abiotic and biotic interactions. The water column was characterized mainly by two contrasting environments, warm O₂-rich surface waters with low dissolved inorganic nitrogen (DIN) and a cold O₂-deficient mesopelagic layer with high concentrations of nitrate (NO₃^–). Thaumarcheotal AOA and bacterial NOB were highly abundant below the deep chlorophyll maximum (DCM) and in the mesopelagic. In the mesopelagic, AOA and NOB represented up to 25 and 3% of the total prokaryotic community, respectively. Interestingly, the AOA community in the mesopelagic was dominated by unclassified genera that may constitute a novel group of AOA highly adapted to the conditions observed at those depths. Several of these unclassified amplicon sequence variants (ASVs) were positively correlated with NO₃^– concentrations and negatively correlated with temperature and O₂, whereas known thaumarcheotal genera exhibited the opposite behavior. Additionally, we found a large network of positive interactions within and between putative nitrifying ASVs and other prokaryotic groups, including 13230 significant correlations and 23 sub-communities of AOA, AOB, NOB, irrespective of their taxonomic classification. This study provides new insights into our understanding of the roles that AOA may play in recycling inorganic nitrogen in the oligotrophic ocean, with potential consequences to primary production in these remote ecosystems.

Comparison and interpretation of characteristics of Rhizosphere microbiomes of three blueberry varieties

Background

Studies on the rhizosphere microbiome of various plants proved that rhizosphere microbiota carries out various vital functions and can regulate the growth and improve the yield of plants. However, the rhizosphere microbiome of commercial blueberry was only reported by a few studies and remains elusive. Comparison and interpretation of the characteristics of the rhizosphere microbiome of blueberry are critical important to maintain its health.

Results

In this study, a total of 20 rhizosphere soil samples, including 15 rhizosphere soil samples from three different blueberry varieties and five bulk soil samples, were sequenced with a high-throughput sequencing strategy. Based on these sequencing datasets, we profiled the taxonomical, functional, and phenotypic compositions of rhizosphere microbial communities for three different blueberry varieties and compared our results with a previous study focused on the rhizosphere microbiome of blueberry varieties. Our results demonstrated significant differences in alpha diversity and beta diversity of rhizosphere microbial communities of different blueberry varieties and bulk soil. The distribution patterns of taxonomical, functional, and phenotypic compositions of rhizosphere microbiome differ across the blueberry varieties. The rhizosphere microbial communities of three different blueberry varieties could be distinctly separated, and 28 discriminative biomarkers were selected to distinguish these three blueberry varieties. Core rhizosphere microbiota for blueberry was identified, and it contained 201 OTUs, which were mainly affiliated with Proteobacteria, Actinobacteria, and Acidobacteria. Moreover, the interactions between OTUs of blueberry rhizosphere microbial communities were explored by a co-occurrence network of OTUs from an ecological perspective.

Conclusions

This pilot study explored the characteristics of blueberry’s rhizosphere microbial community, such as the beneficial microorganisms and core microbiome, and provided an integrative perspective on blueberry’s rhizosphere microbiome, which beneficial to blueberry health and production.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02092-7.

nir gene-based co-occurrence patterns reveal assembly mechanisms of soil denitrifiers in response to fire

Denitrification causes nitrogen losses from terrestrial ecosystems. The magnitude of nitrogen loss depends on the prevalence of denitrifiers, which show ecological differences if they harbour nirS or nirK genes encoding nitrite reductases with the same biological function. Thus, it is relevant to understand the mechanisms of co-existence of denitrifiers, including their response to environmental filters and competition due to niche similarities. We propose a framework to analyse the co-existence of denitrifiers across multiple assemblages by using nir gene-based co-occurrence networks. We applied it in Mediterranean soils before and during 1 year after an experimental fire. Burning did not modify nir community structure, but significantly impacted co-occurrence patterns. Bacteria with the same nir co-occurred in space, and those with different nir excluded each other, reflecting niche requirements: nirS abundance responded to nitrate and salinity, whereas nirK to iron content. Prior to fire, mutual exclusion between bacteria with the same nir suggested competition due to niche similarities. Burning provoked an immediate rise in mineral nitrogen and erased the signals of competition, which emerged again within days as nir abundances peaked. nir co-occurrence patterns can help infer the assembly mechanisms of denitrifying communities, which control nitrogen losses in the face of ecological disturbance.

Co-association of Two nir Denitrifiers Under the Influence of Emergent Macrophytes

Diverse microorganisms perform similar metabolic process in biogeochemical cycles, whereas they are found of highly genomic differentiation. Biotic interactions should be considered in any community survey of these functional groups, as they contribute to community assembly and ultimately alter ecosystem properties. Current knowledge has mainly been achieved based on functional community characterized by a single gene using co-occurrence network analysis. Biotic interactions between functionally equivalent microorganisms, however, have received much less attention. Herein, we propose the nirK- and nirS-type denitrifier communities represented by these two nitrite reductase (nir)-encoding genes, as model communities to investigate the potential interactions of two nir denitrifiers. We evaluated co-occurrence patterns and co-association network structures of nir denitrifier community from an emergent macrophyte-dominated riparian zone of highly active denitrification in Lake Taihu, China. We found a more segregated pattern in combined nir communities than in individual communities. Network analyses revealed a modularized structure of associating nir denitrifiers. An increased proportion of negative associations among combined communities relative to those of individual communities indicated potential interspecific competition between nirK and nirS denitrifiers. pH and NH₄⁺-N were the most important factors driving co-occurrence and mutual exclusion between nirK and nirS denitrifiers. We also showed the topological importance of nirK denitrifiers acting as module hubs for constructing entire association networks. We revealed previously unexplored co-association relationships between nirK and nirS denitrifiers, which were previously neglected in network analyses of individual communities. Using nir denitrifier community as a model, these findings would be helpful for us to understand the biotic interactions and mechanisms underlying how functional groups co-exist in performing biogeochemical cycles.

The limited effects of carbonaceous material amendments on nitrite-oxidizing bacteria in an Alfisol

Carbonaceous materials are soil conditioners that affect nitrogen cycles. However, how carbonaceous materials influence nitrite-oxidizing bacteria (NOB) is yet unclear. In this study, we investigated the NOB community and its potential activities under different treatments (control, biochar, straw, limestone, biochar + limestone, and straw + limestone) in an Alfisol, a type of arable soil depleted in calcium carbonate but enriched in aluminum- and iron-bearing minerals. Treatments with limestone increased soil pH, and straw inputs caused an increment of available potassium (AK). Ammonia (NH₄⁺) was inversely changed under the straw and biochar + limestone amendments. None of the treatments significantly impacted the abundance of Nitrobacter (nxrA) or the potential nitrite oxidation activity (PNO). The abundance of Nitrospira (nxrB) increased in the biochar + limestone-treated samples and was significantly correlated with PNO, pH, and AK. High-throughput sequencing results showed that the α-diversity of NOB did not change in response to the treatments. The dominant Nitrobacter OTUs were affiliated within the Clusters 3, 4, 8, and 9 (a new cluster named in this study), while those of Nitrospira were in the lineage II and Namibian soil cluster 2. The limited compositional variation for Nitrobacter was explained by pH, and that for Nitrospira by pH, TN, and NH₄⁺. Among all available data in this study, the richness of Nitrospira was the most important predictor (73%) for PNO. Therefore, we assumed that the community of nitrite oxidizers (Nitrospira) could be relatively redundant in function, supported by the observation that the carbonaceous inputs did not impact either the potential activity or the α-diversity but did affect the abundance and community composition.

Soil Microbiomes With the Genetic Capacity for Atmospheric Chemosynthesis Are Widespread Across the Poles and Are Associated With Moisture, Carbon, and Nitrogen Limitation

Soil microbiomes within oligotrophic cold deserts are extraordinarily diverse. Increasingly, oligotrophic sites with low levels of phototrophic primary producers are reported, leading researchers to question their carbon and energy sources. A novel microbial carbon fixation process termed atmospheric chemosynthesis recently filled this gap as it was shown to be supporting primary production at two Eastern Antarctic deserts. Atmospheric chemosynthesis uses energy liberated from the oxidation of atmospheric hydrogen to drive the Calvin-Benson-Bassham (CBB) cycle through a new chemotrophic form of ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO), designated IE. Here, we propose that the genetic determinants of this process; RuBisCO type IE (rbcL1E) and high affinity group 1h-[NiFe]-hydrogenase (hhyL) are widespread across cold desert soils and that this process is linked to dry and nutrient-poor environments. We used quantitative PCR (qPCR) to quantify these genes in 122 soil microbiomes across the three poles; spanning the Tibetan Plateau, 10 Antarctic and three high Arctic sites. Both genes were ubiquitous, being present at variable abundances in all 122 soils examined (rbcL1E, 6.25 × 10³–1.66 × 10⁹ copies/g soil; hhyL, 6.84 × 10³–5.07 × 10⁸ copies/g soil). For the Antarctic and Arctic sites, random forest and correlation analysis against 26 measured soil physicochemical parameters revealed that rbcL1E and hhyL genes were associated with lower soil moisture, carbon and nitrogen content. While further studies are required to quantify the rates of trace gas carbon fixation and the organisms involved, we highlight the global potential of desert soil microbiomes to be supported by this new minimalistic mode of carbon fixation, particularly throughout dry oligotrophic environments, which encompass more than 35% of the Earth’s surface.

Deciphering belowground nitrifier assemblages with elevational soil sampling in a subtropical forest ecosystem (Mount Lu, China)

The elevational distribution patterns of microbial functional groups have long been attracting scientific interest. Ammonia-oxidizers (ammonia-oxidizing archaea [AOA] and bacteria [AOB]), complete ammonia oxidation (comammox) Nitrospira and nitrite-oxidizers (e.g. Nitrobacter and Nitrospira) play crucial roles in the nitrogen cycle, yet their activities and abundances in response to elevational gradients in a subtropical forest ecosystem remain unclear. Here, we investigated the distribution of potential functions and abundances of these nitrifiers in forest soils along elevational gradients on Mount Lu, China. Our results showed that AOA and Nitrospira abundance was higher than that of their counterparts. Only AOA, Nitrobacter and comammox Nitrospira abundances followed a hump-backed-model with altitude. Soil potential ammonia-oxidation activity (PAO) and nitrite-oxidation activity (PNO) ranged from 0.003 to 0.084 and 0.34 to 0.53 μg NO2−-N g−1 dry soil h−1, respectively. The biotic (AOA, Nitrobacter, Nitrospira and comammox Nitrospira abundances) and abiotic factors (soil variables) jointly affected PAO, whereas the abiotic factors were mainly responsible for PNO. Variance partitioning analysis showed that contemporary environmental disturbance is the most important driver for the biogeography of nitrifier assemblages. Overall, our findings indicate that forest soil nitrifier assemblages exhibit a biogeographic pattern largely shaped by soil chemistry along an elevational gradient.

The DNRA-Denitrification Dichotomy Differentiates Nitrogen Transformation Pathways in Mountain Lake Benthic Habitats

Effects of nitrogen (N) deposition on microbially-driven processes in oligotrophic freshwater ecosystems are poorly understood. We quantified guilds in the main N-transformation pathways in benthic habitats of 11 mountain lakes along a dissolved inorganic nitrogen gradient. The genes involved in denitrification (nirS, nirK, nosZ), nitrification (archaeal and bacterial amoA), dissimilatory nitrate reduction to ammonium (DNRA, nrfA) and anaerobic ammonium oxidation (anammox, hdh) were quantified, and the bacterial 16S rRNA gene was sequenced. The dominant pathways and associated bacterial communities defined four main N-transforming clusters that differed across habitat types. DNRA dominated in the sediments, except in the upper layers of more productive lakes where nirS denitrifiers prevailed with potential N₂O release. Loss as N₂ was more likely in lithic biofilms, as indicated by the higher hdh and nosZ abundances. Archaeal ammonia oxidisers predominated in the isoetid rhizosphere and rocky littoral sediments, suggesting nitrifying hotspots. Overall, we observed a change in potential for reactive N recycling via DNRA to N losses via denitrification as lake productivity increases in oligotrophic mountain lakes. Thus, N deposition results in a shift in genetic potential from an internal N accumulation to an atmospheric release in the respective lake systems, with increased risk for N₂O emissions from productive lakes.

Geospatial variation in co-occurrence networks of nitrifying microbial guilds

Microbial communities transform nitrogen (N) compounds, thereby regulating the availability of N in soil. The N cycle is defined by interacting microbial functional groups, as inorganic N‐products formed in one process are the substrate in one or several other processes. The nitrification pathway is often a two‐step process in which bacterial or archaeal communities oxidize ammonia to nitrite, and bacterial communities further oxidize nitrite to nitrate. Little is known about the significance of interactions between ammonia‐oxidizing bacteria (AOB) and archaea (AOA) and nitrite‐oxidizing bacterial communities (NOB) in determining the spatial variation of overall nitrifier community structure. We hypothesize that nonrandom associations exist between different AO and NOB lineages that, along with edaphic factors, shape field‐scale spatial patterns of nitrifying communities. To address this, we sequenced and quantified the abundance of AOA, AOB, and Nitrospira and Nitrobacter NOB communities across a 44‐hectare site with agricultural fields. The abundance of Nitrobacter communities was significantly associated only with AOB abundance, while that of Nitrospira was correlated to AOA. Network analysis and geostatistical modelling revealed distinct modules of co‐occurring AO and NOB groups occupying disparate areas, with each module dominated by different lineages and associated with different edaphic factors. Local communities were characterized by a high proportion of module‐connecting versus module‐hub nodes, indicating that nitrifier assemblages in these soils are shaped by fluctuating conditions. Overall, our results demonstrate the utility of network analysis in accounting for potential biotic interactions that define the niche space of nitrifying communities at scales compatible to soil management.