Genotyping and Multivariate Regression Trees Reveal Ecological Diversification within the Microcystis aeruginosa Complex along a Wide Environmental Gradient

Abstract

Addressing the ecological and evolutionary processes underlying biodiversity patterns is essential to identify the mechanisms shaping community structure and function. In bacteria, the formation of new ecologically distinct populations (ecotypes) is proposed as one of the main drivers of diversification. New ecotypes arise when mutations in key functional genes or acquisition of new metabolic pathways by horizontal gene transfer allow the population to exploit new resources, permitting their coexistence with the parental population. We previously reported the presence of microcystin-producing organisms of the Microcystis aeruginosa complex (toxic MAC) through an 800-km environmental gradient ranging from freshwater to estuarine-marine waters in South America. We hypothesize that the success of toxic MAC in such a gradient is due to the existence of very closely related populations that are ecologically distinct (ecotypes), each specialized to a specific arrangement of environmental variables. Here, we analyzed toxic MAC genetic diversity through quantitative PCR (qPCR) and high-resolution melting analysis (HRMA) of a functional gene (mcyJ, microcystin synthetase cluster). We explored the variability of the mcyJ gene along the environmental gradient by multivariate classification and regression trees (mCART). Six groups of mcyJ genotypes were distinguished and associated with different combinations of water temperature, conductivity, and turbidity. We propose that each mcyJ variant associated with a defined environmental condition is an ecotype (or species) whose relative abundances vary according to their fitness in the local environment. This mechanism would explain the success of toxic MAC in such a wide array of environmental conditions. IMPORTANCE Organisms of the Microcystis aeruginosa complex form harmful algal blooms (HABs) in nutrient-rich water bodies worldwide. MAC HABs are difficult to manage owing to the production of potent toxins (microcystins) that resist water treatment. In addition, the role of microcystins in the ecology of MAC organisms is still elusive, meaning that the environmental conditions driving the toxicity of the bloom are not clear. Furthermore, the lack of coherence between morphology-based and genomic-based species classification makes it difficult to draw sound conclusions about when and where each member species of the MAC will dominate the bloom. Here, we propose that the diversification process and success of toxic MAC in a wide range of water bodies involves the generation of ecotypes, each specialized in a particular niche, whose relative abundance varies according to its fitness in the local environment. This knowledge can improve the generation of accurate prediction models of MAC growth and toxicity, helping to prevent human and animal intoxication.