Angel M. Segura

Emergent neutrality drives phytoplankton species coexistence

The mechanisms that drive species coexistence and community dynamics have long puzzled ecologists. Here, we explain species coexistence, size structure and diversity patterns in a phytoplankton community using a combination of four fundamental factors: organism traits, size-based constraints, hydrology and species competition. Using a ‘microscopic’ Lotka–Volterra competition (MLVC) model (i.e. with explicit recipes to compute its parameters), we provide a mechanistic explanation of species coexistence along a niche axis (i.e. organismic volume). We based our model on empirically measured quantities, minimal ecological assumptions and stochastic processes. In nature, we found aggregated patterns of species biovolume (i.e. clumps) along the volume axis and a peak in species richness. Both patterns were reproduced by the MLVC model. Observed clumps corresponded to niche zones (volumes) where species fitness was highest, or where fitness was equal among competing species. The latter implies the action of equalizing processes, which would suggest emergent neutrality as a plausible mechanism to explain community patterns.

Competition Drives Clumpy Species Coexistence in Estuarine Phytoplankton

Understanding the mechanisms that maintain biodiversity is a fundamental problem in ecology. Competition is thought to reduce diversity, but hundreds of microbial aquatic primary producers species coexist and compete for a few essential resources (e.g., nutrients and light). Here, we show that resource competition is a plausible mechanism for explaining clumpy distribution on individual species volume (a proxy for the niche) of estuarine phytoplankton communities ranging from North America to South America and Europe, supporting the Emergent Neutrality hypothesis. Furthermore, such a clumpy distribution was also observed throughout the Holocene in diatoms from a sediment core. A Lotka-Volterra competition model predicted position in the niche axis and functional affiliation of dominant species within and among clumps. Results support the coexistence of functionally equivalent species in ecosystems and indicate that resource competition may be a key process to shape the size structure of estuarine phytoplankton, which in turn drives ecosystem functioning.

Dynamics of toxic genotypes of Microcystis aeruginosa complex (MAC) through a wide freshwater to marine environmental gradient

Bloom-forming species belonging to Microcystis aeruginosa complex (MAC) are the most commonly reported worldwide. MAC blooms are composed by toxic and non-toxic genotypes and the environmental conditions favouring the dominance of toxic genotypes are still a matter of debate among the scientific community. In this study, we evaluated the distribution of toxic MAC genotypes along a seasonal cycle and over an environmental gradient spanning 800km, from a eutrophic freshwater reservoir in Río Uruguay to marine water in the outer limit of Río de la Plata. Abundance of four mcy genes, mcyB, mcyD, mcyE and mcyJ was determined by qPCR and used as a proxy of abundance of toxic MAC genotypes. All the mcy genes were detected through the seasonal cycle at all sampling sites, being systematically higher in the freshwater reservoir and decreasing towards the marine site. The highest toxic genotype abundance was found during the austral summer months. According to generalized linear regressions and random forest models, temperature and conductivity were the most relevant explanatory variables. This suggests that although toxic MAC genotypes grow optimally in freshwater, they are also able to tolerate the high-salinity and low temperature conditions found in estuarine and marine waters. This ability to resist harsh conditions impose a health risk and a management challenge. To our knowledge, this is the first report addressing several mcy genes in a broad gradient that includes a wide array of different environmental conditions.

Morphology-based differences in the thermal response of freshwater phytoplankton

The thermal response of maximum growth rate in morphology-based functional groups (MBFG) of freshwater phytoplankton is analysed. Contrasting an exponential Boltzmann–Arrhenius with a unimodal model, three main features were evaluated: (i) the activation energy of the rise (Er), (ii) the presence of a break in the thermal response and (iii) the activation energy of the fall (Ef). The whole dataset (N = 563) showed an exponential increase (Er ∼ 0.5), a break around 24°C and no temperature dependence after the breakpoint (Ef = 0). Contrasting thermal responses among MBFG were found. All groups showed positive activation energy (Er > 0), four showed no evidence of decline in growth rate (temperature range = 0–35°C) and two presented a breakpoint followed by a sharp decrease in growth rate. Our results evidenced systematic differences between MBFG in the thermal response and a coherent response significantly related to morphological traits other than size (i.e. within MBFG). These results provide relevant information for water quality modelling and climate change predictions.

A multilevel trait-based approach to the ecological performance of Microcystis aeruginosa complex from headwaters to the ocean

The Microcystis aeruginosa complex (MAC) clusters cosmopolitan and conspicuous harmful bloom-forming cyanobacteria able to produce cyanotoxins. It is hypothesized that low temperatures and brackish salinities are the main barriers to MAC proliferation. Here, patterns at multiple levels of organization irrespective of taxonomic identity (i.e. a trait-based approach) were analyzed. MAC responses from the intracellular (e.g. respiratory activity) to the ecosystem level (e.g. blooms) were evaluated in wide environmental gradients. Experimental results on buoyancy and respiratory activity in response to increased salinity (0-35) and a literature review of maximum growth rates under different temperatures and salinities were combined with field sampling from headwaters (800km upstream) to the marine end of the Rio de la Plata estuary (Uruguay-South America). Salinity and temperature were the major variables affecting MAC responses. Experimentally, freshwater MAC cells remained active for 24h in brackish waters (salinity=15) while colonies increased their flotation velocity. At the population level, maximum growth rate decreased with salinity and presented a unimodal exponential response with temperature, showing an optimum at 27.5°C and a rapid decrease thereafter. At the community and ecosystem levels, MAC occurred from fresh to marine waters (salinity 30) with a sustained relative increase of large mucilaginous colonies biovolume with respect to individual cells. Similarly, total biomass and, specific and morphological richness decreased with salinity while blooms were only detected in freshwater both at high (33°C) and low (11°C) temperatures. In brackish waters, large mucilaginous colonies presented advantages under osmotic restrictive conditions. These traits values have also been associated with higher toxicity potential. This suggest salinity or low temperatures would not represent effective barriers for the survival and transport of potentially toxic MAC under likely near future scenarios of increasing human impacts (i.e. eutrophication, dam construction and climate change).

Improved biovolume estimation of Microcystis aeruginosa colonies: A statistical approach

The Microcystis aeruginosa complex (MAC) clusters many of the most common freshwater and brackish bloom-forming cyanobacteria. In monitoring protocols, biovolume estimation is a common approach to determine MAC colonies biomass and useful for prediction purposes. Biovolume (μm3 mL-1) is calculated multiplying organism abundance (orgL-1) by colonial volume (μm3org-1). Colonial volume is estimated based on geometric shapes and requires accurate measurements of dimensions using optical microscopy. A trade-off between easy-to-measure but low-accuracy simple shapes (e.g. sphere) and time costly but high-accuracy complex shapes (e.g. ellipsoid) volume estimation is posed. Overestimations effects in ecological studies and management decisions associated to harmful blooms are significant due to the large sizes of MAC colonies. In this work, we aimed to increase the precision of MAC biovolume estimations by developing a statistical model based on two easy-to-measure dimensions. We analyzed field data from a wide environmental gradient (800 km) spanning freshwater to estuarine and seawater. We measured length, width and depth from ca. 5700 colonies under an inverted microscope and estimated colonial volume using three different recommended geometrical shapes (sphere, prolate spheroid and ellipsoid). Because of the non-spherical shape of MAC the ellipsoid resulted in the most accurate approximation, whereas the sphere overestimated colonial volume (3-80) especially for large colonies (MLD higher than 300 μm). Ellipsoid requires measuring three dimensions and is time-consuming. Therefore, we constructed different statistical models to predict organisms depth based on length and width. Splitting the data into training (2/3) and test (1/3) sets, all models resulted in low training (1.41-1.44%) and testing average error (1.3-2.0%). The models were also evaluated using three other independent datasets. The multiple linear model was finally selected to calculate MAC volume as an ellipsoid based on length and width. This work contributes to achieve a better estimation of MAC volume applicable to monitoring programs as well as to ecological research.