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.     

Lessons learned from microsatellite development for nonmodel organisms using 454 pyrosequencing

Microsatellites, also known as simple sequence repeats (SSRs), are among the most commonly used marker types in evolutionary and ecological studies. Next Generation Sequencing techniques such as 454 pyrosequencing allow the rapid development of microsatellite markers in nonmodel organisms. 454 pyrosequencing is a straightforward approach to develop a high number of microsatellite markers. Therefore, developing microsatellites using 454 pyrosequencing has become the method of choice for marker development. Here, we describe a user friendly way of microsatellite development from 454 pyrosequencing data and analyse data sets of 17 nonmodel species (plants, fungi, invertebrates, birds and a mammal) for microsatellite repeats and flanking regions suitable for primer development. We then compare the numbers of successfully lab‐tested microsatellite markers for the various species and furthermore describe diverse challenges that might arise in different study species, for example, large genome size or nonpure extraction of genomic DNA. Successful primer identification was feasible for all species. We found that in species for which large repeat numbers are uncommon, such as fungi, polymorphic markers can nevertheless be developed from 454 pyrosequencing reads containing small repeat numbers (five to six repeats). Furthermore, the development of microsatellite markers for species with large genomes was also with Next Generation Sequencing techniques more cost and time‐consuming than for species with smaller genomes. In this study, we showed that depending on the species, a different amount of 454 pyrosequencing data might be required for successful identification of a sufficient number of microsatellite markers for ecological genetic studies.     

Microsatellite markers for the endangered root holoparasite Dactylanthus taylorii (Balanophoraceae) from 454 pyrosequencing

? Premise of the study: Microsatellite loci were isolated and developed as polymorphic markers for the New Zealand endemic root holoparasite Dactylanthus taylorii for use in population and conservation genetics studies. ? Methods and Results: Shotgun 454 pyrosequencing was performed on genomic DNA pooled from three individuals of D. taylorii. From 61709 individual sequence reads, primers for 753 microsatellite loci were developed in silico and 72 of these were tested for consistent amplification and variability. Ten microsatellite loci were found to be polymorphic and consistently scorable when screened in 44 individuals from five geographically distant populations. The number of alleles per locus ranged from four to 16 with an average of 9.7, and average observed heterozygosity per locus was between 0.182 and 0.634. ? Conclusions: These polymorphic microsatellite markers establish an important resource for ongoing conservation initiatives and planned population genetic studies of D. taylorii.     

Utilizing next generation sequencing to characterize microsatellite loci in a tropical aquatic plant species Cryptocoryne cordata var. cordata (Araceae)

Cryptocoryne cordata var. cordata (2n = 34) is an aquatic plant species distributed from the southern part of Peninsular Thailand through the Malay Peninsula. It propagates both sexually and asexually via stolons. The current study is aimed at developing nuclear microsatellite markers for the species using next generation sequencing (Roche 454 pyrosequencing) from genomic DNA. A total of 41,653 reads was generated, of which 3636 fragments contained at least one repeat motif. Seventy two primer sets in the flanking region of dinucleotide, trinucleotide and tetranucleotides repeat motifs were designed and tested for efficiency in polymerase chain reaction amplification. Using these primer sets, 11 new microsatellite marker loci were successfully amplified with unambiguous polymorphic alleles exhibited across 30 individuals examined. The number of alleles per locus ranged from 2 to 6, while observed and expected heterozygosity ranged from 0.8190 to 1.0000 and 0.5401 to 0.7548, respectively. Genotype frequencies at all loci departed significantly from Hardy–Weinberg equilibrium. Linkage disequilibrium was not detected between any pair of loci. Cross-species amplification was successful across a panel of ten Cryptocoryne species. The markers described in this study will be useful for evaluating genetic diversity within and between populations, levels of gene flow, and the population dynamics of clones. They will be of further value in the development of effective conservation programs for Cryptocoryne species.     

Thousands of microsatellite loci from the venomous coralsnake Micrurus fulvius and variability of select loci across populations and related species

Studies of population genetics increasingly use next‐generation DNA sequencing to identify microsatellite loci in nonmodel organisms. There are, however, relatively few studies that validate the feasibility of transitioning from marker development to experimental application across populations and species. North American coralsnakes of the Micrurus fulvius species complex occur in the United States and Mexico, and little is known about their population structure and phylogenetic relationships. This absence of information and population genetics markers is particularly concerning because they are highly venomous and have important implications on human health. To alleviate this problem in coralsnakes, we investigated the feasibility of using 454 shotgun sequences for microsatellite marker development. First, a genomic shotgun library from a single individual was sequenced (approximately 7.74 megabases; 26 831 reads) to identify potentially amplifiable microsatellite loci (PALs). We then hierarchically sampled 76 individuals from throughout the geographic distribution of the species complex and examined whether PALs were amplifiable and polymorphic. Approximately half of the loci tested were readily amplifiable from all individuals, and 80% of the loci tested for variation were variable and thus informative as population genetic markers. To evaluate the repetitive landscape characteristics across multiple snakes, we also compared microsatellite content between the coralsnake and two other previously sampled snakes, the venomous copperhead (Agkistrodon contortrix) and Burmese python (Python molurus).     

Isolation and characterization of 20 polynucleotide microsatellite markers in a vulnerable Korean snail, <Emphasis Type="Italic">Ellobium chinense,</Emphasis> using 454 pyrosequencing

A small and air-breathing snail, Ellobium chinense (Ellobiidae), is a vulnerable species by International Union for Conservation of Nature (IUCN). To protect and manage habitat and population of E. chinense, microsatellite markers were developed using 454 pyrosequencing and 20 polymorphic microsatellite markers were identified. A total of 146,704 sequences containing a minimum of four repeat motifs (mean, 631 base pairs) were identified from 499,505 reads. Among 80 loci containing more than nine repeat units, 34 primer sets (42.5 %) produced strong PCR products, of which 20 were polymorphic among 48 samples of E. chinense. All loci exhibited high genetic variability, with an average of 18.9 alleles per locus, and the mean observed and expected heterozygosities were 0.65 and 0.90, respectively. In addition, cross-amplification was tested for all 20 loci in the same family species, Melampus sincaporensis. None of the primer pairs resulted in effective amplification, which might be due to their high mutation rates. Our work demonstrated the utility of next-generation 454 sequencing as a method for the rapid and cost-effective identification of microsatellites. The high degree of polymorphism exhibited by the 20 newly developed microsatellites will be useful in future conservation genetic studies of this species.     

Development of microsatellite markers in Fosterella rusbyi (Bromeliaceae) using 454 pyrosequencing

? Premise of the study: Polymorphic microsatellite markers were developed for Fosterella rusbyi (Bromeliaceae) to evaluate the population genetic structure and genetic diversity of natural populations of F. rusbyi and other Fosterella species in Bolivia. ? Methods and Results: 454 pyrosequencing technology was used to generate 73027 sequence reads from F. rusbyi DNA, which together contained 2796 perfect simple sequence repeats (SSRs). Primer pairs were designed for 30 loci, of which 15 were used to genotype 30 F. rusbyi plants from two geographical areas in Bolivia. All markers were polymorphic, with two to nine alleles in the overall sample. Cross-species amplification was tested in 10 additional Fosterella species. Seven loci showed consistent amplification in six or more species. ? Conclusions: The 15 SSR markers developed for F. rusbyi are promising candidates for population genetic analyses within F. rusbyi and other species of Fosterella.     

Development of 21 polymorphic microsatellite markers for the black-banded sea krait, <Emphasis Type="Italic">Laticauda semifasciata</Emphasis> (Elapidae: Laticaudinae), and cross-species amplification for two other congeneric species

The genus Laticauda (Reptilia: Elapidae), commonly known as sea kraits, is venomous marine amphibious snakes distributed throughout the south and southeast Asian islands and mostly found in coastal waters. To facilitate genetic studies, we have developed microsatellite loci for L. semifasciata using the 454 GS-FLX pyrosequencing technique. A total of 65,680 sequences containing a minimum of five repeat motifs were identified from 451,659 reads. Among 80 loci containing more than nine repeat units, 34 primer sets (42.5%) produced strong PCR products, of which 21 were polymorphic among 36 samples of L. semifasciata. All loci exhibited high genetic variability, with an average of 7.38 alleles per locus, and the mean observed and expected heterozygosities were 0.73 and 0.76, respectively. The cross-species amplification of these loci in two laticaudine species, L. colubrina and L. laticaudata, revealed a high transferability (78.6%) and polymorphism (59.5%) of the loci. Our work demonstrated the utility of next-generation 454 sequencing as the rapid and cost-effective method for development of microsatellite markers. The high level of polymorphism in these microsatellite loci will be useful for the detection of population subdivision and the study of migration, gene flow, relatedness and philopatry of L. semifasciata and other laticaudine species.     

Development of microsatellite markers for the endangered semi‐shrub Chimaphila umbellata (Ericaceae)

Based on the DNA sequencing reads obtained using 454 pyrosequencing, primers amplifying 16 microsatellite loci were developed for the endangered semi‐shrub Chimaphila umbellata, which occurs sporadically in the Japanese Archipelago. These 16 loci were polymorphic in the populations sampled from the Hokkaido and Tohoku Districts; the mean number of alleles was 3.31 and 3.44, and the mean expected heterozygosity was 0.42 and 0.44, respectively. These loci were not linked to each other and contained no null alleles. Amplification using these primers was also tested in the congeneric species C. japonica, but only three of them successfully amplified DNA of the species. These markers will be used to examine genetic diversity and genetic differentiation in populations of C. umbellata.     

Microsatellite development for the endangered Yangtze sturgeon (Acipenser dabryanus Duméril, 1869) using 454 sequencing

The Yangtze sturgeon (Acipenser dabryanus) is an endemic species in China. Using 454 sequencing, eight polymorphic tri‐, tetra‐, penta‐, and hexanucleotide microsatellite loci were isolated in this study. The raw sequence data from a one‐eighth run of 454 sequencings were 38.0 Mbp containing 94 222 reads/sequences. Of 80 microsatellite loci, only eight loci were polymorphic in a population of 30 individuals. The number of alleles per locus ranged from 4 to 14 (mean 7.62), and the observed heterozygosities varied between 0.46 and 0.88 (mean 0.74). Cross amplification was tested in congeneric species Acipenser sturio and Acipenser sinensis. These new microsatellite markers will be useful for further studies on genetic variation, parentage analysis, and conservation management for this critically endangered species.     

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.     

Isolation and characterization of microsatellite markers for Achyranthes bidentata (Amaranthaceae) using next-generation sequencing platform

Achyranthes bidentata Blume (Amaranthaceae) is a perennial herb that is widely distributed in India, Java, China, and Japan. The natural resources of A. bidentata in its geo-authentic product area have rapidly declined in recent years because of the over-collection of its roots. To devise adequate conservation and management strategies for this species, its genetic diversity and population structure should be characterized. Roche 454 pyrosequencing combined with magnetic bead enrichment was used to develop microsatellite markers for A. bidentata. A total of 903 microsatellite loci were identified from 42,004 individual sequence reads. One hundred microsatellite loci were selected to test the primer amplification efficiency across 16 individuals from two A. bidentata populations. Of these tested markers, 8 yielded polymorphic amplification products, 29 yielded single alleles. For polymorphic primer pairs, the number of alleles per locus ranged from 2 to 4, with an average of 2.75. The observed and expected heterozygosities ranged from 0.353 to 0.671 and 0.250 to 0.938, respectively. The inbreeding coefficient varied from −0.692 to 0.627. This set of markers will provide useful tools for examining genetic diversity and population structure, and aid in better understanding of the conservation of A. bidentata.     

MicNeSs: genotyping microsatellite loci from a collection of (NGS) reads

Microsatellites are widely used in population genetics to uncover recent evolutionary events. They are typically genotyped using capillary sequencer, which capacity is usually limited to 9, at most 12 loci for each run, and which analysis is a tedious task that is performed by hand. With the rise of next‐generation sequencing (NGS), a much larger number of loci and individuals are available from sequencing: for example, on a single run of a GS Junior, 28 loci from 96 individuals are sequenced with a 30X cover. We have developed an algorithm to automatically and efficiently genotype microsatellites from a collection of reads sorted by individual (e.g. specific PCR amplifications of a locus or a collection of reads that encompass a locus of interest). As the sequencing and the PCR amplification introduce artefactual insertions or deletions, the set of reads from a single microsatellite allele shows several length variants. The algorithm infers, without alignment, the true unknown allele(s) of each individual from the observed distributions of microsatellites length of all individuals. MicNeSs, a python implementation of the algorithm, can be used to genotype any microsatellite locus from any organism and has been tested on 454 pyrosequencing data of several loci from fruit flies (a model species) and red deers (a nonmodel species). Without any parallelization, it automatically genotypes 22 loci from 441 individuals in 11 hours on a standard computer. The comparison of MicNeSs inferences to the standard method shows an excellent agreement, with some differences illustrating the pros and cons of both methods.     

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.     

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.     

Genome-Wide Microsatellite Identification in the Fungus Anisogramma anomala Using Illumina Sequencing and Genome Assembly

High-throughput sequencing has been dramatically accelerating the discovery of microsatellite markers (also known as Simple Sequence Repeats). Both 454 and Illumina reads have been used directly in microsatellite discovery and primer design (the “Seq-to-SSR” approach). However, constraints of this approach include: 1) many microsatellite-containing reads do not have sufficient flanking sequences to allow primer design, and 2) difficulties in removing microsatellite loci residing in longer, repetitive regions. In the current study, we applied the novel “Seq-Assembly-SSR” approach to overcome these constraints in Anisogramma anomala. In our approach, Illumina reads were first assembled into a draft genome, and the latter was then used in microsatellite discovery. A. anomala is an obligate biotrophic ascomycete that causes eastern filbert blight disease of commercial European hazelnut. Little is known about its population structure or diversity. Approximately 26 M 146 bp Illumina reads were generated from a paired-end library of a fungal strain from Oregon. The reads were assembled into a draft genome of 333 Mb (excluding gaps), with contig N₅₀ of 10,384 bp and scaffold N₅₀ of 32,987 bp. A bioinformatics pipeline identified 46,677 microsatellite motifs at 44,247 loci, including 2,430 compound loci. Primers were successfully designed for 42,923 loci (97%). After removing 2,886 loci close to assembly gaps and 676 loci in repetitive regions, a genome-wide microsatellite database of 39,361 loci was generated for the fungus. In experimental screening of 236 loci using four geographically representative strains, 228 (96.6%) were successfully amplified and 214 (90.7%) produced single PCR products. Twenty-three (9.7%) were found to be perfect polymorphic loci. A small-scale population study using 11 polymorphic loci revealed considerable gene diversity. Clustering analysis grouped isolates of this fungus into two clades in accordance with their geographic origins. Thus, the “Seq-Assembly-SSR” approach has proven to be a successful one for microsatellite discovery.     

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.     

Intracolony variation of bacterial gut microbiota among castes and ages in the fungus-growing termite Macrotermes gilvus

The fungus-growing termites Macrotermes cultivate the obligate ectosymbiontic fungi, Termitomyces. While their relationship has been extesively studied, little is known about the gut bacterial symbionts, which also presumably play a crucial role for the nutrition of the termite host. In this study, we investigated the bacterial gut microbiota in two colonies of Macrotermes gilvus, and compared the diversity and community structure of bacteria among nine termite morphotypes, differing in caste and/or age, using terminal restriction fragment length polymorphism (T-RFLP) and clonal analysis of 16S rRNA. The obtained molecular community profiles clustered by termite morphotype rather than by colony, and the clustering pattern was clearly more related to a difference in age than to caste. Thus, we suggest that the bacterial gut microbiota change in relation to the food of the termite, which comprises fallen leaves and the fungus nodules of Termitomyces in young workers, and leaves degraded by the fungi, in old workers. Despite these intracolony variations in bacterial gut microbiota, their T-RFLP profiles formed a distinct cluster against those of the fungus garden, adjacent soil and guts of sympatric wood-feeding termites, implying a consistency and uniqueness of gut microbiota in M. gilvus. Since many bacterial phylotypes from M. gilvus formed monophyletic clusters with those from distantly related termite species, we suggest that gut bacteria have co-evolved with the termite host and form a microbiota specific to a termite taxonomic and/or feeding group, and furthermore, to caste and age within a termite species.     

Novel Polymorphic Microsatellite Markers for Panulirus ornatus and their Cross-species Primer Amplification in Panulirus homarus

Polymorphic microsatellite loci were isolated for Panulirus ornatus using 454 GS-FLX Titanium pyrosequencing. Fifteen markers containing perfect di-, tri-, tetra-, and penta-nucleotide motifs were consistently co-amplified in five multiplexes in a panel of 91 randomly selected samples. Observed number of alleles varied from 2 to 14 per locus. Observed and expected heterozygosity ranged from 0.090 to 0.79 and 0.08 to 0.87, respectively. Ten loci deviated from Hardy-Weinberg equilibrium after sequential Bonferroni correction. Genetic linkage disequilibrium analysis between all pairs of the loci showed significant departure from the null hypothesis between 11 loci. The microsatellite markers were also amplified successfully in related Panulirus homarus species with adequate level of polymorphism. The successful cross-species primer amplification of the 15 microsatellites indicates the potential of the developed markers to be transferred to other Panulirus species. The 15 novel microsatellite markers reported in this work add to the previously characterized markers by our group, exhibit adequate levels of polymorphism for wide range of future studies investigating population structure, genetic diversity, and evolutionary relationships among Panulirus species.