Susanne Menden-Deuer

The Marine Microbial Eukaryote Transcriptome Sequencing Project (MMETSP): Illuminating the Functional Diversity of Eukaryotic Life in the Oceans through Transcriptome Sequencing.

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

Predator-Induced Fleeing Behaviors in Phytoplankton: A New Mechanism for Harmful Algal Bloom Formation?

In the plankton, heterotrophic microbes encounter and ingest phytoplankton prey, which effectively removes >50% of daily phytoplankton production in the ocean and influences global primary production and biochemical cycling rates. Factors such as size, shape, nutritional value, and presence of chemical deterrents are known to affect predation pressure. Effects of movement behaviors of either predator or prey on predation pressure, and particularly fleeing behaviors in phytoplankton are thus far unknown. Here, we quantified individual 3D movements, population distributions, and survival rates of the toxic phytoplankton species, Heterosigma akashiwo in response to a ciliate predator and predator-derived cues. We observed predator-induced defense behaviors previously unknown for phytoplankton. Modulation of individual phytoplankton movements during and after predator exposure resulted in an effective separation of predator and prey species. The strongest avoidance behaviors were observed when H. akashiwo co-occurred with an actively grazing predator. Predator-induced changes in phytoplankton movements resulted in a reduction in encounter rate and a 3-fold increase in net algal population growth rate. A spatially explicit population model predicted rapid phytoplankton bloom formation only when fleeing behaviors were incorporated. These model predictions reflected field observations of rapid H. akashiwo harmful algal bloom (HAB) formation in the coastal ocean. Our results document a novel behavior in phytoplankton that can significantly reduce predation pressure and suggests a new mechanism for HAB formation. Phytoplankton behaviors that minimize predatory losses, maximize resource acquisition, and alter community composition and distribution patterns could have major implications for our understanding and predictive capacity of marine primary production and biochemical cycling rates.

Broad Salinity Tolerance as a Refuge from Predation in the Harmful Raphidophyte Alga Heterosigma akashiwo (Raphidophyceae)

The ability of harmful algal species to form dense, nearly monospecific blooms remains an ecological and evolutionary puzzle. We hypothesized that predation interacts with estuarine salinity gradients to promote blooms of Heterosigma akashiwo (Y. Hada) Y. Hada ex Y. Hara et M. Chihara, a cosmopolitan toxic raphidophyte. Specifically, H. akashiwos broad salinity tolerance appears to provide a refuge from predation that enhances the net growth of H. akashiwo populations through several mechanisms. (1) Contrasting salinity tolerance of predators and prey. Estuarine H. akashiwo isolates from the west coast of North America grew rapidly at salinities as low as six, and distributed throughout experimental salinity gradients to salinities as low as three. In contrast, survival of most protistan predator species was restricted to salinities >15. (2) H. akashiwo physiological and behavioral plasticity. Acclimation to low salinity enhanced H. akashiwos ability to accumulate and grow in low salinity waters. In addition, the presence of a ciliate predator altered H. akashiwo swimming behavior, promoting accumulation in low‐salinity surface layers inhospitable to the ciliate. (3) Negative effects of low salinity on predation processes. Ciliate predation rates decreased sharply at salinities <25 and, for one species, H. akashiwo toxicity increased at low salinities. Taken together, these behaviors and responses imply that blooms can readily initiate in low salinity waters where H. akashiwo would experience decreased predation pressure while maintaining near‐maximal growth rates. The salinity structure of a typical estuary would provide this HAB species a unique refuge from predation. Broad salinity tolerance in raphidophytes may have evolved in part as a response to selective pressures associated with predation.

Many ways to stay in the game: Individual variability maintains high biodiversity in planktonic microorganisms

In apparent contradiction to competition theory, the number of known, coexisting plankton species far exceeds their explicable biodiversity—a discrepancy termed the Paradox of the Plankton. We introduce a new game-theoretic model for competing microorganisms in which one player consists of all organisms of one species. The stable points for the population dynamics in our model, known as strategic behaviour distributions (SBDs), are probability distributions of behaviours across all organisms which imply a stable population of the species as a whole. We find that intra-specific variability is the key characteristic that ultimately allows coexistence because the outcomes of competitions between individuals with variable competitive abilities are unpredictable. Our simulations based on the theoretical model show that up to 100 species can coexist for at least 10 000 generations, and that even small population sizes or species with inferior competitive ability can survive when there is intra-specific variability. In nature, this variability can be observed as niche differentiation, variability in environmental and ecological factors, and variability of individual behaviours or physiology. Therefore, previous specific explanations of the paradox are consistent with and provide specific examples of our suggestion that individual variability is the mechanism which solves the paradox.

Vibrio cholerae exploits sub-lethal concentrations of a competitor-produced antibiotic to avoid toxic interactions

Vibrio cholerae is a human pathogenic marine bacterium inhabiting coastal regions and is vectored into human food and water supplies via attachment to particles including detritus, phytoplankton, and zooplankton. Particle colonization by the pathogen is inhibited by an antagonistic interaction with the particle-associated Vibrionales bacterium SWAT3, a producer of the antibiotic andrimid. By analyzing the individual movement behaviors of V. cholerae exposed to a gradient of andrimid in a microfluidics device, we show that the pathogen has a concentration dependent avoidance response to sub-lethal concentrations of the pure antibiotic and to the metabolites produced by a growing colony of SWAT3-wild-type. This avoidance behavior includes a 25% increase in swimming speeds, 30% increase in run lengths, and a shift in the direction of the bacteria away from the andrimid source. Consequently, these behavioral shifts at low concentrations of andrimid would lead to higher diffusivity and result in the dispersion of bacteria away from the competitor and source of the antibiotic. Such alterations in motility were not elicited in response to a non-andrimid-producing SWAT3 mutant, suggesting andrimid may be a negative effector of chemotaxis for V. cholerae. The behavioral response of colonizing bacteria to sub-inhibitory concentrations of competitor-produced antibiotics is one mechanism that can influence microbial diversity and interspecific competition on particles, potentially affecting human health in coastal communities and element cycling in the ocean.

Going ballistic in the plankton: Anisotropic swimming behavior of marine protists

Diel vertical migrations (DVMs) of many plankton species, including single-celled protists, are well documented in the field and form a core component of many large-scale numerical models of plankton transport and ecology. However, the sparse quantitative data available describing motility behaviors of individual protists have frequently indicated that motility exhibits only short-term correlation on the order of a few seconds or hundreds of micrometers, resembling diffusive transport at larger scales—a result incompatible with DVM, which requires ballistic (straight-line) motion. We interrogated an extensive set of three-dimensional protistan movement trajectories in an effort to identify spatial and temporal correlation scales. Whereas the horizontal components of movement were diffusive, the vertical component remained highly correlated (i.e., nonrandom) for nearly all species for the duration of observation (up to 120 s and 6.1 mm) and in the absence of any environmental cues besides gravity. These persistent motility patterns may have been obscured in some previous studies due to the use of restrictive containers, dimensionally lumped, isotropic analyses, and/or an observation bias, inherent to observing free-swimming organisms with stationary cameras, which we accounted for in this study. Extrapolated over a 12-h period, conservative estimates of vertical travel ranges for the protists observed here would be 3‐10 m, while diffusive horizontal motion would result in about 10 cm of travel at most. Hence, these extended observations of phylogenetically diverse swimming protists, coupled with a quantitative analysis that accounts for anisotropy in the data, illustrate the small-scale mechanistic underpinnings of DVM.

Inherent High Correlation of Individual Motility Enhances Population Dispersal in a Heterotrophic, Planktonic Protist

Quantitative linkages between individual organism movements and the resulting population distributions are fundamental to understanding a wide range of ecological processes, including rates of reproduction, consumption, and mortality, as well as the spread of diseases and invasions. Typically, quantitative data are collected on either movement behaviors or population distributions, rarely both. This study combines empirical observations and model simulations to gain a mechanistic understanding and predictive ability of the linkages between both individual movement behaviors and population distributions of a single-celled planktonic herbivore. In the laboratory, microscopic 3D movements and macroscopic population distributions were simultaneously quantified in a 1L tank, using automated video- and image-analysis routines. The vertical velocity component of cell movements was extracted from the empirical data and used to motivate a series of correlated random walk models that predicted population distributions. Validation of the model predictions with empirical data was essential to distinguish amongst a number of theoretically plausible model formulations. All model predictions captured the essence of the population redistribution (mean upward drift) but only models assuming long correlation times (minute), captured the variance in population distribution. Models assuming correlation times of 8 minutes predicted the least deviation from the empirical observations. Autocorrelation analysis of the empirical data failed to identify a de-correlation time in the up to 30-second-long swimming trajectories. These minute-scale estimates are considerably greater than previous estimates of second-scale correlation times. Considerable cell-to-cell variation and behavioral heterogeneity were critical to these results. Strongly correlated random walkers were predicted to have significantly greater dispersal distances and more rapid encounters with remote targets (e.g. resource patches, predators) than weakly correlated random walkers. The tendency to disperse rapidly in the absence of aggregative stimuli has important ramifications for the ecology and biogeography of planktonic organisms that perform this kind of random walk.

pCO2 effects on species composition and growth of an estuarine phytoplankton community

The effects of ongoing changes in ocean carbonate chemistry on plankton ecology have important implications for food webs and biogeochemical cycling. However, conflicting results have emerged regarding species-specific responses to pCO2 enrichment and thus community responses have been difficult to predict. To assess community level effects (e.g., production) of altered carbonate chemistry, studies are needed that capitalize on the benefits of controlled experiments but also retain features of intact ecosystems that may exacerbate or ameliorate the effects observed in single-species or single cohort experiments. We performed incubations of natural plankton communities from Narragansett Bay, RI, USA in winter at ambient bay temperatures (5-13 °C), light and nutrient concentrations under three levels of controlled and constant CO2 concentrations, simulating past, present and future conditions at mean pCO2 levels of 224, 361, and 724 μatm respectively. Samples for carbonate analysis, chlorophyll a, plankton size-abundance, and plankton species composition were collected daily and phytoplankton growth rates in three different size fractions (<5, 5-20, and >20 μm) were measured at the end of the 7-day incubation period. Community composition changed during the incubation period with major increases in relative diatom abundance, which were similar across pCO2 treatments. At the end of the experiment, 24-hr growth responses to pCO2 levels varied as a function of cell size. The smallest size fraction (<5 μm) grew faster at the elevated pCO2 level. In contrast, the 5-20 μm size fraction grew fastest in the Present treatment and there were no significant differences in growth rate among treatments in the > 20 μm size fraction. Cell size distribution shifted toward smaller cells in both the Past and Future treatments but remained unchanged in the Present treatment. Similarity in Past and Future treatments for cell size distribution and growth rate (5-20 μm size fraction) illustrate non-monotonic effects of increasing pCO2 on ecological indicators and may be related to opposing physiological effects of high CO2 and low pH both within and among species. Interaction of these effects with other factors (e.g., nutrients, light, temperature, grazing, initial species composition) may explain variability among published studies. The absence of clear treatment-specific effects at the community level suggest that extrapolation of species-specific responses or experiments with only present day and future pCO2 treatments levels would produce misleading predictions of ocean acidification impacts on plankton production.

Persistent Intra-Specific Variation in Genetic and Behavioral Traits in the Raphidophyte, Heterosigma akashiwo

Motility is a key trait that phytoplankton utilize to navigate the heterogeneous marine environment. Quantifying both intra- and inter-specific variability in trait distributions is key to utilizing traits to distinguish groups of organisms and assess their ecological function. Because examinations of intra-specific variability are rare, here we measured three-dimensional movement behaviors and distribution patterns of seven genetically distinct strains of the ichthyotoxic raphidophyte, Heterosigma akashiwo. Strains were collected from different ocean basins but geographic distance between isolates was a poor predictor of genetic relatedness among strains. Observed behaviors were significantly different among all strains examined, with swimming speed and turning rate ranging from 33–115 μm s-1 and 41–110° s-1, respectively. Movement behaviors were consistent over at least 12 h, and in one case identical when measured several years apart. Movement behaviors were not associated with a specific cell size, carbon content, genetic relatedness, or geographic distance. These strain-specific behaviors resulted in algal populations that had distinct vertical distributions in the experimental tank. This study demonstrates that the traits of genetic identity and motility can provide resolution to distinguish strains of species, where variations in size or biomass are insufficient characteristics.

Thalassiosira spp. community composition shifts in response to chemical and physical forcing in the northeast Pacific Ocean.

Diatoms are genetically diverse unicellular photosynthetic eukaryotes that are key primary producers in the ocean. Many of the over 100 extant diatom species in the cosmopolitan genus Thalassiosira are difficult to distinguish in mixed populations using light microscopy. Here, we examine shifts in Thalassiosira spp. composition along a coastal to open ocean transect that encountered a 3-month-old Haida eddy in the northeast Pacific Ocean. To quantify shifts in Thalassiosira species composition, we developed a targeted automated ribosomal intergenic spacer analysis (ARISA) method to identify Thalassiosira spp. in environmental samples. As many specific fragment lengths are indicative of individual Thalassiosira spp., the ARISA method is a useful screening tool to identify changes in the relative abundance and distribution of specific species. The method also enabled us to assess changes in Thalassiosira community composition in response to chemical and physical forcing. Thalassiosira spp. community composition in the core of a 3-month-old Haida eddy remained largely (>80%) similar over a 2-week period, despite moving 24 km southwestward. Shifts in Thalassiosira species correlated with changes in dissolved iron (Fe) and temperature throughout the sampling period. Simultaneously tracking community composition and relative abundance of Thalassiosira species within the physical and chemical context they occurred allowed us to identify quantitative linkages between environmental conditions and community response.