Comparison of random and SSR-enriched shotgun pyrosequencing for microsatellite discovery and single multiplex PCR optimization in Acacia harpophylla F. Muell. Ex Benth

Streamlining the development and genotyping of microsatellites in species for which no genetic information is available represents an important technical challenge to overcome in order to enable mainstream application of state-of-the-art population genetic analysis techniques in nonmodel organisms. Using the example of Acacia harpophylla, an acacia tree endemic of north-eastern Australia, we show that high-throughput shotgun pyrosequencing technology, so-called second-generation sequencing, reduces time and cost of microsatellite marker discovery in nonmodel organisms and of their large-scale typing in natural populations. We found that 0.5% of short sequence reads generated on 454 Genome Sequencer FLX Titanium from random genome sampling and 2.2% of reads generated with prior microsatellite enrichment yielded microsatellite markers with designed polymerase chain reaction (PCR) primers, suggesting that enrichment increases efficiency of pyrosequencing when microsatellite discovery is the primary goal. Using stringent selection criteria to facilitate downstream PCR multiplex design, we identified 1435 microsatellite loci with designed primers from a total of 200,908 short sequence reads. From a subset of 96 loci tested for amplification, 38 were validated for population genetics applications, leading to the optimization of a cost-effective multiplex PCR protocol for the simultaneous typing of nine microsatellites in natural populations of A. harpophylla.     

Development of microsatellite genetic markers in Siberian stone pine (<Emphasis Type="Italic">Pinus sibirica</Emphasis> Du Tour) based on the <Emphasis Type="Italic">de novo</Emphasis> whole genome sequencing

Siberian stone pine, Pinus sibirica Du Tour is one of the most economically and environmentally important forest-forming species of conifers in Russia. To study these forests a large number of highly polymorphic molecular genetic markers, such as microsatellite loci, are required. Prior to the new high-throughput next generation sequencing (NGS) methods, discovery of microsatellite loci and development of micro-satellite markers were very time consuming and laborious. The recently developed draft assembly of the Siberian stone pine genome, sequenced using the NGS methods, allowed us to identify a large number of microsatellite loci in the Siberian stone pine genome and to develop species-specific PCR primers for amplification and genotyping of 70 microsatellite loci. The primers were designed using contigs containing short simple sequence tandem repeats from the Siberian stone pine whole genome draft assembly. Based on the testing of primers for 70 microsatellite loci with tri-, tetra- or pentanucleotide repeats, 18 most promising, reliable and polymorphic loci were selected that can be used further as molecular genetic markers in population genetic studies of Siberian stone pine.     

Isolation and characterization of microsatellite loci from the apple maggot fly Rhagoletis pomonella (Diptera: Tephritidae)

We report the isolation and development of 81 novel primers for amplifying microsatellite loci in the Rhagoletis pomonella sibling species complex, and the sequencing, characterization and analysis of basic population genetic parameters for nine of these genes. We also report the successful cross‐species amplification of several of these loci. The R. pomonella sibling species complex is a textbook example of genetic differentiation in sympatry via host‐plant shifting. Microsatellite markers can be useful for mapping host‐plant‐associated adaptations in Rhagoletis that generate reproductive isolation and facilitate speciation, as well as for resolving the genetic structure and evolutionary history of fly populations.     

Isolation and characterization of new microsatellite markers in the surfclam Mactromeris polynyma

We describe primers and amplification conditions for seven microsatellite loci that were developed and characterized in the commercially harvested surfclam Mactromeris polynyma. Five of these loci were polymorphic: we found 6 to 23 alleles per locus among 100 individuals from one population in Nova Scotia, with some large heterozygote deficits that may reflect unrecognized temporal genetic variation. We are using these markers to investigate temporal and spatial genetic structure in this species.     

Conserved Microsatellites in Ants Enable Population Genetic and Colony Pedigree Studies across a Wide Range of Species

Broadly applicable polymorphic genetic markers are essential tools for population genetics, and different types of markers have been developed for this purpose. Microsatellites have been employed as particularly polymorphic markers for over 20 years. However, PCR primers for microsatellite loci are often not useful outside the species for which they were designed. This implies that a new set of loci has to be identified and primers developed for every new study species. To overcome this constraint, we identified 45 conserved microsatellite loci based on the eight currently available ant genomes and designed primers for PCR amplification. Among these loci, we chose 24 for in-depth study in six species covering six different ant subfamilies. On average, 11.16 of these 24 loci were polymorphic and in Hardy-Weinberg equilibrium in any given species. The average number of alleles for these polymorphic loci within single populations of the different species was 4.59. This set of genetic markers will thus be useful for population genetic and colony pedigree studies across a wide range of ant species, supplementing the markers available for previously studied species and greatly facilitating the study of the many ant species lacking genetic markers. Our study shows that it is possible to develop microsatellite loci that are both conserved over a broad range of taxa, yet polymorphic within species. This should encourage researchers to develop similar tools for other large taxonomic groups.     

De novo discovery and multiplexed amplification of microsatellite markers for black alder (Alnus glutinosa) and related species using SSR-enriched shotgun pyrosequencing

Recent developments in sequencing technologies and bioinformatics analyses provide an unprecedented opportunity for cost and time effective high quality microsatellite marker discovery in nonmodel organisms for which no genomic information is available. Here, we use shotgun pyrosequencing of a microsatellite-enriched library to develop, for the first time, microsatellite markers for Alnus glutinosa, a keystone tree species of European riparian woodland communities. From a total of 17?855 short sequences, we identified 590 perfect microsatellites from which 392 had designed primers. A subset of 48 loci were tested for amplification, 12 of which were polymorphic in A. glutinosa. These 12 loci were successfully coamplified in a single multiplex polymerase chain reaction experiment and validated for population genetics applications. In addition, 10 and 8 of these microsatellites were found to be transferable to the related A. incana and A. cordata species. The developed multiplex of 12 microsatellite markers therefore provides new opportunities for experimental evolutionary and forest genetics research in Alnus.     

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.     

Development of Genomic Microsatellite Markers in Carthamus tinctorius L. (Safflower) Using Next Generation Sequencing and Assessment of Their Cross-Species Transferability and Utility for Diversity Analysis

Background

Safflower (Carthamus tinctorius L.), an Asteraceae member, yields high quality edible oil rich in unsaturated fatty acids and is resilient to dry conditions. The crop holds tremendous potential for improvement through concerted molecular breeding programs due to the availability of significant genetic and phenotypic diversity. Genomic resources that could facilitate such breeding programs remain largely underdeveloped in the crop. The present study was initiated to develop a large set of novel microsatellite markers for safflower using next generation sequencing.

Principal Findings

Low throughput genome sequencing of safflower was performed using Illumina paired end technology providing ~3.5X coverage of the genome. Analysis of sequencing data allowed identification of 23,067 regions harboring perfect microsatellite loci. The safflower genome was found to be rich in dinucleotide repeats followed by tri-, tetra-, penta- and hexa-nucleotides. Primer pairs were designed for 5,716 novel microsatellite sequences with repeat length ≥ 20 bases and optimal flanking regions. A subset of 325 microsatellite loci was tested for amplification, of which 294 loci produced robust amplification. The validated primers were used for assessment of 23 safflower accessions belonging to diverse agro-climatic zones of the world leading to identification of 93 polymorphic primers (31.6%). The numbers of observed alleles at each locus ranged from two to four and mean polymorphism information content was found to be 0.3075. The polymorphic primers were tested for cross-species transferability on nine wild relatives of cultivated safflower. All primers except one showed amplification in at least two wild species while 25 primers amplified across all the nine species. The UPGMA dendrogram clustered C. tinctorius accessions and wild species separately into two major groups. The proposed progenitor species of safflower, C. oxyacantha and C. palaestinus were genetically closer to cultivated safflower and formed a distinct cluster. The cluster analysis also distinguished diploid and tetraploid wild species of safflower.

Conclusion

Next generation sequencing of safflower genome generated a large set of microsatellite markers. The novel markers developed in this study will add to the existing repertoire of markers and can be used for diversity analysis, synteny studies, construction of linkage maps and marker-assisted selection.     

Development of nine polymorphic microsatellite markers for the phytoparasitic nematode Xiphinema index, the vector of the grapevine fanleaf virus

We report isolation, characterization and cross-species amplification of nine microsatellite loci from the phytoparasitic nematode Xiphinema index, the vector of grapevine fanleaf virus. Levels of polymorphism were evaluated in 62 individuals from two X. index populations. The number of alleles varies between two and 10 depending on locus and population. Observed heterozygosity on loci across both populations varied from 0.32 to 0.857 (mean 0.545). The primers were tested for cross-species amplification in three other species of phytoparasitic nematodes of the Xiphinema genus. These nine microsatellite loci constitute valuable markers for population genetics and phylogeographical studies of X. index.     

Development of 13 microsatellite markers for <Emphasis Type="Italic">Castanopsis tribuloides</Emphasis> (Fagaceae) using next-generation sequencing

Catanopsis tribuloides is a climax tree species commonly distributed in evergreen forests and has been used to restore degraded areas in northern Thailand. To aid in study of genetic diversity of the species, microsatellite markers, which are specific to C. tribuloides, were developed using whole genome sequencing by next-generation sequencing technology. The primers for microsatellite were developed and screened for optimal annealing temperature by PCR assay. The loci primers specific with C. tribuloides, 13 polymorphic microsatellite primers were successfully developed. The results from genetic information analyzing showed the number of alleles presented were between 2 and 24. Accordingly, the expected and observed heterozygosity obtained were between 0.298 and 0.920 and 0.364 to 1.000, respectively. Null allele frequency was presented 0.000–0.199. Genetic information was generated 10 loci primers significantly deviated from Hardy–Weinberg Equilibrium. All 13 primer pairs of loci were not significant with linkage disequilibrium. A set of microsatellite markers in this study could be applied to gene flow, genetic structure and population genetic studies in the future.     

High-throughput microsatellite marker development in two sparid species and verification of their transferability in the family Sparidae

Recently, 454 sequencing has emerged as a popular method for isolating microsatellites owing to cost-effectiveness and time saving. In this study, repeat-enriched libraries from two southern African endemic sparids (Pachymetopon blochii and Lithognathus lithognathus) were 454 GS-FLX sequenced. From these, 7370 sequences containing repeats (SCRs) were identified. A brief survey of 23 studies showed a significant difference between the number of SCRs when enrichment was performed first before 454 sequencing. We designed primers for 302 unique fragments containing more than five repeat units and suitable flanking regions. A fraction (<11%) of these loci were characterized with 18 polymorphic microsatellite loci (nine in each of the focal species) being described. Sanger sequencing of alleles confirmed that size variation was because of differences in the number of tandem repeats. However, a case of homoplasy and sequencing errors in the 454 sequencing were identified. These newly developed and four previously isolated loci were successfully used to identify polymorphic markers in nine other economically important species, representative of sparid diversity. The combination of newly developed markers with data from previous sparid cross-species studies showed a significant negative correlation between genetic divergence to focal species and microsatellite transferability. The high level of transferability we described (48% amplification success and 32% polymorphism) suggests that the 302 microsatellite loci identified represent an excellent resource for future studies on sparids. Microsatellite marker development should commonly include tests of transferability to reduce costs and increase feasibility of population genetics studies in nonmodel organisms.     

Utility of pairwise mtDNA genetic distances for predicting cross-species microsatellite amplification and polymorphism success in fishes

In many studies involving microsatellites cross-species amplification, primers designed for one (source) species are used to amplify homologous loci in related (target) species. However, it is not clear how closely related the species must be to attain significant success. Genetic divergence is a clear and easy way to assess similarity between species and provides an accurate measure of their evolutionary distance. Eight Mediterranean target species of the family Serranidae were analysed using twelve primers developed for Serranus cabrilla. Additionally, two mitochondrial genes (12S rRNA and 16S rRNA) were chosen on the basis of their extensive use in phylogenetic and evolutionary analyses to compute genetic divergence between the species. Significant negative correlations were found between genetic divergence and both cross-species amplification and maintained polymorphism of microsatellite markers, which could be generalized by gathering information from different fish studies. The success of obtaining amplifiable and polymorphic microsatellite loci can be a priori approximated knowing the mtDNA genetic divergence between a given source and target species using our inferred regression equations. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Microsatellite markers for Lycium ruthenicum (Solananeae)

We developed microsatellite markers in Lycium ruthenicum, a desert plant widely distributed in northwestern China. In order to investigate its population genetic structure, genetic diversity, and its evolutionary history, we have isolated 11 novel microsatellite loci primers and characterized them in 24 individuals from 3 populations of L. ruthenicum using the combined biotin capture technique. For these microsatellites, one to seven alleles per locus were identified. The observed heterozygosities ranged from 0 to 0.958, meanwhile the expected heterozygosities ranged from 0 to 0.841. These microsatellite markers could be first useful for population level studies like genetic diversity and structure in this species.     

Microsatellite markers from Ceanothus roderickii (Rhamnaceae) using next-generation sequencing technology

? Premise of the study: Ceanothus roderickii is an endangered shrub endemic to California. To investigate the population genetics of this species, including the genetic consequences of population fragmentation and hybridization, 10 microsatellite markers were developed. ? Methods and Results: Using next-generation sequencing (454) data from a single C. roderickii individual, 10 microsatellite markers were developed. A group of 12 individuals representing all of the major C. roderickii populations were analyzed. All loci were found to be polymorphic, with a range from two to 12 alleles per locus. Observed heterozygosity ranged from 0.08 to 0.83 across loci. All 10 loci were also amplifiable in at least one other Ceanothus species. ? Conclusions: Results presented here indicate the utility of our new microsatellite primers in ongoing and future studies concerning population genetics and gene flow in C. roderickii, as well as the potential applicability of these primers in similar studies on other Ceanothus species.     

微卫星跨物种交叉PCR扩增的一个新问题:扩增非同源产物

微卫星已被广泛应用于群体遗传学、生态学和进化生物学研究。然而,一些物种微卫星尚未克隆。为了节省时间和经费,研究人员往往使用一个物种已发表的微卫星引物扩增其近缘物种的微卫星。该研究对属于3个不同科(Clariidae、Heteropneustidae 和Pimelodidae)的7个鲶鱼物种的微卫星跨物种PCR扩增产物进行了序列分析,研究发现扩增非同源（non-orthologous）产物是微卫星跨物种PCR扩增的一个新问题。该研究共采用4对胡子鲶微卫星座位引物对7个鲶鱼物种进行了跨物种PCR扩增。对获得的204个PCR产物的序列分析结果表明,两对微卫星座位引物扩增了所有7个物种的同源特异产物。而其他两个座位的引物扩增了特异但非同源的多态产物,对近缘物种的扩增也获得类似结果。另外,除胡子鲶等位基因大小异源同型(size homoplasy)的特征不明显外,其他物种在3个微卫星座位都具有这一非常明显的特征。这些数据表明,微卫星跨物种间交叉扩增能产生非同源产物;等位基因大小异源同型与微卫星座位本身有关,而与物种间的亲缘关系无明显的相关性。微卫星跨物种扩增产生的非同源产物和等位基因大小异源同型将使系统发育、群体遗传学和进化研究明显复杂化。因此,在应用微卫星跨物种交叉扩增数据以前,最好对跨物种交叉扩增产物进行测序验证。     

Polymorphic microsatellite loci from an indigenous Asian fungus-growing termite, Macrotermes gilvus (Blattodea: Termitidae) and cross amplification in related taxa

The fungus-growing termite, Macrotermes gilvus (Hagen), an indigenous species from Southeast Asia distributed from Myanmar to Indonesia and the Philippines, offers great potential as an ecological model system to elucidate the effects of geography on gene flow within this region. We used next generation sequencing (Roche 454 pyrosequencing) to identify microsatellite markers from the genomic DNA of M. gilvus. A modest sequencing volume generated 34,122 reads, with 1,212 (3.6%) reads contains microsatellites with di-, tri-, tetra-, penta-, and hexa-nucleotide repeat motifs. Thirty-seven loci were selected for primer development and tested for polymorphism across 22 colonies of M. gilvus. Eleven loci were found to be polymorphic with 2-4 alleles per locus. Observed and expected heterozygosities ranged between 0.091-0.727 and 0.090-0.540, respectively. Cross taxa amplification was successful across a panel of four related termite species and four multiplex groups were designed for future population genetic studies. These markers will open new avenues for the study of phylogeography and population genetics of this fungus-growing termite. This study also has effectively demonstrated the use of 454 pyrosequencing for the rapid development of informative microsatellite markers from a termite genome.     

Characterization of microsatellite markers in the threatened sand skink (Neoseps reynoldsi)

We characterized microsatellite loci for the sand skink (Neoseps reynoldsi) for future studies of genetic structure in this threatened taxon. We screened a partial genomic library enriched for microsatellites, designed primers for eight loci and assessed these markers for polymorphism across 11 populations in central Florida. Preliminary analyses indicate deviations from Hardy–Weinberg expectations for most loci, suggesting population genetic structure across the sampled populations; therefore, understanding genetic connectivity is critical for maintaining genetic variation in this species.     

Rapid identification of thousands of copperhead snake (Agkistrodon contortrix) microsatellite loci from modest amounts of 454 shotgun genome sequence

Optimal integration of next-generation sequencing into mainstream research requires re-evaluation of how problems can be reasonably overcome and what questions can be asked. One potential application is the rapid acquisition of genomic information to identify microsatellite loci for evolutionary, population genetic and chromosome linkage mapping research on non-model and not previously sequenced organisms. Here, we report on results using high-throughput sequencing to obtain a large number of microsatellite loci from the venomous snake Agkistrodon contortrix, the copperhead. We used the 454 Genome Sequencer FLX next-generation sequencing platform to sample randomly ∼27 Mbp (128 773 reads) of the copperhead genome, thus sampling about 2% of the genome of this species. We identified microsatellite loci in 11.3% of all reads obtained, with 14 612 microsatellite loci identified in total, 4564 of which had flanking sequences suitable for polymerase chain reaction primer design. The random sequencing-based approach to identify microsatellites was rapid, cost-effective and identified thousands of useful microsatellite loci in a previously unstudied species.     

Microsatellite loci for population and parentage analysis in the Amazon River dolphin (Inia geoffrensis de Blainville, 1817)

We developed specific primers for microsatellite DNA regions for the Amazon River dolphin or boto Inia geoffrensis, for use in population and conservation genetic studies. We also tested their transferability for two other species, Pontoporia blainvillei (sister taxon of I. geoffrensis) and Sotalia guianensis. A total of 12 microsatellite loci were polymorphic for the boto. An additional 25 microsatellite loci previously isolated from other cetacean species were also tested in the boto. The 26 polymorphic microsatellite loci indicate they will be excellent markers for studies of population structure and kinship relations of the boto.     

Transferability of microsatellite primers developed for stingless bees to four other species of the genus Melipona

Microsatellite markers are a useful tool for ecological monitoring of natural and managed populations. A technical limitation is the necessity for investment in the development of primers. Heterologous primers can provide an alternative to searching for new loci. In bees, these markers have been used in populational and intracolonial genetic analyses. The genus Melipona has the largest number of species among bee genera, about 70, occurring throughout the Neotropical region. However, only five species of the genus Melipona have specific microsatellite markers. Given the great diversity of this genus, this number is not representative. We analyzed the transferability of 49 microsatellite loci to four other species of the genus Melipona (M. scutellaris, M. mondury, M. mandacaia, and M. quadrifasciata). Four individuals of each species, from different localities, were used in amplification tests. Primer pairs described for five Melipona species and for Trigona carbonaria were tested. Among the 49 loci, 22 gave amplification products for all four species, while three gave nonspecific bands and five showed no amplification products. The remaining loci varied in the pattern of amplification, according to the species examined. The number of alleles ranged from 1 to 6. The results demonstrate the possibility of using these heterologous markers in other Melipona species, increasing the number of loci that can be analyzed and contributing to further genetic analyses of intra- and intercolonial structure, which is required for conservation measure planning, genetic improvement and resolution of taxonomic problems.