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.     

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.     

A set of microsatellite markers for population genetics of leopard cat (Prionailurus bengalensis) and cross-species amplification in other felids

We describe the isolation and characterization of novel microsatellite loci from the leopard cat, Prionailurus bengalensis Kerr, 1792 (Family Felidae). Using Illumina HiSeq2500 sequencing technology, we sequenced the leopard cat genome and identified 1.5 million loci of simple sequence repeats with di- to deca-nucleotide motifs. We developed twelve polymorphic markers with tetra-nucleotide motif types after screening 35 loci for amplification and polymorphism. The observed and expected heterozygosities of the markers were 0.438 and 0.423, respectively. The number of alleles per locus ranged from 2 to 7, with a mean polymorphism information content of 0.383. Eleven loci were at Hardy-Weinberg equilibrium and no linkage disequilibrium was detected among any pairs of loci. We tested cross-species amplification of these markers across five other felids (Panthera tigris, P. pardus, P. onca, Acinonyx jubatus, and Felis catus). All loci were transferable to at least one other feline species and four amplified all five species. The microsatellite markers developed in this study will be valuable for estimating ecological parameters of populations and to establish conservation and management strategies for feline species.     

Characterization of 13 polymorphic microsatellite loci in the Japanese land leech

We developed 13 polymorphic microsatellite loci of the Japanese land leech (Haemadipsa japonica; Haemadipsidea) using an Illumina MiSeq sequencing approach. A total of 42,064 nuclear DNA contigs were filtered for microsatellite motifs, among which 30,873 simple sequence repeat loci were identified. From these sequences, we selected 30 primer sets, and 13 of these loci were successfully amplified. Polymorphism of the 13 loci was tested using 16 individuals sampled from sixteen populations across Japan. The number of alleles and polymorphism information content varied from 5 to 17 and 0.335 to 0.883, respectively, and observed and expected heterozygosity values ranged from 0.143 to 0.875 and 0.349 to 0.893, respectively, indicating that these loci are polymorphic. Furthermore, we established useful multiplex PCR using these loci. The 13 microsatellite loci described in this paper are the first nuclear microsatellite markers for a land leech species.     

Two highly informative dinucleotide SSR multiplexes for the conifer Larix decidua (European larch)

We have designed two highly polymorphic microsatellite multiplexes for Larix decidua Mill (European larch), a coniferous tree species with a fragmented distribution across Europe. The multiplexes combine microsatellites previously designed for the sister species L. kaempferi and newly identified microsatellites obtained by pyrosequencing of an enriched microsatellite library and subsequent marker candidate selection. As we wanted to target highly polymorphic markers, only microsatellite motifs with a high number of repeats (≥ 12) were selected. An important proportion of the marker candidates presented multiple bands, bad amplification or insufficient polymorphism. Such difficulties were expected owing to the large genome size of the studied species. Our strategy for marker validation followed most recent recommendations for microsatellite development, for example verifying marker quality in terms of polymorphism and accurate allele binning before multiplexing. The most promising loci were combined in two multiplexes, a 7-plex and a 6-plex. These were tested on a sample of 413 individuals from 18 populations distributed across the natural range. The 13 loci had from 9 to 36 alleles. Markers were successfully tested in another laboratory, confirming robustness of the marker protocols. We also tested transferability on six other larch species from Asia and North America. Overall, this study shows that, even in species with large genome size and relatively low overall polymorphism, microsatellites can be successfully developed using next-generation sequencing technologies, provided that some additional precautions are taken compared to species lacking these characteristics.     

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.     

MICROSATELLITE DEVELOPMENT IN RHODOPHYTA USING HIGH‐THROUGHPUT SEQUENCE DATA1

Shotgun genome sequencing is rapidly emerging as the method of choice for the identification of microsatellite loci in nonmodel organisms. However, to the best of our knowledge, this approach has not been applied to marine algae so far. Herein, we report the results of using the 454 next‐generation sequencing (NGS) platform to randomly sample 36.0 and 40.9 Mbp (139,786 and 139,795 reads, respectively) of the genome of two red algae from the northwest Iberian Peninsula [Grateloupia lanceola (J. Agardh) J. Agardh and a still undescribed new member of the family Cruoriaceae]. Using data mining tools, we identified 4,766 and 5,174 perfect microsatellite loci in 4,344 and 4,504 sequences/contigs from G. lanceola and the Cruoriaceae, respectively. After conservative removal of potentially problematic loci (redundant sequences, mobile elements), primer design was possible for 1,371 and 1,366 loci, respectively. A survey of the literature indicates that microsatellite density in our Rhodophyta is at the low end of the values reported for other organisms investigated with the same technology (land plants and animals). A limited number of loci were successfully tested for PCR amplification and polymorphism finding that they may be suitable for population genetic studies. This study demonstrates that random genome sequencing is a rapid, effective alternative to develop useful microsatellite loci in previously unstudied red algae.     

濒危植物羊踯躅全基因组SSR标记开发与鉴定研究

该研究利用基于全基因组限制性酶切位点简化基因组测序技术（RAD seq技术）,开发濒危植物羊踯躅（Rhododendron molle G. Don）全基因组SSR标记,并对3个群体共63份羊踯躅材料进行验证鉴定,为进一步研究羊踯躅的遗传多样性和群体遗传结构以及保护利用提供技术支持。结果显示：（1）羊踯躅基因组测序获得原始数据7.653G bp,过滤后为7.513G bp;经组装发现,羊踯躅171.534 M bp的基因组分布在498 252 contigs中。（2）通过SSR检测,在11 961 SSR位点中获得了11 687对SSR分子标记,并且二核苷酸为基序的重复类型最丰富,达51.76%。（3）随机选取128对SSR标记在6个羊踯躅株系中进行PCR扩增,获得20对高多态性的SSR标记。（4）用所选的20对多态性SSR标记对3个群体共63份羊踯躅材料进行验证分析发现,这些多态性SSR标记位点的等位基因数为4~16个,期望杂合度（H_e）为0.489~0.908。 研究表明,羊踯躅的SSR丰度适中,且二核苷酸为羊踯躅中最丰富的重复序列,该实验进一步证明RAD seq技术是一种经济有效的基因测序方法,实验中开发的SSR引物将有助于进一步研究羊踯躅和其他近缘种的群体结构和多样性。     

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 genomic simple sequence repeat markers for linseed using next-generation sequencing technology

Linseed (Linum usitatissimum L.) is regarded as a cash crop of tomorrow because of the presence of nutraceutically important ??-linolenic acid (ALA) and lignan. However, only limited breeding progress has been made in this crop, mainly due to the lack of sufficient genetic and genomic resources. Among these, simple sequence repeats (SSR) are useful DNA markers for diversity analysis, genetic mapping and tagging traits because of their co-dominant and highly polymorphic nature. In order to develop SSR markers for linseed, we used three microsatellite isolation methods, viz., PCR Isolation of Microsatellite Arrays (PIMA), 5??-anchored PCR method, and Fast Isolation by AFLP of Sequences COntaining repeats (FIASCO). The amplified products from these methods were pooled and sequenced using the 454 GS-FLX platform. A total of 36,332 reads were obtained, which assembled into 2,183 contigs and 2,509 singlets. The contigs and the singlets contained 1,842 microsatellite motifs, with dinucleotide motifs as the most abundant repeat type (54%) followed by trinucleotide motifs (44%). Based on this, 290 SSR markers were designed, 52 of which were evaluated using a panel of 27 diverse linseed genotypes. Among the three enrichment methods, the 5??-anchored PCR method was most efficient for isolation of microsatellites, while FIASCO was most efficient for developing SSR markers. We show the utility of next-generation sequencing technology for efficiently discovering a large number of microsatellite markers in non-model plants.     

Characterization of sequence polymorphisms from microsatellite flanking regions in <Emphasis Type="Italic">Vitis</Emphasis> spp

Single nucleotide polymorphisms or SNPs are the most abundant form of genetic variation in the genome of plants and animals. Microsatellites are hypervariable regions of genome, while their flanking regions are assumed to be as conserved as the average of the genome. In the present study, flanking sequences of 10 microsatellite loci were compared in different cultivars of Vitis to determine the existing polymorphism. For every microsatellite, about 8 homozygous cultivars (regarding the microsatellite genotype) were chosen for sequencing. A total of 45 different varieties of Vitis and 91 sequences were analysed. Sequence polymorphisms were detected for all the microsatellite flanking regions studied, including single nucleotide polymorphisms (SNPs), insertions and deletions. The number of identified changes varied considerably among the loci with a frequency of one polymorphism every 41 nucleotides, being VVMD5 the most polymorphic one. A number of SNPs were used to design SNP markers, which were scored by dideoxy single base primer extension and capillary electrophoresis methodology. These SNP markers were employed to genotype 21 cultivars of Vitis vinifera and 4 varieties of other Vitis species. The utility of the markers developed as well as their utility for varietal identification and pedigree studies is discussed, using a similar study carried out with the 10 microsatellites as a reference.     

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

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

What remains from a 454 run: estimation of success rates of microsatellite loci development in selected newt species (Calotriton asper,Lissotriton helveticus,and Triturus cristatus) and comparison with Illumina‐based approaches

The development of microsatellite loci has become more efficient using next‐generation sequencing (NGS) approaches, and many studies imply that the amount of applicable loci is large. However, few studies have sought to quantify the number of loci that are retained for use out of the thousands of sequence reads initially obtained. We analyzed the success rate of microsatellite loci development for three amphibian species using a 454 NGS approach on tetra‐nucleotide motif‐enriched species‐specific libraries. The number of sequence reads obtained differed strongly between species and ranged from 19,562 for Triturus cristatus to 55,626 for Lissotriton helveticus, with 52,075 reads obtained for Calotriton asper. PHOBOS was used to identify sequences with tetra‐nucleotide repeat motifs with a minimum repeat number of ten and high quality primer binding sites. Of 107 sequences for T. cristatus, 316 for C. asper and 319 for L. helveticus, we tested the amplification success, polymorphism, and degree of heterozygosity for 41 primer combinations each for C. asper and T. cristatus, and 22 for L. helveticus. We found 11 polymorphic loci for T. cristatus, 20 loci for C. asper, and 15 loci for L. helveticus. Extrapolated, the number of potentially amplifiable loci (PALs) resulted in estimated species‐specific success rates of 0.15% (T. cristatus), 0.30% (C. asper), and 0.39% (L. helveticus). Compared with representative Illumina NGS approaches, our applied 454‐sequencing approach on specifically enriched sublibraries proved to be quite competitive in terms of success rates and number of finally applicable loci.     

Identification of the most represented repeated motifs in Arabidopsis thaliana microsatellite loci

The major simple sequence repeats present in the Arabidopsis genome were identified by Southern hybridizations with 49 oligonucleotide probes matching all the possible combinations of motifs up to 4 nucleotides long. The method used allowed us to perform all the hybridizations under the same temperature conditions. A good correlation was observed with the data obtained from database analysis, indicating that the method can be useful for identifying the major classes of microsatellite loci in species for which few or no sequence data are available. AG/CT, AAG/CTT, ATG/CAT and GTG/CAC are the major motifs present in the Arabidopsis genome that can be used as convenient probes to isolate microsatellite loci by screening libraries. AAG/CTT is the more frequent of these motifs, and its relative frequency in Arabidopsis is much higher than averagely found in the plant kingdom. About 8% of the cDNA clones from an immature silique library contains AG/CT, AAG/CTT or ATG/CAT microsatellite loci. Several microsatellite loci were isolated by screening genomic and cDNA libraries. Twenty-six tri-nucleotide loci were PCR amplified from four different ecotypes, and polymorphism was observed for 12 of them; 10 loci showing two alleles and 2 loci showing three alleles.     

Development of microsatellite markers specific for the short arm of rye (Secale cereale L.) chromosome 1

We developed 74 microsatellite marker primer pairs yielding 76 polymorphic loci, specific for the short arm of rye chromosome 1R (1RS) in wheat background. Four libraries enriched for microsatellite motifs AG, AAG, AC and AAC were constructed from DNA of flow-sorted 1RS chromosomes and 1,290 clones were sequenced. Additionally, 2,778 BAC-end-sequences from a 1RS specific BAC library were used for microsatellite screening and marker development. From 724 designed primer pairs, 119 produced 1RS specific bands and 74 of them showed polymorphism in a set of ten rye genotypes. We show that this high attrition rate was due to the highly repetitive nature of the rye genome consisting of a large number of transposable elements. We mapped the 76 polymorphic loci physically into three regions (bins) on 1RS; 29, 30 and 17 loci were assigned to the distal, intercalary and proximal regions of the 1RS arm, respectively. The average polymorphism information content increases with distance from the centromere, which could be due to an increased recombination rate along the chromosome arm toward’s the telomere. Additionally, we demonstrate, using the data of the whole rice genome, that the intra-genomic length variation of microsatellites correlates (r = 0.87) with microsatellite polymorphism. Based on these results we suggest that an analysis of the microsatellite length variation is conducted for each species prior to microsatellite development, provided that sufficient sequence information is available. This will allow to selectively design microsatellite markers for motifs likely to yield a high level of polymorphism. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Development and characterization of 12 polymorphic microsatellite loci in the sea sandwort, <Emphasis Type="Italic">Honckenya peploides</Emphasis>

Codominant marker systems are better suited to analyze population structure and assess the source of an individual in admixture analyses. Currently, there is no codominant marker system using microsatellites developed for the sea sandwort, Honckenya peploides (L.) Ehrh., an early colonizer in island systems. We developed and characterized novel microsatellite loci from H. peploides, using reads collected from whole genome shotgun sequencing on a 454 platform. The combined output from two shotgun runs yielded a total of 62,669 reads, from which 58 loci were screened. We identified 12 polymorphic loci that amplified reliably and exhibited disomic inheritance. Microsatellite data were collected and characterized for the 12 polymorphic loci in two Alaskan populations of H. peploides: Fossil Beach, Kodiak Island (n?=?32) and Egg Bay, Atka Island (n?=?29). The Atka population exhibited a slightly higher average number of alleles (3.9) and observed heterozygosity (0.483) than the Kodiak population (3.3 and 0.347, respectively). The overall probability of identity values for both populations was PID?=?2.892e^?6 and PID_sib?=?3.361e^?3. We also screened the 12 polymorphic loci in Wilhelmsia physodes (Fisch. ex Ser.) McNeill, the most closely related species to H. peploides, and only one locus was polymorphic. These microsatellite markers will allow future investigations into population genetic and colonization patterns of the beach dune ruderal H. peploides on new and recently disturbed islands.