Incubation experiments were conducted on two separate cruises: within regions of the CUZ during July 3–26th 2014 onboard the R/V Melville and along the Line-P transect of the subarctic NE Pacific Ocean during June 7–23rd 2015 onboard the Canadian Coast Guard Ship (CCGS) John P. Tully (Figure 1). The incubated phytoplankton community response was assessed using a combination of physiological measurements and metatranscriptomics to examine the effects of Fe status on diatom physiology and gene expression. Each experiment included three treatments: (1) a 5 nmol L⁻¹ FeCl₃ addition (Fe), (2) a 200 nmol L⁻¹ desferroxamine B (DFB) addition, and (3) an unamended control (Ctl), each sampled at two time points.

During the CUZ cruise, three incubation experiments were performed at separate locations corresponding to distinct Fe and macronutrient regimes (Supplementary Table 1), including sites of high dFe, macronutrients, and phytoplankton biomass (C1: 38°39.30 N, 123°39.87 W), relatively low dFe, high macronutrients and high phytoplankton biomass (C2: 38°15.31 N, 123°57.98 W), and low dFe with high macronutrients and low phytoplankton biomass (C3: 42°40.00 N, 125°0.00 W) (Figure 1). Near-surface seawater was collected from a depth of ~15 m using a trace-metal clean sampling system consisting of a tow-fish sampler attached to Kevlar™ line, PFA Teflon tubing, and a Teflon dual-diaphragm pump that pumped seawater directly into a positive pressure trace-metal clean bubble. The seawater was placed in a large 200 L acid-cleaned HDPE drum for homogenization before being distributed into 10 L flexible acid-cleaned polyethylene cubitainers (Hedwin Corporation). Cleaning protocols for the cubitainers included successive soaks in 1.2 mol L⁻¹ hydrochloric acid (reagent grade) for 3 days, 1.2 mol L⁻¹ hydrochloric acid (trace metal grade) for 1 week and 0.1 mol L⁻¹ acetic acid (trace-metal grade) until use. Prior to filling the cubitainers with seawater, the dilute acetic acid was removed and the cubitainers were rinsed thoroughly three times with ambient seawater from the collection site. The primary objective of these experiments was to elucidate the responses of target diatom genera and the phytoplankton community to variable Fe conditions. Therefore, sites were targeted that would ensure adequate macronutrient concentrations to support phytoplankton growth. However, at C2, 15 μmol L⁻¹ of Si(OH)₄ was added to all cubitainers to support growth of diatoms due to the initially low Si(OH)₄ concentration (<4.7 μmol L⁻¹).

During the Line-P cruise, incubation experiments were conducted at the low NO3- coastal station P4 (48°39 N, 126°40 W; referred to as C4 in this analysis) and at the chronically Fe-limited, HNLC oceanic station P26, also known as Ocean Station Papa (OSP, 50°00 N, 145°00 W; Harrison, 2002; referred to as O5). Seawater was collected at depths corresponding to ~30% of incident irradiance (8–12 m) at both stations using a trace-metal clean sampling system consisting of a Teflon air bellows pump and PTFE lined Kevlar™ tubing attached to a Kevlar™ line. The seawater was pumped directly into 10 L acid-cleaned polyethylene cubitainers placed within an on-deck trace-metal clean positive pressure flowhood. At site C4, 10 μmol L⁻¹ of NO3- was added to all cubitainers to support growth of diatoms due to the initially low NO3- concentration (<1.5 μmol L⁻¹).

At the start of the experiments, ambient seawater was filtered for all initial measurements (T₀). For each incubation experiment, cubitainers were filled to serve as a control (Ctl) or amended with FeCl₃ or DFB just prior to dawn. All cubitainers were placed in on-deck Plexiglass incubators with flow-through seawater to maintain near-ambient surface water temperatures. Incubators were covered with neutral density screening to achieve ~30% of the incident irradiance (Supplementary Figure 1). Following 24–96 h of incubation (Supplementary Table 1; depending on the measured macronutrient drawdown) the cubitainers for a specific time point were removed from the incubators and filtered immediately. The goal for each time point was to achieve measureable drawdowns in macronutrients that would infer stimulation of phytoplankton growth without complete macronutrient depletion. However, for some experiments and time points, depletion of NO3- or other macronutrients occurred. All filtrations were conducted at dawn. Subsamples for dissolved and particulate nutrients, size-fractionated uptake rates of dissolved inorganic carbon (DIC) and NO3-, size-fractionated chlorophyll a, F_v/F_m, and RNA were collected at T₀ and from each cubitainer according to the protocols described below.

For CUZ experiments, dissolved nitrate and nitrite (NO3- + NO2-), phosphate (PO43-), and silicic acid (H₄SiO₄) concentrations were measured onboard using a Lachat Quick Chem 8000 Flow Injection Analysis system (Parsons et al., 1984) with detection limits of 0.05 μM for NO3- + NO2-, 0.03 μM for PO43-, and 0.2 μM for H₄SiO₄ (Bruland et al., 2008). Particles were removed by filtration through a Whatman GF/F filter (25 mm). Reference standards for nutrients in seawater were run for quality control. During Line-P sampling, ~15 mL of seawater was filtered through a Whatman GF/F filter into acid-rinsed polypropylene tubes and frozen at −20°C in aluminum blocks until onshore analysis. Shortly following the cruise, the dissolved NO3- + NO2-, PO43-, and H₄SiO₄ concentrations were determined using an Astoria nutrient analyzer (Barwell-Clarke and Whitney, 1996). Nutrient detection limits were 0.2 μM for NO3- + NO2-, 0.02 μM for PO43-, and 0.5 μM for H₄SiO₄ (Frank Whitney and Mark Belton [IOS], pers. comm.).

For biogenic silica (bSi) measurements, 335 mL (CUZ) or 250 mL (Line P) of seawater was filtered onto polycarbonate filters (1.2 μm pore size for CUZ and 0.6 μm pore size for Line-P, 25 mm), digested with NaOH in Teflon tubes, and measured with the colorimetric ammonium molybdate method (Krause et al., 2013).

Size-fractionated particulate nitrogen (PN), particulate carbon (PC), and NO3- uptake rates were obtained by adding ¹⁵N-NaNO₃ to 618 mL subsample of experimental seawater placed within clear polycarbonate bottles. The concentration of NO3- added was no more than 10% of ambient NO3- level within CUZ incubations, and was 1 μmol L⁻¹ within Line-P incubations (corresponding to NO3- levels of 68% at T₀ and 10% within NO3--amended incubations at C4, and ~10% at O5). DIC uptake within Line-P incubations was measured by additionally spiking subsamples with 120 μmol L⁻¹ NaH¹³CO₃. Bottles were incubated in the same flow-through Plexiglass incubators where cubitainers were kept. Following 8 h of incubation, seawater samples were filtered in series through a polycarbonate filter (5 μm pore size, 47 mm) via gravity filtration, and then through a pre-combusted (450°C for 5 h) GF/F filter by gentle vacuum (<100 mg Hg). Particulates collected on the 5 μm polycarbonate filter were then rinsed onto a separate pre-combusted GF/F filter using an artificial saline solution. Filters were stored at −20°C until onshore analysis. In the laboratory, filters were heated to 60°C for 24 h and pelletized in tin capsules (Elemental Microanalysis) in preparation for analysis of the atom % ¹⁵N, atom % ¹³C (for Line-P), particulate nitrogen (PN), and particulate carbon (PC) using an elemental analyzer paired with an isotope ratio mass spectrometer (EA-IRMS). Biomass-normalized NO3- uptake rates (PN-VNO₃) and DIC uptake rates (PC-VDIC) for the Line-P experiments were obtained by dividing the measured NO3- and DIC biological uptake rates by PN and PC concentrations, respectively.

To quantify VDIC in CUZ incubations, incorporation of ¹⁴C was determined using a protocol adapted from Taylor et al. (2013). Briefly, 60 mL of seawater from each cubitainer was distributed into acid-cleaned light and dark polycarbonate bottles. In each bottle, 1.2 μCi of NaH¹⁴CO₃ was added. Bottles were incubated in the same flow-through Plexiglass incubators where cubitainers were kept for 6.5–8 h. Following incubation, samples were filtered through stacked 47 mm polycarbonate filters (5 and 1 μm) separated with a mesh spacer during filtration. Filters were vacuum dried, placed in 7 mL scintillation vials containing 0.5 mL of 6 M HCl and permitted to degas for 24 h. Disintegrations per minute (DPM) retained on the filters were measured using a Beckman Coulter LS 6500 scintillation counter. Reported values are light bottle DPMs minus dark bottle DPMs. To obtain VDIC, DIC uptake rates were normalized to PC concentrations obtained as part of the NO3- uptake measurements within each incubation and size fraction.

Seawater samples for Fe analysis within the CUZ were acidified at sea with the equivalent of 4 mL 6 N quartz-distilled HCl per L of seawater (pH ~ 1.7) and stored in acid-cleaned LDPE bottles for at least 2 months prior to analysis. Samples were analyzed using an adaption of Biller and Bruland (2012) as described in Parker et al. (2016). Briefly, this method involves preconcentrating the Fe from buffered (pH 6.0) seawater on Nobias-chelate PA1 resin and eluting with 1 N quartz-distilled HNO_3. The eluent was analyzed with a Thermo-Element high resolution XR ICP-MS in counting mode. Line-P dissolved Fe samples were stored in acid-cleaned LDPE bottles, acidified post-cruise with Optima-grade HCl (1 mL 12 N HCl per L of seawater), and allowed to sit for >3 months. Dissolved Fe was measured via ICP-MS by P. Morton at Florida State University following resin preconcentration using the protocol of Milne et al. (2010).

Four hundred milliliters of seawater was gravity-filtered through a polycarbonate filter (5 μm pore size, 47 mm diameter) followed by vacuum filtration through a GF/F filter (0.7 μm nominal pore size, 25 mm diameter) using a series filter cascade for size fractionation. Filters were frozen at −80°C until analysis. Chlorophyll a extraction was performed using 90% acetone at −20°C for 24 h and the extracted Chl a was quantified by fluorometry with a Turner Designs 10-AU fluorometer using the acidification method (Parsons et al., 1984).

Approximately 250 mL of seawater from each CUZ site was filtered through GF/F filters (25 mm) via vacuum pressure (<100 mm Hg) and the filters were frozen at −80°C. The filters were extracted with 20% methanol (MeOH) in water. The mixture was sonicated in an ice bath for 2 min at 30–40 W with a Sonicator 3000, followed by centrifugation (10 min, 1,399 × g). The supernatant was collected and passed through a 0.22 μm syringe filter. Samples were stored at −20°C until analysis. Concentrations with a detection limit of 0.01 μg L⁻¹ were obtained using an enzyme-linked immunosorbent assay (Abraxis, Warminster, PA, USA) following the manufacturer's protocol, including running each sample in duplicate at several dilutions. Final concentrations (pg DA mL extract⁻¹) were calculated using the manufacturer supplied analysis spreadsheet.

The maximum photochemical yield of PSII (F_v/F_m) was measured by fast repetition rate fluorometry (FRRF) using a custom-built fluorescence-induction and relaxation system (Kolber et al., 1998; Gorbunov and Falkowski, 2004). Before each measurement, a 5 mL subsample of seawater from each cubitainer was acclimated to low light for 20 min. A saturating pulse (20,000 μmol photons m⁻² s⁻¹) of blue light (450 nm) was applied to dark-acclimated cells for a duration of 100–200 μs. Measurements were obtained using the single-turnover flash (STF) setting with the average of 50 iterations for the CUZ experiments, and a single iteration for the Line-P experiments. Data were blank corrected using 0.2 μm filtered seawater.

Phytoplankton in seawater samples were filtered onto 0.8 μm Pall Supor filters (142 mm) via peristaltic pumping, immediately flash frozen in liquid nitrogen and stored at −80°C until extraction onshore. The filters were briefly thawed on ice before being extracted individually using the ToTALLY RNA Kit (Ambion). The extraction procedure followed manufacturer protocols with the modified first step of glass bead addition and vortexing to facilitate disruption of cells. Removal of DNA was performed with one round of DNAse I (Ambion) incubation. For the Line P experiments, due to low yields in treatments, RNA from the triplicate cubitainers was pooled prior to sequencing. Within CUZ experiments all triplicate incubation samples were sequenced separately. At the oceanic site O5, RNA yields were too low to successfully sequence metatranscriptomes at the T₁ timepoint, and consequently, transcriptomic analyses were performed using the T₀, T₂ Fe, and T₂ Ctl treatments. Metatranscriptomic library preparation was performed with the Illumina TruSeq Stranded mRNA Library Preparation Kit and HiSeq v4 reagents. Samples were barcoded and run across three lanes of Illumina HiSeq 2000 (125 bp, paired-end) yielding on average 23 million paired-end reads per sample (Supplementary Table 2). The RNA-seq data reported here has been deposited in the National Center for Biotechnology (NCBI) sequence read archive (SRA) under the BioProject accession no. PRJNA320398 and PRJNA388329.

Raw reads were trimmed for quality bases and removal of adapters using Trimmomatic v0.32 (paired-end mode, adaptive quality trim with 40 bp target length, and strictness of 0.6, minimum length of 36 bp; Bolger et al., 2014). Trimmed paired reads were merged into single reads with BBMerge v8.0. For each site, the resulting merged pairs and non-overlapping paired-end reads were assembled using ABySS v1.5.2 with a multi-kmer approach (Birol et al., 2009). The different k-mer assemblies were merged to remove redundant contigs using Trans-ABySS v1.5.3 (Robertson et al., 2010). Read counts were obtained by mapping raw reads to assembled contigs with Bowtie2 v2.2.6 (end-to-end alignment; Langmead and Salzberg, 2012). Alignments were filtered by mapping quality score (MAPQ) of 10 or higher as determined by SAMtools v1.2 (Li et al., 2009). Taxonomic and functional annotations were assigned based on sequence homology to reference databases via BLASTx v2.3.0 with an e-value cutoff of 10⁻³ (Altschul et al., 1990). Functional annotations were assigned according to the top hit using the Kyoto Encyclopedia of Genes and Genomes (KEGG; Release 75), while taxonomic assignments were obtained according to the top hit using MarineRefII (Laboratory of Mary Ann Moran, University of Georgia), a custom-made database comprised of protein sequences of marine prokaryotes and eukaryotes including all sequenced transcriptomes from Marine Microbial Eukaryote Transcriptome Sequencing Project (MMETSP) (Keeling et al., 2014). Taxonomic information was obtained from NCBI's Taxonomy Database (each isolate in MarineRefII is assigned a NCBI taxonomic ID). The information from NCBI was manually curated to ensure proper assignment and use of common phytoplankton taxonomic ranks. For our analysis, we have grouped diatom-associated sequences at the genus level. Therefore, the patterns in gene expression observed could be driven by one dominant species or many equally distributed species belonging to a genus within each site.

All diatom-assigned counts were summed to both the genus taxonomic rank and KEGG Orthology (KO) functional annotation level. For genes of interest without a KO assignments but with an annotated gene definition (i.e., ISIPs and rhodopsin), raw counts corresponding to KEGG gene definitions were summed. EdgeR v3.12.0 was used to calculate Pseudo-nitzschia- or Thalassiosira-specific normalized fold change and counts-per-million (CPM) from pairwise comparisons using the exactTest (Robinson and Smyth, 2008; Robinson and Oshlack, 2010; Robinson et al., 2010; Klingenberg and Meinicke, 2017). Significance (p < 0.05) was calculated with edgeR's estimate of tagwise dispersions across all samples within CUZ sites. Heatmaps were produced with the R package pheatmap v1.0.8, and dendrograms created using Euclidean distance and hierarchical clustering. Assembled contigs, read counts, and functional annotations of contigs are available at marchettilab.web.unc.edu/data.

In order to directly compare transcript abundance across locations for principal component analyses (PCA), the assemblies for all sites were merged with Trans-ABySS. The removal of redundant contigs was verified with GenomeTools v.1.5.1. Counts were obtained by aligning raw reads against this merged metatranscriptome using Salmon v.0.7.3-beta. Normalized counts were then obtained with edgeR v3.12.0. PCA biplots were created using log-transformed normalized counts for genes of interest with ggbiplot v.0.5.

Environmental Pseudo-nitzschia and Thalassiosira contigs functionally annotated as RubisCO (RBCL), rhodopsin (RHO), or superoxide dismutase (SOD) and containing a large number of mapped reads were compared to diatom reference sequences for phylogenetic characterization. Diatom sequences used in reference alignments were obtained through a sequence homology search using BLASTx v2.2.28 with Pseudo-nitzschia RBCL, RHO, and SOD against the database MMETSP using an E-value cutoff of 10⁻⁵ (Altschul et al., 1990). Sequences were aligned using MUSCLE within Geneious v5.6.4 software (Edgar, 2004).

CUZ site C1 (Figure 1) was characterized by high macronutrient and dFe concentrations in the mixed layer supporting a high biomass, nutrient-replete phytoplankton community. The community was dominated by phytoplankton cells in the >5 μm chlorophyll a (chl a) size fraction, constituting 88% of the total chl a concentration (Figure 2B; Supplementary Table 1). Macronutrient concentrations were rapidly consumed during the first 24 h of incubation (T₁), with near complete depletion of the NO3- (≤1 μmol L⁻¹ remaining by 48 h [T₂]; Figure 2A). The initially Fe-replete phytoplankton community (dFe: 3.57 nmol L⁻¹) was mostly unaffected by the additions of Fe or DFB as demonstrated through relatively constant F_v/F_m, phytoplankton biomass, particulate nitrogen (PN)-specific nitrate uptake rates (VNO₃, or nitrate assimilation rates), and particulate carbon (PC)-specific dissolved inorganic carbon uptake rates (VDIC, or carbon assimilation rates) across treatments at each time point (Figures 2B–E). Furthermore, the NO₃:Fe ratio of the initial (T₀) seawater (3.8 μmol:nmol. Supplementary Table 1) was substantially below the predicted threshold ratio for eventual Fe stress of 12 μmol:nmol for phytoplankton in this region as proposed by King and Barbeau (2007), albeit this ratio is subject to variation as a function of phytoplankton iron demands (Bruland et al., 2001), suggesting this phytoplankton community was not likely to be driven into Fe limitation prior to complete NO3- utilization. However, indications of molecular-level responses to Fe and DFB additions were observed; 74 genes were differentially expressed (p < 0.05) in Pseudo-nitzschia between the Fe and DFB treatments (Supplementary Figure 2A). Fe-stress bioindicator genes (FLDA, PETE, and ISIP2A; Whitney et al., 2011; Morrissey et al., 2015; Graff van Creveld et al., 2016) increased in expression following the addition of DFB relative to the added Fe treatment, suggesting the onset of Fe stress following the addition of DFB by the end of the first time point.

CUZ site C2 was located in close geographical proximity to C1 (Figure 1), yet exhibited different mixed layer properties in relation to phytoplankton biomass, silicic acid (Si[OH]₄) and dFe concentrations (0.44 nmol L⁻¹). Nitrate and ortho-phosphate (PO43-) concentrations were similarly high (10.3 and 0.96 μmol L⁻¹, respectively) as found at site C1, although Si(OH)₄ levels were appreciably lower (4.7 μmol L⁻¹) and possibly growth-limiting to certain diatoms (Nelson et al., 1996). Therefore, incubations were amended with 15 μmol L⁻¹ Si(OH)₄ to support potential diatom growth with added Fe (Brzezinski, 1985). Although the chl a concentration in the >5 μm size fraction was initially <1 μg L⁻¹ and biogenic silica (bSi) concentrations were <3 μmol L⁻¹, by 48 h (T₁) the >5 μm chl a fraction reached 5–8 μg L⁻¹, and bSi increased to 10–15 μmol L⁻¹ in all treatments, accompanied by appreciable decreases in NO3-, PO43-, and Si(OH)₄ concentrations (Figures 2A,B). Since this community quickly depleted NO3- concentrations during the experimental period, this site presented an opportunity to couple the physiological indicators of NO3- stress with N-related transport and assimilation genes observed to be elevated in NO₃-starved laboratory diatom cultures (Hildebrand, 2005; Song and Ward, 2007; Bender et al., 2014; Rogato et al., 2015). Apart from F_v/F_m reaching relatively low values in the DFB treatments, indications of Fe stress in bulk physiological measurements across treatments were absent (Figure 2C). However, the initial seawater NO₃:Fe ratio of 23.4 μmol:nmol suggests this community may have been driven into Fe limitation provided sufficient Si(OH)₄ was present. Additionally, a total of 414 Pseudo-nitzschia–associated genes were differentially expressed (p < 0.05) by T₁ between the Fe and DFB treatments (Supplementary Figure 2). This greater number of differentially expressed genes in Pseudo-nitzschia when compared to C1 suggests the C2 diatom community in the DFB treatment experienced a higher degree of Fe stress during the incubation period. The initially low dissolved Si(OH)₄:NO₃ ratio at this site furthermore implies a possible increase in the Si:N ratios of Fe-stressed diatoms (Hutchins and Bruland, 1998; Marchetti and Cassar, 2009; Brzezinski et al., 2015). Interestingly, concentrations of domoic acid (DA), a neurotoxin produced by Pseudo-nitzschia, was 90 pg mL⁻¹ in initial seawater (T₀) and exceeded 3,000 pg mL⁻¹ in the control treatment by T₁ (Supplementary Figure 3). This increase in DA concentration may be linked to both the increase in Pseudo-nitzschia abundance and depletion of Si(OH)₄ resulting in Si-limited cells which has been shown to greatly enhance DA production (Pan et al., 1996).

Site C3 (Figure 1) contained the lowest dFe concentrations (0.31 nmol L⁻¹) among the CUZ sites along with high macronutrient concentrations [17 μmol L⁻¹ NO3-, 19 μmol L⁻¹ Si(OH)₄, and 1.5 μmol L⁻¹ PO43-; Figure 2A]. The corresponding NO₃:Fe ratio of the initial seawater was ~54.9 μmol:nmol (Supplementary Table 1). Following incubation, the chl a, bSi, PN-specific VNO₃, and PC-specific VDIC were all higher in the Fe-amended treatment relative to the unamended control by T₁ (Figures 2B,D,E). By 72 h, NO3- was completely drawn down within the Fe treatment (T₂). Despite the pronounced influence of Fe enrichment on bulk parameters, F_v/F_m values were only slightly higher in the Fe treatment than the control, but they were substantially higher than in the DFB treatment (Figure 2C). This is likely a reflection of the different phytoplankton composition at this location compared to site C2, which did not show indications of an Fe-addition response on the measured bulk parameters, but did demonstrate elevated F_v/F_m values in the added Fe treatment. Site C3 represented the only phytoplankton community in the CUZ that displayed a definite physiological response to Fe addition relative to the control treatment (Supplementary Table 1). The Fe-induced molecular response in diatoms was demonstrated by the differential expression of 458 genes in Pseudo-nitzschia and 1,223 genes in Thalassiosira between the Fe and DFB treatments (Supplementary Figure 2C), and 365 genes in Pseudo-nitzschia and 837 genes in Thalassiosira between the Fe and Ctl treatments (p < 0.05).

Coastal site C4 was located at station P4 of the Line-P transect in the subarctic NE Pacific Ocean (Figure 1). Initial mixed-layer seawater properties were characterized by low concentrations of macronutrients and dFe, which supported a low phytoplankton biomass. Nitrate concentrations were initially 1.5 μmol L⁻¹ (Figure 2A). To facilitate a potential phytoplankton growth response to added Fe, 10 μmol L⁻¹ of NO3- was added to each treatment. Si(OH)₄ concentrations were also initially low (2.2 μmol L⁻¹) and incubation concentrations dropped to <2 μmol L⁻¹ in most treatments by the second time point (T₂; Figure 2A). These low concentrations restricted biomass accumulation as bSi (Figure 2B) and it is likely that the resulting diatom community experienced Si(OH)₄ limitation by the end of the incubation period. Despite its relatively close proximity to land and relatively high dFe concentration (0.64 nmol L⁻¹), there was a pronounced response to Fe addition at C4 as demonstrated through higher F_v/F_m, PN-specific VNO₃, and PC-specific VDIC in the Fe treatment compared to values in the unamended control by T₁ (Figures 2D,E; Supplementary Table 1). The NO₃:Fe ratio following artificial NO3- addition was 18.8 μmol:nmol, sufficiently high to cause Fe stress with phytoplankton growth following an increase in phytoplankton biomass.

Oceanic site O5 was located at Ocean Station Papa (OSP), station P26 of the Line-P transect (Figure 1). This site demonstrated characteristically high macronutrients and low dFe (0.05 nmol L⁻¹), resulting in the highest NO₃:Fe ratio observed across all experimental sites (234 μmol:nmol; Supplementary Table 1). Phytoplankton biomass was initially low, consistent with historical observations from this well-characterized Fe-limited region (Figure 2A; Supplementary Table 1; Boyd and Harrison, 1999). In contrast to most of the coastal sites, the majority of the phytoplankton biomass was dominated by picophytoplankton and other small cells (<5 μm) initially and throughout the incubation period (Supplementary Table 1; Figure 2B). Biogenic Si concentrations only increased after 96 h with similar responses in controls and Fe treatments (Figure 2B). Both large and small chl a size fractions, F_v/F_m, PN-specific VNO₃, and PC-specific VDIC were higher in the Fe treatment than in the unamended control (Ctl), confirming that the phytoplankton community in the initial seawater and in all incubation treatments without added Fe were experiencing Fe limitation (Figures 2B–E).

Metatranscriptomic assembly of sequence data and subsequent taxonomic annotation yielded the relative transcript proportions of phytoplankton functional groups (Figure 3). The CUZ site C1 was predominantly comprised of diatom transcripts at T₀; however, there was a 26% decrease in diatom transcripts in both the Fe and DFB treatments by T₁, accompanied by genus-level shifts within the diatoms. CUZ site C3 contained a phytoplankton community transcript pool almost equally represented by diatoms, prasinophytes, haptophytes, and dinoflagellates with little change in community composition among treatments following incubation. The coastal subarctic Pacific site C4 yielded an initial phytoplankton community transcript pool dominated by dinoflagellate-assigned sequences (24%), although these sequences decreased by ~10% in the Fe treatment, concurrent with a 9% increase in diatom transcripts. At the oceanic site O5, there were initially equal proportions of prasinophyte (22%) and haptophyte (23%) transcripts, with little representation by diatoms (4%). However, diatom-assigned transcripts constituted 9% of the community transcript pool by T₂ in the Fe addition treatment. Pseudo-nitzschia and Thalassiosira were among the top five diatom genera at all sites examined based on relative transcript abundance (Figure 3). These two genera together constituted between 9 and 53% of the transcript proportions in the initial diatom communities, and 25–58% of the Fe-enriched diatom communities.

Gene expression responses among sites were compared using Euclidian distance similarity analyses between Fe and DFB treatments (Fe/DFB, Fe/Ctl for O5) within the diatom genera Pseudo-nitzschia and Thalassiosira (Figure 4). Expression responses within coastal sites clustered together, while the oceanic site O5 displayed distinctly different patterns in both taxa. At site O5, 83 out of 1,334 KEGG Orthology genes (KOs) in Pseudo-nitzschia demonstrated >16-fold higher expression in the added Fe treatment than in the Fe-limited control treatment (Figure 4, Supplemental Table 4). By comparison, 155 out of 1,241 KOs in Thalassiosira showed >16-fold higher expression in the added Fe treatment compared to the low Fe control treatment. The most highly differentially expressed genes in oceanic Pseudo-nitzschia following Fe enrichment were ferritin (FTN, 290-fold), a metal transporter (CNNM, 32-fold), a putative bicarbonate (HCO₃) transporter (ICTB, 133-fold), and an NADPH-dependent glutamate synthase (GLT; 146-fold). In oceanic Thalassiosira, highly differentially expressed genes included ferredoxin-dependent sulfite reductase (Fd-SIR, 74-fold) and ferredoxin-dependent glutamate synthase (Fd-GLT; 416-fold). Fe addition induced both genera to increase the expression of several genes involved in photosynthesis by >16-fold exclusively at this location. Both taxa overexpressed gene products involved in vitamin biosynthesis, including the Fe-dependent vitamin B₇ synthesis protein biotin synthase (BIOB), which increased expression in the Fe enriched treatment by 84- and 49-fold in Pseudo-nitzschia and Thalassiosira, respectively. Furthermore, Pseudo-nitzschia increased expression of the vitamin B₁ (thiamine) biosynthetic gene THIC (by 179-fold) and vitamin B₆ (pyridoxine) biosynthetic genes pyridoxine kinase (PDXK; by 74-fold) and pyridoxine 4-dehydrogenase (PLDH; by 152-fold) following Fe enrichment at the oceanic site.

A number of genes demonstrated higher expression in the Fe-limited control treatment at O5. Forty-eight out of 1,334 genes in Pseudo-nitzschia and 77 out of 1,241 genes in Thalassiosira showed >16-fold higher expression in the Ctl treatment than in the added Fe treatment, patterns that were not found in diatoms from the examined coastal sites (Figure 4). In Thalassiosira, these genes encode proteins such as the copper (Cu)/zinc (Zn) superoxide dismutase (Cu-Zn SOD), an enzyme that removes toxic superoxide radicals by dismuting them into molecular oxygen and hydrogen peroxide, and a divalent metal transporter belonging to the ZIP family (ZIP7) (Marchetti and Maldonado, 2016). In both taxa, ribulose-1,5-bisphosphate carboxylase oxygenase (RubisCO; large subunit; RBCL), which catalyzes C-fixation in the Calvin cycle, had ≥24-fold higher expression in the Ctl treatment at O5.

The expression of genes involved in cellular growth and function, including N and C assimilation, vitamin synthesis, Fe-related metabolism, and trace metal acquisition, were compared in the dominant diatom genera Pseudo-nitzschia and Thalassiosira between the Fe and DFB/Ctl treatments (Figure 5). Genes encoding proteins involved in metal transport were detected at all locations, with expression patterns varying depending on site and taxa. Pseudo-nitzscha increased expression of the Fe transporter ABC.FEV.S by >2-fold under Fe enrichment at all locations where incubated communities showed a physiological Fe effect (C3, C4, O5; Supplementary Table 1). Transcripts for another Fe uptake protein, the high affinity iron permease FTR, were generally more abundant in the DFB/Ctl treatments in Thalassiosira, although the gene was more highly expressed following Fe enrichment in Pseudo-nitzschia at sites C2, C3, and O5 (Figure 5). The putative metal transporter CNNM was 32-fold more highly expressed following Fe enrichment in Pseudo-nitzschia at the oceanic site, but was not detected in Thalassiosira. Conversely, the non-specific metal transporter ZIP7 was 21-fold more highly expressed under Fe-limiting conditions in oceanic Pseudo-nitzschia and similarly not detected in oceanic Thalassiosira. Transcripts for Fe starvation induced proteins (ISIPs), including the recently-identified Fe acquisition protein ISIP2A that binds Fe at the cell surface and is thought to be involved in intracellular Fe transport (Morrissey et al., 2015), were highly abundant in Fe-stressed treatments (e.g., DFB and/or Ctl depending on the site) across all sites and in both taxa (Figure 5). Although their specific functions in diatoms are unclear, other ISIPs were markedly abundant and differentially expressed in the DFB/Ctl treatments, with ISIP1 one of the most differentially expressed genes between Fe-replete and Fe-limited treatments at each experimental site and in both taxa (Supplementary Figure 2).

Other Fe-related metabolic processes similarly varied depending on both site and taxa. Differences in expression patterns between taxa were generally greater for these Fe-related genes than in the N- and C-related genes investigated (Figure 5). At most sites, transcripts for the Fe storage protein ferritin (FTN) were higher in the Fe addition treatments than in the DFB/Ctl treatments. However, at two sites (C2 and C4), FTN transcripts were more abundant in the DFB treatment compared to the Fe addition treatment for one of the two genera (e.g., at site C2, 3.5-fold higher in Pseudo-nitzschia, p = 1 × 10⁻³ and at site C4, 90-fold higher in Thalassiosira). SODs were additionally differentially expressed, but they showed different expression patterns depending on the enzymes' metal cofactor(s) and the diatom genus. Cu-Zn SOD, which contains both Cu and Zn at its active site, showed a >100-fold higher expression in Thalassiosira in the Fe-limited control than in the added Fe treatment at the Fe-limited site O5. In contrast, in the same Fe-limited control treatment at this location, Pseudo-nitzschia demonstrated 2-fold higher expression of Fe-Mn SOD, which contains either Fe or manganese (Mn) as its metal cofactor. Based on the presence of Mn-coordinating amino acids at sites G-77 and Q-146 of the most highly expressed Fe-Mn SOD contigs, this Pseudo-nitzschia SOD was determined to specifically utilize Mn as its metal cofactor (Crowley et al., 2000; Groussman et al., 2015) (Supplementary Figure 4C).

Transcriptional responses of genes encoding Fe-dependent proteins and their functional replacements in photosynthetic electron transport were examined in both diatom genera (Figure 5). Transcripts for the Fe-independent protein flavodoxin (FLDA), which functionally replaces the Fe-protein ferredoxin (PETF) in photosynthetic electron transport, were generally more abundant in the DFB/Ctl treatments than in the Fe treatments in both genera (Figure 5). Conversely, transcripts of PETF were >2-fold higher in the high-Fe treatment only in Thalassiosira and across all sites. In Pseudo-nitzschia, PETF transcripts were either constitutively expressed (C3 and C4), more highly expressed in the DFB treatment (C1), or not present (C2 and O5) (Figure 5). Transcripts of cytochrome c₆ (PETJ) and its functional non-Fe replacement, the copper-protein plastocyanin, also showed differences in gene expression. PETJ transcripts were more abundant in the high Fe treatment at all sites and in both genera, except O5, where it was slightly more abundant in the Fe-limited control treatment (Figure 5). By contrast, transcripts for plastocyanin (PETE) displayed inconsistent expression trends in response to Fe status across sites, being relatively more abundant following Fe enrichment in both genera at C3 (1.4-fold in Pseudo-nitzschia; 1.9-fold in Thalassiosira, p = 5 × 10⁻⁴) and at the initially Fe-limited oceanic site, O5 (1.4-fold in Thalassiosira; Figure 5). At all other locations PETE transcripts were either more abundant under DFB conditions or not detected.

Transcripts for the proton-pumping protein rhodopsin (RHO) furthermore demonstrated differences in expression patterns among genera. This protein can supplement Fe-intensive photosynthesis in the light-driven production of membrane proton gradients and ATP in some diatoms (Marchetti et al., 2015). Rhodopsin was not detected in Thalassiosira at any location while its expression increased in Pseudo-nitzschia by >2-fold in the DFB/Ctl treatments relative to the Fe treatment at the two lowest dFe sites [C3 (p = 0.01) and O5; Figure 5; Supplementary Table 1]. At the other sites RHO expression was constitutive. These rhodopsin contigs were structurally similar to diatom rhodopsins identified within the MMETSP database (≥55% similarity; Supplementary Figure 4B).

Relationships among Fe-related transcript abundance, experimental site and treatment were determined using Principal Components Analysis (PCA) individually for each diatom genus. Principle components P1 and P2 explained 54% of the variation in transcript abundance in Pseudo-nitzschia and 76% in Thalassiosira (Figure 6C). In Pseudo-nitzschia, transcripts for the photosynthetic genes ferredoxin-NADP+ reductase (PETH), PETJ, a cytochrome b₆/f complex protein (PETC), FTN, and Cu-Zn SOD were in higher relative abundance within Fe addition treatments while RHO, ISIPs, FLDA, PETE, and FTR were generally more abundant in the Ctl and/or DFB treatments, as the principle component P1 separated these samples based on Fe treatment. In Thalassiosira, a similar response was observed, although RHO was not detected, and PETF, which was sporadically found and not abundant in Pseudo-nitzschia, strongly co-varied with the other genes highly expressed in the treatments where Fe was added (Figure 6C).

Genes involved in N transport and metabolism were investigated to assess the influence of varying Fe status on N assimilation. Transcripts for genes encoding nitrate (NRT2) and ammonium (AMT) transporters were detected at all locations, with NRT2 increasing in expression by >2-fold in response to Fe addition relative to the DFB/Ctl treatment at the majority of sites in both taxa, while AMT expression varied depending on site (Figure 5). For instance, C4 was the only location with a >2-fold increase in AMT expression in the DFB treatment in both Pseudo-nitzschia and Thalassiosira. Transcripts corresponding to genes encoding components of NO3- assimilation, including nitrate (NR) and nitrite reductases (NIRA, NIRB, NIT-6) were generally more abundant in the treatments with added Fe, although NIRA and NIRB displayed opposite expression patterns in Pseudo-nitzschia and Thalassiosira at site C3 (Figure 5). Furthermore Pseudo-nitzschia increased gene expression of one group of nitrite reductases [NIRB and NIT-6, which use NADPH as the reductant (Brown et al., 2009)] by 11- and 3.6-fold, respectively, following added Fe while Thalassiosira conversely increased NIRB expression by 3.7-fold in the DFB treatment (Figure 5). In addition, Thalassiosira increased gene expression of another form of nitrite reductase (NIRA, which uses ferredoxin/flavodoxin as reductant; Brown et al., 2009) by 8-fold (p = 3 × 10⁻²²) following Fe enrichment while Pseudo-nitzschia constitutively expressed NIRA at this location. Noticeably, transcripts for the genes encoding NIRB and NIT-6 were present in at least one of the two diatom taxa examined at all sites except the oceanic site, O5.

The relationships among transcript abundance for N uptake and assimilation-related genes, experimental sites, treatments and PN-specific VNO₃ measurements within the >5 μm size fraction of the phytoplankton community were examined via PCA bi-plots. Principle components P1 and P2 explained 86% of the variation in N-related transcript abundance in Pseudo-nitzschia and 88% in Thalassiosira (Figure 6A). Sites generally contained high transcript abundances of NRT2 and NR in the added Fe treatment, with the two genes strongly co-varying with one another in both Pseudo-nitzschia and Thalassiosira. Furthermore, the Fe addition treatments at two sites that experienced NO3- depletion following incubation, C2 and C3, clustered together and contained the highest AMT transcript abundance at T₁ and T₂, respectively. Phytoplankton communities within these incubation treatments concomitantly displayed low PN-specific VNO₃ (0.03–0.13 day⁻¹; Figure 6A). The highest PN-specific VNO₃ were observed in the added Fe treatment at site C4 at T₁ and at C1 within the initial (T₀) phytoplankton community (1.4 day⁻¹), which coincided with high abundances of NIRA transcripts in both genera at these locations.

To further gain insight into how variable Fe status influences macronutrient resource utilization and regional biogeochemistry, genes involved in C transport and fixation were examined among sites and between diatom genera. Transcripts corresponding to a carbonic anhydrase belonging to the α-family (α-CA), involved in the carbon concentrating mechanism (CCM) within photosynthetic eukaryotes (Reinfelder, 2011), were either constitutively expressed, not detected, or more highly expressed in the DFB treatment at all locations apart from C1, where expression was 7-fold higher following Fe addition in Thalassiosira (Figure 5).

Members of the solute carrier (SLC) family of bicarbonate transporters (SLC4A-1, -2, and -4), which import bicarbonate ions from the environment also thought to be involved in the CCM (Nakajima et al., 2013), were detected intermittently among sites, though in low transcript abundance (Figure 5). These genes share sequence homology with the P. tricornutum genes PtSLC4-1, -2, and -4 (BLASTP; E <2 × 10⁻⁶⁹) and displayed inconsistent patterns of gene expression with each another, with no clear relationship to carbon assimilation rates. Another putative bicarbonate transporter (ICTB) was detected intermittently across sites and solely in Pseudo-nitzschia, where it was notably more highly expressed by 128-fold following Fe addition at O5. Conversely in Thalassiosira, the gene encoding phosphoenolpyruvate carboxylase (PEPC), which is part of a C₄-CCM in some species of this genus (Reinfelder, 2011), was more highly expressed by 73-fold following Fe addition at O5.

Gene expression of RubisCO (RBCL) was higher by >24-fold in the Fe-limited control treatment in both genera at site O5 while at other sites the gene was either constitutively expressed, increased expression in the added Fe treatment, or not detected (Figure 5). In addition, other genes involved in the Calvin Cycle, including phosphoglycerate kinase (PGK), transketolase (TKL), ribulose-phosphate 3-epimerase (RPE), and phosphoribulokinase (PRK), generally increased in expression following Fe addition compared to the DFB/Ctl treatment at one or more of the three sites experiencing some degree of Fe limitation (C3, C4, and O5; Supplementary Table 1; Figure 5). At the CUZ sites C1 and C2, transcripts for these genes were either not differentially expressed or were more abundant in the DFB treatment within both diatom genera. Fructose-bisphosphate aldolases (FBA), involved in the Calvin Cycle, glycolysis, and gluconeogenesis, demonstrated strong Fe-dependent transcriptional patterns regardless of site and taxa (Figure 5). Transcripts corresponding to class II FBA, likely a metal-dependent aldolase, increased by 1.5 to 69-fold in the added Fe treatment as compared to DFB treatments with the largest fold change attributed to Pseudo-nitzschia from O5. Class II FBA has been previously demonstrated to be abundant under high-Fe conditions in diatoms and is hypothesized to contain Fe²⁺ as a metal cofactor (Horecker et al., 1972; Allen et al., 2012; Lommer et al., 2012). Transcripts corresponding to class I FBA, the metal-independent version of class II FBA, conversely increased by 1.3 to 16-fold in DFB compared to Fe treatments.

The relationships in transcript abundance among C fixation-related genes, experimental sites, incubated treatments and PC-specific VDIC measurements were assessed using PCA bi-plots. Principle components P1 and P2 together explained 80% of the variation in C-related transcript abundance in Pseudo-nitzschia and 78% in Thalassiosira (Figure 6B). Site C4 contained some of the highest PC-specific VDIC measurements within the >5 μm size fraction (0.65–1.6 day⁻¹), and coincided with the highest transcript abundances of PGK, PRK, FBP, TKL, RPE, and GAPDH in Pseudo-nitzschia (Figure 6B). Conversely, Fe-limited treatments from C3 and O5 had the lowest transcript abundances of these genes in both Pseudo-nitzschia and Thalassiosira, with principle component P1 separating these samples from other sites and treatments (Figure 6C). Fe-limited sites C3 and O5 phytoplankton communities additionally displayed some of the lowest PC-specific VDIC observed (0.11–0.17 day⁻¹).

Differences in gene expression patterns in response to Fe status were observed between the coastal (C1-C4) and oceanic sites (O5) examined in this study. This included the >16-fold higher expression of genes in the added Fe treatment relative to the Fe-limited control encoding proteins involved in B₇ synthesis (BIOB) in both taxa, and B₁ (THIC) and B₆ (PDXK, PDLH) synthesis in Pseudo-nitzschia. These increases are consistent with previous field observations demonstrating that Fe enrichment of previously Fe-limited oceanic diatom communities stimulates B-vitamin transcript production (Cohen et al., 2017). Genes encoding an Fe storage protein (ferritin [FTN]) and components of amino acid metabolism (glutamate synthase [GLT] in Pseudo-nitzschia; ferredoxin-dependent glutamate synthase [Fd-GLT] in Thalassiosira) were similarly more highly expressed by >16-fold following Fe addition exclusively at site O5. Conversely, in the Fe-limited control, we observed the >16-fold higher expression of genes encoding the proteins Cu-Zn superoxide dismutase (Cu-Zn SOD) and RubisCO (RBCL), in either one or both taxa investigated. These distinct transcriptomic patterns of genes involved in diverse metabolic processes reflect differences in environmental factors selecting for diatom growth between the chronically Fe-limited open ocean and sporadically Fe-limited coastal regions.

In contrast, many photosynthetic genes were highly expressed following Fe addition regardless of location. A subset of these genes displayed distinct expression responses depending on whether the incubated communities experienced Fe limitation of growth rate (e.g., C3 and O5) or only Fe stress (C4; Supplementary Table 1). One such gene encodes the putative Fe transporter ABC.FEV.S, in which expression increased following Fe addition in Pseudo-nitzschia only at sites C3, C4, and O5. Additional genes include flavodoxin (FLDA) and plastocyanin (PETE), in which transcripts were generally more abundant in the DFB or Fe-limited Ctl treatments, consistent with flavodoxin's role as an Fe-independent photosynthetic electron carrier and plastocyanin's role as a Cu-dependent replacement for Fe-dependent cytochrome c₆. At the Fe-stressed CUZ site C3 however, FLDA was either constitutively expressed or slightly more abundant after Fe addition, depending on the diatom genus. Plastocyanin (PETE) transcripts were similarly more abundant after Fe addition in both diatom genera at C3 and in Thalassiosira at O5. This pattern suggests coastal diatoms from higher-Fe systems tend to temporarily replace Fe-dependent photosynthetic proteins with Fe-independent ones, while certain diatoms in chronically Fe-limited environments may rely exclusively on the Fe-free alternatives (Marchetti et al., 2012).

Pseudo-nitzschia and Thalassiosira demonstrated several distinct responses to changes in Fe status despite co-existing under identical environmental conditions. Ferredoxin (PETF), ferredoxin-dependent glutamate synthase (Fd-GLT), and ferredoxin-dependent sulfite reductase (Fd-SIR) transcripts were more abundant in Thalassiosira at oceanic site O5 following Fe addition with these responses absent in Pseudo-nitzschia. In contrast, ferredoxin-related transcripts in oceanic Pseudo-nitzschia were constitutively expressed or not detected. These patterns may suggest oceanic Thalassiosira strongly utilizes ferredoxin and ferredoxin-dependent proteins following Fe addition while Pseudo-nitzchia relies on Fe-independent machinery. Site O5 was additionally the only location in which Thalassiosira increased gene expression of Cu-Zn SOD under Fe-limitation. This pattern was not evident in oceanic Pseudo-nitzschia, where gene expression of this protein was constitutive, or by either genus at coastal sites, suggesting that the oceanic Thalassiosira species have distinctly evolved to rely on this Cu- and Zn-containing enzyme as the preferred superoxide dismutase in their Fe-limited environment. Pseudo-nitzschia conversely increased expression of Mn SOD following Fe addition, likely as a result of iron-induced increases in photosynthetic rates and photosynthetic production of superoxide radicals (Asada, 2006). These patterns highlight differences in preferred metal cofactors as a function of Fe status and transcriptional tendencies between the two taxa.

Transcripts corresponding to rhodopsin (RHO) increased in abundance within Pseudo-nitzschia in the DFB/Ctl treatments at the two sites experiencing pronounced Fe limitation (C3 and O5), but were not identified in Thalassiosira at any location. This is consistent with rhodopsin being undetected in sequenced Thalassiosira spp. transcriptomes (Marchetti et al., 2015) and supports the notion that Pseudo-nitzschia may have a competitive advantage over non-rhodopsin containing taxa, allowing for an Fe-independent alternative to photosynthesis for ATP generation during times of Fe stress. Ferritin (FTN) gene expression patterns furthermore diverged between the two taxa at the coastal sites C4 (Line-P) and C2 (CUZ). This supports laboratory findings suggesting FTN may exhibit different expression patterns among diverse phytoplankton (Marchetti et al., 2009; Botebol et al., 2015), even between taxa residing in the same location. Lastly, ABC.FEV.S, encoding a membrane Fe transport system protein, displayed divergent expression patterns between the examined genera with only Pseudo-nitzschia increasing ABC.FEV.S expression after Fe addition in all incubations exhibiting signs of iron limitation (C3, C4, and O5).

Taken together, these patterns in gene expression demonstrate that members of the pennate diatom genus Pseudo-nitzschia and the centric diatom genus Thalassiosira restructure their functional metabolisms in response to changes in Fe availability in distinct manners, possibly allowing both species to co-exist in the same environment. Both taxa are equipped with strategies to sustain growth under chronic Fe limitation in the open ocean, as supported by their equal transcript abundance during initial sampling. Following pulse Fe additions however, oceanic Pseudo-nitzschia relies in part on the strategies discussed above to gain a competitive advantage over Thalassiosira and quickly dominates the phytoplankton community. It remains unclear however which combination of environmental factors in the NE Pacific Ocean would select for the preferential growth of Thalassiosira over Pseudo-nitzschia. We conclude that substantial differences in molecular responses to changes in Fe status are observed across taxonomic groups, and patterns in gene expression should not be assumed universal across diverse taxa or environments.

The majority of N transport and assimilation genes investigated increased in expression following Fe addition in both Pseudo-nitzschia and Thalassiosira. Several site- and taxa-specific patterns were identified, with some trends also possibly explained by each site's initial NO3- concentration. For example, most gene copies encoding the NO3- transporter, NRT2, have been demonstrated in laboratory cultures to increase in expression in NO₃-stressed diatoms (Bender et al., 2014; Rogato et al., 2015), and transcripts corresponding to this gene were some of the most abundant in both Pseudo-nitzschia and Thalassiosira at C2—the CUZ site where NO3- concentrations were depleted in all incubations by the first sampling time point. This gene also showed expression trends that correlated with Fe status; NRT2 transcripts were more abundant after Fe addition at all locations, regardless of initial NO3- concentrations. Based on these observations, NRT2 in diatoms also appears to be linked to Fe status and follows the expression of other N-related genes involved in Fe-dependent NO3- assimilation, including those encoding nitrate reductase (NR) and nitrite reductase (NIRA; Marchetti et al., 2012).

Diatoms were perhaps relying on NH₄ in place of NO3- as a source of N based on gene expression patterns at several CUZ sites. Fe-enriched treatments at C2 contained the lowest NO3- after 48 h of incubation (0.06 μmol L⁻¹), and the genes encoding the ammonium transporters AMTs concomitantly increased in expression in the Fe relative to DFB treatment (Figure 6). Furthermore at C3, Fe-enriched communities entered NO3- stress by the end of the incubation period, and AMT expression simultaneously increased in both Pseudo-nitzschia and Thalassiosira. This negative relationship between NO3- concentrations and AMT transcript abundance in natural diatom assemblages is consistent with those in laboratory Pseudo-nitzschia multiseries and Fragilariopsis cylindrus cultures (Bender et al., 2014; Rogato et al., 2015), and is reported here as one of the first observations of this relationship in natural phytoplankton communities.

High AMT transcript abundance at some of these locations may also represent NH₄ rather than NO3- being preferred as an N source by Fe-stressed diatoms conserving their cellular Fe supply, as NO3- assimilation depends on various Fe-dependent processes (Milligan and Harrison, 2000). This is supported by the increased expression of AMT transcripts in both Pseudo-nitzschia and Thalassiosira from the Fe-stressed coastal Line-P incubations at C4. Pseudo-nitzschia from the Fe-limited site O5 also exhibited this pattern whereas Pseudo-nitzschia from C3 and Thalassiosira from both C3 and O5 did not, suggesting other environmental parameters aside from Fe status are influencing whether diatoms utilize NH₄- or NO₃-specific N uptake pathways.

Similar to our Fe-related gene expression results, several N-related genes demonstrated divergent expression responses between Pseudo-nitzschia and Thalassiosira. Expression of the NO2- reductase genes, NIRA and NIRB, displayed opposite patterns between the two genera at the CUZ site where Fe-stress occurred in incubations (C3), with Pseudo-nitzschia highly expressing the gene encoding non-ferredoxin-utilizing NO2- reductase (NIRB) following Fe addition, and Thalassiosira highly expressing the gene encoding the ferredoxin-utilizing nitrite reductase (NIRA). Furthermore at site O5, Pseudo-nitzschia increased expression of AMT and NADPH-dependent glutamate synthase (GLT) following Fe addition while Thalassiosira increased expression of NRT2 and ferredoxin-dependent glutamate synthase (Fd-GLT). These transcriptomic patterns may suggest Pseudo-nitzschia continues to rely on the non-Fe requiring metabolic pathways for assimilating N once Fe becomes available (AMT, NIRB, GLT), whereas Thalassiosira shifts over to Fe-dependent ones (NRT2, NIRA, Fd-GLT) upon Fe resupply.

These expression patterns furthermore support that substantial variations exist between the two diatom taxa in terms of N acquisition and assimilation strategies following changes in Fe supply. Both Pseudo-nitzschia and Thalassiosira are equipped with distinct strategies to compete under a variety of Fe and N conditions, and this may contribute to how multiple diatom species relying upon the same limiting resources in identical environments co-exist (i.e., paradox of the plankton; Hutchinson, 1961). These patterns are consistent with previous reports of resource partitioning among diatoms based on N and phosphate utilization (Alexander et al., 2015). Varying environmental pressure likely maintain populations of diverse diatom genera in the open ocean, with certain species outcompeting others depending on specific sets of external factors, including both macro- and micronutrients (Godhe and Rynearson, 2017).

Genes encoding proteins involved in C uptake and assimilation were surveyed in order to determine the influence of Fe addition or stress on C metabolism. We observed site-specific expression patterns of the diatom RubisCO large subunit protein (RBCL), where gene expression was substantially elevated at site O5 in the Fe-limited control treatment relative to the Fe addition response in both diatom genera. A sequence analysis of RubisCO contigs obtained across experimental sites demonstrates that O5 protein sequences are structurally less similar to known Pseudo-nitzschia and Thalassiosira RubisCO protein sequences within the MMETSP database than those at the four coastal sites (Supplementary Figure 4A; Supplementary Table 3). This distinction in both protein structure and transcriptional expression may indicate a distinct adaptation and utilization of RubisCO in the oceanic diatoms than in those from high-Fe coastal waters. Phylogenetically diverse diatom species have been demonstrated to vary in their RubisCO enzyme kinetics in laboratory cultures, with their RubisCO content inversely linked to the strength of their carbon concentrating mechanism (CCM; Young et al., 2016). The CCM increases CO₂ concentrations in chloroplast stroma in the vicinity of RubisCO and is fueled by the energy (ATP) generated from the Fe-intensive process of photosynthesis (Reinfelder, 2011; Young et al., 2016). We hypothesize that chronically Fe-limited oceanic diatoms are ATP-limited by the scarcity of Fe needed to support photosynthesis, and instead increase their RubisCO protein content to maintain high rates of carbon fixation rather than allocate scarce energy resources to the CCM. Furthermore, the genes encoding a putative bicarbonate transporter (ICBT) and a C₄-CCM component (PEPC; Reinfelder et al., 2000; Sage, 2004; Reinfelder, 2011) were highly expressed following Fe addition in Pseudo-nitzschia and Thalassiosira, respectively, exclusively at O5. This supports that diatoms may be capable of shuffling energy pools into either the CCM or RubisCO production depending on Fe bioavailability. Interestingly, in laboratory-based proteomic analyses with cultures of the coastal diatom T. pseudonana, RubisCO was similarly more highly expressed under Fe limitation, while PEPC protein levels were higher under Fe-replete conditions (Nunn et al., 2013). Consistent with our hypothesis, Hopkinson et al. (2010) attributed increases in biomass following CO₂-enrichment of an Fe-limited phytoplankton community in the HNLC Northeast Pacific Ocean to downregulation of the CCM in order to conserve iron and photosynthetically-produced energy. Laboratory-based RubisCO kinetic work with cultured diatom isolates is needed to confirm whether diatoms from HNLC regions minimize their photosynthetic demand for Fe by synthesizing more RubisCO enzymes rather than allocating scarce energy resources into the CCM.

Other C fixation-related gene expression patterns were largely consistent with C assimilation rates, and generally varied as a function of both Fe status and ocean province. The genes PGK, TKL, RPE, and PRK did not exhibit site-specific expression patterns similar to RBCL, and instead increased in expression following Fe enrichment at sites where Fe addition increased C assimilation rates (C3, C4, and O5). Increased expression of these genes is expected with Fe stimulation of C-fixation and growth. These expression patterns are in agreement with laboratory cultures of the diatom P. tricornutum, which increased expression of genes involved in C fixation during the light portion of their diel cycle, when DIC is being taken up to support photosynthesis (Chauton et al., 2013).

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.