Figures
Abstract
Simple sequence repeat (SSR) markers are widely used tools for inferences about genetic diversity, phylogeography and spatial genetic structure. Their applications assume that variation among alleles is essentially caused by an expansion or contraction of the number of repeats and that, accessorily, mutations in the target sequences follow the stepwise mutation model (SMM). Generally speaking, PCR amplicon sizes are used as direct indicators of the number of SSR repeats composing an allele with the data analysis either ignoring the extent of allele size differences or assuming that there is a direct correlation between differences in amplicon size and evolutionary distance. However, without precisely knowing the kind and distribution of polymorphism within an allele (SSR and the associated flanking region (FR) sequences), it is hard to say what kind of evolutionary message is conveyed by such a synthetic descriptor of polymorphism as DNA amplicon size. In this study, we sequenced several SSR alleles in multiple populations of three divergent tree genera and disentangled the types of polymorphisms contained in each portion of the DNA amplicon containing an SSR. The patterns of diversity provided by amplicon size variation, SSR variation itself, insertions/deletions (indels), and single nucleotide polymorphisms (SNPs) observed in the FRs were compared. Amplicon size variation largely reflected SSR repeat number. The amount of variation was as large in FRs as in the SSR itself. The former contributed significantly to the phylogenetic information and sometimes was the main source of differentiation among individuals and populations contained by FR and SSR regions of SSR markers. The presence of mutations occurring at different rates within a marker’s sequence offers the opportunity to analyse evolutionary events occurring on various timescales, but at the same time calls for caution in the interpretation of SSR marker data when the distribution of within-locus polymorphism is not known.
Citation: Barthe S, Gugerli F, Barkley NA, Maggia L, Cardi C, Scotti I (2012) Always Look on Both Sides: Phylogenetic Information Conveyed by Simple Sequence Repeat Allele Sequences. PLoS ONE 7(7): e40699. https://doi.org/10.1371/journal.pone.0040699
Editor: Shin-Han Shiu, Michigan State University, United States of America
Received: February 7, 2012; Accepted: June 12, 2012; Published: July 13, 2012
This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Funding: Funding provided by PO-FEDER Guyane “ENERGIRAVI” http://www.europe-guyane.eu/index.php?option=com_content&task=view&id=14&Itemid=206. SB’s Doctoral program has been partly funded by the European Social Fund. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Simple sequence repeats (SSRs; see glossary, Table 1) can conveniently be used as genetic markers owing to their polymorphism and relative ease of interpretation [1]. SSRs are widely used for, e.g., reconstructing phylogenetic relationships [2], for analysing spatial genetic structure among and within populations [1], and for detecting and explaining patterns linked to habitat fragmentation and gene flow [3]. In many studies, SSR markers demonstrated recent expansions or bottlenecks in various plant [4], [5], animal [6]–[8] and human populations [9], [10].
SSRs are often expected to mutate following the stepwise mutation model (SMM) [11], whereby mutations alter the length of the repeat either by adding or by deleting a single repeat unit at a fixed rate [12]. More elaborate mutational processes, such as the two-phase model (TPM; [13]) and the generalized stepwise model (GSM; [14], [15]), allow for multi-step mutations. Departures from these models can however occur [16]–[20]. Moreover, SSR alleles are generally scored as the length (in base pairs) of PCR amplicons (see glossary, Table 1) containing the SSR, with differences in amplicon size taken to represent differences in repeat number in the SSR. Amplicon size includes of course SSR repeat number plus the length of the flanking region (FR; see glossary, Table 1). These chunks of sequence can be and actually often are polymorphic and may contain both single nucleotide polymorphisms (SNPs; see glossary, Table 1) and insertions/deletions (indels; see glossary, Table 1) [18], [21]–[24]. Indels which have been demonstrated to occur in the FR sequences of some SSR alleles clearly contribute to total amplicon size, but they are interpreted under the SMM as (false) variation in SSR length. SNPs contribute to the total sequence variance of SSR loci (see glossary, Table 1) that goes largely unnoticed when only amplicon size is recorded (which is typical of most SSR studies). As a consequence, SSR amplicon size data are prone to a particular form of size homoplasy [15], i.e. equally sized alleles may have different sequences and therefore be evolutionarily different. This casts doubt on the evolutionary and population-genetic inferences that can be drawn from estimates of divergence obtained from SSR amplicon size data. It is therefore necessary to examine the kinds of evolutionary information conveyed in SSR alleles, to assess their consequences on fragment size-based evolutionary inferences, and to find ways to take advantage of this so far unexploited source of variability.
Figure 1 succinctly shows how information in FR sequence variation may change one’s view of SSR marker evolutionary information content. Figure 1A shows a “plain” SSR locus, where a mutation (bar on the right branch) changes the number of repeats from n to n+1. The four observed alleles are then linked by a simple genealogy where alleles with the same number of repeats are evolutionarily closer. This is the classical view of SSR allele evolution under repeat-number mutational models. Figure 1B shows four alleles with the same numbers of repeats as in Figure 1A, but in this case information from a SNP in the FR (indicated by the “…A/…C” symbols) is added. Here, alleles with different numbers of repeats share the same SNP. Given that SNP mutation rates are generally lower than those of SSRs, the most-likely genealogy is one in which divergence occurs first at the SNP (mutation in the upper right branch) and then in SSRs within flanking-region haplotypes (mutations in the lower branches). In this case, three mutations have occurred, the genealogy is expected to be older than in 1A, and the alleles are grouped based on their FR sequences and not based on the number of repeats. A simulation study [25] revealed that SSR and SNP sources of variation provide independent evolutionary information, so it is reasonable to study the role of SNP variation in SSR flanking regions. We therefore asked the following questions: How much does amplicon size variation reflect variation in SSR repeat number? How important is the contribution of the various sources of sequence variation in total SSR diversity? Do these extra-SSR polymorphisms carry a different population and phylogenetic differentiation signal than the SSRs themselves? If so, how do the different signals combine? How does one best interpret SSR data, once these departures from the SMM have been acknowledged and assessed? To answer these questions, we sequenced SSR amplicons obtained from three divergent genera of long-lived angiosperm tree species to inspect the distribution of polymorphisms within and around SSRs and to evaluate the impact of each source of sequence variation on the detection of evolutionary and population genetic signals.
SSR alleles are indicated by their number of repeats (n, n+1). The A/C letters indicate a SNP in the flanking regions. Red bars correspond to mutational events in the flanking region sequence or in the number of SSR repeats. (A) No sequence information: deduced genealogy of observed data (third line) groups alleles together according to their number of repeats and involves a single SSR mutation. Genealogy is recent. (B) Consideration of sequence information: deduced genealogy involves a SNP mutation and two SSR mutations (alternative topology will involve two identical and independent substitutions at the same nucleotide site and a unique SSR mutation, which is less likely). Genealogy is ancient and SSR alleles do not group according to their numbers of repeats.
Materials and Methods
Choice of Data Sets
Three data sets were utilized in this study. The structures of the three data sets analysed here varied in the following ways: the Citrus data represent sequences of individual trees each belonging to a different species, variety or cultivar; the Jacaranda data were collected in four different sampling sites considered as four distinct populations; and the oak data came from a single mixed stand of sessile/pedunculate oaks (Quercus robur, Q. petraea). The Citrus, Jacaranda, and Quercus data sets were submitted to partially different types of analyses according to the structure of each data set.
Given the taxonomy and prevailing theory on many Citrus species being derived by natural hybridization from ancestral species, alleles were chosen to minimize known hybrids from the data set of Barkley et al. [26], [27], who examined the genetic diversity of 370 Citrus accessions. The individual identity of each accession (as provided by the University of California, Riverside citrus database, http://www.citrusvariety.ucr.edu/) is given in the Supplementary Table 1 of Barkley et al. [26]. The Jacaranda data were obtained from three sites in French Guiana (Counami, Paracou and Saint-Laurent) and one in Brazil (Tapajos) (Table S1). The Quercus alleles were taken from genetic analyses of Q. petraea and Q. robur in a mixed oak stand in Switzerland [28], [29]. The SSR markers used for this study were: cAGG9, CCT01 and GT03 for Citrus [26], [30]; Jc3H10, Jc3F4 and Jc3H10 for Jacaranda; QrZAG30 for Quercus [31]. Allele sequencing for the Citrus and Quercus data sets was as reported in Barkley et al. [26] and Gugerli et al. [29].
Molecular Methods
Jacaranda markers were obtained following the method developed for SSR markers in tropical crops by Billotte et al. [32]. Ten nanograms of DNA was digested with RsaI, followed by an enrichment step in (GT)n and (GA)n repeats by hybridization of cleaved fragments with biotin-labelled (CA)8 and (CT)8 oligoprobes and capture of the selected sequences with streptavidincoated magnetic beads. Those fragments were ligated into a pGEM-TEasy vector (Promega, Madison, WI), and plasmids were then used to transform competent XL-1 Blue strain Escherichia coli (Stratagene, La Jolla, CA). Positive colonies were tested using colony polymerase chain reaction (PCR) to check the presence of inserts. Amplified fragments were then transferred from an agarose gel onto a N+ Hybond membrane for SSR presence screening by hybridization with labeled γ−32 P (GT)15 and (GA)15 oligoprobes. One hundred and thirteen bacterial clones containing plasmids with inserts that gave strong hybridization signals, with sizes ranging from 500 to 1000bp, were selected and cultured. Plasmids were extracted and sequenced using the universal T7 primer on an automated ABI 377 capillary sequencer (Applied Biosystems, Foster City, CA). After discarding duplicates, hybrid clones and clones with the SSR region too close to either end of the sequence, thirteen sequences were suitable for primer design and allowed a successful design using Oligo 3 software (Molecular Biology Insights Inc, USA). Upon screening for polymorphism and clarity of PCR patterns, three primer pairs were chosen for subsequent analyses. Primer sequences, annealing temperatures and GenBank/EMBL accession numbers are shown in Table S2. Dried cambium discs and leaf samples were flash-frozen in liquid nitrogen and later ground to a powder using a mortar and pestle. DNA extractions were performed by following Colpaert et al. [33]. PCRs for the detection of SSR polymorphisms were carried out in a 12 µL volume containing 6 µL 20-fold diluted DNA, 1× Taq buffer, 0.26 mM of dNTP, 0.03 U/µL Taq DNA polymerase (all products from Invitrogen, Carlsbad, CA) and 0.54 µM of each primer (MWG Biotech, Ebersberg, Germany). For Jc3A10 primers, 0.3% BSA was added. An initial denaturation at 94°C for 5 min was followed by 35 cycles of 94°C for 30 s, annealing temperature for 30 s and 72°C for 30 s, and a final extension at 72°C for 10 min. Genotyping was performed using fluorescently labeled primers (PET, 6FAM and NED) in the previous PCR protocol and fragments were separated on an ABI 3130XL capillary sequencer (Applied Biosystems, Foster City, CA) using ABI POP4 and Applied Biosystems LIZ-500 as internal standard, following the manufacturer’s instructions. DNA was bidirectionally sequenced directly from PCR amplification products. When necessary, gametic phase was determined in heterozygotes by cloning amplicons and sequencing at least one allele per individual following the TA cloning kit protocol from Invitrogen (Carlsbad, CA). PCR products were ligated with the plasmid pCR®2.1 and used to transform competent DH5 cells. After an overnight incubation at 37°C, white colonies were isolated. The presence of an SSR allele was checked by PCR with universal M13 primers. PCRs consisted of 1 µL of cultivated colonies, 1× Taq buffer, 0.125 mM dNTPs, 0.025 U/µL Taq DNA polymerase (all products from Invitrogen, Carlsbad, CA) and 0.83 µM forward and reverse primers. Cycling conditions consisted of 94°C for 5 min; 35 cycles of: 94°C for 1 min, 50°C for 1 min, 72°C for 2 min; and one cycle of 72°C for 5 min. PCR products were checked on a 4% agarose gel in 0.5× TAE. Plasmids were sequenced with the ABI BigDye® Terminator V3.1 kit (Applied Biosystems, Foster City, CA) following the manufacturer’s protocol. Cloning products were diluted 1∶40, purified with ExoSap-IT (USB Corporation, Cleveland, OH), and separated on an ABI 3130XL capillary sequencer (Applied Biosystems, Foster City, CA). All sequences were aligned and edited using CodonCode Aligner V1.6.3 (Codoncode Corporation, Dedham, MA).
Data Analyses
For each data set (Citrus, Jacaranda, and Quercus), each allele was characterised by the following information: (a) amplicon size variation, (b) SSR variation, (c) FR sequence variation and (d) amplicon sequence variation. Each of these portions of information conveyed by SSR marker data was analysed as a separate source of variation. Results were compared across data sets to assess their evolutionary information content. Indels were coded as SNP mutations to represent them as single mutational events.
Levels of polymorphism (number of alleles (A) or haplotypes (h), Nei’s genetic diversity (He), number of SNPs and indels) were recorded for each source of variation. Linkage disequilibrium (LD) was computed for FR sequence variation.
Data were analysed in two ways: (a) for all data sets, at the level of individual alleles; (b) for Jacaranda only, at the population level.
- At the individual-allele level, matrices of pairwise genetic distances were computed between individuals and between alleles based on each source of variation. For amplicon length, distance is represented by difference in length; for SSR length, distance is the difference in number of repeat units; for whole amplicon sequence, distance is the total number of differences, including SNPs, indels and number of repeat units; for flanking region sequences, distance is the total number of SNPs and indels. The correlation between matrices was tested by a Mantel test. LD between all pairs of polymorphic sites in FR sequences was tested with a Markov chain Monte Carlo (MCMC) procedure of 10,000 steps and a burn-in phase of 1,000 steps (default values). Pairwise genetic distances for amplicon and FR sequence variation were computed based on the number of mutational steps between sequences. Pairwise genetic distances for SSR variation and amplicon size variation data were computed as Euclidian distances (genetic distances based on the SSR variation were not computed for Citrus loci CCT01 and GT03, which displayed extremely complex repeat sequences; see Figure S1). Most of the above calculations were performed with ARLEQUIN versions 3.11 and 3.5.1.2 [34]. Genetic distance matrices for SSR variation and amplicon size variation data were computed with an ad-hoc routine in R [35], and correlations between genetic distance matrices were computed with the mantel.test function of the R package “ape”.
- At the population level, for Jacaranda, global and pairwise FST values [36] were computed for all sources of variation; RST values [37] were computed on sources of variability comprising amplicon size variation and SSR variation; NST values [38] were computed on amplicon and FR sequence variations. In NST calculations on amplicon sequence variation, mutations in repeat number were weighted less than substitutions (weight = 0.05) to account for their faster mutation rate. Weight choice was based on the observation that SSRs have, in our Jacaranda sample, approximately ten times more alleles than each SNP (see Results section). Considering that rare alleles might be missing from our sample, we estimated a ratio of 1∶20 in allele richness between SNPs and SSRs and used this ratio to establish the relative weight. A hierarchical analysis of molecular variance (AMOVA) was performed to partition genetic variance into within- and among-population components. The significance of variance components was tested by 1,000 permutations [39]. These analyses were performed with ARLEQUIN 3.11 and 3.5.3.1. Locus-by-locus pairwise FST, RST, and NST matrices were used to construct population-level consensus UPGMA (Unweighted Pair Group Method with Arithmetic Mean) cladograms by averaging genetic distance information of three loci with SplitsTree4 v4.6 [40], [41].
Results
Distribution of Polymorphism
In order to evaluate homoplasy and the variation in phylogenetic signal from different regions of SSR alleles (repeat region and FR), sequence data produced from three divergent tree genera were evaluated. Forty-six, 101, and 47 samples constituted the Citrus, Jacaranda, and Quercus data sets, respectively. Three, three and one SSRs were genotyped and sequenced for the Citrus, Jacaranda and Quercus data sets, respectively. The lengths of the sequenced fragments varied between 123 and 255 bp (Table 2). Whenever possible, multiple alleles of the same fragment size were genotyped and sequenced to estimate the degree of sequence divergence among equally sized alleles (“size homoplasy”). This includes both variation in SSR repeat number (that is taken into account by the SMM) and variation in FR sequences, which is the focus of our study.
Numbers in the top line indicate positions relative to the consensus sequence of (A) 190 bp and (B) 247 bp. Bold nucleotides in brackets indicate SSR motifs and their number of repeats. Dashes indicate gaps, highlighted nucleotides in green (light grey*) indicate indels of one to three bases, highlighted nucleotides in red (dark grey*) mark mutations from one base to another, and yellow (light grey*) boxes indicate groups of insertion/deletions longer than three bases and considered as a single mutational event. * In shades-of-gray printouts.
Levels of polymorphism (number of alleles/haplotypes, Nei’s genetic diversity, number of SNPs and indels) for each amplicon and each source of variation are displayed in Table 3. Large amounts of variability were observed, as expected, in the amplicon sequences. SNPs were found to interrupt the repeat in all three taxa (imperfect SSRs: 33 sequences, data not shown). Sequences with SNPs within the repeat were excluded from subsequent analyses owing to the complexity of the mutation model of SSR repeats containing SNPs. Four loci had two or more repeats (compound SSRs). The highest number of alleles and Nei’s genetic diversity were detected for amplicon sequence variation, then for amplicon size variation and lowest was found for SSR variation (Table 3). FR sequence variation showed the lowest diversity. Indels were observed in six loci out of seven, while SNP density varied between 0.90 and 12.14 SNPs per 100 bp.
Dark cells contain LD values significantly different from zero. An asterisk next to a locus name indicates that the SSR repeat(s) is (are) located between that locus and the next. Colouring (shading) indicates the degree of significance of the test: green (light grey*), P>0.05; orange (grey*), 0.05<P>0.001; red (dark grey*), P<0.001. * In shades-of-gray printouts.
An example of the complexity of the polymorphisms observed is shown for markers Jc3H10 and QrZAG30 (Jacaranda and Quercus data set, Figure 2; the complete alignments for all sequenced amplicons are given in Figure S1). Three indels of several nucleotides (between 12 and 15 bases) were detected in the FR sequences of the compound marker Jc3H10. One was upstream of the first SSR, the second was between the repeat motifs and the last was detected after the SSR repeats. In the oak marker QrZAG30, in addition to SNPs, a large 66-nucleotide indel was observed, with SNPs occurring among alleles carrying the DNA fragment involved in the indel. Several cases of size homoplasy were detected, such as: (i) indels in the FRs compensating differences in number of repeats, (ii) compound SSRs with the same number of repeats, but composed of different numbers of repeats in different motifs, (iii) SSRs with the same amplicon size and number of repeats but with SNPs in the FR sequences or interrupting the repeat (Figure S1). For the three data sets, 50 alleles with the same amplicon sizes out of 68 (74%) showed size homoplasy (Table S3).
Linkage Disequilibrium
LD was tested for all pairs of SNP and indel sites within the FR sequence variation of each marker (Figure 3). For the Citrus data set, only 19 pairs out of 217 (8.8%) showed significant LD. For the three SSR loci, LD was very irregular and there was no relationship between LD and sequence distance (Mantel tests, P>0.05 for all loci). For the Jacaranda data set, ten pairs out of 27 (37%) showed significant LD. For marker Jc3F4, the SNP at position 81 was in disequilibrium with two indels and one SNP, all located upstream of the SSR. LD was strongest for Jc3H10, where all three pairs were significant. The relationship between LD and sequence distance was significant (Mantel test, P = 0.023). Tests of correlation between LD and distance in bp could not be performed for Jc3A10 and Jc3H10 because only three nucleotide sites were polymorphic (three LD values). For the Quercus data set (marker QrZAG30), 146 pairs of SNPs and indels out of 325 (44.9%) showed significant LD. In particular, 13 or 14 contiguous SNPs and indels located downstream of the SSR at positions 84–107 formed a disequilibrium block, itself in weaker disequilibrium with upstream polymorphic sites at positions 23, 45 and 46 and another SNP downstream the SSR (position 60).
SSR variation (A), amplicon size variation (B), flanking region (FR) sequence variation (C) and amplicon sequence variation (D). Each type of data was analysed according to a pair of suitable genetic distance estimators: FST, suitable for loci following the infinite allele model (IAM); RST, for loci following the stepwise mutation model (SMM); and NST, for loci following the infinite site model (ISM) model. There are four geographic populations: Counami (C), Paracou (P) and Saint-Laurent (S) in French Guiana; Tapajos (T) in Brazil. Note that scales are not the same among trees.
Correlation between Genetic Distance Matrices
We tested whether different sources of variation conveyed the same information on relatedness among alleles. To do so, genetic distances were computed for each pair of individuals with each source of variation (SSR variation, amplicon size, FR sequence, and amplicon sequence variation) within each marker, and the correlation between genetic distance matrices was tested by a Mantel test. Among the 36 pairs of matrices, 27 (75%) showed a significant correlation (Figure S2). Genetic distance matrices based on SSR variation and on FR sequence variation were the least correlated (one significant pair out of five). Amplicon sequence variation and amplicon size variation, as well as SSR variation and amplicon size variation, were correlated for all markers. The matrices computed on amplicon sequence variation were correlated to the remaining matrices in 95% of the cases (18 out of 19). Conversely, FR sequence variation alone predicted other matrices in only 53% of the cases.
Population Differentiation
For the Jacaranda data set, we tested population genetic differentiation for each source of variation by computing FST, RST and NST as appropriate (see Methods; Figure 4). All FST values were significant. FR sequence variation showed the largest value (0.42) and SSR variation the smallest (0.05). RST values were very small (<0.02) and non-significant when computed on both amplicon size variation and SSR variation. On the contrary, NST values were large and significant, with values of 0.50 for FR sequence variation and 0.36 for amplicon sequence variation. In the phylogenetic trees obtained from each combination of source of variation/genetic distance measures, FST and NST provided the same topology for all cases they were applied to, indicating a more or less close link between pairs of populations Tapajos/Saint-Laurent and Counami/Paracou (Figure 4). RST grouped Paracou with Tapajos and Counami with Saint-Laurent (for amplicon size variation) or separated the three Guiana shield populations from Tapajos (for SSR variation; Figure 4).
Discussion
Distribution of Polymorphism
SSR polymorphisms are widely used to infer population history and biogeographic patterns. These inferences rely on assumptions about SSR mutational models, and departures from such models are likely to bias divergence and diversity estimates. SSR alleles have been demonstrated to contain other sources of polymorphism in addition to variation in repeat number. In this study we have characterised the contribution of molecular variation occurring outside SSR repeats to SSR marker variability. Whole allele sequencing from SSR markers was analyzed from three tree genera (Citrus, Jacaranda and Quercus). A total of seven (dinucleotide and trinucleotide) SSR markers were targeted, including compound and imperfect markers. More genetic diversity was found for amplicon sequence variation than for amplicon size variation, suggesting that much of data polymorphism is neglected when SSRs are described through amplicon size alone. On average, 7.7 SNPs and 3.4 indels were detected in the flanking sequences of each marker. Genetic diversity (He) ranged between 0.306 and 0.922 for the FRs. LD among SNPs and indels in the FR was generally significant. These results demonstrate that SSR alleles are very often riddled with abundant non-SSR indels, plus vast amounts of SNPs both within and near the repeat itself. Similar genetic variation was previously observed at intra- or inter-specific levels [18], [42]–[44] and does not seem to be exceptional in plants. SSR allele mutations are clearly not restricted to the hypervariable tandem repeat region, and polymorphisms outside SSR motifs are also present in animals [22], [45], [46] and humans [21]. Hidden indels can directly confound the estimation of repeat number from amplicon size (“allele size homoplasy”). While undetected SNPs do not directly influencing repeat number estimates, they do cause haplotypes with different evolutionary histories to be treated as if they were identical. Both forms of homoplasy can have profound consequences on the way data should be interpreted.
Correlation of Polymorphism at the Individual-allele Level
Our results show that pieces of genetic information carried by the different portions of an amplicon are only partially correlated. In particular, blocks of linkage disequilibrium never cover the whole amplicon, and matrices of allele-level pairwise genetic distances built on each source of variation are not always well correlated. In other terms, the different sources of variation are not all equally capable to represent the “true” allele genealogy. Not surprisingly, whole-amplicon sequences were the most parsimonious predictors of relatedness described by the different parts of the amplicons in all our data sets; conversely, FR sequences alone did not adequately synthesize the information from the whole amplicon sequence. Correlation between genetic distances as obtained from FR sequence variation and from SSR variation was globally weak (Figure S2). This recalls the finding of Payseur and Cutter [25] which stated that coalescence times for SSRs and linked SNPs were uncorrelated. Thus, analysing all parts of SSR amplicons increases the amount of available, phylogenetically independent information, with the additional advantage of providing data from DNA regions with different mutation rates. Nevertheless, the correlation analysis shows that the best predictor of amplicon size is SSR length. This indicates that amplicon size is a reliable first-approximation proxy for SSR repeat length, even though extensive homoplasy blurs the correlation.
Correlation of Polymorphism at the Population Level
The consequences – and usefulness – of hidden sequence variation are shown when we compare how sources of variation perform in detecting population divergence. This analysis was performed in the Jacaranda data set comprising four populations. In this example, two alternative patterns were expected for population relatedness: (a) a strictly geographical clustering (the three Guiana shield populations cluster and are separated from the Amazonian population), and (b) a pattern derived from independent results on chloroplast DNA divergence (Counami and Paracou populations form a group, Tapajos and Saint-Laurent form another group (Caroline Scotti-Saintagne, INRA, UMR « Ecologie des forêts de Guyane », article under second round of revision for Journal of Biogeography). Our results support the latter hypothesis, as proven by the convergent topology of trees obtained with FST and NST independently of the source of variation. The purely geographical hypothesis, on the other hand, is only supported by SSR variation in combination with RST (although with a non-significant global RST value). RST estimates have larger variance than FST and NST, particularly when based on small numbers of loci, and RST values obtained here may be unreliable. Nevertheless, the pattern of RST differentiation closely follows geographical distance. This may indicate that repeat number variation has arisen locally from small numbers of founding alleles and reflects recent population divergence, while variation in flanking regions may follow larger phylogeographic patterns and reflect deeper divergence. If RST estimated can be trusted, all sources of variation converge with independent, chloroplast-based results, except SSR variation itself. Therefore, information derived from the “true” nature of SSR data (SSR variation) is actually not representative of the information conveyed by alleles as a whole. This finding strongly hints that other sources of variability, more correctly described by mutation models that apply to DNA sequences, actually provide the bulk of information carried by SSRs. In other words, phylogeographic information carried by amplicon length is not provided by the SSR (which amplicon length is intended to represent) but by variation in flanking regions. Note that this contrasts with analyses carried out at the individual level (see above), where information on genetic distance match when measured as amplicon size and as SSR size. This issue raises the question of the application of the SMM, or other step-based mutational models, to SSRs, as already discussed by Colson and Goldstein [47] and Ellegren [12]. At the same time, these results substantiate the hypothesis of Cornuet et al. [48], who postulated that not all SSRs would evolve under the same mutational model; additional mutational models are required to appropriately handle variation in SSR data produced from different genera. Moreover, variation in the FRs contributed significantly to the phylogenetic signal and sometimes represented the main source of differentiation among individuals and populations, as shown by Jacaranda SSR markers and by the Citrus SSR marker GT03 [27].
Conclusion
Given the complexity of amplicon sequences described here, one may wonder how to properly exploit SSR data. The weight of polymorphisms other than variation in repeat number cannot and should not be overlooked: identity in allele size does not necessarily indicate identity in sequence content or the number of repeats within same sized alleles. When considering individual alleles, amplicon size correlates well with SSR repeat number. Nevertheless, even if repeat number could be obtained directly, it may provide only weak and inconsistent signals of population genetic differentiation, as shown by the Jacaranda data set (but larger numbers of markers may reduce estimation variance and mitigate this problem). Taking sequence variation into account actually adds a significant piece of information to phylogenetic or phylogeographic reconstruction. Population structuring emerges more clearly on the basis of sequence data, and the combination of sequences and SSR variation provides higher resolution. Consequently, genetic distance measures assuming the SMM or related models should be restricted to perfect SSRs with invariable FRs. On the other hand, measures such as FST (for amplicon size data) and NST (for amplicon sequences) seem to be the tool of choice for the analysis of the average, not-so-ideal SSR markers that a population geneticist meets in his or her everyday work.
Supporting Information
Figure S1.
Complete alignments of DNA fragments for all SSR markers studied. Marker names are shown above each alignment.
https://doi.org/10.1371/journal.pone.0040699.s001
(PDF)
Figure S2.
Z-statistics and summary of Mantel tests for the correlations of pairs of genetic distance matrices within species. Colouring (shading) indicates the level of significance of the test: green (light grey*), P>0.05; orange (grey*), 0.05<P>0.001; red (dark grey*), P<0.001. Identity of DNA fragment sections is reported: simple sequence repeat (SSR) variation; amplicon size variation, flanking region (FR) sequence variation, and amplicon sequence variation. Marker names are displayed above each panel and data sets refer to Citrus (C), Jacaranda (J) and Quercus (Q). Sequences with single nucleotide polymorphisms (SNPs) within the repeat were excluded from subsequent analyses because of the complexity to describe the mutation model of SSR repeats bracketing SNPs.
https://doi.org/10.1371/journal.pone.0040699.s002
(PDF)
Table S1.
List and details of the sampled sites for the Jacaranda (J) data set.
https://doi.org/10.1371/journal.pone.0040699.s003
(DOC)
Table S2.
Description, amplification conditions and polymorphism of the nuclear simple sequence repeats (SSRs) analysed in Jacaranda copaia .
https://doi.org/10.1371/journal.pone.0040699.s004
(DOC)
Table S3.
List and details of alleles for the three data sets. Amplicon sizes, frequencies and number of associated haplotypes are indicated for each allele.
https://doi.org/10.1371/journal.pone.0040699.s005
(DOC)
Acknowledgments
We wish to thank Maristerra Lemes for providing a part of the samples, Raphaël Leblois and Caroline Scotti-Saintagne for fruitful discussions and for help with preliminary data analyses, and Valérie Troispoux and Saint-Omer Cazal for technical support.
Author Contributions
Conceived and designed the experiments: IS. Performed the experiments: SB FG NB. Analyzed the data: SB IS. Contributed reagents/materials/analysis tools: LM CC. Wrote the paper: SB FG NB IS.
References
- 1. Tautz D (1989) Hypervariability of simple sequences as a general source for polymorphic DNA markers. Nucleic Acids Res 17: 6465–6471.
- 2. Takezaki N, Nei M (1996) Genetic distances and reconstruction of phylogenetic trees from microsatellite DNA. Genetics 144: 389–399.
- 3. Rodriguez H, Geistlinger J, Berlyn G, Kahl G, Weising K (2000) Characterization of novel microsatellite loci isolated from the tropical dioecious tree Simarouba amara. Mol Ecol 9: 489–504.
- 4. England PR, Usher AV, Whelan RJ, Ayre DJ (2002) Microsatellite diversity and genetic structure of fragmented populations of the rare, fire-dependent shrub Grevillea macleayana. Mol Ecol 11: 967–977.
- 5. Taylor DR, Keller SR (2007) Historical range expansion determines the phylogenetic diversity introduced during contemporary species invasion. Evolution 61: 335–345.
- 6. Estoup A, J Wilson IJ, Sullivan C, Cornuet J-M, Moritz C (2001) Inferring population history from microsatellite and enzyme data in serially introduced cane toads, Bufo marinus. Genetics 159: 1671–1687.
- 7. Randi E, Tabarroni C, Rimondi S, Lucchini V, Sfougaris A (2003) Phylogeography of the rock partridge (Alectoris graeca). Mol Ecol 12: 2201–2214.
- 8. Bonhomme M, Blancher A, Cuartero S, Chikhi L, Crouau-Roy B (2008) Origin and number of founders in an introduced insular primate: estimation from nuclear genetic data. Mol Ecol 17: 1009–1019.
- 9. Reich DE, Goldstein DB (1998) Genetic evidence for a Paleolithic human population expansion in Africa. Proc Nat Acad Sci USA 95: 8119–8123.
- 10. Liu H, Prugnolle F, Manica A, Balloux F (2006) A geographically explicit genetic model of worldwide human-settlement history. Am J Hum Genet 79: 230–237.
- 11. Kimura M, Ohta T (1978) Stepwise mutation model and distribution of allelic frequencies in a finite population. Proc Nat Acad Sci USA 75: 2868–2872.
- 12. Ellegren H (2004) Microsatellites: simple sequences with complex evolution. Nat Rev Genet 5: 435–445.
- 13. Di Rienzo A, Peterson AC, Garza JC, Valdes AM, Slatkin M, et al. (1994) Mutational processes of simple-sequence repeat loci in human populations. Proc Nat Acad Sci USA 91: 3166–3170.
- 14. Zhivotovsky LA, Feldman MW, Grishechkin SA (1997) Biased mutations and microsatellite variation. Mol Biol Evol 14: 926–933.
- 15. Estoup A, Jarne P, Cornuet J-M (2002) Homoplasy and mutation model at microsatellite loci and their consequences for population genetics analysis. Mol Ecol 11: 1591–1604.
- 16. Primmer CR, Ellegren H (1998) Patterns of molecular evolution in avian microsatellites. Mol Biol Evol 15: 997–1008.
- 17. Symonds VV, Lloyd AM (2003) An analysis of microsatellite loci in Arabidopsis thaliana: Mutational dynamics and application. Genetics 165: 1475–1488.
- 18. Curtu A-L, Finkeldey R, Gailing O (2004) Comparative sequencing of a microsatellite locus reveals size homoplasy within and between European oak species (Quercus spp.). Plant Mol Biol Rep 339–346: 22.
- 19. Stajner N, Jakse J, Kozjak P, Javornik B (2005) The isolation and characterisation of microsatellites in hop (Humulus lupulus L.). Plant Sci 168: 213–221.
- 20. Xie H, Sui Y, Chang F-Q, Xu Y, Ma R-C (2006) SSR allelic variation in almond (Prunus dulcis Mill.). Theor Appl Genet 112: 366–372.
- 21. Grimaldi M-C, Crouau-Roy B (1997) Microsatellite allelic homoplasy due to variable flanking sequences. J Mol Evol 44: 336–340.
- 22. Orti G, Pearse DE, Avise JC (1997) Phylogenetic assessment of length variation at a microsatellite locus. Proc Nat Acad Sci USA 94: 10745–10749.
- 23. Mogg R, Batley J, Hanley S, Edwards D, O’Sullivan H, et al. (2002) Characterization of the flanking regions of Zea mays microsatellites reveals a large number of useful sequence polymorphisms. Theor Appl Genet 105: 532–543.
- 24. Lia VV, Bracco M, Gottlieb AM, Poggio L, Confalonieri VA (2007) Complex mutational patterns and size homoplasy at maize microsatellite loci. Theor Appl Genet 115: 981–991.
- 25. Payseur BA, Cutter AD (2006) Intergrating patterns of polymorphism at SNPs and STRs. Trends Genet 22: 424–429.
- 26. Barkley NA, Roose ML, Krueger RR, Federici CT (2006) Assessing genetic diversity and population structure in a citrus germplasm collection utilizing simple sequence repeat markers (SSRs). Theor Appl Genet 112: 1519–1531.
- 27. Barkley NA, Krueger RR, Federici CT, Roose ML (2009) What phylogeny and gene genealogy analyses reveal about homoplasy in citrus microsatellite alleles. Plant Syst Evol 282: 71–86.
- 28. Gugerli F, Walser J-C, Douvani K, Holderegger R, Finkeldey R (2007) Coincidence of small-scale spatial discontinuities in leaf morphology and nuclear microsatellite variation of Quercus petraea and Q. robur in a mixed forest. Ann Bot 99: 713–722.
- 29. Gugerli F, Brodbeck S, Holderegger R (2008) Insertions–deletions in a microsatellite flanking region may be resolved by variation in stuttering patterns. Plant Mol Biol Rep 26: 255–262.
- 30. Kijas JMH, Thomas MR, Fowler JCS, Roose ML (1997) Integration of trinucleotide microsatellites into a linkage map of Citrus. Theor Appl Genet 94: 701–706.
- 31. Kampfer S, Lexer C, Glössl J, Steinkellner H (1998) Characterization of (GA)n microsatellite loci from Quercus robur. Hereditas 129: 183–186.
- 32. Billotte N, Lagoda PJL, Risterucci A-M, Baurens F-C (1999) Microsatellite-enriched libraries: applied methodology for the development of SSR markers in tropical crops. Fruits 54: 277–288.
- 33. Colpaert N, Cavers S, Bandou E, Caron H, Gheysen G, et al. (2005) Sampling tissue for DNA analysis of trees: Trunk cambium as an alternative to canopy leaves. Silvae Genet 54: 265–269.
- 34. Excoffier L, Lischer HEL (2010) Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Res 10: 564–567.
- 35.
R Development Core Team (2011) R: A language and environment for statistical computing. v. 2.14.1 R Foundation for Statistical Computing, Vienna, Austria. Accessed 2012 June 18.
- 36. Weir BS, Cockerham CC (1984) Estimating F-statistics for the analysis of population structure. Evolution 38: 1358–1370.
- 37. Michalakis Y, Excoffier L (1996) A generic estimation of population subdivision using distances between alleles with special reference for microsatellite loci. Genetics 142: 1061–1064.
- 38. Pons O, Petit RJ (1996) Measuring and testing genetic differentiation with ordered versus unordered alleles. Genetics 144: 1237–1245.
- 39. Excoffier L, Smouse PE, Quattro JM (1992) Analysis of molecular variance inferred from metric distances among DNA haplotypes: Application to human mitochondrial DNA restriction data. Genetics 131: 479–491.
- 40. Huson DH (1998) SplitsTree: analysing and visualizing evolutionary data. Bioinformatics 14: 68–73.
- 41. Huson DH, Bryant D (2006) Application of phylogenetic networks in evolutionary studies. Mol Bio Evol 23: 254–267.
- 42. Peakall R, Gilmore S, Keys W, Morgante M, Rafalski A (1998) Cross-species amplification of Soybean (Glycine max) simple sequence repeats (SSRs) within the genus and other legume genera: Implications for the transferability of SSRs in plants. Mol Biol Evol 15: 1275–1287.
- 43. Saha MC, Mian MAR, Eujayl I, Zwonitzer JC, Wang L, et al. (2004) Tall fescue EST-SSR markers with transferability across several grass species. Theor Appl Genet 109: 783–791.
- 44. McDonald MJ, Wang W-C, Huang H-D, Leu J-Y (2011) Clusters of nucleotide substitutions and insertion/deletion mutations are associated with repeat sequences. PLoS Biol 9: e1000622.
- 45. Brohede J, Ellegren H (1999) Microsatellite evolution: polarity of substitutions within repeats and neutrality of flanking sequences. Proc R Soc Lond B 266: 825–833.
- 46. Shepherd LD, Lambert DM (2005) Mutational bias in penguin microsatellite DNA. J Hered 96: 566–571.
- 47. Colson I, Goldstein DB (1999) Evidence for complex mutations at microsatellite loci in drosophila. Genetics 152: 617–627.
- 48. Cornuet J-M, Beaumont MA, Estoup A, Solignac M, Carey PD (2006) Inference on microsatellite mutation processes in the invasive mite, Varroa destructor, using reversible jump Markov chain Monte Carlo. Theor Popul Biol 69: 129–144.