Molecular Forecasting of Domoic Acid during a Pervasive Toxic Diatom Bloom

In 2015, the largest recorded harmful algal bloom (HAB) occurred in the Northeast Pacific, causing nearly 100 million dollars in damages to fisheries and killing many protected marine mammals. Dominated by the toxic diatom Pseudo-nitzschia australis , this bloom produced high levels of the neurotoxin domoic acid (DA). Through molecular and transcriptional characterization of 52 near-weekly phytoplankton net-tow samples collected at a bloom hotspot in Monterey Bay, California, we identified active transcription of known DA biosynthesis ( dab ) genes from the three identified toxigenic species, including P. australis as the primary origin of toxicity. Elevated expression of silicon transporters ( sit1 ) during the bloom supports the previously hypothesized role of dissolved silica (Si) exhaustion in contributing to bloom physiology and toxicity. We find that co-expression of the dabA and sit1 genes serves as a robust predictor of DA one week in advance, potentially enabling the forecasting of DA-producing HABs. We additionally present evidence that low levels of iron could have co-limited the diatom population along with low Si. Iron limitation represents a previously unrecognized driver of both toxin production and ecological success of the low iron adapted Pseudo-nitzschia genus during the 2015 bloom, and increasing pervasiveness of iron limitation may fuel the escalating magnitude and frequency of toxic Pseudo-nitzschia blooms globally. Our results advance understanding of bloom physiology underlying toxin production, bloom prediction, and the impact of global change on toxic blooms.

Significance

Pseudo-nitzschia diatoms form oceanic harmful algal blooms that threaten human health through production of the neurotoxin domoic acid (DA). DA biosynthetic gene expression is hypothesized to control DA production in the environment, yet what regulates expression of these genes is yet to be discovered. In this study, we uncovered expression of DA biosynthesis genes by multiple toxigenic Pseudo-nitzschia species during an economically impactful bloom along the North American West Coast, and identified genes that predict DA in advance of its production. We discovered that iron and silica co-limitation restrained the bloom and likely promoted toxin production. This work suggests that increasing iron limitation due to global change may play a previously unrecognized role in driving bloom frequency and toxicity.

A cross-regional examination of patterns and environmental drivers of Pseudo-nitzschia harmful algal blooms along the California coast

Pseudo-nitzschia species with the ability to produce the neurotoxin domoic acid (DA) are the main cause of harmful algal blooms (HABs) along the U.S. West Coast, with major impacts on ecosystems, fisheries, and human health. While most Pseudo-nitzschia (PN) HAB studies to date have focused on their characteristics at specific sites, few cross-regional comparisons exist, and mechanistic understanding of large-scale HAB drivers remains incomplete. To close these gaps, we compiled a nearly 20-year time series of in situ particulate DA and environmental observations to characterize similarities and differences in PN HAB drivers along the California coast. We focus on three DA hotspots with the greatest data density: Monterey Bay, the Santa Barbara Channel, and the San Pedro Channel. Coastwise, DA outbreaks are strongly correlated with upwelling, chlorophyll-a, and silicic acid limitation relative to other nutrients. Clear differences also exist across the three regions, with contrasting responses to climate regimes across a north to south gradient. In Monterey Bay, PN HAB frequency and intensity increase under relatively nutrient-poor conditions during anomalously low upwelling intensities. In contrast, in the Santa Barbara and San Pedro Channels, PN HABs are favored under cold, nitrogen-rich conditions during more intense upwelling. These emerging patterns provide insights on ecological drivers of PN HABs that are consistent across regions and support the development of predictive capabilities for DA outbreaks along the California coast and beyond.

Coastal upwelling drives ecosystem temporal variability from the surface to the abyssal seafloor

Long-term biological time series that monitor ecosystems across the ocean's full water column are extremely rare. As a result, classic paradigms are yet to be tested. One such paradigm is that variations in coastal upwelling drive changes in marine ecosystems throughout the water column. We examine this hypothesis by using data from three multidecadal time series spanning surface (0 m), midwater (200 to 1,000 m), and benthic (~4,000 m) habitats in the central California Current Upwelling System. Data include microscopic counts of surface plankton, video quantification of midwater animals, and imaging of benthic seafloor invertebrates. Taxon-specific plankton biomass and midwater and benthic animal densities were separately analyzed with principal component analysis. Within each community, the first mode of variability corresponds to most taxa increasing and decreasing over time, capturing seasonal surface blooms and lower-frequency midwater and benthic variability. When compared to local wind-driven upwelling variability, each community correlates to changes in upwelling damped over distinct timescales. This suggests that periods of high upwelling favor increase in organism biomass or density from the surface ocean through the midwater down to the abyssal seafloor. These connections most likely occur directly via changes in primary production and vertical carbon flux, and to a lesser extent indirectly via other oceanic changes. The timescales over which species respond to upwelling are taxon-specific and are likely linked to the longevity of phytoplankton blooms (surface) and of animal life (midwater and benthos), which dictate how long upwelling-driven changes persist within each community.

The founding charter of the Omic Biodiversity Observation Network (Omic BON)

Omic BON is a thematic Biodiversity Observation Network under the Group on Earth Observations Biodiversity Observation Network (GEO BON), focused on coordinating the observation of biomolecules in organisms and the environment. Our founding partners include representatives from national, regional, and global observing systems; standards organizations; and data and sample management infrastructures. By coordinating observing strategies, methods, and data flows, Omic BON will facilitate the co-creation of a global omics meta-observatory to generate actionable knowledge. Here, we present key elements of Omic BON's founding charter and first activities.

Pelagic calcifiers face increased mortality and habitat loss with warming and ocean acidification

Global change is impacting the oceans in an unprecedented way, and multiple lines of evidence suggest that species distributions are changing in space and time. There is increasing evidence that multiple environmental stressors act together to constrain species habitat more than expected from warming alone. Here, we conducted a comprehensive study of how temperature and aragonite saturation state act together to limit Limacina helicina, globally distributed pteropods that are ecologically important pelagic calcifiers and an indicator species for ocean change. We co-validated three different approaches to evaluate the impact of ocean warming and acidification (OWA) on the survival and distribution of this species in the California Current Ecosystem. First, we used colocated physical, chemical, and biological data from three large-scale west coast cruises and regional time series; second, we conducted multifactorial experimental incubations to evaluate how OWA impacts pteropod survival; and third, we validated the relationships we found against global distributions of pteropods and carbonate chemistry. OWA experimental work revealed mortality increases under OWA, while regional habitat suitability indices and global distributions of L. helicina suggest that a multi-stressor framework is essential for understanding pteropod distributions. In California Current Ecosystem habitats, where pteropods are living close to their thermal maximum already, additional warming and acidification through unabated fossil fuel emissions (RCP 8.5) are expected to dramatically reduce habitat suitability.

Environmental DNA reveals the fine-grained and hierarchical spatial structure of kelp forest fish communities

Biodiversity is changing at an accelerating rate at both local and regional scales. Beta diversity, which quantifies species turnover between these two scales, is emerging as a key driver of ecosystem function that can inform spatial conservation. Yet measuring biodiversity remains a major challenge, especially in aquatic ecosystems. Decoding environmental DNA (eDNA) left behind by organisms offers the possibility of detecting species sans direct observation, a Rosetta Stone for biodiversity. While eDNA has proven useful to illuminate diversity in aquatic ecosystems, its utility for measuring beta diversity over spatial scales small enough to be relevant to conservation purposes is poorly known. Here we tested how eDNA performs relative to underwater visual census (UVC) to evaluate beta diversity of marine communities. We paired UVC with 12S eDNA metabarcoding and used a spatially structured hierarchical sampling design to assess key spatial metrics of fish communities on temperate rocky reefs in southern California. eDNA provided a more-detailed picture of the main sources of spatial variation in both taxonomic richness and community turnover, which primarily arose due to strong species filtering within and among rocky reefs. As expected, eDNA detected more taxa at the regional scale (69 vs. 38) which accumulated quickly with space and plateaued at only ~ 11 samples. Conversely, the discovery rate of new taxa was slower with no sign of saturation for UVC. Based on historical records in the region (2000-2018) we found that 6.9 times more UVC samples would be required to detect 50 taxa compared to eDNA. Our results show that eDNA metabarcoding can outperform diver counts to capture the spatial patterns in biodiversity at fine scales with less field effort and more power than traditional methods, supporting the notion that eDNA is a critical scientific tool for detecting biodiversity changes in aquatic ecosystems.

Seasonal and Geographical Transitions in Eukaryotic Phytoplankton Community Structure in the Atlantic and Pacific Oceans

Much is known about how broad eukaryotic phytoplankton groups vary according to nutrient availability in marine ecosystems. However, genus- and species-level dynamics are generally unknown, although important given that adaptation and acclimation processes differentiate at these levels. We examined phytoplankton communities across seasonal cycles in the North Atlantic (BATS) and under different trophic conditions in the eastern North Pacific (ENP), using phylogenetic classification of plastid-encoded 16S rRNA amplicon sequence variants (ASVs) and other methodologies, including flow cytometric cell sorting. Prasinophytes dominated eukaryotic phytoplankton amplicons during the nutrient-rich deep-mixing winter period at BATS. During stratification (‘summer’) uncultured dictyochophytes formed ∼35 ± 10% of all surface plastid amplicons and dominated those from stramenopile algae, whereas diatoms showed only minor, ephemeral contributions over the entire year. Uncultured dictyochophytes also comprised a major fraction of plastid amplicons in the oligotrophic ENP. Phylogenetic reconstructions of near-full length 16S rRNA sequences established 11 uncultured Dictyochophyte Environmental Clades (DEC). DEC-I and DEC-VI dominated surface dictyochophytes under stratification at BATS and in the ENP, and DEC-IV was also important in the latter. Additionally, although less common at BATS, Florenciella-related clades (FC) were prominent at depth in the ENP. In both ecosystems, pelagophytes contributed notably at depth, with PEC-VIII (Pelagophyte Environmental Clade) and (cultured) Pelagomonas calceolata being most important. Q-PCR confirmed the near absence of P. calceolata at the surface of the same oligotrophic sites where it reached ∼1,500 18S rRNA gene copies ml^–1 at the DCM. To further characterize phytoplankton present in our samples, we performed staining and at-sea single-cell sorting experiments. Sequencing results from these indicated several uncultured dictyochophyte clades are comprised of predatory mixotrophs. From an evolutionary perspective, these cells showed both conserved and unique features in the chloroplast genome. In ENP metatranscriptomes we observed high expression of multiple chloroplast genes as well as expression of a selfish element (group II intron) in the psaA gene. Comparative analyses across the Pacific and Atlantic sites support the conclusion that predatory dictyochophytes thrive under low nutrient conditions. The observations that several uncultured dictyochophyte lineages are seemingly capable of photosynthesis and predation, raises questions about potential shifts in phytoplankton trophic roles associated with seasonality and long-term ocean change.

Depth-Differentiation and Seasonality of Planktonic Microbial Assemblages in the Monterey Bay Upwelling System

Coastal upwelling regions are hotspots of biological productivity, supporting diverse communities of microbial life and metabolisms. Monterey Bay (MB), a coastal ocean embayment in central California, experiences seasonal upwelling of cold, nutrient-rich waters that sustain episodes of high phytoplankton production in surface waters. While productivity in surface waters is intimately linked to metabolisms of diverse communities of Archaea and Bacteria, a comprehensive understanding of the microbial community in MB is missing thus far, particularly in relation to the distinct hydrographic seasons characteristic of the MB system. Here we present the results of a 2-year microbial time-series survey in MB, investigating community composition and structure across spatiotemporal gradients. In deciphering these patterns, we used unique sequence variants (SVs) of the 16S rRNA gene (V4–V5 region), complemented with metagenomes and metatranscriptomes representing multiple depth profiles. We found clear depth-differentiation and recurring seasonal abundance patterns within planktonic communities, particularly when analyzed at finer taxonomic levels. Compositional changes were more pronounced in the upper 0–40 m of the water column, whereas deeper depths were characterized by temporally stable populations. In accordance with the dynamic nutrient profiles, the system appears to change from a Bacteroidetes- and Rhodobacterales-dominated upwelling period to an oceanic season dominated by oligotrophic groups such as SAR11 and picocyanobacteria. The cascade of environmental changes brought about by upwelling and relaxation events thus impacts microbial community structure in the bay, with important implications for the temporal variability of nutrient and energy fluxes within the MB ecosystem. Our observations emphasize the need for continued monitoring of planktonic microbial communities in order to predict and manage the behavior of this sensitive marine sanctuary ecosystem, over projected intensification of upwelling in the region.

Environmental DNA reveals seasonal shifts and potential interactions in a marine community

Environmental DNA (eDNA) analysis allows the simultaneous examination of organisms across multiple trophic levels and domains of life, providing critical information about the complex biotic interactions related to ecosystem change. Here we used multilocus amplicon sequencing of eDNA to survey biodiversity from an eighteen-month (2015–2016) time-series of seawater samples from Monterey Bay, California. The resulting dataset encompasses 663 taxonomic groups (at Family or higher taxonomic rank) ranging from microorganisms to mammals. We inferred changes in the composition of communities, revealing putative interactions among taxa and identifying correlations between these communities and environmental properties over time. Community network analysis provided evidence of expected predator-prey relationships, trophic linkages, and seasonal shifts across all domains of life. We conclude that eDNA-based analyses can provide detailed information about marine ecosystem dynamics and identify sensitive biological indicators that can suggest ecosystem changes and inform conservation strategies.

Increasingly, eDNA is being used to infer ecological interactions. Here the authors sample eDNA over 18 months in a marine environment and use co-occurrence network analyses to infer potential interactions among organisms from microbes to mammals, testing how they change over time in response to oceanographic factors.

Humpback whale song occurrence reflects ecosystem variability in feeding and migratory habitat of the northeast Pacific

This study examines the occurrence of humpback whale (Megaptera novaeangliae) song in the northeast Pacific from three years of continuous recordings off central California (36.713°N, 122.186°W). Song is prevalent in this feeding and migratory habitat, spanning nine months of the year (September-May), peaking in winter (November-January), and reaching a maximum of 86% temporal coverage (during November 2017). From the rise of song in fall through the end of peak occurrence in winter, song length increases significantly from month to month. The seasonal peak in song coincides with the seasonal trough in day length and sighting-based evidence of whales leaving Monterey Bay, consistent with seasonal migration. During the seasonal song peak, diel variation shows maximum occurrence at night (69% of the time), decreasing during dawn and dusk (52%), and further decreasing with increasing solar elevation during the day, reaching a minimum near solar noon (30%). Song occurrence increased 44% and 55% between successive years. Sighting data within the acoustic detection range of the hydrophone indicate that variation in local population density was an unlikely cause of this large interannual variation. Hydrographic data and modeling of acoustic transmission indicate that changes in neither habitat occupancy nor acoustic transmission were probable causes. Conversely, the positive interannual trend in song paralleled major ecosystem variations, including similarly large positive trends in wind-driven upwelling, primary productivity, and krill abundance. Further, the lowest song occurrence during the first year coincided with anomalously warm ocean temperatures and an extremely toxic harmful algal bloom that affected whales and other marine mammals in the region. These major ecosystem variations may have influenced the health and behavior of humpback whales during the study period.

Depth distributions of nitrite reductase (nirK) gene variants reveal spatial dynamics of thaumarchaeal ecotype populations in coastal Monterey Bay

Ammonia-oxidizing archaea (AOA) of the phylum Thaumarchaeota are key players in nutrient cycling, yet large gaps remain in our understanding of their ecology and metabolism. Despite multiple lines of evidence pointing to a central role for copper-containing nitrite reductase (NirK) in AOA metabolism, the thaumarchaeal nirK gene is rarely studied in the environment. In this study, we examine the diversity of nirK in the marine pelagic environment, in light of previously described ecological patterns of pelagic thaumarchaeal populations. Phylogenetic analyses show that nirK better resolves diversification patterns of marine Thaumarchaeota, compared to the conventionally used marker gene amoA. Specifically, we demonstrate that the three major phylogenetic clusters of marine nirK correspond to the three 'ecotype' populations of pelagic Thaumarchaeota. In this context, we further examine the relative distributions of the three variant groups in metagenomes and metatranscriptomes representing two depth profiles in coastal Monterey Bay. Our results reveal that nirK effectively tracks the dynamics of thaumarchaeal ecotype populations, particularly finer-scale diversification patterns within major lineages. We also find evidence for multiple copies of nirK per genome in a fraction of thaumarchaeal cells in the water column, which must be taken into account when using it as a molecular marker.

Impacts of urban carbon dioxide emissions on sea-air flux and ocean acidification in nearshore waters

Greatly enhanced atmospheric carbon dioxide (CO₂) levels relative to well-mixed marine air are observed during periods of offshore winds at coastal sensor platforms in Monterey Bay, California, USA. The highest concentrations originate from urban and agricultural areas, are driven by diurnal winds, and peak in the early morning. These enhanced atmospheric levels can be detected across a ~100km wide nearshore area and represent a significant addition to total oceanic CO₂ uptake. A global estimate puts the added sea-air flux of CO₂ from these greatly enhanced atmospheric CO₂ levels at 25 million tonnes, roughly 1% of the ocean’s annual CO₂ uptake. The increased uptake over the 100 km coastal swath is of order 20%, indicating a potentially large impact on ocean acidification in productive coastal waters.

A Sixteen-year Decline in Dissolved Oxygen in the Central California Current

A potential consequence of climate change is global decrease in dissolved oxygen at depth in the oceans due to changes in the balance of ventilation, mixing, respiration, and photosynthesis. We present hydrographic cruise observations of declining dissolved oxygen collected along CalCOFI Line 66.7 (Line 67) off of Monterey Bay, in the Central California Current region, and investigate likely mechanisms. Between 1998 and 2013, dissolved oxygen decreased at the mean rate of 1.92 µmol kg⁻¹ year⁻¹ on σ_θ 26.6–26.8 kg m⁻³ isopycnals (250–400 m), translating to a 40% decline from initial concentrations. Two cores of elevated dissolved oxygen decline at 130 and 240 km from shore, which we suggest are a California Undercurrent and a California Current signal respectively, occurred on σ_θ ranges of 26.0–26.8 kg m⁻³ (100–400 m). A box model suggests that small annual changes in dissolved oxygen in source regions are sufficient to be the primary driver of the mid-depth declines. Variation in dissolved oxygen at the bottom of the surface mixed layer suggests that there is also a signal of increased local remineralization.

Declining oxygen in the global ocean and coastal waters

Oxygen is fundamental to life. Not only is it essential for the survival of individual animals, but it regulates global cycles of major nutrients and carbon. The oxygen content of the open ocean and coastal waters has been declining for at least the past half-century, largely because of human activities that have increased global temperatures and nutrients discharged to coastal waters. These changes have accelerated consumption of oxygen by microbial respiration, reduced solubility of oxygen in water, and reduced the rate of oxygen resupply from the atmosphere to the ocean interior, with a wide range of biological and ecological consequences. Further research is needed to understand and predict long-term, global- and regional-scale oxygen changes and their effects on marine and estuarine fisheries and ecosystems.

Biomonitoring of marine vertebrates in Monterey Bay using eDNA metabarcoding

Molecular analysis of environmental DNA (eDNA) can be used to assess vertebrate biodiversity in aquatic systems, but limited work has applied eDNA technologies to marine waters. Further, there is limited understanding of the spatial distribution of vertebrate eDNA in marine waters. Here, we use an eDNA metabarcoding approach to target and amplify a hypervariable region of the mitochondrial 12S rRNA gene to characterize vertebrate communities at 10 oceanographic stations spanning 45 km within the Monterey Bay National Marine Sanctuary (MBNMS). In this study, we collected three biological replicates of small volume water samples (1 L) at 2 depths at each of the 10 stations. We amplified fish mitochondrial DNA using a universal primer set. We obtained 5,644,299 high quality Illumina sequence reads from the environmental samples. The sequence reads were annotated to the lowest taxonomic assignment using a bioinformatics pipeline. The eDNA survey identified, to the lowest taxonomic rank, 7 families, 3 subfamilies, 10 genera, and 72 species of vertebrates at the study sites. These 92 distinct taxa come from 33 unique marine vertebrate families. We observed significantly different vertebrate community composition between sampling depths (0 m and 20/40 m deep) across all stations and significantly different communities at stations located on the continental shelf (<200 m bottom depth) versus in the deeper waters of the canyons of Monterey Bay (>200 m bottom depth). All but 1 family identified using eDNA metabarcoding is known to occur in MBNMS. The study informs the implementation of eDNA metabarcoding for vertebrate biomonitoring.

Interacting environmental mosaics drive geographic variation in mussel performance and predation vulnerability

Although theory suggests geographic variation in species' performance is determined by multiple niche parameters, little consideration has been given to the spatial structure of interacting stressors that may shape local and regional vulnerability to global change. Here, we use spatially explicit mosaics of carbonate chemistry, food availability and temperature spanning 1280 km of coastline to test whether persistent, overlapping environmental mosaics mediate the growth and predation vulnerability of a critical foundation species, the mussel Mytilus californianus. We find growth was highest and predation vulnerability was lowest in dynamic environments with frequent exposure to low pH seawater and consistent food. In contrast, growth was lowest and predation vulnerability highest when exposure to low pH seawater was decoupled from high food availability, or in exceptionally warm locations. These results illustrate how interactions among multiple drivers can cause unexpected, yet persistent geographic mosaics of species performance, interactions and vulnerability to environmental change.

Interacting environmental mosaics drive geographic variation in mussel performance and predation vulnerability

Although theory suggests geographic variation in species' performance is determined by multiple niche parameters, little consideration has been given to the spatial structure of interacting stressors that may shape local and regional vulnerability to global change. Here, we use spatially explicit mosaics of carbonate chemistry, food availability and temperature spanning 1280 km of coastline to test whether persistent, overlapping environmental mosaics mediate the growth and predation vulnerability of a critical foundation species, the mussel Mytilus californianus. We find growth was highest and predation vulnerability was lowest in dynamic environments with frequent exposure to low pH seawater and consistent food. In contrast, growth was lowest and predation vulnerability highest when exposure to low pH seawater was decoupled from high food availability, or in exceptionally warm locations. These results illustrate how interactions among multiple drivers can cause unexpected, yet persistent geographic mosaics of species performance, interactions and vulnerability to environmental change.

Climate mediates hypoxic stress on fish diversity and nursery function at the land-sea interface

Coastal ecosystems provide numerous important ecological services, including maintenance of biodiversity and nursery grounds for many fish species of ecological and economic importance. However, human population growth has led to increased pollution, ocean warming, hypoxia, and habitat alteration that threaten ecosystem services. In this study, we used long-term datasets of fish abundance, water quality, and climatic factors to assess the threat of hypoxia and the regulating effects of climate on fish diversity and nursery conditions in Elkhorn Slough, a highly eutrophic estuary in central California (United States), which also serves as a biodiversity hot spot and critical nursery grounds for offshore fisheries in a broader region. We found that hypoxic conditions had strong negative effects on extent of suitable fish habitat, fish species richness, and abundance of the two most common flatfish species, English sole (Parophrys vetulus) and speckled sanddab (Citharichthys stigmaeus). The estuary serves as an important nursery ground for English sole, making this species vulnerable to anthropogenic threats. We determined that estuarine hypoxia was associated with significant declines in English sole nursery habitat, with cascading effects on recruitment to the offshore adult population and fishery, indicating that human land use activities can indirectly affect offshore fisheries. Estuarine hypoxic conditions varied spatially and temporally and were alleviated by strengthening of El Niño conditions through indirect pathways, a consistent result in most estuaries across the northeast Pacific. These results demonstrate that changes to coastal land use and climate can fundamentally alter the diversity and functioning of coastal nurseries and their adjacent ocean ecosystems.

Cyanobacterial distributions along a physico-chemical gradient in the Northeastern Pacific Ocean

The cyanobacteria Prochlorococcus and Synechococcus are important marine primary producers. We explored their distributions and covariance along a physico-chemical gradient from coastal to open ocean waters in the Northeastern Pacific Ocean. An inter-annual pattern was delineated in the dynamic transition zone where upwelled and eastern boundary current waters mix, and two new Synechococcus clades, Eastern Pacific Clade (EPC) 1 and EPC2, were identified. By applying state-of-the-art phylogenetic analysis tools to bar-coded 16S amplicon datasets, we observed higher abundance of Prochlorococcus high-light I (HLI) and low-light I (LLI) in years when more oligotrophic water intruded farther inshore, while under stronger upwelling Synechococcus I and IV dominated. However, contributions of some cyanobacterial clades were proportionally relatively constant, e.g. Synechococcus EPC2. In addition to supporting observations that Prochlorococcus LLI thrive at higher irradiances than other LL taxa, the results suggest LLI tolerate lower temperatures than previously reported. The phylogenetic precision of our 16S rRNA gene analytical approach and depth of bar-coded sequencing also facilitated detection of clades at low abundance in unexpected places. These include Prochlorococcus at the coast and Cyanobium-related sequences offshore, although it remains unclear whether these came from resident or potentially advected cells. Our study enhances understanding of cyanobacterial distributions in an ecologically important eastern boundary system.

Ammonium uptake by phytoplankton regulates nitrification in the sunlit ocean

Nitrification, the microbial oxidation of ammonium to nitrate, is a central part of the nitrogen cycle. In the ocean’s surface layer, the process alters the distribution of inorganic nitrogen species available to phytoplankton and produces nitrous oxide. A widely held idea among oceanographers is that nitrification is inhibited by light in the ocean. However, recent evidence that the primary organisms involved in nitrification, the ammonia-oxidizing archaea (AOA), are present and active throughout the surface ocean has challenged this idea. Here we show, through field experiments coupling molecular genetic and biogeochemical approaches, that competition for ammonium with phytoplankton is the strongest regulator of nitrification in the photic zone. During multiday experiments at high irradiance a single ecotype of AOA remained active in the presence of rapidly growing phytoplankton. Over the course of this three day experiment, variability in the intensity of competition with phytoplankton caused nitrification rates to decline from those typical of the lower photic zone (60 nmol L⁻¹ d⁻¹) to those in well-lit layers (<1 nmol L⁻¹ d⁻¹). During another set of experiments, nitrification rates exhibited a diel periodicity throughout much of the photic zone, with the highest rates occurring at night when competition with phytoplankton is lowest. Together, the results of our experiments indicate that nitrification rates in the photic zone are more strongly regulated by competition with phytoplankton for ammonium than they are by light itself. This finding advances our ability to model the impact of nitrification on estimates of new primary production, and emphasizes the need to more strongly consider the effects of organismal interactions on nutrient standing stocks and biogeochemical cycling in the surface of the ocean.

Seasonal change in the abundance of Synechococcus and multiple distinct phylotypes in Monterey Bay determined by rbcL and narB quantitative PCR

Synechococcus is a cosmopolitan marine cyanobacterial genus, and is often the most abundant picocyanobacterial genus in coastal waters. Little is known about Synechococcus seasonal dynamics in coastal zones highly impacted by upwelling. This was investigated by collecting seasonal samples from an upwelling-impacted Monterey Bay (MB) monitoring station M0, in parallel with measurements of oceanographic conditions during 2006-2008. Synechococcus abundances were determined using quantitative PCR (qPCR) assays and flow cytometry (FCM). A new qPCR assay was designed to target dominant Synechococcus in MB using the rbcL gene, while previously designed assays targeted distinct phylotypes (called narB subgroups) with the narB gene. The rbcL qPCR assay successfully tracked abundant Synechococcus in MB, accounting for on average 89% (± 57%) of FCM-based counts. Annual spring upwelling caused decreases in Synechococcus and narB subgroup abundances. Differences in narB subgroup abundance maxima and abundance patterns support the view that subgroups differ in their ecologies, including subgroup D_C1, which seems to specifically thrive in coastal waters. Correlations between narB subgroup abundances and measured environmental variables were similar among the subgroups. Therefore, non-measured environmental factors (e.g. metals, mortality) likely had different influences on subgroups, which led to their distinct abundance patterns at M0.

Differential distributions of synechococcus subgroups across the california current system

Synechococcus is an abundant marine cyanobacterial genus composed of different populations that vary physiologically. Synechococcus narB gene sequences (encoding for nitrate reductase in cyanobacteria) obtained previously from isolates and the environment (e.g., North Pacific Gyre Station ALOHA, Hawaii or Monterey Bay, CA, USA) were used to develop quantitative PCR (qPCR) assays. These qPCR assays were used to quantify populations from specific narB phylogenetic clades across the California Current System (CCS), a region composed of dynamic zones between a coastal-upwelling zone and the oligotrophic Pacific Ocean. Targeted populations (narB subgroups) had different biogeographic patterns across the CCS, which appear to be driven by environmental conditions. Subgroups C_C1, D_C1, and D_C2 were abundant in coastal-upwelling to coastal-transition zone waters with relatively high to intermediate ammonium, nitrate, and chl. a concentrations. Subgroups A_C1 and F_C1 were most abundant in coastal-transition zone waters with intermediate nutrient concentrations. E_O1 and G_O1 were most abundant at different depths of oligotrophic open-ocean waters (either in the upper mixed layer or just below). E_O1, A_C1, and F_C1 distributions differed from other narB subgroups and likely possess unique ecologies enabling them to be most abundant in waters between coastal and open-ocean waters. Different CCS zones possessed distinct Synechococcus communities. Core California current water possessed low numbers of narB subgroups relative to counted Synechococcus cells, and coastal-transition waters contained high abundances of Synechococcus cells and total number of narB subgroups. The presented biogeographic data provides insight on the distributions and ecologies of Synechococcus present in an eastern boundary current system.

Marine primary production in relation to climate variability and change

Marine photosynthetic plankton are responsible for approximately 50 petagrams (10(15)) of carbon per year of net primary production, an amount equivalent to that on land. This primary production supports essentially all life in the oceans and profoundly affects global biogeochemical cycles and climate. This review discusses the general distribution of primary production in the sea, the processes that regulate this distribution, and how marine primary production is sensitive to climate variability and change. Statistical modes of ocean variability and their characteristic interannual to multi-decadal timescales over the last century are described. Recent in situ and satellite time-series of primary production can be clearly linked to interannual ocean variability. Global marine primary production appears to have increased over the past several decades in association with multi-decadal variations. A paleoclimate record extends discussion to the centennial scale, providing contrasting insights into how marine primary production might vary in the future.

Biological consequences of el nino

Observations of the 1982-1983 El Niño make it possible to relate the anomalous ocean conditions to specific biological responses. In October 1982 upwelling ecosystems in the eastern equatorial Pacific began a series of transitions from the normal highly productive condition to greatly reduced productivity. The highly productive condition had returned by July 1983. Nutrients, phytoplankton biomass, and primary productivity are clearly regulated by the physical changes of El Niño. Evidence from 1982 and 1983 also suggests effects on higher organisms such as fish, seabirds, and marine mammals, but several more years of observation are required to accurately determine the magnitude of the consequences on these higher trophic levels.

Dimethylmercury in coastal upwelling waters, Monterey Bay, California

Depth profiles of dimethylmercury (DMHg) concentration were determined at nearshore to offshore sites in Monterey Bay, California. The onset of spring upwelling in the bay was accompanied by increases in DMHg concentrations. Profiles show DMHg increasing gradually with depth in fall and winter from <0.03 pM at the surface to 0.5 pM at 200 m. During the spring, DMHg concentrations increased between 30 and 100 m, first within Monterey Bay, then offshore. This change was accompanied by an increase in DMHg concentrations in the surface water DMHg between fall/winter (<0.03 pM) and spring (0.06-0.29 pM). Microbial activity associated with the remineralization of sinking organic matter produced by the high primary production in the bay may result in the relatively high DMHg in subsurface water in the bay, which when upwelled may facilitate the incorporation of organomercury into biota. As a result, productive coastal upwelling areas may represent an important source of methylated mercury to surface waters, and thus be an important source of mercury to marine ecosystems.

A massive phytoplankton bloom induced by an ecosystem-scale iron fertilization experiment in the equatorial Pacific Ocean

The seeding of an expanse of surface waters in the equatorial Pacific Ocean with low concentrations of dissolved iron triggered a massive phytoplankton bloom which consumed large quantities of carbon dioxide and nitrate that these microscopic plants cannot fully utilize under natural conditions. These and other observations provide unequivocal support for the hypothesis that phytoplankton growth in this oceanic region is limited by iron bioavailability.

Controls on tropical Pacific Ocean productivity revealed through nutrient stress diagnostics

In situ enrichment experiments have shown that the growth of bloom-forming diatoms in the major high-nitrate low-chlorophyll (HNLC) regions of the world's oceans is limited by the availability of iron. Yet even the largest of these manipulative experiments represents only a small fraction of an ocean basin, and the responses observed are strongly influenced by the proliferation of rare species rather than the growth of naturally dominant populations. Here we link unique fluorescence attributes of phytoplankton to specific physiological responses to nutrient stress, and use these relationships to evaluate the factors that constrain phytoplankton growth in the tropical Pacific Ocean on an unprecedented spatial scale. On the basis of fluorescence measurements taken over 12 years, we delineate three major ecophysiological regimes in this region. We find that iron has a key function in regulating phytoplankton growth in both HNLC and oligotrophic waters near the Equator and further south, whereas nitrogen and zooplankton grazing are the primary factors that regulate biomass production in the north. Application of our findings to the interpretation of satellite chlorophyll fields shows that productivity in the tropical Pacific basin may be 1.2-2.5 Pg C yr(-1) lower than previous estimates have suggested, a difference that is comparable to the global change in ocean production that accompanied the largest El Niño to La Niña transition on record.

Southern Ocean Iron Enrichment Experiment: Carbon Cycling in High- and Low-Si Waters

The availability of iron is known to exert a controlling influence on biological productivity in surface waters over large areas of the ocean and may have been an important factor in the variation of the concentration of atmospheric carbon dioxide over glacial cycles. The effect of iron in the Southern Ocean is particularly important because of its large area and abundant nitrate, yet iron-enhanced growth of phytoplankton may be differentially expressed between waters with high silicic acid in the south and low silicic acid in the north, where diatom growth may be limited by both silicic acid and iron. Two mesoscale experiments, designed to investigate the effects of iron enrichment in regions with high and low concentrations of silicic acid, were performed in the Southern Ocean. These experiments demonstrate iron's pivotal role in controlling carbon uptake and regulating atmospheric partial pressure of carbon dioxide.

From anchovies to sardines and back: multidecadal change in the Pacific Ocean

In the Pacific Ocean, air and ocean temperatures, atmospheric carbon dioxide, landings of anchovies and sardines, and the productivity of coastal and open ocean ecosystems have varied over periods of about 50 years. In the mid-1970s, the Pacific changed from a cool "anchovy regime" to a warm "sardine regime." A shift back to an anchovy regime occurred in the middle to late 1990s. These large-scale, naturally occurring variations must be taken into account when considering human-induced climate change and the management of ocean living resources.

Mortality of sea lions along the central California coast linked to a toxic diatom bloom

Over 400 California sea lions (Zalophus californianus) died and many others displayed signs of neurological dysfunction along the central California coast during May and June 1998. A bloom of Pseudo-nitzschia australis (diatom) was observed in the Monterey Bay region during the same period. This bloom was associated with production of domoic acid (DA), a neurotoxin that was also detected in planktivorous fish, including the northern anchovy (Engraulis mordax), and in sea lion body fluids. These and other concurrent observations demonstrate the trophic transfer of DA resulting in marine mammal mortality. In contrast to fish, blue mussels (Mytilus edulus) collected during the DA outbreak contained no DA or only trace amounts. Such findings reveal that monitoring of mussel toxicity alone does not necessarily provide adequate warning of DA entering the food web at levels sufficient to harm marine wildlife and perhaps humans.

Biological and chemical response of the equatorial pacific ocean to the 1997-98 El Nino

During the 1997-98 El Nino, the equatorial Pacific Ocean retained 0. 7 x 10(15) grams of carbon that normally would have been lost to the atmosphere as carbon dioxide. The surface ocean became impoverished in plant nutrients, and chlorophyll concentrations were the lowest on record. A dramatic recovery occurred in mid-1998, the system became highly productive, analogous to coastal environments, and carbon dioxide flux out of the ocean was again high. The spatial extent of the phytoplankton bloom that followed recovery from El Nino was the largest ever observed for the equatorial Pacific. These chemical and ecological perturbations were linked to changes in the upwelling of nutrient-enriched waters. The description and explanation of these dynamic changes would not have been possible without an observing system that combines biological, chemical, and physical sensors on moorings with remote sensing of chlorophyll.

Glacial to Interglacial Fluctuations in Productivity in the Equatorial Pacific as Indicated by Marine Barite

An empirical correlation between marine barite (BaSO4) accumulation rate in core-top sediment samples from two equatorial Pacific transects (at 140&deg;W and 110&deg;W) and the estimated primary productivity of the overlying water column were used to evaluate glacial to interglacial changes in productivity. Fluctuations in barite accumulation rates down-core indicate that during glacial periods of the past 450,000 years, the productivity in the central and eastern equatorial Pacific was about two times that during intervening interglacial periods. This result is consistent with other evidence that productivity was high in the eastern and central equatorial Pacific during the last glacial.

Methyl bromide: ocean sources, ocean sinks, and climate sensitivity

The oceans play an important role in the geochemical cycle of methyl bromide (CH3Br), the major carrier of O3-destroying bromine to the stratosphere. The quantity of CH3Br produced annually in seawater is comparable to the amount entering the atmosphere each year from natural and anthropogenic sources. The production mechanism is unknown but may be biological. Most of this CH3Br is consumed in situ by hydrolysis or reaction with chloride. The size of the fraction which escapes to the atmosphere is poorly constrained; measurements in seawater and the atmosphere have been used to justify both a large oceanic CH3Br flux to the atmosphere and a small net ocean sink. Since the consumption reactions are extremely temperature-sensitive, small temperature variations have large effects on the CH3Br concentration in seawater, and therefore on the exchange between the atmosphere and the ocean. The net CH3Br flux is also sensitive to variations in the rate of CH3Br production. We have quantified these effects using a simple steady state mass balance model. When CH3Br production rates are linearly scaled with seawater chlorophyll content, this model reproduces the latitudinal variations in marine CH3Br concentrations observed in the east Pacific Ocean by Singh et al. [1983] and by Lobert et al. [1995]. The apparent correlation of CH3Br production with primary production explains the discrepancies between the two observational studies, strengthening recent suggestions that the open ocean is a small net sink for atmospheric CH3Br, rather than a large net source. The Southern Ocean is implicated as a possible large net source of CH3Br to the atmosphere. Since our model indicates that both the direction and magnitude of CH3Br exchange between the atmosphere and ocean are extremely sensitive to temperature and marine productivity, and since the rate of CH3Br production in the oceans is comparable to the rate at which this compound is introduced to the atmosphere, even small perturbations to temperature or productivity can modify atmospheric CH3Br. Therefore atmospheric CH3Br should be sensitive to climate conditions. Our modeling indicates that climate-induced CH3Br variations can be larger than those resulting from small (+/- 25%) changes in the anthropogenic source, assuming that this source comprises less than half of all inputs. Future measurements of marine CH3Br, temperature, and primary production should be combined with such models to determine the relationship between marine biological activity and CH3Br production. Better understanding of the biological term is especially important to assess the importance of non-anthropogenic sources to stratospheric ozone loss and the sensitivity of these sources to global climate change.

Ribosomal DNA sequences discriminate among toxic and non-toxic Pseudonitzschia species

Cultured isolates of Pseudonitzschia australis Frenguelli, P. delicatissima (Cleve) Heiden, P. americana (Hasle) Fryxell, P. pungens (Grunow) Hasle, and P. pungens f. multiseries (Hasle) Hasle from Monterey Bay, California, were compared on the basis of their large-subunit ribosomal RNA gene (LsrDNA). Pseudonitzschia australis, P. pungens f. multiseries, and P. delicatissima were previously shown to produce the neurotoxin domoic acid; the remaining isolates are considered non-toxic. For each isolate approximately 800 base pairs of LsrDNA, encompassing both evolutionarily conserved and evolutionarily variable regions of the molecule, were amplified using the polymerase chain reaction (PCR) and sequenced. Phylogenetic trees generated by parsimony analysis of aligned sequences afford a preliminary view of the organisms genetic relationships. Species defined by morphological criteria are also distinguishable by LsrDNA sequence. Organisms known or suspected to produce domoic acid cluster at different termini on the phylogenetic tree. Two genetically distinct strains of P. australis and P. pungens were identified. Development of a restriction fragment length polymorphism (RFLP) assay of the LsrDNA is described. The RFLP assay discriminates each species, including distinguished strains of P. australis and P. pungens. The restriction test provides a rapid and convenient method for screening isolates' LsrDNA, facilitating further tests of the apparent positive correlation between Pseudonitzschia species' ribosomal gene signatures, morphology, and capacity to produce domoic acid.