Genetic architecture of a plant adaptive trait: QTL mapping of intraspecific variation for tolerance to metal pollution in Arabidopsis halleri.

Anthropogenic activities are among the main drivers of global change and result in drastic habitat modifications, which represent strong evolutionary challenges for biological species that can either migrate, adapt, or disappear. In this context, understanding the genetics of adaptive traits is a prerequisite to enable long-term maintenance of populations under strong environmental constraints. To examine these processes, a QTL approach was developed here using the pseudometallophyte Arabidopsis halleri, which displays among-population adaptive divergence for tolerance to metallic pollution in soils. An F2 progeny was obtained by crossing individuals from metallicolous and non-metallicolous populations from Italian Alps, where intense metallurgic activities have created strong landscape heterogeneity. Then, we combined genome de novo assembly and genome resequencing of parental genotypes to obtain single-nucleotide polymorphism markers and achieve high-throughput genotyping of the progeny. QTL analysis was performed using growth parameters and photosynthetic yield to assess zinc tolerance levels. One major QTL was identified for photosynthetic yield. It explained about 27% of the phenotypic variance. Functional annotation of the QTL and gene expression analyses highlighted putative candidate genes. Our study represents a successful approach combining evolutionary genetics and advanced molecular tools, helping to better understand how a species can face new selective pressures of anthropogenic origin.

Evolutionary dynamics of quantitative variation in an adaptive trait at the regional scale: The case of zinc hyperaccumulation in Arabidopsis halleri.

Metal hyperaccumulation in plants is an ecological trait whose biological significance remains debated, in particular because the selective pressures that govern its evolutionary dynamics are complex. One of the possible causes of quantitative variation in hyperaccumulation may be local adaptation to metalliferous soils. Here, we explored the population genetic structure of Arabidopsis halleri at fourteen metalliferous and nonmetalliferous sampling sites in southern Poland. The results were integrated with a quantitative assessment of variation in zinc hyperaccumulation to trace local adaptation. We identified a clear hierarchical structure with two distinct genetic groups at the upper level of clustering. Interestingly, these groups corresponded to different geographic subregions, rather than to ecological types (i.e., metallicolous vs. nonmetallicolous). Also, approximate Bayesian computation analyses suggested that the current distribution of A. halleri in southern Poland could be relictual as a result of habitat fragmentation caused by climatic shifts during the Holocene, rather than due to recent colonization of industrially polluted sites. In addition, we find evidence that some nonmetallicolous lowland populations may have actually derived from metallicolous populations. Meanwhile, the distribution of quantitative variation in zinc hyperaccumulation did separate metallicolous and nonmetallicolous accessions, indicating more recent adaptive evolution and diversifying selection between metalliferous and nonmetalliferous habitats. This suggests that zinc hyperaccumulation evolves both ways-towards higher levels at nonmetalliferous sites and lower levels at metalliferous sites. Our results open a new perspective on possible evolutionary relationships between A. halleri edaphic types that may inspire future genetic studies of quantitative variation in metal hyperaccumulation.

Can zinc pollution promote adaptive evolution in plants? Insights from a one-generation selection experiment.

Human activities generate environmental stresses that can lead plant populations to become extinct. Population survival would require the evolution of adaptive responses that increase tolerance to these stresses. Thus, in pseudometallophyte species that have colonized anthropogenic metalliferous habitats, the evolution of increased metal tolerance is expected in metallicolous populations. However, the mechanisms by which metal tolerance evolves remain unclear. In this study, parent populations were created from non-metallicolous families of Noccaea caerulescens. They were cultivated for one generation in mesocosms and under various levels of zinc (Zn) contamination to assess whether Zn in soil represents a selective pressure. Individual plant fitness estimates were used to create descendant populations, which were cultivated in controlled conditions with moderate Zn contamination to test for adaptive evolution in functional traits. The number of families showing high fitness estimates in mesocosms was progressively reduced with increasing Zn levels in soil, suggesting increasing selection for metal tolerance. In the next generation, adaptive evolution was suggested for some physiological and ecological traits in descendants of the most exposed populations, together with a significant decrease of Zn hyperaccumulation. Our results confirm experimentally that Zn alone can be a significant evolutionary pressure promoting adaptive divergence among populations.

Demographic history of the trace metal hyperaccumulator Noccaea caerulescens (J. Presl and C. Presl) F. K. Mey. in Western Europe.

Noccaea caerulescens (Brassicaceae) is a major pseudometallophyte model for the investigation of the genetics and evolution of metal hyperaccumulation in plants. We studied the population genetics and demographic history of this species to advance the understanding of among-population differences in metal hyperaccumulation and tolerance abilities. Sampling of seven to 30 plants was carried out in 62 sites in Western Europe. Genotyping was carried out using a combination of new chloroplast and nuclear neutral markers. A strong genetic structure was detected, allowing the definition of three genetic subunits. Subunits showed a good geographic coherence. Accordingly, distant metallicolous populations generally belonged to distinct subunits. Approximate Bayesian computation analysis of demographic scenarios among subunits further supported a primary isolation of populations from the southern Massif Central prior to last glacial maximum, whereas northern populations may have derived during postglacial recolonization events. Estimated divergence times among subunits were rather recent in comparison with the species history, but certainly before the establishment of anthropogenic metalliferous sites. Our results suggest that the large-scale genetic structure of N. caerulescens populations pre-existed to the local adaptation to metalliferous sites. The population structure of quantitative variation for metal-related adaptive traits must have established independently in isolated gene pools. However, features of the most divergent genetic unit (e.g. extreme levels of Cd accumulation observed in previous studies) question the putative relationships between adaptive evolution of metal-related traits and subunits isolation. Finally, admixture signals among distant metallicolous populations suggest a putative role of human activities in facilitating long-distance genetic exchanges.

Potential preadaptation to anthropogenic pollution: evidence from a common quantitative trait locus for zinc and cadmium tolerance in metallicolous and nonmetallicolous accessions of Arabidopsis halleri.

As a drastic environmental change, metal pollution may promote the rapid evolution of genetic adaptations contributing to metal tolerance. In Arabidopsis halleri, genetic bases of zinc (Zn) and cadmium (Cd) tolerance have been uncovered only in a metallicolous accession, although tolerance is species-wide. The genetic determinants of Zn and Cd tolerance in a nonmetallicolous accession were thus investigated for the first time. The genetic architecture of tolerance was investigated in a nonmetallicolous population (SK2) by using first backcross progeny obtained from crosses between SK2 and Arabidopsis lyrata petraea, a nonmetallophyte species. Only one significant and common quantitative trait locus (QTL) region was identified explaining 22.6% and 31.2% of the phenotypic variation for Zn and Cd tolerance, respectively. This QTL co-localized with HEAVY METAL ATPASE 4 (AhHMA4), which was previously validated as a determinant of Zn and Cd tolerance in a metallicolous accession. Triplication and high expression of HMA4 were confirmed in SK2. In contrast, gene duplication and high expression of METAL TOLERANT PROTEIN 1A (MTP1A), which was previously associated with Zn tolerance in a metallicolous accession, were not observed in SK2. Overall, the results support the role of HMA4 in tolerance capacities of A. halleri that may have pre-existed in nonmetallicolous populations before colonization of metal-polluted habitats. Preadaptation to metal-contaminated sites is thus discussed.

Recent and ancient signature of balancing selection around the S-locus in Arabidopsis halleri and A. lyrata.

Balancing selection can maintain different alleles over long evolutionary times. Beyond this direct effect on the molecular targets of selection, balancing selection is also expected to increase neutral polymorphism in linked genome regions, in inverse proportion to their genetic map distances from the selected sites. The genes controlling plant self-incompatibility are subject to one of the strongest forms of balancing selection, and they show clear signatures of balancing selection. The genome region containing those genes (the S-locus) is generally described as nonrecombining, and the physical size of the region with low recombination has recently been established in a few species. However, the size of the region showing the indirect footprints of selection due to linkage to the S-locus is only roughly known. Here, we improved estimates of this region by surveying synonymous polymorphism and estimating recombination rates at 12 flanking region loci at known physical distances from the S-locus region boundary, in two closely related self-incompatible plants Arabidopsis halleri and A. lyrata. In addition to studying more loci than previous studies and using known physical distances, we simulated an explicit demographic scenario for the divergence between the two species, to evaluate the extent of the genomic region whose diversity departs significantly from neutral expectations. At the closest flanking loci, we detected signatures of both recent and ancient indirect effects of selection on the S-locus flanking genes, finding ancestral polymorphisms shared by both species, as well as an excess of derived mutations private to either species. However, these effects are detected only in a physically small region, suggesting that recombination in the flanking regions is sufficient to quickly break up linkage disequilibrium with the S-locus. Our approach may be useful for distinguishing cases of ancient versus recently evolved balancing selection in other systems.

Investigation of the geographical scale of adaptive phenological variation and its underlying genetics in Arabidopsis thaliana.

Despite the increasing number of genomic tools, identifying the genetics underlying adaptive complex traits remains challenging in the model species Arabidopsis thaliana. This is due, at least in part, to the lack of data on the geographical scale of adaptive phenotypic variation. The aims of this study were (i) to tease apart the historical roles of adaptive and nonselective processes in shaping phenological variation in A. thaliana in France and (ii) to gain insights into the spatial scale of adaptive variation by identifying the putative selective agents responsible for this selection. Forty-nine natural stands from four climatically contrasted French regions were characterized (i) phenologically for six traits, (ii) genetically using 135 SNP markers and (iii) ecologically for 42 variables. Up to 63% of phenological variation could be explained by neutral genetic diversity. The remaining phenological variation displayed stronger associations with ecological variation within regions than among regions, suggesting the importance of local selective agents in shaping adaptive phenological variation. Although climatic conditions have often been suggested as the main selective agents acting on phenology in A. thaliana, both edaphic conditions and interspecific competition appear to be strong selective agents in some regions. In a first attempt to identify the genetics of phenological variation at different geographical scales, we phenotyped worldwide accessions and local polymorphic populations from the French RegMap in a genome-wide association (GWA) mapping study. The genomic regions associated with phenological variation depended upon the geographical scale considered, stressing the need to account for the scale of adaptive phenotypic variation when choosing accession panels for GWAS.

Nuclear and chloroplast DNA phylogeography reveals vicariance among European populations of the model species for the study of metal tolerance, Arabidopsis halleri (Brassicaceae).

Arabidopsis halleri is a pseudometallophyte involved in numerous molecular studies of the adaptation to anthropogenic metal stress. In order to test the representativeness of genetic accessions commonly used in these studies, we investigated the A. halleri population genetic structure in Europe. Microsatellite and nucleotide polymorphisms from the nuclear and chloroplast genomes, respectively, were used to genotype 65 populations scattered over Europe. The large-scale population structure was characterized by a significant phylogeographic signal between two major genetic units. The localization of the phylogeographic break was assumed to result from vicariance between large populations isolated in southern and central Europe, on either side of ice sheets covering the Alps during the Quaternary ice ages. Genetic isolation was shown to be maintained in western Europe by the high summits of the Alps, whereas admixture was detected in the Carpathians. Considering the phylogeographic literature, our results suggest a distinct phylogeographic pattern for European species occurring in both mountain and lowland habitats. Considering the evolution of metal adaptation in A. halleri, it appears that recent adaptations to anthropogenic metal stress that have occurred within either phylogeographic unit should be regarded as independent events that potentially have involved the evolution of a variety of genetic mechanisms.

Nuclear microsatellite loci for Arabidopsis halleri (Brassicaceae), a model species to study plant adaptation to heavy metals.

PREMISE OF THE STUDY: Arabidopsis halleri is a model species to study the adaptation of plants to soils contaminated by zinc, cadmium, and lead. To provide a neutral genetic background with which adaptive genetic markers could be compared, we developed highly polymorphic neutral microsatellite markers. METHODS AND RESULTS: Using a microsatellite-enriched library method, we identified 120 microsatellite loci for quantitative trait locus (QTL) mapping analysis, of which eight primer pairs were developed in a single multiplex for population genetic studies. Analyses were performed on 508 individuals from 26 populations. All loci were polymorphic with six to 23 alleles per locus. Genetic diversity varied between 0.56 and 0.76. CONCLUSIONS: Our results demonstrated the value of these eight microsatellite markers to investigate neutral population genetic structure in A. halleri. To increase the resolution of population genetic analyses, we suggest adding them to the 11 markers previously developed independently.

Data on the identification of microsatellite markers in Eisenia fetida and Eisenia andrei.

Eisenia fetida and Eisenia andrei are closely related earthworm species that play a crucial part in soil and influence its structure and organic matter cycling. Due to their essential environmental role, they are widely used as model organisms in a vast spectrum of research areas. In this work, we partially sequenced genomes of E. fetida and E. andrei, using Illumina technology (Nano 2 × 250 v2 - MiSeq) and de novo assembly strategy. A total of 3785 and 4258 microsatellite or Simple Sequence Repeat (SSR) markers were identified within E. fetida and E. andrei genomic DNA, respectively. The microsatellite markers will facilitate the analyses of genetic diversity and population genetics studies for the two selected earthworm species and their interspecific hybrids.

Sublethal effects of metal toxicity and the measure of plant fitness in ecotoxicological experiments.

Anthropogenic pollution is a major driver of global environmental change. To be properly addressed, the study of the impact of pollutants must consider both lethal effects and sublethal effects on individual fitness. However, measuring fitness remains challenging. In plants, the total number of seeds produced, i.e. the seed set, is traditionally considered, but is not readily accessible. Instead, performance traits related to survival, e.g., vegetative biomass and reproductive success, can be measured, but their correlation with seed set has rarely been investigated. To develop accurate estimates of seed set, relationships among 15 vegetative and reproductive traits were analyzed. For this purpose, Noccaea caerulescens (Brassicaceae), a model plant to study local adaptation to metal-contaminated environments, was used. To investigate putative variation in trait relationships, sampling included several accessions cultivated in contrasting experimental conditions. To test their applicability, selected estimates were used in the first generation of a Laboratory Natural Selection (LNS) experiment exposing experimentally plants to zinc soil pollution. Principal component analyses revealed statistical independence between vegetative and reproductive traits. Traits showing the strongest positive correlation with seed set were the number of non-aborted silicles, and the product of this number and mean silicle length. They thus appeared the most appropriate to document sublethal or fitness effects of environmental contaminants in plant ecotoxicological studies. The relevance of both estimates was confirmed by using them to assess the fitness of parental plants of the first generation of an LNS experiment: the same families consistently displayed the highest or the lowest performance values in two independent experimental metal-exposed populations. Thus, both these fitness estimates could be used to determine the expected number of offspring and the composition of successive generations in further LNS experiments investigating the impact of multi-generational exposure of a plant species to environmental pollution.

Polymorphic microsatellite markers demonstrate hybridization and interspecific gene flow between lumbricid earthworm species, Eisenia andrei and E. fetida.

The lumbricid earthworms Eisenia andrei (Ea) and E. fetida (Ef) have been used as model organisms for studies on hybridization. Previously they have been identified by species specific sequences of the mitochondrial COI gene of maternal origin ('a' or 'f') and the nuclear 28S gene of maternal/paternal origin ('A' or 'F'). In experimental crosses, these hermaphroditic species produce progeny of genotypes Ea (aAA), Ef (fFF) and hybrids (aAF and fFA) originating by self-fertilization or cross-fertilization. To facilitate studies on new aspects of the breeding biology and hybridization of earthworms, polymorphic microsatellite markers were developed based on 12 Ea and 12 Ef specimens and validated on DNA samples extracted from 24 genotyped specimens (aAA, fFF, aAF and fFA) from three laboratory-raised families and 10 of them were applied in the present study. The results indicate that microsatellite markers are valuable tools for tracking interspecific gene flow between these species.

Biochemical composition of Phleum pratense pollen grains: A review.

The Poaceae family is composed of 12,000 plant species. Some of these species produce highly allergenic anemophilous pollen grains (PGs). Phleum pratense pollen grains (PPPGs) emerged as a model for studies related to grass allergy. The biochemical composition of allergenic PGs has not yet been fully described despite potential health effects of PG constituents other than allergenic proteins. This review brings together the information available in literature aiming at creating a comprehensive picture of the current knowledge about the chemical composition of allergenic PGs from timothy grass. PPPGs have an average diameter between 30-35 µm and the mass of a single PG was reported between 11 and 26 ng. The pollen cytoplasm is filled with two types of pollen cytoplasmic granules (PCGs): the starch granules and the polysaccharide particles (p-particles). Starch granules have a size between 0.6-2.5 µm with an average diameter of 1.1 µm (estimated number of 1000 granules per PG) while p-particles have a size ranging around 0.3 to 0.4 µm (estimated number between 61,000-230,000 p-particles per PG). The rupture of PG induces the release of PCGs and the dispersion of allergens in the inhalable fraction of atmospheric aerosol. PPPGs are composed of sporopollenin, sugars, polysaccharides, starch, glycoproteins (including allergens), amino-acids, lipids, flavonoids (including isorhamnetin), various elements (the more abundant being Si, Mg and Ca), phenolic compounds, phytoprostanoids, carotenoids (pigments) metals and adsorbed pollutants. PPPG contains about a hundred different proteins with molecular masses ranging from 10 to 94 kDa, with isoelectric points from 3.5-10.6. Among these proteins, allergens are classified in eleven groups from 1 to 13 with allergens from groups 1 and 5 being the major contributors to Phl p pollen allergy. Major allergen Phl p 5 was quantified in PPPGs by several studies with concentration ranging from 2.7 and 3.5 µg.mg-1 in unpolluted environment. Values for other allergens are scarce in literature; only one quantitative assessment exists for allergen groups Phl p 1, 2 and 4. The extractible lipid fraction of PPPGs is estimated between 1.7-2.2% of the total PG mass. The main chemical families of lipids reported in PPPGs are: alkanes, alkenes, alcohols, saturated and unsaturated fatty acids, di- and tri-hydroxylated fatty acids, aldehydes and sterols. Several lipid compounds with potential adjuvant effects on allergy have been specifically quantified in PPPGs: E2-like prostaglandin (PGE2), B4-like leukotriene (LTB4), unsaturated fatty acids (linoleic and linolenic acids and their hydroxylated derivatives), adenosine, vitamins and phenolic compounds. Some other biochemical characteristics such as NAD(P)H oxidase, protease activity and pollen microbiome were described in the literature. The bioaccessibility in physiological conditions has not been described for most biochemicals transported by allergenic PPPGs. There is also a considerable lack of knowledge about the potential health effects of pollen constituents other than allergens. The variability of pollen composition remains also largely unknown despite its importance for plant reproduction and allergy in an environment characterized by chemical pollution, climate change and loss of biodiversity.

When population genetics serves genomics: putting adaptation back in a spatial and historical context.

Recent advances in molecular biology have opened new perspectives for the study of plant adaptation, especially at the intraspecific level. Nowadays, scientists employing -omic results in multiple scientific fields can be optimistic of their chances of revealing mechanisms involved in adaptive population divergence. However, the investment required by integrative studies greatly reduces the number of experiments that can be performed. In this context, a comprehensive choice of accessions under study is crucial. We maintain this choice could be appreciably enlightened by population genetics because it helps putting adaptive population divergence in a spatial and historical context. As an example, we highlight the usefulness of knowledge about population genetic structure in the integrative study of metal tolerance in Arabidopsis halleri.

Variability of zinc tolerance among and within populations of the pseudometallophyte species Arabidopsis halleri and possible role of directional selection.

We estimated the level of quantitative polymorphism for zinc (Zn) tolerance in neighboring metallicolous and nonmetallicolous populations of Arabidopsis halleri and tested the hypothesis that divergent selection has shaped this polymorphism. A short-term hydroponic test was used to capture the quantitative polymorphism present between edaphic types, among and within populations. We measured six morphological and physiological traits on shoots and roots to estimate the response of A. halleri to Zn. In order to assess the adaptive value of Zn tolerance polymorphism, we compared differentiation of quantitative traits with that of molecular markers. Zinc tolerance of metallicolous populations was, on average, higher than that of nonmetallicolous populations according to the morphological and physiological traits measured. Phenotypic variability within edaphic types was very high and mainly explained by polymorphism among individuals within populations. Genetic differentiation for photosystem II yield of leaves (a measure of photosynthetic efficiency) was greater than the differentiation for microsatellite and thus, probably shaped by divergent selection. Overall, these results suggest that, in the sampled populations, Zn tolerance has been increased in metallicolous populations through selection on standing genetic variation within local nonmetallicolous ancestral populations.

Impacts of metallic trace elements on an earthworm community in an urban wasteland: Emphasis on the bioaccumulation and genetic characteristics in Lumbricus castaneus.

Metallic trace elements (MTEs) soil pollution has become a worldwide concern, particularly regarding its impact on earthworms. Earthworms, which constitute the dominant taxon of soil macrofauna in temperate regions and are crucial ecosystem engineers, are in direct contact with MTEs. The impacts of MTE exposure on earthworms, however, vary by species, with some able to cope with high levels of contamination. We combined different approaches to study the effects of MTEs at different levels of biological organisation of an earthworm community, in a contaminated urban wasteland. Our work is based on field collection of soil and earthworm samples, with a total of 891 adult earthworms from 8 species collected, over 87 quadrats across the study plot. We found that MTE concentrations are highly structured at the plot scale and that some elements, such as Pb, Zn, and Cu, are highly correlated. Comparing species assemblage to MTE concentrations, we found that the juvenile and adult abundances, and community composition, were significantly affected by pollution. Along the pollution gradient, as species richness decreased, Lumbricus castaneus became more dominant. We thus investigated the physiological response of this species to a set of specific elements (Pb, Zn, Cu, and Cd) and studied the impacts of MTE concentrations at the plot scale on its population genetic. These analyses revealed that L. castaneus is able to bioaccumulate high quantities of Cd and Zn, but not of Cu and Pb. The population genetic analysis, based on the genotyping of 175 individuals using 8 microsatellite markers, provided no evidence of the role of the heterogeneity in MTE concentrations as a barrier to gene flow. The multidisciplinary approach we used enabled us to reveal the comparatively high tolerance of L. castaneus to MTE concentrations, suggesting that this is a promising model to study the molecular bases of MTE tolerance.

Differential retention of transposable element-derived sequences in outcrossing Arabidopsis genomes.

BACKGROUND: Transposable elements (TEs) are genomic parasites with major impacts on host genome architecture and host adaptation. A proper evaluation of their evolutionary significance has been hampered by the paucity of short scale phylogenetic comparisons between closely related species. Here, we characterized the dynamics of TE accumulation at the micro-evolutionary scale by comparing two closely related plant species, Arabidopsis lyrata and A. halleri. RESULTS: Joint genome annotation in these two outcrossing species confirmed that both contain two distinct populations of TEs with either 'recent' or 'old' insertion histories. Identification of rare segregating insertions suggests that diverse TE families contribute to the ongoing dynamics of TE accumulation in the two species. Orthologous TE fragments (i.e. those that have been maintained in both species), tend to be located closer to genes than those that are retained in one species only. Compared to non-orthologous TE insertions, those that are orthologous tend to produce fewer short interfering RNAs, are less heavily methylated when found within or adjacent to genes and these tend to have lower expression levels. These findings suggest that long-term retention of TE insertions reflects their frequent acquisition of adaptive roles and/or the deleterious effects of removing nearly neutral TE insertions when they are close to genes. CONCLUSION: Our results indicate a rapid evolutionary dynamics of the TE landscape in these two outcrossing species, with an important input of a diverse set of new insertions with variable propensity to resist deletion.

Phylogeography of Arabidopsis halleri (Brassicaceae) in mountain regions of Central Europe inferred from cpDNA variation and ecological niche modelling.

The present study aimed to investigate phylogeographical patterns present within A. halleri in Central Europe. 1,281 accessions sampled from 52 populations within the investigated area were used in the study of genetic variation based on chloroplast DNA. Over 500 high-quality species occurrence records were used in ecological niche modelling experiments. We evidenced the presence of a clear phylogeographic structure within A. halleri in Central Europe. Our results showed that two genetically different groups of populations are present in western and eastern part of the Carpathians. The hypothesis of the existence of a glacial refugium in the Western Carpathians adn the Bohemian Forest cannot be rejected from our data. It seems, however, that the evidence collected during the present study is not conclusive. The area of Sudetes was colonised after LGM probably by migrants from the Bohemian Forest.

Characterization and multiplexing of 21 microsatellite markers for the herb Noccaea caerulescens (Brassicaceae).

PREMISE OF THE STUDY: Multiplexed microsatellite markers were developed for population genetic studies in the pseudometallophyte Noccaea caerulescens (Brassicaceae), a model species to investigate metal tolerance and hyperaccumulation in higher plants. METHODS AND RESULTS: Microsatellite loci were isolated through pyrosequencing of an enriched DNA library. Three multiplexes combining four previously published and 17 newly designed markers were developed. The new markers were screened in metallicolous and nonmetallicolous populations from southern France. The total number of alleles per locus ranged from five to 18. The observed heterozygosity per locus and per population ranged from 0 to 0.83, and expected heterozygosity ranged from 0 to 0.89. CONCLUSIONS: The investigated loci showed reasonable to high levels of polymorphism at the regional scale. The multiplex set should be helpful in investigating genetic diversity, population structure, and demographic history in N. caerulescens at various spatial scales.

Using biomarkers in an evolutionary context: lessons from the analysis of biological responses of oligochaete annelids to metal exposure.

Anthropogenic activities may lead to the accumulation of inorganic and organic compounds in topsoils. Biota living in close contact with contaminated soils may experience stress at different levels of biological organization throughout the continuum from molecular to community level. Biological responses observed at the individual or infra-individual level of biological organization led to the development of biomarkers. The development of biomarkers consists often in evidencing biological modifications following a contaminant stress in laboratory conditions, using naïve organisms and it is sometime proposed to use the biological state of individuals from sentinel species collected in the field to evaluate the level of environmental exposure. However, considering the possibility of local adaptation following long-term exposure, organisms response sampled in the field may substantially differ from laboratory specimens. In this review, we discuss this point focusing on the definition and validity of molecular biomarkers of metal pollution using earthworms of the Lumbricidae family.