Homoplasy and mutation model at microsatellite loci and their consequences for population genetics analysis

Homoplasy has recently attracted the attention of population geneticists, as a consequence of the popularity of highly variable stepwise mutating markers such as microsatellites. Microsatellite alleles generally refer to DNA fragments of different size (electromorphs). Electromorphs are identical in state (i.e. have identical size), but are not necessarily identical by descent due to convergent mutation(s). Homoplasy occurring at microsatellites is thus referred to as size homoplasy. Using new analytical developments and computer simulations, we first evaluate the effect of the mutation rate, the mutation model, the effective population size and the time of divergence between populations on size homoplasy at the within and between population levels. We then review the few experimental studies that used various molecular techniques to detect size homoplasious events at some microsatellite loci. The relationship between this molecularly accessible size homoplasy size and the actual amount of size homoplasy is not trivial, the former being considerably influenced by the molecular structure of microsatellite core sequences. In a third section, we show that homoplasy at microsatellite electromorphs does not represent a significant problem for many types of population genetics analyses realized by molecular ecologists, the large amount of variability at microsatellite loci often compensating for their homoplasious evolution. The situations where size homoplasy may be more problematic involve high mutation rates and large population sizes together with strong allele size constraints.     

Comparison of single-nucleotide polymorphisms and microsatellites in inference of population structure

Single-nucleotide polymorphisms (SNPs) are a class of attractive genetic markers for population genetic studies and for identifying genetic variations underlying complex traits. However, the usefulness and efficiency of SNPs in comparison to microsatellites in different scientific contexts, e.g., population structure inference or association analysis, still must be systematically evaluated through large empirical studies. In this article, we use the Collaborative Studies on Genetics of Alcoholism (COGA) data from Genetic Analysis Workshop 14 (GAW14) to compare the performance of microsatellites and SNPs in the whole human genome in the context of population structure inference. A total of 328 microsatellites and 15,840 SNPs are used to infer population structure in 236 unrelated individuals. We find that, on average, the informativeness of random microsatellites is four to twelve times that of random SNPs for various population comparisons, which is consistent with previous studies. Our results also indicate that for the combined set of microsatellites and SNPs, SNPs constitute the majority among the most informative markers and the use of these SNPs leads to better inference of population structure than the use of microsatellites. We also find that the inclusion of less informative markers may add noise and worsen the results.     

On the Structural Differences Between Markers and Genomic AC Microsatellites

AC microsatellites have proved particularly useful as genetic markers. For some purposes, such as in population biology, the inferences drawn depend on the quantitative values of their mutation rates. This, together with intrinsic biological interest, has led to widespread study of microsatellite mutational mechanisms. Now, however, inconsistencies are appearing in the results of marker-based versus non-marker-based studies of mutational mechanisms. The reasons for this have not been investigated, but one possibility, pursued here, is that the differences result from structural differences between markers and genomic microsatellites. Here we report a comparison between the CEPH AC marker microsatellites and the global population of AC microsatellites in the human genome. AC marker microsatellites are longer than the global average. Controlling for length, marker microsatellites contain on average fewer interruptions, and have longer segments, than their genomic counterparts. Related to this, marker microsatellites show a greater tendency to concentrate the majority of their repeats into one segment. These differences plausibly result from scientists selecting markers for their high polymorphism. In addition to the structural differences, there are differences in the base composition of flanking sequences, marker flanking regions being richer in C and G and poorer in A and T. Our results indicate that there are profound differences between marker and genomic microsatellites that almost certainly affect their mutation rates. There is a need for a unified model of mutational mechanisms that accounts for both marker-derived and genomic observations. A suggestion is made as to how this might be done.Reviewing Editor: Dr. Magnto Nordborg     

Conserved flanking microsatellite sequences (ReFS) differentiate between Lepidoptera species,and provide insight into microsatellite evolution

High rates of mutation and homoplasy mean that microsatellites generally are not considered to be useful molecular markers for inferring systematic relationships between species. However, an earlier pilot study suggested that conserved flanking microsatellite sequences, also known as repetitive flanking sequences (ReFS), may form a basis for a dominant marker that can differentiate between species of Lepidoptera. We present data that demonstrate that ReFS are quick and easy to use, and generate highly repeatable banding patterns from a range of Lepidoptera species. Sequence data from a subset of ReFS‐amplified bands revealed microsatellite families with flanking sequences that are more conserved within than among species: this is probably attributable to recombination‐mediated events, transposition of mobile elements or a combination of the two. Our data support the use of ReFS as dominant interspecific molecular markers, and add to the growing literature on the evolution of microsatellites in Lepidoptera.     

The utility of fast evolving molecular markers for studying speciation in the Antarctic benthos

The Southern Ocean is surprisingly rich in species that coexist in one of the most extreme environments on Earth yet the processes leading to speciation in this ecosystem are not well understood. To remedy this, tools that measure the genetic connectedness within a species are needed. Although useful for phylogenetic purposes, the readily available mitochondrial markers (e.g. 16S, COI) suffer from numerous shortcomings for population genetics. Therefore, molecular markers are needed that are sufficiently variable, unlinked, biparentally inherited, and distributed over the whole genome. We argue that microsatellites are suitable markers that have not been widely used in exploratory studies due to their difficult initial set-up. Working with the Ceratoserolis trilobitoides species complex (Isopoda), we demonstrate that using a novel protocol many microsatellites can be identified quickly. An increased availability of these highly sensitive markers will be useful for studies addressing the origin of species in the Southern Ocean and their response to future climate change. Christoph Held and Florian Leese contributed equally to this paper.     

Challenges in analysis and interpretation of microsatellite data for population genetic studies

Advancing technologies have facilitated the ever‐widening application of genetic markers such as microsatellites into new systems and research questions in biology. In light of the data and experience accumulated from several years of using microsatellites, we present here a literature review that synthesizes the limitations of microsatellites in population genetic studies. With a focus on population structure, we review the widely used fixation (F_ST) statistics and Bayesian clustering algorithms and find that the former can be confusing and problematic for microsatellites and that the latter may be confounded by complex population models and lack power in certain cases. Clustering, multivariate analyses, and diversity‐based statistics are increasingly being applied to infer population structure, but in some instances these methods lack formalization with microsatellites. Migration‐specific methods perform well only under narrow constraints. We also examine the use of microsatellites for inferring effective population size, changes in population size, and deeper demographic history, and find that these methods are untested and/or highly context‐dependent. Overall, each method possesses important weaknesses for use with microsatellites, and there are significant constraints on inferences commonly made using microsatellite markers in the areas of population structure, admixture, and effective population size. To ameliorate and better understand these constraints, researchers are encouraged to analyze simulated datasets both prior to and following data collection and analysis, the latter of which is formalized within the approximate Bayesian computation framework. We also examine trends in the literature and show that microsatellites continue to be widely used, especially in non‐human subject areas. This review assists with study design and molecular marker selection, facilitates sound interpretation of microsatellite data while fostering respect for their practical limitations, and identifies lessons that could be applied toward emerging markers and high‐throughput technologies in population genetics.     

Genetic diversity and differentiation in Eryngium alpinum L. (Apiaceae): comparison of AFLP and microsatellite markers

Genetic diversity and structure of 12 populations of Eryngium alpinum L. were investigated using 63 dominant amplified fragment length polymorphism (AFLP) and seven codominant microsatellite (48 alleles) markers. Within-population diversity estimates obtained with both markers were not correlated, but the microsatellite-based fixation index Fis was correlated with both AFLP diversity indices (number of polymorphic bands and Nei's expected heterozygosity). Only AFLP diversity indices increased with the size of populations, although they did not significantly differ among them (Kruskall-Wallis test). The discrepancy between AFLPs and microsatellites may be explained by a better coverage of the genome with numerous AFLPs, the higher mutation rates of microsatellites or the absence of significant difference among within-population diversity estimates. Genetic differentiation was higher with AFLPs (theta=0.40) than with microsatellites (theta=0.23), probably due to the higher polymorphism of microsatellites. Thus, we considered global qualitative patterns rather than absolute estimates to compare the performance of both types of markers. On a large geographic scale, the Mantel test and multivariate analysis showed that genetic patterns were more congruent with the spatial arrangement of populations when inferred from microsatellites than from AFLPs, suggesting higher homoplasy of AFLP markers. On a small spatial scale, AFLPs managed to discriminate individuals from neighboring populations whereas microsatellites did not (multivariate analysis), and the percentage of individuals correctly assigned to their population of origin was higher with AFLPs than with microsatellites. However, dominant AFLPs cannot be used to study heterozygosity-related topics. Thus, distinct molecular markers should be used depending on the biological question and the geographical scale investigated.     

Isolation and characterization of dinucleotide repeat microsatellites in Drosophila ananassae

Drosophila ananassae is a cosmopolitan species with a geographic range throughout most of the tropical and subtropical regions of the world. Previous studies of DNA sequence polymorphism in three genes has shown evidence of selection affecting broad expanses of the genome in regions with low rates of recombination in geographically local populations in and around India. The studies suggest that extensive physical and genetic maps based on molecular markers, and detailed studies of population structure may provide insight into the degree to which natural selection affects DNA sequence polymorphism across broad regions of chromosomes. We have isolated 85 dinucleotide repeat microsatellite sequences and developed assay conditions for genotyping using PCR. The dinucleotide repeats we isolated are shorter, on average, than those isolated in many other Drosophila species. Levels of genetic variation are high, comparable to Drosophila melanogaster. The levels of variation indicate the effective population size of an Indonesian population of D. ananassae is 58,692 (infinite allele model) and 217,284 (stepwise mutation model), similar to estimates of effective population size for D. melanogaster calculated using dinucleotide repeat microsatellites. The data also show that the Indonesian population is in a rapid expansion phase. Cross-species amplification of the microsatellites in 11 species from the Ananassae, Elegans, Eugracilis and Ficusphila subgroups indicates that the loci may be useful for studies of the sister species, D. pallidosa, but will have limited use for more distantly related species.     

Mutational Dynamics of Microsatellites

Microsatellites are a ubiquitous class of simple repetitive DNA sequences, which are widespread in both eukaryotic and prokaryotic genomes. The use of microsatellites as polymorphic DNA markers has considerably increased both in the number of studies and in the number of organisms, primarily for genetic mapping, studying genomic instability in cancer, population genetics, forensics, conservation biology, molecular anthropology and in the studies of human evolutionary history. Although simple sequence repeats have been extensively used in studies encompassing varied areas of genetics, the mutation dynamics of these genome regions is still not well understood. The present review focuses on the mutational dynamics of microsatellite DNA with special reference to mutational mechanisms and their role in microsatellite evolution.     

Heterosis, marker mutational processes and population inbreeding history.

Genotype-fitness correlations (GFC) have previously been studied using allozyme markers and have often focused on short-term processes such as recent inbreeding. Thus, models of GFC usually neglect marker mutation and only use heterozygosity as a genotypic index. Recently, GFC have also been reported (i) with DNA markers such as microsatellites, characterized by high mutation rates and specific mutational processes and (ii) using new individual genotypic indices assumed to be more precise than heterozygosity. The aim of this article is to evaluate the theoretical impact of marker mutation on GFC. We model GFC due to short-term processes generated by the current breeding system (partial selfing) and to long-term processes generated by past population history (hybridization). Various mutation rates and mutation models corresponding to different kinds of molecular markers are considered. Heterozygosity is compared to other genotypic indices designed for specific marker types. Highly mutable markers (such as microsatellites) are particularly suitable for the detection of GFC that evolve in relation to short-term processes, whereas GFC due to long-term processes are best observed with intermediate mutation rates. Irrespective of the marker type and population scenario, heterozygosity usually provides higher correlations than other genotypic indices under most biologically plausible conditions.     

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.     

An integrated chromosome map of microsatellite markers and inversion breakpoints for an Asian malaria mosquito, Anopheles stephensi

Anopheles stephensi is one of the major vectors of malaria in the Middle East and Indo-Pakistan subcontinent. Understanding the population genetic structure of malaria mosquitoes is important for developing adequate and successful vector control strategies. Commonly used markers for inferring anopheline taxonomic and population status include microsatellites and chromosomal inversions. Knowledge about chromosomal locations of microsatellite markers with respect to polymorphic inversions could be useful for better understanding a genetic structure of natural populations. However, fragments with microsatellites used in population genetic studies are usually too short for successful labeling and hybridization with chromosomes. We designed new primers for amplification of microsatellite loci identified in the A. stephensi genome sequenced with next-generation technologies. Twelve microsatellites were mapped to polytene chromosomes from ovarian nurse cells of A. stephensi using fluorescent in situ hybridization. All microsatellites hybridized to unique locations on autosomes, and 7 of them localized to the largest arm 2R. Ten microsatellites were mapped inside the previously described polymorphic chromosomal inversions, including 4 loci located inside the widespread inversion 2Rb. We analyzed microsatellite-based population genetic data available for A. stephensi in light of our mapping results. This study demonstrates that the chromosomal position of microsatellites may affect estimates of population genetic parameters and highlights the importance of developing physical maps for nonmodel organisms.     

Waste not,want not: Microsatellites remain an economical and informative technology for conservation genetics

Comparisons of microsatellites and single‐nucleotide polymorphisms (SNPs) have found that SNPs outperform microsatellites in population genetic analyses, questioning the continued utility of microsatellites in population and landscape genetics. Yet, highly polymorphic markers may be of value in species that have reduced genetic variation. This study repeated previous analyses that used microsatellites with SNPs developed from ddRAD sequencing in the black‐capped vireo source‐sink system. SNPs provided greater resolution of genetic diversity, population differentiation, and migrant detection but could not reconstruct parentage relationships due to insufficient heterozygosities. The biological inferences made by both sets of markers were similar: asymmetrical gene flow from source sites to the remaining sink sites. With the landscape genetic analyses, we found different results between the two molecular markers, but associations of the top environmental features (riparian, open habitat, agriculture, and human development) with dispersal estimates were shared between marker types. Despite the higher precision of SNPs, we find that microsatellites effectively uncover population processes and patterns and are superior for parentage analyses in this species with reduced genetic diversity. This study illustrates the continued applicability and relevance of microsatellites in population genetic research.     

Conservation and evolution of microsatellite loci in primate taxa

Microsatellites are promising genetic markers for the study of demographic structure and phylogenetic history in populations. However, little information exists on the molecular nature of the repeats and their flanking sequences of a same microsatellite in a large range of species. In this study, we report polymorphism and consensus sequences of eight microsatellite loci using human primers in 20 primate species. The results show size polymorphism in almost all species and microsatellites. These loci are therefore useful markers for population genetic studies between populations of the same species. Insertion/deletion events are frequent in the flanking regions, the majority concerning several contiguous bases. This is in contrast with the more usual single base pair events in non-coding regions. The ranges of allele lengths in non-human primates often show no overlap with that of human, usually due to the deletion/insertion events in the flanking sequences, producing smaller allele lengths rather than smaller numbers of repeats. The use of length of PCR product will bias the inter-species interpretation reducing the number of observable alleles and treating as the same allele very divergent molecular sequences. Caution should be used when employing microsatellites in cross-species comparisons in which the species under study are separated by significant amounts of evolutionary time: in such cases allele comparison cannot be based on lengths alone.     

G'ST and D do not replace FST

The genetic differentiation among populations is affected by mutation as well as by migration, drift and selection. For loci with high mutation rates, such as microsatellites, the amount of mutation can influence the values of indices of differentiation such as G(ST) and F(ST). For many purposes, this effect is undesirable, and as a result, new indices such as G'(ST) and D have been proposed to measure population differentiation. This paper shows that these new indices are not effective measures of the causes or consequences of population structure. Both G'(ST) and D depend heavily on mutation rate, but both are insensitive to any population genetic process when the mutation rate is high relative to the migration rate. Furthermore, D is specific to the locus being measured, and so little can be inferred about the population demography from D. However, at equilibrium, D may provide an index of whether a particular marker is more strongly affected by mutation than by migration. I argue that F(ST) is a more important summary of the effects of population structure than D and that R(ST) or other measures that explicitly account for the mutation process are much better than G(ST), G'(ST), or D for highly mutable markers. Markers with lower mutation rates will often be easier to interpret.     

Frequent non-reciprocal exchange in microsatellite-containing-DNA-regions of vertebrates

Microsatellites are DNA-fragments containing short repetitive motifs with 2–10 bp. They are highly variable in most species and distributed throughout the whole genome. It is broadly accepted that their high degree of variability is closely associated with mispairing of DNA-strands during the replication phase, termed slippage, although recombination is also observed. The aim of this study is to demonstrate evidence that non-reciprocal recombination processes changing the total genomic structure are common in microsatellites and flanking regions. We sequenced DNA fragments from birds in which microsatellites are located, and analyzed the structure of the microsatellites and their flanking regions. Additionally, other data and those from literature of three microsatellite regions of primates coding for the Ataxin-2, the Huntingtin and the TATA-box binding protein were analyzed. The structures of seven avian and three primate microsatellites support the hypothesis that non-reciprocal recombination is a common process that may also contribute considerably to the variation at microsatellite loci. We conclude that results of population genetic studies that are analyzed statistically with methods based on stepwise mutation models should be interpreted with caution if no detailed information on the allelic variation of microsatellites is available.     

A comparison of RAPD versus microsatellite DNA markers in population studies of the massasauga rattlesnake

We compared genetic differentiation among populations of the threatened massasauga rattlesnake (Sistrurus c. catenatus) using two types of nuclear molecular markers: randomly amplified polymorphic DNA (RAPD) markers and microsatellites. Analyses of molecular variance (AMOVA) and G(ST) and F(ST) analyses indicated that levels of among-population differentiation between regional populations (>100 km) were comparable for both markers. However, microsatellites were superior in population assignment tests and at discerning fine-scale genetic differentiation between subpopulations separated by tens of kilometers. These results argue that both types of markers are suitable for defining broad-scale genetic structures in snake populations and can provide important inputs into conservation initiatives of focal taxa. However, our analyses suggest that microsatellites 3re better for detecting structure at limited spatial scales.     

Differentiating salmon populations at broad and fine geographical scales with microsatellites and single nucleotide polymorphisms

Single nucleotide polymorphisms (SNPs) are appealing genetic markers due to several beneficial attributes, but uncertainty remains about how many of these bi-allelic markers are necessary to have sufficient power to differentiate populations, a task now generally accomplished with highly polymorphic microsatellite markers. In this study, we tested the utility of 37 SNPs and 13 microsatellites for differentiating 29 broadly distributed populations of Chinook salmon ( n  = 2783). Information content of all loci was determined by In and     , and the top 12 markers ranked by In were microsatellites, but the 6 highest, and 7 of the top 10     ranked markers, were SNPs. The mean ratio of random SNPs to random microsatellites ranged from 3.9 to 4.1, but this ratio was consistently reduced when only the most informative loci were included. Individual assignment test accuracy was higher for microsatellites (73.1%) than SNPs (66.6%), and pooling all 50 markers provided the highest accuracy (83.2%). When marker types were combined, as few as 15 of the top ranked loci provided higher assignment accuracy than either microsatellites or SNPs alone. Neighbour-joining dendrograms revealed similar clustering patterns and pairwise tests of population differentiation had nearly identical results with each suite of markers. Statistical tests and simulations indicated that closely related populations were better differentiated by microsatellites than SNPs. Our results indicate that both types of markers are likely to be useful in population genetics studies and that, in some cases, a combination of SNPs and microsatellites may be the most effective suite of loci.     

A test for concordance between the multilocus genealogies of genes and microsatellites in the pathogenic fungus Coccidioides immitis

Uncovering the correct phylogeny of closely related species requires analysis of multiple gene genealogies or, alternatively, genealogies inferred from the multiple alleles found at highly polymorphic loci, such as microsatellites. However, a concern in using microsatellites is that constraints on allele sizes may occur, resulting in homoplasious distributions of alleles, leading to incorrect phylogenies. Seven microsatellites from the pathogenic fungus Coccidioides immitis were sequenced for 20 clinical isolates chosen to represent the known genetic diversity of the pathogen. An organismal phylogeny for C. immitis was inferred from microsatellite-flanking sequence polymorphisms and other restriction fragment length polymorphism-containing loci. Two microsatellite genetic distances were then used to determine phylogenies for C. immitis, and the trees found by these three methods were compared. Congruence between the organismal and microsatellite phylogenies occurred when microsatellite distances were based on simple allele frequency data. However, complex mutation events at some loci made distances based on stepwise mutation models unreliable. Estimates of times of divergence for the two species of C. immitis based on microsatellites were significantly lower than those calculated from flanking sequence, most likely due to constraints on microsatellite allele sizes. Flanking-sequence insertions/deletions significantly decreased the accuracy of genealogical information inferred from microsatellite loci and caused interspecific length homoplasies at one of the seven loci. Our analysis shows that microsatellites are useful phylogenetic markers, although care should be taken to choose loci with appropriate flanking sequences when they are intended for use in evolutionary studies.     

Long tomato microsatellites are predominantly associated with centromeric regions.

Microsatellites as genetic markers are used in many crop plants. Major criteria for their usability as molecular markers include that they are highly polymorphic and evenly spread throughout a genome. In tomato, it has been reported that long arrays of tetranucleotide microsatellites containing the motif GATA are highly clustered around the centromeres of all chromosomes. In this study, we have isolated tomato microsatellites containing long arrays (> 20 repeats) of the dinucleotide motifs GA, GT, AT, as well as GATA, assessed their variability within Lycopersicon esculentum varieties and mapped them onto a genetic map of tomato. The investigated microsatellite markers exhibited between 1 and 5 alleles in a diverse set of L. esculentum lines. Mapping of the microsatellites onto the genetic map of tomato demonstrates that, as previously shown, GATA microsatellites are highly clustered in the regions of the tomato centromeres. Interestingly, the same centromeric location was now found for long dinucleotide microsatellite markers. Because of this uneven distribution, genetic mapping of the entire tomato genome using long dinucleotide microsatellites will be very difficult to achieve and microsatellite markers with shorter arrays of microsatellites could be more suitable for mapping experiments albeit their lower level of polymorphism. Some microsatellite markers described in this study might provide a useful tool to study the molecular structure of tomato centromeric regions and for variety identification.