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

MYB59 transcription factor behaves differently in metallicolous and non-metallicolous populations of Arabidopsis halleri

In Arabidopsis thaliana (L.) Heynh., MYB59 transcription factor participates in regulating Ca homeostasis and signal transduction and is induced by Cd excess. To investigate its role in the facultative metallophyte Arabidopsis halleri ssp. halleri (L.) O'Kane and Al-Shehbaz, MYB59 expression was investigated under Cd treatment or Ca depletion in three populations belonging to distinct phylogeographic units (metallicolous PL22 and I16 and non-metallicolous I29), and compared with the expression in A. thaliana. In control conditions, MYB59 transcription in A. thaliana and the non-metallicolous population I29 follow a comparable trend with higher expression in roots than shoots, whereas in metallicolous populations I16 and PL22 its expression is similar in roots and shoots, suggesting a convergent evolution associated with adaptation to metalliferous environments. After 6 h of Ca depletion, MYB59 transcript levels were very high in I16 and PL22 populations, indicating that the adaptation to metalliferous environments requires tightly regulated Ca homeostasis and signalling. Cd treatment caused variability in MYB59 expression. In I29, MYB59 expression, as in A. thaliana, is likely associated to stress response, whereas its modulation in the two metallicolous populations reflects the different strategies for Cd tolerance and accumulation. In conclusion, MYB59 regulation in A. halleri is part of the network linking mineral nutrition and Cd tolerance/accumulation.

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.

Local adaptation (mostly) remains local: reassessing environmental associations of climate-related candidate SNPs in Arabidopsis halleri

Numerous landscape genomic studies have identified single-nucleotide polymorphisms (SNPs) and genes potentially involved in local adaptation. Rarely, it has been explicitly evaluated whether these environmental associations also hold true beyond the populations studied. We tested whether putatively adaptive SNPs in Arabidopsis halleri (Brassicaceae), characterized in a previous study investigating local adaptation to a highly heterogeneous environment, show the same environmental associations in an independent, geographically enlarged set of 18 populations. We analysed new SNP data of 444 plants with the same methodology (partial Mantel tests, PMTs) as in the original study and additionally with a latent factor mixed model (LFMM) approach. Of the 74 candidate SNPs, 41% (PMTs) and 51% (LFMM) were associated with environmental factors in the independent data set. However, only 5% (PMTs) and 15% (LFMM) of the associations showed the same environment-allele relationships as in the original study. In total, we found 11 genes (31%) containing the same association in the original and independent data set. These can be considered prime candidate genes for environmental adaptation at a broader geographical scale. Our results suggest that selection pressures in highly heterogeneous alpine environments vary locally and signatures of selection are likely to be population-specific. Thus, genotype-by-environment interactions underlying adaptation are more heterogeneous and complex than is often assumed, which might represent a problem when testing for adaptation at specific loci.

Estimating genomic diversity and population differentiation - an empirical comparison of microsatellite and SNP variation in Arabidopsis halleri

BACKGROUND

Microsatellite markers are widely used for estimating genetic diversity within and differentiation among populations. However, it has rarely been tested whether such estimates are useful proxies for genome-wide patterns of variation and differentiation. Here, we compared microsatellite variation with genome-wide single nucleotide polymorphisms (SNPs) to assess and quantify potential marker-specific biases and derive recommendations for future studies. Overall, we genotyped 180 Arabidopsis halleri individuals from nine populations using 20 microsatellite markers. Twelve of these markers were originally developed for Arabidopsis thaliana (cross-species markers) and eight for A. halleri (species-specific markers). We further characterized 2 million SNPs across the genome with a pooled whole-genome re-sequencing approach (Pool-Seq).

RESULTS

Our analyses revealed that estimates of genetic diversity and differentiation derived from cross-species and species-specific microsatellites differed substantially and that expected microsatellite heterozygosity (SSR-H e) was not significantly correlated with genome-wide SNP diversity estimates (SNP-H e and θ Watterson) in A. halleri. Instead, microsatellite allelic richness (A r) was a better proxy for genome-wide SNP diversity. Estimates of genetic differentiation among populations (F ST) based on both marker types were correlated, but microsatellite-based estimates were significantly larger than those from SNPs. Possible causes include the limited number of microsatellite markers used, marker ascertainment bias, as well as the high variance in microsatellite-derived estimates. In contrast, genome-wide SNP data provided unbiased estimates of genetic diversity independent of whether genome- or only exome-wide SNPs were used. Further, we inferred that a few thousand random SNPs are sufficient to reliably estimate genome-wide diversity and to distinguish among populations differing in genetic variation.

CONCLUSIONS

We recommend that future analyses of genetic diversity within and differentiation among populations use randomly selected high-throughput sequencing-based SNP data to draw conclusions on genome-wide diversity patterns. In species comparable to A. halleri, a few thousand SNPs are sufficient to achieve this goal.

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.

Plant Defensin type 1 (PDF1): protein promiscuity and expression variation within the Arabidopsis genus shed light on zinc tolerance acquisition in Arabidopsis halleri

Plant defensins are recognized for their antifungal properties. However, a few type 1 defensins (PDF1s) were identified for their cellular zinc (Zn) tolerance properties after a study of the metal extremophile Arabidopsis halleri. In order to investigate whether different paralogues would display specialized functions, the A. halleri PDF1 family was characterized at the functional and genomic levels. Eleven PDF1s were isolated from A. halleri. Their ability to provide Zn tolerance in yeast cells, their activity against Fusarium oxysporum f. sp. melonii, and their level of expression in planta were compared with those of the seven A. thaliana PDF1s. The genomic organization of the PDF1 family was comparatively analysed within the Arabidopsis genus. AhPDF1s and AtPDF1s were able to confer Zn tolerance and AhPDF1s also displayed antifungal activity. PDF1 transcripts were constitutively more abundant in A. halleri than in A. thaliana. Within the Arabidopsis genus, the PDF1 family is evolutionarily dynamic, in terms of gain and loss of gene copy. Arabidopsis halleri PDF1s display no superior abilities to provide Zn tolerance. A constitutive increase in AhPDF1 transcript accumulation is proposed to be an evolutionary innovation co-opting the promiscuous PDF1 protein for its contribution to Zn tolerance in A. halleri.

The Influence of Cadmium Stress on the Content of Mineral Nutrients and Metal-Binding Proteins in Arabidopsis halleri

We investigated the influence of cadmium stress on zinc hyperaccumulation, mineral nutrient uptake, and the content of metal-binding proteins in Arabidopsis halleri. The experiments were carried out using plants subjected to long-term cadmium exposure (40 days) in the concentrations of 45 and 225 μM Cd(2+). Inductively coupled plasma-mass spectrometry, size exclusion chromatography coupled with plasma-mass spectrometry, and laser ablation inductively coupled plasma-mass spectrometry used for ablation of polyacylamide gels were employed to assess the content of investigated elements in plants as well as to identify metal-binding proteins. We found that A. halleri is able to translocate cadmium to the aerial parts in high amounts (translocation index >1). We showed that Zn content in plants decreased significantly with the increase of cadmium content in the growth medium. Different positive and negative correlations between Cd content and mineral nutrients were evidenced by our study. We identified more than ten low-molecular-weight (<100 kDa) Cd-binding proteins in Cd-treated plants. These proteins are unlikely to be phytochelatins or metallothioneins. We hypothesize that low-molecular-weight Cd-binding proteins can be involved in cadmium resistance in A. halleri. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11270-012-1292-4) contains supplementary material, which is available to authorized users.

Allelic heterogeneity and trade-off shape natural variation for response to soil micronutrient

As sessile organisms, plants have to cope with diverse environmental constraints that may vary through time and space, eventually leading to changes in the phenotype of populations through fixation of adaptive genetic variation. To fully comprehend the mechanisms of evolution and make sense of the extensive genotypic diversity currently revealed by new sequencing technologies, we are challenged with identifying the molecular basis of such adaptive variation. Here, we have identified a new variant of a molybdenum (Mo) transporter, MOT1, which is causal for fitness changes under artificial conditions of both Mo-deficiency and Mo-toxicity and in which allelic variation among West-Asian populations is strictly correlated with the concentration of available Mo in native soils. In addition, this association is accompanied at different scales with patterns of polymorphisms that are not consistent with neutral evolution and show signs of diversifying selection. Resolving such a case of allelic heterogeneity helps explain species-wide phenotypic variation for Mo homeostasis and potentially reveals trade-off effects, a finding still rarely linked to fitness.

Plants are studied for their ability to adapt to their environment and especially to the physical constraints to which they are subjected. It is expected that they evolve in promoting genetic variants favorable under their native conditions, which could lead to negative consequences in other conditions. One approach to study the mechanisms and dynamics of these adaptations is to discover genetic variants that control potentially adaptive traits, and to study directly these variants in wild populations to try to reveal their evolutionary trajectory. We have identified a new polymorphism in a gene coding for a transporter of molybdenum (an essential micronutrient for the plant) in Arabidopsis; we show that this variant has strong phenotypic consequences at the level of plant growth and reproductive value in specific conditions, and that it explains a lot of the species diversity for these traits. Especially, the variant is associated with a clear negative effect under molybdenum-deficient conditions (caused by soil acidity) and with a subtle positive effect under molybdenum-plethoric conditions. Interestingly, the landscape distribution of the variant is not random among Asian populations and correlates well with the availability of molybdenum in the soil at the precise location where the plants are growing in the wild.