Differential association of key bacterial groups with diatoms and Phaeocystis spp. during spring blooms in the Southern Ocean.

Interactions between phytoplankton and heterotrophic bacteria significantly influence the cycling of organic carbon in the ocean, with many of these interactions occurring at the micrometer scale. We explored potential associations between specific phytoplankton and bacteria in two size fractions, 0.8-3 µm and larger than 3 µm, at three naturally iron-fertilized stations and one high nutrient low chlorophyll station in the Southern Ocean. The composition of phytoplankton and bacterial communities was determined by sequencing the rbcL gene and 16S rRNA gene from DNA and RNA extracts, which represent presence and potential activity, respectively. Diatoms, particularly Thalassiosira, contributed significantly to the DNA sequences in the larger size fractions, while haptophytes were dominant in the smaller size fraction. Correlation analysis between the most abundant phytoplankton and bacterial operational taxonomic units revealed strong correlations between Phaeocystis and picoeukaryotes with SAR11, SAR116, Magnetospira, and Planktomarina. In contrast, most Thalassiosira operational taxonomic units showed the highest correlations with Polaribacter, Sulfitobacteria, Erythrobacter, and Sphingobium, while Fragilariopsis, Haslea, and Thalassionema were correlated with OM60, Fluviicola, and Ulvibacter. Our in-situ observations suggest distinct associations between phytoplankton and bacterial taxa, which could play crucial roles in nutrient cycling in the Southern Ocean.

Giant polyketide synthase enzymes in the biosynthesis of giant marine polyether toxins.

Prymnesium parvum are harmful haptophyte algae that cause massive environmental fish kills. Their polyketide polyether toxins, the prymnesins, are among the largest nonpolymeric compounds in nature and have biosynthetic origins that have remained enigmatic for more than 40 years. In this work, we report the "PKZILLAs," massive P. parvum polyketide synthase (PKS) genes that have evaded previous detection. PKZILLA-1 and -2 encode giant protein products of 4.7 and 3.2 megadaltons that have 140 and 99 enzyme domains. Their predicted polyene product matches the proposed pre-prymnesin precursor of the 90-carbon-backbone A-type prymnesins. We further characterize the variant PKZILLA-B1, which is responsible for the shorter B-type analog prymnesin-B1, from P. parvum RCC3426 and thus establish a general model of haptophyte polyether biosynthetic logic. This work expands expectations of genetic and enzymatic size limits in biology.

Phaeocystis globosa and diatom blooms promote distinct bacterial communities and associations in a coastal ecosystem.

Phytoplankton and bacteria form the foundation of marine food webs. While most studies on phytoplankton bloom influence on bacteria dynamics focus on diatom-dominated blooms due to their global ecological significance, it is unclear if similar patterns extend to other species that compete with diatoms like Phaeocystis spp. This study aimed to contribute to the understanding of associations between phytoplankton and bacteria in a temperate ecosystem. For this, we studied the dynamics of phytoplankton and bacteria, combining 16S metabarcoding, microscopy, and flow cytometry over 4 years (282 samples). Phytoplankton and bacterial communities were studied throughout the year, particularly during contrasting phytoplankton blooms dominated by the Haptophyte Phaeocystis globosa or diatoms. We applied extended local similarity analysis (eLSA) to construct networks during blooming and non-blooming periods. Overall, the importance of seasonal and species-specific interactions between phytoplankton and bacteria is highlighted. In winter, mixed diatom communities were interconnected with bacteria, indicating a synergistic degradation of diverse phytoplankton-derived substrates. In spring, despite the intensity variations of P. globosa blooms, the composition of bacterial communities remained consistent over several years, suggesting establishing a stable-state environment for bacterial communities. Specific associations between monospecific diatom blooms and bacteria were evidenced in summer.

Insights into the nature of ichthyotoxins from the Chrysochromulina leadbeateri blooms in Northern Norwegian fjords.

In May-June 2019, the microalga Chrysochromulina leadbeateri caused a massive fish-killing event in several fjords in Northern Norway, resulting in the largest direct impact ever on aquaculture in northern Europe due to toxic algae. Motivated by the fact that no algal toxins have previously been described from C. leadbeateri, we set out to investigate the chemical nature and toxicity of secondary metabolites in extracts of two strains (UIO 393, UIO 394) isolated from the 2019 bloom, as well as one older strain (UIO 035) isolated during a bloom in Northern Norway in 1991. Initial LC-DAD-MS/MS-based molecular networking analysis of the crude MeOH extracts of the cultivated strains showed that their profiles of small organic molecules, including a large number of known lipids, were very similar, suggesting that the same class of toxin(s) were likely the causative agents of the two harmful algal bloom (HAB) events. Next, bioassay-guided fractionation using the RTgill-W1 cell line and metabolomics analysis pointed to a major compound affording [M + H]+ ions at m/z 1399.8333 as a possible toxin, corresponding to a compound with the formula C67H127ClO27. Moreover, our study unveiled a series of minor analogues exhibiting distinct patterns of chlorination and sulfation, together defining a new family of compounds, which we propose to name leadbeaterins. Remarkably, these suspected toxins were detected in situ in samples collected during the 2019 bloom close to Tromsø, thereby consistent with a role in fish kills. The elemental compositions of the putative C. leadbeateri ichthyotoxins strongly indicate them to be long linear polyhydroxylated polyketides, structurally similar to sterolysins reported from a number of dinoflagellates.

Cellular morphological trait dataset for extant coccolithophores from the Atlantic Ocean.

Calcification and biomass production by planktonic marine organisms influences the global carbon cycle and fuels marine ecosystems. The major calcifying plankton group coccolithophores are highly diverse, comprising ca. 250-300 extant species. However, coccolithophore size (a key functional trait) and degree of calcification are poorly quantified, as most of our understanding of this group comes from a small number of species. We generated a novel reference dataset of coccolithophore morphological traits, including cell-specific data for coccosphere and cell size, coccolith size, number of coccoliths per cell, and cellular calcite content. This dataset includes observations from 1074 individual cells and represents 61 species from 25 genera spanning equatorial to temperate coccolithophore populations that were sampled during the Atlantic Meridional Transect (AMT) 14 cruise in 2004. This unique dataset can be used to explore relationships between morphological traits (cell size and cell calcite) and environmental conditions, investigate species-specific and community contributions to pelagic carbonate production, export and plankton biomass, and inform and validate coccolithophore representation in marine ecosystem and biogeochemical models.

Phagotrophy in the nitrogen-fixing haptophyte Braarudosphaera bigelowii.

Biological nitrogen fixation provides fixed nitrogen for microbes living in the oligotrophic open ocean. UCYN-A2, the previously known symbiont of Braarudosphaera bigelowii, now believed to be an early-stage B. bigelowii organelle that exchanges fixed nitrogen for fixed carbon, is globally distributed. Indirect evidence suggested that B. bigelowii might be a mixotrophic (phagotrophic) phototrophic flagellate. The goal of this study was to determine if B. bigelowii can graze on bacteria using several independent approaches. The results showed that B. bigelowii grazed on co-occurring bacteria at a rate of 5-7 cells/h/B. bigelowii and that the overall grazing rate was significantly higher at nighttime than at daytime. Bacterial abundance changes, assessed with 16S rRNA gene amplicon sequencing analysis, may have indicated preferential grazing by B. bigelowii on specific bacterial genotypes. In addition, Lysotracker™ staining of B. bigelowii suggested digestive activity inside B. bigelowii. Carbon and nitrogen fixation measurements revealed that the carbon demand of B. bigelowii could not be fulfilled by photosynthesis alone, implying supplementation by heterotrophy. These independent lines of evidence together revealed that B. bigelowii engages in phagotrophy, which, beyond serving as a supplementary source of carbon and energy, may also facilitate the indirect assimilation of inorganic nutrients.

Unpredicted ecosystem response to compound human impacts in a European river.

Climate change elevates the threat of compound heat and drought events, with their ecological and socioeconomic impacts exacerbated by human ecosystem alterations such as eutrophication, salinization, and river engineering. Here, we study how multiple stressors produced an environmental disaster in a large European river, the Oder River, where a toxic bloom of the brackish-water planktonic haptophyte Prymnesium parvum (the "golden algae") killed approximately 1000 metric tons of fish and most mussels and snails. We uncovered the complexity of this event using hydroclimatic data, remote sensing, cell counts, hydrochemical and toxin analyses, and genetics. After incubation in impounded upstream channels with drastically elevated concentrations of salts and nutrients, only a critical combination of chronic salt and nutrient pollution, acute high water temperatures, and low river discharge during a heatwave enabled the riverine mass proliferation of B-type P. parvum along a 500 km river section. The dramatic losses of large filter feeders and the spreading of vegetative cells and resting stages make the system more susceptible to new harmful algal blooms. Our findings show that global warming, water use intensification, and chronic ecosystem pollution could increase likelihood and severity of such compound ecoclimatic events, necessitating consideration in future impact models.

Algal methylated compounds shorten the lag phase of Phaeobacter inhibens bacteria.

The lag phase is key in resuming bacterial growth, but it remains underexplored particularly in environmental bacteria. Here we use transcriptomics and 13C-labelled metabolomics to show that the lag phase of the model marine bacterium Phaeobacter inhibens is shortened by methylated compounds produced by the microalgal partner, Emiliania huxleyi. Methylated compounds are abundantly produced and released by microalgae, and we show that their methyl groups can be collected by bacteria and assimilated through the methionine cycle. Our findings underscore the significance of methyl groups as a limiting factor during the lag phase and highlight the adjustability of this growth phase. In addition, we show that methylated compounds, typical of photosynthetic organisms, prompt diverse reductions in lag times in bacteria associated with algae and plants, potentially favouring early growth in some bacteria. These findings suggest ways to accelerate bacterial growth and underscore the significance of studying bacteria within an environmental context.

From genes to toxins: Profiling Prymnesium parvum during a riverine harmful algal bloom.

Blooms of Prymnesium parvum, a unicellular alga globally distributed in marine and brackish environments, frequently result in massive fish kills due to the production of toxins called prymnesins by this haptophyte. In August 2022, a harmful algal bloom (HAB) of this species occurred in the lower Oder River (Poland and Germany), which caused mass mortalities of fish and other organisms. This HAB was linked to low discharge of the Oder and mining activities that caused a significant increase in salinity. In this context, we report on the molecular detection and screening of this haptophyte and its toxins in environmental samples and clonal cultures derived thereof. Both conventional PCR and droplet digital PCR assays reliably detected P. parvum in environmental samples. eDNA metabarcoding using the V4 region of the 18S rRNA gene revealed a single Prymnesium sequence variant, but failed to identify it to species level. Four clonal cultures established from environmental samples were unambiguously identified as P. parvum by molecular phylogenetics (near full-length 18S rRNA gene) and light microscopy. Phylogenetic analysis (ITS1-5.8S-ITS2 marker region) placed the cultured phylotype within a clade containing other P. parvum strains known to produce B-type prymnesins. Toxin-screening of the cultures using liquid chromatography-electrospray ionization - time of flight mass spectrometry identified B-type prymnesins, which were also detected in extracts of filter residues from water samples of the Oder collected during the HAB. Overall, our investigation provides a detailed characterization of P. parvum, including their prymnesins, during this HAB in the Oder River, contributing valuable insights into this ecological disaster. In addition, the droplet digital PCR assay established here will be useful for future monitoring of low levels of P. parvum on the Oder River or any other salt-impacted and brackish water bodies.

A rapid and multi-endpoint ecotoxicological test using Mychonastes afer for efficient screening of metals and herbicides.

Microalgal growth-based tests are international standards for ecotoxicity assessment; however, their long exposure times, large sample volumes, and reliance on a single growth-endpoint make them inadequate for rapid toxicity screening. Here, we aimed to develop a rapid and simple ecotoxicological test using the fast-growing green alga Mychonastes afer, with multiple endpoints-growth, lipid content, and photosynthesis. We exposed M. afer to two metals-silver and copper-and two herbicides-atrazine and diuron-for 24 h and identified the most sensitive and reliable endpoints for each toxicant: the maximum electron transport rate (ETRmax) for Ag, Cu and atrazine, and the lipid content for diuron. Lipid content was found to be both a sensitive and reliable biomarker, meeting the effluent limit guidelines in both the Republic of Korea and the USA. The sensitivity of M. afer to Ag and atrazine also closely matched the HC5 values derived from the species sensitivity distribution approach, confirming its reliability for setting regulatory concentrations of these contaminants. Our calculated predicted no-effect concentration (PNEC) values were similar to established European Union PNECs for Ag, Cu, atrazine, and diuron, underlining the utility of these biological endpoints for ecological risk assessment and regulatory decision making. This method required lower sample volume (2 mL vs 100 mL) and exposure time (24 h vs 72-120 h) than conventional green algal tests, and eliminated the need for labour-intensive cell counting, expensive equipment, and chlorophyll fluorescence measurement expertise. Overall, this M. afer test can be a valuable tool for the rapid screening of wastewater for metals and herbicides, contributing to environmental protection and management practices.

Hydrogen isotope fractionation is controlled by CO2 in coccolithophore lipids.

Hydrogen isotope ratios (δ2H) represent an important natural tracer of metabolic processes, but quantitative models of processes controlling H-fractionation in aquatic photosynthetic organisms are lacking. Here, we elucidate the underlying physiological controls of 2H/1H fractionation in algal lipids by systematically manipulating temperature, light, and CO2(aq) in continuous cultures of the haptophyte Gephyrocapsa oceanica. We analyze the hydrogen isotope fractionation in alkenones (αalkenone), a class of acyl lipids specific to this species and other haptophyte algae. We find a strong decrease in the αalkenone with increasing CO2(aq) and confirm αalkenone correlates with temperature and light. Based on the known biosynthesis pathways, we develop a cellular model of the δ2H of algal acyl lipids to evaluate processes contributing to these controls on fractionation. Simulations show that longer residence times of NADPH in the chloroplast favor a greater exchange of NADPH with 2H-richer intracellular water, increasing αalkenone. Higher chloroplast CO2(aq) and temperature shorten NADPH residence time by enhancing the carbon fixation and lipid synthesis rates. The inverse correlation of αalkenone to CO2(aq) in our cultures suggests that carbon concentrating mechanisms (CCM) do not achieve a constant saturation of CO2 at the Rubisco site, but rather that chloroplast CO2 varies with external CO2(aq). The pervasive inverse correlation of αalkenone with CO2(aq) in the modern and preindustrial ocean also suggests that natural populations may not attain a constant saturation of Rubisco with the CCM. Rather than reconstructing growth water, αalkenone may be a powerful tool to elucidate the carbon limitation of photosynthesis.

Dynamic patterns of carbohydrate metabolism genes in bacterioplankton during marine algal blooms.

Carbohydrates play a pivotal role in nutrient recycling and regulation of algal-bacterial interactions. Despite their ecological significance, the intricate molecular mechanisms governing regulation of phycosphere carbohydrates by bacterial taxa linked with natural algal bloom have yet to be fully elucidated. Here, a comprehensive temporal metagenomic analysis was conducted to explore the carbohydrate-active enzyme (CAZyme) genes in two discrete algal bloom microorganisms (Gymnodinium catenatum and Phaeocystis globosa) across three distinct bloom stages: pre-bloom, peak bloom, and post-bloom. Elevated levels of extracellular carbohydrates, primarily rhamnose, galactose, glucose, and arabinose, were observed during the initial and post-peak stages. The prominent CAZyme families identified-glycoside hydrolases (GH) and carbohydrate-binding modules (CBMs)-were present in both algal bloom occurrences. In the G. catenatum bloom, GH23/24 and CBM13/14 were prevalent during the pre-bloom and peak bloom stages, whereas GH2/3/30 and CBM12/24 exhibited increased prevalence during the post-bloom phase. In contrast, the P. globosa bloom had a dominance of GH13/23 and CBM19 in the initial phase, and this was succeeded by GH3/19/24/30 and CBM54 in the later stages. This gene pool variation-observed distinctly in specific genera-highlighted the dynamic structural shifts in functional resources driven by temporal alterations in available substrates. Additionally, ecological linkage analysis underscored a correlation between carbohydrates (or their related genes) and phycospheric bacteria, hinting at a pattern of bottom-up control. These findings contribute to understanding of the dynamic nature of CAZymes, emphasizing the substantial influence of substrate availability on the metabolic capabilities of algal symbiotic bacteria, especially in terms of carbohydrates.

The impact of artificial light pollution at night on the life history parameters of rotifer Brachionus plicatilis with different food experiences.

Artificial light at night (ALAN) may pose threat to rotifer Brachionus plicatilis. Additionally, the food of rotifer, i.e. algal community composition, often fluctuates. Thus, we selected five wavelengths of ALAN (purple, blue, green, red, white) and a three-colored light flashing mode (3-Flash) to test their impacts on life history traits of B. plicatilis with different food experiences, including those feeding Chlorella vulgaris (RC) or Phaeocystis globosa (RP). Results indicated purple ALAN promoted RC development, white ALAN inhibited RC development, while 3-Flash and white ALAN promoted RP development. Under red and white ALAN, RP increased fecundity but decreased lifespan. High-quality food enhanced rotifer's resistance to the impact of ALAN on lifespan. ALAN and food experience interacted on B. plicatilis. The effect of blue ALAN has less negative effects on B. plicatilis, based on hierarchical cluster analysis. Such findings are helpful to evaluate the potential impact of ALAN on marine zooplankton.

Dynamics and scale variations of blooms of Phaeocystis globosa "giant colony" ecotype and key environmental factors in the Beibu Gulf, China.

The Beibu Gulf has experienced blooms of Phaeocystis globosa "giant colony" ecotype (PGGCE), with noticeable variations in bloom scale across years. However, driving environmental factors and their roles remain poorly understood. In this study, we quantified dynamics of PGGCE cells in 2016-2017 and 2018-2019, and analyzed their correlations with environment factors. The results revealed that PGGCE blooms primarily occurred in Guangxi coast and western waters of Leizhou Peninsula during winter months, exhibiting distinct developmental processes. Bloom intensity, duration, and distribution differed significantly between two bloom events. In 2016-2017, peak PGGCE density exceeded 2.0 × 105 cells L-1 nearly double that of 2018-2019. Furthermore, bloom sustained five months during 2016-2017, compared to three months during 2018-2019. Prolonged period of low temperatures and elevated nitrate concentrations favored PGGCE growth and colony formation, resulting in a larger scale bloom during winter 2016 as opposed to winter 2018.

Diversity and seasonal occurrence of haptophyta in northern South China Sea through size-fractionated metabarcoding.

Haptophyta plays a key role in marine pico-nanoeukaryote communities but information on their diversity and ecology is extremely limited. A total of 103 water samples were collected in northern South China Sea to assess the diversity of haptophyta through metabarcoding targeting 18S V4 rDNA. Furthermore, we investigated the potential genetic differentiation among seasonal occurring Phaeocystis globosa using the high resolution molecular marker pgcp1. 18S V4 rDNA metabarcoding dataset revealed 41 species of haptophytes, with 16 of them as the first record in this region. Notably, six harmful species were detected, including Chrysochromulina leadbeateri, Phaeocystis globosa, and Prymnesium parvum. The pgcp1 marker revealed two clades of Phaeocystis globosa and both of them were present around the year. Clade I was found to predominate in warm season, while Clade III tended to bloom in cold waters. Our results highlight the risk potential of harmful haptophytes in the northern South China Sea.

The micro-ecological feature of colonies is a potential strategy for Phaeocystis globosa bloom formation.

Phaeocystis globosa is among the dominant microalgae associated with harmful algal blooms. P. globosa has a polymorphic life cycle and its ecological success has been attributed to algal colony formation, however, few studies have assessed differences in microbial communities and their functional profiles between intra- and extra-colonies during P. globosa blooms. To address this, environmental and metagenomics tools were used to conduct a time-series analysis of the bacterial composition and metabolic characteristics of intra- and extra-colonies during a natural P. globosa bloom. The results show that bacterial composition, biodiversity, and network interactions differed significantly between intra- and extra-colonies. Dominant extra-colonial bacteria were Bacteroidia and Saccharimonadis, while dominant intra-colonial bacteria included Alphaproteobacteria and Gammaproteobacteria. Despite the lower richness and diversity observed in the intra-colonial bacterial community, relative to extra-colonies, the complexity and interconnectedness of the intra-colonial networks were higher. Regarding bacterial function, more functional genes were enriched in substance metabolism (polysaccharides, iron element and dimethylsulfoniopropionate) and signal communication (quorum sensing, indoleacetic acid-IAA) pathways in intra- than in extra-colonies. Conceptual model construction showed that microbial cooperative synthesis of ammonium, vitamin B12, IAA, and siderophores were strongly related to the P. globosa bloom, particularly in the intra-colonial environment. Overall, our data highlight the differences in bacterial structure and functions within and outside the colony during P. globosa blooms. These findings represent fundamental information indicating that phenotypic heterogeneity is a selective strategy that improves microbial population competitiveness and environmental adaptation, benefiting P. globosa bloom formation and persistence.

Impact of environmental factors changes induced by marine heatwaves and heavy precipitation on antibiotic toxicity to Isochrysis galbana: Implications for climate change adaptation.

Isochrysis galbana, a crucial primary producer and food source in aquatic ecosystems, faces increasing challenges from climate change and emerging contaminants like antibiotics. This study investigates the combined effects of sudden temperature increase (representing marine heatwaves) and rapid salinity change (representing extreme precipitation events) on the toxicity of tetracycline (TC) and oxytetracycline (OTC) to I. galbana. Short-term experiments reveal heightened antibiotic toxicity at 31 °C or salinities of 18 PSU, surpassing algal tolerance limits. Long-term tests show decreased inhibition of algal growth on day 9, indicating algal adaptation to the environment. Analyses of photosynthesis II efficiency, pigment content, and macromolecular composition support this, suggesting adaptation mechanism activation. While algae acclimate to the environment during long-term antibiotic exposure, extreme weather conditions may compromise this adaptation. These findings have implications for managing antibiotics in aquatic environments under climate change.

Enantioselective transformation of phytoplankton-derived dihydroxypropanesulfonate by marine bacteria.

Chirality, a fundamental property of matter, is often overlooked in the studies of marine organic matter cycles. Dihydroxypropanesulfonate (DHPS), a globally abundant organosulfur compound, serves as an ecologically important currency for nutrient and energy transfer from phytoplankton to bacteria in the ocean. However, the chirality of DHPS in nature and its transformation remain unclear. Here, we developed a novel approach using chiral phosphorus-reagent labeling to separate DHPS enantiomers. Our findings demonstrated that at least one enantiomer of DHPS is present in marine diatoms and coccolithophores, and that both enantiomers are widespread in marine environments. A novel chiral-selective DHPS catabolic pathway was identified in marine Roseobacteraceae strains, where HpsO and HpsP dehydrogenases at the gateway to DHPS catabolism act specifically on R-DHPS and S-DHPS, respectively. R-DHPS is also a substrate for the dehydrogenase HpsN. All three dehydrogenases generate stable hydrogen bonds between the chirality-center hydroxyls of DHPS and highly conserved residues, and HpsP also form coordinate-covalent bonds between the chirality-center hydroxyls and Zn2+, which determines the mechanistic basis of strict stereoselectivity. We further illustrated the role of enzymatic promiscuity in the evolution of DHPS metabolism in Roseobacteraceae and SAR11. This study provides the first evidence of chirality's involvement in phytoplankton-bacteria metabolic currencies, opening a new avenue for understanding the ocean organosulfur cycle.

Effects of micro- and nano-plastics on growth, antioxidant system, DMS, and DMSP production in Emiliania huxleyi.

Due to the potential impacts of microplastics (MPs) and nanoplastics (NPs) on algal growth and thereby affect the climate-relevant substances, dimethylsulfoniopropionate (DMSP) and dimethyl sulfide (DMS), we studied the polystyrene (PS) MPs and NPs of 1 µm and 80 nm impacts on the growth, chlorophyll content, reactive oxygen species (ROS), antioxidant enzyme activity, and DMS/DMSP production in Emiliania huxleyi. E. huxleyi is a prominent oceanic alga that plays a key role in DMS and DMSP production. The results revealed that high concentrations of MPs and NPs inhibited the growth, carotenoid (Car), and Chl a concentrations of E. huxleyi. However, short-time exposure to low concentrations of PS MPs and NPs stimulated the growth of E. huxleyi. Furthermore, high concentrations of MPs and NPs resulted in an increase in the superoxide anion radical (O2.-) production rate and a decrease in the malondialdehyde (MDA) content compared with the low concentrations. Exposure to MPs and NPs at 5 mg L-1 induced superoxide dismutase (SOD) activity as a response to scavenging ROS. High concentrations of MPs and NPs significantly inhibited the production of DMSP and DMS. The findings of this study support the potential ecotoxicological impacts of MPs and NPs on algal growth, antioxidant system, and dimethylated sulfur compounds production, which maybe potentially impact the global climate.

Modelling phytoplankton-virus interactions: phytoplankton blooms and lytic virus transmission.

A dynamic reaction-diffusion model of four variables is proposed to describe the spread of lytic viruses among phytoplankton in a poorly mixed aquatic environment. The basic ecological reproductive index for phytoplankton invasion and the basic reproduction number for virus transmission are derived to characterize the phytoplankton growth and virus transmission dynamics. The theoretical and numerical results from the model show that the spread of lytic viruses effectively controls phytoplankton blooms. This validates the observations and experimental results of Emiliana huxleyi-lytic virus interactions. The studies also indicate that the lytic virus transmission cannot occur in a low-light or oligotrophic aquatic environment.