Rise of the machines--recommendations for ecologists when using next generation sequencing for microsatellite development

Next generation sequencing is revolutionizing molecular ecology by simplifying the development of molecular genetic markers, including microsatellites. Here, we summarize the results of the large-scale development of microsatellites for 54 nonmodel species using next generation sequencing and show that there are clear differences amongst plants, invertebrates and vertebrates for the number and proportion of motif types recovered that are able to be utilized as markers. We highlight that the heterogeneity within each group is very large. Despite this variation, we provide an indication of what number of sequences and consequent proportion of a 454 run are required for the development of 40 designable, unique microsatellite loci for a typical molecular ecological study. Finally, to address the challenges of choosing loci from the vast array of microsatellite loci typically available from partial genome runs (average for this study, 2341 loci), we provide a microsatellite development flowchart as a procedural guide for application once the results of a partial genome run are obtained.     

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 Interruptions Stabilize Primate Genomes and Exist as Population-Specific Single Nucleotide Polymorphisms within Individual Human Genomes

Interruptions of microsatellite sequences impact genome evolution and can alter disease manifestation. However, human polymorphism levels at interrupted microsatellites (iMSs) are not known at a genome-wide scale, and the pathways for gaining interruptions are poorly understood. Using the 1000 Genomes Phase-1 variant call set, we interrogated mono-, di-, tri-, and tetranucleotide repeats up to 10 units in length. We detected ∼26,000–40,000 iMSs within each of four human population groups (African, European, East Asian, and American). We identified population-specific iMSs within exonic regions, and discovered that known disease-associated iMSs contain alleles present at differing frequencies among the populations. By analyzing longer microsatellites in primate genomes, we demonstrate that single interruptions result in a genome-wide average two- to six-fold reduction in microsatellite mutability, as compared with perfect microsatellites. Centrally located interruptions lowered mutability dramatically, by two to three orders of magnitude. Using a biochemical approach, we tested directly whether the mutability of a specific iMS is lower because of decreased DNA polymerase strand slippage errors. Modeling the adenomatous polyposis coli tumor suppressor gene sequence, we observed that a single base substitution interruption reduced strand slippage error rates five- to 50-fold, relative to a perfect repeat, during synthesis by DNA polymerases α, β, or η. Computationally, we demonstrate that iMSs arise primarily by base substitution mutations within individual human genomes. Our biochemical survey of human DNA polymerase α, β, δ, κ, and η error rates within certain microsatellites suggests that interruptions are created most frequently by low fidelity polymerases. Our combined computational and biochemical results demonstrate that iMSs are abundant in human genomes and are sources of population-specific genetic variation that may affect genome stability. The genome-wide identification of iMSs in human populations presented here has important implications for current models describing the impact of microsatellite polymorphisms on gene expression.     

In silico whole-genome EST analysis reveals 2322 novel microsatellites for the silver-lipped pearl oyster, Pinctada maxima

Molecular stock improvement techniques such as marker assisted selection have great potential in accelerating selective breeding programmes for animal production industries. However, the discovery and application of trait/marker associations usually requires a large number of genome-wide polymorphic loci. Here, we present 2322 unique microsatellites for the silver-lipped pearl oyster, Pinctada maxima, a species of aquaculture importance throughout the Indo-Australian Archipelago for production of the highly valued South Sea pearl. More than 1.2 million Roche 454 expressed sequence tag (EST) reads were screened for microsatellite repeat motifs. A total of 12,604 sequences contained either a di, tri, tetra, penta or hexa microsatellite repeat motif (n ≥ 6), with 6435 of these sequences having sufficient flanking regions for primer development. All identified microsatellites with designed primers were condensed into 2322 unique clusters (i.e., unique loci) of which 360 were shown to be polymorphic based on multiple sequence reads with different repeat motifs. Genotyping of five microsatellite loci demonstrated that in silico evaluation of polymorphism levels was a very useful method for identification of polymorphic loci, with the variation uncovered being a lower bound. Gene Ontology annotations of sequences containing microsatellites suggest that most are derived from a diverse array of unique genes. This EST derived microsatellite database will be a valuable resource for future studies in genetic map construction, diversity analysis, quantitative trait loci analysis, association mapping and marker assisted selection, not only for P. maxima, but also closely related species within the genus Pinctada.     

MICROSATELLITE EVOLUTION IN VERTEBRATES: INFERENCE FROM AC DINUCLEOTIDE REPEATS

Abstract We analyze published data from 592 AC microsatellite loci from 98 species in five vertebrate classes including fish, reptiles, amphibians, birds, and mammals. We use these data to address nine major questions about microsatellite evolution. First, we find that larger genomes do not have more microsatellite loci and therefore reject the hypothesis that microsatellites function primarily to package DNA into chromosomes. Second, we confirm that microsatellite loci are relatively rare in avian genomes, but reject the hypothesis that this is due to physical constraints imposed by flight. Third, we find that microsatellite variation differs among species within classes, possibly relating to population dynamics. Fourth, we reject the hypothesis that microsatellite structure (length, number of alleles, allele dispersion, range in allele sizes) differs between poikilotherms and homeotherms. The difference is found only in fish, which have longer microsatellites and more alleles than the other classes. Fifth, we find that the range in microsatellite allele size at a locus is largely due to the number of alleles and secondarily to allele dispersion. Sixth, length is a major factor influencing mutation rate. Seventh, there is a directional mutation toward an increase in microsatellite length. Eighth, at the species level, microsatellite and allozyme heterozygosity covary and therefore inferences based on large-scale studies of allozyme variation may also reflect microsatellite genetic diversity. Finally, published microsatellite loci (isolated using conventional hybridization methods) provide a biased estimate of the actual mean repeat length of microsatellites in the genome.     

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.     

Negative covariance suggests mutation bias in a two-locus microsatellite system in the fish Sparus aurata.

Constraints on microsatellite length appear to vary in a species-specific manner. We know very little about the nature of these constraints and why they should vary among species. While surveying microsatellite variation in the Mediterranean gilthead sea bream, Sparus aurata, we discovered an unusual pattern of covariation between two closely linked microsatellite loci. One- and two-locus haplotypes were scored from PCR amplification products of each locus separately and both loci together. In a sample of 211 fish, there was a strong negative covariance in repeat number between the two loci, which suggests a mechanism that maintains the combined length below a constrained size. In addition, there were two clusters of the same combined haplotype length, one consisting of a long repeat array at one locus and a short array at the other and vice versa. We demonstrate that several models of biased mutation or natural selection, in theory, could generate this pattern of covariance. The common feature of all the models is the idea that tightly linked microsatellites do not evolve in complete independence, and that whatever size dependence there is to the process, it appears to "read" the combined size of the two loci.     

Development of a multiplex PCR assay for fine-scale population genetic analysis of the Komodo monitor Varanus komodoensis based on 18 polymorphic microsatellite loci

Multiplex PCR assays for the coamplification of microsatellite loci allow rapid and cost-effective genetic analyses and the production of efficient screening protocols for international breeding programs. We constructed a partial genomic library enriched for di-nucleotide repeats and characterized 14 new microsatellite loci for the Komodo monitor (or Komodo dragon, Varanus komodoensis). Using these novel microsatellites and four previously described loci, we developed multiplex PCR assays that may be loaded on a genetic analyser in three separate panels. We tested the novel set of microsatellites for polymorphism using 69 individuals from three island populations and evaluated the resolving power of the entire panel of 18 loci by conducting (i) a preliminary assignment test to determine population(s) of origin and (ii) a parentage analysis for 43 captive Komodo monitors. This panel of polymorphic loci proved useful for both purposes and thus can be exploited for fine-scale population genetic analyses and as part of international captive breeding programs directed at maintaining genetically viable ex situ populations and reintroductions.     

Causes of Size Homoplasy Among Chloroplast Microsatellites in Closely Related <Emphasis Type="Italic">Clusia</Emphasis> Species

Chloroplast DNA sequences and microsatellites are useful tools for phylogenetic as well as population genetic analyses of plants. Chloroplast microsatellites tend to be less variable than nuclear microsatellites and therefore they may not be as powerful as nuclear microsatellites for within-species population analysis. However, chloroplast microsatellites may be useful for phylogenetic analysis between closely related taxa when more conventional loci, such as ITS or chloroplast sequence data, are not variable enough to resolve phylogenetic relationships in all clades. To determine the limits of chloroplast microsatellites as tools in phylogenetic analyses, we need to understand their evolution. Thus, we examined and compared phylogenetic relationships of species within the genus Clusia, using both chloroplast sequence data and variation at seven chloroplast microsatellite loci. Neither ITS nor chloroplast sequences were variable enough to resolve relationships within some sections of the genus, yet chloroplast microsatellite loci were too variable to provide any useful phylogenetic information. Size homoplasy was apparent, caused by base substitutions within the microsatellite, base substitutions in the flanking regions, indels in the flanking regions, multiple microsatellites within a fragment, and forward/reverse mutations of repeat length resulting in microsatellites of identical base composition that were not identical by descent.     

Evidence for complex mutations at microsatellite loci in Drosophila.

Fifteen lines each of Drosophila melanogaster, D. simulans, and D. sechellia were scored for 19 microsatellite loci. One to four alleles of each locus in each species were sequenced, and microsatellite variability was compared with sequence structure. Only 7 loci had their size variation among species consistent with the occurrence of strictly stepwise mutations in the repeat array, the others showing extensive variability in the flanking region compared to that within the microsatellite itself. Polymorphisms apparently resulting from complex nonstepwise mutations involving the microsatellite were also observed, both within and between species. Maximum number of perfect repeats and variance of repeat count were found to be strongly correlated in microsatellites showing an apparently stepwise mutation pattern. These data indicate that many microsatellite mutation events are more complex than represented even by generalized stepwise mutation models. Care should therefore be taken in inferring population or phylogenetic relationships from microsatellite size data alone. The analysis also indicates, however, that evaluation of sequence structure may allow selection of microsatellites that more closely match the assumptions of stepwise models.     

Evolution of microsatellites in Arabis petraea and Arabis lyrata, outcrossing relatives of Arabidopsis thaliana

We examined microsatellite variation in two diploid, outcrossing relatives of Arabidopsis thaliana, Arabis petraea and Arabis lyrata. The primer sequences were derived from A. thaliana. About 50% (14 loci) of the A. thaliana primers could successfully amplify microsatellites in the related species. Analysis of microsatellite structure in the related species showed that there had been large changes in the microsatellites: there were large differences in repeat numbers and many of the A. thaliana simple repeats were shorter in the related species. For the loci we compared, the related species had a much lower level of variability at the microsatellites than Japanese wild populations of A. thaliana. This is presumably related to the different microsatellite structures, because allozyme data showed that the outcrossing relatives were highly polymorphic compared to other outcrossing herbaceous species. Use of microsatellites in assessing variability or phylogenetic relationships between different species requires caution, because changes in microsatellite structure may alter evolutionary rates.      

Microsatellite evolution--a reciprocal study of repeat lengths at homologous loci in cattle and sheep

The application of microsatellites in evolutionary studies requires an understanding of the patterns governing their evolution in different species. The finding that homologous microsatellite loci are longer, i.e., containing more repeat units, in human and in other primates has been taken as evidence for directional microsatellite evolution and for a difference in the rate of evolution between species. However, it has been argued that this finding is an inevitable consequence of biased selection of longer-than-average microsatellites in human, because cloning procedures are adopted to generate polymorphic and, hence, long markers. As a test of this hypothesis, we conducted a reciprocal comparison of the lengths of microsatellite loci in cattle and sheep using markers derived from the bovine genome as well as the ovine genome. In both cases, amplification products were longer in the focal species, and loci were also more polymorphic in the species from which they were originally cloned. The crossing pattern that we found suggests that interspecific length differences detected at homologous microsatellite loci are the result of biased selection of loci associated with cloning procedures. Hence, comparisons of microsatellite evolution between species are flawed unless they are based on reciprocal analyses or on genuinely random selection of loci with respect to repeat length.      

Heterogeneous Nature and Distribution of Interruptions in Dinucleotides May Indicate the Existence of Biased Substitutions Underlying Microsatellite Evolution

Some aspects of microsatellite evolution, such as the role of base substitutions, are far from being fully understood. To examine the significance of base substitutions underlying the evolution of microsatellites we explored the nature and the distribution of interruptions in dinucleotide repeats from the human genome. The frequencies that we inferred in the repetitive sequences were statistically different from the frequencies observed in other noncoding sequences. Additionally, we detected that the interruptions tended to be towards the ends of the microsatellites and 5'-3' asymmetry. In all the estimates nucleotides forming the same repetitive motif seem to be affected by different base substitution rates in AC and AG. This tendency itself could generate patterning and similarity in flanking sequences and reconcile these phenomena with the high mutation rate found in flanking sequences without invoking convergent evolution. Nevertheless, our data suggest that there is a regional bias in the substitution pattern of microsatellites. The accumulation of random substitutions alone cannot explain the heterogeneity and the asymmetry of interruptions found in this study or the relative frequency of different compound microsatellites in the human genome. Therefore, we cannot rule out the possibility of a mutational bias leading to convergent or parallel evolution in flanking sequences.     

Characterization of AT-rich microsatellites in common bean (<Emphasis Type="Italic">Phaseolus vulgaris</Emphasis> L.)

Polymorphism of microsatellite markers is often associated with the simple sequence repeat motif targeted. AT-rich microsatellites tend to be highly variable and this appears to be notable, especially in legume genomes. To analyze the value of AT-rich microsatellites for common bean (Phaseolus vulgaris L.), we developed a total of 85 new microsatellite markers, 74 of which targeted ATA or other AT-rich motif loci and 11 of which were made for GA, CA or CAC motif loci. We evaluated the loci for the level of allelic diversity in comparison to previously characterized microsatellites using a panel of 18 standard genotypes and genetically mapped any loci polymorphic in the DOR364 × G19833 population. The majority of the microsatellites produced single bands and detected single loci, however, 15 of the AT-rich microsatellites produced multiple or double banding patterns; while only one of the GA or CA-rich microsatellites did. The polymorphism information content (PIC) values averaged 0.892 and 0.600 for the AT and ATA motif microsatellites, respectively, but only 0.140 for the CA-rich microsatellites. GA microsatellites, which had a large average number of repeats, had high to intermediate PIC, averaging 0.706. A total of 45 loci could be genetically mapped and distribution of the loci across the genome was skewed towards non-distal locations with a greater prevalence of loci on linkage groups b02, b09 and b11. AT-rich microsatellites were found to be a useful source of polymorphic markers for mapping and diversity assessment in common bean that appears to uncover higher diversity than other types of simple sequence repeat markers.     

Comparative allozyme and microsatellite population structure in a narrow endemic plant species, Centaurea corymbosa Pourret (Asteraceae)

Centaurea corymbosa Pourret (Asteraceae) is a narrow endemic species known only from six populations located in a 3-km2 area in the south of France. Earlier field experiments have suggested that pollen and seed dispersal were highly restricted within and among populations. Consistent with the field results, populations were highly differentiated for five allozyme loci and among-population variation fitted an isolation-by-distance model. In the present study, we investigated the genetic structure of C. corymbosa using six microsatellite loci. As with allozymes, microsatellites revealed no within-population structure and a large differentiation among populations. However, allozyme loci were less powerful than microsatellites in detecting the extent of gene flow assessed by assignment tests. The patterns of structuration greatly varied among loci for both types of marker; we suggest that differences in single-locus pattern could mainly be an effect of stochastic variation for allozymes and an effect of variation in mutation rate for microsatellites. In contrast to the multilocus results, the two most polymorphic microsatellite loci did not show any isolation-by-distance pattern. Our results suggest that highly variable loci might not always be the best suited markers to quantify levels of gene flow among populations.     

High degree of conservation of nuclear microsatellite loci in the genus Clusia.

In plants, microsatellites and their flanking DNA are rarely conserved across a whole genus, let alone other genera in the same family. Therefore, the possibility of using microsatellite primers developed for a species across a large number of plant species in the same genus is often limited. Remarkably, dinucleotide nuclear microsatellites developed for Clusia minor and for Clusia nemorosa amplified homologous microsatellites in species across the whole genus Clusia. In this present study, we report on the DNA sequence variation across the genus of 3 microsatellite loci with varying levels of variation. Compared over the species, there was a correlation between the lengths of the microsatellite loci. Interrupts occurred multiple times and did not seem to lead to the death of the microsatellite. These highly conserved markers will be useful for studying the variable reproductive systems in the genus Clusia.     

Size homoplasy and mutational processes of interrupted microsatellites in two bee species, Apis mellifera and Bombus terrestris (Apidae)

Similar microsatellite electromorphs (PCR products of the same size) can arise from independent mutational events. Such alleles are not identical by descent. This phenomenon, termed size homoplasy, was studied by sequencing electromorphs of two microsatellite loci in which the stretch of basic repeats is interrupted by different short (1-2 bp) DNA motifs. The number and position of these interruptions were established for electromorphs from closely and distantly related populations of honeybees and bumblebees. No sequence difference was found when electromorphs came from the same subspecies or from closely related subspecies, suggesting that they were probably identical by descent. In contrast, sequence differences were often detected in distantly related subspecies, showing that size homoplasy frequently occurs at this level of population differentiation. Size homoplasy is increased by limits to free length variation of alleles, a phenomenon that seems to act on interrupted microsatellites when comparing distantly related taxa, that is, honeybee subspecies from different evolutionary lineages. Electromorph sequences suggest that, within the scope of these limits, large mutation events have occurred frequently at both interrupted loci studied. In good agreement with the molecular data, computations based on the observed heterozygosity and number of electromorphs and simulation studies showed that neither locus fits the one-step stepwise mutant model (SMM). We speculate that interrupted microsatellites in general could be characterized by a higher variance in repeat number and consequently a lower homoplasy rate than pure ones. Hence, interrupted microsatellites should be most appropriate for investigating population differentiation and evolutionary relationship between relatively distant populations.      

A simple method for isolation of microsatellites from the Japanese squirrel,Sciurus lis,without constructing a genomic library

We developed a simple and easy method to isolate microsatellites without screening genomic libraries by hybridization. The method requires only three basic techniques: polymerase chain reaction, DNA cloning and sequencing. We applied this method to develop microsatellite markers for the Japanese squirrel and isolated 45 clones that contained repetitive sequences. Among the 22 clones that we tested further, we found 11 diagnostic microsatellite loci that are applicable to the molecular ecological study of Japanese squirrels.     

Conservation of microsatellites in three species of salmonids

The conservation of flanking regions of ten microsatellites has been investigated in three salmonid species. The proportion of conserved microsatellites and the level of polymorphism are high in the three species. This confirms the potential interest of microsatellite loci for interspecific comparisons.     

Design and implementation of degenerate microsatellite primers for the mammalian clade

Microsatellites are popular genetic markers in molecular ecology, genetic mapping and forensics. Unfortunately, despite recent advances, the isolation of de novo polymorphic microsatellite loci often requires expensive and intensive groundwork. Primers developed for a focal species are commonly tested in a related, non-focal species of interest for the amplification of orthologous polymorphic loci; when successful, this approach significantly reduces cost and time of microsatellite development. However, transferability of polymorphic microsatellite loci decreases rapidly with increasing evolutionary distance, and this approach has shown its limits. Whole genome sequences represent an under-exploited resource to develop cross-species primers for microsatellites. Here we describe a three-step method that combines a novel in silico pipeline that we use to (1) identify conserved microsatellite loci from a multiple genome alignments, (2) design degenerate primer pairs, with (3) a simple PCR protocol used to implement these primers across species. Using this approach we developed a set of primers for the mammalian clade. We found 126,306 human microsatellites conserved in mammalian aligned sequences, and isolated 5,596 loci using criteria based on wide conservation. From a random subset of ~1000 dinucleotide repeats, we designed degenerate primer pairs for 19 loci, of which five produced polymorphic fragments in up to 18 mammalian species, including the distinctly related marsupials and monotremes, groups that diverged from other mammals 120-160 million years ago. Using our method, many more cross-clade microsatellite loci can be harvested from the currently available genomic data, and this ability is set to improve exponentially as further genomes are sequenced.