Table of Contents

• Coverage
• Dataset Description
  • Methods & Sampling
  • Data Processing Description
  • BCO-DMO Curation Notes
• Data Files
• Related Publications
• Related Datasets
• Parameters
• Instruments
• Deployments
• Project Information
• Funding

This dataset is one of many generated from an iron incubation experiment conducted as part of the PUPCYCLE II cruise in May and June 2023. The different analyses are listed here, with links to other datasets from this study in the Related Datasets section.

Analyses include:

• Physiology (this dataset 994607)
  • Dissolved inorganic nutrients (nitrate, phosphate, silicic acid)
  • Size-fractionated chlorophyll a
  • Size-fractionated isotope uptake rates (nitrate and DIC)
  • Flow cytometry
  • FlowCAM analysis
  • FLP experiments for mixotroph cell-specific grazing rates
• RNA taxonomy data for normalized reads (dataset 1003525)
• Differential expression for +Fe versus DFB treatments at T1 and T2: Pfam
  • Mixotroph Pfam (dataset 1004659)
  • Diatom Pfam (dataset 1004674)
• Differential expression for +Fe versus DFB treatments at T1 and T2 for mixotrophs and diatoms: KEGG gene (dataset 1004689)
• Differential expression for +Fe versus DFB treatments at T1 and T2 for mixotrophs and diatoms: KEGG pathways (this dataset 1004704)
• Dissolved iron concentrations (dataset 1004007)

See Related Datasets section below for links to above mentioned datasets.

Data collection took place on the R/V Sally Ride from May 29th to June 10th, 2023. To simulate upwelling conditions under different iron treatments, an onboard incubation experiment was conducted. Seawater for the incubation experiment was collected within the northern CCS, off the southern coast of Oregon at 43°02'42.7"N, 124°33'07.2"W. The collected seawater was deemed as freshly upwelled water and taken from a depth of 55 meters – which corresponds to a depth slightly below the euphotic zone receiving less than 1% irradiance – using trace-metal clean techniques on May 30th, 2023, 13:30 GMT. The seawater was pumped and homogenized using trace metal clean techniques, then transferred into a total of thirty 20 L low-density polyethylene cubitainers. Three cubitainers were immediately harvested for the initial timepoint (T0). The remaining twenty-seven cubitainers were assigned treatments: nine were unamended (labeled Ctrl), nine were amended with 5 nM FeCl₂ (labeled +Fe), and nine were amended with 200 nM desferrioxamine B, a strong iron chelator that inhibits dissolved iron uptake (labeled DFB). The cubitainers were placed in an on-deck incubator covered with two layers of neutral-density screening to simulate 26% of incident irradiance and supplied with flow-through surface seawater to maintain ambient surface temperature. Three cubitainers from each treatment were harvested for each of the three subsequent timepoints: 48 hours (T1), 168 hours (T2), and 264 hours (T3) after incubation.

Dissolved inorganic nutrients. Dissolved inorganic nutrient measurements were collected by filtering 30 mL of seawater through a GF/F filter into an acid-rinsed polypropylene scintillation vial. The filtrate was immediately frozen and then sent to Wetland Biogeochemistry Analytic Services at Louisiana State University and analyzed via their OI Analytical Flow Solutions IV auto analyzer for nitrate + nitrite, nitrite, phosphate (PO₄^3-), and silicic acid (Si(OH)₄). Detection limits were 0.09 μM for nitrate + nitrite, 0.02 μM for PO₄^3-, and 0.02 μM for Si(OH)₄. Nitrate (NO₃^-) concentrations were derived by subtracting the concentration of nitrite from nitrate + nitrite. NO₃^- concentrations were also used in calculations for NO₃^- uptake rates.

Size-fractionated chlorophyll a. Chlorophyll a (chl a) measurements were collected by gravity filtering 250 mL of seawater through 5 μm Isopore membrane filters (47 mm) to collect large cells (≥ 5 μm), then vacuum filtering the remaining seawater onto GF/F (25 mm) filters under 100 mmHg of vacuum pressure to collect small cells (< 5 μm). Filters were then rinsed with filtered seawater and stored at -20℃. Chl a was extracted on the ship using 90% acetone solution at -20℃ for 24 hours, then measured on a 10-AU fluorometer (Turner Designs, San Jose, CA) using the acidification method (Parsons et al., 1984). 

Size-fractionated isotope uptake rates. Isotope uptake rates of dissolved inorganic carbon (DIC) and NO₃^- were assessed by collecting seawater into 1 L polycarbonate bottles (filled to the top, resulting in a total volume of 1.09 L). Each bottle was immediately spiked with NaH¹³CO₃ and Na¹⁵NO₃, with spiking concentrations varying depending on the ambient DIC and NO₃^- concentrations. The ambient DIC concentrations were estimated at 2000 μM based on literature values (Fassbender et al., 2011), while the NO₃^- ambient concentrations were estimated using the Submersible Ultraviolet Nitrate Analyzer (SUNA) (Johnson & Coletti, 2002). These estimates resulted in adding 200 μM of ¹³C to every sample as well as 2, 2, 0.50, and 0.20 μM of ¹⁵N at T0, T1, T2, and T3 respectively, to achieve approximate additions of 10% of the ambient DIC and NO₃^- concentrations.

Samples were then placed in the ship incubator for 6 hours, simulated under 25% of incident irradiance. After incubation, the seawater from each bottle was gravity filtered onto 5 μm Isopore membrane filters (47 mm) to collect large cells and vacuum filtered onto pre-combusted (450℃ for 5 hours) GF/F filters to collect small cells. Large cells were then rinsed off the Isopore filter using 0.2 mm filtered seawater and vacuum filtered onto pre-combusted (450℃ for 5 hours) GF/F filters. The standard filtering volume was 1.09 L, but due to accumulation of biomass hindering filtering at the later timepoints, all samples in T2 and T3 had a filtering volume of 0.5 L. Filters were stored in small petri dishes at -20℃ until lab preparation.

To prepare for lab analysis, filters were dried at 60℃ for 24 hours, encapsulated in tin, and pelletized. Pelletized samples were sent to the UC Davis Stable Isotope Facility, where they were analyzed for particulate organic carbon (POC), particulate organic nitrogen (PON), and atom percentages of ¹³C and ¹⁵N using an isotope ratio mass spectrometer (EA-IRMS). POC and PON concentrations (μM) were calculated by dividing the mass of POC/PON (μg) by the respective atomic mass of carbon and nitrogen over the volume filtered. Absolute uptake rates (⍴) of DIC and NO₃^- were calculated using the constant transport model and biomass-normalized uptake rates (V) of DIC and NO₃^- were calculated using the specific uptake model; these calculations utilize actual dissolved NO₃^- concentrations (Dugdale and Wilkerson, 1986).

Flow Cytometry. Flow cytometry was performed on a Guava easyCyte HT flow cytometer live on the cruise. The Guava easyCyte HT was outfitted with a blue 488 nm laser. Gain settings were forward scatter (FSC) 4, side scatter (SSC) 1, green fluorescence 2.95, yellow fluorescence 98.7, and red fluorescence 4. Communities of interest were gated based on community recommendations (Thyssen et al., 2022). Synechococcus were detected based on high yellow fluorescence and low FSC, picoeukaryotes (~1-3 µm) were detected based on medium yellow and red fluorescence, and nanoeukaryotes (~3-20 µm) were detected based on high red fluorescence. Samples for bacteria abundance were preserved with 0.25% glutaraldehyde and frozen at -80 ºC. Samples were then stained SYBR Green I (Thermo Fisher; 10,000x stock) diluted to a final concentration of 1x and measured on the Guava EasyCyte HT flow cytometer.

FlowCAM Analysis. FlowCAM imaging was performed for T2 samples of the iron incubation experiment based on the protocol provided by Pierce et al. 2023 using a Portable Series FlowCam connected to a peristaltic pump. A 10x FlowCell objective was used on autoimage mode for samples preserved in Lugols. The sets of images produced were analyzed using the Yokogawa Fluid Imaging VisualSpreadsheet software. Any images not containing particles (such as bubbles or blurry images) and/or less than 15 µm were deleted from the dataset. For each sample, the first 300 identifiable diatom cells were counted and categorized based on cell morphology and assigned to taxa at the genus level when possible. Final categories included: Asterionellopsis spp., Chaetoceros spp., Pseudo-nitzschia spp., other centric diatoms, and other pennate diatoms. To calculate concentrations of each diatom taxonomy within each sample (cells/L), the following equation was used:

(# of cells counted * image factor) / volume run through FlowCAM

The image factor was determined by taking the total number of images for the sample (after quality control) divided the number of images counted for all taxonomy. The volume processed through the FlowCAM was 25 mL for all samples.

FLP Experiments for Mixotroph Grazing. Fluorescently labeled particle (FLP) incubations were performed and analyzed for the iron incubation experiment based on the protocol provided and used by Cook et al. (2025). 500 mL of seawater was pre-screened through 200 µm mesh to remove mesozooplankton. Fluorescently labeled surrogate microspheres (0.5 µm diameter, Fluoresbrite, Polysciences) were spiked into the volume at concentrations of 10^⁵ particles mL^-1. Immediately after surrogate prey addition, a time zero (T_initial) sample was taken and preserved. For preservation, 10-30 mL were fixed with 4% glutaraldehyde, incubated at 4 ºC for 30 minutes, filtered onto 3 µm polycarbonate filters, mounted using DAPI Vectashield (Vector Laboratories), and frozen at -20 ºC. Incubations were carried out for one hour in a clear plexiglass flow-through incubator with mesh bags emulating in situ light availability. After an hour incubation, another sample (T_final) was taken and preserved in the same manner as T_initial. 

FLP incubation analysis was performed using an Olympus CKX53 inverted microscope at 400x total magnification. The field-of-view counting method with a calibrated reticle was used. At least 300 chloroplast-containing cells were enumerated on every slide to ensure statistical reliability. Cell specific grazing rate (CSGR, bacteria mixotroph^-1 hr^-1), was calculated using the following equation:

CSGR = (((C_ingested / C_photo)_Tfinal - (C_ingested / C_photo)_Tinitial) / t) * (C_bacteria / C_FLP)

C_ingested is the count of ingested beads divided by the count of all pigmented phototrophic cells (C_photo), with T_initial values subtracted from T_final values to account for non-ingestion adsorption. Values were divided by time of incubation in hours (t) then multiplied by the Bacteria:FLP ratio. C_bacteria is the natural concentration of heterotrophic bacteria measured via flow cytometry and C_FLP is the concentration of FLPs spiked into the sample (10⁵ particles mL^-1). 

The Yokogawa Fluid Imaging VisualSpreadsheet software was used to analyze FlowCAM samples. Excel and R v4.3.2 were used to input sample information and calculations. 

CURATION ACTIONS PERFORMED ON DATA
- Imported source file "BCO_DMO_PUPCYCLEII_Physiology.csv" with NA as missing data identifier replaced with empty cells.
- Converted the Datetime column into ISO 8601 DateTime format
- Added columns for latitude and longitude of sampling station
- Set data types as string for the first 5 columns, set time point hours as integer, date as datetime, and all others as numeric
- Exported the final file as "994607_v1_physiology_pupcycle2.csv"

CURATION ACTIONS PERFORMED ON METADATA
- Added summary for the broader study showing connections between this dataset and related datasets
- Changed in-line citations to have format of (author, year) replacing number references.

ISSUES POTENTIALLY IMPACTING REUSE
- N/A

NSF Award Abstract:
Upwelling zones are hotspots of photosynthesis that are very dynamic in space and time. Microscopic algae, known as phytoplankton, bloom when deep, nutrient-rich waters are upwelled into sunlit surface layers of the ocean, providing nourishment that supports productive food webs and draws down carbon dioxide (CO2) from the atmosphere to the deep ocean. Photosynthetic microbes in these regions must constantly adapt to changes in their chemical and physical environments. For example, subsurface populations respond to changes in light as they approach the surface. When upwelled waters move offshore, cells sink out of the illuminated zone, establishing seed populations that remain inactive until the next upwelling event. This process is called the upwelling conveyor belt cycle (UCBC). How phytoplankton respond to these changes in environmental conditions and how they may influence their nutrient requirements remains unknown. With future ocean changes predicted to alter seawater chemistry, including ocean acidification and decreased iron availability, some phytoplankton groups may be more vulnerable than others. Accompanying educational activities provide learning experiences to enhance understanding and awareness of marine microbes. The development of a research hub at UNC aims to provide infrastructure and support for scientists and students conducting research on environmental genomics. A laboratory component for an upper-level undergraduate course focused on marine phytoplankton is being developed. Educational outreach activities to broader communities include creation of a lesson plan on phytoplankton in upwelling zones and a virtual research cruise experience for middle-school students, as well as a hands-on lab activity for a local museum focused on marine phytoplankton and the important roles they play in shaping our planet.

The project examines how phytoplankton respond at the molecular and physiological level to the different UCBC stages, which seed populations (i.e., surface versus subsurface) contribute most to phytoplankton blooms during upwelling events of varying intensity, how phytoplankton elemental compositions are altered throughout UCBC stages, and how future predicted ocean conditions will affect the phytoplankton responses to UCBC conditions. This project contains both laboratory and fieldwork. In the laboratory, phytoplankton isolates recently obtained from upwelling regions are exposed to simulated UCBC conditions to examine changes in gene expression, growth and photosynthetic characteristics and elemental composition. Cultures are subjected to both current and future ocean conditions, including reduced iron availability and higher CO2. In the field, research cruises within upwelling regions study the dynamics of natural phytoplankton communities (both surface and subsurface) experiencing upwelling and relaxation and within simulated upwelling incubation experiments. Knowledge of how phytoplankton are affected by UCBC conditions at an integrated molecular, physiological and elemental level under both current and future scenarios is imperative for the proper conservation and management of these critically important ecosystems.