- 1Marine and Environmental Biology, University of Southern California, Los Angeles, CA, United States
- 2Earth System Science, Stanford University, Stanford, CA, United States
Macroalgae, commonly known as seaweed, are foundational species in coastal ecosystems and contribute significantly to coastal primary production globally. However, the impact of macroalgal decomposition on benthic biological nitrogen fixation (BNF) after deposition to the seafloor remains largely unexplored. In this study, we measure BNF rates at three different sites at the Big Fisherman's Cove on Santa Catalina Island, CA, USA, which is representative of globally distributed rocky bottom macroalgal habitats. Unamended BNF rates varied among sites (0.001–0.05 nmol N g−1 h −1) and were generally within the lower end of previously reported ranges. We hypothesized that the differences in BNF between sites were linked to the availability of organic matter. Indeed, additions of glucose, a labile carbon source, resulted in 2–3 orders of magnitude stimulation of BNF rates in bottle incubations of sediment from all sites. To assess the impact of complex, autochthonous organic matter, we simulated macroalgal deposition and remineralization with additions of brown (i.e., Macrocystis pyrifera and Dictyopteris), green (i.e., Codium fragile), and red (i.e., Asparagopsis taxiformis) macroalgae. While brown and green macroalgal amendments resulted in 53- to 520-fold stimulation of BNF rates—comparable to the labile carbon addition—red alga was found to significantly inhibit BNF rates. Finally, we employed nifH sequencing to characterize the diazotrophic community associated with macroalgal decomposition. We observed a distinct community shift in potential diazotrophs from primarily Gammaproteobacteria in the early stages of remineralization to a community dominated by Deltaproteobacteria (e.g., sulfate reducers), Bacteroidia, and Spirochaeta toward the latter phase of decomposition of brown, green, and red macroalgae. Notably, the nifH-containing community associated with red macroalgal detritus was distinct from that of brown and green macroalgae. Our study suggests coastal benthic diazotrophs are limited by organic carbon and demonstrates a significant and phylum-specific effect of macroalgal loading on benthic microbial communities.
Introduction
Nitrogen (N) is an important macronutrient that governs microbial community structure and is often considered to be limiting in many regions of the ocean (Zehr and Capone, 2020). Diazotrophs or N2 fixers are specialized prokaryotes that utilize the nitrogenase enzyme to carry out biological N fixation (BNF), the process of converting dinitrogen (N2) gas into a bioavailable source of N. These microorganisms have been well-investigated in open-ocean oligotrophic regions (Zehr and Capone, 2023). In comparison, benthic coastal ecosystems, which are considered to be eutrophic and areas of net N2 loss, are often under-surveyed for BNF (Fulweiler, 2023). However, numerous studies report BNF by an assortment of diazotrophic phylotypes in illuminated coastal habitats with different benthic primary producers such as photosynthetic microbial mats (Zehr et al., 1995; Paerl et al., 1996), seagrass (Capone, 1982; Welsh, 2000; Mohr et al., 2021), and salt marsh plants (Nedwell and Azni bin Abdul Aziz, 1980; Gandy and Yoch, 1988). More recently, there has been a growing body of evidence suggesting that benthic BNF is ubiquitous throughout global coastlines (Andersson et al., 2014; Hou et al., 2018; Jabir et al., 2018, 2021).
Notably, benthic macroalgae contribute significantly to global coastal primary production and create important coastal habitats (Duarte and Cebrián, 1996; Krause-Jensen and Duarte, 2016). There have also been many reports of diazotrophic macroalgal associations (DMAs) with living macroalgae (Head and Carpenter, 1975; Capone et al., 1977; Johnson et al., 2023). In recent years, there has been a major uptick in macroalgal blooms, often linked to anthropogenic nutrient loading, affecting coastal marine ecosystems, often with detrimental ecological impacts throughout the world (Lyons et al., 2014; Lapointe et al., 2018; Wang et al., 2020). Despite their global occurrence and potential to harbor diazotrophs, coastal macroalgal habitats and the underlying marine sediments remain largely under-investigated for BNF. Therefore, one major objective of this study is to investigate the prevalence of benthic BNF along a coastal rocky bottom habitat, i.e., Big Fisherman's Cove (BFC), which often house dense kelp forests, diverse macroalgae, and other macrophytes such as eelgrass (e.g., Zostera) beds (Tanner et al., 2019; Looby and Ginsburg, 2021; Ginsburg and Huang, 2022).
In many parts of the oligotrophic ocean (e.g., Pacific Ocean and S. Atlantic), iron, an important constituent of the nitrogenase metalloenzyme, often constrains diazotrophic biogeography (Sohm et al., 2011). In other regions of the ocean (e.g., N. Atlantic), phosphorous has been reported to limit diazotrophic activity (Sohm et al., 2008). Numerous studies have explored how similar nutrients (e.g., iron) impact diazotrophic communities in benthic ecosystems (Howarth et al., 1988; Yao et al., 2023). While carbon is rarely thought to limit photosynthetic diazotrophs inhabiting the euphotic zone, the availability of organic matter, ranging from simple carbon substrates to whale falls, has routinely been reported to enhance heterotrophic microbial activity, including BNF, in diverse marine sediments (Meyer-Reil, 1987; Bertics et al., 2010; Gier et al., 2016; Polyak et al., 2017; Dekas et al., 2018). In this study, we similarly expect the addition of labile organic substrates (i.e., glucose and macroalgal detritus) to stimulate benthic BNF rates in shallow coastal marine sediments.
Macroalgae tend to be carbon-rich, often with much higher carbon to N ratios (>15:1) than typical Redfield ratios observed with other phytodetritus such as phytoplankton (Duarte, 1992; Duarte and Cebrián, 1996). During macroalgal decomposition, the N pool can be depleted by the microbial community much faster than carbon, and accumulation of non-labile N substrates can result in N limitation (Tenore et al., 1979; Rice, 1982; Rice and Hanson, 1984; Raut et al., 2018). This provides an ecological advantage for heterotrophic diazotrophs which can readily supply their own N while tapping into the residual carbon pool of the macroalgal detritus.
In previous studies, relatively low BNF rates were associated with healthy macroalgal biomass at the beginning of long-term (16–32 days) decomposition experiments (Hamersley et al., 2015; Raut and Capone, 2021). There was a substantial increase in BNF rates following a few days (3–7) of decomposition that were likely fueled by leaching of labile carbon substrates from the macroalgal detritus and sustained for an intermediary time period (Hamersley et al., 2015; Raut and Capone, 2021). As the macroalgal detritus became more recalcitrant, nitrogenase activity declined significantly toward the latter stages of decomposition (Raut et al., 2018; Raut and Capone, 2021). Interestingly, labile carbon amendments throughout macroalgal remineralization were found to stimulate nitrogenase activity only at the beginning and terminal phases of decomposition. This suggests that rapid microbial degradation of macroalgae in the intermediate phase of decomposition can alleviate carbon limitation for the associated diazotrophic community (Raut et al., 2018; Raut and Capone, 2021).
A significant portion of macroalgal biomass undergoes burial and remineralization in the surrounding sediments of the macroalgal beds and can also be exported to deep sea sediments (Duarte and Cebrián, 1996; Dierssen et al., 2009; Krause-Jensen and Duarte, 2016). Currently, the impact of macroalgal deposition and decomposition on the benthic microbial community remains largely unknown. Therefore, another major objective of this study is to investigate the effects of macroalgal sedimentation and remineralization on benthic BNF rates in shallow coastal marine sediments. We employ acetylene reduction assays (ARAs) to measure BNF rates with sediment slurry incubations at the Wrigley Institute for Environment and Sustainability (WIES), adjacent to BFC. To assess the impact of labile carbon and autochthonous organic matter on benthic BNF rates, we carry out glucose (10 mM) additions and macroalgal amendment experiments using a variety of species across all three major clades of brown, green, and red seaweeds in parallel to the unamended sediment slurry incubations.
Sulfate-reducing bacteria are often found to be the major constituents of benthic diazotrophic communities (Capone, 1982; Gandy and Yoch, 1988; Jabir et al., 2021; Liesirova et al., 2023). Indeed, parallel benthic sulfate reduction and BNF rates further support their importance in shallow marine sediments (Bertics et al., 2010, 2013). In this study, we use additions of sodium molybdate, an inhibitor of sulfate-reducing bacteria, to assess the contribution of sulfate-reducing diazotrophs to benthic BNF at BFC. Interestingly, similar inhibition assays have been used to identify sulfate-reducing bacteria to be a major constituent of the diazotrophic community throughout macroalgal decomposition (Raut et al., 2018; Raut and Capone, 2021). However, there have been molecular surveys, typically focusing on the nifH gene, that have revealed a greater diversity of diazotrophs in coastal marine sediments (Jabir et al., 2018) and associated with macroalgal remineralization (Zhang et al., 2015; Aires et al., 2019). In this study, we leverage a similar metabarcoding approach with nifH sequencing to characterize the macroalgae-associated diazotrophic diversity and identify community shifts during the decomposition of several macroalgae typically found within the study site.
Materials and methods
Site description and sediment sampling
We surveyed three coastal marine sites for nitrogenase activity using the acetylene reduction assay (ARA) in the Big Fisherman's Cove (BFC), Santa Catalina Island, CA, USA, between spring and fall of 2017 (Figure 1A). The first location is a shallow, intertidal rocky site (~3–5 feet water depth) with mixed macroalgal assemblages consisting of various red (e.g., Asparagopsis taxiformis), green (e.g., Codium fragile), and brown (e.g., Dictyopteris) seaweed but primarily dominated by A. taxiformis (Figure 1B). The second location is a sandy bottom site (~10–15 feet water depth) situated adjacent to a stable eelgrass (Zostera) bed and shallow, rocky subtidal zone with mixed macroalgal assemblage (Figures 1C, D). The third sampling location is a sandy bottom site situated under a dock at an average water depth of ~30 feet with a low abundance of canopy-forming giant kelp (Macrocystis pyrifera) and a consistent aggregation of common Southern California fish. Site 3 is also adjacent to a rocky substratum housing many different fleshy brown macroalgae (Figures 1E, F). To investigate the role of autochthonous carbon sources on BNF rates at BFC, we carried out parallel incubations (experiments 2–4) with and without macroalgal amendments with sediments from the three sites. We measured the response of nitrogenase activity with the addition of brown (i.e., M pyrifera and Dictyopteris,), green (i.e., C. fragile), and red (i.e., A. taxiformis) macroalgae. These are commonly found at BFC (collected as healthy living seaweed at the time of sediment sampling) and along rocky bottom habitats all along the California coastline (Looby and Ginsburg, 2021; Ginsburg and Huang, 2022).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1312843/fmicb-14-1312843-HTML/image_m/fmicb-14-1312843-g001.jpg)
Figure 1. (A) Aerial view of Big Fisherman's Cove, Santa Catalina Island, California, USA depicting the three sites of survey for benthic BNF rates. (B) Site 1 (~3–5 feet water depth) is a more intertidal rocky zone populated by a mixture of macroalgae but often dominated by Asparagopsis taxiformis. (C, D) Site 2 (~10–15 feet water depth) is closer to the cliff wall where you might have terrigenous input and there's a stable population of eel grass (Zostera) and macroalgal assemblages of various brown and red macroalgae. (E, F) Site 3 is situated around the dock (~30 feet water depth) and populated by some giant kelp and steady aggregation of fish. Image credits: Yubin Raut—(A), David W. Ginsburg—(B), Maurice Roper—(C), Jelani A. Williams—(D, F), and Jasmine Henderson—(E).
Sediment core samples were collected by hand on SCUBA using impact-resistant polycarbonate round tubes (3.5-inch OD, 3.25-inch ID x 1 ft, McMaster-Carr) and chemical-resistant tapered rubber plugs (McMaster-Carr). Duplicate cores were collected at each site, with an average depth of ~20 cm in May 2017. After collection, sediment cores were stored at 15°C in a constant temperature room at the adjacent WIES and subsampled for unamended BNF rate measurements representative of in situ benthic diazotrophic activity (experiment 1).
Additionally, ~25 L of sediment was also collected by hand on SCUBA at each site. Here, 5-gallon buckets were filled with material from the top 35 cm of the sediments, capped, and brought to the surface via a lift bag. To ensure the removal of both macrofauna and large sediment components (e.g., shells, rocks, and detritus), sediments were sieved through 500-μm copper mesh immediately upon collection. As a result of the sieving, sediment grains <500 μm were homogenized, disrupting any established geochemical gradients. Sieved materials were stored in 5-gallon buckets in the WIES 15°C room (similar to in situ water temperature), with aerated, UV-treated, 0.2 μm-filtered seawater being constantly pumped across each sediment collection (~20 L/day) to replenish oxygen levels within the overlying waters. The constant flow of filtered and UV-treated seawater prevents stagnation of the overlying waters and the removal of metabolic byproducts, creating a mesocosm mimicking that of the natural environment (Lam et al., 2019). The mesocosm established in May 2017 was subsequently sampled to conduct parallel unamended and amended sediment slurry incubations for BNF rate measurements in June and July 2017 (experiments 2 and 3). An additional mesocosm was freshly collected and set up in October 2017 which was subsequently sampled to conduct parallel unamended and amended sediment slurry incubations (experiment 4).
Experimental design
There were a total of four separate experiments conducted where BNF rates were measured for a variety of unamended and amended sediment slurry incubations sampled from each of the three sites. Experiment 1 consisted of unamended BNF rate measurements for sediment sampled from intact sediment cores collected in May 2017 and are representative of in situ benthic diazotrophic activity at the three sites. For experiment 2 (June 2017), we subsampled from the initially established mesocosm (May 2017) to conduct parallel unamended and amended sediment slurry incubations. Here, we employed the following treatments: (1) unamended, (2) M. pyrifera additions, and (3) M. pyrifera + 20 mM sodium molybdate additions. For experiment 3 (July 2017), we similarly subsampled from the same mesocosm (May 2017) to conduct another set of sediment slurry incubations consisting of (1) unamended treatment and (2) A. taxiformis additions. For experiment 4 (October 2017), a new mesocosm was established and subsampled to conduct a final set of sediment slurry incubations with the following treatments: (1) unamended, (2) C. fragile additions, (3) Dictyopteris additions, (4) M. pyrifera additions, and (5) 10 mM glucose additions. All macroalgae used for the amendments were collected on the day of incubation and consisted of healthy, living tissue that was then either autoclaved or freeze-dried before being added to the incubation. Similarly, the glucose and sodium molybdate solutions were freshly made on the day of the incubation.
Nitrogenase activity determined using the acetylene reduction assay (ARA)
All nitrogenase activity measurements for sediment slurry incubations were carried out in 27 mL serum vials flushed with N2 using the acetylene reduction assay (ARA) as detailed by Capone (1993) with some modifications. The volume of sediment sampled, using a 5-mL syringe, and the volume of fresh, 0.2 μM-filtered seawater added to the slurry incubations were 9 mL of sediment: 9 mL of seawater. For any of the macroalgal amendments, an additional ~0.5 g of biomass (autoclaved or freeze-dried living macroalgal tissue) was added to each vial. For the glucose and sodium molybdate amendments, a final concentration of 10 mM glucose and 20 mM sodium molybdate, respectively, was added in the 0.2-μM filtered seawater solution. In the remaining 9 mL of head space, 2 mL of acetylene gas (C2H2)—produced by reacting water with calcium carbide—was added with a gas-tight syringe resulting in a final volume of C2H2 (~7% of total volume or ~22% of the gas phase) that sufficiently saturated the nitrogenase enzyme (Flett et al., 1976). All ARAs were carried out in triplicate alongside a triplicate set of sediment slurry incubations without C2H2 additions that were monitored in parallel and showed no background ethylene (C2H4) production.
The vials were all incubated in the dark inside the same 15°C room as where the mesocosms were stored at WIES. We subsampled 100 μL of the headspace of each vial at different time points throughout the incubation and analyzed it on a Shimadzu Mini 2 gas chromatograph equipped with a flame ionization detector and coupled to a PeakSimple Chromatography Data System that was routinely calibrated at each time point with 100 ppm C2H4 standards to accurately quantify C2H4 peak heights. The initial time point consisted of a measurement at the time of C2H2 injection to establish a baseline of C2H4 production. Subsequently, we conducted time points every 6–15 h over the course of the ~5-day incubations, with increased monitoring during the intermediate period of incubations that excluded the lag and plateau in nitrogenase activity observed at the beginning and end, respectively, of incubations. It is important to note that ARA can have pronounced impacts on the sediment microbial community, including N2 fixers, sulfate reducers, and denitrifiers, in a relatively short time span (~7 h), and thus warrant careful interpretation of the results from longer incubations used in this study (Fulweiler et al., 2015). Furthermore, ARA can impact the physiology of the organisms, e.g., inhibition of sulfate-reducing bacteria (Dekas et al., 2009), which might result in under-estimation of total nitrogenase activity as it might inhibit some diazotrophic members of any community.
Macroalgal litter bag decomposition
To observe potential shifts in the diazotrophic community over time, we sampled macroalgal detritus at multiple time points throughout separate long-term (16–32 days) decomposition experiments and sequenced the nifH genes of their associated microbial communities. The nifH gene encodes the essential iron protein subunit of nitrogenase, the enzyme that catalyzes BNF. This analysis was conducted separately from the benthic ARA incubations but included three of the four macroalgal species (i.e., Dictyopteris, C. fragile, and A. taxiformis) used in the sediment slurry incubations.
Representative macroalgae from each of the three clades (brown—Dictyopteris, green—C. fragile, and red—A. taxiformis) were separately collected between 2017 and 2019 at sites surrounding Santa Catalina Island, CA, USA, for monitoring DMAs throughout long-term macroalgal decomposition experiments as reported in Raut and Capone (2021). The freshly collected macroalgal biomass was divided into multiple individual 200 μM mesh litter bags, and a different bag was subsampled throughout the long-term decomposition of the macroalgae in flow-through seawater tanks at WIES. The different macroalgae were subsampled on the following days and immediately stored in −80°C freezers for downstream molecular work: Dictyopteris (2017—16-day decomposition experiment: days 0, 8, 10, and 14), C. fragile (2019—32-day decomposition experiment: days 0, 2, 4, 7, 10, 13, 17, 22, 26, and 32), and A. taxiformis (2017—18-day decomposition experiment: days 0, 5, 12, 15, and 18). The full nature of the long-term decomposition experimental setup is detailed in Raut and Capone (2021).
DNA extraction, nifH gene amplification, and sequencing
Bulk community DNA extractions were done with macroalgal biomass subsampled on the previously specified days using the DNeasy PowerSoil Kit (Qiagen, Hilden, Germany, Cat. No. 12888-100). The following modifications were made to the quick-start protocol: (1) ~0.5 g of wet macroalgal biomass used, (2) samples vortexed for 20 min between moderate and high speeds during the initial bead beating step in the PowerBead tubes, and (3) centrifuged at 10,000 × g for 90 s during the initial supernatant collection. DNA concentrations were quantified using the protocol recommended with the Qubit 1X dsDNA HS Assay Kits (Invitrogen, USA, Cat. No. Q33230-Q33231).
The nifH gene was amplified in a two-step PCR process, with the first reaction containing: 1 μL DNA template, 12.5 μL Takara 2x Ex Taq Hot Start Premix (Takara Bio Inc., Kusatsu, Shiga, Japan), 0.5 μL 2.5 mg/mL bovine serum albumin (New England Biolabs, Ipswich, MA, USA), 0.5 μL each of forward and reverse nifH primers (10 μM; Mehta et al., 2003), and 10 μL sterile molecular-grade water. A touchdown PCR was conducted under the conditions reported in Supplementary Table 3 of Dekas et al. (2016). The second reaction was modified to include forward and reverse Illumina overhang adapters (5 μm) and was performed as follows: initial denaturation at 95°C for 120 s, followed by eight cycles of 95°C for 30 s, 55°C for 30 s, and 72°C for 45 s, and a final elongation step at 72°C for 5 min. A negative control (1 μL of molecular-grade water as template) was included for each forward and reverse Illumina barcode pair, and a positive control (various nifH sequences from a mock community) plus two randomly chosen sample replicates were run for each reverse Illumina barcode used. Barcoded PCR products were purified using the Ampure XP bead kit (Beckman Coulter) at 0.7 × bead solution, pooled in equal concentrations, purified again at 0.8 × bead solution, and sent to the UC Davis DNA Technologies Core Facility (UC Davis, CA, USA) for 2 × 250 bp paired-end sequencing on an Illumina MiSeq platform (as described in Kapili et al., 2020).
Denoising and quality filtering nifH sequences
The nifH primer sequences were trimmed from demultiplexed fastq files using cutadapt (v.3.5; Martin, 2011), reads were quality filtered and pooled, chimeras were removed, and amplicon sequence variants (ASVs) were inferred using DADA2 (V1.12.1; Callahan et al., 2016). Further quality control included pruning ASVs that did not align to the nifH target region, and the remaining ASVs were incorporated into the reference alignment and tree using SEPP (v.4.4.0; Mirarab et al., 2012) as described by Kapili and Dekas (2021). Finally, we use the Phylogenetic Placement for Inferring Taxonomy (PPIT v.1.2.0), to interpret the nifH phylogenetic placement and infer the taxonomy of the ASVs (Kapili and Dekas, 2021). For the purpose of this study, ASVs identified as suspected nifH homologs were removed, and only ASVs with a taxonomic inference were used to assess the diazotrophic community.
Data analysis and statistical methods
The data analysis and statistical methods described below were carried out using R (version 4.1.3). Following the modifications described by Raut et al. (2018), a linear regression model—excluding any lags and plateaus in nitrogenase activity at the beginning and end, respectively—of ARA incubations was fit to calculate rates of C2H4 production in nmol by converting peak heights of C2H4 measured throughout the ARA incubations (Capone, 1993; Breitbarth et al., 2004). Then, using a 3:1 ratio of C2H2 reduction:N2 reduction (Capone, 1993), theoretical BNF rates were estimated and multiplied by 2 to express as nmol N × g−1 (dw of sediment) × h−1 for all experiments.
Subsequently, to determine departures from normality for BNF rates, the Shapiro–Wilk test (“shapiro.test”) was used. To test for homogeneity of group variances in BNF rates, the Fligner-Killeen test (“fligner.test”), a non-parametric method based on ranks, was used. Since all experiments incorporated a two-factorial design, i.e., sites (three locations) and amendments (variable depending on the experiment), a two-way analysis of variance (ANOVA; “aov”) was performed, and if a significant p-value (< 0.05) was obtained, multiple pairwise comparisons were conducted between BNF rates of different groups using Tukey's honest significant difference (HSD) test (“TukeyHSD”). For multiple experiments, where the data distribution violated assumptions of normality and homoscedasticity, a non-parametric alternative of the two-way ANOVA (“t2way”) was carried out. For instances with a significant p-value (< 0.05), a post-hoc test similar to Tukey's HSD that conducts multiple pairwise comparisons using medians (“medpb”) was used to compare between independent groups of BNF rates (Wilcox, 2017). The non-parametric alternatives to the two-way ANOVA and Tukey's HSD are functions stored in the file “Rallfun-v40.txt” (January 2023) that is publicly available from https://osf.io/xhe8u/. The p-values from these comparisons that are statistically significant are defined in the different figures displaying BNF rates as follows: p-value < 0.001 (***), p-value < 0.01 (**), and p-value < 0.05 (*).
The Bray–Curtis dissimilarity metric was calculated using the % relative abundance of ASVs found in different macroalgae subsampled throughout long-term decomposition experiments using the “vegdist” function with default specifications. Non-metric multidimensional scaling (NMDS), which uses ranks orders of species and makes fewer assumptions of the dataset, was run using the Bray–Curtis distance matrix with the “metaMDS” function under some modifications from default settings (i.e., autotransform = F, wascores = T) to visualize the nifH-containing communities. A hierarchical cluster analysis was performed on the Bray–Curtis distances with the “hclust” function using the complete linkage method to identify similar clusters (k = 4). An analysis of similarities (ANOSIM) test was performed using the function “anosim” with default specifications for different factors of grouping observations (i.e., macroalgal species, hierarchical cluster, and stage of decomposition). A permutational multivariate analysis of variance (PERMANOVA) test was also performed using the function “Adonis” with default specifications for the previously specified factors of grouping observations. The “vegdist,” “metaMDS,” “anosim,” and “adonis” functions were all part of the vegan package (v.2.6-2).
Results and discussion
BNF rates at Big Fisherman's Cove (BFC) and the impact of labile carbon availability
BNF rates measured at all three sites in BFC were within the range of previously reported BNF rates throughout a diverse array of coastal marine habitats, albeit at the lower end (Table 1). These lower rates could potentially be due to a deviation from the theoretical ratio (3:1) used to estimate BNF from ARA as observed in other systems (Soper et al., 2021) and could be refined in future studies by employing the more sensitive stable isotope tracer (e.g., 15 N2) methodology. Nonetheless, the diazotrophic activity was consistently higher at sites 2 and 3 in comparison with site 1 (Figure 2). We hypothesize that this could be due to greater sedimentation of metabolizable organic matter at sites 2 and 3 since benthic diazotrophic activity has been linked to organic matter availability previously (Capone and Budin, 1982; McGlathery et al., 1998). Additionally, potentially greater variability in environmental factors near the intertidal zone might also influence the lower BNF rates measured at site 1. For instance, increased wave exposure, as observed in this site, has previously been reported to have an overall negative impact on benthic diazotrophic activity across multiple shallow-water coastal sites (Andersson et al., 2014). Furthermore, the decreased availability of labile organic carbon can have another deleterious effect on the broader heterotrophic bacterial communities, including diazotrophs, at sites with greater wave exposure (Sander and Kalff, 1993).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1312843/fmicb-14-1312843-HTML/image_m/fmicb-14-1312843-t001.jpg)
Table 1. Compilation of benthic biological nitrogen fixation (BNF) rates from coastal marine sediments around the world.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1312843/fmicb-14-1312843-HTML/image_m/fmicb-14-1312843-g002.jpg)
Figure 2. Biological nitrogen fixation (BNF) rates measured at the three different sampling sites in Big Fisherman's Cove, Santa Catalina Island, CA, USA during four separate experiments. The different experimental conditions are as follows: experiment 1—subsampled from intact sediment cores in-situ in May 2017, experiments 2 and 3—subsampled sediments from longterm (LT) mesocosms that were 1 (June 2017) and 2 (July 2017) months old, respectively, and experiment 4—established a new, short-term (ST) mesocosm (October 2017) that was subsampled shortly thereafter. The unamended BNF rates measured at each site during those four experiments are depicted in bar and whisker plots and statistically significant differences between sites are highlighted with brackets (***p-value < 0.001).
To test the potential for organic carbon to stimulate diazotroph activity, we conducted bottle incubations (experiment 4) with sediments from all three sites and added a labile carbon substrate. In previous studies, there have been a wide variety of carbon substrates (e.g., glucose, sucrose, and succinate) across a range of concentrations (1–10 mM) used to stimulate benthic BNF rates (Capone and Budin, 1982; McGlathery et al., 1998). Among the different carbon substrates tested, 10 mM glucose yielded the greatest stimulation (Capone and Budin, 1982), and therefore, the same concentration and substrate was selected in this study. Despite the baseline variability in unamended BNF rates among the three sites, the addition of glucose ubiquitously stimulated nitrogenase activity and resulted in comparable BNF rates between the sites (Figure 3). The pronounced 2–3 orders of magnitude increase in nitrogenase activity with 10 mM glucose addition compared to unamended BNF rates (p-value < 0.001) suggests that each site harbors a benthic diazotrophic community poised to respond to labile organic matter input (Figure 3). This highlights the BNF potential at each site, suggesting that the availability of labile organic matter can be a strong limiting factor for the benthic diazotrophic community at BFC. Furthermore, similar stimulation of nitrogenase activity with the addition of labile carbon substrates has been observed previously in many other benthic coastal habitats (Capone and Budin, 1982; McGlathery et al., 1998). The results from this study further support that labile carbon availability can strongly influence benthic diazotrophic activity.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1312843/fmicb-14-1312843-HTML/image_m/fmicb-14-1312843-g003.jpg)
Figure 3. Average biological nitrogen fixation (BNF) rates (+/- standard error) measured at the three different sampling sites in Big Fisherman's Cove, Santa Catalina Island, CA, USA during experiment four in Fall 2017 are depicted by the bar plots. The blue bar represents unamended BNF rates, orange bar represents BNF rates with addition of brown macroalga Dictyopteris, green bar represents BNF rates with addition of green macroalga Codium fragile, brown bar represents BNF rates with addition of brown macroalga Macrocystis pyrifera, and gray bar represents BNF rates with addition of 10 mM glucose. Only the statistically significant difference between unamended (treatment A) and macroalgal and glucose amendments (treatments B - E) are depicted on the plot (***p-value < 0.001).
Role of autochthonous organic carbon on benthic diazotrophic community at BFC
In this study, we simulate the deposition of brown (i.e., Dictyopteris and M. pyrifera), green (i.e., C. fragile), and red (i.e., A. taxiformis) macroalgae by adding autoclaved or freeze-dried living macroalgal tissue to the sediment slurry incubations. Due to the length of ARA incubations (~5 days), there was sufficient time for the healthy macroalgal tissue to undergo remineralization, likely entering the intermediate stages of decomposition when the highest BNF rates are observed (Hamersley et al., 2015; Raut et al., 2018; Raut and Capone, 2021). Indeed, the brown and green macroalgal amendments (experiments 2 and 4) significantly stimulated BNF rates by ~2 orders of magnitude compared to unamended BNF rates (p-value < 0.001) across all three sites (Figures 3, 4). The remineralization of these brown and green macroalgae can release a substantial source of metabolizable carbon substrates, e.g., glucose, galactose, arabinose, mannitol, and laminarin (Zimmerman and Kremer, 1986; Ciancia et al., 2007), likely eliciting a similar microbial metabolic response as observed with the addition of 10 mM glucose (Figure 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1312843/fmicb-14-1312843-HTML/image_m/fmicb-14-1312843-g004.jpg)
Figure 4. Average biological nitrogen fixation (BNF) rates (+/- standard error) measured at the three different sampling sites in Big Fisherman's Cove, Santa Catalina Island, CA, USA during experiment 2 (June 2017) are depicted by bar plots. The unamended BNF rates are in blue, BNF rates with addition of brown macroalga Macrocystis pyrifera in brown, and M. pyrifera plus 20 mM sodium molybdate, a strong inhibitor of sulfate reducers, additions are in gray. Only statistically significant differences are highlighted between the comparisons of (1) unamended incubation and two different amendments and (2) macroalgae only amendment vs. macroalgae + 20 mM sodium molybdate amendments: p-value < 0.001***.
Interestingly, the BNF rates between brown (i.e., Dictyopteris vs. M. pyrifera) and green (i.e., C. fragile) macroalgal amendments were notably different at sites 2 and 3 (p-value < 0.001; Figure 3). This potentially suggests that carbon loading of different macroalgal species, even within the same clade, might have different impacts on benthic diazotrophic activity (Figure 3). Furthermore, the ubiquitous response of nitrogenase activity across all three sites reinforces the notion that the benthic diazotrophic community is very strongly limited by the availability of labile organic matter. These results suggest that the autochthonous organic matter deposition of macroalgal biomass is an important source of labile carbon to the benthic diazotrophic community along coastal macroalgal habitats such as BFC (Figures 1, 3, 4).
In previous studies, the highest BNF rates measured during the intermediate phase of decomposition of several red macroalgae (e.g., Plocamium and A. taxiformis) were ~2 orders of magnitude lower than BNF rates associated with brown and green macroalgae (Raut and Capone, 2021). Therefore, the addition of the red macroalga A. taxiformis was not expected to stimulate benthic BNF rates at any of the three sites. Surprisingly, the A. taxiformis amendments (experiment 3) were found to inhibit diazotrophic activity, yielding significantly lower BNF rates compared to the unamended BNF rates (p-value < 0.001) measured at all three sites (Figure 5). Reduced BNF rates could potentially be linked to the notable differences in stoichiometric ratios of carbon to N between A. taxiformis (<7) and other brown and green macroalgae (>10) (Raut and Capone, 2021). Here, the lower carbon to N content of A. taxiformis might have resulted in insufficient release of labile carbon to stimulate the benthic diazotrophic activity.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1312843/fmicb-14-1312843-HTML/image_m/fmicb-14-1312843-g005.jpg)
Figure 5. Biological nitrogen fixation (BNF) rates measured at the three different sampling sites in Big Fisherman's Cove, Santa Catalina Island, CA, USA during experiment 3 (July 2017). The unamended BNF rates (blue) and BNF rates with addition of red macroalga Asparagopsis taxiformis (magenta) are depicted using a box and whisker plot with dashed line representing the mean BNF rates for each treatment. Statistically significant values: p-value < 0.05*, p-value < 0.001***.
The release of halogenated secondary metabolites (e.g., bromophenol), which have strong antimicrobial activity (Mata et al., 2012), might have also prohibited microbial colonization during remineralization of A. taxiformis resulting in reduced diazotrophic activity. Additionally, A. taxiformis can be a significant source of phenolic compounds (Gao et al., 2022), which have been shown in other anoxic sediment systems, e.g., seagrass rhizosphere, to inhibit heterotrophic microbial consumption of sugars, e.g., sucrose (Sogin et al., 2022). Interestingly, the few microbes identified with the genes to degrade phenols and utilize sucrose within the seagrass rhizosphere were sulfate reducers and N2 fixers (Sogin et al., 2022). The decomposition of A. taxiformis could be selecting for fewer heterotrophic N2 fixers with the ability to degrade phenolic compounds and utilize the available sugars from the macroalgal detritus, resulting in lower bulk BNF rates than observed with remineralization of other macroalgae. Ultimately, this response by the diazotrophic community indicates that not all macroalgal deposition to the benthos will necessarily stimulate BNF rates.
Taken together, this suggests that the species composition of macroalgal sedimentation will play an important role in dictating the response by the benthic diazotrophic community. In fact, the increased prevalence of A. taxiformis throughout site 1 might further explain the consistently lower unamended BNF rates observed here in comparison with sites 2 and 3 (p-value < 0.001) in all four experiments (Figure 2). The chemical composition of the autochthonous source of carbon, e.g., macroalgal detritus, likely plays a significant role in determining the response by the benthic diazotrophic community to organic matter loading occurring along coastal marine habitats. Since significantly different BNF rates have been observed throughout the three stages of macroalgal decomposition (Hamersley et al., 2015; Raut et al., 2018; Raut and Capone, 2021), future studies should also explore how long-term macroalgal burial and remineralization in the marine sediment might impact benthic BNF rates. It is likely that the quality and lability of macroalgal detritus will influence benthic BNF rates, with the addition of more recalcitrant macroalgal detritus likely resulting in lower BNF rates compared to additions of healthier macroalgal tissue as used in this study.
Diversity of macroalgae-associated diazotrophic community during decomposition
The beta diversity of the diazotrophic communities was visualized on a non-metric multidimensional scaling (NMDS) plot using the Bray–Curtis dissimilarity indices and shows a clear distinction of diazotrophic communities associated with the different macroalgal species (Figure 6). An ANOSIM test was used to determine that the diazotrophic ASV assemblages were highly dissimilar between the three different species of macroalgae (R-value = 0.75, p-value = 0.001). A hierarchical cluster analysis was also used to identify four major clusters of samples (Figure 6), and an ANOSIM test similarly determined the diazotrophic assemblages to be highly dissimilar between the four clusters (R-value = 0.84, p-value = 0.001). A PERMANOVA test identified the hierarchical clusters (R2 = 0.14, p-value = 0.005) to explain a lower percentage of variation (14%) in diazotrophic communities across samples than macroalgal species composition (R2 = 0.38, p-value = 0.001, 38%; Figure 6).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1312843/fmicb-14-1312843-HTML/image_m/fmicb-14-1312843-g006.jpg)
Figure 6. Non-metric dimensional scaling (NMDS) plot of Bray-Curtis dissimilarity indices of diazotrophic communities associated with different macroalgal species (Asparagopsis taxiformis—red macroalga; Codium fragile—green macroalga; Dictyopteris—brown macroalga) subsampled at different stages of decomposition: early, middle, and end. The stage of decomposition is denoted in the label next to each sample along with the exact day of sub-sampling from longterm litter bag decomposition. Hierarchical cluster analysis was also used to identify four major clusters denoted by the different shapes. The analysis of similarities (ANOSIM) test was used to determine how similar (closer to 0) or dissimilar (closer to 1) diazotrophic communities were to each other based on 3 different grouping variables: (1) macroalgal species composition, (2) hierarchical clusters, and (3) stage of decomposition of the macroalgal detritus. The permutational multivariate analysis of variance (PERMANOVA) test was used to further determine the percentage (R2 value) of variation explained by the 3 corresponding grouping variables.
Notably, there is an emerging feature with macroalgal samples where the diazotrophic community associated with the earlier stages of decomposition is distinct and partitions separately from samples in the intermediate and latter stages of decomposition (Figure 6). Indeed, categorical assignment of the three stages of decomposition, i.e., early, middle, and end, based on observation of the physical and structural integrity of each distinct macroalgal sample at the time of collection revealed that diazotrophic communities were very similar within each stage of decomposition (ANOSIM: R-value = 0.18, p-value = 0.001). The phase of decomposition was found to explain 20% (PERMANOVA: R2 = 0.20, p-value = 0.01) of the variation in diazotrophic communities across samples. Collectively, these observations suggest that the species composition and stage of decomposition of the macroalgal substrate likely play an important role in structuring the associated diazotrophic community.
For each macroalgal sample, microbial lineages inferred to host nifH that collectively contributed >0.5 % relative abundance were aggregated at the class or phylum level based on the inferred taxonomy using PPIT (Kapili and Dekas, 2021) and visualized in a stacked bar plot (Figures 7–9). The macroalgal-associated nifH-containing assemblages form distinct communities throughout the decomposition of brown, green, and red macroalgae, notably on the initial day of sampling where a higher abundance of cyanobacterial-nifH ASVs within Oscillatoriophycideae are observed (Figures 7–9). This is followed by a shift to nifH sequences associated with heterotrophic bacteria (Figure 8). This includes the ubiquitous presence of nifH sequences inferred to be from Gammaproteobacteria and Deltaproteobacteria (Figures 7–9).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1312843/fmicb-14-1312843-HTML/image_m/fmicb-14-1312843-g007.jpg)
Figure 7. Stacked bar plots depicting the taxonomic composition of major nifH-containing amplicon sequence variants (ASVs) aggregated at the class or phylum level associated with macroalgal detritus sub-sampled on different days (i.e., 0, 8, 10, 14) throughout a long-term (16 days) litter bag decomposition of the brown macroalga, Dictyopteris. The remaining ASVs contributing <0.5% relative abundance of the total community per sample are aggregated together into a distinct group and can consist of a wide array of taxonomic diversity.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1312843/fmicb-14-1312843-HTML/image_m/fmicb-14-1312843-g008.jpg)
Figure 8. Stacked bar plots depicting the taxonomic composition of major nifH-containing amplicon sequence variants (ASVs) aggregated at the class or phylum level associated with macroalgal detritus sub-sampled on different days [i.e., 0, 2, 4, 7, 10 (technical sequencing replicate of day 10–10.2), 13, 17, 22, 26, 32] throughout a long-term (32 days) litter bag decomposition of the green macroalga, Codium fragile. The remaining ASVs contributing <0.5% relative abundance of the total community per sample are aggregated together into a distinct group and can consist of a wide array of taxonomic diversity.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1312843/fmicb-14-1312843-HTML/image_m/fmicb-14-1312843-g009.jpg)
Figure 9. Stacked bar plots depicting the taxonomic composition of major nifH-containing amplicon sequence variants (ASVs) aggregated at the class or phylum level associated with macroalgal detritus sub-sampled on different days (i.e., 0, 5, 12, 15, 18) throughout a long-term (18 days) litter bag decomposition of the red macroalga, Asparagopsis taxiformis. The remaining ASVs contributing <0.5% relative abundance of the total community per sample are aggregated together into a distinct group and can consist of a wide array of taxonomic diversity.
For C. fragile, Dictyopteris, and A. taxiformis, the most prominent nifH-containing Gammaproteobacteria was from the genus Vibrio, where at times, a single nifH ASV could contribute up to 29, 13, and 19%, respectively, of the total nifH-containing community composition (Day 0). Furthermore, there were also ASVs with putative hosts within a wide array of orders including Vibrionales, Alteromonadales, Aeromonadales, Methylococcales, Oceanospirillales, Thiotrichales, and Chromatiales within the Gammaproteobacteria class. The most prevalent inferred nifH hosts within the Deltaproteobacteria were from the order Desulfovibrionales or Desulfuromonadales, with the most prevalent ASV contributing up to 57 (Day 26), 25 (Day 10), and 12% (Day 15) of the nifH-containing community composition associated with C. fragile, Dictyopteris, and A. taxiformis, respectively. There was also a wide array of ASVs from different families within the Deltaproteobacteria class including Desulfovibrionaceae, Desulfuromonadaceae, Desulfobulbaceae, Desulfobacteraceae, and Geobacteraceae. Similar diazotrophic members have also previously been observed with nifH surveys in coastal sediments (Jabir et al., 2018), including ecosystems populated by other macrophytes such as mangroves (Zhang et al., 2017) and seagrasses (Sun et al., 2015).
Collectively, there is a successional shift within the nifH-containing community from Gammaproteobacteria-dominated ASVs at the earlier stages of decomposition toward a more anaerobic assemblage of ASVs identified as from Deltaproteobacteria, Spirochaeta, and Bacteroidia in the intermediate and latter stages of remineralization (Figures 7–9). Members of the class Spirochaetia, including those within the genus Spirochaeta, are also inferred to host nifH sequences in these samples and are particularly prevalent during the intermediate and latter stages of macroalgal decomposition of all three brown, red, and green seaweeds (Figures 7–9). Spirochete bacteria, including N2 fixers, are often observed in anoxic aquatic environments, e.g., coastal marine habitats (Shivani et al., 2015), oxygen minimum zones (Loescher et al., 2014), and surface marine sediments (Jabir et al., 2021). After deposition onto marine sediments, it is likely that macroalgal decomposition fosters the development of anoxic conditions that would be favorable for similar Spirochaetia ASVs as those observed with macroalgal-associated detritus throughout long-term litter bag decomposition (Figures 7–9). There is also a prevalence of Bacteroidia associated with the brown macroalga, Dictyopteris, and the ASVs within this class consist of different families including Prolixibacteraceae and Marinilabiliaceae (Figure 7). Facultatively anaerobic Draconibacterium strains isolated from marine sediments (Du et al., 2014) and N2 fixing strains of Mangrovibacterium isolated from mangrove sediment (Huang et al., 2014) are both within the family Prolixibacteraceae which are the prominent Bacteroidia ASVs observed in this study. Furthermore, other facultatively anaerobic N2 fixing organisms within the Bacteroidia class have also been isolated from marine sediments with macrophyte (e.g., cordgrass and mangrove) presence (Huang et al., 2020).
Notably, the nifH-containing community associated with A. taxiformis was distinct from the community associated with the brown and green macroalgae, Dictyopteris and C. fragile (Figure 9). During the community succession throughout A. taxiformis remineralization, there is also a shift from Gammaproteobacteria during the earlier stages of decomposition to a Deltaproteobacteria-dominated nifH-containing community in the intermediate and latter stages. However, this transition is not marked by a homogenous shift toward a Deltaproteobacteria community but rather also consists of equal or greater contribution by Bacteroidia and Clostridia nifH-containing ASVs (Figure 9). Additional anaerobic fermentative Firmicutes species within the Clostridia class, including several representatives from the families Defluviitaleaceae, Lachnospiraceae, and Clostridiaceae, were also observed to be in high abundance toward the latter stages of A. taxiformis degradation (Figure 9). These findings suggest that while Gammaproteobacteria and Deltaproteobacteria are the dominant members of the macroalgal-associated nifH-containing communities, they may also be accompanied by a more diverse assemblage of Spirochaeta, Bacteroidia, and Clostridia N2 fixers (Figures 7–9).
Since nifH sequences inferred to be from Gammaproteobacteria were observed throughout the earlier (<5 days) stage of macroalgal decomposition (Figures 7–9), we hypothesize that macroalgal additions during the ~5-day ARA of sediment slurry incubations might stimulate the resident Gammaproteobacteria members of the benthic diazotrophic community. In this study, the shorter incubation and remineralization time period would likely not provide sufficient time for the succession to a Deltaproteobacteria-dominated community consisting of sulfate-reducing phylotypes (Jabir et al., 2021). Consistent with this, incubations amended with a macroalgal carbon source and 20 mM sodium molybdate, a potent inhibitor of sulfate reduction, only resulted in a significant reduction of BNF rates (p-value < 0.001) compared to incubations with only the macroalgal carbon source at site 1 (Figure 4). We hypothesize that this might be due to an existing difference in benthic diazotrophic assemblages between the three sites, perhaps with greater sulfate reducers present at site 1. In the future, nifH sequencing of marine sediment from the different sites could be helpful in characterizing nifH-containing community and deconvoluting potential differences in the microbial community between sites. Additionally, a longer time series could be used to directly monitor the succession of benthic diazotrophic assemblages throughout macroalgal remineralization, especially to investigate whether a similar shift from Gammaproteobacteria toward a Deltaproteobacteria-dominated community occurs as observed with the macroalgal-associated nifH-containing community.
Conclusion
Conventionally, BNF is often thought to occur where N is limiting and thus is often under-surveyed in coastal marine sediments, especially macroalgal ecosystems where high macroalgal productivity (e.g., bloom events) could serve as an important source of organic matter. This study reevaluates that paradigm using biogeochemical measurements concurrently with molecular tools to disentangle the systematic impact that remineralization of macroalgae has on benthic BNF rates and characterizes the nifH-containing microbes driving this process in coastal marine ecosystems. Our study suggests coastal benthic diazotrophs are limited by organic carbon and demonstrates a significant and phylum-specific effect of macroalgal loading on benthic microbial communities.
Data availability statement
The original contributions presented in the study are deposited in the European Nucleotide Archive (ENA) repository under accession number PRJEB71000 (secondary accession: ERP155889).
Ethics statement
Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author contributions
YR: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Visualization, Writing—original draft, Writing—review & editing. CRB: Conceptualization, Investigation, Methodology, Writing—review & editing. ERP: Data curation, Formal analysis, Methodology, Writing—review & editing. BJK: Data curation, Formal analysis, Methodology, Software, Writing—review & editing. AED: Funding acquisition, Resources, Supervision, Writing—review & editing. DGC: Funding acquisition, Resources, Supervision, Writing—review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. YR was supported to do this study at WIES through the summer Wrigley/Victoria J. Bertics Fellowships (2017-2019). This work was partially supported by Stanford University through research and graduate student funding to AED.
Acknowledgments
The authors would like to thank the USC Wrigley Institute for Environment and Sustainability (WIES) and their staff, especially Lauren Czarnecki Oudin and Kellie Spafford, for their accommodations and support for this research. The authors thank previous Research Experiences for Undergraduates (NSF Award OCE-1559941) participants for their help: Taylor Dillon, Calyn Crawford, and Sarah Ortiz. The authors thank Jasmine Henderson for logistical support and help with the fieldwork. The authors also appreciate the two reviewers' feedback and suggestions to help improve the manuscript.
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.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Aires, T., Muyzer, G., Serrão, E. A., and Engelen, A. H. (2019). Seaweed loads cause stronger bacterial community shifts in coastal lagoon sediments than nutrient loads. Front. Microbiol. 10, 3283. doi: 10.3389/fmicb.2018.03283
Andersson, B., Sundbäck, K., Hellman, M., Hallin, S., and Alsterberg, C. (2014). Nitrogen fixation in shallow-water sediments: spatial distribution and controlling factors. Limnol. Oceanogr. 59, 1932–1944. doi: 10.4319/lo.2014.59.6.1932
Bertics, V. J., Löscher, C. R., Salonen, I., Dale, A. W., Gier, J., Schmitz, R. A., et al. (2013). Occurrence of benthic microbial nitrogen fixation coupled to sulfate reduction in the seasonally hypoxic Eckernförde Bay, Baltic Sea. Biogeosciences 10, 1243–1258. doi: 10.5194/bg-10-1243-2013
Bertics, V. J., Sohm, J. A., Treude, T., Chow, C. E. T., Capone, D. G., Fuhrman, J. A., et al. (2010). Burrowing deeper into benthic nitrogen cycling: the impact of bioturbation on nitrogen fixation coupled to sulfate reduction. Mar. Ecol. Prog. Ser. 409, 1–15. doi: 10.3354/meps08639
Breitbarth, E., Mills, M. M., Friedrichs, G., and Laroche, J. (2004). The Bunsen gas solubility coefficient of ethylene as a function of temperature and salinity and its importance for nitrogen fixation assays. Limnol. Oceanogr. Methods 2, 282–288. doi: 10.4319/lom.2004.2.282
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869
Capone, D. G. (1982). Nitrogen fixation(acetylene reduction) by rhizosphere sediments of the eelgrass Zostera marina. Mar. Ecol. Prog. Ser. 10, 67–75. doi: 10.3354/meps010067
Capone, D. G. (1993). “Determination of nitrogenase activity in aquatic samples using the acetylene reduction procedure,” in Handbook of Methods in Aquatic Microbial Ecology (Boca Raton, FL: Lewis), 621–631.
Capone, D. G., and Budin, J. M. (1982). Nitrogen fixation associated with rinsed roots and rhizomes of the eelgrass Zostera marina. Plant Physiol. 70, 1601–1604. doi: 10.1104/pp.70.6.1601
Capone, D. G., Taylor, D. L., and Taylor, B. F. (1977). Nitrogen fixation (acetylene reduction) associated with macroalgae in a coral-reef community in the Bahamas. Mar. Biol. 40, 29–32. doi: 10.1007/BF00390624
Ciancia, M., Quintana, I., Vizcargüénaga, M. I., Kasulin, L., de Dios, A., Estevez, J. M., et al. (2007). Polysaccharides from the green seaweeds Codium fragile and C. vermilara with controversial effects on hemostasis. Int. J. Biol. Macromol. 41, 641–649. doi: 10.1016/j.ijbiomac.2007.08.007
Dekas, A. E., Connon, S. A., Chadwick, G. L., Trembath-Reichert, E., and Orphan, V. J. (2016). Activity and interactions of methane seep microorganisms assessed by parallel transcription and FISH-NanoSIMS analyses. ISME J. 10, 678–692. doi: 10.1038/ismej.2015.145
Dekas, A. E., Fike, D. A., Chadwick, G. L., Green-Saxena, A., Fortney, J., Connon, S. A., et al. (2018). Widespread nitrogen fixation in sediments from diverse deep-sea sites of elevated carbon loading. Environ. Microbiol. 20, 4281–4296. doi: 10.1111/1462-2920.14342
Dekas, A. E., Poretsky, R. S., and Orphan, V. J. (2009). Deep-sea archaea fix and share nitrogen in methane-consuming microbial consortia. Science 326, 422–426. doi: 10.1126/science.1178223
Dierssen, H. M., Zimmerman, R. C., Drake, L. A., and Burdige, D. J. (2009). Potential export of unattached benthic macroalgae to the deep sea through wind-driven Langmuir circulation. Geophys. Res. Lett. 36, 188. doi: 10.1029/2008GL036188
Du, Z. J., Wang, Y., Dunlap, C., Rooney, A. P., and Chen, G. J. (2014). Draconibacterium orientale gen. nov., sp. nov., isolated from two distinct marine environments, and proposal of Draconibacteriaceae fam. nov. Int. J. Syst. Evol. Microbiol. 64, 1690–1696. doi: 10.1099/ijs.0.056812-0
Duarte, C. M. (1992). Nutrient concentration of aquatic plants: patterns across species. Limnol. Oceanogr. 37, 882–889. doi: 10.4319/lo.1992.37.4.0882
Duarte, C. M., and Cebrián, J. (1996). The fate of marine autotrophic production. Limnol. Oceanogr. 41, 1758–1766. doi: 10.4319/lo.1996.41.8.1758
Flett, R. J., Hamilton, R. D., and Campbell, N. E. R. (1976). Aquatic acetylene-reduction techniques: solutions to several problems. Can. J. Microbiol. 22, 43–51. doi: 10.1139/m76-006
Fulweiler, R. W. (2023). More foxes than hedgehogs: the case for nitrogen fixation in coastal marine sediments. Glob. Biogeochem. Cycl. 37, 7777. doi: 10.1029/2023GB007777
Fulweiler, R. W., Heiss, E. M., Rogener, M. K., Newell, S. E., LeCleir, G. R., Kortebein, S. M., et al. (2015). Examining the impact of acetylene on N-fixation and the active sediment microbial community. Nat. Inst. Health 6, 1–9. doi: 10.3389/fmicb.2015.00418
Gandy, E. L., and Yoch, D. C. (1988). Relationship between nitrogen-fixing sulfate reducers and fermenters in salt marsh sediments and roots of Spartina alterniflora. Appl. Environ. Microbiol. 54, 2031–2036. doi: 10.1128/aem.54.8.2031-2036.1988
Gao, H., Wang, Y., Guo, Z., Liu, Y., Wu, Q., and Xiao, J. (2022). Optimization of ultrasound-assisted extraction of phenolics from Asparagopsis taxiformis with deep eutectic solvent and their characterization by ultra-high-performance liquid chromatography-mass spectrometry. Front. Nutr. 9, 1036436. doi: 10.3389/fnut.2022.1036436
Gier, J., Sommer, S., Löscher, C. R., Dale, A. W., Schmitz, R. A., and Treude, T. (2016). Nitrogen fixation in sediments along a depth transect through the Peruvian oxygen minimum zone. Biogeosciences 13, 4065–4080. doi: 10.5194/bg-13-4065-2016
Ginsburg, D. W., and Huang, A. H. (2022). Over, under, sideways and down: patterns of marine species richness in nearshore habitats off Santa Catalina Island, California. Diversity 14, d14050366. doi: 10.3390/d14050366
Hamersley, M. R., Sohm, J. A., Burns, J. A., and Capone, D. G. (2015). Nitrogen fixation associated with the decomposition of the giant kelp Macrocystis pyrifera. Aquat. Bot. 125, 57–63. doi: 10.1016/j.aquabot.2015.05.003
Head, W. D., and Carpenter, E. J. (1975). Nitrogen fixation associated with the marine macroalga Codium fragile. Limnol. Oceanogr. 20, 815–823. doi: 10.4319/lo.1975.20.5.0815
Hernández-López, J., Camacho-Ibar, V. F., Macías-Tapia, A., McGlathery, K. J., Daesslé, L. W., and Sandoval-Gil, J. M. (2017). Benthic nitrogen fixation in Zostera marina meadows in an upwelling-influenced coastal lagoon. Ciencias Marinas 43, 35–53. doi: 10.7773/cm.v43i1.2700
Hou, L., Wang, R., Yin, G., Liu, M., and Zheng, Y. (2018). Nitrogen fixation in the intertidal sediments of the Yangtze Estuary: occurrence and environmental implications. J. Geophys. Res. Biogeosci. 123, 936–944. doi: 10.1002/2018JG004418
Howarth, R. W., Marino, R., and Cole, J. J. (1988). Nitrogen fixation in freshwater, estuarine, and marine ecosystems. 2. Biogeochemical controls. Limnol. Oceanogr. 33, 688–701. doi: 10.4319/lo.1988.33.4_part_2.0688
Huang, X. F., Liu, Y. J., Dong, J., De Qu, L. Y., Zhang, Y. Y., Wang, F. Z., et al. (2014). Mangrovibacterium diazotrophicum gen. nov., sp. nov., a nitrogen-fixing bacterium isolated from a mangrove sediment, and proposal of Prolixibacteraceae fam. nov. Int. J. Syst. Evol. Microbiol. 64, 875–881. doi: 10.1099/ijs.0.052779-0
Huang, Z., Hu, Y., Lai, Q., and Guo, Y. (2020). Description of Maribellus sediminis sp. nov., a marine nitrogen-fixing bacterium isolated from sediment of cordgrass and mangrove. Syst. Appl. Microbiol. 43, 126099. doi: 10.1016/j.syapm.2020.126099
Jabir, T., Jesmi, Y., Vipindas, P. V., and Mohamed Hatha, A. (2018). Diversity of nitrogen fixing bacterial communities in the coastal sediments of southeastern Arabian Sea (SEAS). Deep Sea Res 2 Top Stud. Oceanogr. 156, 51–59. doi: 10.1016/j.dsr2.2018.09.010
Jabir, T., Vipindas, P. V., Jesmi, Y., Divya, P. S., Adarsh, B. M., Nafeesathul Miziriya, H. S., et al. (2021). Influence of environmental factors on benthic nitrogen fixation and role of sulfur reducing diazotrophs in a eutrophic tropical estuary. Mar. Pollut. Bull. 165, 112126. doi: 10.1016/j.marpolbul.2021.112126
Johnson, C., Dubbs, L. L., and Piehler, M. (2023). Reframing the contribution of pelagic Sargassum epiphytic N2 fixation. PLoS ONE 18, 289485. doi: 10.1371/journal.pone.0289485
Kapili, B. J., Barnett, S. E., Buckley, D. H., and Dekas, A. E. (2020). Evidence for phylogenetically and catabolically diverse active diazotrophs in deep-sea sediment. ISME J. 14, 971–983. doi: 10.1038/s41396-019-0584-8
Kapili, B. J., and Dekas, A. E. (2021). PPIT: an R package for inferring microbial taxonomy from nifH sequences. Bioinformatics 37, 2289–2298. doi: 10.1093/bioinformatics/btab100
Krause-Jensen, D., and Duarte, C. M. (2016). Substantial role of macroalgae in marine carbon sequestration. Nat. Geosci. 9, 737–742. doi: 10.1038/ngeo2790
Lam, B. R., Barr, C. R., Rowe, A. R., and Nealson, K. H. (2019). Differences in applied redox potential on cathodes enrich for diverse electrochemically active microbial isolates from a marine sediment. Front. Microbiol. 10, 1979. doi: 10.3389/fmicb.2019.01979
Lapointe, B. E., Burkholder, J. M., and Van Alstyne, K. L. (2018). “Harmful macroalgal blooms in a changing world: causes, impacts, and management,” in Harmful Algal Blooms (New York, NY: Wiley), 515–560. doi: 10.1002/9781118994672.ch15
Liesirova, T., Aarenstrup-Launbjerg, T., Hallstrøm, S., Bittner, M. J., Riemann, L., and Voss, M. (2023). Nitrogen-fixing sulfate reducing bacteria in shallow coastal sediments under simulated resuspension. Estuar. Coast Shelf Sci. 280, 108165. doi: 10.1016/j.ecss.2022.108165
Loescher, C. R., Großkopf, T., Desai, F. D., Gill, D., Schunck, H., Croot, P. L., et al. (2014). Facets of diazotrophy in the oxygen minimum zone waters off Peru. ISME J. 8, 2180–2192. doi: 10.1038/ismej.2014.71
Looby, A., and Ginsburg, D. W. (2021). Nearshore species biodiversity of a marine protected area off Santa Catalina Island, California. West N Am. Nat. 81, 113–130. doi: 10.3398/064.081.0110
Luo, Z., Zhong, Q., Han, X., Hu, R., Liu, X., Xu, W., et al. (2021). Depth-dependent variability of biological nitrogen fixation and diazotrophic communities in mangrove sediments. Microbiome 9, 1164. doi: 10.1186/s40168-021-01164-0
Lyons, D. A., Arvanitidis, C., Blight, A. J., Chatzinikolaou, E., Guy-Haim, T., Kotta, J., et al. (2014). Macroalgal blooms alter community structure and primary productivity in marine ecosystems. Glob. Chang Biol. 20, 2712–2724. doi: 10.1111/gcb.12644
Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 17, 10–12. doi: 10.14806/ej.17.1.200
Mata, L., Gaspar, H., and Santos, R. (2012). Carbon/nutrient balance in relation to biomass production and halogenated compound content in the red alga asparagopsis taxiformis (bonnemaisoniaceae). J. Phycol. 48, 248–253. doi: 10.1111/j.1529-8817.2011.01083.x
McGlathery, K. J., Risgaard-Petersen, N., and Christensen, P. B. (1998). Temporal and spatial variation in nitrogen fixation activity in the eelgrass Zostera marina rhizosphere. Mar. Ecol. Prog. Ser. 168, 245–258. doi: 10.3354/meps168245
Mehta, M. P., Butterfield, D. A., and Baross, J. A. (2003). Phylogenetic diversity of nitrogenase (nifH) genes in deep-sea and hydrothermal vent environments of the Juan de Fuca Ridge. Appl. Environ. Microbiol. 69, 960–970. doi: 10.1128/AEM.69.2.960-970.2003
Meyer-Reil, L. A. (1987). Seasonal and spatial distribution of extracellular enzymatic activities and microbial incorporation of dissolved organic substrates in marine sediments. Appl. Environ. Microbiol. 53, 1748–1755. doi: 10.1128/aem.53.8.1748-1755.1987
Mirarab, S., Nguyen, N., and Warnow, T. (2012). SEPP: SATé-enabled phylogenetic placement. Biocomputing 2012, 247–258. doi: 10.1142/9789814366496_0024
Mohr, W., Lehnen, N., Ahmerkamp, S., Marchant, H. K., Graf, J. S., Tschitschko, B., et al. (2021). Terrestrial-type nitrogen-fixing symbiosis between seagrass and a marine bacterium. Nature 600, 105–109. doi: 10.1038/s41586-021-04063-4
Nedwell, D. B., and Azni bin Abdul Aziz, S. (1980). Heterotrophic nitrogen fixation in an intertidal saltmarsh sediment. Estuar. Coast. Mar. Sci. 10, 699–702. doi: 10.1016/S0302-3524(80)80097-1
Paerl, H. W., Fitzpatrick, M., and Bebout, B. M. (1996). Seasonal nitrogen fixation dynamics in a marine microbial mat: potential roles of cyanobacteria and microheterotrophs. Limnol. Oceanogr. 41, 419–427. doi: 10.4319/lo.1996.41.3.0419
Polyak, Y., Shigaeva, T., Gubelit, Y., Bakina, L., Kudryavtseva, V., and Polyak, M. (2017). Sediment microbial activity and its relation to environmental variables along the eastern Gulf of Finland coastline. J. Mar. Syst. 171, 101–110. doi: 10.1016/j.jmarsys.2016.11.017
Raut, Y., and Capone, D. G. (2021). Macroalgal detrital systems: an overlooked ecological niche for heterotrophic nitrogen fixation. Environ. Microbiol. 23, 4372–4388. doi: 10.1111/1462-2920.15622
Raut, Y., Morando, M., and Capone, D. G. (2018). Diazotrophic macroalgal associations with living and decomposing sargassum. Front. Microbiol. 9, 1–14. doi: 10.3389/fmicb.2018.03127
Rice, D. L. (1982). The detritus nitrogen problem: new observations and perspectives from organic geochemistry. Mar. Ecol. Prog. Ser. 9, 153–162. doi: 10.3354/meps009153
Rice, D. L., and Hanson, R. B. (1984). A kinetic model for detritus nitrogen: role of the associated bacteria in nitrogen accumulation. Bull. Mar. Sci. 35, 326–340.
Sander, B. C., and Kalff, J. (1993). Factors controlling bacterial production in marine and freshwater sediments. Microb. Ecol. 26, 79–99. doi: 10.1007/BF00177045
Shivani, Y., Subhash, Y., Tushar, L., Sasikala, C., and Ramana, C. V. (2015). Spirochaeta lutea sp. nov., isolated from marine habitats and emended description of the genus Spirochaeta. Syst. Appl. Microbiol. 38, 110–114. doi: 10.1016/j.syapm.2014.11.002
Sogin, E. M., Michellod, D., Gruber-Vodicka, H. R., Bourceau, P., Geier, B., Meier, D. V., et al. (2022). Sugars dominate the seagrass rhizosphere. Nat. Ecol. Evol. 6, 866–877. doi: 10.1038/s41559-022-01740-z
Sohm, J. A., Mahaffey, C., and Capone, D. G. (2008). Assessment of relative phosphorus limitation of Trichodesmium spp. in the North Pacific, North Atlantic, and the north coast of Australia. Limnol. Oceanogr. 53, 2495–2502. doi: 10.4319/lo.2008.53.6.2495
Sohm, J. A., Webb, E. A., and Capone, D. G. (2011). Emerging patterns of marine nitrogen fixation. Nat. Rev. Microbiol. 9, 499–508. doi: 10.1038/nrmicro2594
Soper, F. M., Simon, C., and Jauss, V. (2021). Measuring nitrogen fixation by the acetylene reduction assay (ARA): is 3 the magic ratio? Biogeochemistry 152, 345–351. doi: 10.1007/s10533-021-00761-3
Spinette, R. F., Brown, S. M., Ehrlich, A. L., Puggioni, G., Deacutis, C., and Jenkins, B. D. (2019). Diazotroph activity in surface Narragansett Bay sediments in summer is stimulated by hypoxia and organic matter delivery. Mar. Ecol. Progr. Ser. 614, 35–50. doi: 10.3354/meps12901
Sun, F., Zhang, X., Zhang, Q., Liu, F., Zhang, J., and Gong, J. (2015). Seagrass (Zostera marina) colonization promotes the accumulation of diazotrophic bacteria and alters the relative abundances of specific bacterial lineages involved in benthic carbon and sulfur cycling. Appl. Environ. Microbiol. 81, 6901–6914. doi: 10.1128/AEM.01382-15
Tanner, R. L., Obaza, A. K., and Ginsburg, D. W. (2019). Secondary Production of Kelp Bass Paralabrax clathratus in Relation to Coastal Eelgrass Zostera marina Habitat in a Southern California Marine Protected Area. Available online at: http://meridian.allenpress.com/scasbulletin/article-pdf/118/3/158/2967572/i2162-4534-118-3-158.pdf (accessed August 9, 2023).
Tenore, K. R., Hanson, R. B., Dornseif, B. E., and Wiederhold, C. N. (1979). The effect of organic nitrogen supplement on the utilization of different sources of detritus. Limnol. Oceanogr. 24, 350–355. doi: 10.4319/lo.1979.24.2.0350
Wang, H., Wang, G., and Gu, W. (2020). Macroalgal blooms caused by marine nutrient changes resulting from human activities. J. Appl. Ecol. 57, 766–776. doi: 10.1111/1365-2664.13587
Welsh, D. T. (2000). Nitrogen fixation in seagrass meadows: regulation, plant-bacteria interactions and significance to primary productivity. Ecol. Lett. 3, 58–71. doi: 10.1046/j.1461-0248.2000.00111.x
Wilcox, R. (2017). Modern Statistics for the Social and Behavioral Sciences, 2nd Edn. New York, NY: Chapman and Hall/CRC.
Yao, C., Zhang, Q., Lin, X., Zhao, J., and Zhang, X. (2023). Seasonal hypoxia enhances benthic nitrogen fixation and shapes specific diazotrophic community in the eutrophic marine ranch. Processes 11, 10138. doi: 10.3390/pr11010138
Zehr, J. P., and Capone, D. G. (2020). Changing perspectives in marine nitrogen fixation. Science 368, aay9514. doi: 10.1126/science.aay9514
Zehr, J. P., and Capone, D. G. (2023). Unsolved mysteries in marine nitrogen fixation. Trends Microbiol. 8, 4. doi: 10.1016/j.tim.2023.08.004
Zehr, J. P., Mellon, M., Braun, S., Litaker, W., Steppe, T., and Paerl, H. W. (1995). Diversity of heterotrophic nitrogen fixation genes in a marine cyanobacterial mat. Appl. Environ. Microbiol. 61, 2527–2532. doi: 10.1128/aem.61.7.2527-2532.1995
Zhang, X., Song, Y., Liu, D., Keesing, J. K., and Gong, J. (2015). Macroalgal blooms favor heterotrophic diazotrophic bacteria in nitrogen-rich and phosphorus-limited coastal surface waters in the Yellow Sea. Estuar. Coast. Shelf Sci. 163, 75–81. doi: 10.1016/j.ecss.2014.12.015
Zhang, Y., Yang, Q., Ling, J., Van Nostrand, J. D., Shi, Z., Zhou, J., et al. (2017). Diversity and structure of diazotrophic communities in mangrove rhizosphere, revealed by high-throughput sequencing. Front. Microbiol. 8, 2032. doi: 10.3389/fmicb.2017.02032
Keywords: benthic N2 fixation, macroalgal decomposition, diazotroph community, coastal marine sediments, nifH amplicon sequencing
Citation: Raut Y, Barr CR, Paris ER, Kapili BJ, Dekas AE and Capone DG (2024) Autochthonous carbon loading of macroalgae stimulates benthic biological nitrogen fixation rates in shallow coastal marine sediments. Front. Microbiol. 14:1312843. doi: 10.3389/fmicb.2023.1312843
Received: 10 October 2023; Accepted: 23 November 2023;
Published: 05 January 2024.
Edited by:
Xianbiao Lin, Ocean University of China, ChinaReviewed by:
Xiaoli Zhang, Chinese Academy of Sciences (CAS), ChinaFenfang Wang, Jimei University, China
Copyright © 2024 Raut, Barr, Paris, Kapili, Dekas and Capone. 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: Yubin Raut, yubinrau@usc.edu