Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Research article

Contrasting patterns of genetic divergence in two sympatric pseudo-metallophytes: Rumex acetosa L. and Commelina communis L.

M Ye13, B Liao1, JT Li1, A Mengoni2, M Hu1, WC Luo1 and WS Shu1*

Author Affiliations

1 School of Life Sciences and State Key Laboratory of Biocontrol, Sun Yat-sen University, Guangzhou, 510275, People's Republic of China

2 Department of Evolutionary Biology, University of Firenze, via Romana 17, I-50125, Florence, Italy

3 Guangdong Provincial Academy of Environmental Science, Guangzhou, 510045, People's Republic of China

For all author emails, please log on.

BMC Evolutionary Biology 2012, 12:84  doi:10.1186/1471-2148-12-84

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/12/84


Received:5 December 2011
Accepted:18 May 2012
Published:13 June 2012

© 2012 Ye et al.; licensee BioMed Central Ltd

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Patterns of genetic divergence between populations of facultative metallophytes have been investigated extensively. However, most previous investigations have focused on a single plant species making it unclear if genetic divergence shows common patterns or, conversely, is species-specific. The herbs Rumex acetosa L. and Commelina communis L. are two pseudo-metallophytes thriving in both normal and cupriferous soils along the middle and lower reaches of the Yangtze River in China. Their non-metallicolous and metallicolous populations are often sympatric thus providing an ideal opportunity for comparative estimation of genetic structures and divergence under the selective pressure derived from copper toxicity.

Results

In the present study, patterns of genetic divergence of R. acetosa and C. communis , including metal tolerance, genetic structure and genetic relationships between populations, were investigated and compared using hydroponic experiments, AFLP, ISSR and chloroplast genetic markers. Our results show a significant reduction in genetic diversity in metallicolous populations of C. communis but not in R. acetosa . Moreover, genetic differentiation is less in R. acetosa than in C. communis , the latter species also shows a clustering of its metallicolous populations.

Conclusions

We propose that the genetic divergences apparent in R. acetosa and C. communis , and the contrasting responses of the two species to copper contamination, might be attributed to the differences in their intrinsic physiological and ecological properties. No simple and generalised conclusions on genetic divergence in pseudo-metallophytes can thus be drawn.

Background

Understanding the genetic basis of how organisms adapt to extreme and novel environments remains as one of the major challenges facing evolutionary biologists. Intraspecific differentiation of metallophytes is one of the most striking adaptive changes that result from natural outcrops of heavy metal-rich rock or mining activities [1-3]. Metallophytes refer to plant species that have developed metabolic mechanisms to resist, tolerate, or even thrive on toxic metalliferous soils [4-6]. Populations of metallophytes growing in metalliferous soils are often genetically distinct from those living in nearby non-metalliferous soils [7]. Metalliferous habitats are always characterized by their patchy distribution and therefore have been considered as ecological islands, which provide ideal opportunities to investigate the population differentiation in plants under severe edaphic pressure [8]. In the course of evolution, plants have adapted to widely differing metal availabilities in soils and therefore represent an important source of natural variation of metal homeostasis networks [9]. In the past 50 years, genetic differentiation within and among populations of metallophytes has attracted uninterrupted interest. Yet, the evolution of metallophytes remains heatedly debated [10-13]. Previous studies on genetic divergence of metallophytes have mainly focused on the following two aspects: (1) the effects of heavy metal pollution on genetic diversity, and (2) the origin and evolutionary history of metallicolous populations. In general, substantial founder and bottle-neck effects have been presumed to occur during the colonization of metallophytes in a heavy metal contaminated area. Therefore, a reduction in genetic diversity in metallicolous populations would be theoretically expected. This idea has been well-supported by several studies [14,15]. Conflicting results, however, have also been obtained in some cases where recently established tolerant populations maintain a high level of genetic variation comparable to that in non-metallicolous populations [16,17]. On the other hand, Schat et al. [18] have suggested that geographically distant conspecific metallicolous populations could have evolved independently and therefore constitute a polyphyletic group. The hypothesis of multiple and independent origins has been validated in metallicolous populations of several metallophytes such as Noccaea caerulescens (formerly known as Thlaspi caerulescens ), Silene paradoxa Arabidopsis halleri and A. lyrata [17,19-21]. Beyond all doubt, both mosaic environments and edaphic conditions should contribute to the multiple origins of metallicolous populations. However, it is still not clear whether or not edaphic conditions played a more important role than geographical factors in this event. The conflicting results associated with such an issue might be due to: (1) differences in biological traits, ecological properties, population size, and history of individual populations of different species; and/or (2) various molecular markers (i .e . RAPD, RFLP, ISSR, and AFLP) employed in different studies. Moreover, it should also be noted that until now, most previous studies have focused on a single plant species only, which inevitably influences the comparability of results from different studies. Therefore, choosing two or more sympatric metallophytes with similar geographic distributions and population history, and then comparing their characteristics of genetic divergence might shed light on the evolutionary processes of plants under heavy metal pressure.

The copper (Cu) belts distributing along the middle and lower reaches of the Yangtze River are one of the largest areas of Cu mining in China. In this area, there are considerable spoil heaps and mineral deposits, some of which are currently in practice. A national ruin of an ancient Cu mine built over 1000 years ago (Da-Gong Mountain Mining Ruin) is located in this area, suggesting a long history of mining activities. Extensive Cu mining activities in this area have resulted in many waste heaps with elevated concentrations of Cu in soils, and the plant species growing in this area have experienced long-term selection under Cu toxicity derived from mining activities dating back to at least 1000 years. Although these waste heaps are often far away from each other, they harbor similar assemblages of metallophytes, which are distinctly different from the surrounding communities in non-metalliferous soils [22-24]. Both common sorrel Rumex acetosa L. (Polygonaceae) and Asiatic dayflower Commelina communis L. (Commelinaceae ) have consistently been documented as dominant plants on most of these cupriferous habitats and are also present in surrounding non-metalliferous sites [21,22]. Besides, R. acetosa is well-known metallophyte thriving on metalliferous soils of NW & SW Europe. Many edaphic ecotypes of R. acetosa from a wide range of soil types have been reported, some of which exhibit high potential for phytoremediation of contaminated soils [2,25,26].

Rumex acetosa is a perennial herb with an erect stem, about 40 to 120 cm high. Its flowering season ranges from March to May and the fruiting season ranges from April to June. Unisexual and dioecious flowers are arranged in acrogenous panicle inflorescences. Plants are capable of producing large numbers of small achenes encircled with an aliform membrane and therefore can be easily dispersed by wind. The species is commonly found along roadsides, in forests, and in habitats along rice fields in most regions of China [27]. In contrast, C . communis is an annual multi-branched herb with erect stems in the upper part and creeping stems in the lower part. It is able to grow clonally by elongation of the creeping stem, which bears a number of erect ascending shoots. It flowers from May to September and fructifies from June to November. Bisexual flowers are arranged in a cyme. Commelina communis often occupies habitats similar to those of R. acetosa and is mainly distributed in southeastern China [28].

Despite their similarities, the two pseudo-metallophtes R. acetosa and C. communis have different life histories (perennial vs annual) and reproductive systems (obligately outcrossing vs potentially selfing), which result in them being in different functional groups. Therefore, their metallicolous and non-metallicolous populations provide an ideal opportunity to study genetic divergence under selective pressure derived from Cu toxicity between taxa of two distinct functional groups.

Previous studies have shown that many factors such as life history and reproductive system may significantly influence population genetic divergence through their effects on genetic drift and gene flow [8,15-17]. We therefore hypothesized that the two sympatric pseudo-metallophytes R. acetosa and C. communis would differ greatly in population genetic divergence. To test our hypothesis, 12 R. acetosa populations and 13 C. communis populations collected from different sites along the Yangtze River of Eastern China were used to investigate comparative genetic differentiation patterns. AFLP (Amplified Fragment Length Polymorphism) and ISSR (Inter-Simple Sequence Repeats) were employed to assess genetic diversity and population structures. Six R. acetosa and 10 C. communis populations were selected to determine Cu tolerance of the metallicolous and non-metallicolous populations. In addition, two cpDNA regions (psbJ-petA and 3′rps16-5′trnK ) were sequenced and analyzed to infer uniparental patterns of colonization.

Results

Cu concentrations in soils

Sampling locations for 12 R. acetosa and 13 C. communis populations were investigated at areas along the Yangtze River (Figure 1). The soil Cu concentrations recorded at these different sampling sites are shown in Table 1. In general, contaminated sites contained high concentrations of total and DTPA-extractable Cu, about 100 times greater than those in uncontaminated sites. Large differences in Cu concentrations were detected between cupriferous sites; the highest total and extractable Cu concentrations were 8587 and 2194 mg kg−1 at TLS, whilst the lowest values were 2787 and 359 mg kg−1 at FHS, respectively.

thumbnailFigure 1 . The sampling locations of 12Rumex acetosa and 13Commelina communis populations. · NM populations of R. acetosa ; ▪ NM populations of C. communis ; ▴ M populations of R. acetosa ; ▾ M populations of C. communis .(1, 3: AJS; 2: JHS; 4, 5: JSCK; 6: TJ; 7, 8: DGSCK; 9, 10: DGS; 11, 12: FHS 13, 14: TLCK;15, 18: SZS; 16, 17: CZCK; 19, 20: SCCK; 21, 22: CZ; 23: HZCK; 24: TLS; 25: TSK).

Table 1. Locations of populations ofRumex acetosa andCommelina communis , and the total and extractable Cu concentrations in the substrates associated with plants sampled

Cu tolerance

Tolerance indices based upon growth quotients of relative root length under the four Cu treatments for the different populations are shown in Figure 2. Overall, Cu tolerance indices for the 6 R. acetosa populations and 10 C. communis populations under four Cu treatments were significantly different between the metallicolous (M) and non-metallicolous (NM) populations; tolerance indices of M populations were significantly higher than those of NM populations for both species. Interestingly, tolerance indices of NM populations of both species were also similar. Commelina communis M populations exhibited higher tolerances than those of R. acetosa in the 160 μM and 320 μM Cu treatments. Shoot elongation was, with few exceptions, considerably depressed in all NM populations investigated when seedlings were treated with 40 μM Cu or above. Moreover, the treatment with 320 μM Cu which killed all NM plants revealed a higher tolerance in C. communis M populations than R. acetosa M populations (Additional file 1: Figure S1 and S2).

Additional file 1 . This file provides the tolerance indices for the 6 Rumex acetosa populations and 10 Commelina communis populations in the four CuSO4 treatments.

Format: DOC Size: 248KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

thumbnailFigure 2 . Tolerance indices (%) based on relative root length of different populations ofRumex acetosa andCommelina communis under four Cu treatments in hydroculture (0.25, 40, 160, and 320 μmol L-1CuSO4, respectively; M: metallicolous; NM: non-metallicolous). Different letters on the error bars indicate a significant difference at P  < 0.05 (LSD test) under the same treatment.

Genetic diversity of Rumex acetosa and Commelina communis populations

The 6 primer combinations used in the AFLP analysis yielded 212 and 238 bands in total and the percentage of polymorphic loci was 84.4% and 86.5% for R. acetosa and C. communis, respectively. ISSR markers scored were 65 and 77 in total and the percentage of polymorphic loci was 66.2% and 74% for R. acetosa and C. communis, respectively. Nei’s gene diversity and Shannon indices of each population based on AFLP and ISSR data are given in Tables 2 and 3.

Table 2. Genetic diversity statistics presented as Nei′s gene diversity (H) and Shannon index (H′) inRumex acetosa andCommelina communis based on AFLP data (for each population separately and also divided in two types according to edaphic conditions)

Table 3. Genetic diversity statistics presented as Nei′s gene diversity (H) and Shannon index (H′) inRumex acetosa andCommelina communis based on ISSR data (for each population separately and also divided in two types according to edaphic conditions)

For R. acetosa no significant difference in Nei’s gene diversity between M populations and NM populations was detected. Further, there was no significant correlation between Nei’s gene diversity and the DTPA-extractable Cu concentrations in associated soils within each population (Pearson’s correlation coefficient rS = 0.119, P  = 0.712, and rS = −0.045, P  = 0.890, for AFLP data and ISSR data, respectively). In contrast, for C. communis significant (P  < 0.001) differences in genetic diversity were detected between the M and NM populations, the latter having higher values than the former. Furthermore, significantly negative correlations between genetic diversity and DTPA-extractable Cu concentrations were detected (Pearson’s correlation coefficient rS =–0.793, P  = 0.001; rS =–0.652, P  = 0.016 for AFLP and ISSR, respectively).

Genetic structure analysis

The genetic structure of R. acetosa and C. communis populations was investigated by an Analysis of Molecular Variance (AMOVA). Relatively high genetic differentiation among populations of both C. communis , (F ST  = 0.129, P  < 0.0001) and R. acetosa (F ST  = 0.11, P  < 0.0001) was detected from AFLP data. Conversely, for R. acetosa ISSRs were not able to detect significant differentiation among populations (F ST  = 0.003, P  > 0.05), whilst a low differentiation was detected for C. communis (F ST  = 0.029, P  < 0.0001). Mantel’s tests carried out to investigate the existence of any correlations between geographical and genetic distances (Additional file 2: Table S1) did not detect significant correlations for AFLP data in either plant species. In contrast, differentiation at ISSR loci did show significant correlation with geographical distance especially for R. acetosa (r = 0.88, P  < 0.0001).

Additional file 2 . This file provides some supporting information for the article, including the result of the Mantel’s test, the results of Slatkin’s linearized pairwise FST, the results of AMOVA, GenBank accession numbers of two cpDNA non-coding regions, the lists of primers for AFLP and ISSR analysis.

Format: DOC Size: 234KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Genetic relationships among populations were then investigated by computing NJ dendrograms from Slatkin’s linearized pairwise F ST (Additional file 2: Table S2). Results are reported in Figures 3 and 4 for C. communis and R. acetosa, respectively. In general, only dendrograms based on AFLP were supported by significant P-values for pairwise F ST for all pairwise comparisons. For C. communis (Figure 3), in both AFLP and ISSR data, M populations were more differentiated than NM populations, these latter being grouped in a single cluster (except for CZ population). Conversely, the Rumex acetosa populations were clustered into several groups which approximately reflected their geographic locations, supporting the results of the Mantel’s test. For example, the four populations sampled from Jiangsu Province (JSCK, TJ, AJS, and JHS) were clustered in one main branch in both AFLP and ISSR analyses (Figure 4).

thumbnailFigure 3 . NJ dendrogram ofC. communis populations. A, AFLP data; b ,ISSR data. In bold M populations. Scale bar, Slatkin’s linearized pairwise F ST .

thumbnailFigure 4 . NJ dendrogram ofR . acetosa populations. a, AFLP data; b ,ISSR data. In bold M populations. Scale bar, Slatkin’s linearized pairwise F ST .

Hierarchical partition was then computed at three levels (groups, populations and individuals). Groupings were related to both edaphic and geographic differences (Additional file 2: Tables S3 and S4 for C. communis and R. acetosa, respectively). Edaphic differentiation was strongly supported for C. communis (18.51% and 6.64% of total variance for ISSR and AFLP data respectively), whereas lower levels were found for R. acetosa (3.68 and 4.04% for ISSR and AFLP data). However, for C. communis, a geographical structuring of populations was also strongly supported by the ISSR data (11.01% of total variance).

Chloroplast markers analysis

Two cpDNA non-coding regions (psbJ -petA and 3′rps16-5′trnK ) of each population investigated were sequenced (GenBank accession numbers HM041054-HM041103, Additional file 2: Table S5). For C. communis , the aligned concatemers formed linearly combining the sequences of psbJ-petA and 3′rps16-5′trnK regions resulted 1492 bp long. The number of variable sites was 60, 37 of which were parsimony-informative. For R. acetosa, the alignment was 1688 bp long, with 22 variable sites, only 3 of which were parsimony-informative. However, due to the low sequence variability, haplotype network was not informative on R. acetosa , while for C. communis populations some clustering of haplotypes was recognized (Figure 5). In particular for C. communis , four out of five metallicolous populations (TSK, CZ, TLS, SZS, DGS, and FHS) were included in one single group, suggesting a common chloroplastic (i.e. maternal) origin.

thumbnailFigure 5 . Haplotype network ofCommelina communis populations based on combined sequences of psbJ-petA and 3′rps16-5′trnK regions.

Discussion

Divergence of metallophytes under heavy metal selection has been characterized by other researchers, with a focus on only a few species, notably N. caerulescens [29-33], A. halleri [1,16,34-36], Cistus ladanifer [37] and species in the genus Silene [15,18,33,38,39]. To our knowledge, the present study is the first report on a comparison of the genetic divergence of two sympatric metallophytes.

Heavy metal tolerance in higher plants is a well-documented example of micro-evolution and is a basic necessity for survival in metal-contaminated soils [9,40]. Based on the assumption that genotypes with a sufficiently high level of metal tolerance are rare in non-metallicolous populations, the founder effect was hypothesized to be a result of strong selection occurred during the colonization of contaminated areas by metallophytes [8]. A comparison of heavy metal tolerance among populations from metalliferous and non-metalliferous sites would provide an estimation of the levels of stress occurred during the colonization of heavy metal-rich areas. Differences in heavy metal tolerance among populations of metallophytes have been reported widely (see [41-44] for examples). In general, populations or ecotypes growing in contaminated soils exhibit higher levels of heavy metal tolerance than those from uncontaminated sites.

Tolerance index is a measure that can be used to assess the relative degree of tolerance of various plant varieties [45]. In this sense, it is important to test the Cu tolerance of non-tolerant species (or non-metallicolous populations in our case) as a control. On the other hand, in most of the relevant literature [4,45], heavy metal tolerances of non-metallicolous species (populations) have been compared to those of metallicolous ones. We therefore believe that this kind of comparison is suitable for a relative assessment within our two species. Yet, it should be noted that we cannot establish the presence of any baseline (constitutive) tolerance above that of other local species.

In the present study, significant differences in Cu tolerance at the population level were detected. Further, differences in Cu tolerance between the two species under investigation were also detected. Under the similar level of Cu stress, C. communis exhibited higher Cu tolerance than R. acetosa , especially when Cu concentration reached the maximum of 320 μmol L−1, Cu tolerances of the 4 C. communis populations from cupriferous sites were significantly higher than those of the 4 R. acetosa populations from the same sites (P  < 0.05). Interestingly, there was no significant difference in Cu tolerance of non-metallicolous populations between the two species, suggesting a possible similar level of selective pressure on the ancestral non-tolerant colonizers of metalliferous sites. Although the concentrations of Cu in the contaminated sites we sampled were extremely high when compared with those of the non-metalliferous sites, distinct differences in Cu concentration were also observed among contaminated sites. Therefore, heterogeneity in the metalliferous sites could partly explain the difference in Cu tolerance among metallicolous populations. However, considering that the Cu concentrations in contaminated sites were generally 100 times greater than those in uncontaminated sites, the significant differences in Cu tolerance should be mostly attributed to the intense and long-term selection due to heavy metal stresses and to the differences in biological traits of the two metallophytes. Because they were sympatric in most sampling sites, and so subjected to the similar edaphic conditions. On the other hand, it is possible that the heterogeneous nature of the metalliferous sites could also partly explain the difference of Cu tolerance between the two species. Because there were some significant differences in Cu concentration among the metalliferous sites (Table 1). However, it should also be noted that we cannot exclude the potential contribution of maternal environmental effects on Cu tolerance of different populations of the two species, since field collected seeds were used and in the present study measurements were taken on young seedlings.

In order to determine more thoroughly whether the different phenotypes observed are linked to population histories and are reflected in the genetic diversity of populations, three different molecular techniques (AFLP, ISSR, and cpDNA) were employed to analyze the genetic structures and genetic relationships among different populations of the two metallophytes. Firstly, we found that levels of genetic diversity in the two co-occurring species were distinctly different. For R. acetosa , no significant difference in genetic diversity between metallicolous and non-metallicolous populations was found, whereas significantly lower levels of genetic diversity were detected in metallicolous populations than in non-metallicolous populations of C. communis . Similarly, Mengoni et al. [15] reported a significant decrease in genetic diversity of populations in S. paradoxa collected from Cu deposits. However, contrasting results were obtained in studies on other species. For example, the metallicolous and non-metallicolous populations of A. halleri [16], N. caerulescens [31] and Onosma echioides [46] maintained similar levels of genetic diversity. In a more recent study on N. caerulescens , genetic differentiation linked to heavy metal concentrations in soil was detected, and the gene flow observed at some nuclear loci was shown to be significantly reduced between plants encountering different levels of heavy metal contamination in the soil, suggesting that natural selection limits gene flow between metalliferous and non- metalliferous locations [32]. In the present study, the reduction of genetic diversity in metallicolous populations of C. communis might be due to a strong selective pressure during the colonization of Cu-contaminated sites by this species, coupled with a limited gene flow from surrounding non-metallicolous populations (as suggested by the relatively high values of pairwise F ST ). In contrast, the similar levels of genetic diversity in metallicolous and non-metallicolous populations of R. acetosa might be due to substantial gene flow between populations, as indicated by the low values of pairwise F ST . High levels of gene flow between metallicolous and non-metallicolous populations of metallophytes have previously been reported in some other comparable studies [16,17,31]. In addition, non-metallicolous populations of R. acetosa showed higher levels of background tolerance to Cu (especially for DGSCK) compared to those of C. communis. Considering that metallicolous populations are presumed to evolve from non-metallicolous populations, this evidence may imply that selection for Cu tolerance in C. communis has been greater than in R. acetosa , although present-day levels of Cu tolerance in metallicolous populations are lower in R. acetosa .

Concerning that the habitats of metallicolous populations of metallophytes are often fragmented and disjunct, it is unlikely that dispersal from a single tolerant ancestral population could have produced the wide ranges of geographic distribution today. Hence, the hypothesis of polyphyletic origin is proposed and has been subsequently corroborated by considerable data [15,17-19,31]. In the case of R. acetosa , metallicolous populations branched in separate positions in the NJ dendrograms obtained by AFLP and ISSR data, supporting the hypothesis of independent origins from nearby normal populations. Furthermore, we also found that there was a direct correlation between genetic differentiation and geographical distances between pairs of R. acetosa populations (in particular for ISSR markers, Additional file 2: Table S1). For C. communis , we also speculate that the metallicolous populations evolved independently, since the site-to-site distances of the metallicolous populations were too far to be overcome by natural dispersal. However, metallicolous and non-metallicolous populations of C. communis had a distinct branching pattern in the NJ dendrograms with AFLP and ISSR data. The non-metallicolous populations were included in one group, whilst the metallicolous populations dispersed into single long branches. Additionally, components of variance of differentiation between populations of C. communis were also higher than those of R. acetosa . These results suggested a low gene flow in C. communis not only between metallicolous and non-metallicolous populations, but also among metallicolous populations, supporting higher levels (and possibly rates) of genetic differentiation.

The contrasting pattern of genetic differentiation between C. communis and R. acetosa was also confirmed by the analysis of chloroplast markers. Though this analysis was performed on only one individual per population, results for C. communis were similar with those obtained with AFLP, suggesting the presence of one main cluster of metallicolous populations. On this basis, it could be reasonable to hypothesize that metallicolous populations of C. communis may derive from seeds of one single ancestral population. Interestingly, we observed a very low variability of chloroplast sequences of R. acetosa , which seemed to be consistent with the idea of a recent spread and low differentiation between metallicolous and non-metallicolous populations.

The contrasting results of the two metallophytes also implied that populations of C. communis were more remarkably influenced by heavy metal stress than those of R. acetosa , although the genetic structures of the two species were influenced by geographic isolation and heavy metal contamination.

Although ecogeographic isolation has long been viewed as the most important reproductive barrier in plants [47], barriers to gene flow were also demonstrated between non-metallicolous and closely adjacent metallicolous populations of some species such as the grass Anthoxanthum odoratum [48]. Additionally, it has been reported that heavy metal contamination could have a greater impact on the population structure of the hyperaccumulator Sedum alfredii than geographic distance [49]. Based on our results, we can therefore speculate that both heavy metals and geographic distance play a significant role in determining the population structure of R. acetosa . In contrast, Cu contamination seemed to play a more important role in determining the population structure of C. communis than geographic distance.

In previous studies, it has been shown that the genetic divergence and evolutionary processes of plants could be affected by various factors, such as edaphic conditions, stress intensity and duration, geographic isolation, bottleneck and founder effects, life history, reproductive system, and so on [16,29]. Considering that R. acetosa and C. communis populations in the present study almost share the same edaphic conditions (mainly Cu toxicity) and geographic distribution pattern, the different characteristics of genetic divergence might result largely from their different life history and reproductive system. On the one hand, as a perennial, unisexual, and dioecious species, R. acetosa may have more opportunities to experience gene exchange among individuals and gene flow between populations than the bisexual annual plant C. communis . In addition, the winged seeds of R. acetosa are presumed to have higher dispersal capacity than those of C. communis , which might enhance the success of colonization by R. acetosa [27,28].

Conclusions

Contrasting patterns of genetic divergence in two pseudo-metallophytes that experience similar selective pressures provides insight into evolutionary processes under heavy metal stress. Our findings indicate that genetic divergence of metallophytes is the result of the interactions between biological properties of a species (e.g. seed dispersal, background tolerance) and geographical patterns of colonization. Yet, no simplistic conclusions (e.g. the presence of founder effects etc. ) can be drawn concerning the microevolutionary dynamics of all metallophytes.

Methods

Sampling sites

Nine metalliferous sites along the Yangtze River of Eastern China and six nearby non-metalliferous sites were selected as sampling sites (Table 1 and Figure 1). The sampling was conducted during the period from May to June 2007.

Soil sampling and analysis

From 3 to 10 soil samples were collected at each sampling site (0–20 cm depth). Soil samples were air-dried and passed through a 2-mm sieve. For analysis of total Cu, soil samples were digested with 65% HNO3 + concentrated HClO4 (5:1) [50]; for analysis of extractable Cu, soil samples were extracted with a diethylenetetraminepentaacetic acid (DTPA) solution [51]. Cu concentrations in digested or extracted solutions were then determined by atomic absorption spectrometry (Hitachi-Z-5300, Hitachi, Japan).

Hydroponic experiment

Six R. acetosa populations and 10 C. communis populations were selected for hydoponic experiments to investigate Cu tolerance in each population (Figure 2). A bulk collection of mature seeds was made from 5–10 plants sampled in the field, depending on the number of seeds on each plant. Seeds were germinated on sand and then 3-week old seedlings were transferred to a nutrient solution in a glasshouse equipped with supplementary lighting with fluorescent tubes (44 W/m2; 14-h photoperiod; 24–26°C). The composition of the nutrient solution was as follows: Ca(NO3)2·4H2O 2.00 mmol L−1, KH2PO4 0.10 mmol L−1,MgSO4·7H2O 0.50 mmol L−1,KCl 0.10 mmol L−1,K2SO4 0.70 mmol L−1,H3BO3 0.01 mmol L−1,MnSO4·H2O 0.50 × 10−3 mmol L−1,ZnSO4·7H2O 0.50 × 10−3 mmol L−1,CuSO4·5H2O 0.20 × 10−3 mmol L−1,(NH4)6Mo7O2·4H2O 0.01 × 10−3 mmol L−1,Fe-EDTA 0.10 mmol L−1. Hydroponic culture vessels were positioned randomly in the growth chamber and rearranged once a week. The nutrient solutions were renewed every 5 days. Before Cu treatment, roots of all seedlings were blackened with activated charcoal and rinsed in deionized water to remove the excess powder. Seedlings were then treated with a range of Cu concentrations (0.25, 40, 160, and 320 μM), supplied as copper sulphate (CuSO4·5H2O). Here, 0.25 μM CuSO4 was the control. Because Cu is an essential plant micronutrient and our preliminary tests did suggest that 0.25 μM CuSO4 ensured good growth in all populations. Each treatment was repeated 6 times. After 25 days of Cu treatment, 6 seedlings were harvested from each treatment for Cu tolerance measurement. The presence of new roots (uncoated) visible beyond the charcoal-coated roots were recorded. The roots are in direct contact with the nutrient solution (or the soil solution in metalliferous soil under field conditions) and are therefore the first target for any toxic effect. Hence, the response of root growth is generally considered to be a very appropriate measure of the metal tolerance of plants. The determination of tolerance indices followed the methods described by Baker & Walker [45], and were calculated based on relative root growth (lengths) generated in parallel treatment and control units.

Sampling strategy for populations used in AFLP and ISSR analyses

For AFLP and ISSR analyses, 12 and 13 populations of R. acetosa and C. communis were selected, respectively (Figure 1). For each population, at least 10 individuals were collected in the field at irregular intervals of 10–20 m in order to avoid collecting genetically-uniform individuals. Each population sampled covered a similar area of about 0.2 ha. All the plant materials were dried in silica gel before DNA extraction.

DNA extraction

Genomic DNA was extracted using the cetyltrimethylammonium bromide (CTAB) protocol of Lodhi et al. [52] with minor modifications. Both the concentration and purity of the extracted DNA were checked using UV absorbance spectrophotometry. The resulting DNA samples were stored at −70°C and −20°C for long-term and short-term storage, respectively.

AFLP protocol

AFLP analysis was carried out according to the method proposed by Vos et al. [53] with minor modifications. In brief, genomic DNA samples were digested using the enzyme combination Eco RI/Mse I and synthetic adaptors were ligated to the restriction fragments using T4 DNA ligase. A pre-amplification PCR was performed using AFLP primers having a single selective nucleotide (EcoRI  + A/MseI  + C). PCR products of the pre-amplification reaction were diluted 30-fold in sterile water then used as template for selective amplification using a combination of AFLP primers. In total, six selective primer pair combinations were used to generate fingerprints for the two species (Additional file 2: Table S6). Selective PCR products were then mixed with 2 μL of formamide loading dye. Mixtures were heated for 5 min at 94°C and then quickly cooled on ice. Three μL of each sample were loaded onto 6% denaturing polyacrylamide gel and run for 2 h at 50 W. Band patterns were visualized with a silver staining method [54]. The developed gel plate was then air-dried overnight and the images of AFLP profiles were recorded electronically with a scanner.

ISSR protocol

Two DNA samples (one from the metallicolous population and the other from non-metallicolous population) from each species were selected to screen the ISSR polymorphisms by using 93 primers (UBC Set #9 plus all primers from the ISSR Resource Website http://www.biosci.ohio-state.edu/~awolfe/ISSR.html webcite) in duplicate PCR reactions. Totally, 8 and 9 primers which produced highly readable, reproducible and polymorphic band patterns were selected for R. acetosa and C. communis , respectively (Additional file 2: Table S7). PCR amplification and gel electrophoresis analysis were carried out according to Yakimowski & Eckert [55].

Chloroplast DNA sequencing

One individual from each sampled population was randomly selected for sequencing of two cpDNA regions (psbJ-petA and 3′rps16-5′trnK ). Primers and PCR protocols are as described in Shaw et al. [56]. PCR products were purified using High Pure PCR Product Purification Kit (Roche Diagnostics, USA) and were sequenced using an ABI 3730 sequencer. Both forward and reverse strands of psbJ-petA and 3′rps16-5′trnK cpDNA regions were sequenced with a 100% overlap.

Statistical and bioinformatic analysis

Data from hydroponic experiments were analyzed by using the statistical package SPSS 16.0 for Windows (SPSS Inc., USA) and non-parametric tests were used to assess the statistical significance of data.

Band patterns were scored in AFLP and ISSR fingerprint experiments and binary matrices for analysis were constructed for all populations of each species (POPGENE Version 1.32 [57]). Genetic diversity within populations was evaluated as the percentage of polymorphic loci, Nei’s gene diversity [58] and Shannon information index [59], using the software POPGENE Version 1.32 [57]. An independent-samples t -test as available in the SPSS statistical package was used to detect the significant differences in genetic diversity values between populations from cupriferous sites and those from uncontaminated sites. The same software was used to assess the correlation between the values of genetic diversity of each population and the associated Cu concentrations in soils. The partitioning of the genetic variance within and among populations and groups of populations was obtained with an analysis of molecular variance (AMOVA; [60]) conducted by Arlequin 3.5.1 [61]. For each species, two alternative groupings of populations were examined: (i) a geographically based arrangement of populations with sampling sites located in different provinces of China; and (ii) a strictly edaphic arrangement with Cu-contaminated sites vs. uncontaminated sites. The analyses of three hierarchical levels were conducted in order to partition the genetic variance into components attributable to different hierarchical levels, i.e. between groups, among populations, and within populations. The significance levels of variance components were calculated by 1000 permutations in each analysis. Dendrograms of the populations of the two metallophytes were constructed from Slatkin’s linearized pairwise F ST calculated by Arlequin 3.5.1 [61], using the Neighbor-Joining method implemented in the program MEGA Version 4.1 [62]. Significance levels of pairwise F ST were computed after 1000 random permutations Slatkin’s linearized pairwise F ST values were compared with geographical distances between populations according to the Mantel’s test [63] with the implementation present by Arlequin 3.5.1 [61]; normalized Mantel Z-statistics were calculated after 1000 permutations.

Chloroplast DNA sequences were concatenated and then aligned using ClustalW algorithm resident within BioEdit (version 7.0.4.1; [64]). Haplotype networks were obtained on aligned sequences by NETWORK software ver. 4.6 (Fluxus Technologies Ltd.) by using the Reduced Median method.

Authors’ contributions

MY designed experiments, carried out laboratory and statistical analysis, and drafted the manuscript. BL designed experiments and carried out the molecular marker studies. JTL, MH, WCL performed the analyses, AM constructed sequence alignment and helped to draft the manuscript. WSS designed and coordinated the project and drafted the manuscript. All authors discussed data and drafted the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We thank Professor Alan Baker of the University of Melbourne for his help in improving the quality of the manuscript. We are grateful to the editor and three anonymous referees for their insightful comments on earlier versions of this manuscript. This research was financially supported by the National Natural Science Foundation of China (No. 4093212 and 30970548), the Hi-Tech Research and Development Program (863) of China (2007AA061001), and the National Environmental Public Welfare Research Project of China (No. 200909065).

References

  1. Meyer CL, Vitalis R, Saumitou-Laprade P, Castric V: Genomic pattern of adaptive divergence in Arabidopsis halleri, a model species for tolerance to heavy metal.

    Mol Ecol 2009, 18:2050-2062. PubMed Abstract | Publisher Full Text OpenURL

  2. Słomka A, Sutkowska A, Szczepaniak M, Malec P, Mitka J, Kuta E: Increased genetic diversity of Viola tricolor L. (Violaceae) in metal-polluted environments.

    Chemosphere 2011, 83:435-442. PubMed Abstract | Publisher Full Text OpenURL

  3. Kazakou E, Adamidis GC, Baker AJM, Reeves RD, Godino M, Dimitrakopoulos PG: Species adaptation in serpentine soils in Lesbos Island (Greece): metal hyperaccumulation and tolerance.

    Plant Soil 2010, 332:369-385. Publisher Full Text OpenURL

  4. Antonovics J, Bradshaw AD, Turner RG: Heavy metal tolerance in plants.

    Adv Ecol Res 1971, 7:1-85. OpenURL

  5. Whiting SN, Reeves RD, Richards D, Johnson MS, Cooke JA, Malaisse F, Paton A, Smith JAC, Angle JS, Chaney RL, Ginocchio R, Jaffré T, Johns R, McIntyre T, Purvis OW, Salt D, Schat H, Zhao FJ, Baker AJM: Research priorities for conservation of metallophyte biodiversity and their potential for restoration and site remediation.

    Restor Ecol 2004, 12:106-116. Publisher Full Text OpenURL

  6. Baker AJM, Ernst WHO, Van der Ent A, Malaisse F, Ginocchio R: Metallophytes: the unique biological resource, its ecology and conservational status in Europe,Central Africa and Latin America. Edited by Batty LC, Hallberg K. London: Cambridge University Press; 2010:7-40. OpenURL

  7. Macnair MR, Tilstone GH, Smith SE: The genetics of metal tolerance and accumulation in higher plants. Edited by Terry N, Bañuelos G. Boca Raton, Florida: CRC Press; 2000:235-250. OpenURL

  8. Lefèbvre C, Vernet P: Microevolutionary processes on contaminated deposits. Edited by Shaw J. Boca Raton, Florida: CRC Press; 1990:286-297. OpenURL

  9. Hanikenne M, Nouet C: Metal hyperaccumulation and hypertolerance: a model for plant evolutionary genomics.

    Curr Opin Plant Biol 2011, 14:252-259. PubMed Abstract | Publisher Full Text OpenURL

  10. Bradshaw AD: Populations of Agrostis tenuis resistant to lead and zinc poisoning.

    Nature 1952, 169:1098. PubMed Abstract | Publisher Full Text OpenURL

  11. Bradshaw AD, McNeilly T: Evolution and Pollution. London: Edward Arnold (Publishers) Limited; 1981. OpenURL

  12. Macnair MR: The genetics of metal tolerance in vascular plants.

    New Phytol 1993, 124:541-559. Publisher Full Text OpenURL

  13. Macnair MR: The evolution of plants in metal-contaminated environments. Edited by Bijlsma R, Loeschcke V. Switzerland: Verlag, Birkhäuser; 1997:1-24. OpenURL

  14. Nordal I, Haraldsen KB, Ergon A, Eriksen AB: Copper resistance and genetic diversity in Lychnis alpina (Caryophyllaceae) populations on mining sites.

    Folia Geobot 1999, 34:471-481. Publisher Full Text OpenURL

  15. Mengoni A, Barabesi C, Gonnelli C, Galardi F, Gabbrielli R: Genetic diversity of heavy metal-tolerance populations in Silene paradoxa L. (Caryophyllaceae): a chloroplast microsatellite analysis.

    Mol Ecol 2001, 10:1909-1916. PubMed Abstract | Publisher Full Text OpenURL

  16. Muller LAH, Lambaerts M, Vangronsveld J, Colpaert JV: AFLP-based assessment of the effects of environmental heavy metal pollution on the genetic structure of pioneer populations of Suillus luteus.

    New Phytol 2004, 164:297-303. Publisher Full Text OpenURL

  17. Pauwels M, Saumitou-Laprade P, Holl AC, Petit D, Bonnin I: Multiple origin of metalliferous populations of the pseudometallophyte Arabidopsis halleri (Brassicaceae) in Central Europe: the cpDNA testimony.

    Mol Ecol 2005, 14:4403-4414. PubMed Abstract | Publisher Full Text OpenURL

  18. Schat H, Vooijs R, Kuiper E: Identical major gene loci for heavy metal tolerances that have independently evolved in different local populations and subspecies of Silene vulgaris.

    Evolution 1996, 50:1888-1895. Publisher Full Text OpenURL

  19. Koch M, Mummenhoff K, Hurka H: Systematics and evolutionary history of heavy metal tolerant Thlaspi caerulescens in Western Europe: evidence from genetic studies based on isozyme analysis.

    Biochem Syst Ecol 1998, 26:823-838. Publisher Full Text OpenURL

  20. Mengoni A, Gonnelli C, Hakvoort HW, Galardi F, Bazzicalupo M, Gabbrielli R, Schat H: Evolution of copper-tolerance and increased expression of a 2b-type metallothionein gene in Silene paradoxa L. populations.

    Plant Soil 2003, 257:451-457. OpenURL

  21. Turner TL, Bourne EC, Von Wettberg EJ, Hu TT, Nuzhdin SV: Population resequencing reveals local adaptation of Arabidopsis lyrata to serpentine soils.

    Nat Genet 2010, 42:260-263. PubMed Abstract | Publisher Full Text OpenURL

  22. Tang SR, Wilke BM, Huang CY: The uptake of copper by plants dominantly growing on copper mining spoils along the Yangtze River, the People’s Republic of China.

    Plant Soil 1999, 209:225-232. Publisher Full Text OpenURL

  23. Shu WS, Yang KY, Zhang ZQ, Yang B, Lan CY: Flora and heavy metals in dominant plants growing on an ancient Cu spoil heap on Tonglushan in Hubei province, China.

    Chin J Appl Environ Biol 2001, 7:7-12. OpenURL

  24. Ye M, Li JT, Tian SN, Hu M, Yi S, Liao B: Biogeochemical studies of metallophytes from four copper-enriched sites along the Yangtze River, China.

    Environ Geol 2009, 56:1313-1322. Publisher Full Text OpenURL

  25. Barrutia O, Epelde L, Garcia-Plazaola JI, Garbisu C, Becerril JM: Phytoextraction potential of two Rumex acetosa L. accessions collected from metalliferous and non-metalliferous sites: effect of fertilization.

    Chemosphere 2009, 74:259-264. PubMed Abstract | Publisher Full Text OpenURL

  26. Barrutia O, Garbisu C, Hernandez-Allica J, Garcia-Plazaola JI, Becerril JM: Differences in EDTA-assisted metal phytoextraction between metallicolous and non-metallicolous accessions of Rumex acetosa L.

    Environ Pollut 2010, 158:1710-1715. PubMed Abstract | Publisher Full Text OpenURL

  27. Li AR, Gao ZJ, Mao ZM, Liu YL: Polygonaceae. In Flora Reipublicae Popularis Sinicae, Tomus 25 (1). Edited by Delectis Florae Reipublicae Popularis Sinicae Agendae Academiae Sinicae Edita. Beijing: Science Press; 1998:150-151. OpenURL

  28. Wu GF, Ma WL, Hong DY, Ji ZH: Commelinaceae. In Flora Reipublicae Popularis Sinicae, Tomus 13. Editied by Delectis Florae Reipublicae Popularis Sinicae Agendae Academiae Sinicae Edita. Beijing: Science Press; 1997:125-127. OpenURL

  29. Dubois S, Cheptou PO, Petit C: Genetic structure and mating systems of metalliferous and nonmetalliferous populations of Thlaspi caerulescens.

    New Phytol 2003, 157:633-641. Publisher Full Text OpenURL

  30. Molitor M, Dechamps C, Gruber W, Meerts P: Thlaspi caerulescens on nonmetalliferous soil in Luxembourg: ecological niche and genetic variation in mineral element composition.

    New Phytol 2005, 165:503-512. PubMed Abstract | Publisher Full Text OpenURL

  31. Jiménez-Ambriz G, Petit C, Bourrié I, Dubois S, Olivieri I, Ronce O: Life history variation in the heavy metal tolerant plant Thlaspi caerulescens growing in a network of contaminated and noncontaminated sites in southern France: role of gene flow, selection and phenotypic plasticity.

    New Phytol 2007, 173:199-215. PubMed Abstract | Publisher Full Text OpenURL

  32. Besnard G, Basic N, Chiristin P, Savova-Bianchi D, Galland N: Thlaspi caerulescens (Brassicaceae) population genetics in western Switzerland: is the genetic structure affected by natural variation of soil heavy metal concentrations?

    New Phytol 2009, 181:974-984. PubMed Abstract | Publisher Full Text OpenURL

  33. Mohtadi A, Ghaderian SM, Schat H: A comparison of lead accumulation and tolerance among heavy metal hyperaccumulating and non-hyperaccumulating metallophytes.

    Plant Soil 2012, 352:267-276. Publisher Full Text OpenURL

  34. Bert V, Macnair MR, de Laguerie P, Saumitou-Laprade P, Petit D: Zinc tolerance and accumulation in metalliferous and non-metalliferous populations of Arabidopsis halleri (Brassicaceae).

    New Phytol 2000, 146:225-233. Publisher Full Text OpenURL

  35. Van Rossum F, Bonnin I, Fenart S, Pauwels M, Petit D, Saumitou-laprade P: Spatial genetic structure within a metalliferous population of Arabidopsis halleri, a clonal, self-incompatible and heavy-metal-tolerant species.

    Mol Ecol 2004, 13:2959-2967. PubMed Abstract | Publisher Full Text OpenURL

  36. Pauwels M, Frérot H, Bonnin I, Saumitou-Laprade P: A broad-scale analysis of population differentiation for Zn tolerance in an emerging model species for tolerance study: Arabidopsis halleri (Brassicaceae).

    J Evol Biol 2006, 19:1838-1850. PubMed Abstract | Publisher Full Text OpenURL

  37. Quintela-Sabarís C, Ribeiro MM, Poncet B, Costa R, Castro-Fernández D, Fraga MI: AFLP analysis of the pseudometallophyte Cistus ladanifer: comparison with cpSSRs and exploratory genome scan to investigate loci associated to soil variables.

    Plant Soil 2012. Publisher Full Text OpenURL

  38. Mengoni A, Gonnelli C, Galardi F, Gabbrielli R, Bazzicalupo M: Genetic diversity and heavy metal tolerance in populations of Silene paradoxa L. (Caryophyllaceae): a random amplified polymorphic DNA analysis.

    Mol Ecol 2000, 9:1319-1324. PubMed Abstract | Publisher Full Text OpenURL

  39. Bratteler M, Lexer C, Widmer A: Genetic architecture of traits associated with serpentine adaptation of Silene vulgaris.

    J Evol Biol 2006, 19:1149-1156. PubMed Abstract | Publisher Full Text OpenURL

  40. Antonovics J, Bradshaw AD: Evolution in closely adjacent plant populations. VIII. Clinal patterns at a mine boundary.

    Heredity 1970, 25:349-362. Publisher Full Text OpenURL

  41. Shu WS, Ye ZH, Lan CY, Zhang ZQ, Wong MH: Lead, zinc and copper accumulation and tolerance in populations of Paspalum distichum and Cynodon dactylon.

    Environ Pollut 2002, 120:445-453. PubMed Abstract | Publisher Full Text OpenURL

  42. Yang XE, Long XX, Ye HB, He ZL, Calvert DV, Stoffella PJ: Cadmium tolerance and hyperaccumulation in a new Zn-hyperaccumulating plant species (Sedum alfredii Hance).

    Plant Soil 2004, 259:181-189. OpenURL

  43. Galardi F, Corrales I, Mengoni A, Pucci S, Barletti L, Arnetoli M, Gabbrielli R, Gonnelli C: Intra-specific differences in nickel tolerance and accumulation in the Ni-hyperaccumulator Alyssum bertolonii.

    Environ Exp Bot 2007, 60:377-384. Publisher Full Text OpenURL

  44. Meyer CL, Kostecka AA, Saumitou-Laprade P, Créach A, Castric V, Pauwels M, Frérot H: Variability of zinc tolerance among and within populations of the pseudometallophyte species Arabidopsis halleri and possible role of directional selection.

    New Phytol 2010, 185:130-142. PubMed Abstract | Publisher Full Text OpenURL

  45. Baker AJM, Walker PL: Physiological responses of plants to heavy metals and the quantification of tolerance and toxicity.

    Chem Spec Bioavail 1989, 1:7-17. OpenURL

  46. Mengoni A, Selvi F, Cusimano N, Galardi F, Gonnelli C: Genetic diversity inferred from AFLP fingerprinting in populations of Onosma echioides (Boraginaceae) from serpentine and calcareous soils.

    Plant Biosyst 2006, 140:211-219. Publisher Full Text OpenURL

  47. Rieseberg LH, Willis JH: Plant speciation.

    Science 2007, 317:910-914. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Antonovics J: Evolution in closely adjacent plant populations X: long-term persistence of prereproductive isolation at a mine boundary.

    Heredity 2006, 97:33-37. PubMed Abstract | Publisher Full Text OpenURL

  49. Deng JC, Liao B, Ye M, Deng DM, Lan CY, Shu WS: The effects of heavy metal pollution on genetic diversity in zinc/cadmium hyperaccumulator Sedum alfredii populations.

    Plant Soil 2007, 297:83-92. Publisher Full Text OpenURL

  50. Allen SE: Chemical Analysis of Ecological Materials (2nd ed). Oxford: Blackwell Scientific Publications; 1989. OpenURL

  51. Lindsay WL, Norvell WA: Development of a DTPA test for zinc, iron, manganese, and copper.

    Soil Sci Soc Am J 1978, 42:421-428. Publisher Full Text OpenURL

  52. Lodhi MA, Ye GN, Weeden NF, Daly MJ, Reisch BI: A simple and efficient method for DNA extraction from grapevine cultivars, Vitis species and Ampelopsis.

    Plant Mol Biol Rep 1994, 12:6-13. Publisher Full Text OpenURL

  53. Vos P, Hogers R, Bleeker M, Reijans M, Zabeau M: AFLP: a new technique for DNA fingerprinting.

    Nucl Aci Res 1995, 23:4407-4414. Publisher Full Text OpenURL

  54. Bassam BJ, Caetano-Anollés G, Gresshoff PM: Fast and sensitive silver staining of DNA in polyacrylamide gels.

    Anal Biochem 1991, 196:80-83. PubMed Abstract | Publisher Full Text OpenURL

  55. Yakimowski SB, Eckert CG: Populations do not become less genetically diverse or more differentiated towards the northern limit of the geographical range in clonal Vaccinium stamineum (Ericaceae).

    New Phytol 2008, 180:534-544. PubMed Abstract | Publisher Full Text OpenURL

  56. Shaw J, Lickey EB, Schilling EE, Small RL: Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: the tortoise and the hare III.

    Am J Bot 2007, 94:275-288. PubMed Abstract | Publisher Full Text OpenURL

  57. Yeh FC, Yang RC, Boyle T, Ye ZH, Mao JX: POPGENE, the User-friendly Shareware for Population Genetic Analysis. Alberta: Molecular Biology and Biotechnology Centre, University of Alberta; 1997. OpenURL

  58. Nei M: Analysis of gene diversity in subdivided populations.

    Proc Natl Acad Sci USA 1973, 70:3321-3323. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  59. Shannon CE, Weaver W: The Mathematical Theory of Communication. Urbana: University of Illinois Press; 1949. OpenURL

  60. Excoffier L, Smouse PE, Quattro M: Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction site data.

    Genetics 1992, 131:479-491. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  61. Excoffier L, Laval G, Schneider S: Arlequin ver. 3.0: An integrated software package for population genetics data analysis.

    Evol Bioinform Online 2005, 1:47-50. OpenURL

  62. Kumar S, Tamura K, Nei M: MEGA3: Integrated software for molecular evolutionary genetics analysis and sequence alignment.

    Brief Bioinform 2004, 5:150-163. PubMed Abstract | Publisher Full Text OpenURL

  63. Mantel N: The detection of disease clustering and a generalized regression approach.

    Cancer Res 1967, 27:209-220. PubMed Abstract | Publisher Full Text OpenURL

  64. Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT.

    Nucl Aci Symp Ser 1999, 41:95-98. OpenURL