Genetic diversity within and among Atlantic cod (Gadus morhua) farmed in marine cages: a proof-of-concept study for the identification of escapees

This study presents a molecular genetic characterization of Atlantic cod reared in commercial marine farms. Samples consisted of approximately 47 fish collected from nine cages located on four farms throughout Norway. In addition, 28 farmed escapees were recaptured in the sea (443 fish in total). Nine microsatellite loci and the Pan I gene were analysed, revealing a total of 181 alleles. Each sample contained 43–63% of total allelic variation. Comparing variation with published data for wild cod indicates that lower genetic variation exists within single cages than in wild populations. Significant linkage disequilibrium was observed amongst pairs of loci in all samples, suggesting a low number of contributing parental fish. Global F_ST was 0.049, and the highest pairwise F_ST value (pooled loci) was 0.085. For single loci, the Pan I gene was the most diagnostic, displaying a global F_ST of 0.203. Simulations amongst the samples collected on farms revealed an overall correct self-assignment percentage of 75%, demonstrating a high probability of identifying individuals to their farm of origin. Identification of the 28 escapees revealed a single cage as the most likely source of origin for half of the escapees, whilst the remaining fish were assigned to a mixture of samples, suggesting more than one source of escapees.     

Population structure of the malaria vector Anopheles sinensis (Diptera: Culicidae) in China: two gene pools inferred by microsatellites

Background

Anopheles sinensis is a competent malaria vector in China. An understanding of vector population structure is important to the vector-based malaria control programs. However, there is no adequate data of A. sinensis population genetics available yet.

Methodology/Principal Findings

This study used 5 microsatellite loci to estimate population genetic diversity, genetic differentiation and demographic history of A. sinensis from 14 representative localities in China. All 5 microsatellite loci were highly polymorphic across populations, with high allelic richness and heterozygosity. Hardy–Weinberg disequilibrium was found in 12 populations associated with heterozygote deficits, which was likely caused by the presence of null allele and the Wahlund effect. Bayesian clustering analysis revealed two gene pools, grouping samples into two population clusters; one includes six and the other includes eight populations. Out of 14 samples, six samples were mixed with individuals from both gene pools, indicating the coexistence of two genetic units in the areas sampled. The overall differentiation between two genetic pools was moderate (F _ST = 0.156). Pairwise differentiation between populations were lower within clusters (F _ST = 0.008–0.028 in cluster I and F _ST = 0.004–0.048 in cluster II) than between clusters (F _ST = 0.120–0.201). A reduced gene flow (Nm = 1–1.7) was detected between clusters. No evidence of isolation by distance was detected among populations neither within nor between the two clusters. There are differences in effective population size (Ne = 14.3-infinite) across sampled populations.

Conclusions/Significance

Two genetic pools with moderate genetic differentiation were identified in the A. sinensis populations in China. The population divergence was not correlated with geographic distance or barrier in the range. Variable effective population size and other demographic effects of historical population perturbations could be the factors affecting the population differentiation. The structured populations may limit the migration of genes under pressures/selections, such as insecticides and immune genes against malaria.     

Global phylogeography with mixed-marker analysis reveals male-mediated dispersal in the endangered scalloped hammerhead shark (Sphyrna lewini)

Background

The scalloped hammerhead shark, Sphyrna lewini, is a large endangered predator with a circumglobal distribution, observed in the open ocean but linked ontogenetically to coastal embayments for parturition and juvenile development. A previous survey of maternal (mtDNA) markers demonstrated strong genetic partitioning overall (global Φ_ST = 0.749) and significant population separations across oceans and between discontinuous continental coastlines.

Methodology/Principal Findings

We surveyed the same global range with increased sample coverage (N = 403) and 13 microsatellite loci to assess the male contribution to dispersal and population structure. Biparentally inherited microsatellites reveal low or absent genetic structure across ocean basins and global genetic differentiation (F _ST = 0.035) over an order of magnitude lower than the corresponding measures for maternal mtDNA lineages (Φ_ST = 0.749). Nuclear allelic richness and heterozygosity are high throughout the Indo-Pacific, while genetic structure is low. In contrast, allelic diversity is low while population structure is higher for populations at the ends of the range in the West Atlantic and East Pacific.

Conclusions/Significance

These data are consistent with the proposed Indo-Pacific center of origin for S. lewini, and indicate that females are philopatric or adhere to coastal habitats while males facilitate gene flow across oceanic expanses. This study includes the largest sampling effort and the most molecular loci ever used to survey the complete range of a large oceanic predator, and findings emphasize the importance of incorporating mixed-marker analysis into stock assessments of threatened and endangered shark species.     

Genetic Diversity in Introduced Golden Mussel Populations Corresponds to Vector Activity

We explored possible links between vector activity and genetic diversity in introduced populations of Limnoperna fortunei by characterizing the genetic structure in native and introduced ranges in Asia and South America. We surveyed 24 populations: ten in Asia and 14 in South America using the mitochondrial cytochrome c oxidase subunit I (COI) gene, as well as eight polymorphic microsatellite markers. We performed population genetics and phylogenetic analyses to investigate population genetic structure across native and introduced regions. Introduced populations in Asia exhibit higher genetic diversity (H _E = 0.667–0.746) than those in South America (H _E = 0.519–0.575), suggesting higher introduction effort for the former populations. We observed pronounced geographical structuring in introduced regions, as indicated by both mitochondrial and nuclear markers based on multiple genetic analyses including pairwise Ф_ST, F _ST, Bayesian clustering method, and three-dimensional factorial correspondence analyses. Pairwise F _ST values within both Asia (F _ST = 0.017–0.126, P = 0.000–0.009) and South America (F _ST = 0.004–0.107, P = 0.000–0.721) were lower than those between continents (F _ST = 0.180–0.319, P = 0.000). Fine-scale genetic structuring was also apparent among introduced populations in both Asia and South America, suggesting either multiple introductions of distinct propagules or strong post-introduction selection and demographic stochasticity. Higher genetic diversity in Asia as compared to South America is likely due to more frequent propagule transfers associated with higher shipping activities between source and donor regions within Asia. This study suggests that the intensity of human-mediated introduction vectors influences patterns of genetic diversity in non-indigenous species.     

Genetic screening of farmed Atlantic salmon escapees demonstrates that triploid fish display reduced migration to freshwater

Each year, hundreds of thousands of farmed Atlantic salmon escape from fish farms into the wild. Some of these escapees enter freshwater, and manage to interbreed with native populations. To hinder further genetic introgression in native populations, the use of sterile triploid salmon within commercial aquaculture is being examined. However, if triploid escapees migrate into freshwater, they may still have ecological impacts on local populations. In the present study, we used microsatellite DNA genotyping to determine the ploidy of 3794 farmed escapees captured in 17 Norwegian rivers in the period 2007–2014. Although a previous study has reported an average of 2 % triploids in Norwegian fish farms during this exact period, here, we only observed 7 (0.18 %) triploids among the escapees captured in freshwater. In addition, we identified three trisomic escapees. For the triploids where the within-river capture location was determined, they were only observed in the lower reaches and not on the spawning grounds. It is concluded that propensity for triploid Atlantic salmon to migrate into freshwater following escape from a fish farm is significantly lower than for normal diploid salmon escapees. Therefore, commercial production of triploids should not only be seen as an effective way of stopping genetic introgression, it will also significantly reduce the numbers of escapees entering rivers, which in turn limits ecological interactions and potential disease transmission.     

Genetic Variability of Wild Cherry (Prunus avium L.) Seed Stands in Slovenia as Revealed by Nuclear Microsatellite Loci

Microsatellite markers were used to describe the genetic variability of four seed stands of wild cherry (Prunus avium L.). One hundred and thirty one individuals were genotyped at ten nuclear microsatellite loci. Total genetic diversity was high (H _E = 0.704), while differences between stands were small but significant (F _ST = 0.053, G′ _ST = 0.234). There was a significant amount of clonal reproduction in one stand, with only 11 genotypes identified among 36 trees. One stand showed a significant excess (F _IS = −0.044) of heterozygosity, and one showed a deficit (F _IS = 0.044). Our results demonstrate the importance of taking into account the biological and genetic characteristics of species in forest management, especially when determining a new seed stand. The small genetic differences found between seed stands indicate that a large number of stands are not required. However, they should be carefully selected and should possess adequate genetic variability to ensure low relatedness between seed trees.     

Microsatellites reveal a high population structure in Triatoma infestans from Chuquisaca, Bolivia

Background

For Chagas disease, the most serious infectious disease in the Americas, effective disease control depends on elimination of vectors through spraying with insecticides. Molecular genetic research can help vector control programs by identifying and characterizing vector populations and then developing effective intervention strategies.

Methods and Findings

The population genetic structure of Triatoma infestans (Hemiptera: Reduviidae), the main vector of Chagas disease in Bolivia, was investigated using a hierarchical sampling strategy. A total of 230 adults and nymphs from 23 localities throughout the department of Chuquisaca in Southern Bolivia were analyzed at ten microsatellite loci. Population structure, estimated using analysis of molecular variance (AMOVA) to estimate F_ST (infinite alleles model) and R_ST (stepwise mutation model), was significant between western and eastern regions within Chuquisaca and between insects collected in domestic and peri-domestic habitats. Genetic differentiation at three different hierarchical geographic levels was significant, even in the case of adjacent households within a single locality (R _ST = 0.14, F _ST = 0.07). On the largest geographic scale, among five communities up to 100 km apart, R _ST = 0.12 and F _ST = 0.06. Cluster analysis combined with assignment tests identified five clusters within the five communities.

Conclusions

Some houses are colonized by insects from several genetic clusters after spraying, whereas other households are colonized predominately by insects from a single cluster. Significant population structure, measured by both R _ST and F _ST, supports the hypothesis of poor dispersal ability and/or reduced migration of T. infestans. The high degree of genetic structure at small geographic scales, inferences from cluster analysis and assignment tests, and demographic data suggest reinfesting vectors are coming from nearby and from recrudescence (hatching of eggs that were laid before insecticide spraying). Suggestions for using these results in vector control strategies are made.     

Geographic and temporal genetic structure in <Emphasis Type="Italic">Lepeophtheirus salmonis</Emphasis> from four salmon farms along the northwest and west coasts of Ireland: results from a microsatellite analysis

The marine ectoparasitic copepod of salmonids, Lepeophtheirus salmonis (Krøyer), is a major pest of farmed Atlantic salmon (Salmo salar L.) causing great economic impact. The spatial scales over which L. salmonis populations in different salmon farms are typically connected, and the temporal scales over which L. salmonis from the same farm typically undergo genetic change are largely unknown. These questions were posed in a small-scale geographic study of population structure in L. salmonis from four salmon farms, along the northwest and west coasts of Ireland, using two outgroups from Norway and Canada. The temporal stability of genetic composition was also studied in samples collected quarterly during one year from one salmon farm in Ireland. Genetic composition in L. salmonis was characterised using four nuclear microsatellites. Significant but low genetic differentiation was observed between all sites (F _ST = 0.08), with no evidence that differentiation was correlated with geographic distance. Temporal genetic differentiation was also evident (F _ST = 0.07). An analysis of all L. salmonis samples except the ones from Norway detected two separate clusters. Each cluster contained both geographical and temporal samples. These results are consistent with a population model in which L. salmonis in salmon farms along the northwest and west coasts of Ireland are not isolated, but are potentially subject to (i) localised ecological factors at the particular farm sites or (ii) selection post-settlement or a combination thereof.     

Genetic diversity and population history of a critically endangered primate, the northern muriqui (Brachyteles hypoxanthus)

Social, ecological, and historical processes affect the genetic structure of primate populations, and therefore have key implications for the conservation of endangered species. The northern muriqui (Brachyteles hypoxanthus) is a critically endangered New World monkey and a flagship species for the conservation of the Atlantic Forest hotspot. Yet, like other neotropical primates, little is known about its population history and the genetic structure of remnant populations. We analyzed the mitochondrial DNA control region of 152 northern muriquis, or 17.6% of the 864 northern muriquis from 8 of the 12 known extant populations and found no evidence of phylogeographic partitions or past population shrinkage/expansion. Bayesian and classic analyses show that this finding may be attributed to the joint contribution of female-biased dispersal, demographic stability, and a relatively large historic population size. Past population stability is consistent with a central Atlantic Forest Pleistocene refuge. In addition, the best scenario supported by an Approximate Bayesian Computation analysis, significant fixation indices (Φ_ST = 0.49, Φ_CT = 0.24), and population-specific haplotypes, coupled with the extirpation of intermediate populations, are indicative of a recent geographic structuring of genetic diversity during the Holocene. Genetic diversity is higher in populations living in larger areas (>2,000 hectares), but it is remarkably low in the species overall (θ = 0.018). Three populations occurring in protected reserves and one fragmented population inhabiting private lands harbor 22 out of 23 haplotypes, most of which are population-exclusive, and therefore represent patchy repositories of the species'' genetic diversity. We suggest that these populations be treated as discrete units for conservation management purposes.     

Invasion Pathway of the Ctenophore Mnemiopsis leidyi in the Mediterranean Sea

Gelatinous zooplankton outbreaks have increased globally owing to a number of human-mediated factors, including food web alterations and species introductions. The invasive ctenophore Mnemiopsis leidyi entered the Black Sea in the early 1980s. The invasion was followed by the Azov, Caspian, Baltic and North Seas, and, most recently, the Mediterranean Sea. Previous studies identified two distinct invasion pathways of M. leidyi from its native range in the western Atlantic Ocean to Eurasia. However, the source of newly established populations in the Mediterranean Sea remains unclear. Here we build upon our previous study and investigate sequence variation in both mitochondrial (Cytochrome c Oxidase subunit I) and nuclear (Internal Transcribed Spacer) markers in M. leidyi, encompassing five native and 11 introduced populations, including four from the Mediterranean Sea. Extant genetic diversity in Mediterranean populations (n = 8, N _a = 10) preclude the occurrence of a severe genetic bottleneck or founder effects in the initial colonizing population. Our mitochondrial and nuclear marker surveys revealed two possible pathways of introduction into Mediterranean Sea. In total, 17 haplotypes and 18 alleles were recovered from all surveyed populations. Haplotype and allelic diversity of Mediterranean populations were comparable to populations from which they were likely drawn. The distribution of genetic diversity and pattern of genetic differentiation suggest initial colonization of the Mediterranean from the Black-Azov Seas (pairwise F _ST = 0.001–0.028). However, some haplotypes and alleles from the Mediterranean Sea were not detected from the well-sampled Black Sea, although they were found in Gulf of Mexico populations that were also genetically similar to those in the Mediterranean Sea (pairwise F _ST = 0.010–0.032), raising the possibility of multiple invasion sources. Multiple introductions from a combination of Black Sea and native region sources could be facilitated by intense local and transcontinental shipping activity, respectively.     

Uniform selection as a primary force reducing population genetic differentiation of cavitation resistance across a species range

Background

Cavitation resistance to water stress-induced embolism determines plant survival during drought. This adaptive trait has been described as highly variable in a wide range of tree species, but little is known about the extent of genetic and phenotypic variability within species. This information is essential to our understanding of the evolutionary forces that have shaped this trait, and for evaluation of its inclusion in breeding programs.

Methodology

We assessed cavitation resistance (P ₅₀), growth and carbon isotope composition in six Pinus pinaster populations in a provenance and progeny trial. We estimated the heritability of cavitation resistance and compared the distribution of neutral markers (F _ST) and quantitative genetic differentiation (Q _ST), for retrospective identification of the evolutionary forces acting on these traits.

Results/Discussion

In contrast to growth and carbon isotope composition, no population differentiation was found for cavitation resistance. Heritability was higher than for the other traits, with a low additive genetic variance (h² _ns = 0.43±0.18, CV_A = 4.4%). Q _ST was significantly lower than F _ST, indicating uniform selection for P ₅₀, rather than genetic drift. Putative mechanisms underlying Q_STST are discussed.     

Genetic Diversity and Population Structure of Sitodiplosis mosellana in Northern China

The wheat midge, Sitodiplosis mosellana, is an important pest in Northern China. We tested the hypothesis that the population structure of this species arises during a range expansion over the past 30 years. This study used microsatellite and mitochondrial loci to conduct population genetic analysis of S. mosellana across its distribution range in China. We found strong genetic structure among the 16 studied populations, including two genetically distinct groups (the eastern and western groups), broadly consistent with the geography and habitat fragmentation. These results underline the importance of natural barriers in impeding dispersal and gene flow of S. mosellana populations. Low to moderate genetic diversity among the populations and moderate genetic differentiation (F _ST = 0.117) between the two groups were also found. The populations in the western group had lower genetic diversity, higher genetic differentiation and lower gene flow (F _ST = 0.116, Nm = 1.89) than those in the eastern group (F _ST = 0.049, Nm = 4.91). Genetic distance between populations was positively and significantly correlated with geographic distance (r = 0.56, P<0.001). The population history of this species provided no evidence for population expansion or bottlenecks in any of these populations. Our data suggest that the distribution of genetic diversity, genetic differentiation and population structure of S. mosellana have resulted from a historical event, reflecting its adaptation to diverse habitats and forming two different gene pools. These results may be the outcome of a combination of restricted gene flow due to geographical and environmental factors, population history, random processes of genetic drift and individual dispersal patterns. Given the current risk status of this species in China, this study can offer useful information for forecasting outbreaks and designing effective pest management programs.     

Genome-wide analysis of single nucleotide polymorphisms uncovers population structure in Northern Europe

Background

Genome-wide data provide a powerful tool for inferring patterns of genetic variation and structure of human populations.

Principal Findings

In this study, we analysed almost 250,000 SNPs from a total of 945 samples from Eastern and Western Finland, Sweden, Northern Germany and Great Britain complemented with HapMap data. Small but statistically significant differences were observed between the European populations (F_ST = 0.0040, p<10⁻⁴), also between Eastern and Western Finland (F_ST = 0.0032, p<10⁻³). The latter indicated the existence of a relatively strong autosomal substructure within the country, similar to that observed earlier with smaller numbers of markers. The Germans and British were less differentiated than the Swedes, Western Finns and especially the Eastern Finns who also showed other signs of genetic drift. This is likely caused by the later founding of the northern populations, together with subsequent founder and bottleneck effects, and a smaller population size. Furthermore, our data suggest a small eastern contribution among the Finns, consistent with the historical and linguistic background of the population.

Significance

Our results warn against a priori assumptions of homogeneity among Finns and other seemingly isolated populations. Thus, in association studies in such populations, additional caution for population structure may be necessary. Our results illustrate that population history is often important for patterns of genetic variation, and that the analysis of hundreds of thousands of SNPs provides high resolution also for population genetics.     

Finding Markers That Make a Difference: DNA Pooling and SNP-Arrays Identify Population Informative Markers for Genetic Stock Identification

Genetic stock identification (GSI) using molecular markers is an important tool for management of migratory species. Here, we tested a cost-effective alternative to individual genotyping, known as allelotyping, for identification of highly informative SNPs for accurate genetic stock identification. We estimated allele frequencies of 2880 SNPs from DNA pools of 23 Atlantic salmon populations using Illumina SNP-chip. We evaluated the performance of four common strategies (global F _ST, pairwise F _ST, Delta and outlier approach) for selection of the most informative set of SNPs and tested their effectiveness for GSI compared to random sets of SNP and microsatellite markers. For the majority of cases, SNPs selected using the outlier approach performed best followed by pairwise F _ST and Delta methods. Overall, the selection procedure reduced the number of SNPs required for accurate GSI by up to 53% compared with randomly chosen SNPs. However, GSI accuracy was more affected by populations in the ascertainment group rather than the ranking method itself. We demonstrated for the first time the compatibility of different large-scale SNP datasets by compiling the largest population genetic dataset for Atlantic salmon to date. Finally, we showed an excellent performance of our top SNPs on an independent set of populations covering the main European distribution range of Atlantic salmon. Taken together, we demonstrate how combination of DNA pooling and SNP arrays can be applied for conservation and management of salmonids as well as other species.     

Estimates of Genetic Differentiation Measured by F(ST) Do Not Necessarily Require Large Sample Sizes When Using Many SNP Markers

Population genetic studies provide insights into the evolutionary processes that influence the distribution of sequence variants within and among wild populations. F_ST is among the most widely used measures for genetic differentiation and plays a central role in ecological and evolutionary genetic studies. It is commonly thought that large sample sizes are required in order to precisely infer F_ST and that small sample sizes lead to overestimation of genetic differentiation. Until recently, studies in ecological model organisms incorporated a limited number of genetic markers, but since the emergence of next generation sequencing, the panel size of genetic markers available even in non-reference organisms has rapidly increased. In this study we examine whether a large number of genetic markers can substitute for small sample sizes when estimating F_ST. We tested the behavior of three different estimators that infer F_ST and that are commonly used in population genetic studies. By simulating populations, we assessed the effects of sample size and the number of markers on the various estimates of genetic differentiation. Furthermore, we tested the effect of ascertainment bias on these estimates. We show that the population sample size can be significantly reduced (as small as n = 4–6) when using an appropriate estimator and a large number of bi-allelic genetic markers (k>1,000). Therefore, conservation genetic studies can now obtain almost the same statistical power as studies performed on model organisms using markers developed with next-generation sequencing.     

Tertiary Origin and Pleistocene Diversification of Dragon Blood Tree (Dracaena cambodiana-Asparagaceae) Populations in the Asian Tropical Forests

Background

The origin of extraordinarily rich biodiversity in tropical forests is often attributed to evolution under stable climatic conditions over a long period or to climatic fluctuations during the recent Quaternary period. Here, we test these two hypotheses using Dracaena cambodiana, a plant species distributed in paleotropical forests.

Methods

We analyzed nucleotide sequence data of two chloroplast DNA (cpDNA: atpB-rbcL and trnD-trnT) regions and genotype data of six nuclear microsatellites from 15 populations (140 and 363 individuals, respectively) distributed in Indochina Peninsular and Hainan Island to infer the patterns of genetic diversity and phylogeographic structure. The population bottleneck and genetic drift were estimated based upon nuclear microsatellites data using the software programs BOTTLENECK and 2MOD. The lineage divergence times and past population dynamics based on cpDNA data were estimated using coalescent-based isolation-with-migration (IMa) and BEAST software programs.

Results

A significant phylogeographic structure (N _ST = 0.876, G _ST = 0.796, F _ST-SSR =0.329, R _ST =0.449; N _ST>G _ST, R _ST>F _ST-SSR, P<0.05) and genetic differentiation among populations were detected. Bottleneck analyses and Bayesian skyline plot suggested recent population reduction. The cpDNA haplotype network revealed the ancestral populations from the southern Indochina region expanded to northward. The most recent ancestor divergence time of D. cambodiana dated back to the Tertiary era and rapid diversification of terminal lineages corresponded to the Quaternary period.

Conclusions

The results indicated that the present distribution of genetic diversity in D. cambodiana was an outcome of Tertiary dispersal and rapid divergence during the Quaternary period under limited gene flow influenced by the uplift of Himalayan-Tibetan Plateau and Quaternary climatic fluctuations respectively. Evolutionary processes, such as extinction-recolonization during the Pleistocene may have contributed to the fast diversification in D. cambodiana.     

Tracking the emergence of a new breed using 49,034 SNP in sheep

Domestic animals are unique in that they have been organised into managed populations called breeds. The strength of genetic divergence between breeds may vary dependent on the age of the breed, the scenario under which it emerged and the strength of reproductive isolation it has from other breeds. In this study, we investigated the Gulf Coast Native breed of sheep to determine if it contains lines of animals that are sufficiently divergent to be considered separate breeds. Allele sharing and principal component analysis (PCA) using nearly 50,000 SNP loci revealed a clear genetic division that corresponded with membership of either the Florida or Louisiana Native lines. Subsequent analysis aimed to determine if the strength of the divergence exceeded that found between recognised breed pairs. Genotypes from 14 breeds sampled from Europe and Asia were used to obtain estimates of pair-wise population divergence measured as F _ST. The divergence separating the Florida and Louisiana Native (F _ST = 6.2%) was approximately 50% higher than the average divergence separating breeds developed within the same region of Europe (F _ST = 4.2%). This strongly indicated that the two Gulf Coast Native lines are sufficiently different to be considered separate breeds. PCA using small SNP sets successfully distinguished between the Florida and Louisiana Native animals, suggesting that allele frequency differences have accumulated across the genome. This is consistent with a population history involving geographic separation and genetic drift. Suggestive evidence was detected for divergence at the poll locus on sheep chromosome 10; however drift at neutral markers has been the largest contributor to the genetic separation observed. These results document the emergence of populations that can be considered separate breeds, with practical consequences for bio-conservation priorities, animal registration and the establishment of separate breed societies.     

Environmental and climatic determinants of molecular diversity and genetic population structure in a coenagrionid damselfly

Identifying environmental factors that structure intraspecific genetic diversity is of interest for both habitat preservation and biodiversity conservation. Recent advances in statistical and geographical genetics make it possible to investigate how environmental factors affect geographic organisation and population structure of molecular genetic diversity within species. Here we present a study on a common and wide ranging insect, the blue tailed damselfly Ischnuraelegans, which has been the target of many ecological and evolutionary studies. We addressed the following questions: (i) Is the population structure affected by longitudinal or latitudinal gradients?; (ii) Do geographic boundaries limit gene flow?; (iii) Does geographic distance affect connectivity and is there a signature of past bottlenecks?; (iv) Is there evidence of a recent range expansion and (vi) what is the effect of geography and climatic factors on population structure? We found low to moderate genetic sub-structuring between populations (mean F_ST = 0.06, D_est = 0.12), and an effect of longitude, but not latitude, on genetic diversity. No significant effects of geographic boundaries (e.g. water bodies) were found. F_ST-and D_est-values increased with geographic distance; however, there was no evidence for recent bottlenecks. Finally, we did not detect any molecular signatures of range expansions or an effect of geographic suitability, although local precipitation had a strong effect on genetic differentiation. The population structure of this small insect has probably been shaped by ecological factors that are correlated with longitudinal gradients, geographic distances, and local precipitation. The relatively weak global population structure and high degree of genetic variation within populations suggest that I. elegans has high dispersal ability, which is consistent with this species being an effective and early coloniser of new habitats.     

Environmental versus Anthropogenic Effects on Population Adaptive Divergence in the Freshwater Snail Lymnaea stagnalis

Repeated pesticide contaminations of lentic freshwater systems located within agricultural landscapes may affect population evolution in non-target organisms, especially in species with a fully aquatic life cycle and low dispersal ability. The issue of evolutionary impact of pollutants is therefore conceptually important for ecotoxicologists. The impact of historical exposure to pesticides on genetic divergence was investigated in the freshwater gastropod Lymnaea stagnalis, using a set of 14 populations from contrasted environments in terms of pesticide and other anthropogenic pressures. The hypothesis of population adaptive divergence was tested on 11 life-history traits, using Q _ST -F _ST comparisons. Despite strong neutral differentiation (mean F _ST = 0.291), five adult traits or parameters were found to be under divergent selection. Conversely, two early expressed traits showed a pattern consistent with uniform selection or trait canalization, and four adult traits appeared to evolve neutrally. Divergent selection patterns were mostly consistent with a habitat effect, opposing pond to ditch and channel populations. Comparatively, pesticide and other human pressures had little correspondence with evolutionary patterns, despite hatching rate impairment associated with global anthropogenic pressure. Globally, analyses revealed high genetic variation both at neutral markers and fitness-related traits in a species used as model in ecotoxicology, providing empirical support for the need to account for genetic and evolutionary components of population response in ecological risk assessment.