Full Text

Turn on search term navigation

Introduction

Biological dinitrogen (N₂) fixation is an important source of new N in N-limited ocean gyres [1]. N₂ fixation, the energetically expensive process that converts N₂ into biologically available ammonia, is carried out by diverse Bacteria and Archaea called diazotrophs. N₂ fixation accounts for up to 70% of new N in the oligotrophic gyres [2], making diazotrophs critical components of open ocean biogeochemical cycles. In contrast, the magnitude and importance of N₂ fixation in temperate coastal environments is less well understood, having been historically ignored primarily because high concentrations of dissolved inorganic N can inhibit N₂ fixation [3]. However, there have been numerous recent reports of diazotrophs and N₂ fixation in various temperate [4–14] and polar [15–17] coastal environments. In one of the most well-studied temperate coastal regions, the Western North Atlantic continental shelf, N₂ fixation can support up to 50% of net community production [6, 7].

Marine diazotrophic cyanobacteria have diverse morphologies. The unicellular cyanobacterial group A (UCYN-A), which are obligate symbionts of single-celled haptophyte alga related to Braarudosphaera bigelowii [18], are emerging as important marine diazotrophs. The basis of the symbiosis is the transfer of photosynthetically-fixed C from the host in exchange for fixed N [18], but other metabolic interdependencies may also be essential [19], since UCYN-A lacks several important metabolic pathways [20, 21]. As such, it differs from the well-studied free-living marine diazotrophs; the biology of the host cell must play a significant role in their biogeography and activity. The impact of the UCYN-A symbiosis on the oceanic N budget is potentially significant due to its broad distribution [22–26], high cell-specific N₂ fixation and growth rates [27–29], and potential for transfer of fixed N into the food web through grazers [30–32]. UCYN-A symbioses have been reported in regions not typically assumed important for N₂ fixation, including temperate waters [4, 9, 10, 33], polar seas [16, 17, 26, 34], upwelling and neritic regions [35–37], and the California Current (CC) in the presence of dissolved inorganic N [8, 29]. Despite documentation of UCYN-A symbioses in coastal regions, their N₂ fixation activity and contribution to the biogeochemistry of these waters is unknown.

This study evaluated N₂ fixation along nearshore to offshore transects in the Southern California Current System (SCCS), which has been largely ignored as a region with respect to N₂ fixation, despite net primary production being generally N-limited at the regional scale [38]. Building upon observations of UCYN-A symbioses in the central California Current System [8, 37], this study sought to estimate their contribution to bulk community N₂ fixation along the continental shelf of the Baja California Sur Peninsula by identifying and quantifying the abundance of different diazotrophs, measuring bulk N₂ fixation rates (NFRs) and UCYN-A symbiosis single-cell NFRs.

Materials and methods

Study location and sample collection

Two cruises were conducted aboard the R/V Robert Gordon Sproul from May 3–10, 2017 (SP1714) and October 4–11, 2017 (SP1727) in the waters between San Diego (32^o 50.68' N, 117^o 31.85' W) and Cedros Island (28^o 17.34' N, 118^o 12.66' W) in Sebastián Vizcaíno Bay, Baja California Sur, Mexico. Each cruise followed a similar track, with three transects extending from nearshore to offshore waters (Fig. 1). Transect 1 (T1) was along California Cooperative Oceanic Fisheries Investigations (CalCOFI) program line 93, while Transects 2 (T2) and 3 (T3) were along lines 107 and 117, respectively, from the Investigaciones Mexicanas de la Corriente de California (IMECOCAL) program. The October cruise sampled an additional transect north of T1 (T4). Seawater was collected using a rosette of Niskin® bottles equipped with a conductivity, temperature, depth (CTD) package including fluorescence and transmissivity sensors (Seabird, Bellevue, WA), and a photosynthetically active radiation sensor (PAR; Biospherical Instruments, San Diego, CA). At each station, samples from ca. 8 depths between 0 and 200 m were collected for DNA extraction, and dissolved nutrient and chlorophyll a (Chl a) concentrations. Water for rate incubations were collected from 2 m, 10 m, and the deep chlorophyll maximum.

Fig. 1 [Images not available. See PDF.]

Map of the study region.

Schematic of the major currents and regions in the California Current System (A). Combined station maps for May 2017 (SP1714) and October 2017 (SP1727) cruises (B). Stations sampled during both cruises are solid circles, stations sampled only during SP1714 are gray circles, and stations only sampled in SP1727 are white circles.

Dissolved nutrient and chlorophyll a concentration

Samples for the measurement of nitrate plus nitrite (NO₃⁻+NO₂⁻) and phosphate (PO₄³⁻) concentrations were filtered through precombusted (450 °C for 4 h) 25 mm Whatman^TM GF/F filters (MilliporeSigma, Burlington, MA) and stored in acid-cleaned Falcon^TM tubes (Thermo Fisher Scientific, Waltham, MA) at −20 °C until analysis using standard techniques [39] on a Lachat QuikChem 8000 Flow Injection Analyzer. The limit of detection (LOD) and limit of quantitation (LOQ), respectively, were 0.01 and 0.04 µmol l⁻¹ for NO₃⁻+NO₂⁻, 0.01 and 0.02 µmol l⁻¹ for PO₄³⁻. Three separate Chl a samples (ca. 150 ml) from each depth were filtered onto GF/F filters, extracted in the dark at 3 °C in 90% acetone for 24 h and measured fluorometrically using a Turner Fluorometer TD-700 (Turner Designs, Inc., San Jose, CA) as described in Welschmeyer et al. [40]. The LOD, calculated as 3× the standard deviation of the blank (90% acetone), was 0.002 µg Chl a l⁻¹. With the exception of 200 m samples at five stations, the coefficient of variation for the averaged values was <10%.

DNA collection and extraction

Seawater was sampled into acid-cleaned polycarbonate 2 l bottles and filtered through Sterivex^TM (MilliporeSigma) filters using gentle peristaltic pumping. Sterivex^TM filters were capped, flash-frozen in liquid N₂, and stored at −80 °C. DNA was extracted using the DNeasy Plant Kit (Qiagen, Germantown, MD) using modifications to the manufacturer’s guidelines described in detail in Moisander et al. [41]. On-column steps were automated using a QIACube (Qiagen). DNA was quantified using the Picogreen® dsDNA Quantitation kit (Molecular Probes, Eugene, OR).

Diazotroph community composition and UCYN-A nifH oligotyping

Diazotroph community composition was characterized using nifH amplicon sequencing using a custom pipeline described in detail by Cabello et al. [37]. Operational taxonomic units (OTUs; herein identified as “denovo”) were defined at 97% nucleotide identity in QIIME [42] using Usearch6.1 [43, 44]. A total of 4,288,726 nifH sequences remained after removing low quality sequences and OTUs with <100 sequences, ranging from 1 to 54,726 per sample. The resulting OTU table was subsampled to 2312 sequences, which removed 33 samples primarily from SP1714. Taxonomy was assigned via BLASTX using full length nifH amino acid sequences (publicly available at jzehrlab.com/nifh) with nifH cluster designations based on Zehr et al. [45]. Raw sequences are available in the NCBI Sequence Read Archive (BioProject PRJNA695866).

UCYN-A nifH sequence types were resolved using oligotyping [46], using the exact entropy positions and oligotyping arguments described by Turk-Kubo et al. [47]. To reevaluate global patterns of UCYN-A oligotype distributions, SCCS data were combined with global survey data [47] and other recently published datasets from the Noumea Lagoon in New Caledonia [48], Bering and Chukchi Seas [17], North Pacific Subtropical Gyre (NPSG) [49], and Monterey Bay, CA [37]. This compiled dataset was rarefied to 1029 sequences prior to ordination analyses using Phyloseq [50] in R (r-project.org). The dissimilarity between samples was determined using Bray–Curtis ecological indices and Principal Coordinate Analysis (PCoA) on the resulting distance matrices to visualize dissimilarity between samples and co-occurring oligotypes.

Targeted diazotroph abundance estimates via qPCR

Gene-based abundance estimates of UCYN-A1 [51], UCYN-A2 [52], Crocosphaera (UCYN-B; [24]), Trichodesmium [51], Richelia associated with Hemiaulus (Het-2; [53]), and gamma A (γ-24774A11; [41]) were determined using Taqman® qPCR assays. Protocols used for all aspects of qPCR analysis, including reaction conditions, the use of linearized plasmids and inhibition reactions, and calculation of unknowns follow those described in detail by Goebel et al. [54], apart from a 64 °C annealing temperature for the UCYN-A2 assay. The LOD and LOQ for all assays ranged between 25-31 and 200-250 nifH copies l⁻¹, respectively. Targets with nifH copies >LOD and <LOQ are detected not quantified (DNQ).

The UCYN-A1 qPCR assay has high specificity, but the UCYN-A2 assay cross-reacts with UCYN-A3 and UCYN-A4 [55]. SCCS UCYN-A assemblages were dominated by UCYN-A1 and UCYN-A2; UCYN-A4 was present at low relative abundance, while UCYN-A3 was rarely detected. Therefore, nifH-based abundance using the UCYN-A2 assay may also include UCYN-A4.

Bulk community N₂ fixation rate measurements

Seawater was sampled directly from Niskin® bottles into acid-washed 1.2 l polycarbonate bottles through 210 µm Nitex® mesh (Wildco, Yulee, FL) to remove large grazers. Incubation bottles received 100 ml of ¹⁵N₂-enriched seawater. ¹⁵N₂-enriched seawater [56] was generated and atom% enrichment was measured according to procedures described in detail by Mills et al. [29]. The ¹⁵N₂-enriched seawater atom% enrichment ranged from 2.0–6.1% for SP1714 and 5.1–24.7% for SP1727. Bottles were incubated (24 h) under simulated in situ light using neutral density screening and maintained at surface seawater temperatures in flow-through on-deck incubators. Samples for atom% ¹⁵N of the ambient particulate matter were taken from corresponding depths at T₀. At the termination of the incubation, samples for the analysis of ¹⁵N enrichment into particulate organic matter (ca. 1000 ml) were processed and measured, and NFRs were calculated, as detailed in Mills et al. [29]. LOD and minimum quantifiable rates (MQRs) were calculated as in Montoya et al. [57] and Gradoville et al. [58] in accordance with recommendations by White et al. [59]. LODs ranged from 0.2–3.0 nmol N l⁻¹ d⁻¹ (May) and 0.3–5.4 nmol N l⁻¹ d⁻¹ (October).

UCYN-A symbioses single-cell N₂ fixation rates

Subsamples (95 ml) taken from the incubation bottles were fixed with sterile-filtered formaldehyde (MilliporeSigma) at a final concentration of 1.85% (v/v) for >1 h at 4 °C, then concentrated with 0.6 µm pore-size polycarbonate filters (MilliporeSigma) under gentle vacuum, air-dried and stored at −80 °C. UCYN-A1 and UCYN-A2 symbioses were targeted using 5′-horseradish peroxidase-labeled oligonucleotide probes (Biomers, Inc., Ulm/Donau, Germany), using helper and competitor probes for both symbionts and hosts (Biomers) as described in detail by Cornejo-Castillo et al. [60]. Protocols for CARD-FISH hybridizations followed procedures described in detail by Cabello et al. [61].

Samples were visualized, transferred, and mapped to facilitate nanoSIMS analyses according to protocols detailed in Mills et al. [29]. Individual symbioses were analyzed on a Cameca nanoSIMS 50 L at the Stanford Nano Shared Facilities (Stanford, CA). Once targets were located using the charged-coupled device camera and the secondary electron image, image fields were rastered with a 16 keV Cesium primary ion beam (~5 pA) focused into ca. 120 nm spot diameter (256 × 256 pixels, dwell time 1 ms per pixel). Images of ¹²C⁻, ¹³C⁻, ¹²C¹⁴N⁻ and ¹²C¹⁵N⁻ were measured over 30–100 planes with a mass resolving power of ca. 8000. Regions of interest were defined around UCYN-A and host cells using Look@nanoSIMS [62]. Isotope ratios of UCYN-A and haptophyte cells were calculated as described in Mills et al. [29]. Single-cell NFRs were determined as follows $N_{2} f i x a t i o n r a t e (f m o l {c e l l}^{- 1} d^{- 1}) = \frac{A_{P N}^{f i n a l} - A_{P N}^{t = 0}}{(A_{N_{2}} - A_{P N}^{t = 0})} * \frac{[{P N}_{c e l l}]}{Δ t}$ where ρ equals the absolute uptake rate per cell, $A_{P N}^{f i n a l}$ and $A_{P N}^{t = 0}$ equal the atom% ¹⁵N of the enriched (final) or unenriched (t = 0) UCYN-A symbiosis, $A_{N_{2}}$ is the atom% enrichment of the N₂ source pool, Δt is the incubation time, and $[{P N}_{c e l l}]$ is the per cell N quota. Cell biomass estimates utilized biovolumes according to Krupke et al. [63] and were converted to per cell N quotas using C:N estimates from Martinez-Perez et al. [28]. As with bulk NFR, the LOD and MQR were calculated as in Montoya et al. [57] and Gradoville et al. [58]. NFRs for associated UCYN-A and haptophyte cells were calculated individually and then summed to obtain a single-cell NFR for the entire symbiosis.

Contributions of UCYN-A symbioses to bulk NFR were estimated using single-cell NFRs and nifH-based abundance assuming 1 nifH copy cell⁻¹ for the UCYN-A1 symbiosis and 10 nifH copy cell⁻¹ for the UCYN-A2 symbiosis, which is at the high end of the estimated nifH copies cell⁻¹ range previously reported [52]. At present, this is the best estimate that can be justified, but may result in an underestimation of UCYN-A2 symbiosis contribution. Maximum single-cell NFRs were determined using a dilution factor of 75% to correct for isotope dilution during sample processing [64].

Statistical analysis

Normality was assessed using the Shapiro–Wilk test. The association between diazotroph abundance, NFR, and measured environmental variables was evaluated using the non-parametric Spearman’s rank correlation (ρ), to account for the non-normal distribution of some variables. General linear models were also developed (see Supplemental methods). Both approaches agreed overall; Spearman’s ρ results are emphasized in the Results and Discussion sections. Comparison of seasonal means for environmental parameters was evaluated using the Mann–Whitney U test or t-test for non-normally and normally distributed data, respectively.

Results

Seasonal differences in oceanographic conditions

Surface water (0–50 m) temperatures were lower in May (14.4 ± 2.2 °C) than in October (17.5 ± 3.2 °C) (Fig. 2A, F, Tables S1 and S2) while salinity differences were small (0.1 ± 0.2) between expeditions (Fig. 2B, G). The influence of tropical surface water in the study region was evident in October where surface waters had a potential density anomaly (γ_θ) < 24 kg m⁻³, a signature of poleward transport of tropical surface water [(65]; Fig. S1). During May, the depth of the mixed layer increased with distance from shore (Tables S1 and S2). Along T1 and T2, mixed layer depths averaged 15 ± 17.5 and 18 ± 12.2 m, respectively, with the deepest mixed layers at stations furthest from the coast. Along T3, mixed layer depths were shallower (10 ± 6.7 m), but also deepened offshore. In October, mixed layer depths were shallower and less variable, averaging 11 ± 1.9 m, 10 ± 1.3, and 8 ± 3.9 m along T1, T2, T3, respectively.

Fig. 2 [Images not available. See PDF.]

Physical, chemical and biological conditions in the SCCS.

Section plots of temperature (A, F), salinity (B, G), NO₃⁻ (C, H), PO₄³⁻ (D, I), and Chl a (E, J) along T1–T4 from May 2017 (SP1714 cruise; A–E) and October 2017 (SP1727 cruise; F–J). Transects and stations along each transect are noted at top.

The Bakun Upwelling Index [66] indicated upwelling just prior to the May cruise along all transects (Fig. S2). Upwelling decreased by the start of the cruise and remained relatively weak during T2 and T3. A large upwelling event took place on May 7 along T1, which preceded the re-sampling at Stn. 1 (Stn 1.1) and Stns. 2–5. In October there were no distinct upwelling events prior to or during sampling.

NO₃⁻+NO₂⁻ concentrations in waters ≤10 m were significantly lower (Mann–Whitney U test, n_May = 24, n_Oct = 28, U = 90, p < 0.01) during October (<0.03 µmol l⁻¹) than in May (0.57 ± 1.10 µmol l⁻¹) (Fig. 2C, H). The pattern was similar at depths between 11 and 50 m, with NO₃⁻+NO₂⁻ concentrations averaging 8.01 ± 7.10 µmol l⁻¹ in May and 2.43 ± 3.65 µmol l⁻¹ in October (Mann–Whitney U test, n_May = 30, n_Oct = 34, U = 257, p < 0.01). Below 50 m, where waters were influenced by Equatorial Subsurface water (Fig. S1), NO₃⁻+NO₂⁻ concentrations in May (20.22 ± 8.42 µmol l⁻¹) and October (21.30 ± 8.29 µmol l⁻¹) were similar (Mann–Whitney U test, n_May = 35, n_Oct = 47, U = 717, p > 0.05). Like NO₃⁻, average surface water (≤10 m) PO₄³⁻ concentrations were also higher in May (0.29 ± 0.79 µmol l⁻¹) than in October (0.05 ± 0.02 µmol l⁻¹) (Mann–Whitney U test, n_May = 24, n_Oct = 26, U = 81, p < 0.01, Fig. 2D, I) and no seasonal difference was detected below 50 m (May—1.28 ± 0.58 µmol l⁻¹, October—1.44 ± 0.53 µmol l⁻¹, Mann–Whitney U test, n_May = 35, n_Oct = 47, U = 647, p > 0.05). Both nitraclines and phosphoclines generally deepened with distance from the coast in May and along T3 in October, but varied little with depth across T1, T2, and T4 in October.

P* (P* = PO₄³⁻ − (NO₃⁻ + NO₂⁻)/16), the amount of dissolved PO₄³⁻ in the environment relative to what is expected if N and P uptake and remineralization proceed according to Redfield proportions [67], was on average slightly positive across the study region (Tables S1 and S2), suggestive of conditions conducive to N₂ fixation, and possibly reflecting a signature of ODZ-derived waters that are transported into the Southern California Bight by the California Undercurrent [68]. In May, P* ranged from −0.81 to 3.91 µmol l⁻¹, while in October the range was smaller (−0.20 to 0.52 µmol l⁻¹). P* varied little with depth or distance from the coast.

At the surface (<10 m), May Chl a concentrations were higher than in October, 0.8 ± 1.17 vs. 0.1 ± 0.04 mg Chl a l⁻¹ (Fig. 2E, J). Chl a concentrations were also higher in May at depths between 11 and 50 m (0.6 ± 1.22 vs. 0.3 ± 0.16 mg Chl a l⁻¹), but deep Chl a maxima were more prominent in October, evidenced by the peak in Chl a concentration between 11–50 m (Figs. 2J and S1). Chl a concentrations generally decreased with distance from coast.

Diazotroph community composition

Partial nifH genes were amplified from all samples (76 from May and 105 from October); 148 remained after resampling (49 from May and 99 from October). The rarefied dataset is represented by 2178 OTUs. The majority of sequences affiliated with UCYN-A sublineages (82.6%), cluster 1G (putative γ-proteobacteria; 9.4%), and cluster III (putative δ- proteobacteria and other anaerobes; 4.8%) (Fig. 3). Relative abundances of the cyanobacterial diazotrophs Trichodesmium, Crocosphaera, Cyanothece, and Richelia associated with Rhizosolenia (Het-1) and Hemiaulus (Het-2) were low (0.2%).

Fig. 3 [Images not available. See PDF.]

SCCS diazotroph community composition.

Relative abundance of diazotrophs based on the sequencing of partial nifH genes during May 2017 (A) and October 2017 (B) cruises. Columns represent diazotroph taxonomic groups, and the size of the bubble scales to relative abundances for each group at each station/depth (rows). NCD taxonomic groups are identified using shaded columns, and gray lines delineate stations. Cluster designations are as described in Zehr et al. [45].

Non-cyanobacterial diazotrophs (NCDs) accounted for 17.2% of the sequences, had high relative abundances in May, and were recovered from surface waters and below the deep chlorophyll max (Fig. 3A). Cluster 1G sequences were predominantly affiliated with gamma A [69] which had high relative abundances in October along T1 and T4 (Fig. 3B). Cluster III sequences were found in deeper samples (>50 m) and were dominated by several OTUs (Fig. 3). The most highly recovered cluster III OTU, denovo5, had 99% nucleotide identity to a South Pacific sequence type, cIII-ETSP [70]. The second most abundant cluster III OTU, denovo9, was not similar to any previously reported nifH sequence. Notably, two additional NCD OTUs have previously been reported: denovo13 which is similar to a Klebsiella-like OTU (OTU0009; [11]) and denovo30, which is identical to a Western Pacific sequence type, Alpha-MH144511 [71].

UCYN-A nifH oligotyping

UCYN-A sublineages are comprised of multiple oligotypes but global sequence libraries are dominated by oligo1 (UCYN-A1), oligo2 (UCYN-A3), and oligo3 (UCYN-A2) [47]. SCCS UCYN-A assemblages contained 89 oligotypes; 34 were previously unreported but together accounted for less than 1% of UCYN-A sequences. The five oligotypes with the highest relative abundance were oligo1 (71.5%), oligo3 (20.7%), oligo4 (affiliated with UCYN-A4; 2.5%), oligo46 (UCYN-A2; 1.0%), and oligo45 (UCYN-A1; 0.5%) (Figs. S3 and S4). UCYN-A4 oligotypes (dominated by oligo4) had high relative abundances in some nearshore stations at depths >25 m (Figs. 3, S3, and S4), particularly May Stns. 1.1 and 7 and along the two northernmost transects in October (Figs. S3 and S4).

PCoA using the Bray–Curtis dissimilarity index indicates that oligo1 co-occurred with minor UCYN-A1 oligotypes oligo8, oligo9, and oligo11, and most UCYN-A2 oligotypes co-occurred and cluster separately from UCYN-A1 (Fig. 4). Oligo4 co-occurred with UCYN-A2 along the PCoA axis 1. Notably, oligo46 (UCYN-A2) and oligo45 (UCYN-A1) formed a separate cluster (along with oligo13 and oligo40) that deviated from the predominant co-occurrence of oligotypes within a given sublineage. Oligo46 and oligo45 were found at Stn. 2 in May samples and throughout the region in October (Figs. S3 and S4).

Fig. 4 [Images not available. See PDF.]

Patterns of UCYN-A nifH oligotype biogeography. Dissimilarity between samples is based on Principal Coordinate Analysis using the Bray–Curtis ecological index. Co-occurrence of SCCS UCYN-A oligotypes in the SCCS (A). Oligotypes with the highest relative abundances are labeled. Patterns in UCYN-A biogeography using samples collected in nearshore and offshore environments throughout the global oceans (B). Dissimilarity patterns for UCYN-A nifH oligotypes (B) and samples (C) indicate that SCCS samples fall between nearshore (circle) and offshore (triangle) endpoints, which are characterized by UCYN-A2 and UCYN-A1 symbioses, respectively. SCCS data were combined with other recent datasets [17, 37, 47–49].

Abundance and distribution of targeted diazotrophs

UCYN-A symbioses were the most abundant diazotrophs and detected in nearly all surface samples during both cruises. Average surface abundances of the UCYN-A1 symbiosis in October were significantly (Mann–Whitney U test, n_May = 8, n_Oct = 14, U = 20, p < 0.01) higher than in May, at 1.6 × 10⁶ ± 1.4 × 10⁶ and 3.6 × 10⁵ ± 4.7 × 10⁵nifH l⁻¹, respectively (Fig. 5, Tables S1 and S2), and in October, maximum abundances were generally associated with waters carrying a signature of tropical surface water (Fig. S1). Surface abundances of the UCYN-A2 symbiosis were lower and did not significantly differ between cruises, with May averages of 1.8 × 10⁵ ± 2.5 × 10⁵nifH l⁻¹ and October averages of 1.2 × 10⁵ ± 1.7 × 10⁴nifH l⁻¹. UCYN-A1 and UCYN-A2 symbiosis abundances were positively and significantly associated across both cruises (ρ = 0.71, p < 0.001). In addition, both had significant association with temperature and oxygen (ρ > 1) and depth, NO₃⁻+NO₂⁻, PO₄^3- and σθ (ρ < 1; Table S3).

Fig. 5 [Images not available. See PDF.]

Patterns in qPCR-based abundances of dominant diazotroph taxa.

Distributions of UCYN-A1 symbiosis (A, D), UCYN-A2 symbiosis (B, E), and Het-2 (C, F) in May 2017 (A–C) and October 2017 (D–F). Stations and transects are indicated at the top. Stn. 14 (October, T4) is not included due to missing surface data. Diazotrophs present at low abundances or sporadically detected are not included (Tables S1 and S2).

May UCYN-A symbioses distribution patterns were heterogenous. Peak UCYN-A1 symbiosis abundance along T1 was in warm, surface coastal waters above a shoaling nitracline, while found offshore along T2 and T3 (Figs. 2, 5A, S5a). UCYN-A2 symbiosis distribution patterns were similar, with the exception of T2, where high volumetric and depth-integrated abundances were measured at nearshore Stn. 7 (Figs. 5B and S5C).

In October, both symbioses were distributed throughout surface waters along all four transects and were detected at every station (Figs. 5D, E and S5B, D). Notably, the UCYN-A1 symbiosis was detected at high abundances in nearshore stations along T1–T3 (Fig. 5D). Considering only October data, UCYN-A1 symbiosis abundances were significantly and inversely (ρ < 1) associated with depth, salinity, NO₃⁻+NO₂⁻, PO₄³⁻ and σθ, and positively (ρ > 1) associated with temperature, oxygen and Chl a. In addition, total UCYN-A symbioses abundance was significantly and inversely associated with P* during October (ρ = −0.21, p = 0.03) but not May (ρ = 0.14, p = 0.21; Table S3).

The second most abundant N₂-fixer was Het-2, with peak abundances of 2.1 × 10³ ± 1.3 × 10²nifH l⁻¹ (Stn. 9, 25 m) and 2.2 × 10⁴ ± 6.8 × 10²nifH l⁻¹ (Stn. 16, 20 m) in May and October, respectively (Fig. 5C, F, Tables S1 and S2). In May, Het-2 was undetected or DNQ in most samples, except for Stn. 9. Depth-integrated abundances were highest in October, notably along the northern two transects (Fig. S5E, F). Het-2 abundance was significantly and positively associated with UCYN-A1 (ρ = 0.60, p < 0.001) and UCYN-A2 (ρ = 0.54, p < 0.001) symbioses abundances and temperature, oxygen and PAR, and significantly and inversely associated with depth, salinity, NO₃⁻+NO₂⁻, PO₄³⁻, σθ (Table S3).

Trichodesmium, Crocosphaera (UCYN-B), and gamma A were sporadically detected (Tables S1 and S2). Trichodesmium, present at low abundance in October, had peak abundances at Stn. 14 (3.6 × 10⁴ ± 4.2 × 10¹nifH l⁻¹). Crocosphaera was DNQ at two stations in May (Stns. 4, 10), and at three stations in October (Stns. 2, 7, 13) and only quantifiable in one October sample (Stn. 14, 50 m, 5.1 × 10² ± 5.1 × 10⁰nifH l⁻¹). Gamma A was abundant in Stn. 1 surface waters in May (1.4 × 10⁵ ± 3.4 × 10⁴nifH l⁻¹) and at Stn. 9 (1.5 × 10⁴ ± 7.5 × 10³nifH l⁻¹) but was undetected in all October surface samples. Gamma A sequences were recovered from virtually all October samples (Fig. 3B), but qPCR data indicated they were below quantitation (<25 nifH l⁻¹), consistent with known over-representation in sequence libraries [72].

Bulk community N₂ fixation rates

In May, N₂ fixation was detected at the surface and deep chlorophyll maximum at every station (except Stn. 8 where NFR was not measured), but October NFRs were patchy. May NFRs averaged 8.5 ± 6.5 nmol N l⁻¹d⁻¹ in surface waters and 6.7 ± 3.3 nmol N l⁻¹ d⁻¹ at the deep chlorophyll maximum (Tables 1 and S4). The highest May volumetric rates were measured in nearshore surface waters at Stn 1.1 (14.2 ± 4.5 nmol N l⁻¹ d⁻¹) and Stn. 13 (23.0 ± 3.8 nmol N l⁻¹ d⁻¹). Concurrent with deeper MLDs and higher upwelling indices, N₂ fixation was often detected at the deep chlorophyll maximum in the presence of NO₃⁻+NO₂⁻ concentrations >1 µM. May NFRs were not significantly associated with environmental parameters or diazotroph abundances (Table S3).

Table 1. Compilation of environmental parameters and volumetric euphotic NFR during May 2017 (SP1714) and October 2017 (SP1727).

Cruise	Station	Depth (m)	Temp. (°C)	Salinity	Nitrate + Nitrite (µM)	Phosphate (µM)	P*	Chl a (mg l⁻¹)	NFR (nmol N l⁻¹ d⁻¹)	NFR LOD (nmol N l⁻¹ d⁻¹)
SP1714	1	2	18.22	33.40	0.05	0.09	0.08	0.173 ± 0.003	2.3 ± 1.3	0.4
	1	30^a	11.91	33.41	9.82	0.88	0.27	0.900 ± 0.048	3.1 ± 0.7	0.4
	1.1	2	16.82	33.41	0.01	0.01	0.01	0.213 ± 0.004	14.2 ± 4.5	0.7
	1.1	27^a	12.33	33.36	nm	nm	na	3.012 ± 0.046	11.5 ± 13.1	0.5
	2	2	16.32	33.42	0.10	0.07	0.07	0.508 ± 0.003	6.7 ± 0.7	0.3
	2	29^a	13.63	33.27	9.05	0.64	0.07	0.950 ± 0.033	3.2 ± 2.0	0.6
	3	2	nm	nm	1.65	0.13	0.03	0.943 ± 0.009	4.8 ± 1.9	0.6
	3	25^a	14.37	33.47	nm	nm	na	0.877 ± 0.013	4.0 ± 0.9	0.7
	4	2	nm	nm	nm	nm	na	nm	5.2 ± 1.4	0.9
	4	25^a	14.18	33.32	1.31	0.19	0.11	0.509 ± 0.010	8.2 ± 6.4	1.0
	5	2	14.54	33.23	0.12	0.10	0.09	0.307 ± 0.013	7.1 ± 3.6	0.5
	5	50	14.51	33.23	0.28	0.06	0.05	0.283 ± 0.013	6.4 ± 2.9	0.5
	7	10	14.16	33.31	4.58	0.45	0.16	0.257 ± 0.006	6.1 ± 5.3	1.1
	7	32^a	12.15	33.37	12.73	0.75	−0.04	0.275 ± 0.012	9.5 ± 1.6	0.8
	9	10	17.70	33.43	−0.01	0.05	0.05	0.164 ± 0.016	10.0 ± 2.6	0.5
	9	52^a	13.53	33.26	4.38	0.32	0.05	0.398 ± 0.022	9.5 ± 2.6	0.3
	10	90	14.76	33.31	1.22	0.11	0.04	0.241 ± 0.006	3.2 ± 1.4	0.2
	12	2	16.10	33.51	0.02	0.06	0.06	0.765 ± 0.127	3.6 ± 2.0	0.4
	12	20^a	13.47	33.44	4.97	0.59	0.28	1.884 ± 0.207	3.3 ± 0.9	0.6
	13	2	14.99	33.59	0.25	0.31	0.29	3.994 ± 0.055	23.0 ± 3.8	2.3
	13	10	13.60	33.49	2.84	0.33	0.15	4.268 ± 0.173	DNQ	2.5
	13	25^a	13.54	33.58	4.60	0.30	0.02	6.044 ± 0.916	11.1 ± 7.4	3.0
	14	2	18.03	33.53	0.04	0.14	0.14	0.178 ± 0.009	9.9 ± 4.9	0.8
	14	10	18.00	33.52	0.02	0.16	0.16	0.174 ± 0.013	9.0 ± 4.1	0.7
	14	30^a	14.09	33.18	4.55	0.43	0.15	1.133 ± 0.084	8.0 ± 7.9	1.4
SP1727	1	2	21.76	33.56	0.04	0.02	0.02	0.120 ± 0.007	2.2 ± 0.2	0.7
	1	10	21.35	33.51	0.00	0.02	0.02	0.149 ± 0.010	5.1 ± 0.2	1.8
	1	38^a	14.41	33.32	2.22	0.27	0.13	0.595 ± 0.058	BDL	2.8
	2	2	20.46	33.58	0.00	0.03	0.03	0.115 ± 0.005	3.9 ± 0.2	1.4
	2	10	20.47	33.58	0.00	0.03	0.03	0.110 ± 0.001	5.8 ± 0.8	2.0
	2	35^a	13.71	33.25	1.03	0.19	0.13	0.495 ± 0.017	BDL	3.3
	3	3	18.91	33.51	0.00	0.02	0.02	0.188 ± 0.013	5.7 ± 0.2	1.2
	3	10	18.88	33.49	0.00	0.07	0.07	0.201 ± 0.006	6.3 ± 0.4	1.5
	3	35^a	14.62	33.16	0.21	0.11	0.10	0.390 ± 0.010	BDL	1.0
	4	3	18.58	33.44	0.00	0.07	0.07	0.100 ± 0.001	BDL	0.6
	4	10	18.58	33.44	0.00	0.08	0.08	0.095 ± 0.003	BDL	0.6
	4	50^a	14.07	33.20	0.65	0.16	0.12	0.249 ± 0.002	BDL	0.5
	5	3	18.90	33.43	0.00	0.09	0.09	0.099 ± 0.001	BDL	1.0
	5	10	18.83	33.43	0.00	0.07	0.07	0.105 ± 0.005	BDL	0.7
	5	55^a	14.26	33.34	0.24	0.13	0.12	0.309 ± 0.012	BDL	0.5
	7	3	19.95	33.39	0.00	0.07	0.07	0.161 ± 0.007	6.7 ± 0.8	1.2
	7	10	19.54	33.23	0.00	0.07	0.07	0.196 ± 0.010	16.5 ± 1.8	3.3
	7	27^a	13.71	33.23	0.26	0.22	0.20	0.689 ± 0.044	BDL	5.4
	9	3	20.77	33.51	0.00	0.03	0.03	0.092 ± 0.012	DNQ	0.4
	9	10	20.75	33.50	0.00	0.03	0.03	0.094 ± 0.003	5.8 ± 1.2	0.9
	9	55^a	13.85	33.27	1.46	0.23	0.14	0.271 ± 0.014	DNQ	0.8
	11	3	21.04	33.56	0.00	nm	na	0.088 ± 0.000	5.9 ± 0.9	0.6
	11	10	21.02	33.57	0.00	0.04	0.04	0.091 ± 0.003	5.8 ± 0.38	0.5
	11	67^a	14.11	33.29	0.06	0.18	0.18	0.431 ± 0.015	BDL	0.5
	12	2	20.18	33.44	0.00	0.05	0.05	0.142 ± 0.003	14.1 ± 10.0	0.3
	12	10	19.76	33.34	0.00	0.05	0.05	0.154 ± 0.010	7.0 ± 0.2	0.6
	12	38^a	14.80	33.29	1.73	0.33	0.22	0.678 ± 0.024	BDL	0.9
	13	3	21.14	33.50	0.00	0.04	0.04	0.098 ± 0.003	9.3 ± 1.9	0.5
	13	10	21.14	33.50	0.00	0.05	0.05	0.093 ± 0.001	19.6 ± 1.5	1.1
	13	68^a	13.96	33.30	0.78	0.25	0.20	0.269 ± 0.011	1.4 ± 1.3	1.2
	14	3	21.74	33.57	0.00	0.05	0.05	0.063 ± 0.004	nm	nm
	14	10	21.53	33.56	0.00	0.05	0.05	0.078 ± 0.002	3.4 ± 0.2	1.1
	14	50	18.28	33.33	0.00	0.08	0.08	0.105 ± 0.005	BDL	1.6
	15	3	19.31	33.43	nm	nm	na	0.021 ± 0.006	7.3 ± 0.6	0.6
	15	10	19.26	33.41	nm	nm	na	0.211 ± 0.008	DNQ	0.8
	15	38^a	14.05	33.30	nm	nm	na	0.868 ± 0.039	DNQ	0.6
	17	3	17.67	33.42	0.00	0.05	0.05	0.314 ± 0.129	BDL	1.5
	17	10	17.47	33.42	0.00	0.05	0.05	0.388 ± 0.037	BDL	1.9
	17	27^a	15.95	33.28	1.40	0.18	0.09	1.494 ± 0.402	BDL	2.9

Detection limits and minimum quantifiable rates are detailed in Table S4.

nm not measured, na not applicable; DNQ detected not quantified, BDL below detection limit.

^aDepth in the vicinity of the deep chlorophyll maximum.

NFRs were lower in October, with surface rates averaging 4.3 ± 4.4 nmol N l⁻¹ d⁻¹ and near or below detection limits at the deep chlorophyll maximum. NFRs in October were highest at 10 m at Stns. 13 (19.6 ± 1.5 nmol N l⁻¹ d⁻¹) and 7 (16.5 ± 1.8 nmol N l⁻¹ d⁻¹). October NFRs were significantly and positively associated with temperature, salinity, UCYN-A1 and UCYN-A2 symbioses abundance, and inversely associated with depth, fluorescence, oxygen, σθ, NO₃⁻+NO₂⁻, PO₄³⁻, P* and Chl a (Table S3).

Depth-integrated NFR ranged from 62.0–409.3 µmol N m⁻² d⁻¹ in May (Fig. S5G), and from 0–709.1 µmol N m⁻² d⁻¹ in October (Fig. S5H). In May, high depth-integrated NFRs were measured in both nearshore and offshore stations, while in October, depth-integrated NFR was highest at Stn. 13, where the bottom depth was less than 200 m (Table S2, Fig. S5H).

UCYN-A single-cell N₂ fixation rates

Single-cell NFRs for the UCYN-A1 symbiosis were measured at Stn. 1 (May), Stn. 5 (October) and Stn. 14 (October), and ranged from BDL—30.5 fmol N cell⁻¹ d⁻¹, with an average of 6.6 ± 8.8 fmol N cell⁻¹ d⁻¹ (Fig. 6, Tables 2 and S5). Per cell N quotas ranged from 2.6–27.7 fmol N cell⁻¹ (Table S5). UCYN-A1 symbiosis single-cell NFRs had high variability, most notably at Stn. 5 (October), where the average rate was 13.1 ± 14.6 fmol N cell⁻¹ d⁻¹, and below detection in 3 of 6 associations (Table S5).

Fig. 6 [Images not available. See PDF.]

Compilation of UCYN-A symbioses single-cell N₂ fixation rates.

UCYN-A1 symbiosis (A) and UCYN-A2 symbiosis (B) single-cell NFRs in the SCCS are similar to rates reported in offshore waters. SCCS stations are identified on the x-axis. Note the different scales on the y-axis. Data compiled from the following studies where rates from intact symbioses are reported: “Arctic” [17], “SIO” [29], “Natl” [28], and “NPSG” [49].

Table 2. Estimated contribution of UCYN-A1 and UCYN-A2 symbioses to community NFRs at select near- and offshore stations.

Cruise	Stn	UCYN-A1/UCYN-A2	Number of cells	Measured			Estimated		Corrected estimate (maximum)
Cruise	Stn	UCYN-A1/UCYN-A2	Number of cells	UCYN-A single-cell NFR (fmol N cell⁻¹ d⁻¹)	UCYN-A abundance (nifH copies l⁻¹)	Community NFR (nmol N l⁻¹ d⁻¹)	UCYN-A volumetric NFR (nmol N l⁻¹ d⁻¹)	UCYN-A contribution to community NFR	Maximum UCYN-A single-cell NFR (fmol N cell⁻¹ d⁻¹)	Maximum UCYN-A volumetric NFR (nmol N l⁻¹ d⁻¹)	Maximum UCYN-A contribution to community NFR
SP1714	1	A1	10	3.6 (2.2)	4.1 × 10⁴ (7.3 × 10²)	2.3 (1.3)	0.15 (0.0)	6% (43%)	14.4	0.6	26%
SP1727	1	A2	6	125.8 (41.6)	9.3 × 10⁴ (4.2 × 10³)	2.2 (0.2)	1.7 (0.1)^c	53% (16%)	503.2	4.7	213%
SP1727	5	A1	6	12.8 (13.5)	6.2 × 10⁴ (1.5 × 10³)	0.8 (0.1)^a	0.8 (0.01)	99% (6%)	51.2	3.2	397%
SP1727	7	A2	6	279.9 (53.7)	1.2 × 10⁵ (2.9 × 10³)	6.7 (0.8)	3.6 (0.1)^c	50% (28%)	1119.6	134	201%
SP1727	14	A1	10	5.7 (5.4)	8.8 × 10⁵ (8.1 × 10⁴)	3.4 (0.2)^b	5.0 (2.7)	148% (270%)	22.8	20.1	590%
SP1727	14	A2	5	27.0 (15.4)	2.9 × 10⁴ (1.1 × 10²)	3.4 (0.2)^b	0.8 (0.0)^c	2% (7%)	108	0.3	9%

Standard deviations for measured parameters and propagated errors for estimated parameters are in parenthesis. Maximum single-cell NFR using correction factors defined by [64], and the resulting maximum UCYN-A volumetric NFR and contribution to community are also reported.

^a NFR was BDL, value used is the detection limit.

^b NFR data from 10 m; N content at 3 m was below levels needed for accurate and/or precise isotopic measurement.

^c Potential contributions to volumetric NFR rates for the UCYN-A2 symbiosis assume that there are 10 nifH copies cell⁻¹ [64].

Single-cell NFRs for the UCYN-A2 symbiosis were measured at three stations in October (Stns. 1, 7, and 14) and ranged from 2.2–362.2 fmol N cell⁻¹ d⁻¹, with an average of 151.1 ± 112.7 fmol N cell⁻¹ d⁻¹. Per cell N quotas ranged from 49.3–159.8 fmol N cell⁻¹ (Table S5). UCYN-A2 symbiosis single-cell NFR was the highest nearshore at Stn. 7 (279.9 ± 52.1 fmol N cell⁻¹ d⁻¹) and lowest offshore at Stn. 14 (27.0 ± 16.7 fmol N cell⁻¹ d⁻¹).

Discussion

N₂ fixation in the SCCS

Marine N cycle measurements have not been focused on N₂ fixation in coastal SCCS waters; however, our study demonstrates that it is widespread along the continental shelf between the Southern California Bight and Sebastián Vizcaíno Bay, and in adjacent offshore waters. NFRs in surface waters were high throughout the study region in both May and October and volumetric rates fall into the mid-to-high range of previously reported rates in nearshore environments (Table S6). Together with rates measured in the northern portion of the Southern California Bight, the Eastern Tropical North Pacific, and the Gulf of California [8, 73, 74], our study extends the region of N₂ fixation in the SCCS. The relationships between chemical and physical properties of the water column and bulk N₂ fixation are broadly consistent with the general paradigm of marine N₂ fixation occurring in warm, sunlit, N-deplete ocean waters [75], particularly during the fall oceanic season, where NFRs were positively associated with temperature and inversely associated with depth, NO₃⁻+NO₂⁻, and Chl a (Tables S3, S7, and S8). However, whole community NFRs in this region appear insensitive to excess PO₄³⁻ (P*), which may reflect the dominance of UCYN-A and the insensitivity of N₂ fixation by the UCYN-A symbioses to nitrate concentrations [29, 63].

Despite this overall trend, some of the highest NFRs were measured in May, when upwelling conditions existed 4 days prior to the measurements (Fig. S2). The upwelling was reflected in lower sea surface temperatures and higher NO₃⁻+NO₂⁻, PO₄³⁻, and Chl a concentrations, compared to October. Notably, one station was occupied before (Stn. 1 on May 3, 2017) and after (Stn 1.1 on May 9, 2017) the upwelling event and surface NFRs significantly were higher after the upwelling, at 2.3 ± 1.3 prior to and 14.2 ± 4.5 post-upwelling (unpaired t-test, p = 0.012). At this station, UCYN-A1 and UCYN-A2 symbioses abundance increased by 1 and 2 orders of magnitude post-upwelling, respectively, and Chl a concentrations from the surface to 50 m also increased (Tables 1 and S1). Together these observations suggest a link between post-upwelling conditions and increased NFRs by UCYN-A symbioses. However, a focused field campaign designed to test this directly is needed to reproduce these results and determine whether post-upwelling conditions differentially stimulate NFRs in UCYN-A1 vs. UCYN-A2 symbioses.

Enhanced primary productivity is suspected to stimulate N₂ fixation in UCYN-A symbioses. A link between UCYN-A symbioses presence and/or activity and indicators of primary productivity (namely Chl a) has been reported [6, 24], but the underlying mechanism(s) remain unidentified. However, stimulation of NFRs in the UCYN-A1 symbiosis upon the addition of NO₃⁻, despite the lack of NO₃⁻ assimilation by the haptophyte host, has been experimentally demonstrated and is speculated to result from changes in the productivity or activity of the broader microbial community [29]. The production of organic matter, vitamins (e.g., B₁₂), or Fe-binding compounds by bacterioplankton and phytoplankton is well documented [76–78], and in environments such as the SCCS, which experience upwelling-stimulated net primary production, may be particularly important to the productivity of UCYN-A symbioses.

Averaged across the study region, depth-integrated NFR was markedly similar between the two study periods (May—195 ± 108 µmol N m⁻² d⁻¹, October—195 ± 189 µmol N m⁻² d⁻¹). Assuming Redfield C:N proportions for phytoplankton, this N would support the production of ca. 1.4 mmol C m⁻² d⁻¹. The climatological estimate (1970–2008) of new production in our study region was determined to be 17 and 1.7 mmol C m⁻² d⁻¹ in May and October, respectively [79]. Thus, at the time of this study, N₂ fixation accounted for ~8% of new production in May and 83% in October. The climatological estimate for total production in the region was greater and not substantially different between May (37.5 mmol C m⁻² d⁻¹) and October (27.5 mmol C m⁻² d⁻¹) [79]. As such, N₂ fixation-based production would only contribute to ~4–5% of total production. Future estimates of new production for the region are lower due to community composition shifts (less diatoms and more picophytoplankton) [80]. Such future changes would likely increase the role of N₂ fixation in supporting new production. Improved representation in ecological models of the magnitude and seasonality of nearshore N₂ fixation as well as the activity of coastal diazotrophs in the SCCS and other coastal regions will aid our understanding of and ability to predict changes to the balance between the global ocean’s N sources (N₂ fixation) and sinks (denitrification and anammox), particularly in coastally influenced waters rapidly undergoing climatological changes.

Contribution of UCYN-A symbioses to SCCS N₂ fixation

Our study indicates that UCYN-A symbioses are important N₂-fixers in the SCCS. This is supported by high UCYN-A abundances and low abundances of other diazotrophs, along with the highest single-cell NFR rates reported (UCYN-A2 symbiosis; 362 fmol N cell⁻¹ d⁻¹, Table S5), and average rates comparable to those reported from offshore regions (Fig. 6). Furthermore, UCYN-A symbioses can account for a substantial portion, if not all, of the bulk N₂ fixation at each station, with the exception of Stn. 1 (sampled May 3, prior to the upwelling event). Although our estimates are limited to a subset of stations and subject to large errors, N₂ fixation by the UCYN-A symbioses can account for between 6% and up to 100% of bulk N₂ fixation (Table 2). A substantial percentage of bulk N₂ fixation measured at offshore stations is attributed to UCYN-A1 symbioses, while the UCYN-A2 symbiosis appears to contribute predominantly to nearshore waters, which is a new line of evidence supporting the hypothesis that sublineages are different ecotypes [47]. This is most evident at Stn. 14 where potential contributions from the UCYN-A2 symbiosis were low (2 ± 7%) even though the UCYN-A1 symbiosis could account for bulk NFRs.

Estimates of the UCYN-A symbioses contribution to bulk rates are valuable but require cautious interpretation. The accuracy of cell abundance estimates using qPCR is impacted by numerous factors including DNA extraction efficiency, the quality of qPCR standards, and qPCR assay specificity/efficiency, and can be complicated by polyploidy, which is not well known in UCYN-A. Detecting ¹⁵N assimilation into single cells using nanoSIMS is subject to large uncertainties resulting from error associated with each measured parameter (e.g., $A_{P N}^{f i n a l}$ , $A_{P N}^{t = 0}$ ), the small number of data points, and underestimates of ¹⁵N incorporation due to isotope dilution [64]. More direct measurements of UCYN-A single-cell NFRs in other temperate coastal regions are needed, but our findings confirm their importance to N₂ fixation in coastal ecosystems.

It should be noted that heterocyst-forming Richelia associated with diatoms (diatom/diazotroph associations; DDAs) may contribute to nearshore N₂ fixation. They have been reported in the SCCS [29, 81] and Gulf of California [73, 82], and were detected at low abundances throughout the region in October, However, DDAs have some of the highest per cell NFRs reported [83, 84], thus, even at low abundances they have the potential to be substantive contributors to bulk N₂ fixation [85].

UCYN-A symbioses assemblages reflect mixing of nearshore and offshore ecotypes

SCCS UCYN-A assemblages were dominated by the UCYN-A1 symbiosis, even in nearshore stations, which was surprising given it was undetected in a 2010–2011 seasonal study in the Santa Catalina Bight [52]. Regional warming trends have been hypothesized to increase advection of oligotrophic waters into the Southern California Bight [68] and SST remained elevated in 2017 after the anomalous 2014–2016 warming trend and the 2015 El Niño (0.5–1 °C SST anomalies between 28–32° N; [86]), which have been linked to the persistence of warm Prochlorococcus ecotypes into 2018 [87] and may also explain the high abundances of the UCYN-A1 symbiosis in 2017.

In addition, UCYN-A oligotype assemblages reflect a mixture of nearshore and offshore populations in the SCCS. Ordination analysis on the global ocean dataset shows offshore samples clustered together and apart from nearshore samples (Fig. 4C) driven by the co-occurrence of UCYN-A1/UCYN-A3 in offshore samples and UCYN-A2/UCYN-A4 in nearshore samples (Fig. 4B) and are consistent with previously described patterns [47]. SCCS and Noumea Lagoon samples were found in each cluster, reflecting sample-specific similarities to both offshore and nearshore assemblages.

However, several oligotypes were found primarily in coastal waters influenced by the advection of oligotrophic waters. Many SCCS and Noumea Lagoon samples did not cluster with offshore and nearshore endpoints, and these dissimilarities are influenced by a mixture of UCYN-A1 (oligo45) and UCYN-A2 oligotypes (oligo46 and oligo40), along with oligo4 (Fig. 4B). Morphological and genetic studies of B. bigelowii suggest it is a complex of biologically distinct species, many of which can be found in coastal environments [88]. Thus, minor oligotypes such as oligo45, 46, and 40 may be affiliated with B. bigelowii genotypes distinct from known hosts. However, it is unclear whether all B. bigelowii genotypes form symbioses with UCYN-A and more work is needed to link oligotypes with genotypes.

Summary and conclusions

This work directly demonstrates that UCYN-A symbioses actively fix N₂ in the temperate, coastal SCCS ecosystem, where their cell-specific rates are high enough to account for whole community N₂ fixation. N₂ fixation by UCYN-A1 was more important in offshore waters, while UCYN-A2 contributed a higher share of the total N₂ fixation in nearshore waters, suggesting that even when these sublineages co-occur, their activity may be influenced by different environmental factors. More measurements in other coastal ecosystems are needed to better understand the variability and environmental controls on UCYN-A N₂ fixation, however, these findings lend support to speculation about its importance to N₂ fixation in other coastal systems [6, 7, 10, 89, 90]. In addition, although whole community N₂ fixation does not account for a large amount of new production in this system, this work establishes that N₂ fixation is a widespread process along the Baja California peninsula, during both upwelling and oceanic seasons and throughout the euphotic zone, with the highest volumetric rates in euphotic waters overlying the continental shelf. These measurements provide an important baseline for understanding how N₂ fixation and N-cycling will be impacted by future changes in climatology of the SCCS.

Acknowledgements

We gratefully acknowledge Sarah Spaugh (Stanford), Mehr Kumar (Stanford), Jenna Spackeen (VIMS), Cristian Hakspiel-Segura and Victor Camacho-Ibar (Universidad Autónoma de Baja California), Jose’ Antonia Jiménez Quiñóez (Armada de México), Meredith McPherson (UCSC), Mary Hogan (UCSC), Ana Maria Cabello Perez (UCSC/Instituto Espanol de Oceanografıa) for support processing samples in the field, as well as Josh Manger, James Convery and the captains and crew of the R/V Robert Gordon Sproul (SIO) for field logistical support. We also gratefully acknowledge David Karl (UH Manoa) for providing access to the MIMS, Mary Rosie Gradoville (UCSC) for helpful discussions on N₂ fixation rate calculations, Pete Raimondi (UCSC) for statistical advice, Kendra Negrey (UCSC), Corrine Gibble (UCSC) and Raphael Kudela (UCSC) for assistance with nutrient measurements, and Farooq Azam (SIO) for providing us with laboratory and field support. Finally, we thank Stanford Nano Shared Facilities and Chuck Hitzman, supported by the National Science Foundation under award ECCS-1542152, for nanoSIMS consultation, Stefan Green and his staff at the DNA Services Facility and the University of Illinois, Chicago for sequencing consultation, and the SOEST Biogeochemical Stable Isotope Facility at the University of Hawai’i, Manoa. This work was supported by grants from the National Science Foundation (OCE-1559165 to JPZ and 1559152 to KRA and 1756524 to STW) and the Simons Foundation (SCOPE #72440 and #329108 to JPZ and #824082 to JPZ).

Author contributions

MMM and KATK conceptualized the research and carried out formal analysis of the data; MMM, KATK, GVD, BAH, and BS collected field samples; MIMS analysis was conducted by STW; funding acquisition and supervision were carried out by KRA and JPZ; the original draft was written by KATK and MMM with input from all authors.

Competing interests

The authors declare no conflicts of interest.

Supplementary information

The online version contains supplementary material available at https://doi.org/10.1038/s43705-021-00039-7.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

1. Karl, D; Letelier, R; Tupas, L; Dore, J; Christian, J; Hebel, D. The role of nitrogen fixation in biogeochemical cycling in the subtropical North Pacific Ocean. Nature; 1997; 388, pp. 533-8.[COI: 1:CAS:528:DyaK2sXlt1art70%3D] [DOI: https://dx.doi.org/10.1038/41474]

2. Jickells, TD; Buitenhuis, E; Altieri, K; Baker, AR; Capone, D; Duce, RA et al. A reevaluation of the magnitude and impacts of anthropogenic atmospheric nitrogen inputs on the ocean. Glob Biogeochem Cycles.; 2017; 31, pp. 289-305.[COI: 1:CAS:528:DC%2BC2sXitlCrurw%3D]

3. Knapp, A. The sensitivity of marine N₂ fixation to dissolved inorganic nitrogen. Front Microbiol.; 2012; 3, 374.[COI: 1:CAS:528:DC%2BC3sXmtlKisbs%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23091472][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3476826]

4. Rees, AP; Gilbert, JA; Kelly-Gerreyn, BA. Nitrogen fixation in the western English Channel (NE Atlantic Ocean). Mar Ecol Prog Ser.; 2009; 374, pp. 7-12.[COI: 1:CAS:528:DC%2BD1MXjsVKgsbo%3D] [DOI: https://dx.doi.org/10.3354/meps07771]

5. Shiozaki, T; Nagata, T; Ijichi, M; Furuya, K. Nitrogen fixation and the diazotroph community in the temperate coastal region of the northwestern North Pacific. Biogeosciences; 2015; 12, pp. 4751-64. [DOI: https://dx.doi.org/10.5194/bg-12-4751-2015]

6. Tang, W; Cerdán-García, E; Berthelot, H; Polyviou, D; Wang, S; Baylay, A et al. New insights into the distributions of nitrogen fixation and diazotrophs revealed by high-resolution sensing and sampling methods. ISME J; 2020; 14, pp. 2514-26.[COI: 1:CAS:528:DC%2BB3cXht1GjsbzE] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32581316][DOI: https://dx.doi.org/10.1038/s41396-020-0703-6] [PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC32581316]

7. Tang, W; Wang, S; Fonseca-Batista, D; Dehairs, F; Gifford, S; Gonzalez, AG et al. Revisiting the distribution of oceanic N₂ fixation and estimating diazotrophic contribution to marine production. Nat Commun.; 2019; 10, pp. 1-10. [DOI: https://dx.doi.org/10.1038/s41467-018-07882-8] [COI: 1:CAS:528:DC%2BC1MXmvVCms78%3D]

8. Hamersley, M; Turk, K; Leinweber, A; Gruber, N; Zehr, J; Gunderson, T; Capone, D. Nitrogen fixation within the water column associated with two hypoxic basins in the Southern California Bight. Aquat Microb Ecol.; 2011; 63, pp. 193-205. [DOI: https://dx.doi.org/10.3354/ame01494]

9. Mulholland, MR; Bernhardt, PW; Blanco-Garcia, JL; Mannino, A; Hyde, K; Mondragon, E et al. Rates of dinitrogen fixation and the abundance of diazotrophs in North American coastal waters between Cape Hatteras and Georges Bank. Limnol Oceanogr.; 2012; 57, pp. 1067-83.[COI: 1:CAS:528:DC%2BC38XhtlCqs7fL] [DOI: https://dx.doi.org/10.4319/lo.2012.57.4.1067]

10. Mulholland, MR; Bernhardt, PW; Widner, BN; Selden, CR; Chappell, PD; Clayton, S et al. High rates of N₂ fixation in temperate, western North Atlantic coastal waters expands the realm of marine diazotrophy. Glob Biogeochem Cycles.; 2019; 33, pp. 826-40.[COI: 1:CAS:528:DC%2BC1MXhtlejt77K] [DOI: https://dx.doi.org/10.1029/2018GB006130]

11. Bentzon-Tilia, M; Traving, SJ; Mantikci, M; Knudsen-Leerbeck, H; Hansen, JL; Markager, S et al. Significant N₂ fixation by heterotrophs, photoheterotrophs and heterocystous cyanobacteria in two temperate estuaries. ISME J; 2015; 9, pp. 273-85.[COI: 1:CAS:528:DC%2BC2MXhsVaiu7c%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25026373][DOI: https://dx.doi.org/10.1038/ismej.2014.119] [PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC25026373]

12. Wen, Z; Lin, W; Shen, R; Hong, H; Kao, SJ; Shi, D. Nitrogen fixation in two coastal upwelling regions of the Taiwan Strait. Sci Rep.; 2017; 7, pp. 1-10. [DOI: https://dx.doi.org/10.1038/s41598-016-0028-x] [COI: 1:CAS:528:DC%2BC1cXhs1CltLvP]

13. Voss, M; Bombar, D; Loick, N; Dippner, JW. Riverine influence on nitrogen fixation in the upwelling region off Vietnam, South China Sea. Geophys Res Lett.; 2006; 33, L07604. [DOI: https://dx.doi.org/10.1029/2005GL025569] [COI: 1:CAS:528:DC%2BD28XlsFChsbg%3D]

14. Shiozaki T, Furuya K, Kodama T, Kitajima S, Takeda S, Takemura T, et al. New estimation of N₂ fixation in the western and central Pacific Ocean and its marginal seas. Glob Biogeochem Cycles. 2010;24:GB1015–n/a.

15. Blais, M; Tremblay, JÉ; Jungblut, AD; Gagnon, J; Martin, J; Thaler, M et al. Nitrogen fixation and identification of potential diazotrophs in the Canadian Arctic. Glob Biogeochem Cycles.; 2012; 26, GB3022. [DOI: https://dx.doi.org/10.1029/2011GB004096] [COI: 1:CAS:528:DC%2BC38XhsFSmt77L]

16. Shiozaki, T; Fujiwara, A; Inomura, K; Hirose, Y; Hashihama, F; Harada, N. Biological nitrogen fixation detected under Antarctic sea ice. Nat Geosci.; 2020; 13, pp. 729-32.[COI: 1:CAS:528:DC%2BB3cXitFOgtLnL] [DOI: https://dx.doi.org/10.1038/s41561-020-00651-7]

17. Harding, K; Turk-Kubo, KA; Sipler, RE; Mills, MM; Bronk, DA; Zehr, JP. Symbiotic unicellular cyanobacteria fix nitrogen in the Arctic Ocean. Proc Natl Acad Sci USA.; 2018; 115, pp. 13371-5.[COI: 1:CAS:528:DC%2BC1cXisF2rsrfO] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30538206][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6310837][DOI: https://dx.doi.org/10.1073/pnas.1813658115]

18. Thompson, AW; Foster, RA; Krupke, A; Carter, BJ; Musat, N; Vaulot, D et al. Unicellular cyanobacterium symbiotic with a single-celled eukaryotic alga. Science; 2012; 337, pp. 1546-50.[COI: 1:CAS:528:DC%2BC38XhtlGrur7O] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22997339][DOI: https://dx.doi.org/10.1126/science.1222700]

19. Zehr, JP; Shilova, IN; Farnelid, HM; del Carmen Muñoz-MarínCarmen, M; Turk-Kubo, KA. Unusual marine unicellular symbiosis with the nitrogen-fixing cyanobacterium UCYN-A. Nat Microbiol.; 2016; 2, 16214. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27996008][DOI: https://dx.doi.org/10.1038/nmicrobiol.2016.214] [COI: 1:CAS:528:DC%2BC2sXkvFyqsLg%3D]

20. Zehr, JP; Bench, SR; Carter, BJ; Hewson, I; Niazi, F; Shi, T et al. Globally distributed uncultivated oceanic N₂-fixing cyanobacteria lack oxygenic photosystem II. Science; 2008; 322, pp. 1110-2.[COI: 1:CAS:528:DC%2BD1cXhtlGgsr7P] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19008448][DOI: https://dx.doi.org/10.1126/science.1165340]

21. Tripp, HJ; Bench, SR; Turk, KA; Foster, RA; Desany, BA; Niazi, F et al. Metabolic streamlining in an open-ocean nitrogen-fixing cyanobacterium. Nature; 2010; 464, pp. 90-4.[COI: 1:CAS:528:DC%2BC3cXit1SmtLw%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20173737][DOI: https://dx.doi.org/10.1038/nature08786] [PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC20173737]

22. Church MJ, Mahaffey C, Letelier RM, Lukas R, Zehr JP, Karl DM. Physical forcing of nitrogen fixation and diazotroph community structure in the North Pacific subtropical gyre. Glob Biogeochem Cycles. 2009;23:GB2020.

23. Langlois, RJ; Hummer, D; LaRoche, J. Abundances and distributions of the dominant nifH phylotypes in the Northern Atlantic Ocean. Appl Environ Microbiol.; 2008; 74, pp. 1922-31.[COI: 1:CAS:528:DC%2BD1cXjs1ags7k%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18245263][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2268318][DOI: https://dx.doi.org/10.1128/AEM.01720-07]

24. Moisander, PH; Beinart, RA; Hewson, I; White, AE; Johnson, KS; Carlson, CA et al. Unicellular cyanobacterial distributions broaden the oceanic N₂ fixation domain. Science; 2010; 327, pp. 1512-4.[COI: 1:CAS:528:DC%2BC3cXjt1Gnu7w%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20185682][DOI: https://dx.doi.org/10.1126/science.1185468] [PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC20185682]

25. Krupke, A; Lavik, G; Halm, H; Fuchs, BM; Amann, RI; Kuypers, MM. Distribution of a consortium between unicellular algae and the N₂ fixing cyanobacterium UCYN-A in the North Atlantic Ocean. Environ Microbiol; 2014; 16, pp. 3153-67.[COI: 1:CAS:528:DC%2BC2cXhslCqtLvI] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24612325][DOI: https://dx.doi.org/10.1111/1462-2920.12431] [PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC24612325]

26. Shiozaki, T; Bombar, D; Riemann, L; Hashihama, F; Takeda, S; Yamaguchi, T et al. Basin scale variability of active diazotrophs and nitrogen fixation in the North Pacific, from the tropics to the subarctic Bering Sea. Glob Biogeochem Cycles; 2017; 31, pp. 996-1009.[COI: 1:CAS:528:DC%2BC2sXhtVKhtL%2FK] [DOI: https://dx.doi.org/10.1002/2017GB005681]

27. Krupke, A; Musat, N; Laroche, J; Mohr, W; Fuchs, BM; Amann, RI et al. In situ identification and N₂ and C fixation rates of uncultivated cyanobacteria populations. Syst Appl Microbiol.; 2013; 36, pp. 259-71.[COI: 1:CAS:528:DC%2BC3sXltVeku7w%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23541027][DOI: https://dx.doi.org/10.1016/j.syapm.2013.02.002] [PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC23541027]

28. Martínez-Pérez, C; Mohr, W; Löscher, CR; Dekaezemacker, J; Littmann, S; Yilmaz, P et al. The small unicellular diazotrophic symbiont, UCYN-A, is a key player in the marine nitrogen cycle. Nat Microbiol.; 2016; 1, 16163. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27617976][DOI: https://dx.doi.org/10.1038/nmicrobiol.2016.163] [COI: 1:CAS:528:DC%2BC2sXkvFyrt74%3D] [PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC27617976]

29. Mills, MM; Turk-Kubo, KA; van Dijken, GL; Henke, BA; Harding, K; Wilson, ST et al. Unusual marine cyanobacteria/haptophyte symbiosis relies on N₂ fixation even in N-rich environments. ISME J; 2020; 14, pp. 2395-406.[COI: 1:CAS:528:DC%2BB3cXhtFels7%2FK] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32523086][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7490277][DOI: https://dx.doi.org/10.1038/s41396-020-0691-6]

30. Scavotto, RE; Dziallas, C; Bentzon-Tilia, M; Riemann, L; Moisander, PH. Nitrogen-fixing bacteria associated with copepods in coastal waters of the North Atlantic Ocean. Environ. Microbiol.; 2015; 17, pp. 3754-65.[COI: 1:CAS:528:DC%2BC2MXhslCltbvM] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25655773][DOI: https://dx.doi.org/10.1111/1462-2920.12777] [PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC25655773]

31. Conroy, BJ; Steinberg, DK; Song, B; Kalmbach, A; Carpenter, EJ; Foster, RA. Mesozooplankton graze on cyanobacteria in the amazon river plume and western tropical North Atlantic. Front Microbiol.; 2017; 8, 1436. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28824569][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5540951][DOI: https://dx.doi.org/10.3389/fmicb.2017.01436]

32. Turk-Kubo, KA; Connell, P; Caron, D; Hogan, ME; Farnelid, HM; Zehr, JP. In situ diazotroph population dynamics under different resource ratios in the North Pacific Subtropical Gyre. Front Microbiol.; 2018; 9, 1616. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30090092][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6068237][DOI: https://dx.doi.org/10.3389/fmicb.2018.01616]

33. Needham, DM; Fuhrman, JA. Pronounced daily succession of phytoplankton, archaea and bacteria following a spring bloom. Nat Microbiol.; 2016; 1, 16005.[COI: 1:CAS:528:DC%2BC2sXkvF2jsb0%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27572439][DOI: https://dx.doi.org/10.1038/nmicrobiol.2016.5] [PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC27572439]

34. Shiozaki, T; Fujiwara, A; Ijichi, M; Harada, N; Nishino, S; Nishi, S et al. Diazotroph community structure and the role of nitrogen fixation in the nitrogen cycle in the Chukchi Sea (western Arctic Ocean). Limnol Oceanogr.; 2018; 63, pp. 2191-205.[COI: 1:CAS:528:DC%2BC1cXitFeis7zE] [DOI: https://dx.doi.org/10.1002/lno.10933]

35. Sohm, JA; Hilton, JA; Noble, AE; Zehr, JP; Saito, MA; Webb, EA. Nitrogen fixation in the South Atlantic Gyre and the Benguela Upwelling system. Geophys Res Lett.; 2011; 38, pp. L16608-n/a. [DOI: https://dx.doi.org/10.1029/2011GL048315] [COI: 1:CAS:528:DC%2BC38XhtFShurfK]

36. Moreira-Coello, V; Mouriño-Carballido, B; Marañón, E; Fernández-Carrera, A; Bode, A; Varela, MM. Biological N₂ fixation in the upwelling region off NW Iberia: magnitude, relevance, and players. Front Mar Sci; 2017; 4, 303. [DOI: https://dx.doi.org/10.3389/fmars.2017.00303]

37. Cabello, AM; Turk-Kubo, KA; Hayashi, K; Jacobs, L; Kudela, RM; Zehr, JP. Unexpected presence of the nitrogen-fixing symbiotic cyanobacterium UCYN-A in Monterey Bay, California. J Phycol.; 2020; 56, pp. 1521-33.[COI: 1:CAS:528:DC%2BB3MXotlGjuw%3D%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32609873][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7754506][DOI: https://dx.doi.org/10.1111/jpy.13045]

38. Deutsch, C; Frenzel, H; McWilliams, JC; Renault, L; Kessouri, F; Howard, E et al. Biogeochemical variability in the California Current System. Prog Oceanogr.; 2021; 196, 102565. [DOI: https://dx.doi.org/10.1016/j.pocean.2021.102565]

39. Grasshoff K, Kremling K, Ehrhardt M, editors. Methods of seawater analysis. 3rd ed. Weinheim: Wiley-VCH; 1999.

40. Welschmeyer, NA. Fluorometric analysis of chlorophyll a in the presence of chlorophyll b and phaeopigments. Limnol Oceanogr.; 1994; 39, pp. 1985-92.[COI: 1:CAS:528:DyaK2MXks1Sru70%3D] [DOI: https://dx.doi.org/10.4319/lo.1994.39.8.1985]

41. Moisander, PH; Beinart, RA; Voss, M; Zehr, JP. Diversity and abundance of diazotrophic microorganisms in the South China Sea during intermonsoon. ISME J; 2008; 2, pp. 954-67.[COI: 1:CAS:528:DC%2BD1cXhtF2qsrnL] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18528417][DOI: https://dx.doi.org/10.1038/ismej.2008.51] [PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC18528417]

42. Caporaso, JG; Kuczynski, J; Stombaugh, J; Bittinger, K; Bushman, FD; Costello, EK et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods; 2010; 7, pp. 335-6.[COI: 1:CAS:528:DC%2BC3cXksFalurg%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20383131][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3156573]

43. Edgar, RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics; 2010; 26, pp. 2460-1.[COI: 1:CAS:528:DC%2BC3cXht1WhtbzM] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20709691][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC20709691][DOI: https://dx.doi.org/10.1093/bioinformatics/btq461]

44. Edgar, RC; Haas, BJ; Clemente, JC; Quince, C; Knight, R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics; 2011; 27, pp. 2194-200.[COI: 1:CAS:528:DC%2BC3MXhtVSiurvL] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21700674][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC21700674][DOI: https://dx.doi.org/10.1093/bioinformatics/btr381]

45. Zehr, JP; Jenkins, BD; Short, SM; Steward, GF. Nitrogenase gene diversity and microbial community structure: a cross-system comparison. Environ Microbiol.; 2003; 5, pp. 539-54.[COI: 1:CAS:528:DC%2BD3sXlvFChs7s%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12823187][DOI: https://dx.doi.org/10.1046/j.1462-2920.2003.00451.x] [PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12823187]

46. Eren, AM; Maignien, L; Sul, WJ; Murphy, LG; Grim, SL; Morrison, HG et al. Oligotyping: differentiating between closely related microbial taxa using 16S rRNA gene data. Methods Ecol Evol; 2013; 4, pp. 1111-9. [PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3864673][DOI: https://dx.doi.org/10.1111/2041-210X.12114] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/3864673]

47. Turk-Kubo, KA; Farnelid, HM; Shilova, IN; Henke, B; Zehr, JP. Distinct ecological niches of marine symbiotic N₂-fixing cyanobacterium Candidatus Atelocyanobacterium thalassa sublineages. J Phycol.; 2017; 53, pp. 451-61.[COI: 1:CAS:528:DC%2BC2sXmtFClu7c%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27992651][DOI: https://dx.doi.org/10.1111/jpy.12505]

48. Henke, BA; Turk-Kubo, KA; Bonnet, S; Zehr, JP. Distributions and abundances of sublineages of the N₂-fixing cyanobacterium Candidatus Atelocyanobacterium thalassa (UCYN-A) in the New Caledonian Coral Lagoon. Front Microbiol.; 2018; 9, 554. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29674998][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5895702][DOI: https://dx.doi.org/10.3389/fmicb.2018.00554]

49. Gradoville, MR; Farnelid, H; White, AE; Turk‐Kubo, KA; Stewart, B; Ribalet, F et al. Latitudinal constraints on the abundance and activity of the cyanobacterium UCYN‐A and other marine diazotrophs in the North Pacific. Limnol Oceanogr.; 2020; 65, pp. 1858-75.[COI: 1:CAS:528:DC%2BB3cXis1Krs7jI] [DOI: https://dx.doi.org/10.1002/lno.11423]

50. McMurdie, PJ; Holmes, S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE; 2013; 8, e61217.[COI: 1:CAS:528:DC%2BC3sXntVWht7w%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/3632530][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3632530][DOI: https://dx.doi.org/10.1371/journal.pone.0061217]

51. Church, M; Jenkins, B; Karl, D; Zehr, J. Vertical distributions of nitrogen-fixing phylotypes at Stn ALOHA in the oligotrophic North Pacific Ocean. Aquat Microb Ecol; 2005; 38, pp. 3-14. [DOI: https://dx.doi.org/10.3354/ame038003]

52. Thompson, A; Carter, BJ; Turk-Kubo, K; Malfatti, F; Azam, F; Zehr, JP. Genetic diversity of the unicellular nitrogen-fixing cyanobacteria UCYN-A and its prymnesiophyte host. Environ Microbiol.; 2014; 16, pp. 3238-49.[COI: 1:CAS:528:DC%2BC2cXhslCqtLbM] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24761991][DOI: https://dx.doi.org/10.1111/1462-2920.12490] [PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC24761991]

53. Foster, RA; Subramaniam, A; Mahaffey, C; Carpenter, EJ; Capone, DG; Zehr, JP. Influence of the Amazon River plume on distributions of free-living and symbiotic cyanobacteria in the western tropical north Atlantic Ocean. Limnol Oceanogr.; 2007; 52, pp. 517-32.[COI: 1:CAS:528:DC%2BD2sXkt1Wgsr4%3D] [DOI: https://dx.doi.org/10.4319/lo.2007.52.2.0517]

54. Goebel, NL; Turk, KA; Achilles, KM; Paerl, R; Hewson, I; Morrison, AE et al. Abundance and distribution of major groups of diazotrophic cyanobacteria and their potential contribution to N₂ fixation in the tropical Atlantic Ocean. Environ Microbiol.; 2010; 12, pp. 3272-89.[COI: 1:CAS:528:DC%2BC3MXltleitA%3D%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20678117][DOI: https://dx.doi.org/10.1111/j.1462-2920.2010.02303.x] [PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC20678117]

55. Farnelid, H; Turk-Kubo, K; Munoz-Marin, MD; Zehr, JP. New insights into the ecology of the globally significant uncultured nitrogen-fixing symbiont UCYN-A. Aquat Microb Ecol.; 2016; 77, pp. 125-38. [DOI: https://dx.doi.org/10.3354/ame01794]

56. Mohr, W; Grosskopf, T; Wallace, DWR; LaRoche, J. Methodological underestimation of oceanic nitrogen fixation rates. PLoS ONE; 2010; 5, e12583. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20838446][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2933240][DOI: https://dx.doi.org/10.1371/journal.pone.0012583] [COI: 1:CAS:528:DC%2BC3cXhtFCksrfJ]

57. Montoya, JP; Voss, M; Kahler, P; Capone, DG. A simple, high-precision, high-sensitivity tracer assay for N₂ fixation. Appl Environ Microbiol.; 1996; 62, pp. 986-93.[COI: 1:CAS:528:DyaK28XhsVWqsb0%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16535283][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1388808][DOI: https://dx.doi.org/10.1128/aem.62.3.986-993.1996]

58. Gradoville, MR; Bombar, D; Crump, BC; Letelier, RM; Zehr, JP; White, AE. Diversity and activity of nitrogen-fixing communities across ocean basins. Limnol Oceanogr; 2017; 62, pp. 1895-909. [DOI: https://dx.doi.org/10.1002/lno.10542]

59. White, AE; Granger, J; Selden, C; Gradoville, MR; Potts, L; Bourbonnais, A et al. A critical review of the ¹⁵N₂ tracer method to measure diazotrophic production in pelagic ecosystems. Limnol Oceanogr Methods.; 2020; 18, pp. 129-47. [DOI: https://dx.doi.org/10.1002/lom3.10353]

60. Cornejo-Castillo, FM; Cabello, AM; Salazar, G; Sánchez-Baracaldo, P; Lima-Mendez, G; Hingamp, P et al. Cyanobacterial symbionts diverged in the late Cretaceous towards lineage-specific nitrogen fixation factories in single-celled phytoplankton. Nat Commun.; 2016; 7,[COI: 1:CAS:528:DC%2BC28Xks1CqtLc%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27002549][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4804200][DOI: https://dx.doi.org/10.1038/ncomms11071]

61. Cabello, AM; Cornejo-Castillo, FM; Raho, N; Blasco, D; Vidal, M; Audic, S et al. Global distribution and vertical patterns of a prymnesiophyte-cyanobacteria obligate symbiosis. ISME J; 2016; 10, pp. 693-706. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26405830][DOI: https://dx.doi.org/10.1038/ismej.2015.147]

62. Polerecky, L; Adam, B; Milucka, J; Musat, N; Vagner, T; Kuypers, MM. Look@ NanoSIMS–a tool for the analysis of nanoSIMS data in environmental microbiology. Environ Microbiol.; 2012; 14, pp. 1009-23.[COI: 1:CAS:528:DC%2BC38XmvFGgtrY%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22221878][DOI: https://dx.doi.org/10.1111/j.1462-2920.2011.02681.x]

63. Krupke, A; Mohr, W; LaRoche, J; Fuchs, BM; Amann, RI; Kuypers, MM. The effect of nutrients on carbon and nitrogen fixation by the UCYN-A-haptophyte symbiosis. ISME J; 2015; 9, pp. 1635-47.[COI: 1:CAS:528:DC%2BC2MXhtFOmt7jF] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25535939][DOI: https://dx.doi.org/10.1038/ismej.2014.253]

64. Meyer, NR; Fortney, J; Dekas, AE. NanoSIMS sample preparation decreases isotope enrichment: magnitude, variability and implications for single-cell rates of microbial activity. Environ Microbiol.; 2020; 23, pp. 81-98. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33000528][DOI: https://dx.doi.org/10.1111/1462-2920.15264] [COI: 1:CAS:528:DC%2BB3MXmsVegurs%3D]

65. Durazo, R. Seasonality of the transitional region of the California Current System off Baja California. J Geophys Res Oceans; 2015; 120, pp. 1173-96. [DOI: https://dx.doi.org/10.1002/2014JC010405]

66. Bakun A. Coastal upwelling indices, west coast of North America, 1946–71.

67. Redfield AC. On the proportions of organic derivatives in sea water and their relation to the composition of plankton. Vol. 1. Liverpool: University Press of Liverpool; 1934. p. 176–92.

68. Bograd, SJ; Schroeder, ID; Jacox, MG. A water mass history of the Southern California current system. Geophys. Res. Lett.; 2019; 46, pp. 6690-8. [DOI: https://dx.doi.org/10.1029/2019GL082685]

69. Langlois, R; Großkopf, T; Mills, M; Takeda, S; LaRoche, J. Widespread distribution and expression of gamma A (UMB), an uncultured, diazotrophic, γ-proteobacterial nifH phylotype. PLoS ONE; 2015; 10, e0128912. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26103055][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4477881][DOI: https://dx.doi.org/10.1371/journal.pone.0128912] [COI: 1:CAS:528:DC%2BC28Xns1Cqur8%3D]

70. Dekaezemacker, J; Bonnet, S; Grosso, O; Moutin, T; Bressac, M; Capone, DG. Evidence of active dinitrogen fixation in surface waters of the eastern tropical South Pacific during El Niño and La Niña events and evaluation of its potential nutrient controls. Glob Biogeochem Cycles; 2013; 27, pp. 768-79.[COI: 1:CAS:528:DC%2BC3sXhs1SgtrjO] [DOI: https://dx.doi.org/10.1002/gbc.20063]

71. Chen, M; Lu, Y; Jiao, N; Tian, J; Kao, SJ; Zhang, Y. Biogeographic drivers of diazotrophs in the western Pacific Ocean. Limnol Oceanogr.; 2019; 64, pp. 1403-21.[COI: 1:CAS:528:DC%2BC1MXhtFWgtLfP] [DOI: https://dx.doi.org/10.1002/lno.11123]

72. Turk, KA; Rees, AP; Zehr, JP; Pereira, N; Swift, P; Shelley, R et al. Nitrogen fixation and nitrogenase (nifH) expression in tropical waters of the eastern North Atlantic. ISME J; 2011; 5, pp. 1201-12.[COI: 1:CAS:528:DC%2BC3MXnvVaqsb8%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21228888][DOI: https://dx.doi.org/10.1038/ismej.2010.205] [PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC21228888]

73. White, AE; Foster, RA; Benitez-Nelson, CR; Masqué, P; Verdeny, E; Popp, BN et al. Nitrogen fixation in the Gulf of California and the Eastern Tropical North Pacific. Prog Oceanogr; 2013; 109, pp. 1-17. [DOI: https://dx.doi.org/10.1016/j.pocean.2012.09.002]

74. Selden, CR; Mulholland, MR; Bernhardt, PW; Widner, B; Macías‐Tapia, A; Ji, Q et al. Dinitrogen fixation across physico-chemical gradients of the eastern tropical North Pacific oxygen deficient zone. Glob Biogeochem Cycles.; 2019; 33, pp. 1187-202.[COI: 1:CAS:528:DC%2BC1MXhslOgtr7K] [DOI: https://dx.doi.org/10.1029/2019GB006242]

75. Sohm, JA; Webb, EA; Capone, DG. Emerging patterns of marine nitrogen fixation. Nat Rev Microbiol.; 2011; 9, pp. 499-508.[COI: 1:CAS:528:DC%2BC3MXnsFWrs7o%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21677685][DOI: https://dx.doi.org/10.1038/nrmicro2594] [PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC21677685]

76. Carlucci, A; Bowes, PM. Production of vitamin B12, thiamine, and biotin by phytoplankton. J Phycol.; 1970; 6, pp. 351-7.[COI: 1:CAS:528:DyaE3MXht1Wmu70%3D]

77. Gledhill, M; Buck, KN. The organic complexation of iron in the marine environment: a review. Front Microbiol.; 2012; 3, 69. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22403574][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3289268]

78. Biddanda, B; Benner, R. Carbon, nitrogen, and carbohydrate fluxes during the production of particulate and dissolved organic matter by marine phytoplankton. Limnol Oceanogr.; 1997; 42, pp. 506-18.[COI: 1:CAS:528:DyaK2sXmsF2ht7o%3D] [DOI: https://dx.doi.org/10.4319/lo.1997.42.3.0506]

79. Hernández de la Torre, B; Gaxiola Castro, G; Álvarez Borrego, S; Gallegos García, A; Aguirre Gómez, R. New organic carbon in front of the Baja California Peninsula: time series and climatology. Hidrobiológica; 2015; 25, pp. 74-85.

80. Xiu, P; Chai, F; Curchitser, EN; Castruccio, FS. Future changes in coastal upwelling ecosystems with global warming: the case of the California Current System. Sci. Rep.; 2018; 8, [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29434297][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5809506][DOI: https://dx.doi.org/10.1038/s41598-018-21247-7] [COI: 1:CAS:528:DC%2BC1cXhs1ClurvJ]

81. Kimor, B; Reid, F; Jordan, J. An unusual occurrence of Hemiaulus membranaceus Cleve (Bacillariophyceae) with Richelia intracelluaris Schmidt (Cyanophyceae) off the coast of Southern California. Phycologia; 1978; 17, pp. 162-6. [DOI: https://dx.doi.org/10.2216/i0031-8884-17-2-162.1]

82. White AE, Prahl FG, Letelier RM, Popp BN. Summer surface waters in the Gulf of California: Prime habitat for biological N₂ fixation. Glob Biogeochem Cycles. 2007;21:GB2017–n/a.

83. Pyle, AE; Johnson, AM; Villareal, TA. Isolation, growth, and nitrogen fixation rates of the Hemiaulus-Richelia (diatom-cyanobacterium) symbiosis in culture. PeerJ.; 2020; 8, e10115. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33083143][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7548074][DOI: https://dx.doi.org/10.7717/peerj.10115] [COI: 1:CAS:528:DC%2BB3MXhs1ehsLzO]

84. Foster, RA; Kuypers, MM; Vagner, T; Paerl, RW; Musat, N; Zehr, JP. Nitrogen fixation and transfer in open ocean diatom–cyanobacterial symbioses. ISME J; 2011; 5, pp. 1484-93.[COI: 1:CAS:528:DC%2BC3MXhtVKjsLjN] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21451586][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3160684][DOI: https://dx.doi.org/10.1038/ismej.2011.26]

85. Caputo, A; Nylander, JAA; Foster, RA. The genetic diversity and evolution of diatom-diazotroph associations highlights traits favoring symbiont integration. FEMS Microbiol Lett.; 2019; 366, fny297.[COI: 1:CAS:528:DC%2BC1MXhvV2qsbvK] [PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6341774][DOI: https://dx.doi.org/10.1093/femsle/fny297]

86. Thompson AR. State of the California Current 2017–18: still not quite normal in the north and getting interesting in the south. California cooperative oceanic fisheries investigations, Data report. 2018.

87. Larkin, AA; Moreno, AR; Fagan, AJ; Fowlds, A; Ruiz, A; Martiny, AC. Persistent El Nino driven shifts in marine cyanobacteria populations. PLoS ONE; 2020; 15, e0238405.[COI: 1:CAS:528:DC%2BB3cXhvFWitbvP] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32936809][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7494125][DOI: https://dx.doi.org/10.1371/journal.pone.0238405]

88. Hagino, K; Takano, Y; Horiguchi, T. Pseudo-cryptic speciation in Braarudosphaera bigelowii (Gran and Braarud) Deflandre. Mar Micropaleontol.; 2009; 72, pp. 210-21. [DOI: https://dx.doi.org/10.1016/j.marmicro.2009.06.001]

89. Selden, CR; Chappell, PD; Clayton, S; Macías‐Tapia, A; Bernhardt, PW; Mulholland, MR. A coastal N₂ fixation hotspot at the Cape Hatteras front: elucidating spatial heterogeneity in diazotroph activity via supervised machine learning. Limnol Oceanogr.; 2021; 66, pp. 1832-49. [DOI: https://dx.doi.org/10.1002/lno.11727]

90. Wang, S; Tang, W; Delage, E; Gifford, S; Whitby, H; González, AG et al. Investigating the microbial ecology of coastal hotspots of marine nitrogen fixation in the western North Atlantic. Sci Rep.; 2021; 11,[COI: 1:CAS:528:DC%2BB3MXmtlCntbc%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33750865][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7943828][DOI: https://dx.doi.org/10.1038/s41598-021-84969-1]

Word count: 10039

Show less

© The Author(s) 2021. This work is published under http://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.

Abstract

Translate

The availability of fixed nitrogen (N) is an important factor limiting biological productivity in the oceans. In coastal waters, high dissolved inorganic N concentrations were historically thought to inhibit dinitrogen (N₂) fixation, however, recent N₂ fixation measurements and the presence of the N₂-fixing UCYN-A/haptophyte symbiosis in nearshore waters challenge this paradigm. We characterized the contribution of UCYN-A symbioses to nearshore N₂ fixation in the Southern California Current System (SCCS) by measuring bulk community and single-cell N₂ fixation rates, as well as diazotroph community composition and abundance. UCYN-A1 and UCYN-A2 symbioses dominated diazotroph communities throughout the region during upwelling and oceanic seasons. Bulk N₂ fixation was detected in most surface samples, with rates up to 23.0 ± 3.8 nmol N l⁻¹ d⁻¹, and was often detected at the deep chlorophyll maximum in the presence of nitrate (>1 µM). UCYN-A2 symbiosis N₂ fixation rates were higher (151.1 ± 112.7 fmol N cell⁻¹ d⁻¹) than the UCYN-A1 symbiosis (6.6 ± 8.8 fmol N cell⁻¹ d⁻¹). N₂ fixation by the UCYN-A1 symbiosis accounted for a majority of the measured bulk rates at two offshore stations, while the UCYN-A2 symbiosis was an important contributor in three nearshore stations. This report of active UCYN-A symbioses and broad mesoscale distribution patterns establishes UCYN-A symbioses as the dominant diazotrophs in the SCCS, where heterocyst-forming and unicellular cyanobacteria are less prevalent, and provides evidence that the two dominant UCYN-A sublineages are separate ecotypes.

Details

Title

UCYN-A/haptophyte symbioses dominate N₂ fixation in the Southern California Current System

Author

Turk-Kubo, Kendra A.¹

; Mills, Matthew M.²

; Arrigo, Kevin R.²

; van Dijken, Gert²

; Henke, Britt A.¹; Stewart, Brittany³; Wilson, Samuel T.⁴

; Zehr, Jonathan P.¹

¹ University of California at Santa Cruz, Ocean Sciences Department, Santa Cruz, USA (GRID:grid.205975.c) (ISNI:0000 0001 0740 6917)
² Stanford University, Earth System Science, Stanford, USA (GRID:grid.168010.e) (ISNI:0000000419368956)
³ University of California at Santa Cruz, Ocean Sciences Department, Santa Cruz, USA (GRID:grid.205975.c) (ISNI:0000 0001 0740 6917); University of Southern California, Department of Biological Sciences, Los Angeles, USA (GRID:grid.42505.36) (ISNI:0000 0001 2156 6853)
⁴ University of Hawai’i at Manoa, Center for Microbial Oceanography: Research and Education, Honolulu, USA (GRID:grid.410445.0) (ISNI:0000 0001 2188 0957)

Publication year

2021

Publication date

Dec 2021

Publisher

Oxford University Press

e-ISSN

27306151

Source type

Scholarly Journal

Language of publication

English

DOI

https://doi.org/10.1038/s43705-021-00039-7

ProQuest document ID

2730348537

UCYN-A/haptophyte symbioses dominate N₂ fixation in the Southern California Current System

Jump to:

Full Text

Abstract

Details

Suggested sources

UCYN-A/haptophyte symbioses dominate N2 fixation in the Southern California Current System

Jump to:

Full Text

Abstract

Details

UCYN-A/haptophyte symbioses dominate N₂ fixation in the Southern California Current System