Combination of multiple microsatellite data sets to investigate genetic diversity and admixture of domestic cattle

Microsatellite markers are commonly used for population genetic analyses of livestock. However, up to now, combinations of microsatellite data sets or comparison of population genetic parameters from different studies and breeds has proven difficult. Often different genotyping methods have been employed, preventing standardization of microsatellite allele calling. In other cases different sets of markers have been genotyped, providing differing estimates of population genetic parameters. Here, we address these issues and illustrate a general two-step regression approach in cattle using three different sets of microsatellite data, to combine population genetics estimates of diversity and admixture. This regression-based method is independent of the loci genotyped but requires common breeds in the data sets. We show that combining microsatellite data sets can provide new insights on the origin and geographical distribution of genetic diversity and admixture in cattle, which will facilitate global management of this livestock species.     

Highly polymorphic microsatellite markers in poultry

Microsatellite markers have been established for a large number of species, but up till now very few polymorphic microsatellite markers have been reported in poultry. We have isolated 34 polymorphic chicken microsatellite markers of the poly (TG) type. The number of repeats varied from 9 up to 33. Often, other repeats such as poly (T) or poly (GAA) were present adjacent to the poly (TG) repeat. Polymerase chain reaction amplification of the microsatellites resulted in detection of three or more alleles in a test panel of five different animals for 75% of the microsatellites. Segregation of five microsatellite markers has been tested in a small family.     

Analysis of microsatellite markers in the genome of the plant pathogen Ceratocystis fimbriata

Ceratocystis fimbriata sensu lato represents a complex of cryptic and commonly plant pathogenic species that are morphologically similar. Species in this complex have been described using morphological characteristics, intersterility tests and phylogenetics. Microsatellite markers have been useful to study the population structure and origin of some species in the complex. In this study we sequenced the genome of C. fimbriata. This provided an opportunity to mine the genome for microsatellites, to develop new microsatellite markers, and map previously developed markers onto the genome. Over 6000 microsatellites were identified in the genome and their abundance and distribution was determined. Ceratocystis fimbriata has a medium level of microsatellite density and slightly smaller genome when compared with other fungi for which similar microsatellite analyses have been performed. This is the first report of a microsatellite analysis conducted on a genome sequence of a fungal species in the order Microascales. Forty-seven microsatellite markers have been published for population genetic studies, of which 35 could be mapped onto the C. fimbriata genome sequence. We developed an additional ten microsatellite markers within putative genes to differentiate between species in the C. fimbriata s.l. complex. These markers were used to distinguish between 12 species in the complex.     

花生微卫星标记的研究进展

近年来花生微卫星标记的开发取得了一定的进展, 初步揭示了花生在DNA水平上的遗传多样性。花生微卫星标记的开发途径主要包括通过构建小片段基因组文库开发基因组SSR标记, 根据花生EST序列开发EST-SSR标记, 根据豆科植物序 列信息和SSR标记开发花生SSR标记, 将SSR标记与其它分子标记结合开发新的DNA标记, 以及基于SSR核心序列开发ISSR标记。花生微卫星标记主要应用于遗传多样性研究、遗传图谱与品种指纹图谱构建以及分子标记辅助育种等领域。本文综述了花生SSR标记开发研究的进展及应用。     

Inference on microsatellite mutation processes in the invasive mite, Varroa destructor, using reversible jump Markov chain Monte Carlo

Varroa destructor is a parasitic mite of the Eastern honeybee Apis cerana. Fifty years ago, two distinct evolutionary lineages (Korean and Japanese) invaded the Western honeybee Apis mellifera. This haplo-diploid parasite species reproduces mainly through brother-sister matings, a system which largely favors the fixation of new mutations. In a worldwide sample of 225 individuals from 21 locations collected on Western honeybees and analyzed at 19 microsatellite loci, a series of de novo mutations was observed. Using historical data concerning the invasion, this original biological system has been exploited to compare three mutation models with allele size constraints for microsatellite markers: stepwise (SMM) and generalized (GSM) mutation models, and a model with mutation rate increasing exponentially with microsatellite length (ESM). Posterior probabilities of the three models have been estimated for each locus individually using reversible jump Markov Chain Monte Carlo. The relative support of each model varies widely among loci, but the GSM is the only model that always receives at least 9% support, whatever the locus. The analysis also provides robust estimates of mutation parameters for each locus and of the divergence time of the two invasive lineages (67,000 generations with a 90% credibility interval of 35,000-174,000). With an average of 10 generations per year, this divergence time fits with the last post-glacial Korea-Japan land separation.     

Population divergence with or without admixture: selecting models using an ABC approach

Genetic data have been widely used to reconstruct the demographic history of populations, including the estimation of migration rates, divergence times and relative admixture contribution from different populations. Recently, increasing interest has been given to the ability of genetic data to distinguish alternative models. One of the issues that has plagued this kind of inference is that ancestral shared polymorphism is often difficult to separate from admixture or gene flow. Here, we applied an approximate Bayesian computation (ABC) approach to select the model that best fits microsatellite data among alternative splitting and admixture models. We performed a simulation study and showed that with reasonably large data sets (20 loci) it is possible to identify with a high level of accuracy the model that generated the data. This suggests that it is possible to distinguish genetic patterns due to past admixture events from those due to shared polymorphism (population split without admixture). We then apply this approach to microsatellite data from an endangered and endemic Iberian freshwater fish species, in which a clustering analysis suggested that one of the populations could be admixed. In contrast, our results suggest that the observed genetic patterns are better explained by a population split model without admixture.     

花生微卫星标记的研究进展

近年来花生微卫星标记的开发取得了一定的进展,初步揭示了花生在DNA水平上的遗传多样性。花生微卫星标记的开发途径主要包括通过构建小片段基因组文库开发基因组SSR标记,根据花生EST序列开发EST-SSR标记,根据豆科植物序列信息和SSR标记开发花生SSR标记,将SSR标记与其它分子标记结合开发新的DNA标记,以及基于SSR核心序列开发ISSR标记。花生微卫星标记主要应用于遗传多样性研究、遗传图谱与品种指纹图谱构建以及分子标记辅助育种等领域。本文综述了花生SSR标记开发研究的进展及应用。     

Polarity of mutations in tumor-associated microsatellite instability

Many factors have been implicated in influencing the rate of microsatellite mutations, including the length and base composition of the repeat motif, number of repeats, base composition of flanking sequences and, perhaps most importantly, degree of perfection of the repeats. The latter is of clinical relevance, since in both spino-cerebellar ataxia and fragile X syndrome, alleles with imperfect repeats appear to be much more stable than perfect ones. As yet, the relative importance of increased replication slippage and decreased mismatch repair efficiency in the preference of mutations to occur within perfect repeats has not been fully determined. D13S308E is an asymmetric trinucleotide repeat microsatellite with the sequence (CAT)₃CAC(CAT)CAC(CAT)₂CAC(CAT)CAC(CAT) ₁₅ , thus containing two parts: an 11-repeat imperfect portion (underlined above) and a 15-repeat perfect one (bold). We sequenced eight new mutant alleles of D13S308E from three human gastric tumors with instability in this and other microsatellites. In all mutations the size variation occurred exclusively in the perfect part of the microsatellite. These results constitute direct evidence that the molecular basis of microsatellite alterations seen in normal cells is similar to those that occur in human tumors with extensive microsatellite instability. The investigation of mechanisms involved in microsatellite mutations has been handicapped by the fact that they are rare events. The microsatellite instability observed in malignant tumors provides us with a useful general system to study these mechanisms. Received: 24 June 1997 / Accepted 7 October 1997     

微卫星位点筛选方法综述

微卫星标记因其丰富的多态性和共显性等特点,已得到了广泛的应用.应用微卫星标记首先需要获得微卫星位点的序列信息,用来设计引物.获得微卫星位点的方法有多种,本文综述了获得和富集微卫星位点的常用方法.最简便、最省时的方法是从公共数据库(如EMBL、Genbank、EST数据库等)或已发表的文献中查找到微卫星位点,但只限于已经有序列数据发布的物种.第二种方法是种间转移扩增,即从相近物种的数据库中查找微卫星位点,或使用已有数据发表的遗传距离相近物种的微卫星标记.第三种方法是从基因组DNA中筛选微卫星位点,其中用于富集微卫星的方法有引物法、磁珠杂交法、尼龙膜杂交法以及RAPD技术法.     

Challenges of microsatellite isolation in fungi

Although they represent powerful genetic markers in many fields of biology, microsatellites have been isolated in few fungal species. The aim of this study was to assess whether obtaining microsatellite markers with an acceptable level of polymorphism is generally harder from fungi than in other organisms. We therefore surveyed the number, nature and polymorphism level of published microsatellite markers in fungi from the literature and from our own data on seventeen fungal microsatellite-enriched libraries, and in five other phylogroups (angiosperms, insects, fishes, birds and mammals). Fungal microsatellites indeed appeared both harder to isolate and to exhibit lower polymorphism than in other organisms. This appeared to be due, at least in part, to genomic specificities, such as scarcity and shortness of fungal microsatellite loci. A correlation was observed between mean repeat number and mean allele number in the published fungal microsatellite loci. The cross-species transferability of fungal microsatellites also appeared lower than in other phylogroups. However, microsatellites have been useful in some fungal species. Thus, the considerable advantages of these markers make their development worthwhile, and this study provides some guidelines for their isolation.     

Polymorphic microsatellite loci for the genetic analysis of Lycoris radiata (Amaryllidaceae) and cross-amplification in other congeneric species

Lycoris radiata is a perennial herb that has been used in traditional Chinese medicine for a long time and has two main medicinal components in its bulb, lycorine and galanthamine. However, the original microsatellite loci have not been developed for any species of Lycoris. Total genomic DNA was extracted from fresh bulbs using a modified CTAB protocol. We isolated 10 microsatellite loci from 21 L. radiata individuals of a natural population from Yellow Mountain in Anhui Province, China. The number of alleles ranged from two to nine. The observed and expected heterozygosities ranged from 0.238 to 0.952 and from 0.455 to 0.784, respectively. One locus significantly deviated from Hardy-Weinberg equilibrium and no significant linkage disequilibrium was found between pairs of loci. Cross-species amplification of these microsatellite loci was characterized in additional five species (L. sprengeri, L. anhuiensis, L. albiflora, L. longituba, and L. chinensis) of Lycoris. The results suggest that these microsatellite markers would contribute to the population genetic studies of L. radiata and other related species.     

Successful development of microsatellite markers in a challenging species: the horizontal borer Austroplatypus incompertus (Coleoptera: Curculionidae)

The analysis of microsatellite loci has allowed significant advances in evolutionary biology and pest management. However, until very recently, the potential benefits have been compromised by the high costs of developing these neutral markers. High-throughput sequencing provides a solution to this problem. We describe the development of 13 microsatellite markers for the eusocial ambrosia beetle, Austroplatypus incompertus, a significant pest of forests in southeast Australia. The frequency of microsatellite repeats in the genome of A. incompertus was determined to be low, and previous attempts at microsatellite isolation using a traditional genomic library were problematic. Here, we utilised two protocols, microsatellite-enriched genomic library construction and high-throughput 454 sequencing and characterised 13 loci which were polymorphic among 32 individuals. Numbers of alleles per locus ranged from 2 to 17, and observed and expected heterozygosities from 0.344 to 0.767 and from 0.507 to 0.860, respectively. These microsatellites have the resolution required to analyse fine-scale colony and population genetic structure. Our work demonstrates the utility of next-generation 454 sequencing as a method for rapid and cost-effective acquisition of microsatellites where other techniques have failed, or for taxa where marker development has historically been both complicated and expensive.     

Androgen receptor CAG repeat polymorphism in males of six non-human primate species

Background Androgen receptor [CAG]_n microsatellite has been linked to human diseases. Methods Six non‐human primates were genotyped for the [CAG]_n microsatellite. Results Marmosets and macaques are monomorphic, while mangabeys, baboons, and chimpanzees are polymorphic. Conclusions Non‐human primates that are polymorphic for the microsatellite are candidate animal models for CAG‐related diseases.     

Abundance, polymorphism and genetic mapping of microsatellites in rice

Dinucleotide microsatellites have been characterized and used as genetic markers in rice. Screening of a rice genomic library with poly(dG-dA)·(dC-dT) and poly(dG-dT)·(dC-dA) probes indicated that (GA)_n repeats occurred, on average, once every 225 kb and (GT)_n repeats once every 480 kb. DNA sequencing of ten randomly selected microsatellites indicated that the numbers of repeats ranged from 12 to 34 and that the patterns of microsatellites in rice were similar to those of humans and other mammals. Primers to these microsatellite loci as well as to four published microsatellite-containing sequences have been designed and degrees of polymorphism has been examined with 20 rice accessions. Multiple alleles, ranging from 5 to 11, have been observed at all the microsatellite loci in 20 rice accessions. Alleles specific to two cultivated subspecies, indica and japonica, were found in some microsatellite loci. Heterozygosity values of all the microsatellite markers were significantly higher than those of RFLP markers, based upon a parallel comparison. Ten microsatellite loci have been genetically mapped to four rice chromosomes. The genomic distribution of microsatellites appears to be random in rice.     

Isolation and characterization of microsatellite loci in Peganum harmala (Peganaceae), an important resist-drought and medicinal plant

Peganum harmala is a herb grows spontaneously in arid and rocky areas. From ancient time, it has been claimed to be an important medicinal plant in rough environment. In this study, we developed 12 microsatellite loci from P. harmala by the combining biotin capture method for the first time. A total of 31 microsatellite sequences were recovered through screening the library and 12 of them are polymorphic. The number of alleles per locus in 36 sampled individuals ranged from 3 to 8, expected heterozygosity and observed heterozygosity ranged from 0.1381 to 0.6821 and from 0.3573 to 0.8739, respectively. In addition, all markers have been crossly checked in the other congeneric species. These microsatellite markers would provide a useful tool for investigating the genetic diversity and study the population genetic structure in detail.     

Microsatellite cross-amplification in coccolithophores: application in population diversity studies

The development and isolation of microsatellites entails a significant input of time and money. Therefore there is an interest in using existing microsatellites on species from which markers have not yet been developed. Conservation of six previously identified microsatellite loci in the marine coccolithophorid species Emiliana huxleyi was found in a survey of two bloom forming coccolithophorid species--Gephyrocapsa oceanica and Coccolithus pelagicus. The number of alleles per locus varied from 1 to 8, and half of the microsatellite loci tested showed 4 or more alleles. The microsatellite markers used in this study may be applied to other coccolithophorid species for population analysis, eliminating the time-consuming, costly development of microsatellite markers for other coccolithophorid species.     

Microsatellite mutation models: insights from a comparison of humans and chimpanzees

Using genomic data from homologous microsatellite loci of pure AC repeats in humans and chimpanzees, several models of microsatellite evolution are tested and compared using likelihood-ratio tests and the Akaike information criterion. A proportional-rate, linear-biased, one-phase model emerges as the best model. A focal length toward which the mutational and/or substitutional process is linearly biased is a crucial feature of microsatellite evolution. We find that two-phase models do not lead to a significantly better fit than their one-phase counterparts. The performance of models based on the fit of their stationary distributions to the empirical distribution of microsatellite lengths in the human genome is consistent with that based on the human-chimp comparison. Microsatellites interrupted by even a single point mutation exhibit a twofold decrease in their mutation rate when compared to pure AC repeats. In general, models that allow chimps to have a larger per-repeat unit slippage rate and/or a shorter focal length compared to humans give a better fit to the human-chimp data as well as the human genomic data.     

A possible explanation for the population size discrepancy in tuna (genus Thunnus) estimated from mitochondrial DNA and microsatellite data

A recent study using both mitochondrial DNA (mtDNA) and microsatellite data reported on a population size discrepancy in the eastern tiger salamander where the effective population size (N_e) estimate of the former exceeded that of the latter. That study suggested, among other hypotheses, that homoplasy of microsatellite alleles is responsible for the discrepancy. In this investigation, we report 10 new cases of a similar discrepancy in five species of tuna. These cases derive from our Bayesian inferences using data from Pacific Bluefin Tuna (Thunnus orientalis) and Yellowfin Tuna (Thunnus albacares), as well as from published estimates of genetic diversity for additional populations of Yellowfin Tuna and three other tuna species. Phylogenetic character analyses of inferred genealogies of Pacific Bluefin and Yellowfin Tuna reveal similar reduced levels of mtDNA and microsatellite homoplasy. Thus, the discrepancy between inferred population sizes from mtDNA and microsatellite data in tuna is most likely not an artifact of the chosen mutation models used in the microsatellite analyses, but may reflect behavioral differences between the sexes such as female-biased philopatry and male-biased dispersal. This explanation now warrants critical testing with more local populations of tuna and with other animal and plant groups that have different life histories.     

Complex mutations in a high proportion of microsatellite loci from the protozoan parasite Plasmodium falciparum

Microsatellite loci are generally assumed to evolve via a stepwise mutational process and a battery of statistical techniques has been developed in recent years based on this or related mutation models. It is therefore important to investigate the appropriateness of these models in a wide variety of taxa. We used two approaches to examine mutation patterns in the malaria parasite Plasmodium falciparum: (i) we examined sequence variation at 12 tri-nucleotide repeat loci; and (ii) we analysed patterns of repeat structure and heterozygosity at 114 loci using data from 12 laboratory parasite lines. The sequencing study revealed complex patterns of mutation in five of the 12 loci studied. Alleles at two loci contain indels of 24 bp and 57 bp in flanking regions, while in the other three loci, blocks of imperfect microsatellites appear to be duplicated or inserted; these loci essentially consist of minisatellite repeats, with each repeat unit containing four to eight microsatellites. The survey of heterozygosity revealed a positive relationship between repeat number and microsatellite variability for both di- and trinucleotides, indicating a higher mutation rate in loci with longer repeat arrays. Comparisons of levels of variation in different repeat types indicate that the mutation rate of dinucleotide-bearing loci is 1.6-2.1 times faster than trinucleotides, consistent with the lower mean number of repeats in trinucleotide-bearing loci. However, despite the evidence that microsatellite arrays themselves are evolving in a manner consistent with stepwise mutation model in P. falciparum, the high frequency of complex mutations precludes the use of analytical tools based on this mutation model for many microsatellite-bearing loci in this protozoan. The results call into question the generality of models based on stepwise mutation for analysing microsatellite data, but also demonstrate the ease with which loci that violate model assumptions can be detected using minimal sequencing effort.     

Likelihood-based estimation of microsatellite mutation rates

Microsatellites are widely used in genetic analyses, many of which require reliable estimates of microsatellite mutation rates, yet the factors determining mutation rates are uncertain. The most straightforward and conclusive method by which to study mutation is direct observation of allele transmissions in parent-child pairs, and studies of this type suggest a positive, possibly exponential, relationship between mutation rate and allele size, together with a bias toward length increase. Except for microsatellites on the Y chromosome, however, previous analyses have not made full use of available data and may have introduced bias: mutations have been identified only where child genotypes could not be generated by transmission from parents' genotypes, so that the probability that a mutation is detected depends on the distribution of allele lengths and varies with allele length. We introduce a likelihood-based approach that has two key advantages over existing methods. First, we can make formal comparisons between competing models of microsatellite evolution; second, we obtain asymptotically unbiased and efficient parameter estimates. Application to data composed of 118,866 parent-offspring transmissions of AC microsatellites supports the hypothesis that mutation rate increases exponentially with microsatellite length, with a suggestion that contractions become more likely than expansions as length increases. This would lead to a stationary distribution for allele length maintained by mutational balance. There is no evidence that contractions and expansions differ in their step size distributions.