The role of selection in driving landscape genomic structure of the waterflea Daphnia magna

Copy number and sequence variation in rDNA of Daphnia pulex from natural populations: insights from whole-genome sequencing

Ribosomal DNA (rDNA) has a vital role in ribosome biogenesis as it contains the genes that encode ribosomal RNA (rRNA) separated by intergenic spacers (IGSs). The rRNA genes occur in hundreds to tens of thousands of copies per haploid genome in eukaryotes and are generally highly conserved with low variation within species. Due to the repetitive nature and large size of rDNA arrays, detecting intraindividual variation can be difficult. In this study, we use whole-genome sequences of 169 Daphnia pulex individuals from 10 natural populations to measure the copy number and sequence variation in rDNA. This revealed that variation in rDNA copy number between individuals spans an order of magnitude. We further observed a substantial level of sequence variation within individual genomes. As expected, single-nucleotide polymorphisms occurred in regions of lower functional constraint such as the IGS and expansion segments of the rRNA genes. The presence of strong linkage disequilibrium among variants facilitated identification of haplotypes within each population. Although there was evidence of recombination among haplotypes from different populations, it is insufficient to eliminate linkage disequilibrium within populations. Estimating copy number and haplotype diversity within individuals revealed that the level of intraindividual sequence variation is not strongly correlated with copy number. The observed patterns of variation highlight a complex evolutionary history of rDNA in D. pulex. Future research should explore the functional implications of rDNA copy number and sequence variation on organismal phenotypes.

Macroecological predictors of evolutionary and plastic potential do not apply at microgeographic scales for a freshwater cladoceran under climate change

Rapid evolutionary adaptation could reduce the negative impacts of climate change if sufficient heritability of key traits exists under future climate conditions. Plastic responses to climate change could also reduce negative impacts. Understanding which populations are likely to respond via evolution or plasticity could therefore improve estimates of extinction risk. A large body of research suggests that the evolutionary and plastic potential of a population can be predicted by the degree of spatial and temporal climatic variation it experiences. However, we know little about the scale at which these relationships apply. Here, we test if spatial and temporal variation in temperature affects genetic variation and plasticity of fitness and a key thermal tolerance trait (critical thermal maximum; CTmax) at microgeographic scales using a metapopulation of Daphnia magna in freshwater rock pools. Specifically, we ask if (a) there is a microgeographic adaptation of CTmax and fitness to differences in temperature among the pools, (b) pools with greater temporal temperature variation have more genetic variation or plasticity in CTmax or fitness, and (c) increases in temperature affect the heritability of CTmax and fitness. Although we observed genetic variation and plasticity in CTmax and fitness, and differences in fitness among pools, we did not find support for the predicted relationships between temperature variation and genetic variation or plasticity. Furthermore, the genetic variation and plasticity we observed in CTmax are unlikely sufficient to reduce the impacts of climate change. CTmax plasticity was minimal and heritability was 72% lower when D. magna developed at the higher temperatures predicted under climate change. In contrast, the heritability of fitness increased by 53% under warmer temperatures, suggesting an increase in overall evolutionary potential unrelated to CTmax under climate change. More research is needed to understand the evolutionary and plastic potential under climate change and how that potential will be altered in future climates.

Population genetic differentiation of Daphnia sinensis in a lasting high-phosphorus Chinese lake, Lake Chaohu

Ecological shifts (e.g., eutrophication) can affect the genetic differentiation of zooplankton populations in lakes. However, the role of environmental change in a lasting high-phosphorus lake driving the genetic differentiation of zooplankton population structure over time is poorly understood. In this paper, the changes of the genetic diversity and differentiation of Daphnia sinensis population were studied by using the mitochondrial COI gene and microsatellite markers on modern groups (from January to June 2016) and historic groups (obtained from resting eggs in the sediments) in Lake Chaohu. Based on the microsatellite markers, six modern groups were clustered into two clusters (the WG cluster and SG cluster) during the seasonal dynamics, whereas the genetic differentiation of the five historic groups showed a wave-like pattern and had evolved into four clusters. Moreover, the haplotype network showed that six modern groups had one origin center whereas five historic groups had two origin centers based on the mitochondrial COI gene marker. Fu’s Fs neutral test and Tajima’s test indicated that the five historic groups deviated from neutral evolution and showed a bottleneck effect in the history process. Water temperature and total dissolved phosphorus were obviously associated with the seasonal genetic differentiation of D. sinensis, whereas nitrogen content of the sediments was significantly related to the long-term microevolution of D. sinensis in the high-phosphorus environment. Therefore, the changing pattern of D. sinensis population genetic structure was one of the environmental selections probably combined with co-evolutionary, where rapid-increasing nitrogen level had a large impact on D. sinensis population genetic structure in lasting high phosphorus environment in Lake Chaohu.

Spatial and temporal variation of genetic diversity and genetic differentiation in Daphnia galeata populations in four large reservoirs in southern China

Daphnia galeata is a common and dominant species in warmer waters, and has a strong top-down effect on both phytoplankton and bacteria. The knowledge of its temporal and spatial patterns of genetic diversity is fundamental in understanding its population dynamics and potential ecological function in ecosystems. Its population genetics have been investigated at regional scales but few within regions or at smaller spatial scales. Here, we examined the fine-scale spatial genetic variation of D. galeata within four large, deep reservoirs in wet and dry seasons and the six-year variation of genetic diversity in one of the reservoirs by using cytochrome c oxidase subunit I and microsatellites (simple sequence repeat). Our study shows that fine-scale spatial genetic variation commonly occurred within the reservoirs, indicating strong environmental selection at least in the two of reservoirs with strong longitudinal gradients. Since the environmental gradients established in the dry season was largely reduced in the wet season, the fine-scale spatial genetic variation was much higher in the dry season. The dynamics of local genetic diversity did not follow the theoretical pattern of rapid erosion but peaked in mid or mid-late growth season. The local genetic diversity of D. galeata appears to be shaped and maintained not only by recruitment from resting egg banks but also by gene flow within reservoirs. The temporal and fine-scale genetic variation within a water body suggests that it is necessary to pay attention to sampling periods and locations of a given water body in regional studies.

Evolutionary Genomics of a Subdivided Species

The ways in which genetic variation is distributed within and among populations is a key determinant of the evolutionary features of a species. However, most comprehensive studies of these features have been restricted to studies of subdivision in settings known to have been driven by local adaptation, leaving our understanding of the natural dispersion of allelic variation less than ideal. Here, we present a geographic population-genomic analysis of 10 populations of the freshwater microcrustacean Daphnia pulex, an emerging model system in evolutionary genomics. These populations exhibit a pattern of moderate isolation-by-distance, with an average migration rate of 0.6 individuals per generation, and average effective population sizes of ∼650,000 individuals. Most populations contain numerous private alleles, and genomic scans highlight the presence of islands of excessively high population subdivision for more common alleles. A large fraction of such islands of population divergence likely reflect historical neutral changes, including rare stochastic migration and hybridization events. The data do point to local adaptive divergence, although the precise nature of the relevant variation is diffuse and cannot be associated with particular loci, despite the very large sample sizes involved in this study. In contrast, an analysis of between-species divergence highlights positive selection operating on a large set of genes with functions nearly nonoverlapping with those involved in local adaptation, in particular ribosome structure, mitochondrial bioenergetics, light reception and response, detoxification, and gene regulation. These results set the stage for using D. pulex as a model for understanding the relationship between molecular and cellular evolution in the context of natural environments.

Low genome-wide divergence between two lizard populations with high adaptive phenotypic differentiation

Usually, adaptive phenotypic differentiation is paralleled by genetic divergence between locally adapted populations. However, adaptation can also happen in a scenario of nonsignificant genetic divergence due to intense gene flow and/or recent differentiation. While this phenomenon is rarely published, findings on incipient ecologically driven divergence or isolation by adaptation are relatively common, which could confound our understanding about the frequency at which they actually occur in nature. Here, we explore genome-wide traces of divergence between two populations of the lacertid lizard Psammodromus algirus separated by a 600 m elevational gradient. These populations seem to be differentially adapted to their environments despite showing low levels of genetic differentiation (according to previously studies of mtDNA and microsatellite data). We performed a search for outliers (i.e., loci subject to selection) trying to identify specific loci with FST statistics significantly higher than those expected on the basis of overall, genome-wide estimates of genetic divergence. We find that local phenotypic adaptation (in terms of a wide diversity of characters) was not accompanied by genome-wide differentiation, even when we maximized the chances of unveiling such differentiation at particular loci with FST-based outlier detection tests. Instead, our analyses confirmed the lack of genome-wide differentiation on the basis of more than 70,000 SNPs, which is concordant with a scenario of local adaptation without isolation by environment. Our results add evidence to previous studies in which local adaptation does not lead to any kind of isolation (or early stages of ecological speciation), but maintains phenotypic divergence despite the lack of a differentiated genomic background.

Impact of Limited Dispersion Capacity and Natural Barriers on the Population Structure of the Grasshopper Ommexecha virens (Orthoptera: Ommexechidae)

The grasshopper Ommexecha virens Serville has low dispersion capacity, and it is regarded as a specialist, only being found in sandy, dry environments with high incidence of sunlight. Considering these aspects, we evaluated the diversity and genetic structure of O. virens natural populations using ISSR (Inter Simple Sequence Repeat) markers. The data pointed to low expected heterozygosity for some populations (H_E = 0.06-0.09), probably a consequence of positive inbreeding, which is typical of species showing low or null dispersion indices. Moreover, significant genetic differentiation was observed (F_ST = 0.50 and G_ST = 0.51), as well as low number of migrants (N_m = 0.47), indicating that the populations are genetically differentiated. This is likely related to the limitation in dispersing and fragmentation of suitable environment localities colonized by O. virens. The populations of O. virens were structured in three genetic groups associated to different landscapes, revealing the presence of a secondary contact zone, possibly arisen from isolation followed by genetic divergence among populations and subsequent gene flow of divergent individuals of O. virens. At last, we found positive isolation by distance (IBD; r: 0.427; P: 0.025) which is an important factor, since it may be adding to the emergence of reproductive barriers among individuals of O. virens that have been experiencing isolation.

Evidence of local adaptation despite strong drift in a Neotropical patchily distributed bromeliad

Both genetic drift and divergent selection are predicted to be drivers of population differentiation across patchy habitats, but the extent to which these forces act on natural populations to shape traits is strongly affected by species' ecological features. In this study, we infer the genomic structure of Pitcairnia lanuginosa, a widespread herbaceous perennial plant with a patchy distribution. We sampled populations in the Brazilian Cerrado and the Central Andean Yungas and discovered and genotyped SNP markers using double-digest restriction-site associated DNA sequencing. In addition, we analyzed ecophysiological traits obtained from a common garden experiment and compared patterns of phenotypic and genetic divergence (P_ST-F_ST comparisons) in a subset of populations from the Cerrado. Our results from molecular analyses pointed to extremely low genetic diversity and a remarkable population differentiation, supporting a major role of genetic drift. Approximately 0.3% of genotyped SNPs were flagged as differentiation outliers by at least two distinct methods, and Bayesian generalized linear mixed models revealed a signature of isolation by environment in addition to isolation by distance for high-differentiation outlier SNPs among the Cerrado populations. P_ST-F_ST comparisons suggested divergent selection on two ecophysiological traits linked to drought tolerance. We showed that these traits vary among populations, although without any particular macro-spatial pattern, suggesting local adaptation to differences in micro-habitats. Our study shows that selection might be a relevant force, particularly for traits involved in drought stress, even for populations experiencing strong drift, which improves our knowledge on eco-evolutionary processes acting on non-continuously distributed species.

Regarding the F-word: The effects of data filtering on inferred genotype-environment associations

Genotype-environment association (GEA) methods have become part of the standard landscape genomics toolkit, yet, we know little about how to best filter genotype-by-sequencing data to provide robust inferences for environmental adaptation. In many cases, default filtering thresholds for minor allele frequency and missing data are applied regardless of sample size, having unknown impacts on the results, negatively affecting management strategies. Here, we investigate the effects of filtering on GEA results and the potential implications for assessment of adaptation to environment. We use empirical and simulated data sets derived from two widespread tree species to assess the effects of filtering on GEA outputs. Critically, we find that the level of filtering of missing data and minor allele frequency affect the identification of true positives. Even slight adjustments to these thresholds can change the rate of true positive detection. Using conservative thresholds for missing data and minor allele frequency substantially reduces the size of the data set, lessening the power to detect adaptive variants (i.e., simulated true positives) with strong and weak strengths of selection. Regardless, strength of selection was a good predictor for GEA detection, but even some SNPs under strong selection went undetected. False positive rates varied depending on the species and GEA method, and filtering significantly impacted the predictions of adaptive capacity in downstream analyses. We make several recommendations regarding filtering for GEA methods. Ultimately, there is no filtering panacea, but some choices are better than others, depending on the study system, availability of genomic resources, and desired objectives.

Effects of Colonization, Geography and Environment on Genetic Divergence in the Intermediate Leaf-Nosed Bat, Hipposideros larvatus

Determining the evolutionary history and population drivers, such as past large-scale climatic oscillations, stochastic processes and ecological adaptations, represents one of the aims of evolutionary biology. Hipposideros larvatus is a common bat species in Southern China, including Hainan Island. We examined genetic variation in H. larvatus using mitochondrial DNA and nuclear microsatellites. We found a population structure on both markers with a geographic pattern that corresponds well with the structure on mainland China and Hainan Island. To understand the contributions of geography, the environment and colonization history to the observed population structure, we tested isolation by distance (IBD), isolation by adaptation (IBA) and isolation by colonization (IBC) using serial Mantel tests and RDA analysis. The results showed significant impacts of IBD, IBA and IBC on neutral genetic variation, suggesting that genetic variation in H. larvatus is greatly affected by neutral processes, environmental adaptation and colonization history. This study enriches our understanding of the complex evolutionary forces that shape the distribution of genetic variation in bats.

Population connectivity, dispersal, and swimming behavior in Daphnia

The water flea Daphnia has the capacity to respond rapidly to environmental stressors, to disperse over large geographical scales, and to preserve its genetic material by forming egg banks in the sediment. Spatial and temporal distributions of D. magna have been extensively studied over the last decades using behavioral or genetic tools, although the correlation between the two has rarely been the focus. In the present study, we therefore investigated the population genetic structure and behavioral response to a lethal threat, ultraviolet radiation (UVR), among individuals from two different water bodies. Our results show two genetic populations with moderate gene flow, highly correlated with geographical location and with inheritable traits through generations. However, despite the strong genetic differences between populations, we show homogeneous refuge demand between populations when exposed to the lethal threat solar UVR.

Freshwater zooplankton metapopulations and metacommunities respond differently to environmental and spatial variation

Theory predicts that population genetic structure and metacommunity structure are linked by the common processes of drift and migration, but how population genetic structure and metacommunity structure are related in nature is still unknown. Deeper understanding of the processes influencing both genetic and community diversity is vital for better predicting how environmental change will impact biodiversity patterns. We examined how crustacean zooplankton and rotifer species' metapopulation genetic structure and metacommunities respond to environmental and spatial variation both within and across four regions of boreal Canada. Metapopulation and metacommunity variation partitioning results were compared within and across the four regions. Metapopulations and metacommunities responded differently to environmental variation and spatial structure both within and across regions, as metapopulations were influenced by different environmental variables compared to metacommunities. At larger spatial scales both metapopulations and metacommunities exhibited greater spatial and environmental structuring, again responding to a different subset of environmental variables. Our findings suggest that even though both genetic and species diversity are linked by the same processes, regional variation in environmental characteristics and spatial structure influence resulting biodiversity patterns differently. To date, no other empirical research has explored relationships between entire metapopulation and metacommunity assemblages at large regional spatial scales.

Elucidating the drivers of genetic differentiation in Malaysian torrent frogs (Anura: Ranidae: Amolops): a landscape genomics approach

The interplay between environmental attributes and evolutionary processes can provide valuable insights into how biodiversity is generated, partitioned and distributed. This study investigates the role of spatial, environmental and historical factors that could potentially drive diversification and shape genetic variation in Malaysian torrent frogs. Torrent frogs are ecologically conserved, and we hypothesize that this could impose tight constraints on dispersal routes, gene flow and consequently genetic structure. Moreover, levels of gene flow were shown to vary among populations from separate mountain ranges, indicating that genetic differentiation could be influenced by landscape features. Using genome-wide single nucleotide polymorphisms, in conjunction with landscape variables derived from Geographic Information Systems, we performed distance-based redundancy analyses and variance partitioning to disentangle the effects of isolation-by-distance (IBD), isolation-by-resistance (IBR) and isolation-by-colonization (IBC). Our results demonstrated that IBR contributed minimally to genetic variation. Intraspecific population structure can be largely attributed to IBD, whereas interspecific diversification was primarily driven by IBC. We also detected two distinct population bottlenecks, indicating that speciation events were likely driven by vicariance or founder events.

Population Structure of Colletotrichum tanaceti in Australian Pyrethrum Reveals High Evolutionary Potential

Colletotrichum tanaceti, the causal agent of anthracnose, is an emerging pathogen of commercially grown pyrethrum (Tanacetum cinerariifolium) in Australia. A microsatellite marker library was developed to understand the spatio-genetic structure over three sampled years and across two regions where pyrethrum is cultivated in Australia. Results indicated that C. tanaceti was highly diverse with a mixed reproductive mode; comprising both sexual and clonal reproduction. Sexual reproduction of C. tanaceti was more prevalent in Tasmania than in Victoria. Little differentiation was observed among field populations likely due to isolation by colonization but most of the genetic variation was occurring within populations. C. tanaceti was likely to have had a long-distance gene and genotype flow among distant populations within a state and between states. Anthropogenic transmission of propagules and wind dispersal of ascospores are the most probable mechanisms of long-distance dispersal of C. tanaceti. Evaluation of putative population histories suggested that C. tanaceti most likely originated in Tasmania and expanded from an unidentified host onto pyrethrum. Victoria was later invaded by the Tasmanian population. With the mixed mode of reproduction and possible long-distance gene flow, C. tanaceti is likely to have a high evolutionary potential and thereby has ability to adapt to management practices in the future.

Unexpected endemism in the Daphnia longispina complex (Crustacea: Cladocera) in Southern Siberia

The biological significance of regional cladoceran morphotypes in the montane regions of the central Palearctic remains poorly understood. In the Holarctic Daphnia longispina complex (Cladocera: Daphniidae), several variants, lineages and species have been proposed as endemic for Southern Siberia. Daphnia turbinata Sars, for example, named after its unusual head shape, is known only from Southern Siberia. Here we sequence DNA of Daphnia from three mitochondrial genes (12S rRNA, 16S rRNA, and NADH dehydrogenase subunit 2, ND2) from 57 localities in Russia and Mongolia (the majority being from Southern Siberia) and place them in evolutionary context with existing data. Our aim was to examine regional endemism of the Daphnia longispina complex in Southern Siberian; to improve the phylogenetic understanding with improved taxonomic and regional sampling, and to better understand the influence of Pleistocene glaciation on the biogeography of these lineages. At least three lineages showed genetic evidence for endemism in Southern Siberia. There was strong support for D. turbinata as a sister lineage to to D. longispina/D. dentifera. Another endemic, Siberian D. cf. longispina, is a sister group to the longispina group in general. Within D. longispina s. str. there was an endemic Siberian clade with a western range boundary near the Yenisei River Basin. Gene flow estimates among populations (based on F_ST values) were very low for clades of D. longispina on a regional (the original 12S dataset), and on a pan-Eurasian (the extended 12S dataset) scale. Negative values of Fu’s F_S and Tajima’s D tests prevailed for the species examined with significant values found for two D. longispina clades, D. dentifera, D. galeata and D. cristata. Our results support the notion that Southern Siberia is an important biogeographic region for cladocerans as it contained unexpected diversity of endemics (such as D. turbinata, D. cf. longispina and lineages of D. umbra and D. longsipina s.str.) and from being the geographic meeting place of expanding postglacial lineages from eastern and western refugia.

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.

Adaptive genetic differentiation in Pterocarya stenoptera (Juglandaceae) driven by multiple environmental variables were revealed by landscape genomics

BACKGROUND

The investigation of the genetic basis of local adaptation in non-model species is an interesting focus of evolutionary biologists and molecular ecologists. Identifying these adaptive genetic variabilities on the genome responsible can provide insight into the genetic mechanism of local adaptation.

RESULTS

We investigated the spatial distribution of genetic variation in 22 natural populations of Pterocarya stenoptera across its distribution area in China to provide insights into the complex interplay between multiple environmental variables and adaptive genetic differentiation. The Bayesian analysis of population structure showed that the 22 populations of P. stenoptera were subdivided into two groups. Redundancy analysis demonstrated that this genetic differentiation was caused by the divergent selection of environmental difference. A total of 44 outlier loci were mutually identified by Arlequin and BayeScan, 43 of which were environment-associated loci (EAL). The results of latent factor mixed model analysis showed that solar radiation in June (Sr6), minimum temperature of the coldest month (Bio6), temperature seasonality (Bio4), and water vapor pressure in January (Wvp1) were associated with the highest numbers of EAL. Sr6 was associated with the ecological habitat of "prefered light", and Bio6 and Wvp1 were associated with the ecological habitat of "warm and humid environment".

CONCLUSIONS

Our results provided empirical evidence that environmental variables related to the ecological habitats of species play key roles in driving adaptive differentiation of species genome.

Genome-wide SNP analyses reveal high gene flow and signatures of local adaptation among the scalloped spiny lobster (Panulirus homarus) along the Omani coastline

BACKGROUND

The scalloped spiny lobster (Panulirus homarus) is a popular seafood commodity worldwide and an important export item from Oman. Annual catches in commercial fisheries are in serious decline, which has resulted in calls for the development of an integrated stock management approach. In Oman, the scalloped spiny lobster is currently treated as a single management unit (MU) or stock and there is an absence of information on the genetic population structure of the species that can inform management decisions, particularly at a fine-scale level. This work is the first to identify genome-wide single nucleotide polymorphisms (SNPs) for P. homarus using Diversity Arrays Technology sequencing (DArT-seq) and to elucidate any stock structure in the species.

RESULTS

After stringent filtering, 7988 high utility SNPs were discovered and used to assess the genetic diversity, connectivity and structure of P. homarus populations from Al Ashkharah, Masirah Island, Duqm, Ras Madrakah, Haitam, Ashuwaymiyah, Mirbat and Dhalkut landing sites. Pairwise FST estimates revealed low differentiation among populations (pairwise FST range = - 0.0008 - 0.0021). Analysis of genetic variation using putatively directional FST outliers (504 SNPs) revealed higher and significant pairwise differentiation (p < 0.01) for all locations, with Ashuwaymiyah being the most diverged population (Ashuwaymiyah pairwise FST range = 0.0288-0.0736). Analysis of population structure using Discriminant Analysis of Principal Components (DAPC) revealed a broad admixture among P. homarus, however, Ashuwaymiyah stock appeared to be potentially under local adaptive pressures. Fine scale analysis using Netview R provided further support for the general admixture of P. homarus.

CONCLUSIONS

Findings here suggested that stocks of P. homarus along the Omani coastline are admixed. Yet, fishery managers need to treat the lobster stock from Ashuwaymiyah with caution as it might be subject to local adaptive pressures. We emphasize further study with larger number of samples to confirm the genetic status of the Ashuwaymiyah stock. The approach utilised in this study has high transferability in conservation and management of other marine stocks with similar biological and ecological attributes.

Cracking the Code of Biodiversity Responses to Past Climate Change

How individual species and entire ecosystems will respond to future climate change are among the most pressing questions facing ecologists. Past biodiversity dynamics recorded in the paleoecological archives show a broad array of responses, yet significant knowledge gaps remain. In particular, the relative roles of evolutionary adaptation, phenotypic plasticity, and dispersal in promoting survival during times of climate change have yet to be clarified. Investigating the paleo-archives offers great opportunities to understand biodiversity responses to future climate change. In this review we discuss the mechanisms by which biodiversity responds to environmental change, and identify gaps of knowledge on the role of range shifts and tolerance. We also outline approaches at the intersection of paleoecology, genomics, experiments, and predictive models that will elucidate the processes by which species have survived past climatic changes and enhance predictions of future changes in biological diversity.

Subsampling reveals that unbalanced sampling affects STRUCTURE results in a multi-species dataset

Studying the genetic population structure of species can reveal important insights into several key evolutionary, historical, demographic, and anthropogenic processes. One of the most important statistical tools for inferring genetic clusters is the program STRUCTURE. Recently, several papers have pointed out that STRUCTURE may show a bias when the sampling design is unbalanced, resulting in spurious joining of underrepresented populations and spurious separation of overrepresented populations. Suggestions to overcome this bias include subsampling and changing the ancestry model, but the performance of these two methods has not yet been tested on actual data. Here, I use a data set of 12 high-alpine plant species to test whether unbalanced sampling affects the STRUCTURE inference of population differentiation between the European Alps and the Carpathians. For four of the 12 species, subsampling of the Alpine populations-to match the sample size between the Alps and the Carpathians-resulted in a drastically different clustering than the full data set. On the other hand, STRUCTURE results with the alternative ancestry model were indistinguishable from the results with the default model. Based on these results, the subsampling strategy seems a more viable approach to overcome the bias than the alternative ancestry model. However, subsampling is only possible when there is an a priori expectation of what constitute the main clusters. Though these results do not mean that the use of STRUCTURE should be discarded, it does indicate that users of the software should be cautious about the interpretation of the results when sampling is unbalanced.

Weak population structure and no genetic erosion in Pilosocereus aureispinus: A microendemic and threatened cactus species from eastern Brazil

Succulent cacti (Cactaceae) are among the most threatened taxonomic groups assessed to date. Here we evaluated the genetic diversity and population structure of a narrow endemic columnar cactus Pilosocereus aureispinus. This species is only found in a small area of c. 300 km² of rocky savanna from eastern Brazil and it is currently classified as vulnerable (VU) on the International Union for Conservation of Nature (IUCN) red list. Eight microsatellite loci were genotyped for 91 individuals from four localities of the known P. aureispinus range. In contrast with expectations for narrow endemic species, we found relatively high levels of genetic diversity (e.g., H_E = 0.390 to 0.525; H_O = 0.394 to 0.572) and very low population structure based on the variation of six loci. All the analyzed individuals were clustered in one unique genetic group in assignment tests. We also generated the sequences of two plastid markers (trnT-trnL and psbD-trnT) and found no variation on a subsample of 39 individuals. We used Landsat 8 images and Normalized Difference Vegetation Index to estimate a potential extent of occurrence of c. 750 km² for this species. Our results showed that P. aureispinus is not suffering from erosion of nuclear genetic variability due to its narrow distribution. However, we advocate that because of the extremely limited extent of occurrence, the ongoing anthropogenic disturbances in its habitat, and phylogenetic distinctiveness of P. aureispinus, this species should be classified as endangered (EN) on the IUCN Red List.

A planktonic diatom displays genetic structure over small spatial scales

Marine planktonic microalgae have potentially global dispersal, yet reduced gene flow has been confirmed repeatedly for several species. Over larger distances (>200 km) geographic isolation and restricted oceanographic connectivity have been recognized as instrumental in driving population divergence. Here we investigated whether similar patterns, that is, structured populations governed by geographic isolation and/or oceanographic connectivity, can be observed at smaller (6-152 km) geographic scales. To test this we established 425 clonal cultures of the planktonic diatom Skeletonema marinoi collected from 11 locations in the Archipelago Sea (northern Baltic Sea). The region is characterized by a complex topography, entailing several mixing regions of which four were included in the sampling area. Using eight microsatellite markers and conventional F-statistics, significant genetic differentiation was observed between several sites. Moreover, Bayesian cluster analysis revealed the co-occurrence of two genetic groups spread throughout the area. However, geographic isolation and oceanographic connectivity could not explain the genetic patterns observed. Our data reveal hierarchical genetic structuring whereby despite high dispersal potential, significantly diverged populations have developed over small spatial scales. Our results suggest that biological characteristics and historical events may be more important in generating barriers to gene flow than physical barriers at small spatial scales.

Phylogeography of Daphnia magna Straus (Crustacea: Cladocera) in Northern Eurasia: Evidence for a deep longitudinal split between mitochondrial lineages

Species with a large geographic distributions present a challenge for phylogeographic studies due to logistic difficulties of obtaining adequate sampling. For instance, in most species with a Holarctic distribution, the majority of studies has concentrated on the European or North American part of the distribution, with the Eastern Palearctic region being notably understudied. Here, we study the phylogeography of the freshwater cladoceran Daphnia magna Straus, 1820 (Crustacea: Cladocera), based on partial mitochondrial COI sequences and using specimens from populations spread longitudinally from westernmost Europe to easternmost Asia, with many samples from previously strongly understudied regions in Siberia and Eastern Asia. The results confirm the previously suspected deep split between Eastern and Western mitochondrial haplotype super-clades. We find a narrow contact zone between these two super-clades in the eastern part of Western Siberia, with proven co-occurrence in a single lake in the Novosibirsk region. However, at present there is no evidence suggesting that the two mitochondrial super-clades represent cryptic species. Rather, they may be explained by secondary contact after expansion from different refugia. Interestingly, Central Siberia has previously been found to be an important contact zone also in other cladoceran species, and may thus be a crucial area for understanding the Eurasian phylogeography of freshwater invertebrates. Together, our study provides an unprecedented complete, while still not global, picture of the phylogeography of this important model species.

Spiroides shrubs on Qinghai-Tibetan Plateau: Multilocus phylogeography and palaeodistributional reconstruction of Spiraea alpina and S. Mongolica (Rosaceae)

A common hypothesis for the rich biodiversity found in mountains is uplift-driven diversification. Using a multilocus approach, here we assessed the influence of Qinghai-Tibetan Plateau (QTP) uplift and fluctuating regional climate on genetic diversity of two sister spiroides shrubs, Spiraea alpina and S. mongolica. Combined with palaeodistributional reconstruction modelling, we investigated the current and past-predicted distribution of these species under different climatic episodes. The study demonstrated that continuous pulses of retreat and expansion during last glacial-interglacial episodes, combined with the uplifting of QTP shaped the current distribution of these species. All the populations showed high level of genetic diversity based on both cpDNA and SSR markers. The average gene diversity within populations based on cpDNA markers was 0.383 ± 0.052 for S. alpina and 0.477 ± 0.048 for S. mongolica. The observed and expected heterozygosities based on SSR for both Spiraea alpina and S. mongolicawere H_E(0.72-0.90)/H_O(0.35-0.78) and H_E(0.77-0.92)/H_O(0.47-0.77) respectively. Palaeodistributional reconstruction indicated species' preferences at southeastern edge of the plateau during last glacial maximum, at higher altitude areas of QTP and range expansion to central plateau during the interglacial episodes. Assignment tests in STRUCTURE, discriminant analysis of principal coordinates and Immigrants analysis in GENECLASS based on nuclear SSR markers did not support the hypothesis of gene flow between both the species. However, maximum likelihood approach based on cpDNA showed sharing of haplotypes between both species.

Ten Years of Landscape Genomics: Challenges and Opportunities

Landscape genomics is a relatively new discipline that aims to reveal the relationship between adaptive genetic imprints in genomes and environmental heterogeneity among natural populations. Although the interest in landscape genomics has increased since this term was coined, studies on this topic remain scarce. Landscape genomics has become a powerful method to scan and determine the genes responsible for the complex adaptive evolution of species at population (mostly) and individual (more rarely) level. This review outlines the sampling strategies, molecular marker types and research categories in 37 articles published during the first 10 years of this field (i.e., 2007-2016). We also address major challenges and future directions for landscape genomics. This review aims to promote interest in conducting additional studies in landscape genomics.

Simultaneous evaluation of the effects of geographic, environmental and temporal isolation in ecotypic populations of Solidago virgaurea

Early stages of ecological speciation can create populations with an ecology and reproduction timing distinct from those of related populations. Landscape genetic models incorporating environmental heterogeneity and population-specific reproductive traits enable the processes of population genetic differentiation to be inferred. We investigated genome-wide genetic variation in ecotypic populations of Solidago virgaurea sensu lato, a herbaceous plant inhabiting a wide range of habitats (woodlands, serpentine barrens and alpine grasslands) and displaying remarkable variation in flowering time. Simultaneous evaluation of environmental factors revealed an overwhelming effect of soil type differences on neutral genetic differentiation, compared with elevational differences. This result probably reflects the abrupt environmental changes generated by geological boundaries, whereas mountain slopes exhibit clinal changes, facilitating gene exchange between neighbouring populations. Temporal isolation was positively associated with genetic differentiation, with some early-flowering serpentine populations having allele frequencies distinct from adjacent nonserpentine populations. Overall, this study highlights the importance of ecological processes and of evolution of flowering time to promote genetic differentiation of S. virgaurea populations in a complex landscape.

Evolutionary epidemiology of schistosomiasis: linking parasite genetics with disease phenotype in humans

Here we assess the role of parasite genetic variation in host disease phenotype in human schistosomiasis by implementing concepts and techniques from environmental association analysis in evolutionary epidemiology. Schistosomiasis is a tropical disease that affects more than 200 million people worldwide and is caused by parasitic flatworms belonging to the genus Schistosoma. While the role of host genetics has been extensively studied and demonstrated, nothing is yet known on the contribution of parasite genetic variation to host disease phenotype in human schistosomiasis. In this study microsatellite genotypes of 1561 Schistosoma mansoni larvae collected from 44 human hosts in Senegal were linked to host characteristics such as age, gender, infection intensity, liver and bladder morbidity by means of multivariate regression methods (on each parasite locus separately). This revealed a highly significant association between allelic variation at the parasite locus L46951 and host infection intensity and bladder morbidity. Locus L46951 is located in the 3' untranslated region of the cGMP-dependent protein kinase gene that is expressed in reproductive organs of adult schistosome worms and appears to be linked to egg production. This putative link between parasite genetic variation and schistosomiasis disease phenotype sets the stage for further functional research.

Local selection in the presence of high levels of gene flow: Evidence of heterogeneous insecticide selection pressure across Ugandan Culex quinquefasciatus populations

BACKGROUND

Culex quinquefasciatus collected in Uganda, where no vector control interventions directly targeting this species have been conducted, was used as a model to determine if it is possible to detect heterogeneities in selection pressure driven by insecticide application targeting other insect species.

METHODOLOGY/PRINCIPAL FINDINGS

Population genetic structure was assessed through microsatellite analysis, and the impact of insecticide pressure by genotyping two target-site mutations, Vgsc-1014F of the voltage-gated sodium channel target of pyrethroid and DDT insecticides, and Ace1-119S of the acetylcholinesterase gene, target of carbamate and organophosphate insecticides. No significant differences in genetic diversity were observed among populations by microsatellite markers with HE ranging from 0.597 to 0.612 and low, but significant, genetic differentiation among populations (FST = 0.019, P = 0.001). By contrast, the insecticide-resistance markers display heterogeneous allelic distributions with significant differences detected between Central Ugandan (urban) populations relative to Eastern and Southwestern (rural) populations. In the central region, a frequency of 62% for Vgsc-1014F, and 32% for the Ace1-119S resistant allele were observed. Conversely, in both Eastern and Southwestern regions the Vgsc-1014F alleles were close to fixation, whilst Ace1-119S allele frequency was 12% (although frequencies may be underestimated due to copy number variation at both loci).

CONCLUSIONS/SIGNIFICANCE

Taken together, the microsatellite and both insecticide resistance target-site markers provide evidence that in the face of intense gene flow among populations, disjunction in resistance frequencies arise due to intense local selection pressures despite an absence of insecticidal control interventions targeting Culex.

Considering reefscape configuration and composition in biophysical models advance seascape genetics

Previous seascape genetics studies have emphasized the role of ocean currents and geographic distances to explain the genetic structure of marine species, but the role of benthic habitat has been more rarely considered. Here, we compared the population genetic structure observed in West Pacific giant clam populations against model simulations that accounted habitat composition and configuration, geographical distance, and oceanic currents. Dispersal determined by geographical distance provided a modelled genetic structure in better agreement with the observations than dispersal by oceanic currents, possibly due to insufficient spatial resolution of available oceanographic and coastal circulation models. Considering both habitat composition and configuration significantly improved the match between simulated and observed genetic structures. This study emphasizes the importance of a reefscape genetics approach to population ecology, evolution and conservation in the sea.

Isolation mediates persistent founder effects on zooplankton colonisation in new temporary ponds

Understanding the colonisation process in zooplankton is crucial for successful restoration of aquatic ecosystems. Here, we analyzed the clonal and genetic structure of the cyclical parthenogenetic rotifer Brachionus plicatilis by following populations established in new temporary ponds during the first three hydroperiods. Rotifer populations established rapidly after first flooding, although colonisation was ongoing throughout the study. Multilocus genotypes from 7 microsatellite loci suggested that most populations (10 of 14) were founded by few clones. The exception was one of the four populations that persisted throughout the studied hydroperiods, where high genetic diversity in the first hydroperiod suggested colonisation from a historical egg bank, and no increase in allelic diversity was detected with time. In contrast, in another of these four populations, we observed a progressive increase of allelic diversity. This population became less differentiated from the other populations suggesting effective gene flow soon after its foundation. Allelic diversity and richness remained low in the remaining two, more isolated, populations, suggesting little gene flow. Our results highlight the complexity of colonisation dynamics, with evidence for persistent founder effects in some ponds, but not in others, and with early immigration both from external source populations, and from residual, historical diapausing egg banks.

Thermal tolerance in the keystone species Daphnia magna-a candidate gene and an outlier analysis approach

Changes in temperature have occurred throughout Earth's history. However, current warming trends exacerbated by human activities impose severe and rapid loss of biodiversity. Although understanding the mechanisms orchestrating organismal response to climate change is important, remarkably few studies document their role in nature. This is because only few systems enable the combined analysis of genetic and plastic responses to environmental change over long time spans. Here, we characterize genetic and plastic responses to temperature increase in the aquatic keystone grazer Daphnia magna combining a candidate gene and an outlier analysis approach. We capitalize on the short generation time of our species, facilitating experimental evolution, and the production of dormant eggs enabling the analysis of long-term response to environmental change through a resurrection ecology approach. We quantify plasticity in the expression of 35 candidate genes in D. magna populations resurrected from a lake that experienced changes in average temperature over the past century and from experimental populations differing in thermal tolerance isolated from a selection experiment. By measuring expression in multiple genotypes from each of these populations in control and heat treatments, we assess plastic responses to extreme temperature events. By measuring evolutionary changes in gene expression between warm- and cold-adapted populations, we assess evolutionary response to temperature changes. Evolutionary response to temperature increase is also assessed via an outlier analysis using EST-linked microsatellite loci. This study provides the first insights into the role of plasticity and genetic adaptation in orchestrating adaptive responses to environmental change in D. magna.

Paleogenetic records of Daphnia pulicaria in two North American lakes reveal the impact of cultural eutrophication

Understanding the evolutionary consequences of the green revolution, particularly in wild populations, is an important frontier in contemporary biology. Because human impacts have occurred at varying magnitudes or time periods depending on the study ecosystem, evolutionary histories may vary considerably among populations. Paleogenetics in conjunction with paleolimnology enable us to associate microevolutionary dynamics with detailed information on environmental change. We used this approach to reconstruct changes in the temporal population genetic structure of the keystone zooplankton grazer, Daphnia pulicaria, using dormant eggs extracted from sediments in two Minnesota lakes (South Center, Hill). The extent of agriculture and human population density in the catchment of these lakes has differed markedly since European settlement in the late 19th century and is reflected in their environmental histories reconstructed here. The reconstructed environments of these two lakes differed strongly in terms of environmental stability and their associated patterns of Daphnia population structure. We detected long periods of stability in population structure and environmental conditions in South Center Lake that were followed by a dramatic temporal shift in population genetic structure after the onset of European settlement and industrialized agriculture in its watershed. In particular, we noted a 24.3-fold increase in phosphorus (P) flux between pre-European and modern sediment P accumulation rates (AR) in this lake. In contrast, no such shifts were detected in Hill Lake, where the watershed was not as impacted by European settlement and rates of change were less directional with a much smaller increase in sediment P AR (2.3-fold). We identify direct and indirect effects of eutrophication proxies on genetic structure in these lake populations and demonstrate the power of using this approach in understanding the consequences of anthropogenic environmental change on natural populations throughout historic time periods.

Temporal genetic stability in natural populations of the waterflea Daphnia magna in response to strong selection pressure

Studies monitoring changes in genetic diversity and composition through time allow a unique understanding of evolutionary dynamics and persistence of natural populations. However, such studies are often limited to species with short generation times that can be propagated in the laboratory or few exceptional cases in the wild. Species that produce dormant stages provide powerful models for the reconstruction of evolutionary dynamics in the natural environment. A remaining open question is to what extent dormant egg banks are an unbiased representation of populations and hence of the species' evolutionary potential, especially in the presence of strong environmental selection. We address this key question using the water flea Daphnia magna, which produces dormant stages that accumulate in biological archives over time. We assess temporal genetic stability in three biological archives, previously used in resurrection ecology studies showing adaptive evolutionary responses to rapid environmental change. We show that neutral genetic diversity does not decline with the age of the population and it is maintained in the presence of strong selection. In addition, by comparing temporal genetic stability in hatched and unhatched populations from the same biological archive, we show that dormant egg banks can be consulted to obtain a reliable measure of genetic diversity over time, at least in the multidecadal time frame studied here. The stability of neutral genetic diversity through time is likely mediated by the buffering effect of the resting egg bank.

A portrait of a sucker using landscape genetics: how colonization and life history undermine the idealized dendritic metapopulation

Dendritic metapopulations have been attributed unique properties by in silico studies, including an elevated genetic diversity relative to a panmictic population of equal total size. These predictions have not been rigorously tested in nature, nor has there been full consideration of the interacting effects among contemporary landscape features, colonization history and life history traits of the target species. We tested for the effects of dendritic structure as well as the relative importance of life history, environmental barriers and historical colonization on the neutral genetic structure of a longnose sucker (Catostomus catostomus) metapopulation in the Kogaluk watershed of northern Labrador, Canada. Samples were collected from eight lakes, genotyped with 17 microsatellites, and aged using opercula. Lakes varied in differentiation, historical and contemporary connectivity, and life history traits. Isolation by distance was detected only by removing two highly genetically differentiated lakes, suggesting a lack of migration-drift equilibrium and the lingering influence of historical factors on genetic structure. Bayesian analyses supported colonization via the Kogaluk's headwaters. The historical concentration of genetic diversity in headwaters inferred by this result was supported by high historical and contemporary effective sizes of the headwater lake, T-Bone. Alternatively, reduced allelic richness in headwaters confirmed the dendritic structure's influence on gene flow, but this did not translate to an elevated metapopulation effective size. A lack of equilibrium and upstream migration may have dampened the effects of dendritic structure. We suggest that interacting historical and contemporary factors prevent the achievement of the idealized traits of a dendritic metapopulation in nature.

Applying landscape genetics to the microbial world

Landscape genetics, which explicitly quantifies landscape effects on gene flow and adaptation, has largely focused on macroorganisms, with little attention given to microorganisms. This is despite overwhelming evidence that microorganisms exhibit spatial genetic structuring in relation to environmental variables. The increasing accessibility of genomic data has opened up the opportunity for landscape genetics to embrace the world of microorganisms, which may be thought of as 'the invisible regulators' of the macroecological world. Recent developments in bioinformatics and increased data accessibility have accelerated our ability to identify microbial taxa and characterize their genetic diversity. However, the influence of the landscape matrix and dynamic environmental factors on microorganism genetic dispersal and adaptation has been little explored. Also, because many microorganisms coinhabit or codisperse with macroorganisms, landscape genomic approaches may improve insights into how micro- and macroorganisms reciprocally interact to create spatial genetic structure. Conducting landscape genetic analyses on microorganisms requires that we accommodate shifts in spatial and temporal scales, presenting new conceptual and methodological challenges not yet explored in 'macro'-landscape genetics. We argue that there is much value to be gained for microbial ecologists from embracing landscape genetic approaches. We provide a case for integrating landscape genetic methods into microecological studies and discuss specific considerations associated with the novel challenges this brings. We anticipate that microorganism landscape genetic studies will provide new insights into both micro- and macroecological processes and expand our knowledge of species' distributions, adaptive mechanisms and species' interactions in changing environments.

Colonization history and clonal richness of asexual Daphnia in periglacial habitats of contrasting age in West Greenland

Due to climate change, Arctic ice sheets are retreating. This leads to the formation of numerous new periglacial ponds and lakes, which are being colonized by planktonic organisms such as the water flea Daphnia. This system provides unique opportunities to test genotype colonization dynamics and the genetic assemblage of populations. Here, we studied clonal richness of the Daphnia pulex species complex in novel periglacial habitats created by glacial retreat in the Jakobshavn Isbrae area of western Greenland. Along a 10 km transect, we surveyed 73 periglacial habitats out of which 61 were colonized by Daphnia pulex. Hence, for our analysis, we used 21 ponds and 40 lakes in two clusters of habitats differing in age (estimated <50 years vs. >150 years). We tested the expectation that genetic diversity would be low in recently formed (i.e. young), small habitats, but would increase with increasing age and size. We identified a total of 42 genetically distinct clones belonging to two obligately asexual species of the D. pulex species complex: D. middendorffiana and the much more abundant D. pulicaria. While regional clonal richness was high, most clones were rare: 16 clones were restricted to a single habitat and the five most widespread clones accounted for 68% of all individuals sampled. On average, 3·2 clones (range: 1-12) coexisted in a given pond or lake. There was no relationship between clonal richness and habitat size when we controlled for habitat age. Whereas clonal richness was statistically higher in the cluster of older habitats when compared with the cluster of younger ponds and lakes, most young habitats were colonized by multiple genotypes. Our data suggest that newly formed (periglacial) ponds and lakes are colonized within decades by multiple genotypes via multiple colonization events, even in the smallest of our study systems (4 m(2) ).

Molecular diversity of some species belonging to the genus Daphnia O. F. Müller, 1785 (Crustacea: Cladocera) in Turkey

Daphnia is a freshwater zooplankton species with controversial taxonomy due to its high morphological variation linked to environmental factors and inter-specific hybridization and polyploidy in some groups. The aim of the present study is to examine molecular diversity of some Daphnia species in Turkey and to establish DNA barcodes of Turkish Daphnia species. Sequence analysis was performed using 540 bp region of cytochrome oxidase subunit I gene of mitochondrial DNA. A total of 34 haplotypes have been identified for Turkey. Daphnia pulex complex was divided into two clades with 16.1% sequence divergence according to molecular taxonomy based on Kimura 2-parameter. The clade which was molecularly diverged from Daphnia pulex with 16.1% sequence divergence was found to show 99% similarity with Daphnia cf. pulicaria (sensu Alonso 1996) instead of Daphnia pulicaria Forbes, 1893. Furthermore, this study has contributed to Turkish zoogeography by demonstrating the distribution of Daphnia species in Turkey.

The role of parasite-driven selection in shaping landscape genomic structure in red grouse (Lagopus lagopus scotica)

Landscape genomics promises to provide novel insights into how neutral and adaptive processes shape genome-wide variation within and among populations. However, there has been little emphasis on examining whether individual-based phenotype-genotype relationships derived from approaches such as genome-wide association (GWAS) manifest themselves as a population-level signature of selection in a landscape context. The two may prove irreconcilable as individual-level patterns become diluted by high levels of gene flow and complex phenotypic or environmental heterogeneity. We illustrate this issue with a case study that examines the role of the highly prevalent gastrointestinal nematode Trichostrongylus tenuis in shaping genomic signatures of selection in red grouse (Lagopus lagopus scotica). Individual-level GWAS involving 384 SNPs has previously identified five SNPs that explain variation in T. tenuis burden. Here, we examine whether these same SNPs display population-level relationships between T. tenuis burden and genetic structure across a small-scale landscape of 21 sites with heterogeneous parasite pressure. Moreover, we identify adaptive SNPs showing signatures of directional selection using F(ST) outlier analysis and relate population- and individual-level patterns of multilocus neutral and adaptive genetic structure to T. tenuis burden. The five candidate SNPs for parasite-driven selection were neither associated with T. tenuis burden on a population level, nor under directional selection. Similarly, there was no evidence of parasite-driven selection in SNPs identified as candidates for directional selection. We discuss these results in the context of red grouse ecology and highlight the broader consequences for the utility of landscape genomics approaches for identifying signatures of selection.

Rapid Adaptation of a Daphnia magna Population to Metal Stress Is Associated with Heterozygote Excess

Although natural populations can harbor evolutionary potential to adapt genetically to chemical stress, it is often thought that natural selection leads to a general reduction of genetic diversity and involves costs. Here, a 10 week microevolution experiment was conducted with a genetically diverse and representative sample of one natural Daphnia magna population that was exposed to copper and zinc. Both Cu- and Zn-selected populations developed a significantly higher metal tolerance (i.e., genetic adaptation), indicated by higher reproduction probabilities of clonal lines in Cu and Zn exposures than observed for the original and control populations. The complete recovery of the population densities after 10 weeks of Zn selection (following an initial decrease of 74%) illustrates an example of evolutionary rescue. Microsatellite genotyping revealed a decrease in clonal diversity but no change in allelic richness, and showed an excess in heterozygosity in the Cu- and Zn-selected populations compared to the control and original populations. The excess heterozygosity in metal-selected populations that we observed has important consequences for risk assessment, as it contributes to the maintenance of a higher allelic diversity under multigenerational chemical exposure. This study is, to our knowledge, the first report of an increase in heterozygosity following multigenerational exposure to metal stress, despite a decline in clonal diversity. In a follow-up study with the Zn-selected populations, we observed no effect of Zn selection on the tolerance to heat and cyanobacteria. However, we observed higher tolerance to Cd in the Zn-selected than in the original and control populations if the 20% effective concentration of Cd was considered (cross-tolerance). Our results suggest only limited costs of adaptation but future research is needed to evaluate the adaptive potential of metal-selected populations to novel stressors and to determine to what extent increased heterozygosity is preserved after genetic recombination following periods of sexual reproduction.