Original Research ARTICLE
Mud, Microbes, and Macrofauna: Seasonal Dynamics of the Iron Biogeochemical Cycle in an Intertidal Mudflat
- 1Bigelow Laboratory for Ocean Sciences, Boothbay, ME, United States
- 2San Francisco Bay Regional Water Quality Control Board, Oakland, CA, United States
- 3Department of Earth and Planetary Sciences, Harvard University, Cambridge, MA, United States
- 4Department of Organismal and Evolutionary Biology, Harvard University, Cambridge, MA, United States
Microorganisms and burrowing animals exert a pronounced impact on the cycling of redox sensitive metals in coastal sediments. We hypothesized that the iron biogeochemical cycle and associated sedimentary microbial community will respond to seasonal changes in a bioturbated intertidal mudflat. In this study, we monitored the spatiotemporal dynamics of porewater and highly reactive solid phase iron with the corresponding prokaryotic and eukaryotic sedimentary microbial communities over one annual cycle from November 2015 to November 2016. Continuous and seasonally variable pools of both porewater Fe(II) and highly reactive iron (FeHR) were observed throughout the seasons with significant increases of Fe(II) and FeHR in response to increased sediment temperature in summer months. Maximum concentrations of Fe(II) and FeHR were predominantly confined to the upper 5 cm of sediment throughout the year. Iron-oxidizing and -reducing microorganisms were present and stable temporally, and exhibited a depth-dependent stratification likely due to availability of Fe(II) and FeHR pools, respectively. Zetaproteobacteria, presumptive lithotrophic iron-oxidizing bacteria, were present at abundances around 0.5–1% in the top 5 cm of sediment with decline abundance with depth. As a whole the microbial community was relatively stable across the seasons, and showed strongest separation with depth, probably driven by changes in oxygen availability and organic matter. The Deltaproteobacteria, principally taxa known to be associated with sulfur and iron cycling, were prevalent, especially at >5 cm depth. Gammaproteobacteria and Bacteroidetes were also abundant, with putatively aerobic members especially prevalent in the cm of the sediment. The relative abundance of diatoms, estimated from abundance of 18S rRNA gene counts, showed evidence of a seasonal signal possibly tied to spring and fall blooms. Overall, analysis of phytoplankton found significant abundance at depth, likely due to the feeding and bio-mixing activity of marine worms. Macro-, and meiofauna, consistent with expected taxa, were detected throughout the year via 18S gene counts, and showed some seasonal variations that may influence sedimentary iron transformations by active microbial grazing. In summary, this analysis revealed relatively consistent temporal and spatial trends in iron geochemistry and microbial and macrobial community composition, while also indicating a complex dynamic of microbial and macrobial interactions are responsible for maintaining these processes.
Coastal marine sediments are active and dynamic environments that connect the flux of nutrients and energy from weathered land masses to marine ecosystems. Continentally-derived energy and nutrient sources are utilized and recycled by pelagic and benthic organisms, which act as key biogeochemical agents in coastal ecosystems. Biogeochemical cycles in coastal sedimentary ecosystems can fluctuate on hourly to seasonal time scales, and represent dynamic gradients in space and time influenced by a myriad of biogeophysical processes such as riverine discharge, tidal cycles, seasonal stratification, animal behavior, and microbial activity (Kristensen and Kostka, 2005; Meysman et al., 2006; Bouma et al., 2009). In coastal sediments, sediment-dwelling animals have played a primary top-down control on biogeochemical processes over the last ∼500 million years with important repercussions for the entire earth system (Thayer, 1979; Meysman et al., 2006; Tarhan et al., 2015). These ecosystem-engineering animals actively burrow in sediments and irrigate their burrows with the overlaying bottom water, supplying oxygen and nutrients deep into oxygen-free (i.e., reduced) sediments. In addition to bio-irrigation, burrowing animals, principally marine worms in this case, physically mix sediments, which can act to transport organic matter and minerals within sediment layers. Their burrow walls are hotspots for microbial redox cycling where steep geochemical gradients drive the rapid cycling between reduced and oxidized chemical species (Fenchel, 1969; Reichardt, 1986; Aller, 1994; Fenchel, 1996; Brune et al., 2000; Kristensen and Kostka, 2005; Bertics and Ziebis, 2009). In addition, their complex burrowing and irrigating and bio-mixing behavior essentially increases the exchange of solutes between benthic and pelagic systems (Aller, 1982), and are thus of primary importance to benthic-pelagic coupling and ecosystem function (Griffiths et al., 2017).
Iron is an essential micronutrient for life, and where concentrations of redox active forms or iron are high, such as sedimentary environments, it is an important electron donor and acceptor in microbial metabolism (Emerson et al., 2010; Melton et al., 2014). Numerous, interdependent biological, chemical, and physical processes control the redox cycling of iron in sedimentary environments. Iron is transported to coastal ecosystems by rivers where it settles in sediments and is intensely [re]cycled (Canfield, 1989; Canfield et al., 1993). Sediment-dwelling animals act in concert with diverse microbiota in coastal sediments to regulate iron biogeochemical transformations (Thibault de Chanvalon et al., 2016; Beam et al., 2018; van de Velde et al., 2020). These macrobiota-microbiota relationships result in labile pools of dissolved and highly reactive iron that contribute to important ecosystem processes such as formation of bioavailable iron and subsequent transport to the water column (Elrod et al., 2004; Severmann et al., 2010) and inhibition of porewater hydrogen sulfide accumulation (Antler et al., 2019; van de Velde et al., 2020). These labile pools of iron can also provide a continuous replenishment of reduced and oxidized iron for iron-oxidizing and iron-reducing microorganisms, respectively, (Canfield, 1989; Canfield et al., 1993; Beam et al., 2018). In addition to these direct biological factors, the sedimentary iron biogeochemical cycle may also respond to fluctuations in temperature and organic matter transport from primary productivity in the water column. Based on this knowledge, we have a reasonable understanding of some of the major drivers of coastal iron biogeochemistry. However, less is known about how depth-resolved seasonal changes sedimentary iron biogeochemical cycle changes impact the sedimentary, particularly with regard to the composition and function of sedimentary microbial communities. In many ecological systems, sub-annual seasonal dynamics are as important or can even outweigh year-to-year variation (Staudinger et al., 2019). Thus, it is essential to understand the seasonal dynamics of these interconnected processes, in order to predict how they might respond to effects of anthropogenic climate change in coastal sedimentary ecosystems such as increased temperature, declines in oxygen, and eutrophication (Doney, 2010).
In this study, we tracked the spatial and temporal dynamics of the iron biogeochemical cycle over one seasonal cycle (2015–2016) in a bioturbated intertidal mudflat to determine changes in the reduced porewater Fe(II) and highly reactive Fe (FeHR) mineral pools. The shifts in the sedimentary prokaryotic and eukaryotic microbial communities were also spatiotemporally tracked with 16S and 18S rRNA gene sequencing, respectively. These data suggested a relatively continuous and dynamic pool of both Fe(II) and FeHR over the year, and statistically significant increases of both pools with increasing temperature. The microbial communities were relatively stable over the annual cycle with some pronounced seasonal and spatial variability, possibly due to interactions with water column processes such as phytoplankton blooms. Seasonal and spatial shifts in the sedimentary pools of porewater and solid phase iron are controlled by the intricate interaction between sediment-dwelling animals, sedimentary microbial communities, and water column processes that are a result of complex interconnected processes.
Materials and Methods
An intertidal mudflat (the Eddy, 43.994827, and −69.648632) located within the Sheepscot River Estuary Maine, United States of America was the subject of the current study (Figure 1; Stickney, 1959). As defined by Stickney the Eddy is classified as an upper estuary mudflat with seasonal variations in salinity and temperature that receives a significant influence from brackish and/or freshwater riverine input (Stickney, 1959). From November 2015 to November 2016 sediment cores were sampled (n = 11); for porewater and solid phase geochemistry in combination with microbial community structure (n = 10) at 1 cm depth intervals up to a depth of 10 cm. Due to circumstances beyond our control cores were not collected in August and October 2016. We took care not to take monthly samples in areas disturbed by our previous sampling events. This site was not impacted by local worm and shellfish harvesters as evinced by a lack of visible digging tracks in the surface sediments (Figure 1A), or direct observations of their presence. A wooden post and a washed-up bundle of lobster traps served as a reference point for sampling throughout the study (Figure 1B). At high tide, this area of the mudflat was covered with approximately one meter of water, dependent on tidal amplitude. At low tide, we successively took samples in the same horizontal distance from shore to reduce possible longitudinal effects on the mudflat macrofaunal populations. This mudflat lies adjacent to a tidally flushed salt marsh to the east, and faces south-southwest to the main stem of the Sheepscot River. The salinity and temperature varied throughout the year, but the overlying water was never frozen in the winter months of 2015–2016. Salinity, temperature and spatial/temporal data are presented in Supplementary Table 1.
Figure 1. A photomontage of images over different time and spatial scales from seasonal sampling at “The Eddy” (Maine, United States; 43.994827, −69.648632). A landscape view of the mudflat at low tide (A; 26 May 2016) with dashed arrow line depicting the progression of sampling throughout the annual cycle (tracks created from sampling are visible at the position distal to the arrowhead). (B,C) A representative core (18 November 2015) taken from “The Eddy” at the site of the white arrow showing ∼20 cm of total sediment showing deep mixing by macrofauna up to ∼15–17 cm in depth by visible iron oxides, then the sediments become bluish-gray to black. A one-centimeter core slice from 5 cm depth (D) showing characteristic iron oxide lined worm burrows most likely created by the common polychaete Nereis diversicolor and/or the intertidal inhabiting hemichordate Saccoglossus kowalevskii. Close up image of a dissected iron oxide lined burrow (E; 28 December 2015) and light microscope image of burrow lining at 400x magnification (F). Fluorescence microscope image of DNA stained microorganisms inhabiting the iron oxide burrow lining (G; 8 October 2015). Filamentous Beggiatoa spp. (thick filament) And putative “cable bacteria” (thin filaments) were commonly observed throughout the year.
Porewater Fe(II) Geochemistry
Sediment cores (approximately 15–20 cm in depth) from the Eddy were taken with a Plexiglas core barrel (inner diameter, 7.5 cm) that had predrilled 4 mm holes drilled (sealed with electrical tape before sampling) at 1 cm intervals for Rhizon CSS samplers with a 5 cm length and 0.15 μm pore size (Rhizosphere Research Products, Wageningen, Netherlands; Seeberg-Elverfeldt et al., 2005). Porewaters were extracted by inserting Rhizons into the predrilled holes starting from the top of the core and working down to the final depth interval (9–10 cm). Sterile 10 mL syringes were used to pull negative pressure on the Rhizons, and the plunger was held in place for 1–2 h with a wooden block. After ∼10 mL of porewater had been extracted, the Rhizons were removed, and the porewater was dispensed into a sterile 15 mL tube in a glove bag containing an atmosphere of N2/H2 (95/5%) or under ambient atmospheric conditions. Ferrous iron [Fe(II)] was immediately preserved by pipetting 250 μL of pore water into 250 μL of Ferrozine buffer (10 mM Ferrozine in 50 mM HEPES buffer, pH = 7; Stookey, 1970; McBeth et al., 2011). We did not observe substantial differences in porewater Fe(II) concentrations when rhizon-derived porewater was rapidly preserved, <10 s from rhizon collection to Ferrozine treatment under ambient air conditions, compared to working in the glove bag, thus after the April sampling all samples were collected under ambient atmospheric conditions. The Ferrozine-porewater solution was measured at 562 nm on a Multiskan MCC plate reader (Thermo Fisher) and absorbance values were calculated to concentrations by plotting them on a standard curve of known Fe(II) concentration solutions ranging from 0 to 1,000 μM prepared from a 10 mM solution of FeCl2 in 0.1 N HCl.
Solid Phase Highly Reactive Fe (FeHR) Geochemistry
Sediment cores used for porewater extraction were extruded and sliced into 1 cm sections with a brass wire and dried in an oven at 70°C for 24 h. Highly reactive iron oxides (ferrihydrite and lepidocrocite as well as iron monosulfides that were oxidized during this process) were extracted for 48 h in a 50 mL falcon tube with 10 mL of hydroxylamine-HCl (1 M) in 25% acetic acid (w/v; Poulton and Canfield, 2005). The total extracted Fe was measured with a Ferrozine solution as described above, except the samples were diluted 1:100 in distilled water due to high absorbance values of undiluted samples.
Sediment DNA Extraction
Sub-cores were sampled simultaneously with porewater extraction using 1 mL syringes with cut off tips. The 1 mL syringe corers were inserted into pre-drilled 8 mm holes at a right angle to the Rhizon sampling holes, and gently pressed into the side of the core while drawing back on the plunger to remove approximately 0.2–0.5 cc of sediment. The 1 mL sub-cores were then frozen at −80°C until DNA extraction. DNA was extracted from ∼0.1–0.2 g of wet sediment using the DNeasy PowerSoil Kit (Qiagen) with a modification of the removal of 200 μL of the bead solution and addition of 200 μL of 24:25:1 phenol:chloroform:isoamyl alcohol (Sigma Aldrich); the manufacturers protocol was then followed for the rest of the extraction, with the exception of the DNA column elution step, which was eluted with 100 μL of molecular grade water instead of the elution buffer from the kit. The DNA was quantified using a Qubit 3.0 fluorometer with the Qubit dsDNA quantification kit and protocol (Thermo Fisher Scientific).
16S and 18 rRNA Gene Sequencing and Analysis
The 16S rRNA gene was amplified from the 1 cm interval DNA extractions with the 515F/926R primer pair and the 18S rRNA gene was amplified with the E572F/E1009R primer pair at the Integrated Microbiome Resource (Dalhousie University, Halifax, NS, Canada). Paired end reads were assembled and classified using a mothur (Schloss et al., 2009) pipeline from the mothur MiSeq standard operating procedure (Kozich et al., 2013; Scott et al., 2017). Operational taxonomic units (OTUs) were clustered at either 85% or 97%, which generated 12,437 and 351,057 OTUs, respectively. The 85% identity cutoff, which corresponds roughly to the class-order taxonomic level (Yarza et al., 2014), was chosen for comparative analysis to limit the number of clusters generated for comparisons and lend clarity to data representations. While the 97% cut-off did sub-divide some clusters among the top 100 OTUs (based on relative abundance), we did not note any substantial differences that would change the functional inferences made from OTU analysis. Relative abundances were utilized to build a heatmap in RStudio (RStudio Team) of the top OTUs that were greater than 0.5% abundance in each sample.
The 18S rRNA gene data was processed using mothur (v1.43) following the mothur MiSeq SOP data analysis pipeline (Kozich et al., 2013), but referencing the pR2 database (v4.12) to assign eukaryote taxonomies (Guillou et al., 2013). Eukaryote OTUs were determined at a level of 97% sequence similarity (Hu et al., 2015; Pasulka et al., 2019), which is less conservative than some thresholds proposed for rRNA genes (Caron et al., 2009), and more conservative than others (Callahan et al., 2017; Edgar, 2018). Counts of eukaryotic lineages of interest (i.e., diatoms, polychaetes, molluscs, nematodes, and protozoa) were filtered from the taxonomic data and used to calculate relative abundances for all samples across the year.
One way to organize and find signals in complex ecological communities is to cluster them into coherent “ecostates” and then to map these ecostates back onto time and space to highlight patterns (Record et al., 2017). Each ecostate represents a consortium of taxa that clusters together consistently throughout the sampling program, and maintains the same respective relative abundances of taxa (O’Brien et al., 2016). For this analysis, samples were clustered using average linkage hierarchical clustering of the Bray–Curtis dissimilarity as implemented in the R package vegan (Oksanen et al., 2019). The clustering algorithm used bacterial and archaeal OTUs clustered at 85% identity (n = 12,437 OTUs), since, as discussed above, this identity was a compromise between capturing overall community diversity without overwhelming the analysis with a large number of OTUs. Ecostates were then numbered and mapped back onto the spatiotemporal sampling grid in R Studio (R Studio Team).
The pore water and solid phase sedimentary geochemical data can be found in the Biological and Chemical Oceanography Data Management Office (BCO-DMO) database under the DOI:10.1575/1912/bco-dmo.737962.1. The 16S and 18S rRNA gene raw read data can be found under the Bioproject ID PRJEB36098. All sedimentary geochemical and DNA data can also be found in Supplementary Table 1; all geochemical and DNA profile plots in the main text can be reproduced from this data set.
Results and Discussion
Seasonal Porewater Dissolved Ferrous Iron [Fe(II)] Pools
The dissolved ferrous iron [Fe(II)] in porewaters exhibited both spatial and temporal variability over the sampling period, which was approximately every month during the 2015–2016 annual cycle (Figures 2A,B). Porewater Fe(II) concentrations peaked at a maximum average concentration of approximately 120 μmol L–1 at a depth interval between 3–4 cm, which is a consistent observation with many porewater Fe(II) profiles from coastal and continental shelf systems (Aller, 1980; Canfield, 1989; Hines et al., 1991; Thamdrup et al., 1994; Severmann et al., 2010; Antler et al., 2019), which are more often than not, bioturbated. This peak in porewater Fe(II) is consistent with the activity of iron-reducing microorganisms at this depth interval throughout the study. Throughout the year, the porewater Fe(II) concentration maxima remained in the top 5 cm of sediment, and migrated only once to a depth of 6.5 cm in September (Figure 2B). These Fe(II) porewater maxima suggest that the most active iron cycling zone occurs within the upper 5 cm of sediment. Concentrations of Fe(II) decreased to approximately a seasonal average of 100 μmol L–1 at the maximum depth interval sampled between 9–10 cm (Figure 2A). The decline in Fe(II) with depth is likely due to the reaction with hydrogen sulfide produced by sulfate-reducing bacteria in sediments (Jørgensen, 1982) as well as around macrofaunal burrow walls (Bertics and Ziebis, 2010), which results in the precipitation of iron sulfides in sediments (Berner, 1984). There were consistent pools of porewater Fe(II) over the annual cycle, which exhibited large spatiotemporal changes (Figure 2B). The large variability over space and time could be due, in part, to the heterogeneity of natural populations of macrofauna and microorganisms in the mudflat, which we attempted to minimize by sampling cores within a consistent region over the course of the study (Figure 1). Other studies have shown consistent variability in porewater Fe(II) with space and time (Hines et al., 1991; Thamdrup et al., 1994; van de Velde et al., 2020), while another recent study showed contemporaneously collected salt marsh sediment samples show consistency in geochemical profiles despite large geographic variation (Antler et al., 2019). The abundance and activity of macrofaunal populations likely exerts a primary control on porewater Fe(II) concentrations, due to their sediment mixing, irrigating, and feeding behavior (Aller, 1994; van de Velde et al., 2020).
Figure 2. Porewater Fe(II), highly reactive Fe concentrations, and sedimentary DNA profiles from averages across depths and season (A,C,E) and over the annual seasonal cycle from Nov. 2015 to Nov. 2016 (B,D,F). Blue lines represent mean concentrations calculated from the individual seasonal data points (dark gray circles) and light gray shaded area is the 95% confidence interval.
Generally, the porewater Fe(II) concentrations increased with rising temperature (Figure 3A). This might be indicative of elevated sedimentary carbon mineralization in summer months by iron-reducing bacteria (Canfield, 1989; Thamdrup et al., 1994). Sulfide production via sulfate reduction could also lead to iron reduction (Luther and Church, 1988); however, we did not detect hydrogen sulfide in these sediments, indicating anaerobic respiration coupled to iron reduction is, on balance, an important process. The increases in porewater Fe(II) with elevated temperatures between the months of February and July has also been observed in bioturbated sediments in the Great Bay New Hampshire, United States (Hines et al., 1984). The supply of terrigenous organic matter to this mudflat occurs mainly in the spring that corresponds with the thawing land mass and riverine run-off (Stickney, 1959). Additional organic inputs come from settling of detrital material from the water column, and benthic mats of phytoplankton that can be transported into the sediments via the bio-irrigating and bio-mixing activites of macrofauna (Christensen et al., 2000). It is reasonable that this accounts in part for the high mobility and turnover of porewater Fe(II) due to the behavior of irrigating animals influencing microbial iron cycling processes (van de Velde et al., 2020). This relationship likely explains the large variance of the porewater Fe(II) pools over the 13 month long study, which is overall influenced by seasonal variations in sediment temperature (Figure 3A) and perhaps availability of macrofaunal food sources such as phytoplankton (discussed in detail below).
Figure 3. Porewater Fe(II; A), highly-reactive Fe (B), and sedimentary DNA (C) pools as a function of temperature over the annual cycle at “The Eddy.” Individual data points are gray for each sampling event over the year. Black and blue lines are a loess smooth and a generalized linear model fit, respectively, and gray shading represents the 95% confidence interval for each line fit. The linear model statistics and equation are in blue text.
Seasonal Highly Reactive Fe (FeHR) Pools
Similar to measurements of Fe(II), there was a continuous and variable pool of highly reactive Fe (FeHR) with respect to depth over the seasonal cycle (Figures 2C,D). The concentration of FeHR declined slightly with depth (Figure 2C) perhaps due to the interaction with hydrogen sulfide. The most significant change in the FeHR pool was in the warmer months (July and September) into the summer and early fall, with a nearly two-fold increase in the FeHR pools compared to the other months over the year (Figure 2D). We hypothesize that this change is due to several key biogeochemical processes, which are controlled by the master variable, temperature, over the seasonal cycle. Changes in irrigation rates that likely occur throughout the year (Fang et al., 2019) will have profound effects on the exchange of solutes as well as the control of the pools of sedimentary FeHR—when irrigational activity increases, concomitant microbial oxygen consuming process, such as microaerobic bacterial iron oxidation and aerobic carbon mineralization increase. When irrigation activity is low, the pools of sedimentary FeHR are depleted in the winter months. This is likely due to interaction Fe(II) with hydrogen sulfide produced by sulfate reducing microorganisms, resulting in the formation of more iron sulfides. An increase in the pool of FeHR in the summer has important biogeochemical ramifications, including: the potential to enhance the iron flux from the sediment by stimulating iron reduction (Severmann et al., 2010), and by providing fresh and reactive substrates for microbial iron reduction, which may effectively compete with sulfate reduction in bioturbated sediments (Froelich et al., 1979; Canfield, 1989).
Seasonal Sedimentary DNA Pools
The extractable quantity of sedimentary DNA (μg g–1 sediment) decreased from approximately 10 μg g–1 at 1 cm to ∼7.5 μg g–1 at 10 cm (Figure 2E). Assuming that DNA pools are approximate to microbial biomass, these values infer that microbial abundances decrease with depth, which is consistent with observations in both shallow and deep-sea sediments (Kallmeyer et al., 2012; Chen et al., 2017). Higher microbial biomass near the sediment-water interface is consistent with this being the zone where the most labile organic matter is deposited at the sediment surface, while organic matter becomes more refractory with depth. However, macrofauna may confound these interactions by actively exchanging and flushing bottom water into deeper sediments, as well as by bio-mixing, thereby introducing more labile organic matter, either dissolved or particulate, deep into the sediments. Across the year, DNA pools were distinct between fall and winter months versus summer months (Figure 2F). These variations may be caused by fall phytoplankton blooms that are stored in the sediments over the winter months, and decrease in the summer months due to the activity of microbial mineralization of detrital organic matter, including degradation and metabolism of exogenous DNA. This is evident when sedimentary DNA pools are plotted as a function of temperature, which is negatively correlated; as temperature increased, DNA pools decreased (Figure 3C).
Iron Cycling Microorganisms
There were OTUs that could be assigned to the iron biogeochemical cycle. The iron-oxidizing Zetaproteobacteria formed an OTU (Otu00040) that was 0.28% median seasonal abundance (range 0.01–2%), which is a typical relative abundance for this class of bacteria in marine sediments (Rubin-Blum et al., 2014; McAllister et al., 2015; Laufer et al., 2016, 2017; Beam et al., 2018). The Zetaproteobacteria relative abundances changed with respect to depth (Figure 4A), but did not exhibit pronounced seasonal variation except in one sample in September that appeared to correlate with increased FeHR production (Figures 2D, 4B). The relative abundance of Zetaproteobacteria ranges from 0.25–0.75%, in the upper 5 cm of sediment compared to the deeper sediment, with ranges from around 0–0.25%. The highest relative abundances were at 3–4 cm (Figure 4A), which has been previously documented from a depth profile of 16S rRNA genes from sediment and iron oxide macrofaunal burrow walls at this mudflat (Beam et al., 2018). The peak abundance of Zetaproteobacteria also coincides with the peak in porewater Fe(II) and FeHR from these sediments (Figures 2A,C). These Zetaproteobacteria profiles are consistent with them being most abundant in the active iron cycling zone throughout the year. It is also important to note that in this study we analyzed bulk sediment where we would expect lower overall numbers of Zetaproteobacteria, Previous studies, including one from this same site, have shown their abundance increases significantly in iron-oxide coated worm burrow walls, where they are likely involved in Fe(II) oxidation coupled to oxygen reduction and carbon dioxide fixation (McAllister et al., 2015; Beam et al., 2018). Consistent with these observations is that increased dark CO2 fixation with concomitant Fe(II) oxidation has been observed in burrow walls of the polychaete Neries diversicolor (Reichardt, 1986), which was the most abundant polychaete at the Eddy. Thus, dark CO2 fixation by iron-oxidizing Zetaproteobacteria may contribute to the autochthonous organic matter pools in coastal sediments, similar to that observed in autotrophic sulfur-oxidizing Gammaproteobacteria (Dyksma et al., 2016).
Figure 4. Relative abundance of iron cycling operational taxonomic units (OTUs) identified over the season: Iron-oxidizing Zetaproteobacteria (OTU00040) relative abundance as a function of depth (A) and season (B), and putative iron-reducing Desulfuromondaceae (OTU00028) as a function of depth (C) and season (D).
Microorganisms that catalyze the reduction of iron are phylogenetically diverse (Melton et al., 2014) and, particularly in marine sediments, are more difficult to identify with 16S rRNA gene profiling. However, culture-based studies from marine sediments have suggested that members of the Deltaproteobacteria family Desulfuromonadaceae can reduce iron, as well as other some other Deltaproteobacteria lineages (Vandieken et al., 2006), which have been subsequently identified in coastal sediments (Reyes et al., 2017; Otte et al., 2018; Buongiorno et al., 2019). We identified one OTU containing members of the Desulfuromonadaceae family (Otu00028) that was relatively abundant seasonally accounting for a median relative OTU abundance of 0.80% (range = 0.40–1.7%; Figures 4C,D). The Desulfuromonadaceae OTU was also most abundant in the active iron cycling zone (Figure 6C), which appears to mirror the porewater and highly reactive iron profiles (Figures 2A,C). The Desulfuromonadaceae abundances fluctuated temporally (Figure 4D) perhaps in response to organic matter and/or iron oxide availability. Correlations of Desulfuromonadaceae with porewater and highly reactive iron geochemistry were not indicative of a strong relationship (not shown), but this does not preclude a causal relationship. Canonical iron-reducing bacterial genera such as Geobacter and Shewanella-like OTUs were extremely rare or below detection in these sediments, similar to findings from coastal sediments near Denmark and Sweden (Reyes et al., 2017; Otte et al., 2018; Buongiorno et al., 2019).
Seasonal Dynamics of the Sedimentary Microbial Community
The sedimentary microbial communities were dominated by bacteria (total OTU abundance = 97.5%), with minor contributions by archaea (total OTU abundance = 2.5%)—consistent with the observations of Chen et al. (2017) and Kim et al. (2008) that bacteria dominate over archaea in bioturbated sediments. Two of the three most abundant bacterial OTUs throughout the seasons were related to the Deltaproteobacteria, and one was associated with the Gammaproteobacteria (Figure 5). Deltaproteobacteria and Gammaproteobacteria are commonly observed as abundant members of coastal sediment microbial communities (Kim et al., 2008; Lasher et al., 2009; Wang et al., 2012; Chen et al., 2017; Cleary et al., 2017; Qiao et al., 2018). The two Deltaproteobacteria OTUs were related to the families Desulfobulbaceae (Otu00001; median seasonal abundance = 6.4%; range = 3.0–8.8%) and the genus Desulfococcus in the family Desulfobacteraceae (Otu00002; median seasonal abundance = 5.7%; range = 1.6–9.8%), which are likely involved in sulfur oxidation and sulfate reduction, respectively, in these sediments. The Desulfobulbaceae is a diverse family within the Deltaproteobacteria that contains the so-called “cable bacteria” that transport electrons from sulfur oxidation over centimeter scale distances and catalyze oxygen reduction at sediment-water interfaces (Nielsen et al., 2010; Meysman et al., 2015; Nielsen and Risgaard-Petersen, 2015). We have observed on many occasions the presence of cable bacteria-like filaments from this field site (Figure 1F), and hypothesize that they are functional in these bioturbated sediments, which has been recently demonstrated (Aller et al., 2019), and may impact iron dynamics (Sulu-Gambari et al., 2016). Members of the Desulfobulbaceae are also involved in sulfur disproportionation, which could also be an important process in these sediments, and would impact the iron cycle (Finster, 2011).
Figure 5. A heatmap of the relative abundance of bacterial and archaeal taxonomic units (OTUs) clustered at 85% identity that were at least 0.5% relative abundance in one sample (n = 76 OTUs). The different samples were clustered using average linkage hierarchical clustering of the Bray-Curtis dissimilarity as implemented in the R package vegan (Oksanen et al., 2019). Solid colors depict the different ecostates that were selected from the clustering analysis (see legend in Figure 6). OTU identifiers can be referenced with specific taxonomic groups in Supplementary Table 1.
To gain a more holistic understanding of the microbial community dynamics, each discrete spatiotemporal sample (n = 100) was clustered based on the relative abundance of OTUs at 85% 16S rRNA gene identity. The composition of the sedimentary microbial community did not exhibit large fluctuations with the changing seasons in this bioturbated mudflat (Figures 5, 6). Although some individual OTUs exhibit spatial and temporal variability, the community was generally comprised of approximately 20–25 dominant OTUs (Figure 5). Seven main clusters predominated the various communities with respect to time and space (Figure 5). These clusters comprise different “ecostates” that were identified at various depths and times during the season. Ecostates are operationally defined as a collection of inferred taxonomic groups from distinct ecological samples, and they can be used to track changes in ecological communities through space and time (O’Brien et al., 2016). There is a strong depth stratification of ecostates that persists temporally (Figure 6). Ecostate 2 typically dominates all depth intervals <5 cm at all months. We speculate that this ecostate, which has greater abundances of putatively aerobic microbes (Supplementary Table 3), experiences more oxic conditions, partly due to diffusion of O2 from the sediment air/water interface, but primarily driven by biodiffusive transport by animals that create sporadic “oxygen oasis” in these otherwise anoxic sediments (Glud, 2008; Volkenborn et al., 2010). Ecostate 1 was usually confined to deeper sediments >5 cm where lower oxygen conditions would be prevalent (Glud, 2008; Figure 6). The presumed presence of labile organic matter in the upper 5 cm of sediment also likely contributes to the separation of ecostates 1 and 2 as has been observed between other surface and subsurface coastal sediments (Cleary et al., 2017; Qiao et al., 2018). These ecostates are relatively stable throughout the seasons (Figure 6), and appear to be unaffected by observed seasonal changes in environmental conditions and other stochastic events likely occurring in the intertidal zone. More than half of the samples at the sediment-water interface were classified as ecostate 3, which suggests some unique properties of this interface that contributes to the observed community signal (Figures 5, 6). This consistent pattern of stratification was only disrupted once, in the April profile. This profile differed markedly, with three unique ecostates (4, 5, and 6) that were not observed in any other month. These ecostates contained a cluster of high relative abundances of OTUs, described below, that only comprise a minor community component of the other ecostates (Figure 6). The occurrence of unique ecostates in April could reflect a community response to ecosystem disturbance or a stochastic assembly of horizontal spatial heterogeneity, which would both imply resilience of the system.
Figure 6. Ecostates (n = 7) of bacterial and archaeal microbial communities across one seasonal cycle. The seven different ecostates were identified from the Bray-Curtis dissimilarity matrix of the relative abundance of all OTUs (n = 12,437) by average linkage hierarchical clustering. Taxonomic assignments of all taxa with a relative abundance ≥1.0% for each ecostate are shown in Supplementary Table 2.
Approximately 50% of the taxa assigned to each of the different ecostates, either at the order, family, or genus level, could be associated to known physiological representatives that are presumed to be either predominately anaerobic or aerobic. For example, all known relatives of the family Desulfobulbaceae, or genus Desulfococcus, predominant in Ecostates 1 and 2 are known to be anaerobes; while known relatives of the family Flavobacteriaceae, and family Marcinicellaceae (both abundant in ecostate 3) are obligate aerobes. From this we surmise that the overall dominant ecostates 1 and 2 had relatively high abundances of putative anaerobes, compared to putative aerobes (Supplementary Table 2), while ecostate 3, found primarily in the top 1 cm (Figure 5) had higher relative abundances of putative aerobes, and corroborates the availability of oxygen declining with depth (Glud, 2008). In general, we found members of the phyla Bacteroidetes more prevalent in ecostates 2, 3, and 6, all associated the upper 5 cm of the sediments than in deeper ecostates 1 and 7. The Bacteroidetes are known for their involvement in the degradation of more complex organic matter, either aerobically or anaerobically (Kirchman, 2002; Hahnke et al., 2016). This is consistent with the most labile organic matter being present in the upper layers of sediment, in part driven by spring and fall phytoplankton blooms. Ecostates 4, 5, and 6 were only present in April, and OTUs that increased in abundance were Flavobacteriaceae (Otu00006), Rhodobacteraceae (Otu00022), Vibrionales (Otu00051), Helicobacteraceae (Otu00052), unclassified Bacteroidetes (Otu00009), Clostridiales (Otu00066), and Alteromonadales (Otu00012) OTUs (Supplementary Table 1). These OTUs may be involved in the degradation of phytoplankton-derived organic matter—relatives of some of these OTUs have been observed to change abundance in response to phytoplankton-derived organic matter as well as plant-derived organic matter (Raulf et al., 2015). OTUs that were elevated in April were similar to the elevated OTUs in ecostate 3, which further supports the hypothesis that these OTUs are potentially involved in organic matter-specific responses such as carbon mineralization. Although these microbial responses to phytoplankton blooms are tentative, it would be worth exploring in other coastal environments with respect to seasonal variation.
Seasonal Dynamics of Eukaryotic Phytoplankton in Sediments
An important ancillary question related to our focus on iron-cycling communities, is the role for pelagic-derived phytoplankton, such as diatoms, to be transported to these intertidal sediments by sedimentation and/or bio-mixing activity of suspension feeding macrofauna. These phytoplankton may act as an important source of labile organic matter in sediments at different times of the year (Kristensen, 2001; Kristensen and Mikkelsen, 2003). To address this, we used 18S rRNA gene sequencing to determine the presence and relative abundance of phytoplankton and other eukaryotes in these sediments over the seasonal cycle. We recovered 18S rRNA gene amplicons from 99/100 samples over the different seasons from depths of 1–10 cm. However, the total number of quality filtered sequences from each sample varied significantly with a median of 2,636 sequences per sample and a range of 1–24,484 sequences. One explanation for this variability is that it reflects biological reality, in that there are simply few eukaryotes present in those low-yield samples, thus yielding few 18S rRNA gene sequences. Pasulka et al. (2019) utilized a similar approach to investigate benthic microbial eukaryotes at Guaymas Basin hydrothermal vents, and also recovered a highly variable numbers of sequences among samples. Further sequencing of samples from marine benthic habitats will be needed to understand the relationship between sequence yields, the biomass of microbial eukaryotes or meio-fauna in unprocessed samples, DNA template concentrations in extracted samples, and possible PCR inhibition in these complex and heterogeneous environmental DNA samples. It cannot be ruled out that some sequences are sourced from relict, extracellular DNA bound to sediment particles, and are not indicative of the conditions that prevailed during the course of this study. However, a recent paper looking at marine sediments found little influence of extracellular DNA on amplicon studies similar to these samples (Ramírez et al., 2018).
With these caveats in mind, we found a total of 119 eukaryotic OTUs present at a relative abundance of ≥0.1% of total sequences (range = 0.10–3.87%; Supplementary Table 4). Of these OTUs, 34 were identified as belonging to the two most abundant classes of eukaryotic marine phytoplankton, including the diatoms (Bacillariophyta = 18 OTUs) and dinoflagellates (n = 16 OTUs). For illustrative purposes, the spatiotemporal distribution of diatoms are shown in Figure 7. While there may be temporal patterns to the diatom distribution, further sampling and correlation to additional parameters like quantification of chlorophyll and cell counts in the overlaying water will be necessary to make these associations. The dinoflagellates (data in Supplementary Table 4) showed similar patterns. The fact that significant diatom abundance was found in the depth range of 2–5 cm, i.e., below the water-sediment interface, Figure 7, indicates the transport of diatoms, or diatom-DNA, into the subsurface sediments presumably via bio-mixing by marine worms. This would make an important source of labile carbon that would help sustain biogeochemical cycling of iron and sulfur, and could help explain the profiles of Fe(II) and FeHR discussed above.
Figure 7. Spatiotemporal relative abundance of diatoms (Bacillariophyta) in sediments identified by 18S rRNA gene sequencing.
Many of the most abundant phytoplankton OTUs were associated with benthic habitats (Supplementary Table 4). Diatoms and dinoflagellates were identified as four of the ten most abundant OTUs and included a dinoflagellate from the order Peridiniales (OTU0004, 3.87%), a raphid-pennate diatom (OTU0006, 3.24%), the diatom Psammodictyon sp. (OTU0009, 2.85%) and the dinoflagellate Islandinium minutum (OTU0010, 2.72%). The Peridiniales include a diverse group of dinoflagellates rendering the taxonomic designation of OTU0004 somewhat ambiguous. Nonetheless, OTU0004 is notable with respect to its prevalence across depth, with greater proportions found deeper in the sediments (Supplementary Table 4). Raphid-pennate diatoms include both pelagic and benthic species, but the motile benthic species are thought to be substantially more diverse than pelagic types (Mann et al., 2017), and are likely the type associated with OTU0006. The diatom Psammodictyon (OTU0009) is also associated with benthic habitats, living at the sediment water interface (Round et al., 1990). Islandinium minutum (OTU0010) is primarily associated with a particular “spiny” morphology of dinoflagellate cysts that are common in high latitude sediments (Potvin et al., 2013). The extent to which phytoplankton distributions in the sediments can be used as indicators of subsurface biomixing and ventilation remains to be tested under more controlled experimental conditions, but the possibility exists that these organisms could serve as useful biomarkers for this important biogeochemical process.
Sedimentary Macro-, and Meio-Fauna
The 18S rRNA sequences also revealed the presence of Nereis (Hediste) diversicolor, which was the presumed dominant polychaete worm in this mudflat. This was the most abundant polychaete 18S rRNA gene sequence, accounting for 3.7% of the total 18S rRNA gene sequences (n = 14,149 sequences). Mollusca 18S rRNA gene sequences were not a very abundant component of any samples, although we detected 0.04% of the total 18S rRNA gene sequences (n = 148 sequences) related to the small, intertidal-dwelling clam, Macoma balthica, which is common to the Sheepscot Estuary and to North Atlantic intertidal sediments (Stickney, 1959). We have observed small clams similar to M. balthica from the Eddy on numerous occasions.
Nematodes can comprise a large fraction of the biomass of marine sediments (Heip et al., 1985) and are bewilderingly diverse (Bik et al., 2010). The 18S rRNA gene libraries revealed that nematodes comprised a median relative abundance of 21.9% with a range of 0–67.2% of the total 18S rRNA gene sequences. Their relative abundance was highest near the sediment-water interface, and abundance declined with depth (Figure 8). This finding corroborates observations primarily based from direct species counts and identification (Montagna et al., 1989). Nematode abundance and grazing activity have also been observed to increase bacterial activity and diversity in terrestrial sediments (Jiang et al., 2017). Nematodes play an important role in grazing of bacteria and other microbial eukaryotes in sediments (Pascal et al., 2008), and potentially play an important role in controlling iron biogeochemical dynamics.
Figure 8. Spatiotemporal relative abundance of nematodes (Nematoda) identified by 18S rRNA gene sequencing over the seasonal cycle.
We observed spatial and temporal heterogeneity in both porewater dissolved Fe(II) and highly reactive Fe (FeHR) sedimentary pools over the seasonal cycle that are primarily controlled by the master variable, temperature. We infer that changes in temperature influence both animal feeding and irrigational behavior (Esselink and Zwarts, 1989; Kristensen, 2001; Fang et al., 2019) as well as the direct influence of temperature on increased microbial iron oxidation and reduction. The sedimentary microbial communities were relatively stable throughout the seasonal cycle with a major demarcation between surface (0–5 cm) and subsurface sediments (5–10 cm). Observed perturbations to the microbial community as assessed by ecostates tended to revert to the steady-state ecostate conditions. We speculate that sediment mixing and irrigational behavior of bioturbating animals may be a contributing factor to this microbial community resilience. The presence and quantification of abundances of phytoplankton, especially at depth, in these bioturbated sediments may provide clues about changes in macrofaunal feeding behavior and potential impact on the iron biogeochemical cycle across the seasons. Our observations also suggest that specific benthic micro- and meiofauna associations or interactions likely play an, as yet poorly understood, role in controlling the sedimentary iron biogeochemical cycle. In the rapidly changing coastal ocean system, it is imperative to conduct both temporal and spatial monitoring of sedimentary ecosystems, as environmental perturbations could be reflected in the microbial communities and the essential biogeochemical transformations they catalyze.
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.
JB and DE conceived study with input from PG and DJ. SG assisted with field and laboratory analysis. PC helped analyze eukaryotic OTUs. NR helped with ecostates analysis. JB, DE, and PC wrote manuscript with input from all authors.
Primary funding for this research was provided by a National Science Foundation grant OCE-1459600 to DE, PG, and DJ. Additional support came from NSF OIA-1849227 to DE, NR, and PC, and from NASA grant NNX16AG59G, and Bigelow Laboratory institutional funds.
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.
JB thanks Dr. Peter Larsen for many helpful and insightful discussions about benthic macrofaunal behavior and identification, and Dr. Jarrod Scott for inspiring the microbial community sampling design and insightful discussions on microbial community analysis across spatial and time series data. We also thank Dr. Alex Michaud for insightful comments on iron and sulfur sedimentary biogeochemistry and Dr. Jeremy Rich for helpful comments on the manuscript. We appreciate the assistance of organizing and shipping DNA samples for sequencing by Jennifer Delany at Harvard University.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2020.562617/full#supplementary-material
Supplementary Table 1 | Seasonal data for each sampling.
Supplementary Table 2 | The top 76 OTUs for each ecostate in order from the highest relative abundance to the least relative abundance.
Supplementary Table 3 | This table is derived from the different ecostate taxa counts from Supplementary Table 2. For each ecostate, all taxa of ≥1% relative abundance are represented. Taxa are classified to the highest taxonomic resolution against the Silva database, along with either their phylum, or in the case of the Proteobacteria, their corresponding class, abbreviated as Delta, Gamma, Alpha, or Epsilon. The colors denote whether taxa were judged to be presumptive anaerobes (red), possible anaerobes (purple), or presumptive aerobes (green). These physiological categories were assigned based on multiple cultured representatives of family, order, or class being known anaerobes or aerobes. For approximately half the taxa there were not enough cultured relatives known to make any judgment. For each ecostate, the % of Gammaproteobacteria, Deltaproteobacteria, and Bacteroidetes are given, since these were the three most abundant taxa within all the ecostates. The % anaerobes and % aerobes are the sums of the relative abundances of presumptive and possible anaerobes, or presumptive aerobes for each ecostate.
Supplementary Table 4 | Eukaryotic operational taxonomic units (OTUs) and taxonomic classifications from all seasonal samples.
Aller, R. C. (1982). “The effects of macrobenthos on chemical properties of marine sediment and overlying water,” in Animal-Sediment Relations. Topics in Geobiology, Vol. 100, eds P. L. McCall and M. J. S. Tevesz (Boston, MA: Springer). doi: 10.1007/978-1-4757-1317-6_2
Aller, R. C., Aller, J. Y., Zhu, Q., Heilbrun, C., Klingensmith, I., and Kaushik, A. (2019). Worm tubes as conduits for the electrogenic microbial grid in marine sediments. Sci. Adv. 5:eaaw3651. doi: 10.1126/sciadv.aaw3651
Antler, G., Mills, J. V., Hutchings, A. M., Redeker, K. R., and Turchyn, A. V. (2019). The sedimentary carbon-sulfur-iron interplay – a lesson from East Anglian salt marsh sediments. Front. Earth Sci. 7:140. doi: 10.3389/feart.2019.00140
Beam, J. P., Scott, J. J., McAllister, S. M., Chan, C. S., McManus, J., Meysman, F. J. R., et al. (2018). Biological rejuvenation of iron oxides in bioturbated marine sediments. ISME J. 12, 1389–1394. doi: 10.1038/s41396-017-0032-6
Bertics, V. J., and Ziebis, W. (2009). Biodiversity of benthic microbial communities in bioturbated coastal sediments is controlled by geochemical microniches. ISME J. 3, 1269–1285. doi: 10.1038/ismej.2009.62
Bertics, V. J., and Ziebis, W. (2010). Bioturbation and the role of microniches for sulfate reduction in coastal marine sediments. Environ. Microbiol. 12, 3022–3034. doi: 10.1111/j.1462-2920.2010.02279.x
Bik, H. M., Lambshead, P. J. D., Thomas, W. K., and Lunt, D. H. (2010). Moving towards a complete molecular framework of the Nematoda: a focus on the Enoplida and early-branching clades. BMC Evol. Biol. 10:353. doi: 10.1186/1471-2148-10-353
Bouma, T. J., Olenin, S., Reise, K., and Ysebaert, T. (2009). Ecosystem engineering and biodiversity in coastal sediments: posing hypotheses. Helgol. Mar. Res. 63, 95–106. doi: 10.1007/s10152-009-0146-y
Buongiorno, J., Herbert, L. C., Wehrmann, L. M., Michaud, A. B., Laufer, K., Røy, H., et al. (2019). Complex microbial communities drive iron and sulfur cycling in arctic fjord sediments. Appl. Environ. Microbiol. 85:e00949-19. doi: 10.1128/AEM.00949-19
Callahan, B. J., McMurdie, P. J., and Holmes, S. P. (2017). Exact sequence variants should replace operational taxonomic units in marker-gene data analysis. ISME J. 11, 2639–2643. doi: 10.1038/ismej.2017.119
Canfield, D. E., Thamdrup, B., and Hansen, J. W. (1993). The anaerobic degradation of organic matter in Danish coastal sediments: iron reduction, manganese reduction, and sulfate reduction. Geochim. Cosmochim. Acta 57, 3867–3883. doi: 10.1016/0016-7037(93)90340-3
Caron, D. A., Countway, P. D., Savai, P., Gast, R. J., Schnetzer, A., Moorthi, S. D., et al. (2009). Defining DNA-based operational taxonomic units for microbial-eukaryote ecology. Appl. Environ. Microbiol. 75, 5797–5808. doi: 10.1128/AEM.00298-09
Chen, X., Andersen, T. J., Morono, Y., Inagaki, F., Jørgensen, B. B., and Lever, M. A. (2017). Bioturbation as a key driver behind the dominance of Bacteria over Archaea in near-surface sediment. Sci. Rep. 7:2400. doi: 10.1038/s41598-017-02295-x
Christensen, B., Vedel, A., and Kristensen, E. (2000). Carbon and nitrogen fluxes in sediment inhabited by suspension-feeding (Nereis diversicolor) and non-suspension-feeding (N. virens) polychaetes. Mar. Ecol. Prog. Ser. 192, 203–217. doi: 10.3354/meps192203
Cleary, D. F. R., Coelho, F. J. R. C., Oliveira, V., Gomes, N. C. M., and Polónia, A. R. M. (2017). Sediment depth and habitat as predictors of the diversity and composition of sediment bacterial communities in an inter-tidal estuarine environment. Mar. Ecol. 38:e12411. doi: 10.1111/maec.12411
Dyksma, S., Bischof, K., Fuchs, B. M., Hoffmann, K., Meier, D., Meyerdierks, A., et al. (2016). Ubiquitous Gammaproteobacteria dominate dark carbon fixation in coastal sediments. ISME J. 10, 1939–1953. doi: 10.1038/ismej.2015.257
Elrod, V. A., Berelson, W. M., Coale, K. H., and Johnson, K. S. (2004). The flux of iron from continental shelf sediments: a missing source for global budgets. Geophys. Res. Lett. 31, 2–5. doi: 10.1029/2004GL020216
Fang, X., Mestdagh, S., Ysebaert, T., Moens, T., Soetaert, K., and Van Colen, C. (2019). Spatio-temporal variation in sediment ecosystem processes and roles of key biota in the Scheldt estuary. Estuar. Coast. Shelf Sci. 222, 21–31. doi: 10.1016/j.ecss.2019.04.001
Fenchel, T. (1969). The ecology of marine microbenthos IV. Structure and function of the benthic ecosystem, its chemical and physical factors and the microfauna commuities with special reference to the ciliated protozoa. Ophelia 6, 1–182. doi: 10.1080/00785326.1969.10409647
Froelich, P. N., Klinkhammer, G. P., Bender, M. L., Luedtke, N. A., Heath, G. R., Cullen, D., et al. (1979). Early oxidation of organic matter in pelagic sediments of the eastern equatorial Atlantic: suboxic diagenesis. Geochim. Cosmochim. Acta 43, 1075–1090. doi: 10.1016/0016-7037(79)90095-4
Griffiths, J. R., Kadin, M., Nascimento, F. J. A., Tamelander, T., Törnroos, A., Bonaglia, S., et al. (2017). The importance of benthic–pelagic coupling for marine ecosystem functioning in a changing world. Glob. Chang. Biol. 23, 2179–2196. doi: 10.1111/gcb.13642
Guillou, L., Bachar, D., Audic, S., Bass, D., Berney, C., Bittner, L., et al. (2013). The Protist Ribosomal Reference database (PR2): a catalog of unicellular eukaryote Small Sub-Unit rRNA sequences with curated taxonomy. Nucleic Acids Res. 41, 597–604. doi: 10.1093/nar/gks1160
Hahnke, R. L., Meier-Kolthoff, J. P., García-López, M., Mukherjee, S., Huntemann, M., and Ivanova, N. N. (2016). Genome-based taxonomic classification of bacteroidetes. Front. Microbiol. 7:2003. doi: 10.3389/fmicb.2016.02003
Hines, M. E., Bazylinski, D. A., Tugel, J. B., and Berry Lyons, W. (1991). Anaerobic microbial biogeochemistry in sediments from two Basins in the Gulf of Maine: evidence for iron and manganese reduction. Estuar. Coast. Shelf Sci. 32, 313–324. doi: 10.1016/0272-7714(91)90046-E
Hines, M. E., Berry Lyons, W., Armstrong, P. B., Orem, W. H., Spencer, M. J., Gaudette, H. E., et al. (1984). Seasonal metal remobilization in the sediments of Great Bay, New Hampshire. Mar. Chem. 15, 173–187. doi: 10.1016/0304-4203(84)90014-8
Hu, S. K., Liu, Z., Lie, A. A. Y., Countway, P. D., Kim, D. Y., Jones, A. C., et al. (2015). Estimating protistan diversity using high-throughput sequencing. J. Eukaryot. Microbiol. 62, 688–693. doi: 10.1111/jeu.12217
Jiang, Y., Liu, M., Zhang, J., Chen, Y., Chen, X., Chen, L., et al. (2017). Nematode grazing promotes bacterial community dynamics in soil at the aggregate level. ISME J. 11, 2705–2717. doi: 10.1038/ismej.2017.120
Kallmeyer, J., Pockalny, R., Adhikari, R. R., Smith, D. C., and D’Hondt, S. (2012). From the Cover: global distribution of microbial abundance and biomass in subseafloor sediment. Proc. Natl. Acad. Sci. U.S.A. 109, 16213–16216. doi: 10.1073/pnas.1203849109
Kim, B.-S., Kim, B. K., Lee, J.-K., Kim, M., Lim, Y. W., and Chun, J. (2008). Rapid phylogenetic dissection of prokaryotic community structure in tidal flat using pyrosequencing. J. Microbiology. 46, 357–363. doi: 10.1007/s12275-008-0071-9
Kozich, J. J., Westcott, S. L., Baxter, N. T., Highlander, S. K., and Schloss, P. D. (2013). Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the miseq illumina sequencing platform. Appl. Environ. Microbiol. 79, 5112–5120. doi: 10.1128/AEM.01043-13
Kristensen, E., and Kostka, J. E. (2005). “Macrofaunal burrows and irrigation in marine sediment: microbiological and biogeochemical interactions,” in Interactions between Macro- and Microorganisms in Marine Sediments, eds E. Kristensen, J. E. Kostka, and R. Haese (Washington, DC: American Geophysical Union), 125–157. doi: 10.1029/CE060p0125
Kristensen, E., and Mikkelsen, O. L. (2003). Impact of the burrow-dwelling polychaete Nereis diversicolor on the degradation of fresh and aged macroalgal detritus in a coastal marine sediment. Mar. Ecol. Prog. Ser. 265, 141–153. doi: 10.3354/meps265141
Lasher, C., Dyszynski, G., Everett, K., Edmonds, J., Ye, W., Sheldon, W., et al. (2009). The diverse bacterial community in intertidal, anaerobic sediments at Sapelo Island, Georgia. Microb. Ecol. 58, 244–261. doi: 10.1007/s00248-008-9481-9
Laufer, K., Nordhoff, M., Halama, M., Martinez, R., Obst, M., Nowak, M., et al. (2017). Microaerophilic Fe(II)-oxidizing Zetaproteobacteria isolated from low-Fe marine coastal sediments – physiology and characterization of their twisted stalks. Appl. Environ. Microbiol. 83:e03118-16. doi: 10.1128/AEM.03118-16
Laufer, K., Nordhoff, M., Schmidt, C., Behrens, S., Jørgensen, B. B., and Kappler, A. (2016). Co-existence of microaerophilic, nitrate-reducing, and phototrophic Fe(II)-oxidizers and Fe(III)-reducers in coastal marine sediment. Appl. Environ. Microbiol. 82, 1433–1447. doi: 10.1128/AEM.03527-15
McAllister, S. M., Barnett, J. M., Heiss, J. W., Findlay, A. J., MacDonald, D. J., Dow, C. L., et al. (2015). Dynamic hydrologic and biogeochemical processes drive microbially enhanced iron and sulfur cycling within the intertidal mixing zone of a beach aquifer. Limnol. Oceanogr. 60, 329–345. doi: 10.1111/lno.10029
McBeth, J. M., Little, B. J., Ray, R. I., Farrar, K. M., and Emerson, D. (2011). Neutrophilic iron-oxidizing “Zetaproteobacteria” and mild steel corrosion in nearshore marine environments. Appl. Environ. Microbiol. 77, 1405–1412. doi: 10.1128/AEM.02095-10
Melton, E. D., Swanner, E. D., Behrens, S., Schmidt, C., and Kappler, A. (2014). The interplay of microbially mediated and abiotic reactions in the biogeochemical Fe cycle. Nat. Rev. Microbiol. 12, 797–809. doi: 10.1038/nrmicro3347
Meysman, F. J. R., Risgaard-Petersen, N., Malkin, S. Y., and Nielsen, L. P. (2015). The geochemical fingerprint of microbial long-distance electron transport in the seafloor. Geochim. Cosmochim. Acta 152, 122–142. doi: 10.1016/j.gca.2014.12.014
Montagna, P. A., Bauer, J. E., Hardin, D. H., and Spies, R. B. (1989). Vertical distribution of microbial and meiofaunal populations in sediments of a natural coastal hydrocarbon seep. J. Mar. Res. 47, 657–680. doi: 10.1357/002224089785076226
Nielsen, L. P., and Risgaard-Petersen, N. (2015). Rethinking sediment biogeochemistry after the discovery of electric currents. Annu. Rev. Mar. Sci. 7, 425–442. doi: 10.1146/annurev-marine-010814-015708
Nielsen, L. P., Risgaard-Petersen, N., Fossing, H., Christensen, P. B., and Sayama, M. (2010). Electric currents couple spatially separated biogeochemical processes in marine sediment. Nature 463, 1071–1074. doi: 10.1038/nature08790
Otte, J. M., Harter, J., Laufer, K., Blackwell, N., Straub, D., Kappler, A., et al. (2018). The distribution of active iron-cycling bacteria in marine and freshwater sediments is decoupled from geochemical gradients. Environ. Microbiol. 20, 2483–2499. doi: 10.1111/1462-2920.14260
Pascal, P. Y., Dupuy, C., Richard, P., Rzeznik-Orignac, J., and Niquil, N. (2008). Bacterivory of a mudflat nematode community under different environmental conditions. Mar. Biol. 154, 671–682. doi: 10.1007/s00227-008-0960-9
Pasulka, A., Hu, S. K., Countway, P. D., Coyne, K. J., Cary, S. C., Heidelberg, K. B., et al. (2019). SSU-rRNA gene sequencing survey of benthic microbial eukaryotes from guaymas basin hydrothermal vent. J. Eukaryot. Microbiol. 66, 637–653. doi: 10.1111/jeu.12711
Potvin, É, Rochon, A., and Lovejoy, C. (2013). Cyst-theca relationship of the arctic dinoflagellate cyst Islandinium minutum (Dinophyceae) and phylogenetic position based on SSU rDNA and LSU rDNA. J. Phycol. 49, 848–866. doi: 10.1111/jpy.12089
Poulton, S. W., and Canfield, D. E. (2005). Development of a sequential extraction procedure for iron: implications for iron partitioning in continentally derived particulates. Chem. Geol. 214, 209–221. doi: 10.1016/j.chemgeo.2004.09.003
Qiao, Y., Liu, J., Zhao, M., and Zhang, X. H. (2018). Sediment depth-dependent spatial variations of bacterial communities in mud deposits of the Eastern China marginal seas. Front. Microbiol. 9:1128. doi: 10.3389/fmicb.2018.01128
Ramírez, G. A., Jørgensen, S. L., Zhao, R., and D’Hondt, S. (2018). Minimal influence of extracellular DNA on molecular surveys of marine sedimentary communities. Front. Microbiol. 9:2969. doi: 10.3389/fmicb.2018.02969
Raulf, F. F., Fabricius, K., Uthicke, S., de Beer, D., Abed, R. M. M., and Ramette, A. (2015). Changes in microbial communities in coastal sediments along natural CO2 gradients at a volcanic vent in Papua New Guinea. Environ. Microbiol. 17, 3678–3691. doi: 10.1111/1462-2920.12729
Record, N. R., O’Brien, J. D., Stamieszkin, K., and Runge, J. A. (2017). Omic-Style statistical clustering reveals old and new patterns in the Gulf of Maine ecosystem. Can. J. Fish. Aquat. Sci. 74, 1061–1076. doi: 10.1139/cjfas-2016-0151
Reyes, C., Schneider, D., Thürmer, A., Kulkarni, A., Lipka, M., Sztejrenszus, S. Y., et al. (2017). Potentially active iron, sulfur, and sulfate reducing bacteria in Skagerrak and Bothnian bay sediments. Geomicrobiol. J. 34, 840–850. doi: 10.1080/01490451.2017.1281360
Rubin-Blum, M., Antler, G., Tsadok, R., Shemesh, E., Austin, J. A., Coleman, D. F., et al. (2014). First evidence for the presence of iron oxidizing zetaproteobacteria at the levantine continental margins. PLoS One 9:e0091456. doi: 10.1371/journal.pone.0091456
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
Scott, J. J., Glazer, B. T., and Emerson, D. (2017). Bringing microbial diversity into focus: high-resolution analysis of iron mats from the Lō‘ihi Seamount. Environ. Microbiol. 19, 301–316. doi: 10.1111/1462-2920.13607
Seeberg-Elverfeldt, J., Schluter, M., Feseker, T., and Kolling, M. (2005). Rhizon sampling of porewaters near the sediment-water interface of aquatic systems. Limnol. Oceanogr. 3, 361–371. doi: 10.1016/S0012-821x(02)01064-6
Severmann, S., McManus, J., Berelson, W. M., and Hammond, D. E. (2010). The continental shelf benthic iron flux and its isotope composition. Geochim. Cosmochim. Acta 74, 3984–4004. doi: 10.1016/j.gca.2010.04.022
Staudinger, M. D., Mills, K. E., Stamieszkin, K., Record, N. R., Hudak, C. A., Allyn, A., et al. (2019). It’s about time: a synthesis of changing phenology in the Gulf of Maine ecosystem. Fish. Oceanogr. 28, 532–566. doi: 10.1111/fog.12429
Sulu-Gambari, F., Seitaj, D., Behrends, T., Banerjee, D., Meysman, F. J. R., and Slomp, C. P. (2016). Impact of cable bacteria on sedimentary iron and manganese dynamics in a seasonally-hypoxic marine basin. Geochim. Cosmochim. Acta 192, 49–69. doi: 10.1016/j.gca.2016.07.028
Thamdrup, B., Fossing, H., and Jorgensen, B. B. (1994). Manganese, iron, and sulfur cycling in a coastal marine sediment, Aarhus Bay, Denmark. Geochim. Cosmochim. Acta 58, 5115–5129. doi: 10.1016/0016-7037(94)90298-4
Thibault de Chanvalon, A., Metzger, E., Mouret, A., Knoery, J., Geslin, E., and Meysman, F. J. R. (2016). Two dimensional mapping of iron release in marine sediments at submillimetre scale. Mar. Chem. 191, 34–49. doi: 10.1016/j.marchem.2016.04.003
van de Velde, S. J., Hidalgo-Martinez, S., Callebaut, I., Antler, G., James, R. K., Leermakers, M., et al. (2020). Burrowing fauna mediate alternative stable states in the redox cycling of salt marsh sediments. Geochim. Cosmochim. Acta 276, 31–49. doi: 10.1016/j.gca.2020.02.021
Vandieken, V., Finke, N., and Jørgensen, B. B. (2006). Pathways of carbon oxidation in an Arctic fjord sediment (Svalbard) and isolation of psychrophilic and psychrotolerant Fe(III)-reducing bacteria. Mar. Ecol. Prog. Ser. 322, 29–41. doi: 10.3354/meps322029
Volkenborn, N., Polerecky, L., Wethey, D. S., and Woodin, S. A. (2010). Oscillatory porewater bioadvection in marine sediments induced by hydraulic activities of Arenicola marina. Limnol. Oceanogr. 55, 1231–1247. doi: 10.4319/lo.2010.55.3.1231
Wang, Y., Sheng, H. F., He, Y., Wu, J. Y., Jiang, Y. X., Tam, N. F. Y., et al. (2012). Comparison of the levels of bacterial diversity in freshwater, intertidal wetland, and marine sediments by using millions of illumina tags. Appl. Environ. Microbiol. 78, 8264–8271. doi: 10.1128/AEM.01821-12
Yarza, P., Yilmaz, P., Pruesse, E., Glöckner, F. O., Ludwig, W., Schleifer, K. H., et al. (2014). Uniting the classification of cultured and uncultured bacteria and archaea using 16S rRNA gene sequences. Nat. Rev. Microbiol. 12, 635–645. doi: 10.1038/nrmicro3330
Keywords: ferrous iron, marine sediment, iron cycle, bioturbation, Zetaproteobacteria
Citation: Beam JP, George S, Record NR, Countway PD, Johnston DT, Girguis PR and Emerson D (2020) Mud, Microbes, and Macrofauna: Seasonal Dynamics of the Iron Biogeochemical Cycle in an Intertidal Mudflat. Front. Mar. Sci. 7:562617. doi: 10.3389/fmars.2020.562617
Received: 15 May 2020; Accepted: 09 October 2020;
Published: 16 November 2020.
Edited by:Gordon T. Taylor, Stony Brook University, United States
Reviewed by:Eric Boschker, Delft University of Technology, Netherlands
Sebastiaan van de Velde, Université libre de Bruxelles, Belgium
Copyright © 2020 Beam, George, Record, Countway, Johnston, Girguis and Emerson. 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: David Emerson, email@example.com