A discrete, stochastic model of colonial phytoplankton population size structure: Development and application to in situ imaging-in-flow cytometer observations of Dinobryon

The scientific community lacks models for the dynamic changes in population size structure that occur in colonial phytoplankton. This is surprising, as size is a key trait affecting many aspects of phytoplankton ecology, and colonial forms are very common. We aim to fill this gap with a new discrete, stochastic model of dynamic changes in phytoplankton colonies' population size structure. We use the colonial phytoplankton Dinobryon as a proof-of-concept organism. The model includes four stochastic functions-division, stomatocyst production, colony breakage, and colony loss-to determine Dinobryon population size structure and populations counts. Although the functions presented here are tailored to Dinobryon, the model is readily adaptable to represent other colonial taxa. We demonstrate how fitting our model to in situ observations of colony population size structure can provide a powerful approach to explore colony size dynamics. Here, we have (1) collected high-frequency in situ observations of Dinobryon in Lac (Lake) Montjoie (Quebec, Canada) in 2013 with a moored Imaging FlowCytobot (IFCB) and (2) fit the model to those observations with a genetic algorithm solver that extracts parameter estimates for each of the four stochastic functions. As an example of the power of this model-data integration, we also highlight ecological insights into Dinobryon colony size and stomatocyst production. The Dinobryon population was enriched in larger, flagellate-rich colonies near bloom initiation and shifted to smaller and emptier colonies toward bloom decline.

Dinoflagellate vertical migration fuels an intense red tide

Harmful algal blooms (HABs) are increasing globally, causing economic, human health, and ecosystem harm. In spite of the frequent occurrence of HABs, the mechanisms responsible for their exceptionally high biomass remain imperfectly understood. A 50-y-old hypothesis posits that some dense blooms derive from dinoflagellate motility: organisms swim upward during the day to photosynthesize and downward at night to access deep nutrients. This allows dinoflagellates to outgrow their nonmotile competitors. We tested this hypothesis with in situ data from an autonomous, ocean-wave-powered vertical profiling system. We showed that the dinoflagellate Lingulodinium polyedra's vertical migration led to depletion of deep nitrate during a 2020 red tide HAB event. Downward migration began at dusk, with the maximum migration depth determined by local nitrate concentrations. Losses of nitrate at depth were balanced by proportional increases in phytoplankton chlorophyll concentrations and suspended particle load, conclusively linking vertical migration to the access and assimilation of deep nitrate in the ocean environment. Vertical migration during the red tide created anomalous biogeochemical conditions compared to 70 y of climatological data, demonstrating the capacity of these events to temporarily reshape the coastal ocean's ecosystem and biogeochemistry. Advances in the understanding of the physiological, behavioral, and metabolic dynamics of HAB-forming organisms from cutting-edge observational techniques will improve our ability to forecast HABs and mitigate their consequences in the future.

Temperature dependence of parasitoid infection and abundance of a diatom revealed by automated imaging and classification

Diatoms are a group of phytoplankton that contribute disproportionately to global primary production. Traditional paradigms that suggest diatoms are consumed primarily by larger zooplankton are challenged by sporadic parasitic "epidemics" within diatom populations. However, our understanding of diatom parasitism is limited by difficulties in quantifying these interactions. Here, we observe the dynamics of Cryothecomonas aestivalis (a protist) infection of an important diatom on the Northeast U.S. Shelf (NES), Guinardia delicatula, with a combination of automated imaging-in-flow cytometry and a convolutional neural network image classifier. Application of the classifier to >1 billion images from a nearshore time series and >20 survey cruises across the broader NES reveals the spatiotemporal gradients and temperature dependence of G. delicatula abundance and infection dynamics. Suppression of parasitoid infection at temperatures <4 °C drives annual cycles in both G. delicatula infection and abundance, with an annual maximum in infection observed in the fall-winter preceding an annual maximum in host abundance in the winter-spring. This annual cycle likely varies spatially across the NES in response to variable annual cycles in water temperature. We show that infection remains suppressed for ~2 mo following cold periods, possibly due to temperature-induced local extinctions of the C. aestivalis strain(s) that infect G. delicatula. These findings have implications for predicting impacts of a warming NES surface ocean on G. delicatula abundance and infection dynamics and demonstrate the potential of automated plankton imaging and classification to quantify phytoplankton parasitism in nature across unprecedented spatiotemporal scales.

Predator switching strength controls stability in diamond-shaped food web models

In food web models that include more than one prey type for a single predator, it is common for the predator's functional response to include some form of switching-preferential consumption of more abundant prey types. Predator switching promotes coexistence among competing prey types and increases diversity in the prey community. Here, we show how the dynamics of a diamond-shaped food web model of a marine plankton community are sensitive to a parameter that sets the strength of predator switching. Stronger switching destabilizes the model's coexistence equilibrium and leads to the appearance of limit cycles. Stronger switching also increases the evenness of the asymptotic prey community and promotes synchrony in the dynamics of disparate prey types. Given the dependence of model behavior on the strength of predator switching, it is important that modelers carefully consider the parameterization of functional responses that include switching.

Marine Pelagic Ecosystem Responses to Climate Variability and Change

The marine coastal region makes up just 10% of the total area of the global ocean but contributes nearly 20% of its total primary production and over 80% of fisheries landings. Unicellular phytoplankton dominate primary production. Climate variability has had impacts on various marine ecosystems, but most sites are just approaching the age at which ecological responses to longer term, unidirectional climate trends might be distinguished. All five marine pelagic sites in the US Long Term Ecological Research (LTER) network are experiencing warming trends in surface air temperature. The marine physical system is responding at all sites with increasing mixed layer temperatures and decreasing depth and with declining sea ice cover at the two polar sites. Their ecological responses are more varied. Some sites show multiple population or ecosystem changes, whereas, at others, changes have not been detected, either because more time is needed or because they are not being measured.

Machine learning techniques to characterize functional traits of plankton from image data

Plankton imaging systems supported by automated classification and analysis have improved ecologists' ability to observe aquatic ecosystems. Today, we are on the cusp of reliably tracking plankton populations with a suite of lab-based and in situ tools, collecting imaging data at unprecedentedly fine spatial and temporal scales. But these data have potential well beyond examining the abundances of different taxa; the individual images themselves contain a wealth of information on functional traits. Here, we outline traits that could be measured from image data, suggest machine learning and computer vision approaches to extract functional trait information from the images, and discuss promising avenues for novel studies. The approaches we discuss are data agnostic and are broadly applicable to imagery of other aquatic or terrestrial organisms.

Seasonal environmental variability drives microdiversity within a coastal Synechococcus population

Marine microbes often show a high degree of physiological or ecological diversity below the species level. This microdiversity raises questions about the processes that drive diversification and permit coexistence of diverse yet closely related marine microbes, especially given the theoretical efficiency of competitive exclusion. Here, we provide insight with an 8‐year time series of diversity within Synechococcus, a widespread and important marine picophytoplankter. The population of Synechococcus on the Northeast U.S. Shelf is comprised of six main types, each of which displays a distinct and consistent seasonal pattern. With compositional data analysis, we show that these patterns can be reproduced with a simple model that couples differential responses to temperature and light with the seasonal cycle of the physical environment. These observations support the hypothesis that temporal variability in environmental factors can maintain microdiversity in marine microbial populations. We also identify how seasonal diversity patterns directly determine overarching Synechococcus population abundance features.

Dynamics and functional diversity of the smallest phytoplankton on the Northeast US Shelf

Picophytoplankton are the most abundant primary producers in the ocean. Knowledge of their community dynamics is key to understanding their role in marine food webs and global biogeochemical cycles. To this end, we analyzed a 16-y time series of observations of a phytoplankton community at a nearshore site on the Northeast US Shelf. We used a size-structured population model to estimate in situ division rates for the picoeukaryote assemblage and compared the dynamics with those of the picocyanobacteria Synechococcus at the same location. We found that the picoeukaryotes divide at roughly twice the rate of the more abundant Synechococcus and are subject to greater loss rates (likely from viral lysis and zooplankton grazing). We describe the dynamics of these groups across short and long timescales and conclude that, despite their taxonomic differences, their populations respond similarly to changes in the biotic and abiotic environment. Both groups appear to be temperature limited in the spring and light limited in the fall and to experience greater mortality during the day than at night. Compared with Synechococcus, the picoeukaryotes are subject to greater top-down control and contribute more to the region's primary productivity than their standing stocks suggest.

Seasons of Syn

Synechococcus is a widespread and important marine primary producer. Time series provide critical information for identifying and understanding the factors that determine abundance patterns. Here, we present the results of analysis of a 16-yr hourly time series of Synechococcus at the Martha's Vineyard Coastal Observatory, obtained with an automated, in situ flow cytometer. We focus on understanding seasonal abundance patterns by examining relationships between cell division rate, loss rate, cellular properties (e.g., cell volume, phycoerythrin fluorescence), and environmental variables (e.g., temperature, light). We find that the drivers of cell division vary with season; cells are temperature-limited in winter and spring, but light-limited in the fall. Losses to the population also vary with season. Our results lead to testable hypotheses about Synechococcus ecophysiology and a working framework for understanding the seasonal controls of Synechococcus cell abundance in a temperate coastal system.

An Ocean-Colour Time Series for Use in Climate Studies: The Experience of the Ocean-Colour Climate Change Initiative (OC-CCI)

Ocean colour is recognised as an Essential Climate Variable (ECV) by the Global Climate Observing System (GCOS); and spectrally-resolved water-leaving radiances (or remote-sensing reflectances) in the visible domain, and chlorophyll-a concentration are identified as required ECV products. Time series of the products at the global scale and at high spatial resolution, derived from ocean-colour data, are key to studying the dynamics of phytoplankton at seasonal and inter-annual scales; their role in marine biogeochemistry; the global carbon cycle; the modulation of how phytoplankton distribute solar-induced heat in the upper layers of the ocean; and the response of the marine ecosystem to climate variability and change. However, generating a long time series of these products from ocean-colour data is not a trivial task: algorithms that are best suited for climate studies have to be selected from a number that are available for atmospheric correction of the satellite signal and for retrieval of chlorophyll-a concentration; since satellites have a finite life span, data from multiple sensors have to be merged to create a single time series, and any uncorrected inter-sensor biases could introduce artefacts in the series, e.g., different sensors monitor radiances at different wavebands such that producing a consistent time series of reflectances is not straightforward. Another requirement is that the products have to be validated against in situ observations. Furthermore, the uncertainties in the products have to be quantified, ideally on a pixel-by-pixel basis, to facilitate applications and interpretations that are consistent with the quality of the data. This paper outlines an approach that was adopted for generating an ocean-colour time series for climate studies, using data from the MERIS (MEdium spectral Resolution Imaging Spectrometer) sensor of the European Space Agency; the SeaWiFS (Sea-viewing Wide-Field-of-view Sensor) and MODIS-Aqua (Moderate-resolution Imaging Spectroradiometer-Aqua) sensors from the National Aeronautics and Space Administration (USA); and VIIRS (Visible and Infrared Imaging Radiometer Suite) from the National Oceanic and Atmospheric Administration (USA). The time series now covers the period from late 1997 to end of 2018. To ensure that the products meet, as well as possible, the requirements of the user community, marine-ecosystem modellers, and remote-sensing scientists were consulted at the outset on their immediate and longer-term requirements as well as on their expectations of ocean-colour data for use in climate research. Taking the user requirements into account, a series of objective criteria were established, against which available algorithms for processing ocean-colour data were evaluated and ranked. The algorithms that performed best with respect to the climate user requirements were selected to process data from the satellite sensors. Remote-sensing reflectance data from MODIS-Aqua, MERIS, and VIIRS were band-shifted to match the wavebands of SeaWiFS. Overlapping data were used to correct for mean biases between sensors at every pixel. The remote-sensing reflectance data derived from the sensors were merged, and the selected in-water algorithm was applied to the merged data to generate maps of chlorophyll concentration, inherent optical properties at SeaWiFS wavelengths, and the diffuse attenuation coefficient at 490 nm. The merged products were validated against in situ observations. The uncertainties established on the basis of comparisons with in situ data were combined with an optical classification of the remote-sensing reflectance data using a fuzzy-logic approach, and were used to generate uncertainties (root mean square difference and bias) for each product at each pixel.

Resonance control of acoustic focusing systems through an environmental reference table and impedance spectroscopy

Acoustic standing waves can precisely focus flowing particles or cells into tightly positioned streams for interrogation or downstream separations. The efficiency of an acoustic standing wave device is dependent upon operating at a resonance frequency. Small changes in a system's temperature and sample salinity can shift the device's resonance condition, leading to poor focusing. Practical implementation of an acoustic standing wave system requires an automated resonance control system to adjust the standing wave frequency in response to environmental changes. Here we have developed a rigorous approach for quantifying the optimal acoustic focusing frequency at any given environmental condition. We have demonstrated our approach across a wide range of temperature and salinity conditions to provide a robust characterization of how the optimal acoustic focusing resonance frequency shifts across these conditions. To generalize these results, two microfluidic bulk acoustic standing wave systems (a steel capillary and an etched silicon wafer) were examined. Models of these temperature and salinity effects suggest that it is the speed of sound within the liquid sample that dominates the resonance frequency shift. Using these results, a simple reference table can be generated to predict the optimal resonance condition as a function of temperature and salinity. Additionally, we show that there is a local impedance minimum associated with the optimal system resonance. The integration of the environmental results for coarse frequency tuning followed by a local impedance characterization for fine frequency adjustments, yields a highly accurate method of resonance control. Such an approach works across a wide range of environmental conditions, is easy to automate, and could have a significant impact across a wide range of microfluidic acoustic standing wave systems.

Satellite sensor requirements for monitoring essential biodiversity variables of coastal ecosystems

The biodiversity and high productivity of coastal terrestrial and aquatic habitats are the foundation for important benefits to human societies around the world. These globally distributed habitats need frequent and broad systematic assessments, but field surveys only cover a small fraction of these areas. Satellite-based sensors can repeatedly record the visible and near-infrared reflectance spectra that contain the absorption, scattering, and fluorescence signatures of functional phytoplankton groups, colored dissolved matter, and particulate matter near the surface ocean, and of biologically structured habitats (floating and emergent vegetation, benthic habitats like coral, seagrass, and algae). These measures can be incorporated into Essential Biodiversity Variables (EBVs), including the distribution, abundance, and traits of groups of species populations, and used to evaluate habitat fragmentation. However, current and planned satellites are not designed to observe the EBVs that change rapidly with extreme tides, salinity, temperatures, storms, pollution, or physical habitat destruction over scales relevant to human activity. Making these observations requires a new generation of satellite sensors able to sample with these combined characteristics: (1) spatial resolution on the order of 30 to 100-m pixels or smaller; (2) spectral resolution on the order of 5 nm in the visible and 10 nm in the short-wave infrared spectrum (or at least two or more bands at 1,030, 1,240, 1,630, 2,125, and/or 2,260 nm) for atmospheric correction and aquatic and vegetation assessments; (3) radiometric quality with signal to noise ratios (SNR) above 800 (relative to signal levels typical of the open ocean), 14-bit digitization, absolute radiometric calibration <2%, relative calibration of 0.2%, polarization sensitivity <1%, high radiometric stability and linearity, and operations designed to minimize sunglint; and (4) temporal resolution of hours to days. We refer to these combined specifications as H4 imaging. Enabling H4 imaging is vital for the conservation and management of global biodiversity and ecosystem services, including food provisioning and water security. An agile satellite in a 3-d repeat low-Earth orbit could sample 30-km swath images of several hundred coastal habitats daily. Nine H4 satellites would provide weekly coverage of global coastal zones. Such satellite constellations are now feasible and are used in various applications.

Physiological and ecological drivers of early spring blooms of a coastal phytoplankter

Climate affects the timing and magnitude of phytoplankton blooms that fuel marine food webs and influence global biogeochemical cycles. Changes in bloom timing have been detected in some cases, but the underlying mechanisms remain elusive, contributing to uncertainty in long-term predictions of climate change impacts. Here we describe a 13-year hourly time series from the New England shelf of data on the coastal phytoplankter Synechococcus, during which the timing of its spring bloom varied by 4 weeks. We show that multiyear trends are due to temperature-induced changes in cell division rate, with earlier blooms driven by warmer spring water temperatures. Synechococcus loss rates shift in tandem with division rates, suggesting a balance between growth and loss that has persisted despite phenological shifts and environmental change.

Diversity of Synechococcus at the Martha's Vineyard Coastal Observatory: Insights from Culture Isolations, Clone Libraries, and Flow Cytometry

The cyanobacterium Synechococcus is a ubiquitous, important phytoplankter across the world's oceans. A high degree of genetic diversity exists within the marine group, which likely contributes to its global success. Over 20 clades with different distribution patterns have been identified. However, we do not fully understand the environmental factors that control clade distributions. These factors are likely to change seasonally, especially in dynamic coastal systems. To investigate how coastal Synechococcus assemblages change temporally, we assessed the diversity of Synechococcus at the Martha's Vineyard Coastal Observatory (MVCO) over three annual cycles with culture-dependent and independent approaches. We further investigated the abundance of both phycoerythrin (PE)-containing and phycocyanin (PC)-only Synechococcus with a flow cytometric setup that distinguishes PC-only Synechococcus from picoeukaryotes. We found that the Synechococcus assemblage at MVCO is diverse (13 different clades identified), but dominated by clade I representatives. Many clades were only isolated during late summer and fall, suggesting more favorable conditions for isolation at this time. PC-only strains from four different clades were isolated, but these cells were only detected by flow cytometry in a few samples over the time series, suggesting they are rare at this site. Within clade I, we identified four distinct subclades. The relative abundances of each subclade varied over the seasonal cycle, and the high Synechococcus cell concentration at MVCO may be maintained by the diversity found within this clade. This study highlights the need to understand how temporal aspects of the environment affect Synechococcus community structure and cell abundance.

Rapid growth and concerted sexual transitions by a bloom of the harmful dinoflagellate Alexandrium fundyense (Dinophyceae)

Transitions between life cycle stages by the harmful dinoflagellate Alexandrium fundyense are critical for the initiation and termination of its blooms. To quantify these transitions in a single population, an Imaging FlowCytobot (IFCB), was deployed in Salt Pond (Eastham, Massachusetts), a small, tidally flushed kettle pond that hosts near annual, localized A. fundyense blooms. Machine-based image classifiers differentiating A. fundyense life cycle stages were developed and results were compared to manually corrected IFCB samples, manual microscopy-based estimates of A. fundyense abundance, previously published data describing prevalence of the parasite Amoebophrya, and a continuous culture of A. fundyense infected with Amoebophrya. In Salt Pond, a development phase of sustained vegetative division lasted approximately 3 weeks and was followed by a rapid and near complete conversion to small, gamete cells. The gametic period (∼3 d) coincided with a spike in the frequency of fusing gametes (up to 5% of A. fundyense images) and was followed by a zygotic phase (∼4 d) during which cell sizes returned to their normal range but cell division and diel vertical migration ceased. Cell division during bloom development was strongly phased, enabling estimation of daily rates of division, which were more than twice those predicted from batch cultures grown at similar temperatures in replete medium. Data from the Salt Pond deployment provide the first continuous record of an A. fundyense population through its complete bloom cycle and demonstrate growth and sexual induction rates much higher than are typically observed in culture.

Inversion of spectral absorption coefficients to infer phytoplankton size classes, chlorophyll concentration, and detrital matter

Measured spectral absorption coefficients were inverted to infer phytoplankton concentration in three size classes (picoplankton, nanoplankton, and microplankton), chlorophyll concentration [Chl], and both magnitude and spectral shape of absorption by colored detrital matter (CDM). Our algorithm allowed us to solve for the nonlinear factor of CDM absorption slope separately from the other linear factors, thus fully utilizing the additive characteristic inherent in absorption coefficients. We validated the inversion with three datasets: two spatially distributed global datasets, the Laboratoire d'Océanographie de Villefranche dataset and the NASA bio-Optical Marine Algorithm Dataset, and a time series coastal dataset, the Martha's Vineyard Coastal Observatory dataset. Comparison with high performance liquid chromatography analyses showed that the phytoplankton size classes can be retrieved with correlation coefficients (r)>0.7, root mean square errors of 0.2, and median relative errors of 20% in oceanic waters and with similar performance in coastal waters. Much improved agreement was found for the entire phytoplankton population, with r>0.90 for [Chl] and absorption coefficients (a_ph) for all three datasets. The inferred a_CDM(400) and CDM spectral slope agree within ±4% of measurements in both oceanic and coastal waters. The results indicate that the chlorophyll-a specific absorption spectra used as an inversion kernel represent well the global mean states for each of the three phytoplankton size classes. The method can be applied to either bulk or particulate absorption data and is spectrally flexible.

The Marine Microbial Eukaryote Transcriptome Sequencing Project (MMETSP): illuminating the functional diversity of eukaryotic life in the oceans through transcriptome sequencing

Current sampling of genomic sequence data from eukaryotes is relatively poor, biased, and inadequate to address important questions about their biology, evolution, and ecology; this Community Page describes a resource of 700 transcriptomes from marine microbial eukaryotes to help understand their role in the world's oceans.

Complexities of bloom dynamics in the toxic dinoflagellate Alexandrium fundyense revealed through DNA measurements by imaging flow cytometry coupled with species-specific rRNA probes

Measurements of the DNA content of different protist populations can shed light on a variety of processes, including cell division, sex, prey ingestion, and parasite invasion. Here, we modified an Imaging FlowCytobot (IFCB), a custom-built flow cytometer that records images of microplankton, to measure the DNA content of large dinoflagellates and other high-DNA content species. The IFCB was also configured to measure fluorescence from Cy3-labeled rRNA probes, aiding the identification of Alexandrium fundyense (syn. A. tamarense Group I), a photosynthetic dinoflagellate that causes paralytic shellfish poisoning (PSP). The modified IFCB was used to analyze samples from the development, peak and termination phases of an inshore A. fundyense bloom (Salt Pond, Eastham, MA USA), and from a rare A. fundyense 'red tide' that occurred in the western Gulf of Maine, offshore of Portsmouth, NH (USA). Diploid or G2 phase ('2C') A. fundyense cells were frequently enriched at the near-surface, suggesting an important role for aggregation at the air-sea interface during sexual events. Also, our analysis showed that large proportions of A. fundyense cells in both the Salt Pond and red tide blooms were planozygotes during bloom decline, highlighting the importance of sexual fusion to bloom termination. At Salt Pond, bloom decline also coincided with a dramatic rise in infections by the parasite genus Amoebophrya. The samples that were most heavily infected contained many large cells with higher DNA-associated fluorescence than 2C vegetative cells, but these cells' nuclei were also frequently consumed by Amoebophrya trophonts. Neither large cell size nor increased DNA-associated fluorescence could be replicated by infecting an A. fundyense culture of vegetative cells. Therefore we attribute these characteristics of the large Salt Pond cells to planozygote maturation rather than Amoebophrya infection, though an interaction between infection and planozygote maturation may also have contributed. The modified IFCB is a valuable tool for exploring the conditions that promote sexual transitions by dinoflagellate blooms but care is needed when interpreting results from samples in which parasitism is prevalent.

Continuous automated imaging-in-flow cytometry for detection and early warning of Karenia brevis blooms in the Gulf of Mexico

Monitoring programs for harmful algal blooms (HABs) typically rely on time-consuming manual methods for identification and enumeration of phytoplankton, which make it difficult to obtain results with sufficient temporal resolution for early warning. Continuous automated imaging-in-flow by the Imaging FlowCytobot (IFCB) deployed at Port Aransas, TX has provided early warnings of six HAB events. Here we describe the progress in automating this early warning system for blooms of Karenia brevis. In 2009, manual inspection of IFCB images in mid-August 2009 provided early warning for a Karenia bloom that developed in mid-September. Images from 2009 were used to develop an automated classifier that was employed in 2011. Successful implementation of automated file downloading, processing and image classification allowed results to be available within 4 h after collection and to be sent to state agency representatives by email for early warning of HABs. No human illness (neurotoxic shellfish poisoning) has resulted from these events. In contrast to the common assumption that Karenia blooms are near monospecific, post-bloom analysis of the time series revealed that Karenia cells comprised at most 60-75 % of the total microplankton.

Taxonomic classification of phytoplankton with multivariate optical computing, part III: demonstration

We describe the automatic analysis of fluorescence tracks of phytoplankton recorded with a fluorescence imaging photometer. The optical components and construction of the photometer were described in Part I and Part II of this series in this issue. An algorithm first isolates tracks corresponding to a single phytoplankter transit in the nominal focal plane of a flow cell. Then, the fluorescence streaks in the track that correspond to individual optical elements on the filter wheel are identified. The fluorescence intensity of each streak is integrated and used to calculate ratios. This approach was tested using 853 fluorescence measurements of the coccolithophore Emiliania huxleyi and the diatom Thalassiosira pseudonana. Average intensity ratios for the two classes closely follow those predicted in Part I of this series, with a distribution of ratios in each class that is consistent with the signal-to-noise ratio calculations in Part II for single cells. No overlap of the two class ratios was observed, yielding perfect classification.

Massive phytoplankton blooms under Arctic sea ice

Phytoplankton blooms over Arctic Ocean continental shelves are thought to be restricted to waters free of sea ice. Here, we document a massive phytoplankton bloom beneath fully consolidated pack ice far from the ice edge in the Chukchi Sea, where light transmission has increased in recent decades because of thinning ice cover and proliferation of melt ponds. The bloom was characterized by high diatom biomass and rates of growth and primary production. Evidence suggests that under-ice phytoplankton blooms may be more widespread over nutrient-rich Arctic continental shelves and that satellite-based estimates of annual primary production in these waters may be underestimated by up to 10-fold.

PHYLOGENETIC ANALYSIS OF BRACHIDINIUM CAPITATUM (DINOPHYCEAE) FROM THE GULF OF MEXICO INDICATES MEMBERSHIP IN THE KARENIACEAE(1)

Brachidinium capitatum F. J. R. Taylor, typically considered a rare oceanic dinoflagellate, and one which has not been cultured, was observed at elevated abundances (up to 65 cells · mL(-1) ) at a coastal station in the western Gulf of Mexico in the fall of 2007. Continuous data from the Imaging FlowCytobot (IFCB) provided cell images that documented the bloom during 3 weeks in early November. Guided by IFCB observations, field collection permitted phylogenetic analysis and evaluation of the relationship between Brachidinium and Karenia. Sequences from SSU, LSU, internal transcribed spacer (ITS), and cox1 regions for B. capitatum were compared with five other species of Karenia; all B. capitatum sequences were unique but supported its placement within the Kareniaceae. From a total of 71,487 images, data on the timing and frequency of dividing cells was also obtained for B. capitatum, allowing the rate of division for B. capitatum to be estimated. The maximum daily growth rate estimate was 0.22 d(-1) . Images showed a range in morphological variability, with the position of the four major processes highly variable. The combination of morphological features similar to the genus Karenia and a phylogenetic analysis placing B. capitatum in the Karenia clade leads us to propose moving the genus Brachidinium into the Kareniaceae. However, the lack of agreement among individual gene phylogenies suggests that the inclusion of different genes and more members of the genus Karenia are necessary before a final determination regarding the validity of the genus Brachidinium can be made.

Flow cytometric determination of size and complex refractive index for marine particles: comparison with independent and bulk estimates

We advance a method to determine the diameter D and the complex refractive index (n + n'i) of marine particles from flow cytometric measurements of forward scattering, side scattering, and chlorophyll fluorescence combined with Mie theory. To understand better the application of Mie theory with its assumptions to flow cytometry (FCM) measurements of phytoplankton cells, we evaluate our flow cytometric-Mie (FCM-Mie) method by comparing results from a variety of phytoplankton cultures with independent estimates of cell D and with estimates of n and n' from the inversion of bulk measurements. Cell D initially estimated from the FCM-Mie method is lower than independent estimates, and n and n' are generally higher than bulk estimates. These differences reflect lower forward scattering and higher side scattering for single-cell measurements than predicted by Mie theory. The application of empirical scattering corrections improves FCM-Mie estimates of cell size, n, and n'; notably size is determined accurately for cells grown in both high- and low-light conditions, and n' is correlated with intracellular chlorophyll concentration. A comparison of results for phytoplankton and mineral particles suggests that differences in n between these particle types can be determined from FCM measurements. In application to natural mixtures of particles, eukaryotic pico/nanophytoplankton and Synechococcus have minimum mean values of n' in surface waters, and nonphytoplankton particles have higher values of n than phytoplankton at all depths.

Photosynthetic characteristics of marine phytoplankton from pump-during-probe fluorometry of individual cells at sea

BACKGROUND

Active fluorescence techniques are becoming commonly used to monitor the state of the photosynthetic apparatus in natural populations of phytoplankton, but at present these are bulk water measurements that average all the fluorescent material in each sample. Here we describe two instruments that combine individual-cell "pump-during-probe" (PDP) measurements of chlorophyll (Chl) fluorescence induction, on the time scale of 30 to 100 micros, with flow cytometric or visual characterization of each cell.

METHODS

In the PDP flow cytometer, we measure Chl fluorescence yield as a function of time during a 150 micros excitation flash provided by an argon ion laser; each particle is subsequently classified as in a conventional flow cytometer. In the PDP microfluorometer, individual cells in a sample chamber are visually identified, and fluorescence excitation is provided by a blue light-emitting diode that can be configured to provide a saturating flash and also a subsequent series of short flashlets. This sequence allows both saturation and relaxation kinetics to be monitored.

RESULTS

Phytoplankton from natural samples and on-deck iron-enrichment incubation experiments in the Southern Ocean were examined with each PDP instrument, providing estimates of the potential quantum yield of photochemistry and the functional absorption cross section for photosystem 2, for either individuals (for cells larger than a few micrometers) or populations (for smaller cells).

CONCLUSIONS

Results from initial field applications indicate that single-cell PDP measurements can be a powerful tool for investigating the nutritional state of phytoplankton cells and the regulation of phytoplankton growth in the sea.