<title>Introduction</title> <p>Serpentine soil mosaics are a classic context for ecological adaptation ##UREF##0##[1]##, ##UREF##1##[2]##. These widespread soils occur in small patches along fault lines where igneous rocks such as serpentinite are exposed. This environment is characterized by a suite of challenging abiotic factors such as low calcium-to-magnesium ratios, increased heavy metal concentrations, nutrient deficiency, and low moisture retention ##UREF##2##[3]##, leading to sharp transitions in abiotic conditions at the boundaries of serpentine patches ##UREF##3##[4]##. The ecological community is differentiated along these boundaries as well: serpentine soils support an considerable array of endemic plant species, with serpentine specialists comprising 12.5% of native California plants, despite these soils comprising only 1% of the land area in the region ##UREF##4##[5]##. Some widespread species are also found on serpentine, with several documented cases of local adaptation to soil conditions (e.g. <italic>Collinsia</italic>
##UREF##5##[6]##, <italic>Cerastium</italic>
##UREF##6##[7]##, and <italic>Pinus</italic>
##UREF##7##[8]##). Populations of <italic>Arabidopsis lyrata</italic>, a perennial self incompatible crucifer, grow on serpentine outcroppings, interspersed with populations on other soil types such as granitic outcroppings and sand dunes ##UREF##8##[9]##. Though information regarding local adaptation is lacking for <italic>A. lyrata</italic>, its proximity to the genetic model organism <italic>A. thaliana</italic> provides an opportunity to locate polymorphisms which are associated with the serpentine soil habitat. The role of these polymorphisms in serpentine adaptation, if any, can then be experimentally investigated.</p> <p>By hybridizing genomic DNA from <italic>A. lyrata</italic> to <italic>A. thaliana</italic> Affymetrix tiling arrays, we can measure genetic differentiation between soil types at 2,853,369 probes throughout the <italic>A. lyrata</italic> genome. When DNA is hybridized to the array, probes which overlap a polymorphic SNP or indel will hybridize poorly in individuals with the mismatched allele ##REF##15687292##[10]##–##REF##16527929##[14]##; when hybridization intensity is significantly different between populations, a differentiated polymorphism can be mapped to a specific location in the genome ##REF##16076241##[15]##–##REF##18493034##[17]##. This method facilitates the discovery of adaptive variation in several complementary ways. First, we can discover candidate genes for adaptation to important environmental conditions, such as the low Ca∶Mg ratio which is characteristic of serpentine soil. In <italic>A. thaliana</italic> (which cannot normally survive on serpentine soils), an induced loss of function mutation in the calcium-proton antiporter <italic>cax1</italic> enhances survival on soils with a low Ca∶Mg ratio ##REF##15948832##[18]##: differentiation of natural variation at related genes in <italic>A. lyrata</italic> would provide clear candidate loci for adaptation to low Ca∶Mg ratio. Second, we can use the natural distribution of genetic variation between environments to investigate the functions of unannotated genes and non-coding features. If an uncharacterized locus consistently assorts with soil ion content, for example, then it can be hypothesized to interact with this environmental variable to determine fitness. Third, we can use genetic differentiation at genes with known function to form hypotheses about other environmental differences which may be important in nature. Although abiotic factors are thought to be the most important drivers of serpentine adaptation, if differentiation is also found in genes coding for anti-microbial compounds, anti-predatory compounds, or proteins which mediate competitive interactions, then parasitism, predation, or competition can be inferred to be important.</p> <p>In the current work, we have hybridized genomic DNA from two serpentine and two granitic populations of <italic>A. lyrata</italic> to the <italic>A. thaliana</italic> tiling array. This has allowed us to locate many polymorphisms which are differentiated between soil types, including excellent candidates for adaptation to soil conditions. We do not mean to imply that selection is the only force which leads to differentiation between populations, as stochastic demographic forces may also lead to correlated distributions of genetic polymorphisms and environmental conditions, especially in this small sample of populations ##REF##1213275##[19]##–##REF##17473869##[21]##. Indeed, there has long been a debate as to the relative importance of selection and demography in creating population differentiation ##REF##1213275##[19]##, ##REF##4711903##[22]##, ##REF##5803535##[23]##. This debate has only intensified in the genomic era, as it is now clear that selection can be rampant, and lead to adaptive change at loci across the genome ##REF##18302709##[24]##–##REF##18282109##[26]##. Our paper is not intended to resolve this debate–rather, we provide a genomic portrait of differentiation which is necessary but not sufficient to determine which polymorphisms are influenced by spatially varying selection. This data can then inform functional investigation of candidate genes and processes.</p>
<title>Materials and Methods</title> <title>Sample collection</title> <p>Seeds were collected at two serpentine sites and two granitic sites in close proximity (##FIG##0##figure 1##). Although historically there were 26 serpentinic outcroppings in an arc from Virginia to Staten Island, New York, nearly half have been destroyed by landuse changes in greater Philadelphia and Baltimore; of the remaining outcroppings, two are large (1000 acres or greater) and largely protected, while the others are small (2–10 acres) and altered by successional changes ##UREF##10##[37]##. Although a few of the smaller barrens do have small <italic>A. lyrata</italic> populations (frequently fewer than 100 individuals, EVW unpublished data), we restricted our analyses to the two patches which we estimate contained over 10,000 individuals. The two non-serpentinic populations sampled are the only ones approximately equidistant to the two serpentinic populations that we could find based on herbaria records (R. Latham, unpublished) and our own search; the other five sites with herbaria collections on non-serpentinic substrates we visited had disappeared due to ecological succession or human alteration. The serpentinic sites are part of the Nature Conservancy's state line serpentine barren complex in Chester County, Pa, and Soldier's Delight Natural Environment Area, Baltimore County, Md. Granitic sites were located at the Lock 12 recreation area on the Susquehanna river (York County, Pa), and Lancaster County Park (Lancaster County, Pa).</p> <title>Microarray hybridization</title> <p>DNA was extracted from three plants from each of the four populations using Qiagen plant mini-kits and amplified using rolling-circle amplification with the Qiagen Repli-G kit. After phenol-chloroform clean up and ethanol precipitation, 10 µg of DNA from each sample was fragmented with DNase following Turner <italic>et al.</italic> 2008. To fragment DNA, we created a master mix of DNase I (Promega), One-Phor-All buffer (Amersham Biosciences), and Acetylated BSA (Invitrogen); the amounts added per sample were 4 µl 10X One-Phor-All, 0.14 µl Acetylated BSA, and approximately 0.085 µl of DNase per ug of DNA. The precise amount of DNase added per sample depends on the batch of enzyme: our goal was to calibrate the reaction to obtain fragmented DNA of approximately 50 bp, with low variance in fragment size. Fragmentation was performed in a MJ Research PTC-200 thermocycler at 37 C for 16 min, 99 C for 15 min, and 12 C for 15 min. Samples with similar intensities and fragment sizes (approximate fragment size = 50 bp) were labelled with a master mix of Biotin-N6-ddATP (Enzo) and RTdT enzyme (Promega) to each sample. RTdT was diluted from 30 U/µl to 15 U/µl enzyme before use by mixing RTdT enzyme, RTdT 5X buffer, and water at a ratio of 5∶1∶4. For each sample to be labelled, 2 µl of Biotin and 3 µl of RTdT were added to the master mix. Labelling was accomplished in a thermocycler at 37 C for 90 min, 99 C for 15 min, and 12 C for 5 min. Microarrays were hybridized at the Affymetrix core facility in the UC Davis Genome Center using standard conditions for this array.</p> <title>Microarray analysis</title> <p>Using NCBI megablast, we determined that 2,853,369 array probes have a single perfect match in the TAIR7 <italic>A. thaliana</italic> reference genome and retained these probes for analysis (that is, we exclude probes with more than one exact match in the reference genome). We spatially normalized array intensities at these probes by dividing each array into 1600 subarrays of 64 by 64 probes, and divided the log(intensity) of each probe by the median log(intensity) of the local 64×64 probe window. We further normalized these values using quantile normalization in R ##REF##14960456##[38]##. Of these 2.85 million probes, 371,642 (13%) had a perfect match in a preliminary 4× draft of the <italic>A. lyrata</italic> genome (kindly provided by DOE-JG's Community Sequencing Program [proposal coordinated by Detlef Weigel, MPI Tübingen]).</p> <p>To find differentiated polymorphisms, <italic>t</italic>-tests were computed for each probe. We first analyzed only the probes with a match to <italic>A. lyrata</italic>, and corrected these <italic>p</italic>-values for 371,642 multiple tests. We used three significance thresholds in order to have inclusive and conservative lists of differentiated probes: a bonferroni-corrected <italic>p</italic>&lt;0.05, and two false discovery rate (FDR) thresholds: FDR&lt;0.01 and FDR&lt;0.05 ##REF##12883005##[39]##. Next, we computed <italic>t</italic>-test <italic>p</italic>-values for the remaining 2,481,727 probes, correcting these values for the greater number of multiple tests, and discovered additional probes significant at the three thresholds. We combined these two sets of significant probes to generate a list of differentiated loci at the three levels of stringency, which are available in ##SUPPL##0##data set S1##. Note that FDR calculations assume a uniform distribution of truly null features (<italic>p</italic>-value * number of multiple tests / <italic>p</italic>-value rank), and are therefore estimates. These estimates are based on technical, and not biological, false discovery: a probe with low specificity may truly detect differentiation, but not at the expected location.</p> <p>To increase our power to find differentiated loci, we used a sliding window analysis to determine if probes with low FDR are clustered in the genome. For this analysis, we only used probes with a single perfect match in the draft <italic>A. lyrata</italic> genome. There were 5,338 probes with FDR&lt;0.30 (<italic>p</italic>&lt;0.003): these probes have high individual rates of false discovery, but can be used as markers to find regions of the genome that are enriched for differentiated probes. Using permuted chromosomes to generate a null distribution, we mapped windows of 5, 10, and 20 probes that are significantly enriched for candidate markers. Each window recorded the average probe FDR, and compared this value to the permuted chromosomes. Probes with FDR&gt;0.30 were considered to have an FDR = 1.00, as this was previously found to make the analysis more conservative ##REF##18493034##[17]##. Because FDR estimates for windows were found to be less accurate than FDR estimates for individual probes ##REF##18493034##[17]##, we considered windows significant only if they had FDR&lt;0.001 when compared to the permuted chromosomes.</p> <title>Differentiation of functional categories</title> <p>Using the TAIR7 annotation of the <italic>A. lyrata</italic> reference genome, we determined which genes overlapped our differentiated probes and windows; a gene was considered to overlap a differentiated locus if any portion of its transcript overlaps a significant probe or window. Significance of gene ontology (GO) terms among our significant genes was determined in two complementary ways: a binomial sampling test, and a permutation test (GO terms are a controlled vocabulary which classifies genes based on all known functional associations determined computationally or experimentally). First, we compared the observed proportion of significant genes in each GO term to an expected number, which is the binomial sampling probability of sampling an equal or larger number of genes in each category given the number of significant genes in the genome. For the permutation test, we simulated our discovery procedure by randomly sampling probes and genomic regions of the same size and number as our significant probes and regions. For example, in the primary analysis presented below, there are 751 significant probes and 71 significant regions. To generate a null distribution of genes for this analysis, we sampled this same number of probes and regions (of an equal size distribution), and determined which genes overlap the random sample. We considered a GO term to be significantly overrepresented if the observed sample has more genes in the given term than 5% of 500 simulated data sets, and less than 5% probability of occurring by chance in the binomial sampling test.</p> <title>Copy number variants</title> <p>If serpentine and granitic populations are equidistant to the reference genome, then any given differentiated probe should have higher intensity on serpentine half the time, and higher intensity on granite half the time. In some significant regions, however, all probes are more intense on one soil type, consistent with differentiated copy number variants (CNVs). Differentiated regions are considered candidate CNVs if the ratio of mean hybridization intensity for that region is more extreme than expected. We extracted the raw (unnormalized) intensity values for all probes in each significant region, and computed a mean intensity for serpentine and granitic populations. These means were minimally normalized by dividing by the median intensity for all probes for each soil type, and the ratio of the normalized means was used as the test metric. An expected distribution of ratios was created by randomly sampling the number of probes in the test region from the set of probes from all significant regions. Regions were considered candidate CNVs if their ratio was more extreme than 9,990/10,000 random sets (<italic>p</italic>&lt;0.001). This permutation controls for many possible confounding effects. For example, if all derived SNPs were found on serpentine, then granitic populations would have higher average intensities in significant regions; because our permutation samples only from significant regions, the candidate CNV are significantly beyond any such effect.</p> <title>Sequencing and PCR</title> <p>Three loci were amplified via long PCR using Phusion polymerase (Finnzymes), and two were directly sequenced. Population genetic analysis of sequenced loci was done in DNAsp ##REF##14668244##[40]##; see ##TAB##0##table 1## for summary statistics. The third locus, <italic>AT5G27100</italic>, is polymorphically deleted. The presence/absence of this gene was investigated using primers within the deletion as mapped by the arrays. These primers amplified a fragment in 11/12 and 9/12 individuals in the serpentine1 and serpentine2 populations, but only 3/12 and 5/12 individuals in the granitic1 and granitic2 populations; 24 individuals were typed twice with complete reproducibility. Deletion frequency was calculated assuming Hardy-Weinburg equilibrium, as this is a dominant marker. Primers are available from the authors upon request.</p>
<title>Results and Discussion</title> <title>Genomic analysis</title> <p>To map differentiated polymorphisms between serpentine and granitic <italic>A. lyrata</italic> localities, DNA from three individuals from each of the four localities (##FIG##0##figure 1##) was fragmented, labelled, and hybridized to an <italic>A. thaliana</italic> tiling array (12 total arrays). On the array, each perfect match (PM) probe that matches the genome is adjacent to a mismatch (MM) probe, which has a mismatched base at the middle base pair. Comparison of normalized Σ(PM<sub>i</sub>) and Σ(MM<sub>i</sub>) intensities over the <italic>i</italic> chips for each probe indicates that the probes matching the draft <italic>A lyrata</italic> genome are the most sensitive markers of DNA differentiation, as expected (##FIG##1##figure 2##). For these 371,642 probes, the Σ(PM<sub>i</sub>)−Σ(MM<sub>i</sub>) is positive for 98.8% of probes. For the other 2,481,727 probes on the array, the Σ(PM<sub>i</sub>)−Σ(MM<sub>i</sub>) is positive at 65% of probes. This indicates that probes matching the draft genome are the most sensitive markers for detecting DNA differentiation, but also that the 2.48 million other unique probes contain useful information.</p> <p>To determine which probes overlapped differentiated polymorphisms, we computed a <italic>t</italic>-test <italic>p</italic>-value for between soil types at each array probe; the distribution of <italic>p</italic>-values for probes with a perfect match in the draft <italic>A. lyrata</italic> genome is shown in ##FIG##2##figure 3##. There is a considerable excess of probes with low <italic>p</italic>-values: 2402 probes have <italic>p</italic>&lt;0.001, whereas only 371 are expected by chance. For comparison, <italic>t</italic>-tests were computed between the other two possible combinations of our four collections. Dividing the samples along these additional axes revealed a small excess of probes with low <italic>p</italic>-values in both cases (##FIG##2##figure 3##), but substantially fewer than in the serpentine vs. granitic comparison, indicating that most differentiation between these samples assorts with soil type.</p> <p>After correcting for multiple tests, 39 probes are significantly differentiated between soil types at a bonferroni <italic>p</italic>&lt;0.05 (for this and all further analyses, the full 2.48 million probes were used; see <xref ref-type="sec" rid="s3">Materials and Methods</xref> for details). Using a less conservative false discovery (FDR) rate criterion, 168 probes are significant at a FDR&lt;0.01, and 751 probes are significant at a FDR&lt;0.05. The 39 bonferroni significant probes overlap 33 genes, including two genes with multiple significant probes. When the FDR&lt;0.01 probes are included, 112 genes overlap at least one significant probe; adding the FDR&lt;0.05 set of probes increases this number to 455 genes including 81 genes with 2–12 significant probes each (##SUPPL##0##Data Set S1##). To locate additional differentiated genes, we used a sliding window analysis to determine if probes with low, but not individually significant FDR are clustered in the genome. Using permuted data to generate a null distribution, we mapped windows of 5, 10, and 20 probes that are significantly enriched for probes with FDR&lt;0.30. Seventy-one small genomic regions were significantly differentiated between serpentine and granitic samples (windows with a permutation-based FDR&lt;0.001 were considered significant, see <xref ref-type="sec" rid="s3">methods</xref>). These regions are mostly a few kb in size (mean = 6.6 kb), with 11 regions larger than 10 kb and 3 larger than 30 kb. These 71 regions overlap 184 total genes, including 52 genes already significant from the analysis of individual probes.</p> <p>The combined list of significantly differentiated loci (probes with FDR&lt;0.05 and windows with FDR&lt;0.001) includes 586 genes, 545 of which are annotated with at least one function in <italic>A. thaliana</italic> (based on gene ontology assignments, hereafter referred to as GO terms). When the probes with 0.01&lt;FDR&lt;0.05 are excluded to generate a more conservative list, the number of GO associated genes drops to 263. As noted in the introduction, additional data is required in order to determine how many of these genes are differentiated due to spatially varying selection. We can, however, use the annotations of differentiated loci to prioritize specific genes and molecular functions for further investigation. We find that the genes which are differentiated between soil conditions are associated with a non-random subset of biological functions (##SUPPL##0##Data Set S1##). Membrane proteins (<italic>p</italic>&lt;0.001) and transporters (<italic>p</italic>&lt;0.001) are very overrepresented among differentiated genes. GO terms which are obvious candidates for adaptation to soil salt and nutrient content include potassium ion transport (<italic>p</italic> = 0.006), cellular calcium ion homeostasis (<italic>p</italic> = 0.008), and cation transport (<italic>p</italic> = 0.03). These genes provide a starting point for understanding the molecular basis of adaptation to known environmental differences, such as low Ca∶Mg ratios, and also may indicate additional unappreciated environmental variables. An example of this second category is the differentiated gene <italic>AT4G10380</italic>, which responds to boron limitation in <italic>A. thaliana</italic>
##REF##16679457##[27]##. Differentiation of this gene motivates the investigation of the role of this important micronutrient in soil adaptation.</p> <p>We also discovered suites of related genes that would not be obvious <italic>a priori</italic> candidates for soil adaptation: one of the most overrepresented GO terms is microtubule-based movement (<italic>p</italic>&lt;0.001). Six genes with this function are differentiated, including a putative kinesin heavy chain and four genes which share a kinesin motor protein domain. It is possible that these loci interact with spatially varying environmental variables either directly (gene-by-environment interactions), or to compensate for deleterious side effects of adaptive change at other loci (gene-by-gene-by-environment interactions). An alternative, neutralist, explanation for the overrepresentation of GO terms among differentiated loci is that genes differ in polymorphism rates, and genes with more polymorphism are more likely to have some polymorphisms which are differentiated due to stochastic forces. It should be noted that the overrepresentation of GO terms is not specific to the granitic vs. serpentine population comparison. We also conducted GO tests between the other two combinations of our four populations (serp1+granite1 vs. serp2+granite2 and serp1+granite2 vs. serp2+granite1). In order to have the same power for all tests, the GO test of serpentine vs. granitic soils was re-computed using only the 751 significant probes, and compared to the probes with the lowest 751 <italic>t</italic>-test <italic>p</italic>-vales from the other two combinations of the data. For the significant probes, 33 GO terms are significant, while in both other combinations of the data, many different GO terms are significant (20 and 46 terms, respectively). Whether this is due to spatially varying selection on different process along different geographic axes, or simply variation in neutral differentiation rates remains to be determined.</p> <title>Genetic analysis</title> <p>An induced null mutation in the calcium-proton antiporter <italic>cax1</italic> in <italic>A. thaliana</italic> is the only characterized mutation which enhances survival on soils with low Ca∶Mg ratios ##REF##15948832##[18]##. A related gene, <italic>cax7</italic> (<italic>AT5G17860</italic>), contains a significantly differentiated probe in our analysis. This gene, though functionally distinct from <italic>cax1</italic>
##REF##17086450##[28]##, is therefore an excellent candidate for adaptation to low Ca∶Mg ratios in <italic>A. lyrata</italic>. We sequenced 1.3 Kb of this gene in a larger sample of plants from the four <italic>A. lyrata</italic> collections and found a small region of high differentiation towards the 5′ end of the coding sequence (##FIG##3##figure 4##, ##TAB##0##table 1##). Patterns of polymorphisms at this locus are quite different on the different soils: individuals from granitic soils have high nucleotide diversity, many intermediate frequency amino acid polymorphisms, and significantly positive values of Tajima's D, Fu and Li's D, and Fu and Li's F (<italic>p</italic>&lt;0.05; ##REF##2513255##[29]##, ##REF##8454210##[30]##. Plants from serpentine soils, on the contrary, have low levels of polymorphism and negative (though non-significant) values of the above statistics; Fay and Wu's H is −3.31 for serpentine plants, indicating many high frequency derived mutations on these soils. One 70 bp region of <italic>cax7</italic> exemplifies these patterns: in this region, 8 of 23 amino acids are polymorphic in granitic collections, perhaps due to ongoing gene flow from serpentine plants, while all are fixed on serpentine, consistent with strong purifying selection against such gene flow. Fine-scale studies of this gene involving many individuals and environmental conditions can now be conducted to determine the environmental and population-genetic forces influencing this variation, and the phenotypic consequences of the polymorphic alleles.</p> <p>We also sequenced 2.2 kb of <italic>AT4G35290</italic>, one of three putative ligand-gated ion channels with a significant probe, contributing to the significant overrepresentation of genes involved in cellular calcium ion homeostasis. Unlike <italic>cax7</italic>, we found similar levels of polymorphism at this locus in the four populations, with no skews in the allele frequency spectrum (θ<sub>π</sub> = 0.0006 and 0.0012 for granitic populations, 0.0007 and 0.0008 for serpentine populations). Although this gene contained a significant (FDR = 0.04) probe, this probe did not overlap a differentiated polymorphism and is therefore a false discovery. This gene was in our least significant category, 0.05&gt;FDR&gt;0.01, which is expected to include some false discoveries. Including these probes also increases the sensitivity of the analysis, however: the significant probe in <italic>cax7</italic> had a FDR = 0.011, and would have also been excluded had we only used a threshold FDR of 0.01. Despite reduced power when probes with 0.01&lt;FDR&lt;0.05 are excluded, 20 GO terms remain significant, including cellular calcium ion homeostasis. The conservative and comprehensive lists of significant probes, regions, genes and GO terms are available as supplementary data (##SUPPL##0##Data Set S1##).</p> <title>Copy number variation</title> <p>Many of our significant regions appear to be large differentiated duplications and deletions (copy number variants, or CNVs). By comparing the intensity of probes in each significant region to an expected distribution generated using probes from all significant regions, we established that 15 of the 71 differentiated regions are likely CNVs (at a threshold <italic>p</italic>&lt;0.001, see <xref ref-type="sec" rid="s3">methods</xref>); this includes 7/11 regions larger than 10 kb, and all three regions larger than 30 kb. As an example, a 33.9 kb significantly differentiated region on chromosome 4 spans 11 genes (map coordinates refer to the <italic>A. thaliana</italic> assembly). The average intensity of the 99 probes in this region in granitic plants is 1.2 (median normalized, where the median intensity = 1.0), while the serpentine plants have an average intensity of 2.6, which appears to indicate that there is one copy of this region in most granitic individuals and two or three copies in most serpentine individuals. Other cases appear to be differentiated deletions: there are 21 probes which overlap the gene <italic>AT5G27100</italic>, one of three significantly differentiated ligand-gated ion channels. The average intensity of these probes is 0.986 on serpentine, near the median of all probes, and 0.167 on granite. PCR amplification of this gene from an expanded sample of 48 individuals validated the presence of a deletion with higher frequency in granitic populations: deletion frequencies are 0.79 and 0.65 for the granitic collections and 0.24 and 0.38 for the serpentine collections. The existence of CNVs whose frequencies are associated with environmental conditions, including large duplications of many genes, provides an excellent opportunity to study the influence of selection in the initial stages of gene duplication ##UREF##9##[31]##. A similar study of clinal adaptation in <italic>Drosophila melanogaster</italic> also found a considerable contribution of copy number variation to population differentiation, including large variants that affected multiple genes ##REF##18493034##[17]##. These two experiments, together with adaptive CNVs in humans ##REF##17828263##[32]##, indicate that copy number variation is a prolific source of adaptive polymorphism in natural populations. Interestingly, changes in gene copy number have also recently been shown to provide heavy metal tolerance in the closely related species <italic>A. halleri</italic>
##REF##18425111##[33]##.</p> <title>Conclusions</title> <p>Between the four populations used in this study, many more loci are differentiated between serpentine and granitic populations than along other geographic axes. The neutralist interpretation of this result is that these populations are more distantly related due to recent or historical migration patterns. However, recent genomic studies of polymorphism, differentiation, and divergence clearly indicate the pervasive effects of selection ##REF##18246066##[25]##, ##REF##18282109##[26]##, ##REF##17988176##[34]##–##REF##15718463##[36]##. It is therefore very difficult to infer whether the observed patterns of differentiation are due to patterns of migration, rampant selection, or both effects ##REF##18302709##[24]##. We believe that the most powerful approach will be combining “from the genes up” approaches like the current work with more traditional approaches “from the traits down”. For example, reciprocal transplants of individuals between populations can directly measure fitness, which can then be associated with the differentiated polymorphisms detected here.</p>
<title>Results and Discussion</title> <title>Genomic analysis</title> <p>To map differentiated polymorphisms between serpentine and granitic <italic>A. lyrata</italic> localities, DNA from three individuals from each of the four localities (##FIG##0##figure 1##) was fragmented, labelled, and hybridized to an <italic>A. thaliana</italic> tiling array (12 total arrays). On the array, each perfect match (PM) probe that matches the genome is adjacent to a mismatch (MM) probe, which has a mismatched base at the middle base pair. Comparison of normalized Σ(PM<sub>i</sub>) and Σ(MM<sub>i</sub>) intensities over the <italic>i</italic> chips for each probe indicates that the probes matching the draft <italic>A lyrata</italic> genome are the most sensitive markers of DNA differentiation, as expected (##FIG##1##figure 2##). For these 371,642 probes, the Σ(PM<sub>i</sub>)−Σ(MM<sub>i</sub>) is positive for 98.8% of probes. For the other 2,481,727 probes on the array, the Σ(PM<sub>i</sub>)−Σ(MM<sub>i</sub>) is positive at 65% of probes. This indicates that probes matching the draft genome are the most sensitive markers for detecting DNA differentiation, but also that the 2.48 million other unique probes contain useful information.</p> <p>To determine which probes overlapped differentiated polymorphisms, we computed a <italic>t</italic>-test <italic>p</italic>-value for between soil types at each array probe; the distribution of <italic>p</italic>-values for probes with a perfect match in the draft <italic>A. lyrata</italic> genome is shown in ##FIG##2##figure 3##. There is a considerable excess of probes with low <italic>p</italic>-values: 2402 probes have <italic>p</italic>&lt;0.001, whereas only 371 are expected by chance. For comparison, <italic>t</italic>-tests were computed between the other two possible combinations of our four collections. Dividing the samples along these additional axes revealed a small excess of probes with low <italic>p</italic>-values in both cases (##FIG##2##figure 3##), but substantially fewer than in the serpentine vs. granitic comparison, indicating that most differentiation between these samples assorts with soil type.</p> <p>After correcting for multiple tests, 39 probes are significantly differentiated between soil types at a bonferroni <italic>p</italic>&lt;0.05 (for this and all further analyses, the full 2.48 million probes were used; see <xref ref-type="sec" rid="s3">Materials and Methods</xref> for details). Using a less conservative false discovery (FDR) rate criterion, 168 probes are significant at a FDR&lt;0.01, and 751 probes are significant at a FDR&lt;0.05. The 39 bonferroni significant probes overlap 33 genes, including two genes with multiple significant probes. When the FDR&lt;0.01 probes are included, 112 genes overlap at least one significant probe; adding the FDR&lt;0.05 set of probes increases this number to 455 genes including 81 genes with 2–12 significant probes each (##SUPPL##0##Data Set S1##). To locate additional differentiated genes, we used a sliding window analysis to determine if probes with low, but not individually significant FDR are clustered in the genome. Using permuted data to generate a null distribution, we mapped windows of 5, 10, and 20 probes that are significantly enriched for probes with FDR&lt;0.30. Seventy-one small genomic regions were significantly differentiated between serpentine and granitic samples (windows with a permutation-based FDR&lt;0.001 were considered significant, see <xref ref-type="sec" rid="s3">methods</xref>). These regions are mostly a few kb in size (mean = 6.6 kb), with 11 regions larger than 10 kb and 3 larger than 30 kb. These 71 regions overlap 184 total genes, including 52 genes already significant from the analysis of individual probes.</p> <p>The combined list of significantly differentiated loci (probes with FDR&lt;0.05 and windows with FDR&lt;0.001) includes 586 genes, 545 of which are annotated with at least one function in <italic>A. thaliana</italic> (based on gene ontology assignments, hereafter referred to as GO terms). When the probes with 0.01&lt;FDR&lt;0.05 are excluded to generate a more conservative list, the number of GO associated genes drops to 263. As noted in the introduction, additional data is required in order to determine how many of these genes are differentiated due to spatially varying selection. We can, however, use the annotations of differentiated loci to prioritize specific genes and molecular functions for further investigation. We find that the genes which are differentiated between soil conditions are associated with a non-random subset of biological functions (##SUPPL##0##Data Set S1##). Membrane proteins (<italic>p</italic>&lt;0.001) and transporters (<italic>p</italic>&lt;0.001) are very overrepresented among differentiated genes. GO terms which are obvious candidates for adaptation to soil salt and nutrient content include potassium ion transport (<italic>p</italic> = 0.006), cellular calcium ion homeostasis (<italic>p</italic> = 0.008), and cation transport (<italic>p</italic> = 0.03). These genes provide a starting point for understanding the molecular basis of adaptation to known environmental differences, such as low Ca∶Mg ratios, and also may indicate additional unappreciated environmental variables. An example of this second category is the differentiated gene <italic>AT4G10380</italic>, which responds to boron limitation in <italic>A. thaliana</italic>
##REF##16679457##[27]##. Differentiation of this gene motivates the investigation of the role of this important micronutrient in soil adaptation.</p> <p>We also discovered suites of related genes that would not be obvious <italic>a priori</italic> candidates for soil adaptation: one of the most overrepresented GO terms is microtubule-based movement (<italic>p</italic>&lt;0.001). Six genes with this function are differentiated, including a putative kinesin heavy chain and four genes which share a kinesin motor protein domain. It is possible that these loci interact with spatially varying environmental variables either directly (gene-by-environment interactions), or to compensate for deleterious side effects of adaptive change at other loci (gene-by-gene-by-environment interactions). An alternative, neutralist, explanation for the overrepresentation of GO terms among differentiated loci is that genes differ in polymorphism rates, and genes with more polymorphism are more likely to have some polymorphisms which are differentiated due to stochastic forces. It should be noted that the overrepresentation of GO terms is not specific to the granitic vs. serpentine population comparison. We also conducted GO tests between the other two combinations of our four populations (serp1+granite1 vs. serp2+granite2 and serp1+granite2 vs. serp2+granite1). In order to have the same power for all tests, the GO test of serpentine vs. granitic soils was re-computed using only the 751 significant probes, and compared to the probes with the lowest 751 <italic>t</italic>-test <italic>p</italic>-vales from the other two combinations of the data. For the significant probes, 33 GO terms are significant, while in both other combinations of the data, many different GO terms are significant (20 and 46 terms, respectively). Whether this is due to spatially varying selection on different process along different geographic axes, or simply variation in neutral differentiation rates remains to be determined.</p> <title>Genetic analysis</title> <p>An induced null mutation in the calcium-proton antiporter <italic>cax1</italic> in <italic>A. thaliana</italic> is the only characterized mutation which enhances survival on soils with low Ca∶Mg ratios ##REF##15948832##[18]##. A related gene, <italic>cax7</italic> (<italic>AT5G17860</italic>), contains a significantly differentiated probe in our analysis. This gene, though functionally distinct from <italic>cax1</italic>
##REF##17086450##[28]##, is therefore an excellent candidate for adaptation to low Ca∶Mg ratios in <italic>A. lyrata</italic>. We sequenced 1.3 Kb of this gene in a larger sample of plants from the four <italic>A. lyrata</italic> collections and found a small region of high differentiation towards the 5′ end of the coding sequence (##FIG##3##figure 4##, ##TAB##0##table 1##). Patterns of polymorphisms at this locus are quite different on the different soils: individuals from granitic soils have high nucleotide diversity, many intermediate frequency amino acid polymorphisms, and significantly positive values of Tajima's D, Fu and Li's D, and Fu and Li's F (<italic>p</italic>&lt;0.05; ##REF##2513255##[29]##, ##REF##8454210##[30]##. Plants from serpentine soils, on the contrary, have low levels of polymorphism and negative (though non-significant) values of the above statistics; Fay and Wu's H is −3.31 for serpentine plants, indicating many high frequency derived mutations on these soils. One 70 bp region of <italic>cax7</italic> exemplifies these patterns: in this region, 8 of 23 amino acids are polymorphic in granitic collections, perhaps due to ongoing gene flow from serpentine plants, while all are fixed on serpentine, consistent with strong purifying selection against such gene flow. Fine-scale studies of this gene involving many individuals and environmental conditions can now be conducted to determine the environmental and population-genetic forces influencing this variation, and the phenotypic consequences of the polymorphic alleles.</p> <p>We also sequenced 2.2 kb of <italic>AT4G35290</italic>, one of three putative ligand-gated ion channels with a significant probe, contributing to the significant overrepresentation of genes involved in cellular calcium ion homeostasis. Unlike <italic>cax7</italic>, we found similar levels of polymorphism at this locus in the four populations, with no skews in the allele frequency spectrum (θ<sub>π</sub> = 0.0006 and 0.0012 for granitic populations, 0.0007 and 0.0008 for serpentine populations). Although this gene contained a significant (FDR = 0.04) probe, this probe did not overlap a differentiated polymorphism and is therefore a false discovery. This gene was in our least significant category, 0.05&gt;FDR&gt;0.01, which is expected to include some false discoveries. Including these probes also increases the sensitivity of the analysis, however: the significant probe in <italic>cax7</italic> had a FDR = 0.011, and would have also been excluded had we only used a threshold FDR of 0.01. Despite reduced power when probes with 0.01&lt;FDR&lt;0.05 are excluded, 20 GO terms remain significant, including cellular calcium ion homeostasis. The conservative and comprehensive lists of significant probes, regions, genes and GO terms are available as supplementary data (##SUPPL##0##Data Set S1##).</p> <title>Copy number variation</title> <p>Many of our significant regions appear to be large differentiated duplications and deletions (copy number variants, or CNVs). By comparing the intensity of probes in each significant region to an expected distribution generated using probes from all significant regions, we established that 15 of the 71 differentiated regions are likely CNVs (at a threshold <italic>p</italic>&lt;0.001, see <xref ref-type="sec" rid="s3">methods</xref>); this includes 7/11 regions larger than 10 kb, and all three regions larger than 30 kb. As an example, a 33.9 kb significantly differentiated region on chromosome 4 spans 11 genes (map coordinates refer to the <italic>A. thaliana</italic> assembly). The average intensity of the 99 probes in this region in granitic plants is 1.2 (median normalized, where the median intensity = 1.0), while the serpentine plants have an average intensity of 2.6, which appears to indicate that there is one copy of this region in most granitic individuals and two or three copies in most serpentine individuals. Other cases appear to be differentiated deletions: there are 21 probes which overlap the gene <italic>AT5G27100</italic>, one of three significantly differentiated ligand-gated ion channels. The average intensity of these probes is 0.986 on serpentine, near the median of all probes, and 0.167 on granite. PCR amplification of this gene from an expanded sample of 48 individuals validated the presence of a deletion with higher frequency in granitic populations: deletion frequencies are 0.79 and 0.65 for the granitic collections and 0.24 and 0.38 for the serpentine collections. The existence of CNVs whose frequencies are associated with environmental conditions, including large duplications of many genes, provides an excellent opportunity to study the influence of selection in the initial stages of gene duplication ##UREF##9##[31]##. A similar study of clinal adaptation in <italic>Drosophila melanogaster</italic> also found a considerable contribution of copy number variation to population differentiation, including large variants that affected multiple genes ##REF##18493034##[17]##. These two experiments, together with adaptive CNVs in humans ##REF##17828263##[32]##, indicate that copy number variation is a prolific source of adaptive polymorphism in natural populations. Interestingly, changes in gene copy number have also recently been shown to provide heavy metal tolerance in the closely related species <italic>A. halleri</italic>
##REF##18425111##[33]##.</p> <title>Conclusions</title> <p>Between the four populations used in this study, many more loci are differentiated between serpentine and granitic populations than along other geographic axes. The neutralist interpretation of this result is that these populations are more distantly related due to recent or historical migration patterns. However, recent genomic studies of polymorphism, differentiation, and divergence clearly indicate the pervasive effects of selection ##REF##18246066##[25]##, ##REF##18282109##[26]##, ##REF##17988176##[34]##–##REF##15718463##[36]##. It is therefore very difficult to infer whether the observed patterns of differentiation are due to patterns of migration, rampant selection, or both effects ##REF##18302709##[24]##. We believe that the most powerful approach will be combining “from the genes up” approaches like the current work with more traditional approaches “from the traits down”. For example, reciprocal transplants of individuals between populations can directly measure fitness, which can then be associated with the differentiated polymorphisms detected here.</p>
<title>Conclusions</title> <p>Between the four populations used in this study, many more loci are differentiated between serpentine and granitic populations than along other geographic axes. The neutralist interpretation of this result is that these populations are more distantly related due to recent or historical migration patterns. However, recent genomic studies of polymorphism, differentiation, and divergence clearly indicate the pervasive effects of selection ##REF##18246066##[25]##, ##REF##18282109##[26]##, ##REF##17988176##[34]##–##REF##15718463##[36]##. It is therefore very difficult to infer whether the observed patterns of differentiation are due to patterns of migration, rampant selection, or both effects ##REF##18302709##[24]##. We believe that the most powerful approach will be combining “from the genes up” approaches like the current work with more traditional approaches “from the traits down”. For example, reciprocal transplants of individuals between populations can directly measure fitness, which can then be associated with the differentiated polymorphisms detected here.</p>
<p>Conceived and designed the experiments: TLT EJvW SVN. Performed the experiments: TLT EJvW SVN. Analyzed the data: TLT EJvW SVN. Wrote the paper: TLT EJvW SVN.</p> <p>Serpentine soil, which is naturally high in heavy metal content and has low calcium to magnesium ratios, comprises a difficult environment for most plants. An impressive number of species are endemic to serpentine, and a wide range of non-endemic plant taxa have been shown to be locally adapted to these soils. Locating genomic polymorphisms which are differentiated between serpentine and non-serpentine populations would provide candidate loci for serpentine adaptation. We have used the <italic>Arabidopsis thaliana</italic> tiling array, which has 2.85 million probes throughout the genome, to measure genetic differentiation between populations of <italic>Arabidopsis lyrata</italic> growing on granitic soils and those growing on serpentinic soils. The significant overrepresentation of genes involved in ion transport and other functions provides a starting point for investigating the molecular basis of adaptation to soil ion content, water retention, and other ecologically and economically important variables. One gene in particular, <italic>calcium-exchanger 7</italic>, appears to be an excellent candidate gene for adaptation to low Ca∶Mg ratio in <italic>A. lyrata</italic>.</p>
<title>Supporting Information</title>
<p>We are grateful to Roger Latham for assistance with locating populations, and Maureen Stanton, Dave Remington, Lindy McBride and Michael Turelli for helpful discussion.</p>
<fig id="pone-0003183-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003183.g001</object-id><label>Figure 1</label><caption><title>Map of collection locations.</title></caption></fig> <fig id="pone-0003183-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003183.g002</object-id><label>Figure 2</label><caption><title>Distributions of hybridization intensities.</title><p>Probes with a perfect match in the draft <italic>A. lyrata</italic> genome have the highest mean intensity (red), and the corresponding mismatches for these probes have the next highest mean (orange). For probes without a match in the draft genome, intensities are still higher at perfect match (grey) than mismatch (purple) probes. Though some of these probes might have a matching DNA sequence in our populations, but not in the reference genome, many of these probes are probably comparisons between a single mismatch and a double mismatch, yet are still providing valuable information. Note that extreme values are not shown for both axes.</p></caption></fig> <fig id="pone-0003183-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003183.g003</object-id><label>Figure 3</label><caption><title>Histogram of <italic>t</italic>-test <italic>p</italic>-values between collection localities (bin size = 0.001).</title><p>Three combinations of our four populations were tested: test 1 is the soil type comparison (blue line) = serp1+serp2 vs. granite1+granite2. For comparison, we also computed the other two combinations of the four localities: test 2 (red line) = serp1+granite1 vs. serp2+granite2; test 3 (green line) = serp2+granite1 vs. serp1+granite2).</p></caption></fig> <fig id="pone-0003183-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003183.g004</object-id><label>Figure 4</label><caption><title>Sequencing of cax7 revealed a small region of high differentiation between soil types: 200 bp sliding windows of FST are shown in black.</title><p>Serpentine populations have much less polymorphism in this region that granitic populations (200 bp sliding windows of heterozygosity are shown in blue for serpentine, red for granite). Vertical lines delineate a region of 70 bp where 8 of 23 amino acids are polymorphic in granitic populations and fixed in serpentine populations.</p></caption></fig>
<table-wrap id="pone-0003183-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003183.t001</object-id><label>Table 1</label><caption><title>Analysis of sequence data, with individuals grouped by soil type.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">bp</td><td align="left" rowspan="1" colspan="1">n</td><td align="left" rowspan="1" colspan="1">F<sub>ST</sub>
</td><td align="left" rowspan="1" colspan="1">θπ</td><td align="left" rowspan="1" colspan="1">Taj D</td><td align="left" rowspan="1" colspan="1">Fu&amp;Li D</td><td align="left" rowspan="1" colspan="1">Fu&amp;Li F</td><td align="left" rowspan="1" colspan="1">Fay&amp;Wu H</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">AT5G17860 (cax7)</td><td align="left" rowspan="1" colspan="1">granite</td><td align="left" rowspan="1" colspan="1">1326</td><td align="left" rowspan="1" colspan="1">24</td><td align="left" rowspan="1" colspan="1">0.45</td><td align="left" rowspan="1" colspan="1">0.0144</td><td align="left" rowspan="1" colspan="1">2.239</td><td align="left" rowspan="1" colspan="1">2.045</td><td align="left" rowspan="1" colspan="1">2.611</td><td align="left" rowspan="1" colspan="1">4.587</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">serpentine</td><td align="left" rowspan="1" colspan="1">1326</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.0036</td><td align="left" rowspan="1" colspan="1">−0.966</td><td align="left" rowspan="1" colspan="1">−0.670</td><td align="left" rowspan="1" colspan="1">−0.970</td><td align="left" rowspan="1" colspan="1">−3.305</td></tr><tr><td align="left" rowspan="1" colspan="1">AT4G35290</td><td align="left" rowspan="1" colspan="1">granite</td><td align="left" rowspan="1" colspan="1">2167</td><td align="left" rowspan="1" colspan="1">22</td><td align="left" rowspan="1" colspan="1">0.11</td><td align="left" rowspan="1" colspan="1">0.0002</td><td align="left" rowspan="1" colspan="1">−0.312</td><td align="left" rowspan="1" colspan="1">0.413</td><td align="left" rowspan="1" colspan="1">0.227</td><td align="left" rowspan="1" colspan="1">1.750</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">serpentine</td><td align="left" rowspan="1" colspan="1">2167</td><td align="left" rowspan="1" colspan="1">22</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.0002</td><td align="left" rowspan="1" colspan="1">−1.720</td><td align="left" rowspan="1" colspan="1">−1.547</td><td align="left" rowspan="1" colspan="1">−1.862</td><td align="left" rowspan="1" colspan="1">1.403</td></tr></tbody></table></alternatives></table-wrap>
<supplementary-material content-type="local-data" id="pone.0003183.s001"><label>Data Set S1</label><caption><p>Complete lists of differentiated probes, regions, genes, and GO terms at multiple significance threshold. (Note that this is a large data set formatted as an .xls file with multiple spreadsheets. We have converted the initial sheet, of GO terms into a .pdf, but cannot easily submit the entire data set as .pdf.)</p><p>(0.03 MB PDF)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><label/><p>bp = base pairs sequenced, n = chromosomes sequenced, θπ = heterozygosity.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>This work was supported by NSF Predoctoral Fellowship and the Center for Population Biology at UC Davis (TLT), an NIH NRSA Postdoctoral Fellowship (EVW), NIH grant RGM076643 (SVN), and NSF grant DEB0548991 (SVN). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p></fn></fn-group>
<graphic xlink:href="pone.0003183.g001"/> <graphic xlink:href="pone.0003183.g002"/> <graphic xlink:href="pone.0003183.g003"/> <graphic xlink:href="pone.0003183.g004"/> <graphic id="pone-0003183-t001-1" xlink:href="pone.0003183.t001"/>
<media xlink:href="pone.0003183.s001.pdf"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["1"], "element-citation": ["\n"], "surname": ["Kruckeberg"], "given-names": ["AR"], "year": ["1951"], "article-title": ["Intraspecific variability in the response of certain native plant species to serpentine soil."], "source": ["American Journal of Botany"], "volume": ["38(6)"], "fpage": ["408"], "lpage": ["419"]}, {"label": ["2"], "element-citation": ["\n"], "surname": ["Kruckeberg"], "given-names": ["AR"], "year": ["1984"], "source": ["California Serpentines: Flora, Vegetation, Geology, Soils, and Management Problems"], "publisher-loc": ["Berkeley, CA"], "publisher-name": ["University of California Press"]}, {"label": ["3"], "element-citation": ["\n"], "surname": ["Brady", "Kruckeberg", "Bradshaw"], "given-names": ["KU", "AR", "HD"], "year": ["2005"], "article-title": ["Evolutionary ecology of plant adaptation to serpentine soils."], "source": ["Annual Review of Ecology Evolution and Systematics"], "volume": ["36"], "fpage": ["243"], "lpage": ["266"]}, {"label": ["4"], "element-citation": ["\n"], "surname": ["Kruckeberg"], "given-names": ["AR"], "year": ["2002"], "source": ["Geology and Plant Life: The Effects of Landforms and Rock Types on Plants"], "publisher-loc": ["Seattle, WA"], "publisher-name": ["University of Washington Press"]}, {"label": ["5"], "element-citation": ["\n"], "surname": ["Safford", "Viers", "Harrison"], "given-names": ["HD", "JH", "SP"], "year": ["2005"], "article-title": ["Serpentine endemism in the California flora: a database of serpentine affinity."], "source": ["Madrono"], "volume": ["52"], "fpage": ["222"], "lpage": ["257"]}, {"label": ["6"], "element-citation": ["\n"], "surname": ["Wright", "Stanton", "Scherson"], "given-names": ["JW", "ML", "R"], "year": ["2006"], "article-title": ["Local adaptation to serpentine and non-serpentine soils in "], "italic": ["Collinsia sparsiflora"], "source": ["Evolutionary Ecology Research"], "volume": ["8(1)"], "fpage": ["1"], "lpage": ["21"]}, {"label": ["7"], "element-citation": ["\n"], "surname": ["Berglund", "Dahlgren", "Westerbergh"], "given-names": ["ABN", "S", "A"], "year": ["2004"], "article-title": ["Evidence for parallel evolution and site-specific selection of serpentine tolerance in "], "italic": ["Cerastium alpinum"], "source": ["New Phytologist"], "volume": ["161(1)"], "fpage": ["199"], "lpage": ["209"]}, {"label": ["8"], "element-citation": ["\n"], "surname": ["Wright"], "given-names": ["JW"], "year": ["2007"], "article-title": ["Local adaptation to serpentine soils in "], "italic": ["Pinus ponderosa"], "source": ["Plant and Soil"], "volume": ["293(1\u20132)"], "fpage": ["209"], "lpage": ["217"]}, {"label": ["9"], "element-citation": ["\n"], "surname": ["Rhoads", "Klein"], "given-names": ["AF", "WH"], "year": ["1993"], "source": ["The Vascular Flora of Pennsylvania: Annotated Checklist and Atlas"], "publisher-loc": ["Philadelphia, PA"], "publisher-name": ["American Philosophical Society"]}, {"label": ["31"], "element-citation": ["\n"], "surname": ["Lynch"], "given-names": ["M"], "year": ["2007"], "source": ["The Origins of Genome Architecture"], "publisher-loc": ["Sunderland, MA"], "publisher-name": ["Sinauer Associates, Inc"]}, {"label": ["37"], "element-citation": ["\n"], "surname": ["Latham"], "given-names": ["RE"], "year": ["1993"], "article-title": ["The serpentine barrens of temperate eastern North America: critical issues in the management of rare species and communities."], "source": ["Bartonia"], "volume": ["57"], "issue": ["Suppl."], "fpage": ["61"], "lpage": ["74"]}]
PLoS One. 2008 Sep 11; 3(9):e3183