BACKGROUND
genetic diversity is essential for organisms to adapt to changing environmental conditions and is recognised as a key component of biodiversity .
microsatellite markers  are a widely used marker system to estimate genetic diversity in population genetic studies and are often implicitly assumed to reflect the genome-wide diversity of a taxon  <cit> .
the use of microsatellites has increased linearly since their detection in the 1980s  <cit> , and they are nowadays extensively applied, for example in conservation genetics , forensic dna profiling, paternity analyses, and studies of neutral genetic population structure .
however, the challenge of correctly interpreting microsatellite data is often strongly underrated  <cit> , and the question whether a limited number of microsatellite markers accurately reflects genome-wide diversity remains a contentious issue .
single nucleotide polymorphisms  on the basis of traditional dna sequencing  <cit>  have long been known, but in contrast to microsatellites, were relatively rarely used in population genetics until recently because of the difficulties associated with their characterization and genotyping in non-model organisms  <cit> .
moreover, their  bi-allelic state limits the information content per locus compared to the more polymorphic microsatellite markers .
in recent years, the use of snps has been exponentially increasing  <cit> , mainly because newly developed high-throughput sequencing techniques can efficiently be applied to a wide range of organisms.
these techniques allow for the identification of thousands to millions of unbiased snps, and the simultaneous estimation of snp frequencies across the genomes of individuals, populations and species .
microsatellites have unique properties that distinguish them from the rest of the genome, and these should be taken into consideration when analysing and interpreting them  <cit> .
microsatellites are codominant markers and typically consist of simple sequence repeats varying in length between one and six base pairs.
their variability originates from dna polymerase slippage during replication, leading to the formation of shorter or longer alleles .
in plants, microsatellite mutation rates range between 10− <dig> and 10− <dig> per locus and generation , thus varying approximately  <dig> -fold, and are affected by various factors, including repeat type, repeat copy number, marker location in the genome, and taxon  <cit> .
in contrast, spontaneous mutation rates for snps only vary about 100-fold  <cit> .
knowledge of direct estimates of snp mutation rates is limited, but the rate has been accurately estimated e.g.
in arabidopsis thaliana to be 7 × 10− <dig> substitutions per site per generation  <cit> .
microsatellite mutation rates are therefore several orders of magnitude higher and much more variable than those of snps.
in combination with the often small number of markers used, microsatellite-based studies typically sample a narrow fraction of the genome with unusually high mutation rate  <cit> .
this may be aggravated when only the most polymorphic microsatellite markers are selected for further analysis after initial screening of a small subsample of individuals or populations.
estimates of genetic diversity may then suffer from ascertainment bias  <cit> .
additionally, amplification variation of primers  <cit>  and fragment size homoplasy  <cit>  potentially reduce the accuracy of genetic estimates inferred from microsatellite markers.
the use of microsatellite markers may thus lead to estimates of genetic diversity and differentiation that do not well reflect genome-wide patterns of variation.
despite these potential caveats, a large number of studies has relied on microsatellite markers to estimate genetic diversity and genetic differentiation, not only within and among populations, but also among species .
in a conservation context, microsatellites are also used to identify conservation units , whose genetic variation and distinctness is potentially relevant for species survival .
well-known case studies are the florida panther  <cit>  or the african elephant from eritrea  <cit> , for which management decisions were taken based on genetic data derived from few microsatellites.
to date, only few studies have explored in detail to what degree microsatellite variation reflects genetic variation at other nuclear loci, and which genetic diversity estimator for microsatellites provides the most accurate prediction of genome-wide diversity.
positive but sometimes weak correlations between expected microsatellite heterozygosity  and snp diversity in nuclear gene sequences have been reported at the population level in salmon  <cit>  and several carnivore species  <cit> , as well as different rice varieties and sheep breeds  <cit> .
most of these studies, however, have investigated only a limited number of snps .
the outcome of the comparison of snp versus microsatellite diversity in these studies was strongly affected by the number of snp markers used.
studies in which low snp numbers  were compared to microsatellites found that the latter had more power to infer differences in genetic summary statistics  or found similar results when approximately  <dig> snps were used  <cit> .
in contrast, studies using larger numbers of snps  found that snps performed better than microsatellites  <cit> .
many of these studies used existing genotyping arrays for snp detection.
these may, however, cause ascertainment bias as a consequence of the overrepresentation of common snps  <cit> .
to date, no unbiased whole-genome re-sequencing approach has been used for comparison.
studies based on reduced representation libraries , which sample a subset of all snps of the genome, showed that snps have more power than microsatellites, e.g.
to detect heterozygosity–fitness correlations in natural populations of oldfield mice  <cit> .
further, demographic inferences drawn from radseq-derived snps in bumble bees reflected important long-term differences in population size better than microsatellites, which instead signalled either recent demographic changes or mutational processes  <cit> .
because of the widespread application of microsatellite markers both in basic research and practical conservation, it is important to evaluate the tenet that microsatellite variation adequately reflects genome-wide genetic diversity, especially for situations in which only a limited number of markers are used, as is often the case in conservation genetics, where on average only  <dig> microsatellites are used per study  <cit> .
it is further relevant to evaluate the power of next-generation sequencing  based genotyping approaches to infer genome-wide diversity and population structure, e.g.
to estimate the number of snps required to achieve accurate and consistent estimates of genome-wide diversity.
we used two types of microsatellite markers  as well as genome-wide snp variation in the meadow rock cress, arabidopsis halleri , to compare estimates of genetic diversity and differentiation.
overall, we genotyped  <dig> individuals of a. halleri from nine natural populations using  <dig> microsatellite markers, which is above the average number of microsatellites typically used in population and conservation genetic studies  <cit> .
we compared them to a pooled whole-genome re-sequencing approach  and tested whether estimates of genetic variation derived from microsatellite polymorphisms are valid and useful proxies of genome-wide genetic variation and differentiation.
specifically, we tested whether estimates from both marker types were correlated  and had similar absolute values .
further, we used down-sampling to assess how many random and presumably unlinked snps are required to calculate accurate estimates of genome-wide diversity.
methods
study system
arabidopsis halleri is a perennial, insect-pollinated, strictly outcrossing and functionally self-incompatible herb  <cit>  with a wide geographic distribution from central europe to eastern asia  <cit> .
it grows in diverse habitats including mountain slopes, grassy meadows, forest margins and rocky crevices  <cit>  and has been widely used as a model to study heavy metal tolerance .
the species is diploid with 2n =  <dig>  <cit>  and has an estimated genome size of  <dig> mbp  <cit> .
sampling and dna extraction
leaf tissue from  <dig> individuals each of nine populations of a. halleri was sampled in south-eastern switzerland and northern italy .
the selected samples size per population should allow to accurately estimate population-specific genetic diversity and differentiation  <cit> .
a minimal distance of two, but preferably 4 m, was maintained between collected individuals, as genetic structure and diversity may be affected by clonal growth when plants are separated by less than one meter  <cit> .
leaf samples were dried in silica gel, and dna was extracted with the dneasy plant mini kit  according to the manufacturer’s protocol.
dna concentrations were measured with a qubit®  <dig>  fluorometer , and dna quality was examined using a nanodrop  <dig> spectrophotometer  as well as  <dig> % agarose gels stained with gelred .
microsatellite analyses
the  <dig> samples were genotyped using  <dig> microsatellite markers in three multiplex pcrs, each amplifying either six or eight microsatellite markers .
the first two multiplex sets included  <dig> microsatellite loci that were originally developed for a. thaliana  <cit> , hereafter referred to as “cross-species” microsatellites.
some of these primer sequences were adapted to a. halleri by comparing them to our own de-novo assembly of the a. halleri genome  <cit>  using igv  <dig>   <cit> , identifying potential mismatches and changing the primer sequences accordingly .
further, eight microsatellite primer pairs that were specifically developed for a. halleri  <cit>  were combined in a third multiplex set.
these markers are hereafter referred to as “species-specific” microsatellites.
detailed lab protocols can be found in the additional methods .
alleles were called using genemapper  <dig>  .
estimating microsatellite-based genetic diversity
for every marker and population, we assessed the following population genetic parameters.
the inbreeding coefficient f
is and its p-value, which indicate whether markers or populations deviate from hardy–weinberg equilibrium, were calculated with genodive  <dig> b <dig>  <cit>  using the heterozygosity-based g
is statistic with  <dig> permutations and applying bonferroni correction for multiple testing.
null allele frequencies were calculated with freena  <cit> .
pairwise values of genetic differentiation among populations, f
st, were calculated based on allele identity with genepop  <dig> .
<dig>  <cit> , whereas allele frequencies, expected heterozygosity  and mean number of alleles  per locus were quantified with genetix  <dig>   <cit> .
we consistently genotyped  <dig> individuals per population, therefore, a
r did not have to be corrected with a rarefaction approach.
all population parameters were computed for three different marker sets including  all,  the cross-species, and  the species-specific microsatellite markers .
to infer marker bias, we tested for quantitative differences in estimates of ssr-h
e and a
r estimated in each population from cross-species and species-specific markers using a paired t-test  in r  <dig>   <cit> .
for the relative comparison of ssr-h
e derived from cross-species and species-specific markers, we used a pearson’s correlation test  in r. to test whether population-specific estimates of ssr-h
e obtained from the different microsatellite types differ  <cit> , we used a pairwise wilcoxon signed-rank test  in r. by plotting ssr-h
e medians and quantiles for each population, we inferred whether non-significant differences were caused by high variance.
when the absence of significant differences between populations was obviously caused by overly high variances in the genetic diversity estimates computed per microsatellite marker and population, we interpreted this as variance bias.
p-values were adjusted for multiple testing using bonferroni correction.
pool-seq and illumina read processing
pooled next-generation sequencing  has been shown to produce accurate population-specific allele frequencies  <cit> .
for ngs, individually extracted high-quality dna was equimolarly pooled using the same  <dig> individuals from the nine populations as for the microsatellite genotyping presented above.
these nine pools were high-throughput sequenced and further processed as described below and, in more detail, in fischer et al.
<cit> .
for a subset of snps and populations of the present dataset, the accuracy of exactly the same pool-seq approach had been validated  <cit> : differences in estimates of population-specific allele frequencies compared to those from individual genotyping were on average less than 4%.
library preparation  and sequencing on an illumina hiseq <dig>  were performed by gatc biotech  and the quantitative genomics facility .
forward and reverse raw reads were trimmed for tags and adaptors with cutadapt  <cit> .
phred-type quality scores q <dig> were used for quality trimming with the fastx toolkit .
the separately trimmed forward and reverse reads were then re-synchronized to pairs with an in-house perl script.
only paired sequences were used for further analysis  <cit> .
read mapping, snp calling and genome-wide population genetic estimates
to estimate genome-wide genetic diversity and differentiation for all nine populations of a. halleri, reads were mapped to the a. thaliana reference genome  using bwa aln, allowing for 10% mismatch, and sampe  <cit> .
all ambiguously mapped reads were removed and the remaining high-quality reads were sorted with samtools  <dig> .
<dig>  <cit> .
snps were called for the nine populations by producing mpileup files with samtools .
to obtain population-specific genome-wide estimates of genetic diversity, we first calculated watterson’s theta , an estimator that takes into account the number of segregating sites to estimate the population mutation rate.
θ
watterson was calculated on a gene-by-gene basis  for each population using the gene and exon annotation of tair <dig> .
the perl script ‘variance-at-position.pl’ of the software package popoolation  <cit>  was used with the mpileup file of each population and the.gtf annotation file  to calculate exonic θ
watterson.
this approach provides unbiased estimates for pooled samples as it corrects for coverage.
exon based θ
watterson is a conservative estimate of genetic diversity as it infers diversity from genomic regions that are predominantly under purifying selection, hence from slightly less diverse regions than the rest of the genome.
exon-based diversity estimates are of direct relevance for the adaptive potential of a population, because exons harbour functionally relevant polymorphisms that allow populations to adapt to changing environments.
to accurately estimate allele frequencies for estimating θ
watterson, minimum counts for minor alleles were set to two to account for sequencing errors, leading to a minor allele frequency threshold of  <dig> .
the minimum coverage per site within populations was set to 20×, which mimics the number of individuals.
to further correct for potential errors caused by repeated sequences, a maximum coverage of 400× per population was used as threshold for snp identification.
in order to be included in the genome-wide estimates of gene diversity, 50% of all snps within a gene had to reach the above-mentioned thresholds in all nine populations  <cit> .
for all analyses, pool size per population was set to  <dig> because  <dig> diploid genomes were represented in each population pool.
second, we calculated genome-wide snp-based expected heterozygosity , taking all snps into account, not only those located in exons.
mpileup files were synchronized and filtered for base quality  with the perl script ‘mpileup2sync.pl’ of popoolation <dig>  <cit> .
next, major and minor allele frequencies were calculated with the script snp-frequency-diff.pl.
the coverage threshold was the same as mentioned above, except that the minor allele count was set to four, as all nine populations were jointly used to infer minor allele frequencies  <cit> , leading to a more sensitive, but less error-prone minor allele frequency threshold of  <dig> .
we only used bi-allelic snps and calculated the average genome-wide snp-h
e as snp-he=1n∑i=1n2pi where n is the number of snps, and p
i is the minor allele frequency of the ith allele.
this approach assumes hardy–weinberg equilibrium within populations.
to infer potential demographic events that could strongly influence genetic diversity within populations, we calculated exome-wide tajima’s d using the tair <dig> gene annotation and the perl script ‘variance-at-position.pl’ in popoolation  <cit> .
it is suggested to use a coverage threshold of less than three times smaller than the pool size  <cit> , which is in our case 13×.
a negative genome-wide tajima’s d is indicative of an expansion after a bottleneck, whereas a positive d is compatible with a scenario of a decrease in population size  <cit> .
to test whether the average of the resulting distribution of tajima’s d was significantly different from zero, we used t-tests against random normal distributions  with an average of zero and the same standard deviation as observed in the real data of each population.
estimates of pairwise population genetic differentiation  were calculated with ‘fst-sliding.pl’ in popoolation <dig>  <cit> .
average values of pairwise f
st were calculated using the same parameters as mentioned above for the estimates of snp-h
e as explained in detail in fischer et al.
<cit> .
comparisons of genetic diversity and differentiation derived from microsatellites and genome-wide snps
to explore associations between estimates of genetic diversity derived from microsatellites and snps, we performed pearson’s correlations  of microsatellite-based allelic richness  and expected heterozygosity .
estimates of genome-wide snp diversity were derived from exon sequences  and genome-wide snp expected heterozygosity .
further, to account for possible confounding effects due to linkage, associations were also tested for a subset of snps, consisting of every 50th snp.
we performed mantel tests to check for correlations between values of f
st derived from genome-wide snps and  all,  cross-species,  and species-specific microsatellite markers.
the same analysis was used to assess correlations between values of f
st derived from species-specific and cross-species microsatellite markers.
all analyses were performed with  <dig> permutations using ecodist  <dig> .
<dig>  <cit>  in r. finally, we used paired t-tests implemented in r to quantitatively evaluate whether h
e and f
st derived from microsatellite markers and snps significantly differ.
estimating the number of unbiased snps required for accurate estimates of genetic diversity
we used a down-sampling procedure to estimate snp-h
e with the aim to infer the required number of randomly selected and unlinked snp markers to obtain accurate estimates of genetic diversity and to reliably rank populations according to their genetic diversity .
thus, each population was resampled for the same k random snp markers drawn from the pool of more than  <dig> million snps.
for each value of k varying between  <dig> and  <dig> , we created  <dig> random subsamples of k snp markers .
we then computed the mean expected heterozygosity  and 95% confidence intervals for each value of k observed in each of the  <dig> subsamples.
obtained results were used to draw curves representing the variation of the estimated h
e as a function of genotyping effort in each population.
we then identified the number of snps for which the upper and lower confidence intervals for expected heterozygosity  fell below ± <dig> , ± <dig> , and ± <dig> .
RESULTS
microsatellite diversity
twenty microsatellite loci were initially used to characterise  <dig> a. halleri individuals from nine populations.
we excluded marker ah <dig> from further analyses, because it deviated significantly from hardy–weinberg equilibrium and exhibited an estimated null allele frequency of 10% .
the remaining  <dig> microsatellite markers harboured  <dig> alleles and only  <dig> % missing data.
allelic richness per locus  ranged between  <dig>  and  <dig>  per population, with an average of  <dig>  , and expected heterozygosity  ranged from  <dig>  to  <dig>  per microsatellite marker.
further details are given in table  <dig> and additional file 4: table s <dig>  population allele frequency distributions were fairly noisy, see additional file 5: figure s1a.
none of the nine populations showed significant deviation from hardy–weinberg equilibrium after bonferroni correction .table  <dig> population genetic parameters inferred from  <dig> microsatellites and genome-wide snps for nine populations of arabidopsis halleri.
allelic richness , expected heterozygosity  and inbreeding coefficient f
is including its one-sided p-value  are given.
no f
is value was significantly different from zero after bonferroni correction.
θ
watterson was calculated for  <dig>  genes.
expected heterozygosity  was calculated from all snps across the genome.
tajima’s d was calculated for  <dig>  genes, p-values refer to deviations from zero
a
r
e
f
is
f
is
θ
watterson
e
all twelve cross-species microsatellite markers, initially developed for a. thaliana, successfully amplified in a. halleri.
a
r at cross-species microsatellite markers was  <dig>  alleles per marker and population and thus significantly lower than a
r at species-specific markers .
the same pattern  was observed for ssr-h
e, which was  <dig>  and  <dig> , respectively.
no significant correlation was observed among estimates of h
e inferred from cross-species and species-specific microsatellite markers .fig.
<dig> comparison of estimates of genetic diversity derived from cross-species  and species-specific microsatellite markers  for a allelic richness  and b expected microsatellite heterozygosity .
c estimates of ssr-h
e inferred separately from cross-species and species-specific microsatellite markers were not significantly correlated .
dots are labelled with population codes .
d no significantly different estimates of h
e were observed among populations after bonferroni correction .
without correction for multiple testing only population aha <dig> and ahan <dig> showed significantly different estimates of h
e
the variance in the estimates of ssr-h
e among microsatellite markers within and among population was so high that no significant differences in genetic diversity among populations could be inferred after bonferroni correction .
without bonferroni correction, only populations aha <dig> and ahan <dig> differed significantly in their estimates of h
e .
similar results were found when cross-species and species-specific microsatellites were analysed separately and corrections for multiple testing were performed .
average pairwise f
st for the  <dig> microsatellite loci was  <dig>  ; species-specific  and cross-species microsatellites markers  did not deviate significantly from each other .
illumina sequencing and genome-wide diversity
the illumina sequencing yielded a total of  <dig> , <dig>  100-bp paired-end reads corresponding to  <dig> , <dig>  nucleotides.
after quality filtering and trimming,  <dig> , <dig>  paired-end reads were mapped to the a. thaliana reference genome , from which organellar dna was excluded.
the average coverage per site, after filtering with the defined thresholds, was  <dig> × with a range of population-wise coverage of  <dig>  to  <dig>  × .
we detected  <dig> , <dig> snps, which were used for the calculation of pairwise f
st, and  <dig> , <dig> bi-allelic snps, which were used for the calculation of population-specific snp-h
e. all populations had even population-specific allele frequency distributions .
values of pairwise f
st ranged between  <dig>  and  <dig> , and population-specific snp-h
e was between  <dig>  and  <dig> .
overall,  <dig>  and  <dig>  genes fulfilled our thresholds of coverage for calculating θ
watterson and tajima’s d, respectively.
θ
watterson ranged from  <dig>  to  <dig> , and tajima’s d values were all slightly negative, ranging from − <dig>  to − <dig> .
only two of the nine populations showed no significant deviation from zero.
in other words, most populations showed weak demographic changes probably related to a bottleneck with later expansion .
comparisons of genetic diversity estimates derived from microsatellites and genome-wide snps
no significant correlation was observed between population-specific estimates of h
e derived from microsatellites and genome-wide snps , independent of whether all snps or a subset of presumably unlinked snps  were used .
estimates of snp-h
e were overall significantly lower than those of ssr-h
e .
ssr-h
e was also not significantly correlated with θ
watterson .
the correlation coefficient was higher when using only species-specific markers , but lower when only cross-species markers were used , though neither of the two correlations was significant.fig.
<dig> relationships between population genetic parameters for nine populations of arabidopsis halleri estimated on the basis of  <dig> microsatellites and  <dig> million genome-wide single nucleotide polymorphisms .
a relationship between expected genome-wide snp heterozygosity  and expected microsatellite heterozygosity .
the insert illustrates differences in h
e between snps and microsatellites across the nine a. halleri populations.
b relationship between genome-wide diversity estimated by average θ
watterson across  <dig>  genes and expected heterozygosity for all  <dig> microsatellites  and separately for species-specific  and cross-species microsatellites .
c comparisons of allelic richness  with genome-wide snp variation estimated by θ
watterson
in marked contrast to heterozygosity, a
r of microsatellite markers was significantly correlated with θ
watterson .
the correlation based on ranks was still significant , but slightly weaker than non-ranked comparisons.
the snp-based genome-wide diversity estimates, θ
watterson and snp-h
e, were highly correlated , even though values of θ
watterson were derived exclusively from coding regions , and estimates of snp-h
e were calculated from more than two million snps across the entire genome.
in fact, θ
watterson estimates inferred from introns and intergenic regions were highly correlated to estimates of exon-based θ
watterson .fig.
<dig> pearson’s correlation  between exome-wide θ
watterson and genome-wide expected snp heterozygosity  in arabidopsis halleri
comparison of genetic differentiation estimates derived from microsatellite versus genome-wide snp variation
mantel tests revealed highly significant correlations between values of pairwise f
st derived from genome-wide snp data and microsatellite markers .
the best correlation was achieved when using all  <dig> microsatellites .
however, values of f
st derived from microsatellite markers were  <dig> -fold higher than those from snps and were significantly different .
if we split the microsatellite markers into species-specific and cross-species, the correlations were slightly weaker for cross-species microsatellites  and for species-specific microsatellites .
the correlation among species-specific and cross-species microsatellites was high .fig.
<dig> comparison of pairwise population genetic differentiation  among nine populations of arabidopsis halleri estimated with different genetic markers.
a comparison of f
st derived from  <dig> , <dig> genome-wide snps with f
st inferred on the basis of all  <dig> microsatellite markers , and separately for  <dig> cross-species , and seven species-specific microsatellites .
given are correlation coefficients  and p-values from mantel tests.
b estimates of f
st inferred from snps were significantly lower than those based on microsatellites
estimating the number of unbiased snps required for accurate estimates of genetic diversity
snp down-sampling revealed that the upper and lower confidence intervals for expected heterozygosity  fell below ± <dig> , ± <dig> , and ± <dig>  with  <dig>   <dig> and  <dig>  random snps, respectively.
to accurately and consistently rank all a. halleri populations according to their genome-wide diversity ,  <dig>  snps were required .
however, populations aha <dig> and ahan <dig> could not be distinguished, as they have the same genome-wide snp-h
e .fig.
<dig> the effect of the number of randomly selected snp markers on the estimation of h
e represented for each of nine populations of arabidopsis halleri.
light colours indicate the 95% confidence intervals as inferred from  <dig> randomly drawn replicates.
dashed horizontal lines  indicate the numbers of snps for which 95% confidence intervals for snp-h
e are below ± <dig> , ± <dig> , and ± <dig>
discussion
the application of microsatellite markers is widespread in population and conservation genetic studies.
however, ngs-based snp genotyping approaches are rapidly developing and can be applied to a wide diversity of model and non-model organisms.
our comparative analysis of genetic diversity estimates based on microsatellites and genome-wide snps revealed interesting differences.
most importantly, we found no significant correlation between expected microsatellite heterozygosity , an estimator of genetic diversity that is widely used and reported in microsatellite studies, and genome-wide snp diversity .
this finding indicates that ssr-h
e does not adequately reflect genome-wide genetic diversity in the investigated populations of a. halleri.
in contrast, microsatellite allelic richness  was a much better proxy for genome-wide diversity .
further, genetic differentiation in terms of f
st estimated from microsatellite variation correlated reasonably well with that based on genome-wide snp data.
however, absolute values of the different summary statistics inferred from different marker types varied considerably .
our results do not question the usefulness of microsatellites per se, but point to research questions for which snps may be better suited than microsatellite markers, given the availability of robust and cost-effective high-throughput sequencing-based snp genotyping approaches.
thanks to the massive advances in sequencing technology, thousands of snps can efficiently be genotyped in any given organism, and these may improve our ability to adequately estimate genetic diversity and differentiation.
previous studies found that the potential of snps to resolve population genetic structure strongly depends on their number.
indeed, in studies using a low number of assay-based snps, microsatellites performed similarly well or better than snps .
however, in studies using larger snp numbers, especially when they were derived from ngs-based approaches, the relative performance of snps clearly improved  <cit> .
these studies show that a large number of snps compensates for the lower information content of these typically bi-allelic markers compared to more polymorphic microsatellite markers.
our study shows that a few thousand snps are enough to accurately estimate genome-wide diversity in terms of h
e .
evidence available from animal species  <cit>  as well as from this study reveals that there is no or only weak congruence between estimates of heterozygosity derived from microsatellites and snps.
indeed, theory suggests that an association of heterozygosity estimates between microsatellites and genome-wide snps is not expected a priori.
according to ljungqvist et al.
<cit> , this association is shaped by identity disequilibrium, i.e.
the non-random association of diploid genotypes between loci.
simulations indicated that a strong positive correlation only emerges when the studied populations are characterized by substantial identity disequilibrium, as was the case in the studies on salmon and carnivores  <cit> .
the absence of this correlation in a. halleri suggests that identity disequilibrium is weak or absent, which is compatible with the existence of a strong self-incompatibility system in this species and observed non-significant values of f
is.
this finding indicates that heterozygosity estimates based on microsatellites  might not be good surrogates for genome-wide diversity in outbreeding species.
further, we found that the marker-specific variation of ssr-h
e is too high to distinguish populations based on their levels of genetic diversity  and that even among microsatellite markers originally developed for different taxa, important discrepancies can be observed .
this variance bias is especially strong when a low number of microsatellites is used, as is the case in many population and conservation genetic studies .
the larger sampling variance associated with a limited number of microsatellite markers is evidenced by their allele frequency distributions, which are much noisier than those derived from the snp data set .
an alternative estimator to h
e is microsatellite allelic richness , one of the simplest estimators of genetic diversity available.
a
r was significantly correlated with genome-wide snp diversity in our study  and thus appears to be a useful proxy of genome-wide genetic diversity.
congruent results of a
r and snp diversity were also reported in other studies  <cit>  and might be explained by several reasons.
first, ssr-h
e estimates are based on few markers with noisy allele frequency distributions  and represent a proportion, ranging between  <dig> and  <dig>  whereas allelic richness is an infinite count.
accordingly, one additional allele, especially when it is rare and many alleles are already present, does not strongly influence ssr-h
e, but will affect a
r estimates, making the latter a more sensitive estimator of diversity.
moreover, especially for microsatellites with a high number of alleles, accurate estimates of population-specific heterozygosity can be problematic  <cit> , and stochasticity may have a strong impact on estimates of h
e at the lower range of allelic diversity.
finally, a
r is more sensitive to population bottlenecks than h
e  <cit> .
therefore, a
r better reflects the population’s demographic history and hence is a more relevant estimator of genetic diversity to predict the short-term survival of a population.
despite the good performance of a
r as a proxy of genome-wide diversity, it was not sufficient to accurately rank populations according to their genetic diversity .
consequently, the identification of conservation units  or decision taking for conservation actions based on microsatellite-derived rankings of genetic diversity may be misguiding.
nevertheless, using microsatellite-derived a
r rather than ssr-h
e provides more accurate estimates of genome-wide genetic diversity derived from a limited number of microsatellite markers.
these considerations are strengthened by a simulation study  <cit> , which found allelic richness to be two to four times more powerful than h
e for the identification of a temporal genetic decline in a population.
a possible reason for the deviation of estimates of genetic diversity derived from microsatellite and genome-wide snp data could be the influence of very recent demographic changes.
as a consequence of their higher mutation rate  <cit> , microsatellites respond more strongly to recent demographic events than genome-wide snps  <cit> , and snps uncover a different and likely older demographic history.
evidence for this hypothesis has also been presented for bumble bees  <cit>  and may indicate a fruitful application of microsatellites in analyses of populations that may have undergone very recent demographic changes.
however, it is important to note that loci with high mutation rate  may violate demographic model assumptions, such as mutation–migration–drift equilibrium  <cit> .
further, for the long-term survival of populations, genetic diversity in coding and regulatory sequences is arguably more relevant than microsatellite diversity, because most microsatellites are located in non-coding regions and are mostly selectively neutral, hence of less evolutionary importance.
consequently, estimates of genome-wide snp diversity better reflect functionally important and potentially adaptive genetic variation  <cit> , and should therefore be used preferentially, especially in conservation genetics studies.
in a. halleri the long-term demographic history inferred from genome-wide data indicates a bottleneck with later expansion for most populations, as values of tajima’s d were slightly negative and significantly different from zero .
the origin of the microsatellite markers used, i.e.
whether they are species-specific or cross-species markers, may further impact estimates of genetic diversity  <cit> .
in this study, species-specific microsatellites displayed significantly higher a
r and ssr-h
e than cross-species markers  originally developed for a. thaliana  <cit> .
further, species-specific markers resulted in more accurate estimates of genetic diversity , but less accurate estimates of divergence  among populations of a. halleri.
hence, the practicability of microsatellites for population genetic studies is limited and difficult to assess a priori  <cit> .
high estimates of genetic diversity derived from species-specific microsatellite markers may be a consequence of ascertainment bias caused by selecting the most polymorphic markers  <cit> , whereas cross-species microsatellites are mostly chosen based on their amplification success in the study species.
the consequence of this ascertainment bias is evident in this study, as estimates of ssr-h
e for cross-species microsatellites are not significantly correlated with ssr-h
e for species-specific microsatellites .
these differences further emphasize our inference that marker choice may substantially bias estimates of genetic diversity and may invalidate comparisons between populations or species, most notably when different microsatellite loci are assessed  <cit> .
a different pattern emerges for estimates of population genetic differentiation in terms of pairwise f
st. we found a significant positive correlation between values of pairwise f
st derived from microsatellites and genome-wide snps .
similar findings were reported for salmons and threespine sticklebacks  <cit> .
a reason for the better correlation between estimates of genetic differentiation compared to estimates of genetic diversity may be that more values are involved in pairwise comparisons, and that differences in allele frequencies of the common alleles are more important for the accurate estimation of genetic differentiation than those of rare alleles.
importantly, estimates of f
st derived from microsatellites were consistently and substantially higher than those based on genome-wide snps.
this seems counterintuitive, because multi-allelic microsatellite markers with high mutation rates  should cause lower f
st values than low-diversity markers like snps  <cit> .
however, pooled whole-genome re-sequencing studies with high coverage, such as this one, also detect rare variants; these low-frequency snps reduce overall f
st  <cit> .
overall, we consider pairwise population genetic differentiation estimated from microsatellites a useful proxy for genome-wide differentiation, but only in relative and not in absolute terms .
this finding has serious implications, because absolute values of f
st continue to be frequently used to infer indirect estimates of gene flow and migration, even though estimates of gene-flow should not be derived from f
st  <cit> .
further, this marker-specific difference in f
st estimates has a major impact on comparative studies of the divergence of quantitative traits, known as q
st–f
st comparisons, because the inference of the role of natural selection and genetic drift as causes of population genetic differentiation in complex polygenic traits is biased  <cit> .
while our results suggest that microsatellites should not be used for estimating genome-wide heterozygosity, we emphasize that microsatellites remain useful molecular markers for other applications.
for example, microsatellites perform very well in genetic stock identification or paternity analysis owing to their high variability  and may therefore continue to play an important role in molecular ecology.
however, before embarking on a molecular analysis, it remains a key issue to carefully assess the inherent strengths and limitations associated with different molecular markers  <cit> .
only then it is possible to select the most appropriate method for a given ecological or evolutionary question  <cit> .
in contrast to microsatellite-derived data, estimates of genome-wide diversity inferred from whole-genome re-sequencing data, e.g.
exome-wide θ
watterson, intronic and intergenic θ
watterson, or genome-wide snp-h
e, were highly correlated with each other and led to the same ranking of populations .
even though values of θ
watterson were either derived exclusively from coding regions or intronic and intergenic regions, and snp-h
e was calculated from positions across the whole genome, their estimates were highly congruent.
the slight variation observed among θ
watterson and snp-h
e  might be explained by differences in the demographic history among populations , because the demographic history of a population has a stronger influence on θ
watterson  than on snp-h
e or snp nucleotide diversity estimates  <cit> , as rare alleles are more likely to be lost during a bottleneck than common ones.
further, we found in our whole-genome re-sequencing study that the confounding effects of genetic linkage are negligible in the highly outcrossing and self-incompatible a. halleri, most likely because the small-scale linkage effects are compensated by the large numbers of unlinked snps.
thus, our snp-h
e estimates based on a subset of putatively unlinked snps were nearly identical to the estimates inferred for all snps, see fig.
2a and additional file 8: figure s <dig> .
similar to our genome re-sequencing study, approaches that use reduced representation libraries  to sample a subset of genome-wide snps can accurately estimate genome-wide heterozygosity.
as a consequence of the much smaller proportion of the genome surveyed with such approaches, however, care should be taken to avoid confounding effects of linkage, for example by considering only one snp per rad-locus  <cit> .
for example, the inbreeding coefficient of a known pedigree in oldfield mice showed strong concordance with the inferred estimates of heterozygosity obtained from  <dig>  radseq snps  <cit> .
this result indicates that, as long as a sufficiently large number of unbiased ngs-based snps is analysed across the genome, snp estimates accurately reflect genome-wide diversity in natural populations.
therefore, approaches like radseq  <cit> , pool-seq  <cit>  and whole-genome re-sequencing at low coverage  <cit>  are more appropriate than array-based snp approaches, which may be affected by strong ascertainment bias  <cit> .
an important question to consider in many studies may be the number of snps that are needed to estimate genetic diversity.
our down-sampling approach indicated that the number of random snps that are required to resolve genetic diversity difference among populations range from  <dig>  to  <dig>  snps .
this number is in the range of snps that can be inferred with standard radseq protocols also in non-model organisms  <cit> .
however, to differentiate among populations with very similar levels of genetic diversity, we required approximately  <dig>  snps .
thus, large snp datasets that are ideally identified de novo through ngs approaches  are highly suitable to distinguish, for example, between populations differing in genetic diversity and may therefore support decision-making in conservation management.
a further advantage of genome-wide snp data is that they not only allow one to estimate neutral genetic diversity, but also to identify adaptive genetic variation , which is considered essential for delimitating conservation units .
the large technical advances in nucleotide sequencing technology in recent years have not only massively increased the number of nucleotides that can be sequenced per individual or population, but have also led to reduced costs per nucleotide to the extent that screening a handful of microsatellite markers may be as expensive as surveying thousands of snps using latest ngs-based genotyping technologies .
CONCLUSIONS
this case study in the perennial and outcrossing plant a. halleri reveals that genetic diversity estimated from microsatellite markers, notably expected heterozygosity, may not adequately reflect genome-wide genetic diversity estimated from single-nucleotide polymorphisms and may therefore be a poor proxy for genome-wide estimates of genetic diversity.
possible causes include the limited number of microsatellite markers used, marker ascertainment bias, as well as the high variance in microsatellite-derived diversity estimates.
interestingly, microsatellite allelic richness  was found to be a reasonable proxy for genome-wide diversity, but the absolute ranking of populations was still inconsistent.
estimates of genetic differentiation  among populations derived from microsatellites were consistently higher than snp-based estimates but were significantly correlated with the latter.
our results do not question the usefulness of microsatellites per se, but point to research questions for which ngs-derived snps may be better suited than microsatellite markers, given the availability of robust and cost-effective snp genotyping approaches based on high-throughput sequencing.
as a consequence, we recommend using genome-wide analyses of snp diversity when the inference and comparison of genetic diversity within and among populations and species is the goal of a study.
a few thousand ngs-derived snps are sufficient for this purpose and this number of unbiased snps can nowadays easily been obtained also for non-model species, for example by using a reduced representation sequencing approach such as radseq  <cit> .
additional files
additional file 1: table s <dig>  sampling locations of the nine study populations of arabidopsis halleri.
additional file 2: table s <dig>  characteristics of  <dig> microsatellite markers used for arabidopsis halleri.
source refers to the species for which a microsatellite marker was originally developed.
sequences of forward and reverse primers are shown in 5′ - 3′ direction.
the fluorescent dyes used in multiplex pcrs and allele size ranges in base pairs are given for each locus.
additional file 3: additional methods.
additional file 4: table s <dig>  locus-specific summary statistics for all  <dig> microsatellite markers genotyped in arabidopsis halleri.
estimated frequency of null alleles, total number of alleles found among the  <dig> individuals, expected heterozygosity , and inbreeding coefficient f
is with one-sided p-value  are given.
additional file 5: figure s <dig>   population specific allele frequency distributions in nine populations of a. halleri for  <dig> microsatellite markers .
minor allele frequency distributions in the same nine populations across  <dig> , <dig> snps .
histograms are labelled with population codes .
additional file 6: figure s <dig>  box-whisker plot of expected heterozygosity calculated for each microsatellite marker  genotyped in arabidopsis halleri.
for cross-species microsatellites , no significantly different estimates of h
e  were observed after bonferroni correction, as a consequence of the high variance.
without correction for multiple testing, only population aha <dig> had a significantly different h
e estimates from aha <dig> and ahan <dig> .
for species-specific markers , no significant differences in h
e estimates were observed after bonferroni correction.
without multiple corrections, only aha <dig> was significantly different from ahan <dig>
additional file 7: figure s <dig>  observed  versus expected  distributions of tajima’s d for nine populations of arabidopsis halleri.
the expected normal distribution consists of  <dig>  values with an average of zero and the same standard deviation as the real dataset.
p-values indicate whether there was a significant deviation of the average from zero using a t-test.
additional file 8: figure s <dig>  relationships between expected snp heterozygosity  and expected microsatellite heterozygosity  of arabidopsis halleri estimated on the basis of  <dig> microsatellites and  <dig>  genetically unlinked genome-wide single nucleotide polymorphisms .
for this analysis only every 50th snp of the more than  <dig> million snps was used.
the median distance among two neighbouring snps is 1600 bp on the a. thaliana reference genome.
additional file 9: figure s <dig>  pearson’s correlation  between exome-wide θ
watterson and θ
watterson calculated from introns and intergenic regions in arabidopsis halleri.
additional file 10: figure s <dig>  comparison of pairwise population genetic differentiation  among nine populations of arabidopsis halleri based on  <dig> cross-species microsatellite markers and seven species-specific microsatellite markers.
r
mt represents the correlation coefficient of the mantel test and p the significance of the correlation.
abbreviations
armicrosatellite allelic richness
heexpected heterozygosity
ngsnext-generation sequencing
pool-seqpooled next-generation sequencing
radseqrestriction-site associated dna sequencing
snpsingle nucleotide polymorphism
ssrsimple sequence repeat
θwattersonwatterson’s theta
we thank c. beisel for supporting our illumina sequencing efforts at the quantitative genomics facility, d-bsse, eth zürich.
we also thank a. bösch, e. schnyder, j. zimmermann, d. zulliger, n. quebre, and n. aellen for sampling a. halleri populations, and c. michel, r. graf, and the genetic diversity centre zürich  for experimental support, a. tedder for helpful input, as well as s. fior and n. zemp for bioinformatics support.
we acknowledge the valuable comments made by anonymous reviewers on an earlier version of this article.
funding
this study was funded by the swiss national science foundation  and the adaptation to a changing environment  center of eth zurich.
availability of data and materials
the datasets supporting the conclusions of this article are available in the dryad digital repository, , in the european nucleotide archive  and in the sequence read archive .
authors’ contributions
mcf, cr, and aw designed the study; fg, mcf, and aw conducted field-work; mcf, cr, ml, and mr performed the research and analyzed the data; mcf, cr, and aw wrote the manuscript with substantial contributions from ml, mr, fg, and rh.
aw, kks and rh conceived the project and obtained funding.
all authors read and approved the final manuscript.
competing interests
the authors declare that they have no competing interests.
consent for publication
not applicable.
ethics approval and consent to participate
not applicable.

 genetic diversity is essential for organisms to adapt to changing environmental conditions and is recognised as a key component of biodiversity. microsatellite markers are a widely used marker system to estimate genetic diversity in population genetic studies and are often implicitly assumed to reflect the genome-wide diversity of a taxon. the use of microsatellites has increased linearly since their detection in the 1980s, and they are nowadays extensively applied, for example in conservation genetics, forensic dna profiling, paternity analyses, and studies of neutral genetic population structure. however, the challenge of correctly interpreting microsatellite data is often strongly underrated, and the question whether a limited number of microsatellite markers accurately reflects genome-wide diversity remains a contentious issue. single nucleotide polymorphisms on the basis of traditional dna sequencing have long been known, but in contrast to microsatellites, were relatively rarely used in population genetics until recently because of the difficulties associated with their characterization and genotyping in non-model organisms. moreover, their bi-allelic state limits the information content per locus compared to the more polymorphic microsatellite markers. in recent years, the use of snps has been exponentially increasing, mainly because newly developed high-throughput sequencing techniques can efficiently be applied to a wide range of organisms. these techniques allow for the identification of thousands to millions of unbiased snps, and the simultaneous estimation of snp frequencies across the genomes of individuals, populations and species. microsatellites have unique properties that distinguish them from the rest of the genome, and these should be taken into consideration when analysing and interpreting them. microsatellites are codominant markers and typically consist of simple sequence repeats varying in length between one and six base pairs. their variability originates from dna polymerase slippage during replication, leading to the formation of shorter or longer alleles. in plants, microsatellite mutation rates range between 10− and 10− per locus and generation, thus varying approximately -fold, and are affected by various factors, including repeat type, repeat copy number, marker location in the genome, and taxon. in contrast, spontaneous mutation rates for snps only vary about 100-fold. knowledge of direct estimates of snp mutation rates is limited, but the rate has been accurately estimated e.g. in arabidopsis thaliana to be 7 × 10− substitutions per site per generation. microsatellite mutation rates are therefore several orders of magnitude higher and much more variable than those of snps. in combination with the often small number of markers used, microsatellite-based studies typically sample a narrow fraction of the genome with unusually high mutation rate. this may be aggravated when only the most polymorphic microsatellite markers are selected for further analysis after initial screening of a small subsample of individuals or populations. estimates of genetic diversity may then suffer from ascertainment bias. additionally, amplification variation of primers and fragment size homoplasy potentially reduce the accuracy of genetic estimates inferred from microsatellite markers. the use of microsatellite markers may thus lead to estimates of genetic diversity and differentiation that do not well reflect genome-wide patterns of variation. despite these potential caveats, a large number of studies has relied on microsatellite markers to estimate genetic diversity and genetic differentiation, not only within and among populations, but also among species. in a conservation context, microsatellites are also used to identify conservation units, whose genetic variation and distinctness is potentially relevant for species survival. well-known case studies are the florida panther or the african elephant from eritrea, for which management decisions were taken based on genetic data derived from few microsatellites. to date, only few studies have explored in detail to what degree microsatellite variation reflects genetic variation at other nuclear loci, and which genetic diversity estimator for microsatellites provides the most accurate prediction of genome-wide diversity. positive but sometimes weak correlations between expected microsatellite heterozygosity and snp diversity in nuclear gene sequences have been reported at the population level in salmon and several carnivore species, as well as different rice varieties and sheep breeds. most of these studies, however, have investigated only a limited number of snps. the outcome of the comparison of snp versus microsatellite diversity in these studies was strongly affected by the number of snp markers used. studies in which low snp numbers were compared to microsatellites found that the latter had more power to infer differences in genetic summary statistics or found similar results when approximately snps were used. in contrast, studies using larger numbers of snps found that snps performed better than microsatellites. many of these studies used existing genotyping arrays for snp detection. these may, however, cause ascertainment bias as a consequence of the overrepresentation of common snps. to date, no unbiased whole-genome re-sequencing approach has been used for comparison. studies based on reduced representation libraries, which sample a subset of all snps of the genome, showed that snps have more power than microsatellites, e.g. to detect heterozygosity–fitness correlations in natural populations of oldfield mice. further, demographic inferences drawn from radseq-derived snps in bumble bees reflected important long-term differences in population size better than microsatellites, which instead signalled either recent demographic changes or mutational processes. because of the widespread application of microsatellite markers both in basic research and practical conservation, it is important to evaluate the tenet that microsatellite variation adequately reflects genome-wide genetic diversity, especially for situations in which only a limited number of markers are used, as is often the case in conservation genetics, where on average only microsatellites are used per study. it is further relevant to evaluate the power of next-generation sequencing based genotyping approaches to infer genome-wide diversity and population structure, e.g. to estimate the number of snps required to achieve accurate and consistent estimates of genome-wide diversity. we used two types of microsatellite markers as well as genome-wide snp variation in the meadow rock cress, arabidopsis halleri, to compare estimates of genetic diversity and differentiation. overall, we genotyped individuals of a. halleri from nine natural populations using microsatellite markers, which is above the average number of microsatellites typically used in population and conservation genetic studies. we compared them to a pooled whole-genome re-sequencing approach and tested whether estimates of genetic variation derived from microsatellite polymorphisms are valid and useful proxies of genome-wide genetic variation and differentiation. specifically, we tested whether estimates from both marker types were correlated and had similar absolute values. further, we used down-sampling to assess how many random and presumably unlinked snps are required to calculate accurate estimates of genome-wide diversity. methods study system arabidopsis halleri is a perennial, insect-pollinated, strictly outcrossing and functionally self-incompatible herb with a wide geographic distribution from central europe to eastern asia. it grows in diverse habitats including mountain slopes, grassy meadows, forest margins and rocky crevices and has been widely used as a model to study heavy metal tolerance. the species is diploid with 2n =  and has an estimated genome size of mbp. sampling and dna extraction leaf tissue from individuals each of nine populations of a. halleri was sampled in south-eastern switzerland and northern italy. the selected samples size per population should allow to accurately estimate population-specific genetic diversity and differentiation. a minimal distance of two, but preferably 4 m, was maintained between collected individuals, as genetic structure and diversity may be affected by clonal growth when plants are separated by less than one meter. leaf samples were dried in silica gel, and dna was extracted with the dneasy plant mini kit according to the manufacturer’s protocol. dna concentrations were measured with a qubit® fluorometer, and dna quality was examined using a nanodrop spectrophotometer as well as % agarose gels stained with gelred. microsatellite analyses the samples were genotyped using microsatellite markers in three multiplex pcrs, each amplifying either six or eight microsatellite markers. the first two multiplex sets included microsatellite loci that were originally developed for a. thaliana, hereafter referred to as “cross-species” microsatellites. some of these primer sequences were adapted to a. halleri by comparing them to our own de-novo assembly of the a. halleri genome using igv, identifying potential mismatches and changing the primer sequences accordingly. further, eight microsatellite primer pairs that were specifically developed for a. halleri were combined in a third multiplex set. these markers are hereafter referred to as “species-specific” microsatellites. detailed lab protocols can be found in the additional methods. alleles were called using genemapper. estimating microsatellite-based genetic diversity for every marker and population, we assessed the following population genetic parameters. the inbreeding coefficient f is and its p-value, which indicate whether markers or populations deviate from hardy–weinberg equilibrium, were calculated with genodive b using the heterozygosity-based g is statistic with permutations and applying bonferroni correction for multiple testing. null allele frequencies were calculated with freena. pairwise values of genetic differentiation among populations, f st, were calculated based on allele identity with genepop., whereas allele frequencies, expected heterozygosity and mean number of alleles per locus were quantified with genetix. we consistently genotyped individuals per population, therefore, a r did not have to be corrected with a rarefaction approach. all population parameters were computed for three different marker sets including all, the cross-species, and the species-specific microsatellite markers. to infer marker bias, we tested for quantitative differences in estimates of ssr-h e and a r estimated in each population from cross-species and species-specific markers using a paired t-test in r. for the relative comparison of ssr-h e derived from cross-species and species-specific markers, we used a pearson’s correlation test in r. to test whether population-specific estimates of ssr-h e obtained from the different microsatellite types differ, we used a pairwise wilcoxon signed-rank test in r. by plotting ssr-h e medians and quantiles for each population, we inferred whether non-significant differences were caused by high variance. when the absence of significant differences between populations was obviously caused by overly high variances in the genetic diversity estimates computed per microsatellite marker and population, we interpreted this as variance bias. p-values were adjusted for multiple testing using bonferroni correction. pool-seq and illumina read processing pooled next-generation sequencing has been shown to produce accurate population-specific allele frequencies. for ngs, individually extracted high-quality dna was equimolarly pooled using the same individuals from the nine populations as for the microsatellite genotyping presented above. these nine pools were high-throughput sequenced and further processed as described below and, in more detail, in fischer et al.. for a subset of snps and populations of the present dataset, the accuracy of exactly the same pool-seq approach had been validated: differences in estimates of population-specific allele frequencies compared to those from individual genotyping were on average less than 4%. library preparation and sequencing on an illumina hiseq were performed by gatc biotech and the quantitative genomics facility. forward and reverse raw reads were trimmed for tags and adaptors with cutadapt. phred-type quality scores q were used for quality trimming with the fastx toolkit. the separately trimmed forward and reverse reads were then re-synchronized to pairs with an in-house perl script. only paired sequences were used for further analysis. read mapping, snp calling and genome-wide population genetic estimates to estimate genome-wide genetic diversity and differentiation for all nine populations of a. halleri, reads were mapped to the a. thaliana reference genome using bwa aln, allowing for 10% mismatch, and sampe. all ambiguously mapped reads were removed and the remaining high-quality reads were sorted with samtools.. snps were called for the nine populations by producing mpileup files with samtools. to obtain population-specific genome-wide estimates of genetic diversity, we first calculated watterson’s theta, an estimator that takes into account the number of segregating sites to estimate the population mutation rate. θ watterson was calculated on a gene-by-gene basis for each population using the gene and exon annotation of tair. the perl script ‘variance-at-position.pl’ of the software package popoolation was used with the mpileup file of each population and the.gtf annotation file to calculate exonic θ watterson. this approach provides unbiased estimates for pooled samples as it corrects for coverage. exon based θ watterson is a conservative estimate of genetic diversity as it infers diversity from genomic regions that are predominantly under purifying selection, hence from slightly less diverse regions than the rest of the genome. exon-based diversity estimates are of direct relevance for the adaptive potential of a population, because exons harbour functionally relevant polymorphisms that allow populations to adapt to changing environments. to accurately estimate allele frequencies for estimating θ watterson, minimum counts for minor alleles were set to two to account for sequencing errors, leading to a minor allele frequency threshold of. the minimum coverage per site within populations was set to 20×, which mimics the number of individuals. to further correct for potential errors caused by repeated sequences, a maximum coverage of 400× per population was used as threshold for snp identification. in order to be included in the genome-wide estimates of gene diversity, 50% of all snps within a gene had to reach the above-mentioned thresholds in all nine populations. for all analyses, pool size per population was set to because diploid genomes were represented in each population pool. second, we calculated genome-wide snp-based expected heterozygosity, taking all snps into account, not only those located in exons. mpileup files were synchronized and filtered for base quality with the perl script ‘mpileup2sync.pl’ of popoolation. next, major and minor allele frequencies were calculated with the script snp-frequency-diff.pl. the coverage threshold was the same as mentioned above, except that the minor allele count was set to four, as all nine populations were jointly used to infer minor allele frequencies, leading to a more sensitive, but less error-prone minor allele frequency threshold of. we only used bi-allelic snps and calculated the average genome-wide snp-h e as snp-he=1n∑i=1n2pi where n is the number of snps, and p i is the minor allele frequency of the ith allele. this approach assumes hardy–weinberg equilibrium within populations. to infer potential demographic events that could strongly influence genetic diversity within populations, we calculated exome-wide tajima’s d using the tair gene annotation and the perl script ‘variance-at-position.pl’ in popoolation. it is suggested to use a coverage threshold of less than three times smaller than the pool size, which is in our case 13×. a negative genome-wide tajima’s d is indicative of an expansion after a bottleneck, whereas a positive d is compatible with a scenario of a decrease in population size. to test whether the average of the resulting distribution of tajima’s d was significantly different from zero, we used t-tests against random normal distributions with an average of zero and the same standard deviation as observed in the real data of each population. estimates of pairwise population genetic differentiation were calculated with ‘fst-sliding.pl’ in popoolation. average values of pairwise f st were calculated using the same parameters as mentioned above for the estimates of snp-h e as explained in detail in fischer et al.. comparisons of genetic diversity and differentiation derived from microsatellites and genome-wide snps to explore associations between estimates of genetic diversity derived from microsatellites and snps, we performed pearson’s correlations of microsatellite-based allelic richness and expected heterozygosity. estimates of genome-wide snp diversity were derived from exon sequences and genome-wide snp expected heterozygosity. further, to account for possible confounding effects due to linkage, associations were also tested for a subset of snps, consisting of every 50th snp. we performed mantel tests to check for correlations between values of f st derived from genome-wide snps and all, cross-species, and species-specific microsatellite markers. the same analysis was used to assess correlations between values of f st derived from species-specific and cross-species microsatellite markers. all analyses were performed with permutations using ecodist. in r. finally, we used paired t-tests implemented in r to quantitatively evaluate whether h e and f st derived from microsatellite markers and snps significantly differ. estimating the number of unbiased snps required for accurate estimates of genetic diversity we used a down-sampling procedure to estimate snp-h e with the aim to infer the required number of randomly selected and unlinked snp markers to obtain accurate estimates of genetic diversity and to reliably rank populations according to their genetic diversity. thus, each population was resampled for the same k random snp markers drawn from the pool of more than million snps. for each value of k varying between and, we created random subsamples of k snp markers. we then computed the mean expected heterozygosity and 95% confidence intervals for each value of k observed in each of the subsamples. obtained results were used to draw curves representing the variation of the estimated h e as a function of genotyping effort in each population. we then identified the number of snps for which the upper and lower confidence intervals for expected heterozygosity fell below ±, ±, and ±.  microsatellite diversity twenty microsatellite loci were initially used to characterise a. halleri individuals from nine populations. we excluded marker ah from further analyses, because it deviated significantly from hardy–weinberg equilibrium and exhibited an estimated null allele frequency of 10%. the remaining microsatellite markers harboured alleles and only % missing data. allelic richness per locus ranged between and per population, with an average of, and expected heterozygosity ranged from to per microsatellite marker. further details are given in table  and additional file 4: table s population allele frequency distributions were fairly noisy, see additional file 5: figure s1a. none of the nine populations showed significant deviation from hardy–weinberg equilibrium after bonferroni correction.table population genetic parameters inferred from microsatellites and genome-wide snps for nine populations of arabidopsis halleri. allelic richness, expected heterozygosity and inbreeding coefficient f is including its one-sided p-value are given. no f is value was significantly different from zero after bonferroni correction. θ watterson was calculated for genes. expected heterozygosity was calculated from all snps across the genome. tajima’s d was calculated for genes, p-values refer to deviations from zero a r e f is f is θ watterson e all twelve cross-species microsatellite markers, initially developed for a. thaliana, successfully amplified in a. halleri. a r at cross-species microsatellite markers was alleles per marker and population and thus significantly lower than a r at species-specific markers. the same pattern was observed for ssr-h e, which was and, respectively. no significant correlation was observed among estimates of h e inferred from cross-species and species-specific microsatellite markers.fig. comparison of estimates of genetic diversity derived from cross-species and species-specific microsatellite markers for a allelic richness and b expected microsatellite heterozygosity. c estimates of ssr-h e inferred separately from cross-species and species-specific microsatellite markers were not significantly correlated. dots are labelled with population codes. d no significantly different estimates of h e were observed among populations after bonferroni correction. without correction for multiple testing only population aha and ahan showed significantly different estimates of h e the variance in the estimates of ssr-h e among microsatellite markers within and among population was so high that no significant differences in genetic diversity among populations could be inferred after bonferroni correction. without bonferroni correction, only populations aha and ahan differed significantly in their estimates of h e. similar results were found when cross-species and species-specific microsatellites were analysed separately and corrections for multiple testing were performed. average pairwise f st for the microsatellite loci was; species-specific and cross-species microsatellites markers did not deviate significantly from each other. illumina sequencing and genome-wide diversity the illumina sequencing yielded a total of, 100-bp paired-end reads corresponding to, nucleotides. after quality filtering and trimming,, paired-end reads were mapped to the a. thaliana reference genome, from which organellar dna was excluded. the average coverage per site, after filtering with the defined thresholds, was × with a range of population-wise coverage of to  ×. we detected, snps, which were used for the calculation of pairwise f st, and, bi-allelic snps, which were used for the calculation of population-specific snp-h e. all populations had even population-specific allele frequency distributions. values of pairwise f st ranged between and, and population-specific snp-h e was between and. overall, and genes fulfilled our thresholds of coverage for calculating θ watterson and tajima’s d, respectively. θ watterson ranged from to, and tajima’s d values were all slightly negative, ranging from − to −. only two of the nine populations showed no significant deviation from zero. in other words, most populations showed weak demographic changes probably related to a bottleneck with later expansion. comparisons of genetic diversity estimates derived from microsatellites and genome-wide snps no significant correlation was observed between population-specific estimates of h e derived from microsatellites and genome-wide snps, independent of whether all snps or a subset of presumably unlinked snps were used. estimates of snp-h e were overall significantly lower than those of ssr-h e. ssr-h e was also not significantly correlated with θ watterson. the correlation coefficient was higher when using only species-specific markers, but lower when only cross-species markers were used, though neither of the two correlations was significant.fig. relationships between population genetic parameters for nine populations of arabidopsis halleri estimated on the basis of microsatellites and million genome-wide single nucleotide polymorphisms. a relationship between expected genome-wide snp heterozygosity and expected microsatellite heterozygosity. the insert illustrates differences in h e between snps and microsatellites across the nine a. halleri populations. b relationship between genome-wide diversity estimated by average θ watterson across genes and expected heterozygosity for all microsatellites and separately for species-specific and cross-species microsatellites. c comparisons of allelic richness with genome-wide snp variation estimated by θ watterson in marked contrast to heterozygosity, a r of microsatellite markers was significantly correlated with θ watterson. the correlation based on ranks was still significant, but slightly weaker than non-ranked comparisons. the snp-based genome-wide diversity estimates, θ watterson and snp-h e, were highly correlated, even though values of θ watterson were derived exclusively from coding regions, and estimates of snp-h e were calculated from more than two million snps across the entire genome. in fact, θ watterson estimates inferred from introns and intergenic regions were highly correlated to estimates of exon-based θ watterson.fig. pearson’s correlation between exome-wide θ watterson and genome-wide expected snp heterozygosity in arabidopsis halleri comparison of genetic differentiation estimates derived from microsatellite versus genome-wide snp variation mantel tests revealed highly significant correlations between values of pairwise f st derived from genome-wide snp data and microsatellite markers. the best correlation was achieved when using all microsatellites. however, values of f st derived from microsatellite markers were -fold higher than those from snps and were significantly different. if we split the microsatellite markers into species-specific and cross-species, the correlations were slightly weaker for cross-species microsatellites and for species-specific microsatellites. the correlation among species-specific and cross-species microsatellites was high.fig. comparison of pairwise population genetic differentiation among nine populations of arabidopsis halleri estimated with different genetic markers. a comparison of f st derived from, genome-wide snps with f st inferred on the basis of all microsatellite markers, and separately for cross-species, and seven species-specific microsatellites. given are correlation coefficients and p-values from mantel tests. b estimates of f st inferred from snps were significantly lower than those based on microsatellites estimating the number of unbiased snps required for accurate estimates of genetic diversity snp down-sampling revealed that the upper and lower confidence intervals for expected heterozygosity fell below ±, ±, and ± with and random snps, respectively. to accurately and consistently rank all a. halleri populations according to their genome-wide diversity, snps were required. however, populations aha and ahan could not be distinguished, as they have the same genome-wide snp-h e.fig. the effect of the number of randomly selected snp markers on the estimation of h e represented for each of nine populations of arabidopsis halleri. light colours indicate the 95% confidence intervals as inferred from randomly drawn replicates. dashed horizontal lines indicate the numbers of snps for which 95% confidence intervals for snp-h e are below ±, ±, and ± discussion the application of microsatellite markers is widespread in population and conservation genetic studies. however, ngs-based snp genotyping approaches are rapidly developing and can be applied to a wide diversity of model and non-model organisms. our comparative analysis of genetic diversity estimates based on microsatellites and genome-wide snps revealed interesting differences. most importantly, we found no significant correlation between expected microsatellite heterozygosity, an estimator of genetic diversity that is widely used and reported in microsatellite studies, and genome-wide snp diversity. this finding indicates that ssr-h e does not adequately reflect genome-wide genetic diversity in the investigated populations of a. halleri. in contrast, microsatellite allelic richness was a much better proxy for genome-wide diversity. further, genetic differentiation in terms of f st estimated from microsatellite variation correlated reasonably well with that based on genome-wide snp data. however, absolute values of the different summary statistics inferred from different marker types varied considerably. our results do not question the usefulness of microsatellites per se, but point to research questions for which snps may be better suited than microsatellite markers, given the availability of robust and cost-effective high-throughput sequencing-based snp genotyping approaches. thanks to the massive advances in sequencing technology, thousands of snps can efficiently be genotyped in any given organism, and these may improve our ability to adequately estimate genetic diversity and differentiation. previous studies found that the potential of snps to resolve population genetic structure strongly depends on their number. indeed, in studies using a low number of assay-based snps, microsatellites performed similarly well or better than snps. however, in studies using larger snp numbers, especially when they were derived from ngs-based approaches, the relative performance of snps clearly improved. these studies show that a large number of snps compensates for the lower information content of these typically bi-allelic markers compared to more polymorphic microsatellite markers. our study shows that a few thousand snps are enough to accurately estimate genome-wide diversity in terms of h e. evidence available from animal species as well as from this study reveals that there is no or only weak congruence between estimates of heterozygosity derived from microsatellites and snps. indeed, theory suggests that an association of heterozygosity estimates between microsatellites and genome-wide snps is not expected a priori. according to ljungqvist et al., this association is shaped by identity disequilibrium, i.e. the non-random association of diploid genotypes between loci. simulations indicated that a strong positive correlation only emerges when the studied populations are characterized by substantial identity disequilibrium, as was the case in the studies on salmon and carnivores. the absence of this correlation in a. halleri suggests that identity disequilibrium is weak or absent, which is compatible with the existence of a strong self-incompatibility system in this species and observed non-significant values of f is. this finding indicates that heterozygosity estimates based on microsatellites might not be good surrogates for genome-wide diversity in outbreeding species. further, we found that the marker-specific variation of ssr-h e is too high to distinguish populations based on their levels of genetic diversity and that even among microsatellite markers originally developed for different taxa, important discrepancies can be observed. this variance bias is especially strong when a low number of microsatellites is used, as is the case in many population and conservation genetic studies. the larger sampling variance associated with a limited number of microsatellite markers is evidenced by their allele frequency distributions, which are much noisier than those derived from the snp data set. an alternative estimator to h e is microsatellite allelic richness, one of the simplest estimators of genetic diversity available. a r was significantly correlated with genome-wide snp diversity in our study and thus appears to be a useful proxy of genome-wide genetic diversity. congruent results of a r and snp diversity were also reported in other studies and might be explained by several reasons. first, ssr-h e estimates are based on few markers with noisy allele frequency distributions and represent a proportion, ranging between and whereas allelic richness is an infinite count. accordingly, one additional allele, especially when it is rare and many alleles are already present, does not strongly influence ssr-h e, but will affect a r estimates, making the latter a more sensitive estimator of diversity. moreover, especially for microsatellites with a high number of alleles, accurate estimates of population-specific heterozygosity can be problematic, and stochasticity may have a strong impact on estimates of h e at the lower range of allelic diversity. finally, a r is more sensitive to population bottlenecks than h e. therefore, a r better reflects the population’s demographic history and hence is a more relevant estimator of genetic diversity to predict the short-term survival of a population. despite the good performance of a r as a proxy of genome-wide diversity, it was not sufficient to accurately rank populations according to their genetic diversity. consequently, the identification of conservation units or decision taking for conservation actions based on microsatellite-derived rankings of genetic diversity may be misguiding. nevertheless, using microsatellite-derived a r rather than ssr-h e provides more accurate estimates of genome-wide genetic diversity derived from a limited number of microsatellite markers. these considerations are strengthened by a simulation study, which found allelic richness to be two to four times more powerful than h e for the identification of a temporal genetic decline in a population. a possible reason for the deviation of estimates of genetic diversity derived from microsatellite and genome-wide snp data could be the influence of very recent demographic changes. as a consequence of their higher mutation rate, microsatellites respond more strongly to recent demographic events than genome-wide snps, and snps uncover a different and likely older demographic history. evidence for this hypothesis has also been presented for bumble bees and may indicate a fruitful application of microsatellites in analyses of populations that may have undergone very recent demographic changes. however, it is important to note that loci with high mutation rate may violate demographic model assumptions, such as mutation–migration–drift equilibrium. further, for the long-term survival of populations, genetic diversity in coding and regulatory sequences is arguably more relevant than microsatellite diversity, because most microsatellites are located in non-coding regions and are mostly selectively neutral, hence of less evolutionary importance. consequently, estimates of genome-wide snp diversity better reflect functionally important and potentially adaptive genetic variation, and should therefore be used preferentially, especially in conservation genetics studies. in a. halleri the long-term demographic history inferred from genome-wide data indicates a bottleneck with later expansion for most populations, as values of tajima’s d were slightly negative and significantly different from zero. the origin of the microsatellite markers used, i.e. whether they are species-specific or cross-species markers, may further impact estimates of genetic diversity. in this study, species-specific microsatellites displayed significantly higher a r and ssr-h e than cross-species markers originally developed for a. thaliana. further, species-specific markers resulted in more accurate estimates of genetic diversity, but less accurate estimates of divergence among populations of a. halleri. hence, the practicability of microsatellites for population genetic studies is limited and difficult to assess a priori. high estimates of genetic diversity derived from species-specific microsatellite markers may be a consequence of ascertainment bias caused by selecting the most polymorphic markers, whereas cross-species microsatellites are mostly chosen based on their amplification success in the study species. the consequence of this ascertainment bias is evident in this study, as estimates of ssr-h e for cross-species microsatellites are not significantly correlated with ssr-h e for species-specific microsatellites. these differences further emphasize our inference that marker choice may substantially bias estimates of genetic diversity and may invalidate comparisons between populations or species, most notably when different microsatellite loci are assessed. a different pattern emerges for estimates of population genetic differentiation in terms of pairwise f st. we found a significant positive correlation between values of pairwise f st derived from microsatellites and genome-wide snps. similar findings were reported for salmons and threespine sticklebacks. a reason for the better correlation between estimates of genetic differentiation compared to estimates of genetic diversity may be that more values are involved in pairwise comparisons, and that differences in allele frequencies of the common alleles are more important for the accurate estimation of genetic differentiation than those of rare alleles. importantly, estimates of f st derived from microsatellites were consistently and substantially higher than those based on genome-wide snps. this seems counterintuitive, because multi-allelic microsatellite markers with high mutation rates should cause lower f st values than low-diversity markers like snps. however, pooled whole-genome re-sequencing studies with high coverage, such as this one, also detect rare variants; these low-frequency snps reduce overall f st. overall, we consider pairwise population genetic differentiation estimated from microsatellites a useful proxy for genome-wide differentiation, but only in relative and not in absolute terms. this finding has serious implications, because absolute values of f st continue to be frequently used to infer indirect estimates of gene flow and migration, even though estimates of gene-flow should not be derived from f st. further, this marker-specific difference in f st estimates has a major impact on comparative studies of the divergence of quantitative traits, known as q st–f st comparisons, because the inference of the role of natural selection and genetic drift as causes of population genetic differentiation in complex polygenic traits is biased. while our results suggest that microsatellites should not be used for estimating genome-wide heterozygosity, we emphasize that microsatellites remain useful molecular markers for other applications. for example, microsatellites perform very well in genetic stock identification or paternity analysis owing to their high variability and may therefore continue to play an important role in molecular ecology. however, before embarking on a molecular analysis, it remains a key issue to carefully assess the inherent strengths and limitations associated with different molecular markers. only then it is possible to select the most appropriate method for a given ecological or evolutionary question. in contrast to microsatellite-derived data, estimates of genome-wide diversity inferred from whole-genome re-sequencing data, e.g. exome-wide θ watterson, intronic and intergenic θ watterson, or genome-wide snp-h e, were highly correlated with each other and led to the same ranking of populations. even though values of θ watterson were either derived exclusively from coding regions or intronic and intergenic regions, and snp-h e was calculated from positions across the whole genome, their estimates were highly congruent. the slight variation observed among θ watterson and snp-h e might be explained by differences in the demographic history among populations, because the demographic history of a population has a stronger influence on θ watterson than on snp-h e or snp nucleotide diversity estimates, as rare alleles are more likely to be lost during a bottleneck than common ones. further, we found in our whole-genome re-sequencing study that the confounding effects of genetic linkage are negligible in the highly outcrossing and self-incompatible a. halleri, most likely because the small-scale linkage effects are compensated by the large numbers of unlinked snps. thus, our snp-h e estimates based on a subset of putatively unlinked snps were nearly identical to the estimates inferred for all snps, see fig. 2a and additional file 8: figure s. similar to our genome re-sequencing study, approaches that use reduced representation libraries to sample a subset of genome-wide snps can accurately estimate genome-wide heterozygosity. as a consequence of the much smaller proportion of the genome surveyed with such approaches, however, care should be taken to avoid confounding effects of linkage, for example by considering only one snp per rad-locus. for example, the inbreeding coefficient of a known pedigree in oldfield mice showed strong concordance with the inferred estimates of heterozygosity obtained from radseq snps. this result indicates that, as long as a sufficiently large number of unbiased ngs-based snps is analysed across the genome, snp estimates accurately reflect genome-wide diversity in natural populations. therefore, approaches like radseq, pool-seq and whole-genome re-sequencing at low coverage are more appropriate than array-based snp approaches, which may be affected by strong ascertainment bias. an important question to consider in many studies may be the number of snps that are needed to estimate genetic diversity. our down-sampling approach indicated that the number of random snps that are required to resolve genetic diversity difference among populations range from to snps. this number is in the range of snps that can be inferred with standard radseq protocols also in non-model organisms. however, to differentiate among populations with very similar levels of genetic diversity, we required approximately snps. thus, large snp datasets that are ideally identified de novo through ngs approaches are highly suitable to distinguish, for example, between populations differing in genetic diversity and may therefore support decision-making in conservation management. a further advantage of genome-wide snp data is that they not only allow one to estimate neutral genetic diversity, but also to identify adaptive genetic variation, which is considered essential for delimitating conservation units. the large technical advances in nucleotide sequencing technology in recent years have not only massively increased the number of nucleotides that can be sequenced per individual or population, but have also led to reduced costs per nucleotide to the extent that screening a handful of microsatellite markers may be as expensive as surveying thousands of snps using latest ngs-based genotyping technologies.  this case study in the perennial and outcrossing plant a. halleri reveals that genetic diversity estimated from microsatellite markers, notably expected heterozygosity, may not adequately reflect genome-wide genetic diversity estimated from single-nucleotide polymorphisms and may therefore be a poor proxy for genome-wide estimates of genetic diversity. possible causes include the limited number of microsatellite markers used, marker ascertainment bias, as well as the high variance in microsatellite-derived diversity estimates. interestingly, microsatellite allelic richness was found to be a reasonable proxy for genome-wide diversity, but the absolute ranking of populations was still inconsistent. estimates of genetic differentiation among populations derived from microsatellites were consistently higher than snp-based estimates but were significantly correlated with the latter. our results do not question the usefulness of microsatellites per se, but point to research questions for which ngs-derived snps may be better suited than microsatellite markers, given the availability of robust and cost-effective snp genotyping approaches based on high-throughput sequencing. as a consequence, we recommend using genome-wide analyses of snp diversity when the inference and comparison of genetic diversity within and among populations and species is the goal of a study. a few thousand ngs-derived snps are sufficient for this purpose and this number of unbiased snps can nowadays easily been obtained also for non-model species, for example by using a reduced representation sequencing approach such as radseq. additional files additional file 1: table s sampling locations of the nine study populations of arabidopsis halleri. additional file 2: table s characteristics of microsatellite markers used for arabidopsis halleri. source refers to the species for which a microsatellite marker was originally developed. sequences of forward and reverse primers are shown in 5′ - 3′ direction. the fluorescent dyes used in multiplex pcrs and allele size ranges in base pairs are given for each locus. additional file 3: additional methods. additional file 4: table s locus-specific summary statistics for all microsatellite markers genotyped in arabidopsis halleri. estimated frequency of null alleles, total number of alleles found among the individuals, expected heterozygosity, and inbreeding coefficient f is with one-sided p-value are given. additional file 5: figure s population specific allele frequency distributions in nine populations of a. halleri for microsatellite markers. minor allele frequency distributions in the same nine populations across, snps. histograms are labelled with population codes. additional file 6: figure s box-whisker plot of expected heterozygosity calculated for each microsatellite marker genotyped in arabidopsis halleri. for cross-species microsatellites, no significantly different estimates of h e were observed after bonferroni correction, as a consequence of the high variance. without correction for multiple testing, only population aha had a significantly different h e estimates from aha and ahan. for species-specific markers, no significant differences in h e estimates were observed after bonferroni correction. without multiple corrections, only aha was significantly different from ahan additional file 7: figure s observed versus expected distributions of tajima’s d for nine populations of arabidopsis halleri. the expected normal distribution consists of values with an average of zero and the same standard deviation as the real dataset. p-values indicate whether there was a significant deviation of the average from zero using a t-test. additional file 8: figure s relationships between expected snp heterozygosity and expected microsatellite heterozygosity of arabidopsis halleri estimated on the basis of microsatellites and genetically unlinked genome-wide single nucleotide polymorphisms. for this analysis only every 50th snp of the more than million snps was used. the median distance among two neighbouring snps is 1600 bp on the a. thaliana reference genome. additional file 9: figure s pearson’s correlation between exome-wide θ watterson and θ watterson calculated from introns and intergenic regions in arabidopsis halleri. additional file 10: figure s comparison of pairwise population genetic differentiation among nine populations of arabidopsis halleri based on cross-species microsatellite markers and seven species-specific microsatellite markers. r mt represents the correlation coefficient of the mantel test and p the significance of the correlation. abbreviations armicrosatellite allelic richness heexpected heterozygosity ngsnext-generation sequencing pool-seqpooled next-generation sequencing radseqrestriction-site associated dna sequencing snpsingle nucleotide polymorphism ssrsimple sequence repeat θwattersonwatterson’s theta we thank c. beisel for supporting our illumina sequencing efforts at the quantitative genomics facility, d-bsse, eth zürich. we also thank a. bösch, e. schnyder, j. zimmermann, d. zulliger, n. quebre, and n. aellen for sampling a. halleri populations, and c. michel, r. graf, and the genetic diversity centre zürich for experimental support, a. tedder for helpful input, as well as s. fior and n. zemp for bioinformatics support. we acknowledge the valuable comments made by anonymous reviewers on an earlier version of this article. funding this study was funded by the swiss national science foundation and the adaptation to a changing environment center of eth zurich. availability of data and materials the datasets supporting the conclusions of this article are available in the dryad digital repository,, in the european nucleotide archive and in the sequence read archive. authors’ contributions mcf, cr, and aw designed the study; fg, mcf, and aw conducted field-work; mcf, cr, ml, and mr performed the research and analyzed the data; mcf, cr, and aw wrote the manuscript with substantial contributions from ml, mr, fg, and rh. aw, kks and rh conceived the project and obtained funding. all authors read and approved the final manuscript. competing interests the authors declare that they have no competing interests. consent for publication not applicable. ethics approval and consent to participate not applicable.
possible causes include the limited number of microsatellite markers used, marker ascertainment bias, as well as the high variance in microsatellite-derived estimates. twelve of these markers were originally developed for arabidopsis thaliana and eight for a. halleri. our analyses revealed that estimates of genetic diversity and differentiation derived from cross-species and species-specific microsatellites differed substantially and that expected microsatellite heterozygosity was not significantly correlated with genome-wide snp diversity estimates in a. halleri. in species comparable to a. halleri, a few thousand snps are sufficient to achieve this goal. microsatellite markers are widely used for estimating genetic diversity within and differentiation among populations. in contrast, genome-wide snp data provided unbiased estimates of genetic diversity independent of whether genome- or only exome-wide snps were used. we recommend that future analyses of genetic diversity within and differentiation among populations use randomly selected high-throughput sequencing-based snp data to draw conclusions on genome-wide diversity patterns. microsatellitesssrarabidopsis hallerigenetic diversityexpected heterozygositysnpspopulation genomicswhole-genome re-sequencingpool-seqconservation unitshttp://dx.doi.org/ /501100001711schweizerischer nationalfonds zur förderung der wissenschaftlichen forschungcrsi33_127155widmer alex adaptation to a changing environment fellowship, eth zurichissue-copyright-statement© the author 2017 however, it has rarely been tested whether such estimates are useful proxies for genome-wide patterns of variation and differentiation. estimates of genetic differentiation among populations based on both marker types were correlated, but microsatellite-based estimates were significantly larger than those from snps. here, we compared microsatellite variation with genome-wide single nucleotide polymorphisms to assess and quantify potential marker-specific biases and derive recommendations for future studies. further, we inferred that a few thousand random snps are sufficient to reliably estimate genome-wide diversity and to distinguish among populations differing in genetic variation. instead, microsatellite allelic richness was a better proxy for genome-wide snp diversity. overall, we genotyped arabidopsis halleri individuals from nine populations using microsatellite markers.
 microsatellite markers are widely used for estimating genetic diversity within and differentiation among populations. however, it has rarely been tested whether such estimates are useful proxies for genome-wide patterns of variation and differentiation. here, we compared microsatellite variation with genome-wide single nucleotide polymorphisms to assess and quantify potential marker-specific biases and derive recommendations for future studies. overall, we genotyped arabidopsis halleri individuals from nine populations using microsatellite markers. twelve of these markers were originally developed for arabidopsis thaliana and eight for a. halleri. we further characterized million snps across the genome with a pooled whole-genome re-sequencing approach.  our analyses revealed that estimates of genetic diversity and differentiation derived from cross-species and species-specific microsatellites differed substantially and that expected microsatellite heterozygosity was not significantly correlated with genome-wide snp diversity estimates in a. halleri. instead, microsatellite allelic richness was a better proxy for genome-wide snp diversity. estimates of genetic differentiation among populations based on both marker types were correlated, but microsatellite-based estimates were significantly larger than those from snps. possible causes include the limited number of microsatellite markers used, marker ascertainment bias, as well as the high variance in microsatellite-derived estimates. in contrast, genome-wide snp data provided unbiased estimates of genetic diversity independent of whether genome- or only exome-wide snps were used. further, we inferred that a few thousand random snps are sufficient to reliably estimate genome-wide diversity and to distinguish among populations differing in genetic variation.  we recommend that future analyses of genetic diversity within and differentiation among populations use randomly selected high-throughput sequencing-based snp data to draw conclusions on genome-wide diversity patterns. in species comparable to a. halleri, a few thousand snps are sufficient to achieve this goal. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords microsatellitesssrarabidopsis hallerigenetic diversityexpected heterozygositysnpspopulation genomicswhole-genome re-sequencingpool-seqconservation unitshttp://dx.doi.org/ /501100001711schweizerischer nationalfonds zur förderung der wissenschaftlichen forschungcrsi33_127155widmer alex adaptation to a changing environment fellowship, eth zurichissue-copyright-statement© the author 2017 