Widespread use of proton-pumping rhodopsin in Antarctic phytoplankton.

Photosynthetic carbon (C) fixation by phytoplankton in the Southern Ocean (SO) plays a critical role in regulating air-sea exchange of carbon dioxide and thus global climate. In the SO, photosynthesis (PS) is often constrained by low iron, low temperatures, and low but highly variable light intensities. Recently, proton-pumping rhodopsins (PPRs) were identified in marine phytoplankton, providing an alternate iron-free, light-driven source of cellular energy. These proteins pump protons across cellular membranes through light absorption by the chromophore retinal, and the resulting pH energy gradient can then be used for active membrane transport or for synthesis of adenosine triphosphate. Here, we show that PPR is pervasive in Antarctic phytoplankton, especially in iron-limited regions. In a model SO diatom, we found that it was localized to the vacuolar membrane, making the vacuole a putative alternative phototrophic organelle for light-driven production of cellular energy. Unlike photosynthetic C fixation, which decreases substantially at colder temperatures, the proton transport activity of PPR was unaffected by decreasing temperature. Cellular PPR levels in cultured SO diatoms increased with decreasing iron concentrations and energy production from PPR photochemistry could substantially augment that of PS, especially under high light intensities, where PS is often photoinhibited. PPR gene expression and high retinal concentrations in phytoplankton in SO waters support its widespread use in polar environments. PPRs are an important adaptation of SO phytoplankton to growth and survival in their cold, iron-limited, and variable light environment.

A surplus no more? Variation in krill availability impacts reproductive rates of Antarctic baleen whales.

The krill surplus hypothesis of unlimited prey resources available for Antarctic predators due to commercial whaling in the 20th century has remained largely untested since the 1970s. Rapid warming of the Western Antarctic Peninsula (WAP) over the past 50 years has resulted in decreased seasonal ice cover and a reduction of krill. The latter is being exacerbated by a commercial krill fishery in the region. Despite this, humpback whale populations have increased but may be at a threshold for growth based on these human-induced changes. Understanding how climate-mediated variation in prey availability influences humpback whale population dynamics is critical for focused management and conservation actions. Using an 8-year dataset (2013-2020), we show that inter-annual humpback whale pregnancy rates, as determined from skin-blubber biopsy samples (n = 616), are positively correlated with krill availability and fluctuations in ice cover in the previous year. Pregnancy rates showed significant inter-annual variability, between 29% and 86%. Our results indicate that krill availability is in fact limiting and affecting reproductive rates, in contrast to the krill surplus hypothesis. This suggests that this population of humpback whales may be at a threshold for population growth due to prey limitations. As a result, continued warming and increased fishing along the WAP, which continue to reduce krill stocks, will likely impact this humpback whale population and other krill predators in the region. Humpback whales are sentinel species of ecosystem health, and changes in pregnancy rates can provide quantifiable signals of the impact of environmental change at the population level. Our findings must be considered paramount in developing new and more restrictive conservation and management plans for the Antarctic marine ecosystem and minimizing the negative impacts of human activities in the region.

Photosynthetic energy conversion efficiency in the West Antarctic Peninsula.

The West Antarctic Peninsula (WAP) is a highly productive polar ecosystem where phytoplankton dynamics are regulated by intense bottom-up control from light and iron availability. Rapid climate change along the WAP is driving shifts in the mixed layer depth and iron availability. Elucidating the relative role of each of these controls and their interactions is crucial for understanding of how primary productivity will change in coming decades. Using a combination of ultra-high-resolution variable chlorophyll fluorescence together with fluorescence lifetime analyses on the 2017 Palmer Long Term Ecological Research cruise, we mapped the temporal and spatial variability in phytoplankton photophysiology across the WAP. Highest photosynthetic energy conversion efficiencies and lowest fluorescence quantum yields were observed in iron replete coastal regions. Photosynthetic energy conversion efficiencies decreased by ~ 60% with a proportional increase in quantum yields of thermal dissipation and fluorescence on the outer continental shelf and slope. The combined analysis of variable fluorescence and lifetimes revealed that, in addition to the decrease in the fraction of inactive reaction centers, up to 20% of light harvesting chlorophyll-protein antenna complexes were energetically uncoupled from photosystem II reaction centers in iron-limited phytoplankton. These biophysical signatures strongly suggest severe iron limitation of photosynthesis in the surface waters along the continental slope of the WAP.

Towards Quantitative Microbiome Community Profiling Using Internal Standards.

An inherent issue in high-throughput rRNA gene tag sequencing microbiome surveys is that they provide compositional data in relative abundances. This often leads to spurious correlations, making the interpretation of relationships to biogeochemical rates challenging. To overcome this issue, we quantitatively estimated the abundance of microorganisms by spiking in known amounts of internal DNA standards. Using a 3-year sample set of diverse microbial communities from the Western Antarctica Peninsula, we demonstrated that the internal standard method yielded community profiles and taxon cooccurrence patterns substantially different from those derived using relative abundances. We found that the method provided results consistent with the traditional CHEMTAX analysis of pigments and total bacterial counts by flow cytometry. Using the internal standard method, we also showed that chloroplast 16S rRNA gene data in microbial surveys can be used to estimate abundances of certain eukaryotic phototrophs such as cryptophytes and diatoms. In Phaeocystis, scatter in the 16S/18S rRNA gene ratio may be explained by physiological adaptation to environmental conditions. We conclude that the internal standard method, when applied to rRNA gene microbial community profiling, is quantitative and that its application will substantially improve our understanding of microbial ecosystems.IMPORTANCE High-throughput-sequencing-based marine microbiome profiling is rapidly expanding and changing how we study the oceans. Although powerful, the technique is not fully quantitative; it provides taxon counts only in relative abundances. In order to address this issue, we present a method to quantitatively estimate microbial abundances per unit volume of seawater filtered by spiking known amounts of internal DNA standards into each sample. We validated this method by comparing the calculated abundances to other independent estimates, including chemical markers (pigments) and total bacterial cell counts by flow cytometry. The internal standard approach allows us to quantitatively estimate and compare marine microbial community profiles, with important implications for linking environmental microbiomes to quantitative processes such as metabolic and biogeochemical rates.

Changes in the upper ocean mixed layer and phytoplankton productivity along the West Antarctic Peninsula.

The West Antarctic Peninsula (WAP) has experienced significant change over the last 50 years. Using a 24 year spatial time series collected by the Palmer Long Term Ecological Research programme, we assessed long-term patterns in the sea ice, upper mixed layer depth (MLD) and phytoplankton productivity. The number of sea ice days steadily declined from the 1980s until a recent reversal that began in 2008. Results show regional differences between the northern and southern regions sampled during regional ship surveys conducted each austral summer. In the southern WAP, upper ocean MLD has shallowed by a factor of 2. Associated with the shallower mixed layer is enhanced phytoplankton carbon fixation. In the north, significant interannual variability resulted in the mixed layer showing no trended change over time and there was no significant increase in the phytoplankton productivity. Associated with the recent increases in sea ice there has been an increase in the photosynthetic efficiency (chlorophyll a-normalized carbon fixation) in the northern and southern regions of the WAP. We hypothesize the increase in sea ice results in increased micronutrient delivery to the continental shelf which in turn leads to enhanced photosynthetic performance.This article is part of the theme issue 'The marine system of the West Antarctic Peninsula: status and strategy for progress in a region of rapid change'.

Spring-summer net community production, new production, particle export and related water column biogeochemical processes in the marginal sea ice zone of the Western Antarctic Peninsula 2012-2014.

New production (New P, the rate of net primary production (NPP) supported by exogenously supplied limiting nutrients) and net community production (NCP, gross primary production not consumed by community respiration) are closely related but mechanistically distinct processes. They set the carbon balance in the upper ocean and define an upper limit for export from the system. The relationships, relative magnitudes and variability of New P (from 15NO3- uptake), O2 : argon-based NCP and sinking particle export (based on the 238U : 234Th disequilibrium) are increasingly well documented but still not clearly understood. This is especially true in remote regions such as polar marginal ice zones. Here we present a 3-year dataset of simultaneous measurements made at approximately 50 stations along the Western Antarctic Peninsula (WAP) continental shelf in midsummer (January) 2012-2014. Net seasonal-scale changes in water column inventories (0-150 m) of nitrate and iodide were also estimated at the same stations. The average daily rates based on inventory changes exceeded the shorter-term rate measurements. A major uncertainty in the relative magnitude of the inventory estimates is specifying the start of the growing season following sea-ice retreat. New P and NCP(O2) did not differ significantly. New P and NCP(O2) were significantly greater than sinking particle export from thorium-234. We suggest this is a persistent and systematic imbalance and that other processes such as vertical mixing and advection of suspended particles are important export pathways.This article is part of the theme issue 'The marine system of the west Antarctic Peninsula: status and strategy for progress in a region of rapid change'.

Inter-decadal variability of phytoplankton biomass along the coastal West Antarctic Peninsula.

The West Antarctic Peninsula (WAP) is a climatically sensitive region where periods of strong warming have caused significant changes in the marine ecosystem and food-web processes. Tight coupling between phytoplankton and higher trophic levels implies that the coastal WAP is a bottom-up controlled system, where changes in phytoplankton dynamics may largely impact other food-web components. Here, we analysed the inter-decadal time series of year-round chlorophyll-a (Chl) collected from three stations along the coastal WAP: Carlini Station at Potter Cove (PC) on King George Island, Palmer Station on Anvers Island and Rothera Station on Adelaide Island. There were trends towards increased phytoplankton biomass at Carlini Station (PC) and Palmer Station, while phytoplankton biomass declined significantly at Rothera Station over the studied period. The impacts of two relevant climate modes to the WAP, the El Niño-Southern Oscillation and the Southern Annular Mode, on winter and spring phytoplankton biomass appear to be different among the three sampling stations, suggesting an important role of local-scale forcing than large-scale forcing on phytoplankton dynamics at each station. The inter-annual variability of seasonal bloom progression derived from considering all three stations together captured ecologically meaningful, seasonally co-occurring bloom patterns which were primarily constrained by water-column stability strength. Our findings highlight a coupled link between phytoplankton and physical and climate dynamics along the coastal WAP, which may improve our understanding of overall WAP food-web responses to climate change and variability.This article is part of the theme issue 'The marine system of the West Antarctic Peninsula: status and strategy for progress in a region of rapid change'.

Molecular physiology of Antarctic diatom natural assemblages and bloom event reveal insights into strategies contributing to their ecological success.

The continental shelf of the Western Antarctic Peninsula (WAP) is a highly variable system characterized by strong cross-shelf gradients, rapid regional change, and large blooms of phytoplankton, notably diatoms. Rapid environmental changes coincide with shifts in plankton community composition and productivity, food web dynamics, and biogeochemistry. Despite the progress in identifying important environmental factors influencing plankton community composition in the WAP, the molecular basis for their survival in this oceanic region, as well as variations in species abundance, metabolism, and distribution, remains largely unresolved. Across a gradient of physicochemical parameters, we analyzed the metabolic profiles of phytoplankton as assessed through metatranscriptomic sequencing. Distinct phytoplankton communities and metabolisms closely mirrored the strong gradients in oceanographic parameters that existed from coastal to offshore regions. Diatoms were abundant in coastal, southern regions, where colder and fresher waters were conducive to a bloom of the centric diatom, Actinocyclus. Members of this genus invested heavily in growth and energy production; carbohydrate, amino acid, and nucleotide biosynthesis pathways; and coping with oxidative stress, resulting in uniquely expressed metabolic profiles compared to other diatoms. We observed strong molecular evidence for iron limitation in shelf and slope regions of the WAP, where diatoms in these regions employed iron-starvation induced proteins, a geranylgeranyl reductase, aquaporins, and urease, among other strategies, while limiting the use of iron-containing proteins. The metatranscriptomic survey performed here reveals functional differences in diatom communities and provides further insight into the environmental factors influencing the growth of diatoms and their predicted response to changes in ocean conditions. IMPORTANCE: In the Southern Ocean, phytoplankton must cope with harsh environmental conditions such as low light and growth-limiting concentrations of the micronutrient iron. Using metratranscriptomics, we assessed the influence of oceanographic variables on the diversity of the phytoplankton community composition and on the metabolic strategies of diatoms along the Western Antarctic Peninsula, a region undergoing rapid climate change. We found that cross-shelf differences in oceanographic parameters such as temperature and variable nutrient concentrations account for most of the differences in phytoplankton community composition and metabolism. We opportunistically characterized the metabolic underpinnings of a large bloom of the centric diatom Actinocyclus in coastal waters of the WAP. Our results indicate that physicochemical differences from onshore to offshore are stronger than between southern and northern regions of the WAP; however, these trends could change in the future, resulting in poleward shifts in functional differences in diatom communities and phytoplankton blooms.

Depth drives the distribution of microbial ecological functions in the coastal western Antarctic Peninsula.

The Antarctic marine environment is a dynamic ecosystem where microorganisms play an important role in key biogeochemical cycles. Despite the role that microbes play in this ecosystem, little is known about the genetic and metabolic diversity of Antarctic marine microbes. In this study we leveraged DNA samples collected by the Palmer Long Term Ecological Research (LTER) project to sequence shotgun metagenomes of 48 key samples collected across the marine ecosystem of the western Antarctic Peninsula (wAP). We developed an in silico metagenomics pipeline (iMAGine) for processing metagenomic data and constructing metagenome-assembled genomes (MAGs), identifying a diverse genomic repertoire related to the carbon, sulfur, and nitrogen cycles. A novel analytical approach based on gene coverage was used to understand the differences in microbial community functions across depth and region. Our results showed that microbial community functions were partitioned based on depth. Bacterial members harbored diverse genes for carbohydrate transformation, indicating the availability of processes to convert complex carbons into simpler bioavailable forms. We generated 137 dereplicated MAGs giving us a new perspective on the role of prokaryotes in the coastal wAP. In particular, the presence of mixotrophic prokaryotes capable of autotrophic and heterotrophic lifestyles indicated a metabolically flexible community, which we hypothesize enables survival under rapidly changing conditions. Overall, the study identified key microbial community functions and created a valuable sequence library collection for future Antarctic genomics research.

Photosynthetic energy storage efficiency in Chlamydomonas reinhardtii, based on microsecond photoacoustics.

Using a novel, pulsed micro-second time-resolved photoacoustic (PA) instrument, we measured thermal dissipation and energy storage (ES) in the intact cells of wild type (WT) Chlamydomonas reinhardtii, and mutants lacking either PSI or PSII reaction centers (RCs). On this time scale, the kinetic contributions of the thermal expansion component due to heat dissipation of absorbed energy and the negative volume change due to electrostriction induced by charge separation in each of the photosystems could be readily distinguished. Kinetic analysis revealed that PSI and PSII RCs exhibit strikingly different PA signals where PSI is characterized by a strong electrostriction signal and a weak thermal expansion component while PSII has a small electrostriction component and large thermal expansion. The calculated ES efficiencies at ~10 µs were estimated to be 80 ± 5 and 50 ± 13% for PSII-deficient mutants and PSI-deficient mutants, respectively, and 67 ± 2% for WT. The overall ES efficiency was positively correlated with the ratio of PSI to PSI + PSII. Our results suggest that the shallow excitonic trap in PSII limits the efficiency of ES as a result of an evolutionary frozen metabolic framework of two photosystems in all oxygenic photoautotrophs.

The role of nutricline depth in regulating the ocean carbon cycle.

Carbon uptake by marine phytoplankton, and its export as organic matter to the ocean interior (i.e., the "biological pump"), lowers the partial pressure of carbon dioxide (pCO(2)) in the upper ocean and facilitates the diffusive drawdown of atmospheric CO(2). Conversely, precipitation of calcium carbonate by marine planktonic calcifiers such as coccolithophorids increases pCO(2) and promotes its outgassing (i.e., the "alkalinity pump"). Over the past approximately 100 million years, these two carbon fluxes have been modulated by the relative abundance of diatoms and coccolithophores, resulting in biological feedback on atmospheric CO(2) and Earth's climate; yet, the processes determining the relative distribution of these two phytoplankton taxa remain poorly understood. We analyzed phytoplankton community composition in the Atlantic Ocean and show that the distribution of diatoms and coccolithophorids is correlated with the nutricline depth, a proxy of nutrient supply to the upper mixed layer of the ocean. Using this analysis in conjunction with a coupled atmosphere-ocean intermediate complexity model, we predict a dramatic reduction in the nutrient supply to the euphotic layer in the coming century as a result of increased thermal stratification. Our findings indicate that, by altering phytoplankton community composition, this causal relationship may lead to a decreased efficiency of the biological pump in sequestering atmospheric CO(2), implying a positive feedback in the climate system. These results provide a mechanistic basis for understanding the connection between upper ocean dynamics, the calcium carbonate-to-organic C production ratio and atmospheric pCO(2) variations on time scales ranging from seasonal cycles to geological transitions.

Decline in plankton diversity and carbon flux with reduced sea ice extent along the Western Antarctic Peninsula.

Since the middle of the past century, the Western Antarctic Peninsula has warmed rapidly with a significant loss of sea ice but the impacts on plankton biodiversity and carbon cycling remain an open question. Here, using a 5-year dataset of eukaryotic plankton DNA metabarcoding, we assess changes in biodiversity and net community production in this region. Our results show that sea-ice extent is a dominant factor influencing eukaryotic plankton community composition, biodiversity, and net community production. Species richness and evenness decline with an increase in sea surface temperature (SST). In regions with low SST and shallow mixed layers, the community was dominated by a diverse assemblage of diatoms and dinoflagellates. Conversely, less diverse plankton assemblages were observed in waters with higher SST and/or deep mixed layers when sea ice extent was lower. A genetic programming machine-learning model explained up to 80% of the net community production variability at the Western Antarctic Peninsula. Among the biological explanatory variables, the sea-ice environment associated plankton assemblage is the best predictor of net community production. We conclude that eukaryotic plankton diversity and carbon cycling at the Western Antarctic Peninsula are strongly linked to sea-ice conditions.

Interpopulational differences in the nutritional condition of Aequiyoldia eightsii (Protobranchia: Nuculanidae) from the Western Antarctic Peninsula during austral summer.

The Western Antarctic Peninsula (WAP) is a hotspot for environmental change and has a strong environmental gradient from North to South. Here, for the first time we used adult individuals of the bivalve Aequiyoldia eightsii to evaluate large-scale spatial variation in the biochemical composition (measured as lipid, protein and fatty acids) and energy content, as a proxy for nutritional condition, of three populations along the WAP: O'Higgins Research Station in the north (63.3°S), Yelcho Research Station in mid-WAP (64.9°S) and Rothera Research Station further south (67.6°S). The results reveal significantly higher quantities of lipids (L), proteins (P), energy (E) and total fatty acids (FA) in the northern population (O'Higgins) (L: 8.33 ± 1.32%; P: 22.34 ± 3.16%; E: 171.53 ± 17.70 Joules; FA: 16.33 ± 0.98 mg g) than in the mid-WAP population (Yelcho) (L: 6.23 ± 0.84%; P: 18.63 ± 1.17%; E: 136.67 ± 7.08 Joules; FA: 10.93 ± 0.63 mg g) and southern population (Rothera) (L: 4.60 ± 0.51%; P: 13.11 ± 0.98%; E: 98.37 ± 5.67 Joules; FA: 7.58 ± 0.48 mg g). We hypothesize these differences in the nutritional condition could be related to a number of biological and environmental characteristics. Our results can be interpreted as a consequence of differences in phenology at each location; differences in somatic and gametogenic growth rhythms. Contrasting environmental conditions throughout the WAP such as seawater temperature, quantity and quality of food from both planktonic and sediment sources, likely have an effect on the metabolism and nutritional intake of this species.

The evolutionary inheritance of elemental stoichiometry in marine phytoplankton.

Phytoplankton is a nineteenth century ecological construct for a biologically diverse group of pelagic photoautotrophs that share common metabolic functions but not evolutionary histories. In contrast to terrestrial plants, a major schism occurred in the evolution of the eukaryotic phytoplankton that gave rise to two major plastid superfamilies. The green superfamily appropriated chlorophyll b, whereas the red superfamily uses chlorophyll c as an accessory photosynthetic pigment. Fossil evidence suggests that the green superfamily dominated Palaeozoic oceans. However, after the end-Permian extinction, members of the red superfamily rose to ecological prominence. The processes responsible for this shift are obscure. Here we present an analysis of major nutrients and trace elements in 15 species of marine phytoplankton from the two superfamilies. Our results indicate that there are systematic phylogenetic differences in the two plastid types where macronutrient (carbon:nitrogen:phosphorus) stoichiometries primarily reflect ancestral pre-symbiotic host cell phenotypes, but trace element composition reflects differences in the acquired plastids. The compositional differences between the two plastid superfamilies suggest that changes in ocean redox state strongly influenced the evolution and selection of eukaryotic phytoplankton since the Proterozoic era.

The role of functional traits and trade-offs in structuring phytoplankton communities: scaling from cellular to ecosystem level.

Trait-based approaches to community structure are increasingly used in terrestrial ecology. We show that such an approach, augmented by a mechanistic analysis of trade-offs among functional traits, can be successfully used to explain community composition of marine phytoplankton along environmental gradients. Our analysis of literature on major functional traits in phytoplankton, such as parameters of nutrient-dependent growth and uptake, reveals physiological trade-offs in species abilities to acquire and utilize resources. These trade-offs, arising from fundamental relations such as cellular scaling laws and enzyme kinetics, define contrasting ecological strategies of nutrient acquisition. Major groups of marine eukaryotic phytoplankton have adopted distinct strategies with associated traits. These diverse strategies of nutrient utilization can explain the distribution patterns of major functional groups and size classes along nutrient availability gradients.

Rapid shelf-wide cooling response of a stratified coastal ocean to hurricanes.

Large uncertainty in the predicted intensity of tropical cyclones (TCs) persists compared to the steadily improving skill in the predicted TC tracks. This intensity uncertainty has its most significant implications in the coastal zone, where TC impacts to populated shorelines are greatest. Recent studies have demonstrated that rapid ahead-of-eye-center cooling of a stratified coastal ocean can have a significant impact on hurricane intensity forecasts. Using observation-validated, high-resolution ocean modeling, the stratified coastal ocean cooling processes observed in two U.S. Mid-Atlantic hurricanes were investigated: Hurricane Irene (2011)-with an inshore Mid-Atlantic Bight (MAB) track during the late summer stratified coastal ocean season-and Tropical Storm Barry (2007)-with an offshore track during early summer. For both storms, the critical ahead-of-eye-center depth-averaged force balance across the entire MAB shelf included an onshore wind stress balanced by an offshore pressure gradient. This resulted in onshore surface currents opposing offshore bottom currents that enhanced surface to bottom current shear and turbulent mixing across the thermocline, resulting in the rapid cooling of the surface layer ahead-of-eye-center. Because the same baroclinic and mixing processes occurred for two storms on opposite ends of the track and seasonal stratification envelope, the response appears robust. It will be critical to forecast these processes and their implications for a wide range of future storms using realistic 3-D coupled atmosphere-ocean models to lower the uncertainty in predictions of TC intensities and impacts and enable coastal populations to better respond to increasing rapid intensification threats in an era of rising sea levels.

Seasonal time bombs: dominant temperate viruses affect Southern Ocean microbial dynamics.

Rapid warming in the highly productive western Antarctic Peninsula (WAP) region of the Southern Ocean has affected multiple trophic levels, yet viral influences on microbial processes and ecosystem function remain understudied in the Southern Ocean. Here we use cultivation-independent quantitative ecological and metagenomic assays, combined with new comparative bioinformatic techniques, to investigate double-stranded DNA viruses during the WAP spring-summer transition. This study demonstrates that (i) temperate viruses dominate this region, switching from lysogeny to lytic replication as bacterial production increases, and (ii) Southern Ocean viral assemblages are genetically distinct from lower-latitude assemblages, primarily driven by this temperate viral dominance. This new information suggests fundamentally different virus-host interactions in polar environments, where intense seasonal changes in bacterial production select for temperate viruses because of increased fitness imparted by the ability to switch replication strategies in response to resource availability. Further, temperate viral dominance may provide mechanisms (for example, bacterial mortality resulting from prophage induction) that help explain observed temporal delays between, and lower ratios of, bacterial and primary production in polar versus lower-latitude marine ecosystems. Together these results suggest that temperate virus-host interactions are critical to predicting changes in microbial dynamics brought on by warming in polar marine systems.

Adaptive evolution of phytoplankton cell size.

We present a simple nutrient-phytoplankton-zooplankton (NPZ) model that incorporates adaptive evolution and allometric relations to examine the patterns and consequences of adaptive changes in plankton body size. Assuming stable environmental conditions, the model makes the following predictions. First, phytoplankton should evolve toward small sizes typical of picoplankton. Second, in the absence of grazers, nutrient concentration is minimized as phytoplankton reach their fitness maximum. Third, increasing nutrient flux tends to increase phytoplankton cell size in the presence of phytoplankton-zooplankton coevolution but has no effect in the absence of zooplankton. Fourth, phytoplankton reach their fitness maximum in the absence of grazers, and the evolutionary nutrient-phytoplankton system has a stable equilibrium. In contrast, phytoplankton may approach their fitness minimum in the evolutionary NPZ system where phytoplankton and zooplankton are allowed to coevolve, which may result in oscillatory (unstable) dynamics of the evolutionary NPZ system, compared with the otherwise stable nonevolutionary NPZ system. These results suggest that evolutionary interactions between phytoplankton and zooplankton may have contributed to observed changes in phytoplankton sizes and associated biogeochemical cycles over geological time scales.