Quantifying ascertainment bias and species-specific length differences in human and chimpanzee microsatellites using genome sequences

Surveys of variability of homologous microsatellite loci among species reveal an ascertainment bias for microsatellite length where microsatellite loci isolated in one species tend to be longer than homologous loci in related species. Here, we take advantage of the availability of aligned human and chimpanzee genome sequences to compare length difference of homologous microsatellites for loci identified in humans to length difference for loci identified in chimpanzees. We are able to quantify ascertainment bias for a range of motifs and microsatellite lengths. Because ascertainment bias should not exist if a microsatellite selected in one species is as likely to be longer as it is to be shorter than its homologue, we propose that the nature of ascertainment bias can provide evidence for understanding how microsatellites evolve. We show that bias is greater for longer microsatellites but also that many long microsatellites have short homologues. These results are consistent with the notion that growth of long microsatellites is constrained by an upper length boundary that, when reached, sometimes results in large deletions. By evaluating ascertainment bias separately for interrupted and uninterrupted repeats we also show that long microsatellites tend to become interrupted, thereby contributing a second component of ascertainment bias. Having accounted for ascertainment bias, in agreement with results published elsewhere, we find that microsatellites in humans are longer on average than those in chimpanzees. This length difference is similar among repeat motifs but surprisingly comprises two roughly equal components, one associated with the repeats themselves and one with the flanking sequences. The differences we find can only be explained if microsatellites are both evolving directionally under a biased mutation process and are doing so at different rates in different closely related species.     

Microsatellite length differences between humans and chimpanzees at autosomal Loci are not found at equivalent haploid Y chromosomal Loci

When homologous microsatellites are compared between species, significant differences in mean length are often noted. A dominant cause of these length differences is ascertainment bias due to selection for maximum repeat number and repeat purity when the markers are being developed. However, even after ascertainment bias has been allowed for through reciprocal comparisons, significant length differences remain, suggesting that the average microsatellite mutation rate differs between species. Two classes of mechanism have been proposed: rapid evolution of enzymes involved in the generation and repair of slippage products (enzyme evolution model) and heterozygote instability, whereby interchromosomal events at heterozygous sites offer extra opportunities for mutations to occur (heterozygote instability model). To examine which of these hypotheses is most likely, we compared ascertainment bias and species length differences between humans and chimpanzees in autosomal and Y chromosomal microsatellites. We find that levels of ascertainment bias are indistinguishable, but that interspecies length differences are significantly greater for autosomal loci compared with haploid Y chromosomal loci. Such a pattern is consistent with predictions from the heterozygote instability model and is not expected under models of microsatellite evolution that do not include interchromosomal events such as the enzyme evolution model.     

Microsatellite Evolution: Testing the Ascertainment Bias Hypothesis

Previous studies suggest the median allele length of microsatellites is longest in the species from which the markers were derived, suggesting that an ascertainment bias was operating. We have examined whether the size distribution of microsatellite alleles between sheep and cattle is source dependent using a set of 472 microsatellites that can be amplified in both species. For those markers that were polymorphic in both species we report a significantly greater number of markers (P < 0.001) with longer median allele sizes in sheep, regardless of microsatellite origin. This finding suggests that any ascertainment bias operating during microsatellite selection is only a minor contributor to the variation observed. Received: 6 January 1997 / Accepted: 19 May 1997     

Microsatellite variation in Drosophila melanogaster and Drosophila simulans: a reciprocal test of the ascertainment bias hypothesis

Interspecific comparisons of microsatellite loci have repeatedly shown that the loci are longer and more variable in the species from which they are derived (the focal species) than are homologous loci in other (nonfocal) species. There is debate as to whether this is due to directional evolution or to an ascertainment bias during the cloning and locus selection processes. This study tests these hypotheses by performing a reciprocal study. Eighteen perfect dinucleotide microsatellite loci identified from a Drosophila simulans library screen and 18 previously identified in an identical Drosophila melanogaster library screen were used to survey natural populations of each species. No difference between focal and nonfocal species was observed for mean PCR fragment length. However, heterozygosity and number of alleles were significantly higher in the focal species than in the nonfocal species. The most common allele in the Zimbabwe population of both species was sequenced for 31 of the 36 loci. The length of the longest stretch of perfect repeat units is, on average, longer in the focal species than in the non-focal species. There is a positive correlation between the length of the longest stretch of perfect repeats and heterozygosity. The difference in heterozygosity can thus be explained by a reduction in the length of the longest stretch of perfect repeats in the nonfocal species. Furthermore, flanking-sequence length difference was noted between the two species at 58% of the loci sequenced. These data do not support the predictions of the directional-evolution hypothesis; however, consistent with the ascertainment bias hypothesis, the lower variability in nonfocal species is an artifact of the microsatellite cloning and isolation process. Our results also suggest that the magnitude of ascertainment bias for repeat unit length is a function of the microsatellite size distribution in the genomes of different species.      

Factors Influencing Ascertainment Bias of Microsatellite Allele Sizes: Impact on Estimates of Mutation Rates

Microsatellite loci play an important role as markers for identification, disease gene mapping, and evolutionary studies. Mutation rate, which is of fundamental importance, can be obtained from interspecies comparisons, which, however, are subject to ascertainment bias. This bias arises, for example, when a locus is selected on the basis of its large allele size in one species (cognate species 1), in which it is first discovered. This bias is reflected in average allele length in any noncognate species 2 being smaller than that in species 1. This phenomenon was observed in various pairs of species, including comparisons of allele sizes in human and chimpanzee. Various mechanisms were proposed to explain observed differences in mean allele lengths between two species. Here, we examine the framework of a single-step asymmetric and unrestricted stepwise mutation model with genetic drift. Analysis is based on coalescent theory. Analytical results are confirmed by simulations using the simuPOP software. The mechanism of ascertainment bias in this model is a tighter correlation of allele sizes within a cognate species 1 than of allele sizes in two different species 1 and 2. We present computations of the expected average allele size difference, given the mutation rate, population sizes of species 1 and 2, time of separation of species 1 and 2, and the age of the allele. We show that when the past demographic histories of the cognate and noncognate taxa are different, the rate and directionality of mutations affect the allele sizes in the two taxa differently from the simple effect of ascertainment bias. This effect may exaggerate or reverse the effect of difference in mutation rates. We reanalyze literature data, which indicate that despite the bias, the microsatellite mutation rate estimate in the ancestral population is consistently greater than that in either human or chimpanzee and the mutation rate estimate in human exceeds or equals that in chimpanzee with the rate of allele length expansion in human being greater than that in chimpanzee. We also demonstrate that population bottlenecks and expansions in the recent human history have little impact on our conclusions.     

Microsatellite variation among divergent populations of stalk-eyed flies, genus Cyrtodiopsis

Microsatellite primers are often developed in one species and used to assess neutral variability in related species. Such analyses may be confounded by ascertainment bias (i.e. a decline in amplification success and allelic variability with increasing genetic distance from the source of the microsatellites). In addition, other factors, such as the size of the microsatellite, whether it consists of perfect or interrupted tandem repeats, and whether it is autosomal or X-linked, can affect variation. To test the relative importance of these factors on microsatellite variation, we examine patterns of amplification and allelic diversity in 52 microsatellite loci amplified from five individuals in each of six populations of Cyrtodiopsis stalk-eyed flies that range from 2.2 % to 11.2% mitochondrial DNA sequence divergence from the population used for microsatellite development. We find that amplification success and most measures of allelic diversity declined with genetic distance from the source population, in some cases an order of magnitude faster than in birds or mammals. The median and range of the repeat array length did not decline with genetic distance. In addition, for loci on the X chromosome, we find evidence of lower observed heterozygosity compared with loci on autosomes. The differences in variability between X-linked and autosomal loci are not adequately explained by differences in effective population sizes of the chromosomes. We suggest, instead, that periodic selection events associated with X-chromosome meiotic drive, which is present in many of these populations, reduces X-linked variation.     

New methods to identify conserved microsatellite loci and develop primer sets of high cross-species utility - as demonstrated for birds

We have developed a new approach to create microsatellite primer sets that have high utility across a wide range of species. The success of this method was demonstrated using birds. We selected 35 avian EST microsatellite loci that had a high degree of sequence homology between the zebra finch Taeniopygia guttata and the chicken Gallus gallus and designed primer sets in which the primer bind sites were identical in both species. For 33 conserved primer sets, on average, 100% of loci amplified in each of 17 passerine species and 99% of loci in five non-passerine species. The genotyping of four individuals per species revealed that 24-76% (mean 48%) of loci were polymorphic in the passerines and 18-26% (mean 21%) in the non-passerines. When at least 17 individuals were genotyped per species for four Fringillidae finch species, 71-85% of loci were polymorphic, observed heterozygosity was above 0.50 for most loci and no locus deviated significantly from Hardy-Weinberg proportions. This new set of microsatellite markers is of higher cross-species utility than any set previously designed. The loci described are suitable for a range of applications that require polymorphic avian markers, including paternity and population studies. They will facilitate comparisons of bird genome organization, including genome mapping and studies of recombination, and allow comparisons of genetic variability between species whilst avoiding ascertainment bias. The costs and time to develop new loci can now be avoided for many applications in numerous species. Furthermore, our method can be readily used to develop microsatellite markers of high utility across other taxa.     

Effect of microsatellite selection on individual and population genetic inferences: an empirical study using cross‐specific and species‐specific amplifications

Although whole‐genome sequencing is becoming more accessible and feasible for nonmodel organisms, microsatellites have remained the markers of choice for various population and conservation genetic studies. However, the criteria for choosing microsatellites are still controversial due to ascertainment bias that may be introduced into the genetic inference. An empirical study of red deer (Cervus elaphus) populations, in which cross‐specific and species‐specific microsatellites developed through pyrosequencing of enriched libraries, was performed for this study. Two different strategies were used to select the species‐specific panels: randomly vs. highly polymorphic markers. The results suggest that reliable and accurate estimations of genetic diversity can be obtained using random microsatellites distributed throughout the genome. In addition, the results reinforce previous evidence that selecting the most polymorphic markers leads to an ascertainment bias in estimates of genetic diversity, when compared with randomly selected microsatellites. Analyses of population differentiation and clustering seem less influenced by the approach of microsatellite selection, whereas assigning individuals to populations might be affected by a random selection of a small number of microsatellites. Individual multilocus heterozygosity measures produced various discordant results, which in turn had impacts on the heterozygosity‐fitness correlation test. Finally, we argue that picking the appropriate microsatellite set should primarily take into account the ecological and evolutionary questions studied. Selecting the most polymorphic markers will generally overestimate genetic diversity parameters, leading to misinterpretations of the real genetic diversity, which is particularly important in managed and threatened populations.     

Long microsatellite alleles in Drosophila melanogaster have a downward mutation bias and short persistence times, which cause their genome-wide underrepresentation

Microsatellites are short tandemly repeated DNA sequence motifs that are highly variable in most organisms. In contrast to mammals, long microsatellites (>15 repeats) are extremely rare in the Drosophila melanogaster genome. To investigate this paucity of long microsatellites in Drosophila, we studied 19 loci with exceptionally long microsatellite alleles. Inter- and intraspecific analysis showed that long microsatellite alleles arose in D. melanogaster only very recently. This lack of old alleles with many repeats indicated that long microsatellite alleles have short persistence times. The size distribution of microsatellite mutations in mutation-accumulation lines suggests that long alleles have a mutation bias toward a reduction in the number of repeat units. This bias causes the short persistence times of long microsatellite alleles. We propose that species-specific, size-dependent mutation spectra of microsatellite alleles may provide a general mechanism to account for the observed differences in microsatellite length between species.     

Ascertainment bias in spatially structured populations: a case study in the eastern fence lizard

Despite increased interest in applying single nucleotide polymorphism (SNP) data to questions in natural systems, one unresolved issue is to what extent the ascertainment bias induced during the SNP discovery phase will impact available analysis methods. Although most studies addressing ascertainment bias have focused on human populations, it is not clear whether existing methods will work when applied to other species with more complex demographic histories and more significant levels of population structure. Here we present findings from an empirical approach to exploring the effect of population structure on issues of ascertainment bias in the Eastern Fence Lizard, Sceloporus undulatus. We find that frequency spectra and summary statistics were highly sensitive to SNP discovery strategy, necessitating careful selection of the initial ascertainment panel. Randomly selected ascertainment panels performed equally well as ascertainment panels chosen to jointly sample geographic, phenotypic, and genetic diversity. Geographically restricted panels resulted in larger biases. Additionally, we found existing ascertainment bias correction methods, which were not developed for geographically structured data sets, were largely effective at reducing the impact of ascertainment bias. Because bias correction methods performed well even when underlying assumptions were violated, our results suggest tools are currently available to analyze SNP data in structured populations.     

Transpecific microsatellites for hard pines

Microsatellites are difficult to recover from large plant genomes so cross-specific utilisation is an important source of markers. Fifty microsatellites were tested for cross-specific amplification and polymorphism to two New World hard pine species, slash pine (Pinus elliottii var. elliottii) and Caribbean pine (P. caribaea var. hondurensis). Twenty-nine (58%) markers amplified in both hard pine species, and 23 of these 29 were polymorphic. Soft pine (subgenus Strobus) microsatellite markers did amplify, but none were polymorphic. Pinus elliottii var. elliottii and P. caribaea var. hondurensis showed mutational changes in the flanking regions and the repeat motif that were informative for Pinus spp. phylogenetic relationships. Most allele length variation could be attributed to variability in repeat unit number. There was no evidence for ascertainment bias.     

Rapid,economical single‐nucleotide polymorphism and microsatellite discovery based on de novo assembly of a reduced representation genome in a non‐model organism: a case study of Atlantic cod Gadus morhua

By combining next‐generation sequencing technology (454) and reduced representation library (RRL) construction, the rapid and economical isolation of over 25 000 potential single‐nucleotide polymorphisms (SNP) and >6000 putative microsatellite loci from c. 2% of the genome of the non‐model teleost, Atlantic cod Gadus morhua from the Celtic Sea, south of Ireland, was demonstrated. A small‐scale validation of markers indicated that 80% (11 of 14) of SNP loci and 40% (6 of 15) of the microsatellite loci could be amplified and showed variability. The results clearly show that small‐scale next‐generation sequencing of RRL genomes is an economical and rapid approach for simultaneous SNP and microsatellite discovery that is applicable to any species. The low cost and relatively small investment in time allows for positive exploitation of ascertainment bias to design markers applicable to specific populations and study questions.     

Factors affecting germline mutations in a hypervariable microsatellite: a comparative analysis of six species of swallows (Aves: Hirundinidae)

Microsatellites mutate frequently by replication slippage. Empirical evidence shows that the probability of such slippage mutations may increase with the length of the repeat region as well as exposure to environmental mutagens, but the mutation rate can also differ between the male and female germline. It has been hypothesized that more intense sexual selection or sperm competition can also lead to elevated mutation rates, but the empirical evidence is inconclusive. Here, we analyzed the occurrence of germline slippage mutations in the hypervariable pentanucleotide microsatellite locus HrU10 across six species of swallow (Aves: Hirundinidae). These species exhibit marked differences in the length range of the microsatellite, as well as differences in the intensity of sperm competition. We found a strong effect of microsatellite length on the probability of mutation, but no residual effect of species or their level of sperm competition when the length effect was accounted for. Neither could we detect any difference in mutation rate between tree swallows (Tachycineta bicolor) breeding in Hamilton Harbour, Ontario, an industrial site with previous documentation of elevated mutation rates for minisatellite DNA, and a rural reference population. However, our cross-species analysis revealed two significant patterns of sex differences in HrU10 germline mutations: (1) mutations in longer alleles occurred typically in the male germline, those in shorter alleles in the female germline, and (2) male germline mutations were more often expansions than contractions, whereas no directional bias was evident in the female germline. These results indicate some fundamental differences in male and female gametogenesis affecting the probability of slippage mutations. Our study also reflects the value of a comparative, multi-species approach for locus-specific mutation analyses, through which a wider range of influential factors can be assessed than in single-species studies.     

Horse microsatellites and their amenability to comparative equid genetics

We investigated the applicability of microsatellite primers, designed in horses, for use in plains and mountain zebras. Fifteen of the 20 tested horse-isolated primer pairs reliably amplified polymorphic loci in two wild equid species. We used this information to assess whether levels of genetic variation and repeat size differed in species from which microsatellites were isolated and in closely related target species. Target equid species exhibited similar levels of genetic variation to the horse, the species from which primers were originally isolated. We show that ascertainment bias results in lower mean and modal repeat size in target species. The data also provide evidence for a bi-directional mutational constraint in allele size across three equid species.     

Vitis Phylogenomics: Hybridization Intensities from a SNP Array Outperform Genotype Calls

Understanding relationships among species is a fundamental goal of evolutionary biology. Single nucleotide polymorphisms (SNPs) identified through next generation sequencing and related technologies enable phylogeny reconstruction by providing unprecedented numbers of characters for analysis. One approach to SNP-based phylogeny reconstruction is to identify SNPs in a subset of individuals, and then to compile SNPs on an array that can be used to genotype additional samples at hundreds or thousands of sites simultaneously. Although powerful and efficient, this method is subject to ascertainment bias because applying variation discovered in a representative subset to a larger sample favors identification of SNPs with high minor allele frequencies and introduces bias against rare alleles. Here, we demonstrate that the use of hybridization intensity data, rather than genotype calls, reduces the effects of ascertainment bias. Whereas traditional SNP calls assess known variants based on diversity housed in the discovery panel, hybridization intensity data survey variation in the broader sample pool, regardless of whether those variants are present in the initial SNP discovery process. We apply SNP genotype and hybridization intensity data derived from the Vitis9kSNP array developed for grape to show the effects of ascertainment bias and to reconstruct evolutionary relationships among Vitis species. We demonstrate that phylogenies constructed using hybridization intensities suffer less from the distorting effects of ascertainment bias, and are thus more accurate than phylogenies based on genotype calls. Moreover, we reconstruct the phylogeny of the genus Vitis using hybridization data, show that North American subgenus Vitis species are monophyletic, and resolve several previously poorly known relationships among North American species. This study builds on earlier work that applied the Vitis9kSNP array to evolutionary questions within Vitis vinifera and has general implications for addressing ascertainment bias in array-enabled phylogeny reconstruction.     

Genetic diversity analysis of elite European maize (Zea mays L.) inbred lines using AFLP, SSR, and SNP markers reveals ascertainment bias for a subset of SNPs

Recent advances in high-throughput sequencing technologies have triggered a shift toward single-nucleotide polymorphism (SNP) markers. A systematic bias can be introduced if SNPs are ascertained in a small panel of genotypes and then used for characterizing a larger population (ascertainment bias). With the objective of evaluating a potential ascertainment bias of the Illumina MaizeSNP50 array with respect to elite European maize dent and flint inbred lines, we compared the genetic diversity among these materials based on 731 amplified fragment length polymorphisms (AFLPs), 186 simple sequence repeats (SSRs), 41,434 SNPs of the MaizeSNP50 array (SNP-A), and two subsets of it, i.e., 30,068 Panzea (SNP-P) and 11,366 Syngenta markers (SNP-S). We evaluated the bias effects on major allele frequency, allele number, gene diversity, modified Roger’s distance (MRD), and on molecular variance (AMOVA). We revealed ascertainment bias in SNP-A, compared to AFLPs and SSRs. It affected especially European flint lines analyzed with markers (SNP-S) specifically developed to maximize differences among North American dent germplasm. The bias affected all genetic parameters, but did not substantially alter the relative distances between inbred lines within groups. For these reasons, we conclude that the SNP markers of the MaizeSNP50 array can be employed for breeding purposes in the investigated material. However, attention should be paid in case of comparisons between genotypes belonging to different heterotic groups. In this case, it is advisable to prefer a marker subset with potentially low ascertainment bias, like in our case the SNP-P marker set.     

A Phylogenetic Perspective on Sequence Evolution in Microsatellite Loci

We examined the evolution of the repeat regions of three noncoding microsatellite loci in 58 species of the Polistinae, a subfamily of wasps that diverged over 140 million years ago. A phylogenetic approach allows two new kinds of approaches to studying microsatellite evolution: character mapping and comparative analysis. The basic repeat structure of the loci was highly conserved, but was often punctuated with imperfections that appear to be phylogenetically informative. Repeat numbers evolved more rapidly than other changes in the repeat region. Changes in number of repeats among species seem consistent with the stepwise mutation model, which is based on slippage during replication as the main source of mutations. Changes in repeat numbers can occur even when there are very few tandem repeats but longer repeats, especially perfect repeats led to greater rates of evolutionary change. Species phylogenetically closer to the one from which we identified the loci had longer stretches of uninterrupted repeats and more different motifs, but not longer total repeat regions. The number of perfect repeats increased more often than it decreased. However, there was no evidence that some species have consistently greater numbers of repeats across loci than other species have, once ascertainment bias is eliminated. We also found no evidence for a population size effect posited by one form of the directionality hypothesis. Overall, phylogenetic variation in repeat regions can be explained by adding neutral evolution to what is already known about the mutation process. The life cycle of microsatellites appears to reflect a balance between growth by slippage and degradation by an essentially irreversible accumulation of imperfections. Received: 13 April 1999 / Accepted: 8 September 1999     

On the power to detect differences between male and female mutation rates for Duchenne muscular dystrophy, using classical segregation analysis and restriction fragment length polymorphisms

The power to detect departures from the theoretical proportion of new mutants in X-linked lethal disorders has been analyzed for several types of segregation analysis, including methods based on completely linked restriction fragment length polymorphisms. It is shown that all methods require large sample sizes in order to detect even large differences between male and female mutation rates. Ascertainment bias is shown to have a great effect on the outcome of the segregation analysis. All reviewed studies concerning the proportion of new mutants in Duchenne muscular dystrophy, whether they claimed equality or inequality between the male and female mutation rates, give insufficient evidence because of ascertainment bias and a too low power. An ascertainment bias-free method is given, with the advantage that information from many studies can be combined. By doing so, in the long run, even moderate departures from equality in mutation rates (if present) can be detected.     

Evolutionary Dynamics of Microsatellite Distribution in Plants: Insight from the Comparison of Sequenced Brassica,Arabidopsis and Other Angiosperm Species

Despite their ubiquity and functional importance, microsatellites have been largely ignored in comparative genomics, mostly due to the lack of genomic information. In the current study, microsatellite distribution was characterized and compared in the whole genomes and both the coding and non-coding DNA sequences of the sequenced Brassica, Arabidopsis and other angiosperm species to investigate their evolutionary dynamics in plants. The variation in the microsatellite frequencies of these angiosperm species was much smaller than those for their microsatellite numbers and genome sizes, suggesting that microsatellite frequency may be relatively stable in plants. The microsatellite frequencies of these angiosperm species were significantly negatively correlated with both their genome sizes and transposable elements contents. The pattern of microsatellite distribution may differ according to the different genomic regions (such as coding and non-coding sequences). The observed differences in many important microsatellite characteristics (especially the distribution with respect to motif length, type and repeat number) of these angiosperm species were generally accordant with their phylogenetic distance, which suggested that the evolutionary dynamics of microsatellite distribution may be generally consistent with plant divergence/evolution. Importantly, by comparing these microsatellite characteristics (especially the distribution with respect to motif type) the angiosperm species (aside from a few species) all clustered into two obviously different groups that were largely represented by monocots and dicots, suggesting a complex and generally dichotomous evolutionary pattern of microsatellite distribution in angiosperms. Polyploidy may lead to a slight increase in microsatellite frequency in the coding sequences and a significant decrease in microsatellite frequency in the whole genome/non-coding sequences, but have little effect on the microsatellite distribution with respect to motif length, type and repeat number. Interestingly, several microsatellite characteristics seemed to be constant in plant evolution, which can be well explained by the general biological rules.     

Mutational bias in penguin microsatellite DNA

Analysis of nucleotide sequence variation at a microsatellite DNA locus revealed extensive size homoplasy of alleles in Adélie penguins (Pygoscelis adeliae). Variation in the flanking regions at this locus allowed discrimination between mechanisms proposed for length changes in microsatellite DNA alleles. We further examined the structure of alleles for the same microsatellite DNA locus across 11 additional species of penguin (Spheniscidae) by mapping allele sequences onto an independent penguin phylogeny. Our analysis indicated that the repeat motifs appear to have evolved independently on several occasions. We observed sequence instability in the region bordering the repeat tract with a transversional bias predominating. We propose that this bias results from inaccurate DNA replication owing to the sequence context of this repeat tract. Because we show that regions flanking repeat sequences exhibit this mutational bias, this cautions against the use of such regions for phylogeny reconstruction.