- 1School of Ocean Science and Engineering, The University of Southern Mississippi, Ocean Springs, MS, United States
- 2Department of Earth and Environmental Sciences, Wright State University, Dayton, OH, United States
- 3Section for Aquaculture, The North Sea Research Centre, DTU Aqua, Technical University of Denmark, Hirtshals, Denmark
Nitrification is an important biological link between oxidized and reduced forms of nitrogen (N). The efficiency of nitrification plays a key role in mitigating excess N in eutrophic systems, including those with cyanobacterial harmful algal blooms (cyanoHABs), since it can be closely coupled with denitrification and removal of excess N. Recent work suggests that competition for ammonium (NH4+) between ammonia oxidizers and cyanoHABs can help determine microbial community structure. Nitrification rates and ammonia-oxidizing archaeal (AOA) and bacterial (AOB) community composition and gene abundances were quantified in Lake Okeechobee and St. Lucie Estuary in southern Florida (United States). We sampled during cyanobacterial (Microcystis) blooms in July 2016 and August 2017 (2 weeks before Hurricane Irma) and 10 days after Hurricane Irma made landfall. Nitrification rates were low during cyanobacterial blooms in Lake Okeechobee and St. Lucie Estuary, while low bloom conditions in St. Lucie Estuary coincided with greater nitrification rates. Nitrification rates in the lake were correlated (R2 = 0.94; p = 0.006) with AOA amoA abundance. Following the hurricane, nitrification rates increased by an order of magnitude, suggesting that nitrifiers outcompeted cyanobacteria for NH4+ under turbid, poor light conditions. After Irma, AOA and AOB abundances increased in St. Lucie Estuary, while only AOB increased in Lake Okeechobee. AOA sequences clustered into three major lineages: Nitrosopumilales (NP), Nitrososphaerales (NS), and Nitrosotaleales (NT). Many of the lake OTUs placed within the uncultured and uncharacterized NS δ and NT β clades, suggesting that these taxa are ecologically important along this eutrophic, lacustrine to estuarine continuum. After the hurricane, the AOA community shifted toward dominance by freshwater clades in St. Lucie Estuary and terrestrial genera in Lake Okeechobee, likely due to high rainfall and subsequent increased turbidity and freshwater loading from the lake into the estuary. AOB community structure was not affected by the disturbance. AOA communities were consistently more diverse than AOB, despite fewer sequences recovered, including new, unclassified, eutrophic ecotypes, suggesting a wider ecological biogeography than the oligotrophic niche originally posited. These results and other recent reports contradict the early hypothesis that AOB dominate ammonia oxidation in high-nutrient or terrestrial-influenced systems.
Introduction
Canonical nitrification is a two-step process in which ammonium (NH4+) is oxidized to nitrite (NO2–), and NO2– is further oxidized to nitrate (NO3–). Nitrification is a key pathway in the nitrogen (N) cycle as the sole biological link between oxidized and chemically reduced N forms, and it is often closely coupled to removal of excess N via denitrification (An and Joye, 2001). Thus, nitrification is a critical link in the ability of eutrophic ecosystems to mitigate excess N loads via coupling with denitrification.
The first step of nitrification, ammonia oxidation, is performed by ammonia-oxidizing archaea (AOA) and bacteria (AOB). Although AOA and AOB coexist in the environment, they generally exhibit niche separation (Prosser and Nicol, 2012; Hink et al., 2018). Previous studies on nitrifier community structure have posited that trophic status of the ecosystem defines the dynamics between AOA and AOB, based mostly on differences in substrate affinity for NH4+ and half-saturation constants (Km; Martens-Habbena et al., 2009; Bristow et al., 2016). Traditionally, AOB were thought to be more abundant in nutrient rich waters and soils (Jia and Conrad, 2009; Verhamme et al., 2011; Hou et al., 2013), while AOA were assumed to be more abundant in the oligotrophic open ocean (Francis et al., 2005; Könneke et al., 2005; Newell et al., 2011; Beman et al., 2012). However, a comprehensive phylogenetic analysis of >33,000 archaeal amoA sequences showed that AOA are ubiquitous, widespread, and highly diverse, and current knowledge on the physiology of cultured AOA does not represent predominant clades in the environment (Alves et al., 2018). Recently, various AOA clades from Nitrososphaerales, Nitrosotaleales, and Nitrosopumilales orders have shown to be highly abundant in shallow, eutrophic waters affected by agriculture runoff and cyanobacterial blooms (Zeng et al., 2012; Damashek et al., 2015; Hampel et al., 2018; Li et al., 2018), and wastewater treatment plants (Alves et al., 2018), where depth distribution of AOA ecotypes may not be applicable, and NH4+ concentrations often exceed those in the oligotrophic ocean by orders of magnitude. Thus, some AOA clades likely inhabit a more eutrophic niche than previously determined in ocean water column studies (Santoro et al., 2010; Sintes et al., 2013), but more work is needed to understand active ecotypes in shallow, freshwater systems.
Environmental controls on nitrification have been studied extensively in marine and coastal ecosystems (Beman et al., 2012; Horak et al., 2013; Urakawa et al., 2014; Damashek et al., 2016; Santoro et al., 2017), while studies in freshwater systems (Hugoni et al., 2013; Vissers et al., 2013; Bollmann et al., 2014; Hayden and Beman, 2014), particularly eutrophic systems, have received less attention (Hampel et al., 2018) and are often limited to sediment studies (Hou et al., 2013; Yao et al., 2018). Work in marine and coastal environments has shown that oxygen, temperature, salinity, light, turbidity, top-down grazing pressure (Lavrentyev et al., 1997; Jochem et al., 2004), and ambient NH4+ concentrations can all influence nitrification rates (Heiss and Fulweiler, 2016; Damashek and Francis, 2018). However, in eutrophic lakes, competition for NH4+ between nitrifiers and photoautotrophs, including bloom-forming cyanobacteria (cyanoHABs), may be more important than ambient NH4+ concentrations (which can be very high). Nitrifiers and cyanobacteria are both highly competitive for NH4+ and urea, and while cyanobacteria can assimilate N2 and NO3–, nitrifiers cannot. Cyanobacteria can also thrive in high light conditions, which may inhibit nitrifiers (Hayden and Beman, 2014). A recent study showed that this competition can reduce nitrification efficiency (Hampel et al., 2018), potentially impacting the capacity for denitrification and removal of excess N.
We examined environmental controls on nitrification rates and ammonia oxidizer community structure along a freshwater-estuarine continuum during cyanoHABs and before and after passage of a major hurricane. Lake Okeechobee and the St. Lucie Estuary (Florida) provide an opportunity to investigate the impacts of dissolved inorganic N (DIN) concentrations, salinity, cyanoHABs, and physical disturbance (hurricane passage) on nitrification. Lake Okeechobee has experienced toxic cyanoHABs for decades (James et al., 2009), resulting in myriad negative effects on the downstream St. Lucie Estuary. To prevent flooding, nutrient-laden water from Lake Okeechobee is released into St. Lucie Estuary via a canal system, reducing salinity and increasing nutrient concentrations and turbidity (Lapointe et al., 2012; Phlips et al., 2012). Toxic Microcystis blooms in St. Lucie Estuary have been observed on numerous occasions (Phlips et al., 2012), including during a massive bloom in 2016, which led the State of Florida to declare a state of emergency (Kramer et al., 2018).
On September 10th, 2017, Hurricane Irma made landfall in the Florida Keys as a Category 4 hurricane (sustained winds >210 km h–1) and moved through central Florida with sustained winds of ∼65 km h–1 in the Okeechobee region (Hampel et al., 2019) and maximum winds up to ∼100 km h–1 (SFWMD). High winds and precipitation runoff have major impacts on hydrology, microbial and phytoplankton community structure, sediments, and nutrients, including reducing salinity in estuaries, often leading to stratification and bottom-water hypoxia (Paerl et al., 2001). Hurricane winds also produce strong currents and seiches, which resuspend sediments and increase total suspended solids (TSS) and turbidity in the water column (James et al., 2008). Short-term hurricane effects on cyanoHABs include decreased NH4+ uptake rates and microcystin synthetase gene abundance (Hampel et al., 2019). Previous studies showed decreased AOA abundance and diversity after an oil spill (Newell et al., 2014) and increased AOB abundance after storm disturbance (Happel et al., 2018). However, very little is known about storm effects on nitrification and ammonia oxidizer community structure in freshwater and estuary ecosystems (Bouskill et al., 2011; Newell et al., 2014; Happel et al., 2018).
Despite its ecological importance, nitrification in eutrophic, cyanoHABs-dominated systems remains understudied. The goal of this study was to quantify nitrification rates and amoA gene abundance, and investigate ammonia oxidizer community structure, along a freshwater-estuarine continuum in central Florida during cyanoHABs in 2016 and 2017 and before and after Hurricane Irma passage. We hypothesized that nitrification rates would be lower during cyanoHABs due to resource competition and that hurricane disturbance would further decrease nitrification rates and the abundance of ammonia oxidizers. We also speculated that there would be a shift in amoA community structure as an immediate response to hurricane disturbance. Climate change predictions forecast escalation of eutrophication, cyanoHABs, and the intensity of large-scale extreme weather events (Paerl et al., 2019). Thus, understanding the effects of these disturbances on nitrification and ammonia oxidizer communities is important for constructing and validating ecosystem models used to evaluate ecosystem resilience and inform management action in eutrophic waters.
Materials and Methods
Sampling
Lake and estuary water samples were collected on July 25–27, 2016, August 22–24, 2017, and September 20–21, 2017. The July 2016 sampling occurred during a large cyanoHAB in both systems (Kramer et al., 2018). In August 2017, sampling occurred 2 weeks before Hurricane Irma passed over the Okeechobee region as a Category 3 hurricane. The September 2017 sampling occurred 10 days after Hurricane Irma passed through south and central Florida. In July 2016, sampling was conducted at two stations in Lake Okeechobee, L004 (surface and bottom water) and LZ40 (surface and bottom water), and four stations in the St. Lucie Estuary following a freshwater to marine salinity gradient from SLE80, SLE2, SLE4, and SLE8 (surface water only; Figure 1). In August 2017, L004 (surface and bottom water) and a northern lake station (SAV165, surface water only) were sampled, along with SLE80, SLE5, and SLE7 in the estuary. Post-hurricane sampling was hindered by flooding and included stations L004, LZ40, and LOBG (nearshore in the south) in the lake (surface water only), and SLE5 and SLE7 in the estuary.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/554278/fmicb-11-01541-HTML/image_m/fmicb-11-01541-g001.jpg)
Figure 1. Map of sampling locations in Lake Okeechobee and St. Lucie Estuary. Annotations indicate sampling date (J16, July 2016; A17, August 2017; S17, September 2017). Modified from Hampel et al. (2019).
Water for nutrient analyses was filtered in the field to 0.2 μm with pre-rinsed syringe filters (Nylon) into 15 ml polypropylene tubes and frozen upon return to the laboratory. Water for nitrification experiments was collected with a van Dorn water sampler into 5 L cubitainers and returned to the lab for incubations. Physicochemical parameters were measured with a YSI multisensor sonde in July 2016 and August 2017, and a Manta 2 multiparameter sonde in September 2017. Dissolved nutrient analyses included NH4+, NO2–, NO3–, orthophosphate (ortho-P), and urea were analyzed using a Lachat Quickchem 8500 FIA nutrient analyzer. Turbidity and TSS measurements were obtained from the DB Hydro database1 and collected by South Florida Water Management District (SFWMD) according to standard EPA protocols.
DNA Collection and Extraction
Environmental DNA was collected in August and September 2017. Samples were hand filtered in the field with 0.2 μm Sterivex filters (SVGP01015; MilliporeSigma, MA, United States). Samples were preserved in the field with ∼2 mL of Ambion RNAlater (Invitrogen, Carlsbad, CA, United States). The volume of site water pushed through Sterivex filters in August 2017 ranged from 120–240 ml. However, due to increased turbidity after the hurricane, only 45–60 ml of water was filtered in September for stations L004, LZ40, SLE5, and SLE7. An exception was station LOBG, where 300 ml of water were filtered. Preserved filters were frozen at −80°C until analysis.
Filter cartridges were thawed on ice in the lab, and RNAlater in the Sterivex filters was removed by pushing 10 ml of Phosphate Buffer Saline 1× Solution (Fisher BioReagents, United States) through the filter. DNA was extracted using the Gentra PureGene kit (Qiagen Inc., United States). Lysis buffer (0.9 ml) and Proteinase K (10 μl) were added to the filters, and the samples were incubated for 1 h at 55°C and 1 h at 65°C, repeated twice, in a hybridization oven at 90 rpm (Newell et al., 2011; Hampel et al., 2018).
amoA Gene Abundance
Quality and quantity of the extracted DNA were measured spectrophotometrically (NanoDrop 2000, Thermo Scientific). There are limitations to using the NanoDrop at low DNA concentrations, but our samples had consistently high concentrations (>10 ng/μL). Archaeal amoA was quantified with Arch-amoAF and Arch-amoAR primers, and bacterial amoA was quantified with amoAF and amoA2R primers (Supplementary Table S4; Rotthauwe et al., 1997; Francis et al., 2005). Archaeal and bacterial amoA qPCR standards were prepared by cloning the fragment of interest with the TOPO TA Cloning Kit (Invitrogen, United States), inserting it into a competent cell plasmid (One Shot E. coli cells, Invitrogen, United States), and isolating the plasmid containing the amoA gene using the UltraClean Standard Mini Plasmid Prep Kit (Mo Bio Laboratories Inc., Carlsbad, CA, United States). Each qPCR run included three no template controls (NTC), six standards from serial dilution in triplicate, and the environmental DNA samples in triplicate. Each sample and standard received 10 μl of Luna Universal qPCR Master Mix (New England Biolabs Inc., United States), 0.5–1 μl of each 10 μM primer, and 20–35 ng of DNA template.
amoA qPCR protocols followed methods of Bollmann et al.(2014; and Supplementary Table S4), followed by melting curve analysis to ensure a single product. Automatic settings for the thermocycler (Realplex, Eppendorf) were used to determine threshold cycle (Ct values), efficiency (93–98%), and a standard curve with R2 values ≥0.99. Gene copy number was calculated as (ng ∗ number mol –1)/(bp ∗ ng g–1 ∗ g mol–1 of bp) and is reported in amoA gene copies ml–1 of sample water.
amoA Amplicon Sequencing and Analysis
Barcoded amplicons for AOA and AOB were prepared by using a two-step PCR approach (Herbold et al., 2015). In the first PCR, bacterial and archaeal amoA fragments were amplified using the standard procedures and primers listed above, with a decreased number of cycles (24 cycles). In the second PCR, the first-step PCR product was used as a template for eight cycle amplification with a primer consisting of the barcodes and Illumina adapters (Supplementary Table S4). The barcoded amplicons were sequenced on an Illumina MiSeq for paired end sequencing (2 × 300 bp PE). Barcoding, cleaning, quantification, and pooling was conducted at the Ohio State Molecular and Cellular Imaging Center, and the Illumina MiSeq Reagent kit v3, paired end 300 sequencing kit, and the Illumina Nextera XT index kit v2 for indexing the libraries were used. Raw sequence data analysis was conducted using Mothur pipeline (version 1.42.3; Schloss et al., 2009). Sequences were assembled into contigs, and adapters and primers were removed. Sequences longer than 555 bp for AOA, and 475 bp for AOB, and low-quality sequences with more than one mismatch in barcode/primer sequences, were removed in Mothur. Frameshift errors in unique archaeal and bacterial reads were corrected using the FunGene FrameBot tool (Wang et al., 2013), and sequence alignment was conducted using aligned archaeal or bacterial amoA sequences retrieved from the FunGene database (Aalto et al., 2018). Chimeric sequences were removed using Uchime in Mothur (Edgar et al., 2011), and sequences were clustered into OTUs at 95% sequence identity.
The taxonomic assignment of archaeal and bacterial sequences was completed in Mothur using the Alves et al. (2018) comprehensive AOA amoA database and custom AOB amoA reference libraries, which were collected from sequences and their taxonomies in the Nucleotide database (Aalto et al., 2018) at a 95% cutoff (Beman and Francis, 2006; Frank et al., 2016; Aalto et al., 2018). Phylogenetic analysis was conducted in MEGA X, where sequences were aligned using the multiple sequence alignment tool, ClustalW, to align representative and reference sequences. Maximum likelihood phylogenetic trees were constructed with 1000 bootstrap replicates using the Tamura-Nei nucleotide substitution model in MEGA X (Kumar et al., 2018). The phylogenetic tree was edited and annotated using iTOL v.5 (Letunic and Bork, 2019). Raw sequences were submitted to NCBI’s Sequence Read Archive under BioProject #PRJNA592084.
Nitrification Rates
Nitrification rates were measured using the 15NH4+ tracer addition method (Ward, 2008; Newell et al., 2013; Heiss and Fulweiler, 2016). Five hundred ml of water from each station was distributed into acid-washed, 1 L polycarbonate bottles and enriched with 98% 15NH4Cl (Isotec) at final concentrations of 0.25 μM (2016) or 0.5 μM (2017). Amended water was mixed thoroughly by inverting 10 times and distributed into three, acid-washed 125 ml polycarbonate incubation bottles. Unamended control samples for each station were distributed into 125 ml incubation bottles. Initial samples (T0) were filtered to 0.22 μm with syringe filters into 50 ml polycarbonate tubes and frozen until analysis, and final samples (Tf) were collected after incubating for 20 h at near in situ light and temperature in a laboratory incubator.
Accumulation of 15NO3– was measured using the Cd reduction/NaN3 reduction method (McIlvin and Altabet, 2005; Heiss and Fulweiler, 2016). Approximately 25 ml from each sample was filtered (0.2 μm syringe filter) into 50 ml centrifuge tubes. NO3– was reduced to NO2– by addition of 100 mg of MgO, 6.6 g of NaCl, and 0.75–1 g of acidified Cd powder to each sample, followed by a 17 h incubation on a shaker table (McIlvin and Altabet, 2005). Samples were centrifuged at 1000 × g for 15 min, and 7.5 ml of supernatant was carefully transferred into 12 ml Exetainers.
NOx (including cadmium-reduced NO2–) was further reduced to N2O with the NaN3 method (McIlvin and Altabet, 2005; Newell et al., 2011; Hampel et al., 2018). Briefly, each sample was treated (with gastight syringe) with 0.25 ml of a 1:1 (v:v) solution of 2 M NaN3:20% CH3COOH solution (Ar-purged) and incubated for 1 h at 30°C (McIlvin and Altabet, 2005). NO2– in the sample from Cd reduction was transformed chemically to N2O. After a 1 h incubation, the reaction was stopped with an injection of 0.15 ml of 10 M NaOH. Samples were inverted and sent to the University of California Davis Stable Isotope Facility for isotopic analysis of 45/44N2O using a Thermo Finnigan GasBench + PreCon trace gas concentration system interfaced to a Thermo Scientific Delta V Plus isotope-ratio mass spectrometer (Bremen, Germany). Nitrification rates were corrected for NaN3 reduction efficiency (determined from reduced standards), and 15NO3– production was calculated as:
Where α = [15NH4+]/([15NH4+] + [14NH4+]).
Statistical Analysis
Statistical analysis was conducted using RStudio software (Version 1.1.463). First, environmental data, nitrification rates, and amoA abundance were checked for normality using the Shapiro-Wilk test. After determining that the data were not normally distributed, the Kendall correlation method for non-parametric data was used to determine correlations between parameters. Wilcoxon and Kruskal-Wallis tests for non-parametric data were used to determine differences between sampling events. Additionally, a stepwise multiple-regression model for nitrification rates was constructed using the MASS package (R Version 7.3). The best-fitting model was selected based on the minimum Akaike’s information criteria (AIC; Akaike, 1987). All non-normally distributed variables were log(x + 1)-transformed prior to running the model to normalize the data for parametric analysis.
Community structure analysis and visualization was performed using phyloseq (Version 1.28.0; McMurdie and Holmes, 2013), vegan (Version 2.5-2; Oksanen et al., 2018), and ggplot2 (Version 3.2.0; Wickham, 2016) packages in RStudio. Alpha diversity was calculated for bacterial and archaeal data. Bray-Curtis dissimilarity analysis was used to evaluate variability of communities across samples, using non-metric multidimensional scaling (NMDS; Supplementary Figure S1). Differences in the AOA and AOB community structure between sampling events were tested separately with permutational multivariate analysis of variance (PERMANOVA; Anderson, 2001; adonis in vegan package), followed by the homogeneity of dispersion test using the betadisper function. Relationships between sequencing data and environmental variables were evaluated using the Constrained Analysis of Principal Coordinates (CAP) ordinations after removing missing values from environmental data (Supplementary Figure S2).
Results
Physical and Biogeochemical Parameters
Environmental data have been reported previously (Hampel et al., 2019) and are presented here in Supplementary Tables S1, S2. Water temperature in the lake and estuary was uniform across all sampling events. After the hurricane, salinity and NH4+ concentrations decreased in the estuary (p = 0.03, p = 0.08, respectively), while conductivity (p = 0.005), turbidity (p = 0.001), TSS (p = 0.005), and NO3– concentrations (p = 0.001) increased (Supplementary Tables S1, S2). Post-hurricane, lake turbidity (p = 0.04) and NO3– concentration (p = 0.05) increased. In the lake, chlorophyll a (chl a) concentration was higher during the August 2017 bloom than during the 2016 bloom and decreased after the hurricane (Supplementary Table S1). In the estuary, chl a concentrations decreased along the salinity gradient and were lowest following the hurricane (Supplementary Table S2).
Nitrification Rates
Nitrification rates in the lake in July 2016 varied between the eastern (L004; Figure 1) and central lake stations (LZ40; Figure 2A). Nitrification rates were higher in bottom water (2 m) at L004 (98.8 ± 5.11 nmol L–1 d–1; mean ± SE) compared to surface water (∼0.2 m; 10.5 ± 1.6 nmol L–1 d–1; p = 0.03). In contrast, nitrification rates in the bottom water at LZ40 were undetectable, and low rates were measured in surface water (2.04 ± 0.52 nmol L–1 d–1). Nitrification rates in the lake in August 2017 were only measurable in surface water and were slightly lower than those measured in 2016 (0.78 ± 0.82 nmol L–1 d–1 at SAV165 to 2.14 ± 0.02 nmol L–1 d–1 at L004; Figure 2A). However, after the hurricane, nitrification rates increased by 2–3 orders of magnitude (p = 0.003) at L004 (1280 ± 71.3 nmol L–1 d–1) and LZ40 (168 ± 96.1 nmol L–1 d–1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/554278/fmicb-11-01541-HTML/image_m/fmicb-11-01541-g002.jpg)
Figure 2. Nitrification rates ± one standard error in: (A) Lake Okeechobee and (B) St. Lucie Estuary. Note the scale difference between (A) and (B).
In July 2016, nitrification rates in the estuary increased along the salinity gradient, with undetectable rates at the most upstream site (SLE80) and greatest rates at the most saline site (SLE8; 79.4 ± 0.01 nmol L–1 d–1; Figure 2B). In August 2017, nitrification rates in the estuary were higher than in 2016 (p = 0.06). Again, the lowest rates were observed upstream near the St. Lucie canal (SLE80; 79.1 ± 2.22 nmol L–1 d–1), and greatest rates were observed downstream at SLE5 (978 ± 294 nmol L–1 d–1; p = 0.04) and SLE7 (563 ± 0.4 nmol L–1 d–1; p < 0.001; Figure 2B). Following the hurricane, nitrification rates increased at both upstream (SLE5: 3590 ± 98.9 nmol L–1 d–1) and downstream stations (SLE7: 4030 ± 888 nmol L–1 d–1; Figure 2B; p < 0.001). Overall, nitrification rates in the estuary exceeded those in the lake on all occasions. Nitrification rates were negatively correlated with temperature and dissolved oxygen (DO; Kendall, p = 0.03; Table 1) and positively correlated with ambient NO3– concentrations (Kendall, p = 0.001). Additionally, the best-fitting multiple regression model for nitrification revealed that ambient NO3– concentration, turbidity, and salinity were the best predictors for nitrification rates in the lake and estuary (Adj. R2 = 0.99; Supplementary Table S3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/554278/fmicb-11-01541-HTML/image_m/fmicb-11-01541-t001.jpg)
Table 1. Details of Kendall correlation for non-parametric data analysis of environmental parameters, geochemical rates, and amoA gene abundance in Lake Okeechobee and St. Lucie Estuary.
amoA Abundance
In the lake in August 2017, AOA (mean = 8.21 ± 0.70 × 105; range: 8.21 × 104–1.83 × 106) exceeded AOB abundance (mean = 1.54 ± 0.29 × 104; range: 2.31–3.29 × 103) at all stations, with AOA abundance 1–3 orders of magnitude higher (Figure 3A). Following the hurricane, AOA abundance (mean = 3.06 ± 1.32 × 106; range: 4.32 × 104–5.73 × 106) was also higher than AOB (mean = 4.87 ± 0.43 × 104; range: 2.07–9.35 × 104) at all stations, except near the lake shore (LOBG; Figure 3A). AOB abundance increased after Hurricane Irma in the lake (p = 0.05), while AOA abundance did not change. Lake nitrification rates were correlated with AOA amoA copy numbers (Figure 4; R2 = 0.94, p = 0.006), while no pattern was observed between nitrification rates and lake AOB or estuarine ammonia oxidizers (Supplementary Figure S3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/554278/fmicb-11-01541-HTML/image_m/fmicb-11-01541-g003.jpg)
Figure 3. Ammonia oxidizer abundance (amoA gene) ± one standard error in Lake Okeechobee (A) and St. Lucie Estuary (B) before and after Hurricane Irma.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/554278/fmicb-11-01541-HTML/image_m/fmicb-11-01541-g004.jpg)
Figure 4. Pearson correlation between the abundance of ammonia oxidizing archaea amoA functional gene and log-transformed nitrification rates in Lake Okeechobee. Black points indicate samples collected in August 2017; red points indicate samples collected in September 2017.
In the estuary in August 2017, AOA abundance (mean = 5.36 ± 4.39 × 104; range: 7.10 × 103–1.29 × 105) exceeded AOB (mean = 2.07 ± 0.36 × 104; range: 5.25 × 102–9.37 × 103) at all stations except upstream (SLE80). AOA abundance was slightly lower than in the lake, while AOB abundances were comparable to those measured in the lake (Figure 3B). After Hurricane Irma, amoA gene copies for both AOA and AOB increased, and AOA (mean = 2.71 ± 0.63 × 105; range: 7.65 × 104–4.65 × 105) outnumbered AOB (mean = 1.03 ± 0.13 × 105; range: 5.06 × 104–1.56 × 105) at all stations (Figure 3B). AOA abundance in both the lake and estuary were negatively correlated with temperature and ambient NH4+ concentration and positively correlated with TSS (Kendall, p = 0.003; Table 1). AOB abundance was positively correlated with conductivity (Kendall, p = 0.03).
Ammonia Oxidizer Community Structure
Following the quality control steps in Mothur, frameshift in FrameBot, and chimera removal, we retrieved 3,699 unique amoA sequences for AOA and 17,030 for AOB. amoA sequences grouped into 574 and 101 OTUs (95% similarity) for AOA and AOB, respectively. We observed a shift in the AOA community after the hurricane (Figure 5; PERMANOVA, F = 2.763, p = 0.03), while the AOB community changed only marginally (F = 1.592, p = 0.08).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/554278/fmicb-11-01541-HTML/image_m/fmicb-11-01541-g005.jpg)
Figure 5. Relative abundance of ammonia-oxidizing archaea (A) and bacteria (B) in Lake Okeechobee and St. Lucie Estuary before (August) and after (September) hurricane disturbance. AOA represent three major lineages: Nitrosopumilales (NP), Nitrososphaerales (NS), and Nitrosotaelales (NT). Lineages are divided into clades (Greek letters) and numbers (genera).
Ammonia-oxidizing archaeal sequences grouped phylogenetically into the lineages and clades of Nitrososphaerales (NS), Nitrosotaleales (NT), and Nitrosopumilales (NP): NP η (Nitrosotenuis spp.), NP γ 2.1 (Nitrosopumilus spp.), NP γ 2.2 (Nitrosoarchaeum spp.), and NP δ (Figure 6). Sequences from the lake and estuary were mostly phylogenetically distinct, as OTUs present only in the lake were phylogenetically close to NS and NT orders, while the estuarine-only OTUs assigned to the NP γ 2.1 clade.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/554278/fmicb-11-01541-HTML/image_m/fmicb-11-01541-g006.jpg)
Figure 6. Phylogenetic tree of ammonia-oxidizing archaea OTUs classified from amoA reads at 95% similarity with reference sequences (Alves et al., 2018). OTU label colors indicate OTUs found in Lake Okeechobee (green), St. Lucie Estuary (blue), or both systems (black). Numbers in parenthesis indicate the number of OTUs present from each station. Nodes with branch length <0.05. Some nodes were further collapsed for clarity.
Before the hurricane, AOA were dominated by unclassified NP γ 2.2 (Nitrosoarchaeum spp.) in the lake and NP γ 2,1 (Nitrosopumilus spp.) OTUs in the estuary (Figure 5A). After the hurricane, AOA in the lake shifted to NP η 1.1 (Ca. Nitrosotenuis cloacae) OTUs in the central lake (L004, LZ40), while most OTUs from the nearshore station (LOBG) clustered within the NS δ clade. The post-hurricane AOA community in the estuary shifted to higher relative abundances of NP η 1.1 (Ca. N. cloacae) upstream (SLE5) and NP η 1.1 and NP γ 2.2 OTUs downstream (SLE7).
Ammonia-oxidizing bacterial sequences were not as clearly separated by sampling location as AOA sequences and were placed phylogenetically into the Nitrosomonas clusters, with most of the OTUs close to Nitrosomonas ureae, and some estuary-only OTUs close to Nitrosomonas nitrosa (Figure 7). Before the hurricane, AOB communities were dominated by Nitrosomonas oligotropha in the upstream estuary (SLE80 and SLE5) and by unclassified Nitrosomonas groups in the lake (Figure 5B). After the hurricane, both lake and estuary AOB communities shifted to a higher proportion of Nitrosomonas oligotropha.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/554278/fmicb-11-01541-HTML/image_m/fmicb-11-01541-g007.jpg)
Figure 7. Phylogenetic tree of ammonia-oxidizing bacteria OTUs classified from amoA reads at 95% similarity. OTU label colors indicate OTUs found in Lake Okeechobee (green), St. Lucie Estuary (blue), or both systems (black). Numbers in parenthesis indicate the number of OTUs present from each station. Nodes with branch length <0.05. Some nodes were further collapsed for clarity.
Overall, despite an order of magnitude fewer sequences, AOA were more diverse than AOB (p = 0.009; Figure 8). AOB alpha diversity increased after the hurricane in both systems (p = 0.02), while AOA diversity remained unchanged (Figure 8). Constrained Analysis of Principal coordinates (CAP) for AOA showed separation by sampling month and the effects of hurricane disturbance (Supplementary Figure S2A). CAP1 (57.9% variance explained) was driven by salinity (r = −0.95), NH4+ concentration (r = −0.91), and turbidity (r = 0.77), which were also affected by the hurricane, while CAP2 (25.9% variance explained) was driven primarily by NO3– (r = −0.61) and DO concentrations (r = 0.53). For AOB communities (Supplementary Figure S2B), CAP1 (35.2% variance explained) was associated with salinity (r = 0.80), NH4+ concentration (r = 0.76), and nitrification rate (r = −0.63), while CAP2 (20.7% variance) was associated with salinity (r = 0.54) and turbidity (r = −0.51).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/554278/fmicb-11-01541-HTML/image_m/fmicb-11-01541-g008.jpg)
Figure 8. Shannon diversity before and after Hurricane Irma in Lake Okeechobee and St. Lucie Estuary for ammonia oxidizing archaea and ammonia oxidizing bacteria.
Discussion
Nitrification During Cyanobacterial Blooms
Nitrification rates in both the lake and estuary varied between sampling events (Figure 2). Nitrification rates in Lake Okeechobee were low during cyanoHABs in 2016 and 2017 (Figure 2A), similar in magnitude to those from oligotrophic Lake Superior (Small et al., 2013), and lower than those previously measured in spring in Lake Okeechobee (James et al., 2009). Competition for NH4+ between nitrifiers and photoautotrophs, including cyanobacteria, may inhibit nitrification during severe cyanoHABs, which were observed in Lake Okeechobee in 2016 (Kramer et al., 2018; Hampel et al., 2019). NH4+ is the most energetically favorable N form for primary producers, including non-N2 fixing cyanobacteria (e.g., Microcystis), which often outcompete other algal groups (e.g., Blomqvist et al., 1994) and ammonia oxidizers for NH4+. Nitrifiers have a lower substrate affinity and higher Km (Martens-Habbena et al., 2009) for NH4+ than Microcystis (Nicklisch and Kohl, 1983; Baldia et al., 2007). In Lake Okeechobee, nitrification accounted for <1% of total NH4+ uptake and was likely suppressed by competition with Microcystis (Hampel et al., 2019). Suppression of nitrification due to competition for NH4+ between Microcystis and ammonia oxidizers has also been observed in hypereutrophic Lake Taihu (China), where nitrification rates during massive Microcystis blooms were equally as low as reported here in Lake Okeechobee, but increased when the bloom was less pronounced (Hampel et al., 2018).
Water column nitrification measurements in lakes are scarce compared to open ocean and coastal marine systems, but the rates measured here are comparable to other systems. In the estuary, nitrification rates in 2016 ranged from undetectable to 79.4 ± 0.02 nmol L–1 d–1 and increased along the salinity gradient (Figure 2B). These rates are similar to those measured in other coastal marine systems (Horak et al., 2013; Damashek et al., 2016; Heiss and Fulweiler, 2016; Laperriere et al., 2018). Undetectable and low nitrification rates in the upstream estuary are likely due to the presence of cyanoHABs during sampling and increased competition for NH4+. High photoautotrophic NH4+ uptake rates (2.32 μmol L–1 h–1; Hampel et al., 2019), coincident with high Microcystis cell densities in the upstream estuary (Kramer et al., 2018), suggest that nitrifiers were outcompeted for NH4+. In contrast, the bloom was less pronounced downstream (Kramer et al., 2018), and photoautotrophic NH4+ uptake rates were lower (1.11 μmol L–1 h–1; Hampel et al., 2019), which may help explain the increase in nitrification rates along the salinity gradient.
Before the hurricane, flows of lake water into the estuary were low in 2017, and cyanoHABs were mostly constrained to the lake (Hampel et al., 2019). Nitrification rates in the estuary were an order of magnitude higher in August 2017 than in the previous year (Figure 2B). The highest nitrification rates were observed at the more saline SLE5 and SLE7 (563–978 nmol L–1 d–1; Figure 2B), and nitrification accounted for 3–7% of total NH4+ uptake, suggesting a more effective competition for NH4+. Nitrification rates of this magnitude have been observed in other high-nutrient environments, such as the Eastern Tropical South Pacific (Lipschultz et al., 1990), the hypoxic zone in the Gulf of Mexico (Bristow et al., 2015), coastal Georgia (United States; Tolar et al., 2017), and Hood Canal, Puget Sound (United States; Urakawa et al., 2014). These results illustrate indirect downstream consequences of HABs along the hydrologic continuum, where HABs and nutrient loading upstream in the estuary can suppress nitrification (as seen in 2016), therefore potentially reducing N removal via denitrification.
Nitrification Following Hurricane Irma
After the hurricane, nitrification rates in both systems increased by an order of magnitude or more (Figure 2). In the lake, nitrification rates were up to 1280 nmol L–1 d–1 (Figure 2A). Nitrification rates on this scale in natural, freshwater systems have, to our knowledge, only been reported in Lake Taihu (Hampel et al., 2018) and Lake Mendota (Hall, 1986). In the estuary, nitrification rates after the hurricane ranged from 3590–4030 nmol L–1 d–1, much higher than rates reported previously in oligotrophic and some eutrophic estuarine systems (Horak et al., 2013; Bronk et al., 2014; Urakawa et al., 2014; Heiss and Fulweiler, 2016). Rates in this range were previously measured in the Chang Jiang (Yangtze) River plume (Hsiao et al., 2014), eutrophic Elbe estuary (Sanders et al., 2018), Pearl River estuary (Dai et al., 2008), and northern Gulf of Mexico (Carini et al., 2010). The dramatic increase in nitrification following the hurricane may be related to sediment resuspension caused by high winds, introducing benthic nitrifiers and NH4+ accumulated in sediments into the water column. With more favorable oxygen conditions in the water column, a fresh supply of NH4+, and higher turbidity, photoautotrophs may have experienced light limitation (James et al., 2008), thus releasing competitive pressure and allowing nitrification to increase (Wengrove et al., 2015).
A comprehensive review of nitrification in estuaries reported a trend of elevated nitrification rates at Estuary Turbidity Maxima (ETM; Damashek and Francis, 2018), some of which were also associated with turbid waters and strong wind events (Hsiao et al., 2014; Happel et al., 2018; Laperriere et al., 2018; Sanders et al., 2018). NH4+ released from resuspended sediments can thus lead to nitrification “hot spots” (Percuoco et al., 2015; Damashek et al., 2016). Hurricanes are also accompanied by increased rainfall, leading to larger runoff and increased nutrient concentrations in the water column (Havens et al., 2001; James et al., 2008). High turbidity in the water column may limit photoautotrophic access to light, thus favoring nitrifiers, which are chemolithoautotrophic and do not require light. Microcystin synthetase (mcyD) abundance and chl a levels decreased following the hurricane (Hampel et al., 2019), supporting this interpretation. In the estuary, nitrification accounted for 60–90% of total NH4+ uptake after the hurricane, further supporting the idea that turbid conditions favored nitrifiers.
Ammonia-Oxidizer Community Structure and Abundance
For most sampling events, AOA abundance was greater than AOB (Figure 3). While AOA are generally more abundant than AOB in the open ocean, estuarine ammonia-oxidizing communities are more complex and dynamic, with different systems exhibiting AOA dominance (Beman and Francis, 2006; Horak et al., 2013), AOB dominance (Bernhard et al., 2010), or interchanging spatial dominance of AOA and AOB (Santoro et al., 2008; Zhang et al., 2014). Salinity can help shape ammonia oxidizer community structure (Caffrey et al., 2007; Mosier and Francis, 2008; Santoro et al., 2008). AOB often dominate lower salinity regions of estuaries, while AOA are more abundant in more marine areas (Bouskill et al., 2012), but this pattern is not universal (Mosier and Francis, 2008; Santoro et al., 2008). In the estuary, AOB abundance was only greater than AOA at the upstream, lower salinity site (SLE80).
Despite AOA outnumbering AOB in many systems, AOB abundance has more often correlated positively with nitrification rates (Bernhard et al., 2010; Zeng et al., 2012; Damashek et al., 2015; Hampel et al., 2018). However, the only observed relationship between amoA abundance and nitrification rate was between AOA amoA copy number and log-transformed nitrification rates in Lake Okeechobee (R2 = 0.94; Figure 4 and Supplementary Figure S3). A relationship between AOA amoA abundance and potential nitrification rates in sediments was reported in hypereutrophic Lake Taihu (Zeng et al., 2012), but we are not aware of any similar relationship reported in a eutrophic lake water column to date.
This study presents a remarkable phylogenetic diversity of recovered AOA OTUs, with close relatives from three of the four major lineages (Alves et al., 2018). Unlike AOB, AOA rarefaction curves (Supplementary Figure S4) did not plateau, suggesting that much of the diversity for AOA was not characterized. Despite the relatively low number of AOA sequences recovered, AOA were still more abundant and diverse than AOB. The low number of sequences, but high diversity and abundance, justifies not rarefying our sequences in the downstream analysis. Recent studies suggest that normalizing sequences via rarefaction can result in discarding important reads and introduces more bias (McMurdie and Holmes, 2014). Overall, the high diversity (even at comparatively low abundance) of AOA sequences along with strong correlation with nitrification rates emphasizes the possible importance of AOA in eutrophic systems and their contribution to nitrification.
The high abundance and diversity of AOA, and strong relationship with nitrification rates within Lake Okeechobee, is puzzling considering the reported oligotrophic nature of AOA (Martens-Habbena et al., 2009; Bollmann et al., 2014; Kits et al., 2017). We suggest that AOA in these eutrophic ecosystems might have evolved from niche differentiation and could represent a new, eutrophic ecotype, separate from oligotrophic AOA (Beman et al., 2008; Sintes et al., 2013; Santoro et al., 2017). High phylogenetic diversity of AOA in this and other studies (Alves et al., 2018; Cheung et al., 2019; Zhao et al., 2020) suggests that broad classification into shallow and deep marine and terrestrial ecotypes is insufficient, and many amoA clades represent ecotype intra-diversity, perhaps dictated by other environmental conditions beyond depth, salinity, and pH.
The majority of lake-only AOA OTUs clustered within the NS δ clade, the most abundant clade in the Alves et al. (2018) database that still lacks cultivated representatives. Although little is known about the physiological characteristics of the NS δ clade, it is abundant in soils, freshwater water columns and sediments, marine sediments, wastewater treatment plants, and salt lakes, suggesting an exceptional versatility of habitats for this group. A BLAST search of the representative sequences from this clade matched with previous studies from sludges (Xia et al., 2019), soils (Li et al., 2015; Zhong et al., 2016), and streams affected by wastewater treatment plants (Merbt et al., 2015). Other lake OTUs clustered in NT β and NP δ and η. The NT β clade has no isolated representative and is mostly present in freshwaters, soils, and sediments, while NP δ is also uncharacterized, with sequences from freshwater, estuarine, and marine environments (Alves et al., 2018). These lake-only OTUs, particularly from the NS δ clade at LOBG, could belong to a eutrophic AOA ecotype that has adapted to high NH4+ concentrations, but evidence from the literature is lacking on AOA communities in lakes. The NH4+ concentration at this station during sampling was ∼20 μM (Supplementary Table S1), which greatly exceeds typical concentrations in the ocean. This eutrophic AOA ecotype may be able to outcompete AOB in eutrophic freshwater systems, as suggested from our results in Lake Okeechobee. Similar conclusions were reached in the eutrophic Pearl River estuary (Zou et al., 2019), but genomic adaptations of AOA in eutrophic systems and their contributions to nitrification remain an intriguing knowledge gap requiring further investigation.
Few studies have documented changes in ammonia oxidizer community structure following disturbances in aquatic systems. Responses of ammonia oxidizers to various perturbations differ and include increased AOB diversity or abundance (wastewater, Aalto et al., 2018; flow restriction, Bernhard et al., 2015; flooding, Bouskill et al., 2011; storms, Kan, 2018), decreased AOA diversity or abundance (storms, Happel et al., 2018; hurricanes, Newell et al., 2014), or little to no effect (oil spill, Bernhard et al., 2019; oil spill, Marton et al., 2015). Here, we observed substantial hurricane impacts on the abundance, community composition, and diversity of ammonia oxidizers. In the lake, only AOB abundance increased after the hurricane, while in the estuary, we observed the opposite, with AOA abundance increasing. AOB communities were more diverse following the hurricane, while AOA diversity increased in the estuary but decreased in the lake (Figure 8). These results agree with those from a previous study in the Gulf of Mexico (Newell et al., 2014), where benthic AOA communities were less diverse following Hurricanes Ike and Gustav, suggesting a weak resistance to disturbance. Increased AOB abundance in our study also supports disturbance studies in boreal lake sediments (Aalto et al., 2018) and estuaries (Happel et al., 2018), where AOB diversity and/or abundance increased at affected sites, suggesting that AOB are more adaptable to disruptive events (Aalto et al., 2018).
Hurricane disturbance caused a significant shift in AOA community structure in both systems (Figure 5A). In the lake, the community changed from high abundance of low-salinity tolerant NP γ 2.2 (Nitrosoarchaeum spp.) to dominance by NP η 1.1 (Ca. Nitrosotenuis cloacae) in the central lake (L004, LZ40) and NS δ OTUs at the littoral station (LOBG). Ca. Nitrosotenuis cloacae was enriched from a wastewater treatment plant and exhibits intolerance to high salinity, and phylogenetically related strains of NP η have been found in lakes and rivers (Li et al., 2016; Alves et al., 2018). Nitrososphaera are often present in soils and sediments (Alves et al., 2018; Li et al., 2018); thus, their high abundance at the littoral station after the hurricane is likely due to sediment resuspension or soil runoff from heavy rainfall. A significant shift in the community was also observed in the estuary, from the salinity tolerant Nitrosopumilus spp. (NP γ 2.1) to salinity intolerant Nitrosoarchaeum spp. (NP γ 2.2) and Ca. Nitrosotenuis cloacae (NP η 1.1), likely reflecting population transfer via water discharges from the lake into the estuary.
In contrast, AOB communities after the hurricane exhibited no community composition changes. The majority of the sequences from both systems were unclassified Nitrosomonas and Nitrosomonas oligotropha (Figures 5B, 7). A BLAST search of unclassified Nitrosomonas consensus sequences revealed that sequences from this study matched with those from other eutrophic systems, such as the Chang Jiang River basin (Jiang et al., 2017), Lake Taihu (Zeng et al., 2012; Zhao et al., 2013; Liu and Yang, 2019), San Francisco Bay (Damashek et al., 2015), an aquaculture system (Bartelme et al., 2017), and wastewater/sludge bioreactors (Wells et al., 2009; Wu et al., 2013; Zhang et al., 2018). Overall, a high abundance of Nitrosomonas-like AOB is expected due to NH4+-rich, eutrophic conditions (Stein et al., 2007; Prosser et al., 2014; Sedlacek et al., 2019) and have been found in many eutrophic natural environments (Zhang et al., 2014; Damashek et al., 2015) and N-rich engineered systems (Zhang et al., 2014; Bartelme et al., 2017). However, some Nitrosomonas species are adapted to lower NH4+ environments. Representative sequences from this study matched Nitrosomonas sp. Is79, an ammonia oxidizer found in low NH4+ freshwater environments (Bollmann et al., 2013). Additionally, sequences in this study classified as Nitrosomonas oligotropha, which is also adapted to lower trophic conditions (French et al., 2012). In the estuary, we observed an increase in Nitrosomonas oligotropha after the hurricane, possibly due to reduced salinity (Mosier and Francis, 2008), sediment resuspension, and lower nutrient concentrations (NH4+ and ortho-P; Supplementary Table S2).
The increased nitrifier abundance, shift in AOA community structure, and higher nitrification rates post-hurricane were likely caused by a combination of physicochemical factors, resuspension of benthic nitrifiers into the water column, and reduced competition with cyanobacteria. Environmental parameters explaining variance in nitrifier community structure were salinity, turbidity, and NH4+ and NO3– concentrations (Supplementary Figure S2). All of these factors were affected by the hurricane disturbance (Supplementary Tables S1, S2). The presence of primarily soil and sediment nitrifiers in water column samples likely contributed to increased nitrification rates, along with reduced competition with Microcystis. The shift (PERMANOVA; p = 0.003) in AOA community structure, but not AOB, might indicate that AOA have lower resistance to disturbance than AOB (Bowen et al., 2011; Shade et al., 2011, 2012), but a follow-up study would be required to assess resilience.
Conclusion
This study evaluated the effects of cyanoHABs and Hurricane Irma on water column nitrification rates and ammonia oxidizer community structure along a freshwater-estuarine continuum. We show that nitrification rates were lower during cyanoHABs, increasing by an order of magnitude when cyanoHABs were less abundant in the estuary the following year. We also show that several OTUs from Lake Okeechobee belong to the uncultured NS δ clade, which is versatile and appears to be ecologically important in eutrophic systems. Following Hurricane Irma, nitrification rates increased by an order of magnitude or more, coincident with increased turbidity and conductivity and decreased salinity and cyanoHABs. We also observed a substantial shift in AOA community structure immediately after the hurricane, suggesting weaker resistance to disturbance than AOB.
This study contributes to the growing literature showing high abundance of AOA in eutrophic waters relative to AOB and presents a strong relationship between AOA abundance and nitrification rates. Given their relatively recent discovery and large gaps in knowledge of AOA physiology, many open questions and suggestions for future research remain:
(1) AOA appear to be adapted to a wider range of trophic conditions than previously thought. The high abundance and diversity of AOA in eutrophic systems, and low number of cultured species (NS δ, NT β, NP δ), suggest that trophic adaptation and niche differentiation are likely, beyond the oceanic/terrestrial ecotype differentiation;
(2) AOB alpha diversity increased after the hurricane, while the AOA community shifted dominance. However, it is not certain which of these groups was responsible for the increase in nitrification rates in the estuary after the hurricane;
(3) Our results present short-term effects of hurricane and cyanoHAB disturbances on nitrification rates and ammonia oxidizer community structure. Studies focused on longer term effects should explore resilience of the ammonia oxidizer community and temporal scales of community recovery. The limited scope of the present study leaves open whether and/or how quickly the community might return to the pre-hurricane state, when turbidity is either flushed or settles to the sediment surface. Current climate change projections forecast amplification of hurricane events and increased eutrophication, and understanding the role of nitrification in converting chemically reduced N forms to oxidized N, and providing substrate for N removal via denitrification, may present valuable opportunities to constrain system N budgets.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material.
Author Contributions
JH, MM, and SN designed the study and collected samples. JH analyzed the samples. SA helped with sequence analysis. JH, SN, and MM interpreted results. All authors contributed to writing and revising the manuscript.
Funding
This work was funded by Florida Sea Grant (Project PD-16-10), the Wright State Women in Science Giving Circle, and the Eppley Foundation.
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.
Acknowledgments
We thank the South Florida Water Management District, especially Therese East, for providing lab space, boat access, and assistance with incubations. We thank Dr. Annette Bollmann for her insights to sequencing approach. We also thank Captain David Lauer for boat access on Lake Okeechobee after Hurricane Irma, Daniel Hoffman, Justin Myers, and Ashlynn Boedecker for help with sample analyses, and Megan Reed and Donnie Peterson for help with field sampling. We thank Richard Doucett and the UC Davis Stable Isotope Facility for processing nitrification samples. We thank the instructors and teaching assistants at STAMPS 2019 MBL workshop for providing valuable advice and tips regarding microbiome data analysis.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.01541/full#supplementary-material
Footnotes
References
Aalto, S. L., Saarenheimo, J., Mikkonen, A., Rissanen, A. J., and Tiirola, M. (2018). Resistant ammonia-oxidizing archaea endure, but adapting ammonia-oxidizing bacteria thrive in boreal lake sediments receiving nutrient-rich effluents: AOA and AOB in wastewater-influenced sediments. Environ. Microbiol. 20, 3616–3628. doi: 10.1111/1462-2920.14354
Akaike, H. (1987). “Factor analysis and AIC,” in Selected Papers of Hirotugu Akaike, eds E. Parzen, K. Tanabe, and G. Kitagawa (New York, NY: Springer), 371–386. doi: 10.1007/978-1-4612-1694-0_29
Alves, R. J. E., Minh, B. Q., Urich, T., von Haeseler, A., and Schleper, C. (2018). Unifying the global phylogeny and environmental distribution of ammonia-oxidising archaea based on amoA genes. Nat. Commun. 9:1517.
An, S., and Joye, S. B. (2001). Enhancement of coupled nitrification-denitrification by benthic photosynthesis in shallow estuarine sediments. Limnol. Oceanogr. 46, 62–74. doi: 10.4319/lo.2001.46.1.0062
Anderson, M. J. (2001). A new method for non-parametric multivariate analysis of variance. Aust. Ecol. 26, 32–46. doi: 10.1111/j.1442-9993.2001.01070.pp.x
Baldia, S. F., Evangelista, A. D., Aralar, E. V., and Santiago, A. E. (2007). Nitrogen and phosphorus utilization in the cyanobacterium Microcystis aeruginosa isolated from Laguna de Bay, Philippines. J. Appl. Phycol. 19, 607–613. doi: 10.1007/s10811-007-9209-0
Bartelme, R. P., McLellan, S. L., and Newton, R. J. (2017). Freshwater recirculating aquaculture system operations drive biofilter bacterial community shifts around a stable nitrifying consortium of ammonia-oxidizing Archaea and comammox Nitrospira. Front. Microbiol. 8:101. doi: 10.3389/fmicb.2017.00101
Beman, J. M., and Francis, C. A. (2006). Diversity of ammonia-oxidizing archaea and bacteria in the sediments of a hypernutrified subtropical estuary: Bahia del Tobari, Mexico. Appl. Environ. Microbiol. 72, 7767–7777. doi: 10.1128/aem.00946-06
Beman, J. M., Popp, B. N., and Francis, C. A. (2008). Molecular and biogeochemical evidence for ammonia oxidation by marine Crenarchaeota in the Gulf of California. ISME J. 2, 429–441. doi: 10.1038/ismej.2007.118
Beman, M., Popp, B. N., and Alford, S. E. (2012). Quantification of ammonia oxidation rates and ammonia-oxidizing archaea and bacteria at high resolution in the Gulf of California and eastern tropical North Pacific Ocean. Limnol. Oceanogr. 57, 711–726. doi: 10.4319/lo.2012.57.3.0711
Bernhard, A. E., Chelsky, A., Giblin, A. E., and Roberts, B. J. (2019). Influence of local and regional drivers on spatial and temporal variation of ammonia-oxidizing communities in Gulf of Mexico salt marshes. Environ. Microbiol. Rep. 11, 825–834.
Bernhard, A. E., Dwyer, C., Idrizi, A., Bender, G., and Zwick, R. (2015). Long-term impacts of disturbance on nitrogen-cycling bacteria in a New England salt marsh. Front. Microbiol. 6:46. doi: 10.3389/fmicb.2015.00046
Bernhard, A. E., Landry, Z. C., Blevins, A., de la Torre, J. R., Giblin, A. E., and Stahl, D. A. (2010). Abundance of ammonia-oxidizing archaea and bacteria along an estuarine salinity gradient in relation to potential nitrification rates. Appl. Environ. Microbiol. 76, 1285–1289. doi: 10.1128/AEM.02018-09
Blomqvist, P., Petterson, A., and Hyenstrand, P. (1994). Ammonium-nitrogen: a key regulatory factor causing dominance of non-nitrogen-fixing cyanobacteria in aquatic systems. Arch. Hydrobiol. 132, 141–164.
Bollmann, A., Bullerjahn, G. S., and McKay, R. M. (2014). Abundance and diversity of ammonia-oxidizing archaea and bacteria in sediments of trophic end members of the Laurentian Great Lakes, Erie and superior. PLoS One 9:97068. doi: 10.1371/journal.pone.0097068
Bollmann, A., Sedlacek, C., Norton, J., Norton, J. M., Laanbroek, H. J., Suwa, Y., et al. (2013). Complete genome sequence of Nitrosomonas sp. Is79, an ammonia oxidizing bacterium adapted to low ammonium concentrations. Stand. Genomic Sci. 7, 469–482.
Bouskill, N. J., Eveillard, D., Chien, D., Jayakumar, A., and Ward, B. B. (2012). Environmental factors determining ammonia-oxidizing organism distribution and diversity in marine environments: ammonia-oxidizers across aquatic environmental gradients. Environ. Microbiol. 14, 714–729. doi: 10.1111/j.1462-2920.2011.02623.x
Bouskill, N. J., Eveillard, D., O’Mullan, G., Jackson, G. A., and Ward, B. B. (2011). Seasonal and annual reoccurrence in betaproteobacterial ammonia-oxidizing bacterial population structure. Environ. Microbiol. 13, 872–886. doi: 10.1111/j.1462-2920.2010.02362.x
Bowen, J. L., Ward, B. B., Morrison, H. G., Hobbie, J. E., Valiela, I., Deegan, L. A., et al. (2011). Microbial community composition in sediments resists perturbation by nutrient enrichment. ISME J. 5, 1540–1548. doi: 10.1038/ismej.2011.22
Bristow, L. A., Dalsgaard, T., Tiano, L., Mills, D. B., Bertagnolli, A. D., Wright, J. J., et al. (2016). Ammonium and nitrite oxidation at nanomolar oxygen concentrations in oxygen minimum zone waters. Proc. Natl. Acad. Sci. U.S.A. 113, 10601–10606. doi: 10.1073/pnas.1600359113
Bristow, L. A., Sarode, N., Cartee, J., Caro-Quintero, A., Thamdrup, B., and Stewart, F. J. (2015). Biogeochemical and metagenomic analysis of nitrite accumulation in the Gulf of Mexico hypoxic zone: nitrite accumulation in hypoxic zones. Limnol. Oceanogr. 60, 1733–1750. doi: 10.1002/lno.10130
Bronk, D. A., Killberg-Thoreson, L., Sipler, R. E., Mulholland, M. R., Roberts, Q. N., Bernhardt, P. W., et al. (2014). Nitrogen uptake and regeneration (ammonium regeneration, nitrification and photoproduction) in waters of the West Florida Shelf prone to blooms of Karenia brevis. Harmful Algae 38, 50–62. doi: 10.1016/j.hal.2014.04.007
Caffrey, J. M., Bano, N., Kalanetra, K., and Hollibaugh, J. T. (2007). Ammonia oxidation and ammonia-oxidizing bacteria and archaea from estuaries with differing histories of hypoxia. ISME J. 1, 660–662. doi: 10.1038/ismej.2007.79
Carini, S. A., McCarthy, M. J., and Gardner, W. S. (2010). An isotope dilution method to measure nitrification rates in the northern Gulf of Mexico and other eutrophic waters. Cont. Shelf Res. 30, 1795–1801. doi: 10.1016/j.csr.2010.08.001
Cheung, S., Mak, W., Xia, X., Lu, Y., Cheung, Y., and Liu, H. (2019). Overlooked genetic diversity of ammonia oxidizing archaea lineages in the global oceans. J. Geophys. Res. Biogeosci. 124, 1799–1811.
Dai, M., Wang, L., Guo, X., Zhai, W., Li, Q., and He, B. (2008). Nitrification and inorganic nitrogen distribution in a large perturbed river/estuarine system: the Pearl River Estuary, China. Biogeosciences 5, 1227–1244. doi: 10.5194/bg-5-1227-2008
Damashek, J., Casciotti, K. L., and Francis, C. A. (2016). Variable nitrification rates across environmental gradients in turbid, nutrient-rich estuary waters of San Francisco bay. Estuar. Coast. 39, 1050–1071. doi: 10.1007/s12237-016-0071-7
Damashek, J., and Francis, C. A. (2018). Microbial nitrogen cycling in estuaries: from genes to ecosystem processes. Estuar. Coast. 41, 626–660. doi: 10.1007/s12237-017-0306-2
Damashek, J., Smith, J. M., Mosier, A. C., and Francis, C. A. (2015). Benthic ammonia oxidizers differ in community structure and biogeochemical potential across a riverine delta. Front. Microbiol. 5:743. doi: 10.3389/fmicb.2014.00743
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381
Francis, C. A., Roberts, K. J., Beman, J. M., Santoro, A. E., and Oakley, B. B. (2005). Ubiquity and diversity of ammonia-oxidizing archaea in water columns and sediments of the ocean. Proc. Natl. Acad. Sci. U.S.A. 102, 14683–14688. doi: 10.1073/pnas.0506625102
Frank, I. E., Turk-Kubo, K. A., and Zehr, J. P. (2016). Rapid annotation of nifH gene sequences using classification and regression trees facilitates environmental functional gene analysis. Environ. Microbiol. Rep. 8, 905–916. doi: 10.1111/1758-2229.12455
French, E., Kozlowski, J. A., Mukherjee, M., Bullerjahn, G., and Bollmann, A. (2012). Ecophysiological characterization of ammonia-oxidizing archaea and bacteria from freshwater. Appl. Environ. Microbiol. 78, 5773–5780. doi: 10.1128/aem.00432-12
Hall, G. H. (1986). “Nitrification in lakes,” in Nitrification, 1st Edn, ed. J. I. Prosser (Washington, DC: IRL Press), 127–156.
Hampel, J. J., McCarthy, M. J., Gardner, W. S., Zhang, L., Xu, H., Zhu, G., et al. (2018). Nitrification and ammonium dynamics in Taihu Lake, China: seasonal competition for ammonium between nitrifiers and Cyanobacteria. Biogeosciences 15, 733–748. doi: 10.5194/bg-15-733-2018
Hampel, J. J., McCarthy, M. J., Reed, M. H., and Newell, S. E. (2019). Short term effects of Hurricane Irma and cyanobacterial blooms on ammonium cycling along a freshwater–estuarine continuum in south Florida. Front. Mar. Sci. 6:640. doi: 10.3389/fmars.2019.00640
Happel, E., Bartl, I., Voss, M., and Riemann, L. (2018). Extensive nitrification and active ammonia oxidizers in two contrasting coastal systems of the Baltic Sea. Environ. Microbiol. 20, 2913–2926. doi: 10.1111/1462-2920.14293
Havens, K. E., Jin, K.-R., Rodusky, A. J., Sharfstein, B., Brady, M. A., East, T. L., et al. (2001). Hurricane effects on a shallow lake ecosystem and its response to a controlled manipulation of water level. Sci. World J. 1, 44–70. doi: 10.1100/tsw.2001.14
Hayden, C. J., and Beman, J. M. (2014). High abundances of potentially active ammonia-oxidizing bacteria and archaea in oligotrophic, high-altitude lakes of the Sierra Nevada, California, USA. PLoS One 9:e111560. doi: 10.1371/journal.pone.0111560
Heiss, E. M., and Fulweiler, R. W. (2016). Coastal water column ammonium and nitrite oxidation are decoupled in summer. Estuar. Coast. Shelf Sci. 178, 110–119. doi: 10.1016/j.ecss.2016.06.002
Herbold, C. W., Pelikan, C., Kuzyk, O., Hausmann, B., Angel, R., Berry, D., et al. (2015). A flexible and economical barcoding approach for highly multiplexed amplicon sequencing of diverse target genes. Front. Microbiol. 6:731. doi: 10.3389/fmicb.2015.00731
Hink, L., Gubry-Rangin, C., Nicol, G. W., and Prosser, J. I. (2018). The consequences of niche and physiological differentiation of archaeal and bacterial ammonia oxidisers for nitrous oxide emissions. ISME J. 12, 1084–1093. doi: 10.1038/s41396-017-0025-5
Horak, R. E., Qin, W., Schauer, A. J., Armbrust, E. V., Ingalls, A. E., Moffett, J. W., et al. (2013). Ammonia oxidation kinetics and temperature sensitivity of a natural marine community dominated by Archaea. ISME J. 7, 2023–2033.
Hou, J., Song, C., Cao, X., and Zhou, Y. (2013). Shifts between ammonia-oxidizing bacteria and archaea in relation to nitrification potential across trophic gradients in two large Chinese lakes (Lake Taihu and Lake Chaohu). Water Res. 47, 2285–2296. doi: 10.1016/j.watres.2013.01.042
Hsiao, S. S. Y., Hsu, T. C., Liu, J. W., Xie, X., Zhang, Y., Lin, J., et al. (2014). Nitrification and its oxygen consumption along the turbid Chang Jiang River plume. Biogeosciences 11, 2083–2098. doi: 10.5194/bg-11-2083-2014
Hugoni, M., Etien, S., Bourges, A., Lepère, C., Domaizon, I., Mallet, C., et al. (2013). Dynamics of ammonia-oxidizing Archaea and Bacteria in contrasted freshwater ecosystems. Res. Microbiol. 164, 360–370. doi: 10.1016/j.resmic.2013.01.004
James, R. T., Havens, K., Zhu, G., and Qin, B. (2009). Comparative analysis of nutrients, chlorophyll and transparency in two large shallow lakes (Lake Taihu, P.R. China and Lake Okeechobee, USA). Hydrobiologia 627, 211–231. doi: 10.1007/s10750-009-9729-5
James, T. R., Chimney, M. J., Sharfstein, B., Engstrom, D. R., Schottler, S. P., East, T., et al. (2008). Hurricane effects on a shallow lake ecosystem, Lake Okeechobee, Florida (USA). Fundam. Appl. Limnol. 172, 273–287. doi: 10.1127/1863-9135/2008/0172-0273
Jia, Z., and Conrad, R. (2009). Bacteria rather than Archaea dominate microbial ammonia oxidation in an agricultural soil. Environ. Microbiol. 11, 1658–1671. doi: 10.1111/j.1462-2920.2009.01891
Jiang, X., Wu, Y., Liu, G., Liu, W., and Lu, B. (2017). The effects of climate, catchment land use and local factors on the abundance and community structure of sediment ammonia-oxidizing microorganisms in Yangtze lakes. AMB Express 7:173.
Jochem, F. J., Lavrentyev, P. J., and First, M. R. (2004). Growth and grazing rates of bacteria groups with different apparent DNA content in the Gulf of Mexico. Mar. Biol. 145, 1213–1225. doi: 10.1007/s00227-004-1406-7
Kan, J. (2018). Storm events restructured bacterial community and their biogeochemical potentials. J. Geophys. Res. Biogeosci. 123, 2257–2267.
Kits, K. D., Sedlacek, C. J., Lebedeva, E. V., Han, P., Bulaev, A., Pjevac, P., et al. (2017). Kinetic analysis of a complete nitrifier reveals an oligotrophic lifestyle. Nature 549, 269–272. doi: 10.1038/nature23679
Könneke, M., Bernhard, A. E., José, R., Walker, C. B., Waterbury, J. B., and Stahl, D. A. (2005). Isolation of an autotrophic ammonia-oxidizing marine archaeon. Nature 437, 543–546. doi: 10.1038/nature03911
Kramer, B. J., Davis, T. W., Meyer, K. A., Rosen, B. H., Goleski, J. A., Dick, G. J., et al. (2018). Nitrogen limitation, toxin synthesis potential, and toxicity of cyanobacterial populations in Lake Okeechobee and the St. Lucie River Estuary, Florida, during the 2016 state of emergency event. PLoS One 13:e0196278. doi: 10.1371/journal.pone.0196278
Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549. doi: 10.1093/molbev/msy096
Laperriere, S. M., Nidzieko, N. J., Fox, R. J., Fisher, A. W., and Santoro, A. E. (2018). Observations of variable ammonia oxidation and nitrous oxide flux in a eutrophic estuary. Estuar. Coast. 42, 33–44. doi: 10.1007/s12237-018-0441-4
Lapointe, B. E., Herren, L. W., and Bedford, B. J. (2012). Effects of hurricanes, land use, and water management on nutrient and microbial pollution: St. Lucie estuary, Southeast Florida. J. Coast. Res. 285, 1345–1361. doi: 10.2112/jcoastres-d-12-00070.1
Lavrentyev, P. J., Gardner, W. S., and Johnson, J. R. (1997). Cascading trophic effects on aquatic nitrification: experimental evidence and potential implications. Aquat. Microb. Ecol. 13, 161–175. doi: 10.3354/ame013161
Letunic, I., and Bork, P. (2019). Interactive tree of life (iTOL) v4: recent updates and new developments. Nucleic Acids Res. 47, W256–W259.
Li, H., Weng, B. S., Huang, F. Y., Su, J. Q., and Yang, X. R. (2015). pH regulates ammonia-oxidizing bacteria and archaea in paddy soils in Southern China. Appl. Microbiol. Biotechnol. 99, 6113–6123. doi: 10.1007/s00253-015-6488-2
Li, M., Wei, G., Shi, W., Sun, Z., Li, H., Wang, X., et al. (2018). Distinct distribution patterns of ammonia-oxidizing archaea and bacteria in sediment and water column of the Yellow River estuary. Sci. Rep. 8:1584. doi: 10.1038/s41598-018-20044-6
Li, Y., Ding, K., Wen, X., Zhang, B., Shen, B., and Yang, Y. (2016). A novel ammonia-oxidizing archaeon from wastewater treatment plant: its enrichment, physiological and genomic characteristics. Sci. Rep. 6:23747.
Lipschultz, F., Wofsy, S. C., Ward, B. B., Codispoti, L. A., Friedrich, G., and Elkins, J. W. (1990). Bacterial transformations of inorganic nitrogen in the oxygen-deficient waters of the Eastern Tropical South Pacific Ocean. Deep Sea Res. A Oceanogr. Res. Pap. 37, 1513–1541. doi: 10.1016/0198-0149(90)90060-9
Liu, T. T., and Yang, H. (2019). An RNA-based quantitative and compositional study of ammonium-oxidizing bacteria and archaea in Lake Taihu, a eutrophic freshwater lake. FEMS Microbiol. Ecol. 95:fiz117.
Martens-Habbena, W., Berube, P. M., Urakawa, H., de la Torre, J. R., and Stahl, D. A. (2009). Ammonia oxidation kinetics determine niche separation of nitrifying Archaea and Bacteria. Nature 461, 976–979. doi: 10.1038/nature08465
Marton, J. M., Roberts, B. J., Bernhard, A. E., and Giblin, A. E. (2015). Spatial and temporal variability of nitrification potential and ammonia-oxidizer abundances in Louisiana salt marshes. Estuar. Coast. 38, 1824–1837. doi: 10.1007/s12237-015-9943-5
McIlvin, M. R., and Altabet, M. A. (2005). Chemical conversion of nitrate and nitrite to nitrous oxide for nitrogen and oxygen isotopic analysis in freshwater and seawater. Anal. Chem. 77, 5589–5595. doi: 10.1021/ac050528s
McMurdie, P. J., and Holmes, S. (2013). phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. doi: 10.1371/journal.pone.0061217
McMurdie, P. J., and Holmes, S. (2014). Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput. Biol. 10:e1003531. doi: 10.1371/journal.pcbi.1003531
Merbt, S. N., Auguet, J. C., Blesa, A., Martí, E., and Casamayor, E. O. (2015). Wastewater treatment plant effluents change abundance and composition of ammonia-oxidizing microorganisms in Mediterranean urban stream biofilms. Microb. Ecol. 69, 66–74. doi: 10.1007/s00248-014-0464-8
Mosier, A. C., and Francis, C. A. (2008). Relative abundance and diversity of ammonia-oxidizing archaea and bacteria in the San Francisco Bay estuary. Environ. Microbiol. 10, 3002–3016. doi: 10.1111/j.1462-2920.2008.01764.x
Newell, S. E., Babbin, A. R., Jayakumar, A., and Ward, B. B. (2011). Ammonia oxidation rates and nitrification in the Arabian Sea: Arabian Sea ammonia oxidation and nitrification. Glob. Biogeochem. Cycles 25:GB4016. doi: 10.1029/2010GB003940
Newell, S. E., Eveillard, D., McCarthy, M. J., Gardner, W. S., Liu, Z., and Ward, B. B. (2014). A shift in the archaeal nitrifier community in response to natural and anthropogenic disturbances in the northern Gulf of Mexico. Environ. Microbiol. Rep. 6, 106–112. doi: 10.1111/1758-2229.12114
Newell, S. E., Fawcett, S. E., and Ward, B. B. (2013). Depth distribution of ammonia oxidation rates and ammonia-oxidizer community composition in the Sargasso Sea. Limnol. Oceanogr. 58, 1491–1500. doi: 10.4319/lo.2013.58.4.149
Nicklisch, A., and Kohl, J. G. (1983). Growth kinetics of Microcystis aeruginosa (Kütz) Kütz as a basis for modelling its population dynamics. Int. Rev. Gesamten Hydrobiol. Hydrogr. 68, 317–326. doi: 10.1002/iroh.19830680304
Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2018). R Package ‘Vegan’ (Vienna). (accessed October, 2019). doi: 10.1002/iroh.19830680304
Paerl, H. W., Bales, J. D., Ausley, L. W., Buzzelli, C. P., Crowder, L. B., Eby, L. A., et al. (2001). Ecosystem impacts of three sequential hurricanes (Dennis, Floyd, and Irene) on the United States’ largest Lagoonal estuary, Pamlico Sound, NC. Proc. Natl. Acad. Sci. U.S.A. 98, 5655–5660. doi: 10.1073/pnas.101097398
Paerl, H. W., Havens, K. E., Hall, N. S., Otten, T. G., Zhu, M., Xu, H., et al. (2019). Mitigating a global expansion of toxic cyanobacterial blooms: confounding effects and challenges posed by climate change. Mar. Freshw. Res. 71, 579–592. doi: 10.1071/mf18392
Percuoco, V. P., Kalnejais, L. H., and Officer, L. V. (2015). Nutrient release from the sediments of the Great Bay Estuary, NH USA. Estuar. Coast. Shelf Sci. 161, 76–87. doi: 10.1016/j.ecss.2015.04.006
Phlips, E. J., Badylak, S., Hart, J., Haunert, D., Lockwood, J., O’Donnell, K., et al. (2012). Climatic influences on autochthonous and allochthonous phytoplankton blooms in a subtropical estuary, St. Lucie Estuary, Florida, USA. Estuar. Coast. 35, 335–352. doi: 10.1007/s12237-011-9442-2
Prosser, J. I., Head, I. M., and Stein, L. Y. (2014). “The family Nitrosomonadaceae,” in The Prokaryotes: Alphaproteobacteria and Betaproteobacteria, eds E. Rosenberg, E. F. DeLong, S. Lory, E. Stackebrandt, and F. Thompson (Berlin: Springer), 901–918. doi: 10.1007/978-3-642-30197-1_372
Prosser, J. I., and Nicol, G. W. (2012). Archaeal and bacterial ammonia-oxidisers in soil: the quest for niche specialisation and differentiation. Trends Microbiol. 20, 523–531. doi: 10.1016/j.tim.2012.08.001
Rotthauwe, J.-H., Witzel, K.-P., and Liesack, W. (1997). The ammonia monooxygenase structural gene amoA as a functional marker: molecular fine-scale analysis of natural ammonia-oxidizing populations. Appl. Environ. Microbiol. 63, 4704–4712. doi: 10.1128/aem.63.12.4704-4712.1997
Sanders, T., Schöl, A., and Dähnke, K. (2018). Hot spots of nitrification in the Elbe estuary and their impact on nitrate regeneration. Estuar. Coast. 41, 128–138. doi: 10.1007/s12237-017-0264-8
Santoro, A. E., Casciotti, K. L., and Francis, C. A. (2010). Activity, abundance and diversity of nitrifying archaea and bacteria in the central California Current. Environ. Microbiol. 12, 1989–2006. doi: 10.1111/j.1462-2920.2010.02205.x
Santoro, A. E., Francis, C. A., de Sieyes, N. R., and Boehm, A. B. (2008). Shifts in the relative abundance of ammonia-oxidizing bacteria and archaea across physicochemical gradients in a subterranean estuary. Environ. Microbiol. 10, 1068–1079. doi: 10.1111/j.1462-2920.2007.01547.x
Santoro, A. E., Saito, M. A., Goepfert, T. J., Lamborg, C. H., Dupont, C. L., and DiTullio, G. R. (2017). Thaumarchaeal ecotype distributions across the equatorial Pacific Ocean and their potential roles in nitrification and sinking flux attenuation. Limnol. Oceanogr. 62, 1984–2003. doi: 10.1002/lno.10547
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
Sedlacek, C. J., McGowan, B., Suwa, Y., Sayavedra-Soto, L., Laanbroek, H. J., Stein, L. Y., et al. (2019). A physiological and genomic comparison of Nitrosomonas cluster 6a and 7 ammonia-oxidizing bacteria. Microb. Ecol. 78, 985–994. doi: 10.1007/s00248-019-01378-8
Shade, A., Peter, H., Allison, S. D., Baho, D., Berga, M., Bürgmann, H., et al. (2012). Fundamentals of microbial community resistance and resilience. Front. Microbiol. 3:417. doi: 10.3389/fmicb.2012.00417
Shade, A., Read, J. S., Welkie, D. G., Kratz, T. K., Wu, C. H., and McMahon, K. D. (2011). Resistance, resilience and recovery: aquatic bacterial dynamics after water column disturbance. Environ. Microbiol. 13, 2752–2767. doi: 10.1111/j.1462-2920.2011.02546.x
Sintes, E., Bergauer, K., De Corte, D., Yokokawa, T., and Herndl, G. J. (2013). Archaeal amoA gene diversity points to distinct biogeography of ammonia-oxidizing Crenarchaeota in the ocean. Environ. Microbiol. 15, 1647–1658. doi: 10.1111/j.1462-2920.2012.02801.x
Small, G. E., Bullerjahn, G. S., Sterner, R. W., Beall, B. F. N., Brovold, S., Finlay, J. C., et al. (2013). Rates and controls of nitrification in a large oligotrophic lake. Limnol. Oceanogr. 58, 276–286. doi: 10.4319/lo.2013.58.1.0276
Stein, L. Y., Arp, D. J., Berube, P. M., Chain, P. S., Hauser, L., Jetten, M. S., et al. (2007). Whole-genome analysis of the ammonia-oxidizing bacterium, Nitrosomonas eutropha C91: implications for niche adaptation. Environ. Microbiol. 9, 2993–3007. doi: 10.1111/j.1462-2920.2007.01409.x
Tolar, B. B., Wallsgrove, N. J., Popp, B. N., and Hollibaugh, J. T. (2017). Oxidation of urea-derived nitrogen by Thaumarchaeota-dominated marine nitrifying communities: oxidation of urea-derived nitrogen in the sea. Environ. Microbiol. 19, 4838–4850. doi: 10.1111/1462-2920.13457
Urakawa, H., Martens-Habbena, W., Huguet, C., de la Torre, J. R., Ingalls, A. E., Devol, A. H., et al. (2014). Ammonia availability shapes the seasonal distribution and activity of archaeal and bacterial ammonia oxidizers in the Puget Sound Estuary. Limnol. Oceanogr. 59, 1321–1335. doi: 10.4319/lo.2014.59.4.1321
Verhamme, D. T., Prosser, J. I., and Nicol, G. W. (2011). Ammonia concentration determines differential growth of ammonia-oxidising archaea and bacteria in soil microcosms. ISME J. 5, 1067–1071. doi: 10.1038/ismej.2010.191
Vissers, E. W., Anselmetti, F. S., Bodelier, P. L., Muyzer, G., Schleper, C., Tourna, M., et al. (2013). Temporal and spatial coexistence of archaeal and bacterial amoA genes and gene transcripts in Lake Lucerne. Archaea 2013:289478.
Wang, Q., Quensen, J. F., Fish, J. A., Lee, T. K., Sun, Y., Tiedje, J. M., et al. (2013). Ecological patterns of nifH genes in four terrestrial climatic zones explored with targeted metagenomics using FrameBot, a new informatics tool. mBio 4:e00592-13.
Ward, B. B. (2008). “Nitrification in marine systems,” in Nitrogen in the Marine Environment, 2nd Edn, D. Capone, D. Bronk, M. Mulholland, and E. Carpenter (Amsterdam: Elsevier).
Wells, G. F., Park, H. D., Yeung, C. H., Eggleston, B., Francis, C. A., and Criddle, C. S. (2009). Ammonia-oxidizing communities in a highly aerated full-scale activated sludge bioreactor: betaproteobacterial dynamics and low relative abundance of Crenarchaea. Environ. Microbiol. 11, 2310–2328. doi: 10.1111/j.1462-2920.2009.01958.x
Wengrove, M. E., Foster, D. L., Kalnejais, L. H., Percuoco, V., and Lippmann, T. C. (2015). Field and laboratory observations of bed stress and associated nutrient release in a tidal estuary. Estuar. Coast. Shelf Sci. 161, 11–24. doi: 10.1016/j.ecss.2015.04.005
Wu, Y. J., Whang, L. M., Fukushima, T., and Chang, S. H. (2013). Responses of ammonia-oxidizing archaeal and betaproteobacterial populations to wastewater salinity in a full-scale municipal wastewater treatment plant. J. Biosci. Bioeng. 115, 424–432. doi: 10.1016/j.jbiosc.2012.11.005
Xia, H., Wu, Y., Chen, X., Huang, K., and Chen, J. (2019). Effects of antibiotic residuals in dewatered sludge on the behavior of ammonia oxidizers during vermicomposting maturation process. Chemosphere 218, 810–817. doi: 10.1016/j.chemosphere.2018.11.167
Yao, L., Chen, C., Liu, G., and Liu, W. (2018). Sediment nitrogen cycling rates and microbial abundance along a submerged vegetation gradient in a eutrophic lake. Sci. Total Environ. 616, 899–907. doi: 10.1016/j.scitotenv.2017.10.230
Zeng, J., Zhao, D.-Y., Huang, R., and Wu, Q. L. (2012). Abundance and community composition of ammonia-oxidizing archaea and bacteria in two different zones of Lake Taihu. Can. J. Microbiol. 58, 1018–1026. doi: 10.1139/w2012-078
Zhang, X., Zheng, S., Zhang, H., and Duan, S. (2018). Autotrophic and heterotrophic nitrification-anoxic denitrification dominated the anoxic/oxic sewage treatment process during optimization for higher loading rate and energy savings. Bioresour. Technol. 263, 84–93. doi: 10.1016/j.biortech.2018.04.113
Zhang, Y., Xie, X., Jiao, N., Hsiao, S. Y., and Kao, S. J. (2014). Diversity and distribution of amoA-type nitrifying and nirS-type denitrifying microbial communities in the Yangtze River estuary. Biogeosciences 11, 2131–2145. doi: 10.5194/bg-11-2131-2014
Zhao, D., Zeng, J., Wan, W., Liang, H., Huang, R., and Wu, Q. L. (2013). Vertical distribution of ammonia-oxidizing archaea and bacteria in sediments of a eutrophic lake. Curr. Microbiol. 67, 327–332. doi: 10.1007/s00284-013-0369-7
Zhao, R., Dahle, H., Ramírez, G. A., and Jørgensen, S. L. (2020). Indigenous ammonia-oxidizing archaea in oxic subseafloor oceanic crust. Msystems 5:e00758-19.
Zhong, W., Bian, B., Gao, N., Min, J., Shi, W., Lin, X., et al. (2016). Nitrogen fertilization induced changes in ammonia oxidation are attributable mostly to bacteria rather than archaea in greenhouse-based high N input vegetable soil. Soil Biol. Biochem. 93, 150–159. doi: 10.1016/j.soilbio.2015.11.003
Keywords: AOA, AOB, eutrophication, harmful algal blooms, hurricane, nitrification
Citation: Hampel JJ, McCarthy MJ, Aalto SL and Newell SE (2020) Hurricane Disturbance Stimulated Nitrification and Altered Ammonia Oxidizer Community Structure in Lake Okeechobee and St. Lucie Estuary (Florida). Front. Microbiol. 11:1541. doi: 10.3389/fmicb.2020.01541
Received: 21 April 2020; Accepted: 12 June 2020;
Published: 10 July 2020.
Edited by:
Sebastian Lücker, Radboud University Nijmegen, NetherlandsReviewed by:
Huiluo Cao, The University of Hong Kong, Hong KongHidetoshi Urakawa, Florida Gulf Coast University, United States
Copyright © 2020 Hampel, McCarthy, Aalto and Newell. 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: Justyna J. Hampel, justyna.hampel@usm.edu