Selection on life-history traits and genetic population divergence in rotifers

Intraspecific variation in sensitivity to food availability and temperature-induced phenotypic plasticity in the rotifer Keratella cochlearis

Organisms with wide environmentally induced morphological plasticity and cosmopolitan distribution, e.g. the common freshwater rotifer Keratella cochlearis, are ideal models to study the evolution of plastic polymorphisms and the capacity of zooplankton to adapt to local selection conditions. We investigated population-level differences (population-by-environment interaction) in sensitivity to food availability and temperature-induced phenotypic plasticity between two clones of K. cochlearis isolated from neighboring populations in Ruidera Natural Park (Spain) with different trophic statuses: Tinaja lake (mesotrophic) and Cueva Morenilla lake (eutrophic). Using common-garden experiments, each clone proved to have a different sensitivity to food availability, with substantial phenotypic differences between them. When rotifers grew at moderate temperature (15.6°C), low food levels were more efficiently used by the Tinaja versus Cueva Morenilla clone, whereas high food levels were more efficiently used by the Cueva Morenilla versus Tinaja clone. The posterior spine was much longer and the lorica wider in the Tinaja versus Cueva Morenilla clone, with no difference in lorica length. Phylogenetic analysis based on cytochrome c oxidase subunit I (COI) gene sequences showed that the two populations have the same haplotype. This is the first study to show possible local adaptation by a rotifer species to habitats that consistently differ in food availability. We also detected an intriguing deviation from the expected negative relationship between posterior spine length and temperature. Our experimental results indicate that intermediate temperatures may activate the gene responsible for spine elongation in K. cochlearis This suggests that rotifers in nature could use water temperature as proxy signal of a change in predation risk before defense is needed.

Ecological genomics of adaptation to unpredictability in experimental rotifer populations

Elucidating the genetic basis of phenotypic variation in response to different environments is key to understanding how populations evolve. Facultatively sexual rotifers can develop adaptive responses to fluctuating environments. In a previous evolution experiment, diapause-related traits changed rapidly in response to two selective regimes (predictable vs unpredictable) in laboratory populations of the rotifer Brachionus plicatilis. Here, we investigate the genomic basis of adaptation to environmental unpredictability in these experimental populations. We identified and genotyped genome-wide polymorphisms in 169 clones from both selective regimes after seven cycles of selection using genotyping by sequencing (GBS). Additionally, we used GBS data from the 270 field clones from which the laboratory populations were established. This GBS dataset was used to identify candidate SNPs under selection. A total of 76 SNPs showed divergent selection, three of which are candidates for being under selection in the particular unpredictable fluctuation pattern studied. Most of the remaining SNPs showed strong signals of adaptation to laboratory conditions. Furthermore, a genotype-phenotype association approach revealed five SNPs associated with two key life-history traits in the adaptation to unpredictability. Our results contribute to elucidating the genomic basis for adaptation to unpredictable environments and lay the groundwork for future evolution studies in rotifers.

Within species expressed genetic variability and gene expression response to different temperatures in the rotifer Brachionus calyciflorus sensu stricto

Genetic divergence is impacted by many factors, including phylogenetic history, gene flow, genetic drift, and divergent selection. Rotifers are an important component of aquatic ecosystems, and genetic variation is essential to their ongoing adaptive diversification and local adaptation. In addition to coding sequence divergence, variation in gene expression may relate to variable heat tolerance, and can impose ecological barriers within species. Temperature plays a significant role in aquatic ecosystems by affecting species abundance, spatio-temporal distribution, and habitat colonization. Recently described (formerly cryptic) species of the Brachionus calyciflorus complex exhibit different temperature tolerance both in natural and in laboratory studies, and show that B. calyciflorus sensu stricto (s.s.) is a thermotolerant species. Even within B. calyciflorus s.s., there is a tendency for further temperature specializations. Comparison of expressed genes allows us to assess the impact of stressors on both expression and sequence divergence among disparate populations within a single species. Here, we have used RNA-seq to explore expressed genetic diversity in B. calyciflorus s.s. in two mitochondrial DNA lineages with different phylogenetic histories and differences in thermotolerance. We identify a suite of candidate genes that may underlie local adaptation, with a particular focus on the response to sustained high or low temperatures. We do not find adaptive divergence in established candidate genes for thermal adaptation. Rather, we detect divergent selection among our two lineages in genes related to metabolism (lipid metabolism, metabolism of xenobiotics).

Founder effects drive the genetic structure of passively dispersed aquatic invertebrates

Populations of passively dispersed organisms in continental aquatic habitats typically show high levels of neutral genetic differentiation despite their high dispersal capabilities. Several evolutionary factors, including founder events, local adaptation, and life cycle features such as high population growth rates and the presence of propagule banks, have been proposed to be responsible for this paradox. Here, we have modeled the colonization process to assess the impact of migration rate, population growth rate, population size, local adaptation and life-cycle features on the population genetic structure in these organisms. Our simulations show that the strongest effect on population structure are persistent founder effects, resulting from the interaction of a few population founders, high population growth rates, large population sizes and the presence of diapausing egg banks. In contrast, the role of local adaptation, genetic hitchhiking and migration is limited to small populations in these organisms. Our results indicate that local adaptation could have different impact on genetic structure in different groups of zooplankters.

Genomic signatures of local adaptation to the degree of environmental predictability in rotifers

Environmental fluctuations are ubiquitous and thus essential for the study of adaptation. Despite this, genome evolution in response to environmental fluctuations —and more specifically to the degree of environmental predictability– is still unknown. Saline lakes in the Mediterranean region are remarkably diverse in their ecological conditions, which can lead to divergent local adaptation patterns in the inhabiting aquatic organisms. The facultatively sexual rotifer Brachionus plicatilis shows diverging local adaptation in its life-history traits in relation to estimated environmental predictability in its habitats. Here, we used an integrative approach —combining environmental, phenotypic and genomic data for the same populations– to understand the genomic basis of this diverging adaptation. Firstly, a novel draft genome for B. plicatilis was assembled. Then, genome-wide polymorphisms were studied using genotyping by sequencing on 270 clones from nine populations in eastern Spain. As a result, 4,543 high-quality SNPs were identified and genotyped. More than 90 SNPs were found to be putatively under selection with signatures of diversifying and balancing selection. Over 140 SNPs were correlated with environmental or phenotypic variables revealing signatures of local adaptation, including environmental predictability. Putative functions were associated to most of these SNPs, since they were located within annotated genes. Our results reveal associations between genomic variation and the degree of environmental predictability, providing genomic evidence of adaptation to local conditions in natural rotifer populations.

Adaptation in response to environmental unpredictability

Understanding how organisms adaptively respond to environmental fluctuations is a fundamental question in evolutionary biology. The Mediterranean region typically exhibits levels of environmental unpredictability that vary greatly in habitats over small geographical scales. In cyclically parthenogenetic rotifers, clonal proliferation occurs along with occasional bouts of sex. These bouts contribute to the production of diapausing eggs, which allows survival between growing seasons. Here, we studied two diapause-related traits in rotifers using clones from nine Brachionus plicatilis natural populations that vary in the degree of environmental unpredictability. We tested the hypothesis that the level of environmental unpredictability is directly related to the propensity for sex and inversely related to the hatching fraction of diapausing eggs. We found significant levels of genetic variation within populations for both traits. Interestingly, a positive correlation between pond unpredictability-quantified in a previous study from satellite imagery-and the propensity for sex was found. This correlation suggests a conservative, bet-hedging strategy that provides protection against unexpectedly short growing seasons. By contrast, the hatching fraction of diapausing eggs was not related to the level of environmental predictability. Our results highlight the ability of rotifer populations to locally adapt to time-varying environments, providing an evolutionarily relevant step forward in relating life-history traits to a quantitative measure of environmental unpredictability.

Diapause and maintenance of facultative sexual reproductive strategies

Facultative sex combines sexual and asexual reproduction in the same individual (or clone) and allows for a large diversity of life-history patterns regarding the timing, frequency and intensity of sexual episodes. In addition, other life-history traits such as a diapause stage may become linked to sex. Here, we develop a matrix modelling framework for addressing the cost of sex in facultative sexuals, in constant, periodic and stochastically fluctuating environments. The model is parametrized using life-history data from Brachionus calyciflorus, a facultative sexual rotifer in which sex and diapause are linked. Sexual propensity was an important driver of costs in constant environments, in which high costs (always > onefold, and sometimes > twofold) indicated that asexuals should outcompete facultative sexuals. By contrast, stochastic environments with high temporal autocorrelation favoured facultative sex over obligate asex, in particular, if the penalty to fecundity in 'bad' environments was large. In such environments, obligate asexuals were constrained by their life cycle length (i.e. time from birth to last reproductive adult age class), which determined an upper limit to the number of consecutive bad periods they could tolerate. Nevertheless, when facultative asexuals with different sexual propensities competed simultaneously against each other and asex, the lowest sex propensity was the most successful in stochastic environments with positive autocorrelation. Our results suggest that a highly specific mechanism (i.e. diapause linked to sex) can alone stabilize facultative sex in these animals, and protect it from invasion of both asexual and pure sexual strategies.This article is part of the themed issue 'Weird sex: the underappreciated diversity of sexual reproduction'.

Competitive advantage and higher fitness in native populations of genetically structured planktonic diatoms

It has been shown that the planktonic diatom Skeletonema from neighbouring areas are genetically differentiated despite absence of physical dispersal barriers. We revisited two sites, Mariager Fjord and Kattegat, NE Atlantic, and isolated new strains. Microsatellite genotyping and F-statistics revealed that the populations were genetically differentiated. An experiment was designed to investigate if populations are locally adapted and have a native competitive advantage. Ten strains from each location were grown individually in native and foreign water to investigate differences in produced biomass. Additionally, we mixed six pairs, one strain from each site, and let them grow together in native and foreign water. Strains from Mariager Fjord and Kattegat produced higher biomass in native water. In the competition experiment, strains from both sites displayed higher relative abundance and demonstrated competitive advantage in their native water. The cause of the differentiated growth is unknown, but could possibly be attributed to differences in silica concentration or viruses in the two water types. Our data show that dispersal potential does not influence the genetic structure of the populations. We conclude that genetic adaptation has not been overruled by gene flow, but instead the responses to different selection conditions are enforcing the observed genetic structure.

From local adaptation to ecological speciation in copepod populations from neighboring lakes

Continental copepods have been derived from several independent invasive events from the sea, but the subsequent evolutionary processes that account for the current diversity in lacustrine environments are virtually unknown. Salinity is highly variable among lakes and constitutes a source of divergent selection driving potential reproductive isolation. We studied four populations of the calanoid copepod Leptodiaptomus cf. sicilis inhabiting four neighboring lakes with a common history (since the Late Pleistocene) located in the Oriental Basin, Mexico; one lake is shallow and varies in salinity periodically (1.4-10 g L(-1)), while three are deep and permanent, with constant salinity (0.5, 1.1 and 6.5 g L(-1), respectively). We hypothesized that (1) these populations belong to a different species than L. sicilis sensu stricto and (2) are experiencing ecologically based divergence due to salinity differences. We assessed morphological and molecular (mtDNA) COI variation, as well as fitness differences and tests of reproductive isolation. Although relationships of the Mexican populations with L. sicilis s.s. could not be elucidated, we identified a clear pattern of divergent selection driven by salinity conditions. The four populations can still be considered a single biological species (sexual recognition and hybridization are still possible in laboratory conditions), but they have diverged into at least three different phenotypes: two locally adapted, specialized in the lakes of constant salinity (saline vs. freshwater), and an intermediate generalist phenotype inhabiting the temporary lake with fluctuating salinity. The specialized phenotypes are poorly suited as migrants, so prezygotic isolation due to immigrant inviability is highly probable. This implication was supported by molecular evidence that showed restricted gene flow, persistence of founder events, and a pattern of allopatric fragmentation. This study showed how ecologically based divergent selection may explain diversification patterns in lacustrine copepods.

Molecular and quantitative genetic differentiation in Sitobion avenae populations from both sides of the Qinling Mountains

Quantitative trait differences are often assumed to be correlated with molecular variation, but the relationship is not certain, and empirical evidence is still scarce. To address this issue, we sampled six populations of the cereal aphid Sitobion avenae from areas north and south of the Qinling Mountains, and characterized their molecular variation at seven microsatellite loci and quantitative variation at nine life-history traits. Our results demonstrated that southern populations had slightly longer developmental times of nymphs but much higher lifetime fecundity, compared to northern populations. Of the nine tested quantitative characters, eight differed significantly among populations within regions, as well as between northern and southern regions. Genetic differentiation in neutral markers was likely to have been caused by founder events and drift. Increased subdivision for quantitative characters was found in northern populations, but reduced in southern populations. This phenomenon was not found for molecular characters, suggesting the decoupling between molecular and quantitative variation. The pattern of relationships between FST and QST indicated divergent selection and suggested that local adaptation play a role in the differentiation of life-history traits in tested S. avenae populations, particularly in those traits closely related to reproduction. The main role of natural selection over genetic drift was also supported by strong structural differences in G-matrices among S. avenae populations. However, cluster analyses did not result in two groups corresponding to northern and southern regions. Genetic differentiation between northern and southern populations in neutral markers was low, indicating considerable gene flow between them. The relationship between molecular and quantitative variation, as well as its implications for differentiation and evolution of S. avenae populations, was discussed.

Seascape analysis reveals regional gene flow patterns among populations of a marine planktonic diatom

We investigated the gene flow of the common marine diatom, Skeletonema marinoi, in Scandinavian waters and tested the null hypothesis of panmixia. Sediment samples were collected from the Danish Straits, Kattegat and Skagerrak. Individual strains were established from germinated resting stages. A total of 350 individuals were genotyped by eight microsatellite markers. Conventional F-statistics showed significant differentiation between the samples. We therefore investigated whether the genetic structure could be explained using genetic models based on isolation by distance (IBD) or by oceanographic connectivity. Patterns of oceanographic circulation are seasonally dependent and therefore we estimated how well local oceanographic connectivity explains gene flow month by month. We found no significant relationship between genetic differentiation and geographical distance. Instead, the genetic structure of this dominant marine primary producer is best explained by local oceanographic connectivity promoting gene flow in a primarily south to north direction throughout the year. Oceanographic data were consistent with the significant FST values between several pairs of samples. Because even a small amount of genetic exchange prevents the accumulation of genetic differences in F-statistics, we hypothesize that local retention at each sample site, possibly as resting stages, is an important component in explaining the observed genetic structure.

Morphological similarity and ecological overlap in two rotifer species

Co-occurrence of cryptic species raises theoretically relevant questions regarding their coexistence and ecological similarity. Given their great morphological similitude and close phylogenetic relationship (i.e., niche retention), these species will have similar ecological requirements and are expected to have strong competitive interactions. This raises the problem of finding the mechanisms that may explain the coexistence of cryptic species and challenges the conventional view of coexistence based on niche differentiation. The cryptic species complex of the rotifer Brachionus plicatilis is an excellent model to study these questions and to test hypotheses regarding ecological differentiation. Rotifer species within this complex are filtering zooplankters commonly found inhabiting the same ponds across the Iberian Peninsula and exhibit an extremely similar morphology-some of them being even virtually identical. Here, we explore whether subtle differences in body size and morphology translate into ecological differentiation by comparing two extremely morphologically similar species belonging to this complex: B. plicatilis and B. manjavacas. We focus on three key ecological features related to body size: (1) functional response, expressed by clearance rates; (2) tolerance to starvation, measured by growth and reproduction; and (3) vulnerability to copepod predation, measured by the number of preyed upon neonates. No major differences between B. plicatilis and B. manjavacas were found in the response to these features. Our results demonstrate the existence of a substantial niche overlap, suggesting that the subtle size differences between these two cryptic species are not sufficient to explain their coexistence. This lack of evidence for ecological differentiation in the studied biotic niche features is in agreement with the phylogenetic limiting similarity hypothesis but requires a mechanistic explanation of the coexistence of these species not based on differentiation related to biotic niche axes.

Polar lakes may act as ecological islands to aquatic protists

A fundamental question in ecology is whether microorganisms follow the same patterns as multicellular organisms when it comes to population structure and levels of genetic diversity. Enormous population sizes, predominately asexual reproduction and presumably high dispersal because of small body size could have profound implications on their genetic diversity and population structure. Here, we have analysed the population genetic structure in a lake-dwelling microbial eukaryote (dinoflagellate) and tested the hypothesis that there is population genetic differentiation among nearby lake subpopulations. This dinoflagellate occurs in the marine-derived saline lakes of the Vestfold Hills, Antarctica, which are ice-covered most of the year. Clonal strains were isolated from four different lakes and were genotyped using amplified fragment length polymorphism (AFLP). Our results show high genetic differentiation among lake populations despite their close geographic proximity (<9 km). Moreover, genotype diversity was high within populations. Gene flow in this system is clearly limited, either because of physical or biological barriers. Our results discard the null hypothesis that there is free gene flow among protist lake populations. Instead, limnetic protist populations may differentiate genetically, and lakes act as ecological islands even on the microbial scale.

Comparisons between Q(ST) and F(ST) --how wrong have we been?

The comparison between quantitative genetic divergence (Q(ST) ) and neutral genetic divergence (F(ST) ) among populations has become the standard test for historical signatures of selection on quantitative traits. However, when the mutation rate of neutral markers is relatively high in comparison with gene flow, estimates of F(ST) will decrease, resulting in upwardly biased comparisons of Q(ST) vs. F(ST) . Reviewing empirical studies, the difference between Q(ST) and F(ST) is positively related to marker heterozygosity. After refuting alternative explanations for this pattern, we conclude that marker mutation rate indeed has had a biasing effect on published Q(ST) -F(ST) comparisons. Hence, it is no longer clear that populations have commonly diverged in response to divergent selection. We present and discuss potential solutions to this bias. Comparing Q(ST) with recent indices of neutral divergence that statistically correct for marker heterozygosity (Hedrick's G'st and Jost's D) is not advised, because these indices are not theoretically equivalent to Q(ST) . One valid solution is to estimate F(ST) from neutral markers with mutation rates comparable to those of the loci underlying quantitative traits (e.g. SNPs). Q(ST) can also be compared to Φ(ST) (Phi(ST) ) of amova, as long as the genetic distance among allelic variants used to estimate Φ(ST) reflects evolutionary history: in that case, neutral divergence is independent of mutation rate. In contrast to their common usage in comparisons of Q(ST) and F(ST) , microsatellites typically have high mutation rates and do not evolve according to a simple evolutionary model, so are best avoided in Q(ST) -F(ST) comparisons.

Long-term coexistence of rotifer cryptic species

Despite their high morphological similarity, cryptic species often coexist in aquatic habitats presenting a challenge in the framework of niche differentiation theory and coexistence mechanisms. Here we use a rotifer species complex inhabiting highly unpredictable and fluctuating salt lakes to gain insights into the mechanisms involved in stable coexistence in cryptic species. We combined molecular barcoding surveys of planktonic populations and paleogenetic analysis of diapausing eggs to reconstruct the current and historical coexistence dynamics of two highly morphologically similar rotifer species, B. plicatilis and B. manjavacas. In addition, we carried out laboratory experiments using clones isolated from eight lakes where both species coexist to explore their clonal growth responses to salinity, a challenging, highly variable and unpredictable condition in Mediterranean salt lakes. We show that both species have co-occurred in a stable way in one lake, with population fluctuations in which no species was permanently excluded. The seasonal occurrence patterns of the plankton in two lakes agree with laboratory experiments showing that both species differ in their optimal salinity. These results suggest that stable species coexistence is mediated by differential responses to salinity and its fluctuating regime. We discuss the role of fluctuating salinity and a persistent diapausing egg banks as a mechanism for species coexistence in accordance with the 'storage effect'.

Widespread secondary contact and new glacial refugia in the halophilic rotifer Brachionus plicatilis in the Iberian Peninsula

Small aquatic organisms harbour deep phylogeographic patterns and highly structured populations even at local scales. These patterns indicate restricted gene flow, despite these organisms' high dispersal abilities, and have been explained by a combination of (1) strong founder effects due to rapidly growing populations and very large population sizes, and (2) the development of diapausing egg banks and local adaptation, resulting in low effective gene flow, what is known as the Monopolization hypothesis. In this study, we build up on our understanding of the mitochondrial phylogeography of the halophilic rotifer Brachionus plicatilis in the Iberian Peninsula by both increasing the number of sampled ponds in areas where secondary contact is likely and doubling sample sizes. We analyzed partial mitochondrial sequences of 252 individuals. We found two deep mitochondrial DNA lineages differing in both their genetic diversity and the complexity of their phylogeographic structure. Our analyses suggest that several events of secondary contact between clades occurred after their expansion from glacial refugia. We found a pattern of isolation-by-distance, which we interpret as being the result of historical colonization events. We propose the existence of at least one glacial refugium in the SE of the Iberian Peninsula. Our findings challenge predictions of the Monopolization hypothesis, since coexistence (i.e., secondary contact) of divergent lineages in some ponds in the Iberian Peninsula is common. Our results indicate that phylogeographic structures in small organisms can be very complex and that gene flow between diverse lineages after population establishment can indeed occur.

Genome size evolution at the speciation level: the cryptic species complex Brachionus plicatilis (Rotifera)

BACKGROUND

Studies on genome size variation in animals are rarely done at lower taxonomic levels, e.g., slightly above/below the species level. Yet, such variation might provide important clues on the tempo and mode of genome size evolution. In this study we used the flow-cytometry method to study the evolution of genome size in the rotifer Brachionus plicatilis, a cryptic species complex consisting of at least 14 closely related species.

RESULTS

We found an unexpectedly high variation in this species complex, with genome sizes ranging approximately seven-fold (haploid '1C' genome sizes: 0.056-0.416 pg). Most of this variation (67%) could be ascribed to the major clades of the species complex, i.e. clades that are well separated according to most species definitions. However, we also found substantial variation (32%) at lower taxonomic levels--within and among genealogical species--and, interestingly, among species pairs that are not completely reproductively isolated. In one genealogical species, called B. 'Austria', we found greatly enlarged genome sizes that could roughly be approximated as multiples of the genomes of its closest relatives, which suggests that whole-genome duplications have occurred early during separation of this lineage. Overall, genome size was significantly correlated to egg size and body size, even though the latter became non-significant after controlling for phylogenetic non-independence.

CONCLUSIONS

Our study suggests that substantial genome size variation can build up early during speciation, potentially even among isolated populations. An alternative, but not mutually exclusive interpretation might be that reproductive isolation tends to build up unusually slow in this species complex.

Hundred years of genetic structure in a sediment revived diatom population

This paper presents research on the genetic structure and diversity of populations of a common marine protist and their changes over time. The bloom-forming diatom Skeletonema marinoi was used as a model organism. Strains were revived from anoxic discrete layers of a (210)Pb-dated sediment core accumulated over more than 100 y, corresponding to >40,000 diatom mitotic generations. The sediment core was sampled from the highly eutrophic Mariager Fjord in Denmark. The genetic structure of S. marinoi was examined using microsatellite markers, enabling exploration of changes through time and of the effect of environmental fluctuations. The results showed a stable population structure among and within the examined sediment layers, and a similar genetic structure has been maintained over thousands of generations. However, established populations from inside the fjord were highly differentiated from open-sea populations. Despite constant water exchange and influx of potential colonizers into the fjord, the populations do not mix. One fjord population, accumulated in 1980, was significantly differentiated from the other groups of strains isolated from the fjord. This differentiation could have resulted from the status of Mariager Fjord, which was considered hypereutrophic, around 1980. There was no significant genetic difference between pre- and posteutrophication groups of strains. Our data show that dispersal potential and generation time do not have a large impact on the genetic structuring of the populations investigated here. Instead, the environmental conditions, such as the extreme eutrophication of the Mariager Fjord, are deemed more important.