Motif mismatches in microsatellites: insights from genome-wide investigation among 20 insect species

We present a detailed genome-wide comparative study of motif mismatches of microsatellites among 20 insect species representing five taxonomic orders. The results show that varying proportions (∼15–46%) of microsatellites identified in these species are imperfect in motif structure, and that they also vary in chromosomal distribution within genomes. It was observed that the genomic abundance of imperfect repeats is significantly associated with the length and number of motif mismatches of microsatellites. Furthermore, microsatellites with a higher number of mismatches tend to have lower abundance in the genome, suggesting that sequence heterogeneity of repeat motifs is a key determinant of genomic abundance of microsatellites. This relationship seems to be a general feature of microsatellites even in unrelated species such as yeast, roundworm, mouse and human. We provide a mechanistic explanation of the evolutionary link between motif heterogeneity and genomic abundance of microsatellites by examining the patterns of motif mismatches and allele sequences of single-nucleotide polymorphisms identified within microsatellite loci. Using Drosophila Reference Genetic Panel data, we further show that pattern of allelic variation modulates motif heterogeneity of microsatellites, and provide estimates of allele age of specific imperfect microsatellites found within protein-coding genes.     

Triplet-repeat microsatellites shared among hard and soft pines

Vascular plant species have shown a low level of microsatellite conservation compared to many animal species. Finding trans-specific microsatellites for plants may be improved by using a priori knowledge of genome organization. Fifteen triplet-repeat microsatellites from hard pine (Pinus taeda L.) were tested for trans-specific amplification across seven hard pines (P. palustris Mill., P. echinata Mill., P. radiata D. Don., P. patula Schiede et Deppe, P. halepensis Mill., P. kesiya Royle), a soft pine (P. strobus L.), and Picea rubens Sargent. Seven of 15 microsatellites had trans-specific amplification in both hard and soft pine subgenera. Two P. taeda microsatellites had conserved flanking regions and repeat motifs in all seven hard pines, soft pine P. strobus, and P. rubens. Perfect triplet-repeat P. taeda microsatellites appear to be better candidates for trans-specific polymorphism than compound microsatellites. Not all perfect triplet-repeat microsatellites were conserved, but all conserved microsatellites had perfect repeat motifs. Persistent microsatellites PtTX2123 and PtTX3020 had highly conserved flanking regions and a conserved repeat motif composition with variable repeat unit numbers. Using trinucleotide microsatellites improved trans-specific microsatellite recovery among hard and soft pine species.     

In silico mining of microsatellites and analysis of genetic diversity among inter‐ and intra‐generic aphids of the subfamily Aphidinae

Nearly 5 000 aphid species damage crops, either by sucking plant sap or as disease‐transmitting vectors. Microsatellites are used for understanding molecular diversity and eco‐geographical relationships among aphid species. Expressed sequence tag (EST)‐microsatellite motifs were identified through an in silico approach using inbuilt simple sequence repeat mining tools in aphid EST dataset. Microsatellite mining revealed one in every five aphid genes as containing a repeat motif, and out of 9 290 EST microsatellites mined from Aphis gossypii Glover and Acyrthosiphon pisum (Harris) (both Hemiptera: Aphididae), 80% were of A and/or T (AT, ATA, AAT, AATA, and ATTT) motifs, and the rest contained G and/or C motifs. All microsatellite sequences were annotated using BLAST. Primers for EST microsatellites were designed using the Primer 3.0 tool. 106 primer pairs of both dinucleotide repeats (DNRs) and trinucleotide repeats (TNRs), representing open reading frames (ORFs) and untranslated regions (UTRs), were synthesized to amplify 15 aphid species belonging to the subfamily Aphidinae, collected from diverse hosts. Four hundred forty‐five polymorphic alleles were amplified. Fifty TNR and 23 DNR microsatellites amplified across the species studied. Polymorphism information content values of microsatellites ranged from 0.23 to 0.91, amplifying 2–16 alleles. Genetic similarity indices were estimated using the ‘NTSYS‐pc’ software package. Unweighted pair group with arithmetic mean and principal component analysis resolved taxonomic relationships of the aphid species studied. The new aphid microsatellites developed will provide valuable information to researchers to study Indian aphid species diversity and genetic relationships.     

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.     

Characterization and isolation of novel microsatellites from the Drosophila dunni subgroup

We have isolated and characterized 77 novel microsatellites from two species, Drosophila dunni and Drosophila nigrodunni, which are closely related Caribbean-island endemics from the Drosophila cardini species group. These species are very distantly related to all other Drosophila from which microsatellites have previously been characterized. We find that the average length of microsatellites isolated in these species is quite small, with an overall mean length of 9.8 repeat units for dinucleotide microsatellites in the two study species. The nucleotide composition of dinucleotides differs between the two species: D. nigrodunni has a predominance of (AC/GT)n repeats, whereas D. dunni has equal numbers of (AC/GT)n and (AG/CT)n repeats. Tri- and tetranucleotide repeats are not abundant in either species. We assayed the variability of eight microsatellites in a closely related third species, Drosophila arawakana, using wild-caught individuals from the island of Guadeloupe. We found the microsatellites to be extremely variable in this population, with observed heterozygosities ranging from 0.541 to 0.889. DNA amplification trials suggest that these eight microsatellites are widely conserved across the D. cardini group, with five of the eight producing amplification products in every species tested. However, the loci are very poorly conserved over greater phylogenetic distances. DNA amplification of the microsatellite loci was unreliable in members of the closely related Drosophila quinaria, Drosophila calloptera, Drosophila guarani and Drosophila tripunctata species groups. Furthermore, these microsatellites could not be detected in the genome of Drosophila melanogaster, despite the conservation of microsatellite flanking regions at some loci. These data indicate that Drosophila microsatellite loci are quite short lived over evolutionary timescales relative to many other taxa.     

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.     

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.     

Genetic variation and differentiation at microsatellite loci in Drosophila simulans. Evidence for founder effects in new world populations.

Drosophila simulans isofemale lines from Africa, South America, and two locations in North America were surveyed for variation at 16 microsatellite loci on the X, second, and third chromosomes, and 18 microsatellites, which are unmapped. D. simulans is thought to have colonized New World habitats only relatively recently (within the last few hundred years). Consistent with a founder effect occurring as colonizers moved into these New World habitats, we find less microsatellite variability in North and South American D. simulans populations than for an African population. Population subdivision as measured at microsatellites is moderate when averaged across all loci (FST = 0.136), but contrasts sharply with previous studies of allozyme variation, which have showed significantly less differentiation in D. simulans than in D. melanogaster. There are substantially fewer private alleles observed in New World populations of D. simulans than seen in a similar survey of D. melanogaster. In addition to possible differences in population size during their evolutionary histories, varying colonization histories or other demographic events may be necessary to explain discrepancies in the patterns of variation observed at various genetic markers between these closely related species.     

Microsatellite allele sequencing in population analyses of the South American cactophilic species Drosophila antonietae (Diptera: Drosophilidae)

Drosophila antonietae belongs to the Drosophila buzzatii cluster, a cactophilic group of species naturally endemic to South America. Morphological and genetic analyses indicate that its populations are the most homogenous in the cluster and that the diversity observed is mainly a result of variation within populations. Seven polymorphic microsatellite loci were described for this species and used in the present study to investigate the genetic diversity of natural populations of D. antonietae by both length and sequence variation. The study aimed to understand how homoplasy and null alleles affect inferences about the population history of this species and to obtain an accurate interpretation of population inferences where these loci could be applied. The results provide useful information on the interpretation of genetic data derived from the microsatellite loci described for D. antonietae and on evolutionary aspects of cactophilic Drosophila. Importantly, the results indicate that size homoplasy and null alleles do not represent significant problems for the population genetics analyses because the large amount of variability at microsatellite loci compensate the low frequency of these problems in the populations. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 573–584.     

Using Next Generation RAD Sequencing to Isolate Multispecies Microsatellites for Pilosocereus (Cactaceae)

Microsatellite markers (also known as SSRs, Simple Sequence Repeats) are widely used in plant science and are among the most informative molecular markers for population genetic investigations, but the development of such markers presents substantial challenges. In this report, we discuss how next generation sequencing can replace the cloning, Sanger sequencing, identification of polymorphic loci, and testing cross-amplification that were previously required to develop microsatellites. We report the development of a large set of microsatellite markers for five species of the Neotropical cactus genus Pilosocereus using a restriction-site-associated DNA sequencing (RAD-seq) on a Roche 454 platform. We identified an average of 165 microsatellites per individual, with the absolute numbers across individuals proportional to the sequence reads obtained per individual. Frequency distribution of the repeat units was similar in the five species, with shorter motifs such as di- and trinucleotide being the most abundant repeats. In addition, we provide 72 microsatellites that could be potentially amplified in the sampled species and 22 polymorphic microsatellites validated in two populations of the species Pilosocereus machrisii. Although low coverage sequencing among individuals was observed for most of the loci, which we suggest to be more related to the nature of the microsatellite markers and the possible bias inserted by the restriction enzymes than to the genome size, our work demonstrates that an NGS approach is an efficient method to isolate multispecies microsatellites even in non-model organisms.     

Bioinformatic survey of Haliotis midae microsatellites reveals a non-random distribution of repeat motifs

Recent studies have shown the non-random distribution of microsatellite motifs between genomic regions within a particular species. This study investigates such microsatellite distributions in the genome of the economically important abalone Haliotis midae, via a bioinformatic survey. In particular, the association of specific repeat motifs to coding regions and transposable elements is investigated. An understanding of microsatellite genomic distribution will facilitate more efficient use and development of this popular molecular marker. A bias toward di- and tetranucleotide repeats was found in the H. midae genome. CA microsatellite units were the most abundant repeat motif, but were notably underrepresented in genic regions where GAGT repeats predominate. Approximately 17.5% and 21% of the microsatellites showed gene and/or transposable element associations, respectively. This could explain the high genomic frequencies of particular motifs across the genome and may allude to a possible functional role. The data presented in this study are the first to demonstrate such non-random dispersal of microsatellites in abalone and support previous findings arguing in favor of non-random distribution of repeat motifs.     

Development, inheritance and cross-species amplification of microsatellite markers from Acacia mangium

Microsatellite markers were developed in Acacia mangium Willd. to provide highly variable co-dominant markers for linkage mapping and studies of the breeding system. After an enrichment procedure 40% of colonies contained microsatellites in contrast with less than 1% from a non-enriched library. The majority of microsatellite sequences were AC repeats. Co-dominant segregation of alleles in two full-sib crosses of A. mangium was demonstrated at 33 microsatellite loci. The markers were highly variable relative to restriction fragment lengths polymorphisms (RFLPs). In the two pedigrees 53% of microsatellite loci were fully informative compared with 15% of RFLPs. Based on alleles detected among four parental genotypes, the microsatellites consisting of dinucleotide repeats were more polymorphic than those with tri- and tetra-nucleotide repeats. The microsatellite markers were not as transferable across species in the genus Acacia as RFLPs. Two thirds of the primers developed in A. mangium (subgenus Phyllodineae, section Juliflorae) amplified DNA from other species within the same section but failed to amplify in species from the subgenus Acacia. The availability of multiallelic, PCR-based, co-dominant microsatellite loci makes possible efficient studies of gene flow and breeding systems in A. mangium, a species with low allozyme variation. Received: 30 December 1999 / Accepted: 10 May 2000     

Conservation of locus-specific microsatellite variability across species: a comparison of two Drosophila sibling species, D. melanogaster and D. simulans

Fifteen microsatellite loci were studied in Drosophila melanogaster and Drosophila simulans, two closely related sibling species which split 2- 3.5 MYA. Within-species variances in repeat number were found to differ up to 1,000-fold among individual microsatellite loci. A significant correlation of log variances between both species indicated a locus- specific mutation rate of microsatellites. Hence, locus-specific effects are apparently among the major forces influencing microsatellite variation and deserve more consideration in microsatellite analysis.      

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.     

Phylogenetic Patterns of Geographical and Ecological Diversification in the Subgenus Drosophila

Colonisation of new geographic regions and/or of new ecological resources can result in rapid species diversification into the new ecological niches available. Members of the subgenus Drosophila are distributed across the globe and show a large diversity of ecological niches. Furthermore, taxonomic classification of Drosophila includes the rank radiation, which refers to closely related species groups. Nevertheless, it has never been tested if these taxonomic radiations correspond to evolutionary radiations. Here we present a study of the patterns of diversification of Drosophila to test for increased diversification rates in relation to the geographic and ecological diversification processes. For this, we have estimated and dated a phylogeny of 218 species belonging to the major species groups of the subgenus. The obtained phylogenies are largely consistent with previous studies and indicate that the major groups appeared during the Oligocene/Miocene transition or early Miocene, characterized by a trend of climate warming with brief periods of glaciation. Ancestral reconstruction of geographic ranges and ecological resource use suggest at least two dispersals to the Neotropics from the ancestral Asiatic tropical disribution, and several transitions to specialized ecological resource use (mycophagous and cactophilic). Colonisation of new geographic regions and/or of new ecological resources can result in rapid species diversification into the new ecological niches available. However, diversification analyses show no significant support for adaptive radiations as a result of geographic dispersal or ecological resource shift. Also, cactophily has not resulted in an increase in the diversification rate of the repleta and related groups. It is thus concluded that the taxonomic radiations do not correspond to adaptive radiations.     

Genetic differentiation and demographic history in Drosophila pachea from the Sonoran Desert

Genetic variation at six microsatellite DNA loci and a segment of the mitochondrial cytochrome oxidase subunit I (COI) locus was used to estimate gene flow, population structure, and demographic history in the cactophilic Drosophila pachea from the Sonoran Desert of North America, a species that shows a strict association with its senita host cactus (genus Lophocereus). For microsatellite analyses, thirteen populations of D. pachea were sampled, five in mainland Mexico and the southwestern USA, and eight on the Baja California (Baja) peninsula, covering essentially the entire range of the species. Analysis of molecular variance (AMOVA) of microsatellite data revealed that populations from both the mainland and the Baja peninsula generally showed little structure, although there were a few exceptions, suggesting some local differentiation and restriction of gene flow within both regions. Pairwise comparisons of F(ST) among each of the mainland and Baja populations showed evidence of both panmixia and population subdivision. AMOVA performed on grouped populations from both the mainland and Baja, however, revealed significant partitioning of genetic variation among the two regions, but no partitioning among localities within each region. Bayesian skyline analyses of the COI data set, consisting of four mainland and seven peninsular populations, revealed population expansions dating to the Pleistocene or late Pliocene in D. pachea from both regions, although regional differences were seen in the estimated timing of the expansions and in changes in effective population size over time.     

Rapid expansion of microsatellite sequences in pines

Microsatellite persistence time and evolutionary change was studied among five species of pines, which included a pair of closely related species (Pinus sylvestris and Pinus resinosa) in the subgenus Pinus, their relative Pinus radiata, and another closely related species pair (Pinus strobus and Pinus lambertiana) in the subgenus Strobus. The effective population sizes of these species are known to have ranged from the very small bottlenecks of P. resinosa to vast populations of P. sylvestris. This background allowed us to place the microsatellite evolution in a well-defined phylogenetic setting. Of 30 loci originating from P. strobus and P. radiata, we were able to consistently amplify 4 in most of the these pine species. These priming sites had been conserved for over 100 Myr. The four microsatellites were sequenced in the five species. Flanking sequences were compared to establish that the loci were orthologous. All microsatellites had persisted in these species, despite very different population sizes. We found a recent microsatellite duplication: a closely related pair of loci in P. strobus, where the other four species had just one locus. On two independent occasions, the repeat area of this same microsatellite (locus RPS 105a/b) had grown from a very low repeat number to 15 or 17 in the last 10-25 Myr. Other parts of the same compound microsatellite had remained virtually unchanged. Locus PR 4.6 is known to be polymorphic in both P. radiata and P. sylvestris, but the polymorphism in the two species is due to different motifs. The very large pine genomes are highly repetitive, and microsatellite loci also occur as gene families.     

Evolutionary change in seed size among some legume species: The effects of phylogeny

The importance of phylogenetic effects in controlling seed size variation at the macroevolutionary levels was examined using species of a well-defined, monophyletic family, the Leguminosae. A nested ANOVA was used to separate variance components at the various taxonomic levels such as subfamily, genus, and subgenus. Statistical significance was found at most of the taxonomic levels examined, which suggests that phylogeny as shown by the accepted taxonomy of the family, exerted a substantial influence over seed size variation. Thus, there appears to be a characteristic mean seed size for each genus and subgenus. The overall positive correlation between plant height and seed size was interpreted as a scaling of seed size to plant height. When the effect of plant height was controlled, the amount of variance changed to some extent across taxonomic levels, though resulting in no change in statistical significance at various taxonomic levels. This results indicate that phylogenetic effects in seed size variation among congeneric or consubgeneric groupings were largely independent of allometric effects. The patterns of variance components and intraclass correlation coefficients for seed size differed to some extent among subfamilies. The difference in those patterns among taxa will reflect evolutionary change more closely with larger sample sizes. A simple inverse relationship between seed size and number per fruit among related legume species was difficult to generalize. This was perhaps due to phylogenetic constraints on seed size and/or seed number (or ovule number). These results demonstrate that phylogeny and plant height affect seed size variation among related species.     

Evolutionary history of microsatellites in the obscura group of Drosophila

The evolutionary origins of microsatellites are not well understood. Some investigators have suggested that point mutations that expand repeat arrays beyond a threshold size trigger microsatellites to become variable. However, little empirical data has been brought forth on this and related issues. In this study, we examine the evolutionary history of microsatellites in six species within the obscura group of Drosophila, tracing changes in microsatellite alleles using both PCR product size and sequence data. We found little evidence supporting a general role of point mutations triggering initial microsatellite expansion, and no consistent threshold size for expansion was observed. Flanking region length variation was extensive when alleles were sequenced in distantly related species, and some species possessed altogether different repeat arrays between the same primer binding sites. Our results suggest extreme caution in using microsatellite allele sizes for phylogenetic analyses or to infer divergences between populations.     

What phylogeny and gene genealogy analyses reveal about homoplasy in citrus microsatellite alleles

Sixty-five microsatellite alleles amplified from ancestral citrus accessions classified in three separate genera were evaluated for sequence polymorphism to establish the basis of inter- and intra-allelic genetic variation, evaluate the extent of size homoplasy, and determine an appropriate model (stepwise or infinite allele) for analysis of citrus microsatellite alleles. Sequences for each locus were aligned and subsequently used to determine relationships between alleles of different taxa via parsimony. Interallelic size variation at each SSR locus examined was due to changes in repeat copy number with one exception. Sequencing these alleles uncovered new distinct point mutations in the microsatellite region and the region flanking the microsatellite. Several of the point mutations were found to be genus, species, or allele specific, and some mutations were informative about the inferred evolutionary relationships among alleles. Overall, homoplasy was observed in alleles from all three loci, where the core microsatellite repeat was changed causing alleles of the same size class to be identical in state but not identical by descent. Because nearly all changes in allele size (with one exception) were due to expansion or contraction of the repeat motif, this suggests that a stepwise mutation model, which assumes homoplasy may occur, would be the most appropriate for analyzing Citrus SSR data. The collected data indicate that microsatellites can be a useful tool for evaluating Citrus species and two related genera since repeat motifs were reasonably well retained. However, this work also demonstrated that the number of microsatellite alleles is clearly an underestimate of the number of sequence variants present.