Original Research ARTICLE
Global Distribution of Anaerobic Ammonia Oxidation (Anammox) Bacteria – Field Surveys in Wetland, Dryland, Groundwater Aquifer and Snow
- 1Key Laboratory of Drinking Water Science and Technology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing, China
- 2Institute of Environmental Research at Greater Bay, Key Laboratory for Water Quality and Conservation of the Pearl River Delta, Guangzhou University, Guangzhou, China
The discovery of anaerobic ammonia oxidation (anammox) expanded our knowledge on the microbial nitrogen cycle. Previous studies report that anammox bacteria are distributed in a wide range of habitats and plays significant roles in the global nitrogen cycle. However, most studies focus only on individual ecosystems or datasets from public databases. To date, our understanding of how anammox bacteria respond to environmental properties and are distributed in different habitats on a global scale, remain unclear. To explore the global distribution of anammox bacteria, samples were collected from different habitats at different locations globally, including wetlands, drylands, groundwater aquifers and snow from 10 countries across six continents. We then used high-throughput amplicon sequencing targeting the functional gene hydrazine synthase (HZS) and generated community profiles. Results showed that Candidatus Brocadia is detected as the dominant genus on a global scale, accounting for 80.0% to 99.9% of the retrieved sequences in different habitats. The Jettenia-like sequences were the second most abundant group, accounting for no more than 19.9% of the retrieved sequences in all sites. The samples in drylands, wetlands and groundwater aquifers showed similar community composition and diversity, with the snow samples being the most different. Deterministic processes seem stronger in regulating the community composition of anammox bacteria, which is supported by the higher proportion explained by local-scale factors. Groundwater aquifers showed high gene abundance and the most complex co-occurrence network among the four habitat types, suggesting that it might be the preferred habitat of anammox bacteria. There is little competition between anammox bacterial species based on co-occurrence analysis. Hence, we could infer that environmental factors such as anaerobic and stable conditions, instead of substrate limitations, may be vital factors determining the anammox bacteria community. These results provide a better understanding of the global distribution of anammox bacteria and the ecological factors that affect their community structuring in diverse habitats.
The process of anammox (anaerobic ammonium oxidation), referring to the oxidation of ammonium coupled with the reduction of nitrite under anoxic conditions, leads to a loss of nitrogen in form of N2 gas (Strous et al., 1999; Jetten, 2001). The anammox bacteria is a deep-branching monophyletic group within the phylum Planctomycetes (Jetten et al., 2003). The roles of anammox bacteria in the nitrogen cycle have drawn extensive attention since it was first described in the 1990s (Mulder et al., 1995; Strous et al., 1999). A total of five genera have already been identified to date, including Candidatus Brocadia, Kuenenia, Scalindua, Anammoxoglobus, and Jettenia with 16 classified species (Sonthiphand et al., 2014). Anammox bacteria were proposed to be responsible for as much as 50% of the global removal of fixed nitrogen from the oceans (Dalsgaard et al., 2005). In freshwater lakes and agricultural soils, the anammox process was reported to account for up to 40 and 37% of the nitrogen loss, respectively (Yoshinaga et al., 2011; Zhu et al., 2011). These bacteria have also been detected in extreme environments such as hydrothermal vents and cold hydrocarbon-rich seeps (Russ et al., 2013). The distribution of anammox bacteria and their role in nitrogen loss are influenced by local environmental conditions: dissolved oxygen (Woebken et al., 2008; Lam and Kuypers, 2011), organic content (Thamdrup and Dalsgaard, 2002; Kalvelage et al., 2013), NOx– (Risgaard-Petersen et al., 2004), environmental stability (Dalsgaard et al., 2005), and salinity (Dale et al., 2009) have been described as key factors. Succeeding studies showed that anammox based processes are more environmentally friendly compared to traditional nitrogen removal processes, as they mitigate greenhouse gas emissions of N2O and CO2 (Kartal et al., 2010; Wang et al., 2018). Until now, our understanding of the ecological functions and environmental significance of anammox bacteria is continues to grow.
In a recent study, anammox bacteria were resuscitated from a >10,000 years dormant state by the addition of water (Zhu et al., 2019). This finding suggests that the presence and distribution of anammox bacteria may be directly linked to their ecological functions. In light of global climate change, where increasing precipitation is expected, the dormant anammox bacteria present in different environments may also be revived, allowing them to be actively involved in global nitrogen cycles (Zhu et al., 2019). Anammox bacteria are highly habitat-specific and their community structure is largely dependent on prevailing environmental conditions (Dalsgaard et al., 2005; Dale et al., 2009; Zhu et al., 2015). For example, the anammox bacteria Candidatus Scalindua was almost exclusively only detected in marine environments (Schmid et al., 2007; Woebken et al., 2008; Lam and Kuypers, 2011). In contrast, the diversity of anammox bacteria in terrestrial and freshwater ecosystems is much higher, and the factors influencing their community structuring and interactions are presumably more complex. Specifically, organic carbon, ammonium, and nitrite were reported to be key factors in determining anammox bacterial community structures (Hu et al., 2012; Wang et al., 2012; Hou et al., 2013). However, to date, most studies done on annamox bacteria have only focused on individual ecosystems or datasets from public databases (e.g., GenBank). On the global scale, our understanding of how anammox bacteria respond to environmental factors and their distribution in different habitats, remains unclear.
The global distribution of microorganisms was first considered to be random (O’Malley, 2007) but recent evidence suggests otherwise (Martiny et al., 2006, 2011; Nelson et al., 2016). Deterministic and stochastic processes were suggested as the two basic mechanisms influencing microbial spatial community structuring (Martiny et al., 2006). Deterministic processes suggest that species are filtered out by local-scale factors (such as disturbance and nutrients level) and only occur at environmentally suitable locations (Leibold et al., 2004; Aguilar et al., 2014), while stochastic processes are driven by ecological drift and affected by regional- scale factors (such as climate conditions and migration) (Martiny et al., 2011; Zhou et al., 2014a). In natural ecosystems, the two processes concurrently act in structuring the assembly of microbial communities, but many environmental parameters seem to be of particular importance in different systems (Zhou et al., 2014a). For anammox bacteria, how multiple environmental factors interact to determine their distribution on a global scale, remains unknown and is critical for understanding the habitat preference and for predicting their potential ecological functions in the global nitrogen cycle.
In this study, field surveys of anammox bacteria were conducted, covering different habitats and locations, but excluding the marine environment due to its expected low diversity (Schmid et al., 2007; Woebken et al., 2008; Lam and Kuypers, 2011). As one of the most wide-ranging studies on the spatial distribution of anammox bacteria to date, this work attempts to address the following questions: (i) What is the global distribution pattern of anammox bacteria? (ii) How do the habitats affect community structuring of anammox bacteria? (iii) What is the most suitable habitat for anammox bacteria and the underlying mechanism of the habitat preference?
Materials and Methods
Samples were collected from four different systems in the biosphere, which included drylands, wetlands, groundwater, and snow. Drylands and wetlands represent two of the most active and contrasting ecosystems on earth’s surface in terms of biogeochemistry. Groundwater and snow represent habitats beneath and above the earth’s surface, respectively. The sampling sites representing these different ecosystems were distributed over 10 countries and six continents around the world (Figure 1; Supplementary Tables S1, S2). The dryland ecosystem habitat types included grasslands, forestlands, agricultural farmlands, and bare land (unutilized land, covered by sparse grass or no plants in some seasons). The sediment samples from wetland ecosystems included freshwater rivers, lakes, riparian zones, paddy fields, estuaries, natural, and constructed wetlands. Fresh snow was collected at various levels of air quality according to the PM2.5 values (25 to 115 μg/m3) which refers to atmospheric particulate matter that has a diameter of less than 2.5 μm at sampling time. All the samples were collected from October 2015 to January 2016, totaling 89 samples.
Figure 1. Global maps showing the sampling sites of four different ecosystems (A) drylands, (B) wetlands, (C) groundwater aquifers, (D) snow. Colors indicate various habitats where samples were taken from, and n indicates the sample number.
The samples were processed for chemical analysis according to previous studies (Zhu et al., 2013; Zhou et al., 2014b). In brief, for each sampling site, three soil/sediment samples were collected parallel to each other from three plots (10 × 10 m). In each soil plot, samples of the soil organic layer (10 × 10 cm in area, and 0–5 cm depth) were collected at three random points using a sterile blade or auger (Ø 2 cm) and mixed together in a sterile plastic bag as a single sample. For each groundwater aquifer soil/sediment sample, three large soil pits (100 × 100 cm, 9 m apart) were dug in an approximately 100 m2 sized area using PowerProbe and GeoProbe direct push drill rigs, and three parallel samples were obtained at each sampled layer (Zhu et al., 2013, 2018a,b). Soil pH was measured after shaking with water at a soil/water ratio of 1:5 (weight/volume). Soil moisture content (MC) was measured gravimetrically by oven-drying at 105°C for 48 h (Ye et al., 2019). NH4+, NO2–, and NO3– were determined by flow injection analysis (FIA Star 5000, FOSS Tecator, Sweden) after extraction with 2 M KCl solution at a soil/KCl ratio of 1:5 (weight/volume) (Ye et al., 2019). Soil total organic matter (TOM) was measured through LOI550 (Loss on ignition at 550°C) (Dean, 1974). The content of TC, TN, and TS was measured by an Element Analyzer (Vario EL cube, Elementar, Germany) (Wang et al., 2019). For snow samples, we collected the fresh snow using a sterile plastic sheet placed on the ground when the snow began to fall. At each snow sample point (100 × 100 m), the four corners and the center position of the square were sampled as five independent replicates, which were mixed into one sample. Once collected, visible roots and residues were removed prior to homogenizing the soil/sediment fraction of each sample. Part of the soil/sediment samples were freeze-dried, and snow (after melting) samples were melted and immediately filtered through a 0.22 μm filter paper (Millipore, MA, United States) for bacterial DNA. Samples were stored in a cooler and then transported to the laboratory for further processing.
After checking by PCR amplification targeting the anammox bacterial specific gene, 74 samples (dryland = 31, wetland = 31, groundwater aquifer = 6, fresh snow = 6) tested positive for annamox bacterial gene and were then kept for downstream analyses.
DNA Extraction and PCR Amplification
Genomic DNA from soil and sediment samples were extracted from ∼ 0.33 g subsamples after freeze drying. The DNA from snow samples was extracted from 300–500 ml thawed snow water after being filtered through a 0.22 μm filter. DNA was extracted with a Fast DNA SPIN Kit for Soil (MP Biomedicals, United States), quantified and assessed for quality using a Nano Drop 2000 UV-Vis Spectrophotometer (Thermo Fisher Scientific, United States). PCR amplification of the anammox bacteria hydrazine synthase (HZS) was done using the hzsB_396F and hzsB_742R primer pairs (Wang et al., 2012). Amplification was performed in 50 μl reaction mixtures containing 5 μl 10 × buffer, 4 μl dNTP (2.5 mmol/L), 1 μl of each primer (10 mmol/L), 0.5 μl BSA, 0.25 μl Taq polymerase (2.5 U) (Takara, Japan) and 2 μl of the DNA template, and dd H2O made up to 50 μl. The amplification thermal protocol consisted of following steps: 10 min at 95°C, then 35 cycles of 60 s at 95°C, followed by 60 s at 59°C, and 45 s at 72°C, and an extension of 10 min at 72°C. The PCR products were purified in 1% agarose gel after electrophoresis using the Wizard SV Gel & PCR clean-up system (Promega, United States) before multiplexing sequencing using the Hiseq 2500 platform (Illumina, United States).
The abundance of anammox bacteria was determined by quantitative PCR using primer sets of hzsB_396F and hzsB_742R (Wang et al., 2012). Amplification was performed in 20 μl reaction mixtures: 10 μl SYBRs Premix Ex Taq (Takara, Japan), 0.4 μl ROX50, 0.2 μl of each primer (200 nmol/L), 0.4 μl bovine serum albumin (BSA) and 2 μl of the DNA template, and ddH2O added to 20 μl. The ddH2O was used as the template of the negative control. The DNA of soil samples that had been identified to have anammox activity were used as templates of the positive controls (Wang et al., 2019). The thermal cycles included 3 min at 95°C, followed by 40 cycles of 30 s at 95°C, 30 s at 59°C, and 30 s at 72°C. Standard curves were generated from 10-fold serial dilutions of the plasmid DNA. The results with an efficiency and correlation coefficient above 90% and 0.95 were utilized.
High-Throughput Sequence Processing
Raw paired-end reads of nucleic acids were assigned to samples based on their unique barcodes, which were later truncated with the primer sequences. Specifically, chimeric sequences produced during incomplete PCR amplification were deleted using the Uchime algorithm (V1.7.0) (Caporaso et al., 2010). Then low-quality sequences were removed using the QIIME. Redundant nucleic acid sequences were then determined using Mothur (Schloss et al., 2009), followed by translation into amino acid sequences using the BioEdit software (Version 7.2). Sequences demonstrating 90% similarity were assigned to the same operational taxonomic unit (OTU). Representative sequences were assigned with taxonomy using the Basic Local Alignment Search Tool (BLAST) and compared against the hzsB gene sequences from the NCBI database. Raw data of the high throughput sequencing were deposited in the NCBI Sequence Read Archive (SRA)1 under accession number SRX3003873.
Statistical and Network Analysis
The Mann–Whitney U test was used to assess significant differences in anammox bacterial α-diversity indices (i.e., Shannon diversity, Chao1 richness, and Pielou evenness) between and among different habitats (IBM SPSS Statistics 20.0 for Windows), as well as non-metric multidimensional scaling (NMDS) using Bray-Curtis dissimilarity of anammox bacterial OTUs and α-diversity. Variation partitioning analysis (VPA) was then used to quantify the relative influences of environmental factors on the anammox bacterial community in CANOCO 5 software. Analysis of similarity (ANOSIM) was calculated at the OTU level to infer the dissimilarity of the anammox bacterial community between and among habitats using the R platform with the vegan package (Version 3.0.1, permutations = 999). Spearman’s rank correlation between the abundance of anammox bacteria and environmental variables was plotted using the “corrplot” package in R (version 3.0.1).
A correlation-based network analysis was carried out to explore the co-occurrence patterns of anammox bacteria in different habitats. To visualize the associations in networks, a correlation matrix was constructed by calculating for pairwise Spearman’s rank correlations based on the OTU level. The Spearman’s rank correlation coefficient (ρ) absolute value >0.8 and the P-value < 0.05 were considered to be a statistically robust correlation between OTUs. The nodes and edges in the networks correspond to the anammox bacterial taxa (OTUs) and the strong correlations between the nodes, respectively. Topological characteristics of the networks were calculated to describe the complex pattern of interrelationships among anammox bacterial OTUs. The Spearman’s rank correlations were calculated with Python 2.6, and networks were visualized using Gephi software (Bastian et al., 2009).
Global Distribution of Anammox Bacterial Community Composition
After filtering out sequences of incongruent lengths and low quality, a total of 975,093 reads were retained for further analyses (Supplementary Table S3). Blasted results in the GenBank database showed that the majority of the processed sequences were highly similar to the anammox bacterial hzsB gene. Most of them belonged to Candidatus Brocadia, Candidatus Jettenia, Candidatus Kuenenia, and Candidatus Scalindua. Sequences that were distantly related to these genera were defined as “unknown” in a subsequent analysis.
Globally, Candidatus Brocadia was predominant in all four systems, accounting for 80.6%, 77.9%, 80.0%, and 99.9% of the retrieved sequences in drylands, wetlands, groundwater aquifers and snow habitats, respectively (Figure 2A). Within Candidatus Brocadia, Brocadia fulgida, and Brocadia sinica related sequences accounted for as much as 23.8–52.7% and 41.2–55.8% of the total sequences, which were the dominant species. The Candidatus Jettenia was the second most abundant genus which composed 18.4%, 16.8%, and 19.9% of the retrieved sequences in the dryland, wetland and groundwater aquifer systems, respectively. No Jettenia-like sequence was detected in the snow samples. The Kuenenia-like sequences only accounted for a very small fraction in each habitat (0–2.6%). Similarly, the Scalindua-like sequences were only 0–2.8% of the sequences, with the highest abundance in a wetland site.
Figure 2. Heatmap showing the relative abundance of each species or genus in each sampling site (A), the non-metric multidimensional scaling ordination (NMDS) based on Bray-Curtis similarity of anammox bacterial communities at each sampling site. Sites in arrows indicated that the dominant genus was Jettenia or Brocadia (B), variation partitioning analysis showing the relative influence of environmental factors on the structuring of the anammox bacterial community in drylands (C) and wetlands (D).
The community composition of sites in drylands and wetlands, revealed by the hzsB gene, showed considerably high variation. For instance, the majority of sites in both drylands and wetlands clustered together (Figure 2B), and were dominated by Candidatus Brocadia. Another three dryland sites and three wetland sites clustered together, and were dominated by Candidatus Jettenia. In another cluster, which included three wetland sites and two dryland sites, the Brocadia- and Jettenia-like sequences were of comparable proportions. In contrast, a relatively low variation in the community composition was observed on the snow samples sites. Notably, one sample from a riparian sediment of a city river was dominated by Scalindua-like sequences which are frequently detected in marine environments in previous studies.
The role of regional-scale factors (e.g., average annual temperature, precipitation, and sunshine duration) and local-scale factors (e.g., soil/sediment physicochemical properties) in the structuring of the anammox community were quantified by a variation partition analysis (VPA). In drylands, the two type of factors explained 26.7% of the observed variation, while local and regional factors independently explained 21.7% and 5.0% of the variation, respectively (Figure 2C). In wetlands, the local factors explained more variation than the regional factors (14.5% versus 12.7%), leaving 72.8% of the variance unexplained. But effects of the local factors were not as strong as in the drylands (14.5% versus 21.7%) (Figure 2D). Samples in groundwater aquifers and snow were excluded in the VPA due to their relatively low sample numbers.
Abundance of Anammox Bacteria
The gene abundance of anammox bacteria in wetlands revealed by the hzsB gene ranged from 5.42 × 104 to 9.56 × 106 copies/g of dry soil (Figure 3A). Among these habitats, sediments from freshwater rivers showed a low abundance (6.9 × 104 to 8.9 × 105 copies/g of dry soil). In contrast, freshwater lakes showed a higher abundance at 5.4 × 104 to 9.6 × 106 copies/g of dry soil. It was noted that the abundance of anammox bacteria in groundwater aquifers was high, ranging from 5.93 × 105 to 9.12 × 106 copies/g of dry soil. A correlation analysis suggested that the chemical conditions had lower to no effect on the abundance of these samples and that only soil moisture content (P = 0.035) and total carbon (P = 0.048) showed significant correlations with their abundance (Figure 3B). In comparison, the samples from drylands and snow failed to be quantified by the qPCR assay due to their extremely low abundance.
Figure 3. Boxplots showing the abundance of the hzsB gene in each wetland habitat and groundwater aquifer. (A) and a correlation analysis of anammox bacterial abundance and chemical parameters (B). The dryland and snow samples failed to be quantified by qPCR assay due to their extremely low abundance. In this figure mean values are indicated by the open square, the central horizontal bars are the medians, and the box lower and upper limits are the 1st and 3rd quartiles, respectively. Crosses are minimum/maximum values. Sample numbers are indicated after the habitat name.
Alpha-Diversity of Anammox Bacteria
Shannon diversity, Chao1 richness, and Pielou evenness were calculated to determine α-diversity according to the retrieved hzsB sequences. Drylands, wetlands, and groundwater aquifers showed comparable levels of diversity based on these indices. The Shannon index in these three types of habitats were 0.13–3.26, 0.61–2.59, and 0.60–2.03, the Chao1 indices were 44–186, 33–191, and 95–120, respectively (Figure 4). The Pielou indices were 0.04–0.64, 0.15–0.52, and 0.14–0.45 in drylands, wetlands, and groundwater aquifers, respectively. In contrast to their extremely low abundance, the alpha-diversity in drylands was high. However, snow had the lowest diversity among the samples (P < 0.01, one-way ANOVA). Such observations were confirmed by the NMDS, where samples in drylands, wetlands, and groundwater aquifers clustered closer together, while the snow samples were farther distributed.
Figure 4. Alpha-diversity represented by Shannon, Chao1 and Pielou indices in the four ecosystems. The P-values indicate differences between groups based on the Mann–Whitney U test (A), and NMDS analysis of α-diversity based on Euclidean distances of the four diversity indices (Shannon, Simpson, Pielou and Chao1) of the anammox bacterial community (B).
Co-occurrence Patterns in Anammox Bacterial Community
The co-occurrence patterns of anammox bacteria were explored using a network analysis based on strong and significant correlations. The number of positive correlations was much higher than the negative correlations in all four networks (ρ > 0.8, P < 0.05) (Figure 5). Drylands and wetlands showed similar co-occurrence patterns with comparable numbers of positive correlation, negative correlation, and edges. However, the ecological network in groundwater aquifers was markedly different in terms of the number of nodes, edges, and significant correlations, which were considerably higher than the other three ecosystems. Furthermore, the average path length (APL), indicating the average network distance between all pairs of nodes, showed the highest value in groundwater aquifers, while the network in the snow sample was the simplest with the lowest number of nodes, edges, and significant correlations. The average clustering coefficient (avgCC), which showed that the nodes tend to cluster together, ranged from 0.113 to 0.324 in the four networks, and the groundwater aquifers and snow showed relative higher values. A similar trend was also observed in the graph density, with higher values in groundwater aquifers and snow.
Figure 5. Ecological networks of co-occurring anammox bacterial OTUs in drylands (A), wetlands (B), groundwater aquifers (C) and snow (D), based on correlations. A connection corresponds to a strong (Spearman’s ρ > 0.8) and significant (P < 0.05) correlation. The size of each node is proportional to the number of connections (i.e., degree); the thickness of each connection between two nodes (i.e., edge) is proportional to the value of Spearman’s correlation coefficients. A blue edge indicates a positive interaction between two individual nodes, while a red edge indicates a negative interaction.
This study revealed that anammox bacteria were widespread in different ecosystems, covering a wide range of habitats and geographical regions. To our knowledge, this is the largest field investigation in terms of scale on anammox bacteria, focusing on their distribution and community structure, to date. Applying similar procedures of sampling and molecular analysis in one study allowed us to avoid experimental biases to compare results.
Globally, the Candidatus Brocadia was the most dominant genera of anammox bacteria in natural ecosystems, followed by Candidatus Jettenia. It is consistent with previous studies showing that the Brocadia, Kuenenia, and Jettenia are the most common representatives in terrestrial and freshwater systems (Humbert et al., 2010; Zhu et al., 2011; Wang et al., 2012; Bai et al., 2015). The dominance of Candidatus Brocadia is mostly associated with their versatile metabolism (e.g., multiple physiological capabilities and environmental tolerance), allowing them to be more adaptable to various environments (Gori et al., 2011).
On the global scale, local-scale factors associated with soil/sediment characteristics were stronger in explaining the variation in the anammox bacterial community. In contrast, regional-scale factors such as temperature, precipitation, and the duration of sunlight plays a lesser role, implying a stronger role of the deterministic process. Similar conclusions, that the global pattern of bacterial community structures is related to local factors such as the nutrient level (Demir-Hilton et al., 2011) and soil moisture (DeBruyn et al., 2011) have been suggested. This is distinct from the global distribution of plants, in which temperature and precipitation are the key factors in shaping a community structure (Takenaka, 2005). This indicates that the growth of plants and microorganisms are generally limited by different factors. For the bacteria of the anammox group, the micro-environment in the habitat is more important.
Local factors explain more variation in the drylands compared to the wetlands. It is probably because soil is more heterogeneous in terms of moisture content, nutrient level and other factors, which could diversify the species’ coexistence and composition (Reynolds et al., 1997). In aquatic systems, however, water might potentially sweep the local environmental gradients, decreasing the effects of local factors (Xiong et al., 2017). Moreover, the presence of water in wetlands satisfied the basic demands of the anammox bacteria, either by creating an anaerobic environment or by carrying the substrates for anammox activity (Zhou et al., 2014b). Therefore, other local factors related to soil properties in wetlands is not as important as in drylands.
Anammox bacteria was also unexpectedly detected in snow samples, with a distinct community composition and diversity. This suggests that anammox bacteria is a component of an airborne bacterial community which could come from biological aerosols like other microorganisms (Morris et al., 2011). However, it remains inconclusive that the atmosphere represented by snow was another habitat for anammox bacteria, in addition to their well-known habitats (i.e., drylands, wetlands, groundwater aquifers and the ocean), since anammox bacteria have strict requirements for a stable and anaerobic environment (Strous et al., 1998; Jetten, 2001; Dalsgaard et al., 2005; Wenk et al., 2013). Moreover, the doubling time of anammox bacteria in the best growth conditions needs at least 10–12 days (van der Star et al., 2007). However, the life span of atmospheric particulates is less than 6 weeks, even for fine particles (PM2.5) (Steinfeld, 1998), which is neither long nor stable enough to support the proliferation of anammox bacteria. Hence, anammox bacteria may merely attach to soil particles and disperse with air flow.
The long-range dispersal of microorganisms including pathogens, algae, and bacteria along with soil particles have been long observed (Griffin et al., 2001). During transport, dust particles can serve as vessels for their global dispersion (Griffin et al., 2003). It has been estimated approximately 2 × 109 metric tons of soil or sediment (dried river and lakebeds) is transported in the planet’s atmosphere each year (Perkins, 2001). In this scenario, anammox bacteria could disperse over very long distances, up to 1000s km (Steinfeld, 1998), by attaching to atmospheric particulates. Transport via atmospheric particulates is probably one of the many paths facilitating the dispersal and the worldwide distribution of anammox bacteria. The presence of anammox in aerial particles could somewhat explain why anammox bacteria could be detected in some environments using molecular approaches but tend to be inactive by isotope probing (Zhu et al., 2015).
Among the four types of habitats, anammox bacteria preferred groundwater aquifers more, which was evidenced by two aspects: a high abundance and the most complex community network. The samples in groundwater aquifers showed high gene abundance compared to those in the drylands, snow and some sites in wetlands. In the complex bacterial communities, abundance was a key indicator of their functional role (Rivett and Bell, 2018). The high anammox bacterial abundance observed here, may therefore be indicative of high anammox rates, which aligns with previous work where anammox was shown to be the dominant nitrogen loss pathway in unconfined aquifer soils (Wang et al., 2017). It has been identified that anammox is the dominant nitrogen loss pathway in unconfined aquifer soils (Wang et al., 2017). On the other hand, co-occurrence patterns in the groundwater aquifer showed the highest number of nodes, edges, positive correlations and average degree. These topological properties suggest that anammox bacteria in the groundwater aquifer had a complex and closely related community structure, which would increase the stability of interaction networks and enhance the resistance of the community to environmental disturbance (Krause et al., 2003).
A co-occurrence network analysis, which usefully represents various biological interactions in an ecosystem (Layeghifard et al., 2017), was conducted in this study. Notably, positive correlations dominated negative correlations in all four habitat types. In a network, positive interactions suggest that species are associated by commensalism or mutualisms, while negative interactions lead to predation and competition relationships (Chow et al., 2013). In this study, the low number of negative correlations implied that the competition between different species of anammox bacteria was weak in the four ecosystems. Inter-species competition, which is a major factor structuring species’ communities, is often emphasized as facing limited resources (Schoener, 1983; Vitousek and Howarth, 1991). Therefore, the low level of competition suggested that the presence and structure of the anammox community was probably not regulated by the availability of nutrient resources.
Although nutrient concentrations remained elevated at most of our sampling sites, which were greater than 0.17 mg.kg–1 (25.51 μM) for NH4+ and 0.02 mg.kg–1 (1.17 μM) for NO2– (Supplementary Tables S1, S2), we cannot rule out the possibility of nutrient competition between anammox and competitors in some micro-habitats. Anammox bacteria were identified to have high substrate affinities, with half-saturation constants (Ks) as low as 7 μM for NH4+ (Jetten et al., 2005) and 5 μM (Brocadia) and 0.2–0.3 μM (Kuenenia) for NO2– (Kartal et al., 2012). Whereas, the Ks (NH4+/NH3) of ammonium oxidizing bacteria were higher than 10 μM (Nitrosomonas) and the Ks (NO2–) of nitrite-oxidizing bacteria of Nitrospira was 10 μM and 21–135 μM of Nitrobacter (Blackburne et al., 2007). Therefore, anammox bacteria would be a stronger competitor for NH4+ and NO2– compared to most of their competitors, with substrate availability not influencing anammox bacterial community. Instead, environmental stress like oxygen exposure and habitat stability are more likely to be key factors.
It is well acknowledged that groundwater is the world’s largest freshwater resource, and it is also an extremely stable habitat. The relatively low influx in groundwater systems with a residence time of about 1,400 years, which is 70 times longer than that of surface waters (rivers, lakes, or wetlands) (Aeschbach-Hertig and Gleeson, 2012), suggests that the groundwater habitat is generally more stable compared to surface water bodies. Therefore, in addition to sufficient substrates for the anammox process (mostly NO3–), the stable environment, which is essential for anammox bacteria (Dalsgaard et al., 2005), would greatly facilitate its growth in groundwater. This could explain why the preferred habitat for anammox bacteria was the groundwater aquifer.
Anammox bacteria were widely detected in the biosphere, including different habitats on earth’s surface (e.g., grasslands, forestlands, agriculture drylands, freshwater rivers, lakes, paddy fields, and wetlands), as well as beneath and above the earth’s surface (e.g., groundwater aquifers and the snow). Globally, the Candidatus Brocadia dominates the community, followed by Candidatus Jettenia. The habitat systems of drylands, wetlands, and groundwater aquifers showed similar patterns of community composition and diversity of anammox bacteria. However, the snow samples were significantly different from the other three systems. Local-scale factors were stronger in explaining the variation of community composition, implying that the global distribution of anammox bacteria were more likely to follow a deterministic process. Groundwater was identified to be a preferred habitat for anammox bacteria based on their high gene abundance and the complex co-occurrence network. Low competition was observed between species, implying that substrate availability (NH4+ and NO2–) is not a regulator for the anammox community, instead, environmental pressures such as oxygen exposure and unstable habitat conditions may play more vital roles.
Data Availability Statement
The datasets generated for this study can be found in NCBI Sequence Read Archive, SRX3003873.
GZ and YW conceived the research. LX and SW performed the experiments. YW wrote the manuscript. GZ, YW, LX, SW, and FY edited the manuscript. LX and FY contributed to the data analysis and graphics pipeline. All authors reviewed and accepted the manuscript.
This work was supported by the National Natural Science Foundation of China (41671471, 41322012, and 91851204), the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB15020303), the National Key R&D Program (2016YFA0602303), the Local Innovative and Research Teams Project of Guangdong Pearl River Talents Program (2017BT01Z176), special funds from the State Key Joint Laboratory of Environment Simulation and Pollution Control (Research Center for Eco-environmental Sciences, Chinese Academy of Sciences) (18Z02ESPCR), and the Program of the Youth Innovation Promotion Association (CAS).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors thank Dr. Yuanyang Huang and Prof. Maohua Ma from the Chongqing Institute of Green and Intelligent Technology, Chinese Academy of Sciences for their constructive comments on this article.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.02583/full#supplementary-material
Aguilar, M., Fiore-Donno, A. M., Lado, C., and Cavalier-Smith, T. (2014). Using environmental niche models to test the ‘everything is everywhere’ hypothesis for Badhamia. ISME J. 8, 737–745. doi: 10.1038/ismej.2013.183
Bai, R., Chen, X., He, J., Shen, J., and Zhang, L. (2015). Candidatus brocadia and Candidatus kuenenia predominated in anammox bacterial community in selected Chinese paddy soils. J. Soils Sediments 15, 1977–1986. doi: 10.1007/s11368-015-1131-5
Bastian, M., Heymann, S., and Jacomy, M. (2009). “Gephi: an open source software for exploring and manipulating networks,” in Proceedings of the International Aaai Conference On Weblogs And Social Media; Third International Aaai Conference On Weblogs And Social Media, San Jose, CA.
Blackburne, R., Vadivelu, V. M., Yuan, Z., and Keller, J. (2007). Kinetic characterisation of an enriched Nitrospira culture with comparison to Nitrobacter. Water Res. 41, 3033–3042. doi: 10.1016/j.watres.2007.01.043
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336.
Chow, C.-E. T., Kim, D. Y., Sachdeva, R., Caron, D. A., and Fuhrman, J. A. (2013). Top-down controls on bacterial community structure: microbial network analysis of bacteria, T4-like viruses and protists. ISME J. 8, 816–829. doi: 10.1038/ismej.2013.199
Dale, O. R., Tobias, C. R., and Song, B. (2009). Biogeographical distribution of diverse anaerobic ammonium oxidizing (anammox) bacteria in cape fear River Estuary. Environ. Microbiol. 11, 1194–1207. doi: 10.1111/j.1462-2920.2008.01850.x
DeBruyn, J. M., Nixon, L. T., Fawaz, M. N., Johnson, A. M., and Radosevich, M. (2011). Global biogeography and quantitative seasonal dynamics of gemmatimonadetes in soil. Appl. Environ. Microbiol. 77, 6295–6300. doi: 10.1128/AEM.05005-11
Demir-Hilton, E., Sudek, S., Cuvelier, M. L., Gentemann, C. L., Zehr, J. P., and Worden, A. Z. (2011). Global distribution patterns of distinct clades of the photosynthetic picoeukaryote Ostreococcus. ISME J. 5, 1095–1107. doi: 10.1038/ismej.2010.209
Gori, F., Tringe, S. G., Kartal, B., Machiori, E., and Jetten, M. S. M. (2011). The metagenomic basis of anammox metabolism in Candidatus ‘Brocadia fulgida’. Biochem. Soc. Trans. 39, 1799–1804. doi: 10.1042/BST20110707
Griffin, D. W., Kellogg, C. A., Garrison, V. H., Lisle, J. T., Borden, T. C., and Shinn, E. A. (2003). Atmospheric microbiology in the northern Caribbean during African dust events. Aerobiologia 19, 143–157. doi: 10.1023/b:aero.0000006530.32845.8d
Hou, L., Zheng, Y., Liu, M., Gong, J., Zhang, X., Yin, G., et al. (2013). Anaerobic ammonium oxidation (anammox) bacterial diversity, abundance, and activity in marsh sediments of the Yangtze Estuary. J. Geophys. Res. 118, 1237–1246. doi: 10.1002/jgrg.20108
Hu, B. L., Shen, L. D., Zheng, P., Hu, A. H., Chen, T. T., Cai, C., et al. (2012). Distribution and diversity of anaerobic ammonium-oxidizing bacteria in the sediments of the Qiantang River. Environ. Microbiol. Rep. 4, 540–547. doi: 10.1111/j.1758-2229.2012.00360.x
Humbert, S., Tarnawski, S., Fromin, N., Mallet, M. P., Aragno, M., and Zopfi, J. (2010). Molecular detection of anammox bacteria in terrestrial ecosystems: distribution and diversity. ISME J. 4, 450–454. doi: 10.1038/ismej.2009.125
Jetten, M. S. M., Cirpus, I., Kartal, B., Niftrik, L. V., Passchoonen, K. T. V. D., Sliekers, O., et al. (2005). 1994–2004: 10 years of research on the anaerobic oxidation of ammonium. Biochem. Soc. Trans. 33, 119–123. doi: 10.1042/bst0330119
Jetten, M. S. M., Sliekers, O., Kuypers, M., Dalsgaard, T., Niftrik, L. V., Cirpus, I., et al. (2003). Anaerobic ammonium oxidation by marine and freshwater planctomycete-like bacteria. Appl. Microbiol. Biotechnol. 63, 107–114. doi: 10.1007/s00253-003-1422-4
Kalvelage, T., Lavik, G., Lam, P., Contreras, S., Arteaga, L., Löscher, C. R., et al. (2013). Nitrogen cycling driven by organic matter export in the South Pacific oxygen minimum zone. Nat. Geosci. 6, 228–234. doi: 10.1038/ngeo1739
Kartal, B., Van Niftrik, L., Keltjens, J. T., Op Den Camp, H. J. M., and Jetten, M. S. M. (2012). “Chapter 3 - anammox—growth physiology, cell biology, and metabolism,” in Advances in Microbial Physiology, ed. R. K. Poole (Cambridge, MA: Academic Press), 211–262. doi: 10.1016/b978-0-12-398264-3.00003-6
Leibold, M. A., Holyoak, M., Mouquet, N., Amarasekare, P., Chase, J. M., Hoopes, M. F., et al. (2004). The metacommunity concept: a framework for multi-scale community ecology. Ecol. Lett. 7, 601–613. doi: 10.1111/j.1461-0248.2004.00608.x
Martiny, J. B. H., Bohannan, B. J. M., Brown, J. H., Colwell, R. K., Fuhrman, J. A., Green, J. L., et al. (2006). Microbial biogeography: putting microorganisms on the map. Nat. Rev. Microbiol. 4, 102–112. doi: 10.1038/nrmicro1341
Martiny, J. B. H., Eisen, J. A., Penn, K., Allison, S. D., and Horner-Devine, M. C. (2011). Drivers of bacterial beta-diversity depend on spatial scale. Proc. Natl. Acad. Sci. U.S.A. 108, 7850–7854. doi: 10.1073/pnas.1016308108
Morris, C. E., Sands, D. C., Bardin, M., Jaenicke, R., Vogel, B., Leyronas, C., et al. (2011). Microbiology and atmospheric processes: research challenges concerning the impact of airborne micro-organisms on the atmosphere and climate. Biogeosciences 8, 17–25. doi: 10.5194/bg-8-17-2011
Mulder, A., Graaf, A. A. V. D., Robertson, L. A., and Kuenen, J. G. (1995). Anaerobic ammonium oxidation discovered in a denitrifying fluidized bed reactor. FEMS Microbiol. Ecol. 16, 177–184. doi: 10.1111/j.1574-6941.1995.tb00281.x
Nelson, M. B., Martiny, A. C., and Martiny, J. B. H. (2016). Global biogeography of microbial nitrogen-cycling traits in soil. Proc. Natl. Acad. Sci. U.S.A. 113, 8033–8040. doi: 10.1073/pnas.1601070113
Risgaard-Petersen, N., Meyer, R. C., Schmid, M., Jetten, M., Enrich-Prast, A., Rysgaard, S., et al. (2004). Anaerobic ammonia oxidation an estuarine sediment. Aqua. Microb. Ecol. 36, 293–304. doi: 10.3354/ame036293
Russ, L., Kartal, B., Hj, O. D. C., Sollai, M., Le, B. J., Caprais, J. C., et al. (2013). Presence and diversity of anammox bacteria in cold hydrocarbon-rich seeps and hydrothermal vent sediments of the Guaymas Basin. Front. Microbiol. 4:219. doi: 10.3389/fmicb.2013.00219
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Schmid, M. C., Nils, R. P., Jack, V. D. V., Kuypers, M. M. M., Gaute, L., Jan, P., et al. (2007). Anaerobic ammonium-oxidizing bacteria in marine environments: widespread occurrence but low diversity. Environ. Microbiol. 9, 1476–1484. doi: 10.1111/j.1462-2920.2007.01266.x
Strous, M., Heijnen, J. J., Kuenen, J. G., and Jetten, M. S. M. (1998). The sequencing batch reactor as a powerful tool for the study of slowly growing anaerobic ammonium-oxidizing microorganisms. Appl. Microbiol. Biotechnol. 50, 589–596. doi: 10.1007/s002530051340
Takenaka, A. (2005). Local coexistence of tree species and the dynamics of global distribution pattern along an environmental gradient: a simulation study. Ecol. Res. 20, 297–304. doi: 10.1007/s11284-005-0045-x
Thamdrup, B., and Dalsgaard, T. (2002). Production of N2 through anaerobic ammonium oxidation coupled to nitrate reduction in marine sediments. Appl. Environ. Microbiol. 68, 1312–1318. doi: 10.1128/aem.68.3.1312-1318.2002
van der Star, W. R. L., Abma, W. R., Blommers, D., Mulder, J.-W., Tokutomi, T., Strous, M., et al. (2007). Startup of reactors for anoxic ammonium oxidation: experiences from the first full-scale anammox reactor in Rotterdam. Water Res. 41, 4149–4163. doi: 10.1016/j.watres.2007.03.044
Wang, S., Wang, W., Liu, L., Zhuang, L., Zhao, S., Su, Y., et al. (2018). Microbial nitrogen cycle hotspots in the plant-bed/ditch system of a constructed wetland with N2O mitigation. Environ. Sci. Technol. 52, 6226–6236. doi: 10.1021/acs.est.7b04925
Wang, S., Wang, W., Zhao, S., Wang, X., Hefting, M. M., Schwark, L., et al. (2019). Anammox and denitrification separately dominate microbial N-loss in water saturated and unsaturated soils horizons of riparian zones. Water Res. 162, 139–150. doi: 10.1016/j.watres.2019.06.052
Wang, Y., Zhu, G., Harhangi, H. R., Zhu, B., Jetten, M. S. M., Yin, C., et al. (2012). Co-occurrence and distribution of nitrite-dependent anaerobic ammonium and methane-oxidizing bacteria in a paddy soil. FEMS Microbiol. Lett. 336, 79–88. doi: 10.1111/j.1574-6968.2012.02654.x
Wenk, C. B., Blees, J., Zopfi, J., Veronesi, M., Bourbonnais, A., Schubert, C. J., et al. (2013). Anaerobic ammonium oxidation (anammox) bacteria and sulfide-dependent denitrifiers coexist in the water column of a meromictic south-alpine lake. Limnol. Oceanogr. 58, 1–12. doi: 10.4319/lo.2013.58.1.0001
Woebken, D., Lam, P., Kuypers, M. M. M., Naqvi, S. W. A., Kartal, B., Strous, M., et al. (2008). A microdiversity study of anammox bacteria reveals a novel Candidatus scalindua phylotype in marine oxygen minimum zones. Environ. Microbiol. 10, 3106–3119. doi: 10.1111/j.1462-2920.2008.01640.x
Xiong, W., Ni, P., Chen, Y. Y., Gao, Y. C., Shan, B. Q., and Zhan, A. B. (2017). Zooplankton community structure along a pollution gradient at fine geographical scales in river ecosystems: the importance of species sorting over dispersal. Mol. Ecol. 26, 4351–4360. doi: 10.1111/mec.14199
Ye, F., Ma, M. H., Wu, S. J., Jiang, Y., Zhu, G. B., Zhang, H., et al. (2019). Soil properties and distribution in the riparian zone: the effects of fluctuations in water and anthropogenic disturbances. Eur. J. Soil Sci. 70, 664–673. doi: 10.1111/ejss.12756
Yoshinaga, I., Amano, T., Yamagishi, T., Okada, K., Ueda, S., Sako, Y., et al. (2011). Distribution and diversity of anaerobic ammonium oxidation (anammox) bacteria in the sediment of a eutrophic freshwater lake, Lake Kitaura, Japan. Microbes Environ. 26, 189–197. doi: 10.1264/jsme2.me10184
Zhou, J. Z., Deng, Y., Zhang, P., Xue, K., Liang, Y. T., Van Nostrand, J. D., et al. (2014a). Stochasticity, succession, and environmental perturbations in a fluidic ecosystem. Proc. Natl. Acad. Sci. U.S.A. 111, E836–E845. doi: 10.1073/pnas.1324044111
Zhou, L., Wang, Y., Long, X.-E., Guo, J., and Zhu, G. (2014b). High abundance and diversity of nitrite-dependent anaerobic methane-oxidizing bacteria in a paddy field profile. FEMS Microbiol. Lett. 360, 33–41. doi: 10.1111/1574-6968.12567
Zhu, G., Wang, M., Li, Y., Wang, Y., Fei, H., and Wang, S. (2018a). Denitrifying anaerobic methane oxidizing in global upland soil: sporadic and non-continuous distribution with low influence. Soil Biol. Biochem. 119, 90–100. doi: 10.1016/j.soilbio.2018.01.011
Zhu, G., Wang, S., Zhou, L., Wang, Y., Zhao, S., Xia, C., et al. (2015). Ubiquitous anaerobic ammonium oxidation in inland waters of China: an overlooked nitrous oxide mitigation process. Sci. Rep. 5:17306. doi: 10.1038/srep17306
Keywords: anammox, nitrogen cycle, co-occurrence network, groundwater aquifer, global-scale, local-scale
Citation: Wang Y, Xu L, Wang S, Ye F and Zhu G (2019) Global Distribution of Anaerobic Ammonia Oxidation (Anammox) Bacteria – Field Surveys in Wetland, Dryland, Groundwater Aquifer and Snow. Front. Microbiol. 10:2583. doi: 10.3389/fmicb.2019.02583
Received: 01 May 2019; Accepted: 24 October 2019;
Published: 12 November 2019.
Edited by:Jennifer Glass, Georgia Institute of Technology, United States
Reviewed by:Laura Anne Bristow, University of Southern Denmark, Denmark
John Senko, The University of Akron, United States
Copyright © 2019 Wang, Xu, Wang, Ye and Zhu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Guibing Zhu, email@example.com