Rapid microsatellite isolation from a butterfly by de novo transcriptome sequencing: performance and a comparison with AFLP-derived distances

High migratory propensity constitutes a single stock of an exploited cutlassfish species in the Northwest Pacific: A microsatellite approach

Cutlassfishes, also known as hairtails, include multiple predatory fishes of the family Trichiuridae. They constitute a top marine fish commodity globally, yet the knowledge about their composition and intraspecific genetic structures is still limited. Trichiurus japonicus accounts for a major amount in the northwest Pacific fishery. Previous studies based on mitochondrial DNA markers reported incongruences in its population structure, hence prompting the need for high-resolution markers and avoiding possible shortcomings in its management. Here we genotyped ten novel de novo-assembled transcriptome-derived microsatellite markers on a total of 150 samples across five major fishing grounds (encompassing latitudes 22-39°N). These markers presented a high number of alleles and heterozygosity compared to other marine fishes, corresponding to the large effective population size of ~20,000 per location and cohort differentiation. Population structuring analyses suggested T. japonicus to be a homogenous well-mixed population. This configuration is likely attributed to the majority of its effective population migrates across locations, and the absence of oceanographic barriers at the continental shelves. Qingdao with reportedly high ocean productivity could be a genetic pseudosink based on the high heterozygosity and migratory preference. Moreover, the results of sign tests suggest that T. japonicus experienced a recent bottleneck likely concurrent with historical glaciation events. Further, we demonstrated satisfactory cross-amplifications of our markers on several congeners, indicating a great promise to use these markers to study the population genetics of trichiurids. Together, our findings will serve as an essential groundwork for enhancing resource conservation and management of cutlassfishes.

Genetic diversity and structure of Chinese grass shrimp, Palaemonetes sinensis, inferred from transcriptome-derived microsatellite markers

BACKGROUND

The Chinese grass shrimp, Palaemonetes sinensis, is an economically important freshwater shrimp in China, and the study of genetic diversity and structure can positively contribute to the exploration of germplasm resources and assist in the understanding of P. sinensis aquaculture. Microsatellite markers are widely used in research of genetic backgrounds since it is considered an important molecular marker for the analyses of genetic diversity and structure. Hence, the aim of this study was to evaluate the genetic diversity and structure of wild P. sinensis populations in China using the polymorphic microsatellite makers from the transcriptome.

RESULTS

Sixteen polymorphic microsatellite markers were developed for P. sinensis from transcriptome, and analyzed for differences in genetic diversity and structure in multiple wild P. sinensis populations in China. Totally of 319 individual shrimps from seven different populations were genotyped to find that allelic polymorphisms varied in two to thirteen alleles seen in the entire loci. Compared to other populations analyzed, the two populations including LD and SJ showed lower genetic diversity. Both the genetic distance (D) and Wrights fixation index (F_ST) comparing any two populations also indicated that LD and SJ populations differed from the other five populations. An UPGMA tree analysis showed three main clusters containing SJ, LD and other populations which were also confirmed using STRUCTURE analysis.

CONCLUSION

This is the first study where polymorphic microsatellite markers from the transcriptome were used to analyze genetic diversity and structures of different wild P. sinensis populations. All the polymorphic microsatellite makers are believed useful for evaluating the extent of the genetic diversity and population structure of P. sinensis. Compared to the other five populations, the LD and SJ populations exhibited lower genetic diversity, and the genetic structure was differed from the other five populations. Therefore, they needed to be protected against further declines in genetic diversity. The other five populations, LP, LA, LSL, LSY and LSH, are all belonging to Liaohe River Drainage with a relatively high genetic diversity, and hence can be considered as hot spots for in-situ conservation of P. sinensis as well as sources of desirable alleles for breeding values.

New microsatellites revealed strong gene flow among populations of a new outbreak pest, Athetis lepigone (Möschler)

Athetis lepigone (Möschler) (Lepidoptera: Noctuidae) is a new outbreak pest in China. Consequently, it is unclear whether the emergence and spread of the outbreak of this pest are triggered by rapid in situ population size increases in each outbreak area, or by immigrants from a potential source area in China. In order to explore the outbreak process of this pest through a population genetics approach, we developed ten novel polymorphic expressed sequence tags (EST)-derived microsatellites. These new microsatellites had moderately high levels of polymorphism in the tested population. The number of alleles per locus ranged from 3 to 19, with an average of 8.6, and the expected heterozygosity ranged from 0.269 to 0.783. A preliminary population genetic analysis using these new microsatellites revealed a lack of population genetic structure in natural populations of A. lepigone. The estimates of recent migration rate revealed strong gene flow among populations. In conclusion, our study developed the first set of EST-microsatellite markers and shed a new light on the population genetic structure of this pest in China.

Development of Microsatellite Markers for the Nipa Palm Hispid Beetle, Octodonta nipae (Maulik)

The nipa palm hispid beetle, Octodonta nipae (Maulik) (Coleoptera: Chrysomelidae), is an important invasive pest on palm plants in southern China. Based on existing transcriptome data, polymorphism simple sequence repeat (SSR) loci were identified. In total, 1274 SSR loci were identified from 49919 unigenes. The majority of them contained mononucleotide, dinucleotide, and trinucleotide motifs (43.56%, 26.14%, and 28.18%), in which A/T (41.21%) and AT/TA (15.86%) were the most abundant motifs. 104 pairs of the SSR primers produced amplification bands of expected sizes in O. nipae, 80 pairs of SSR primers were tested randomly for polymorphism, 9 loci of them were validated to be polymorphic markers, and the number of alleles ranged from 2 to 3, with an average of 2.56 per locus. The population of Zhangzhou and Fuzhou was analyzed by the 9 loci (On1-On9). These SSR transcriptome data can provide invaluable resource for SSR development, population genetics research, invasion and expansion mechanism, paternity testing, and other research on O. nipae and its related species.

Transcriptome sequencing and de novo characterization of Korean endemic land snail, Koreanohadra kurodana for functional transcripts and SSR markers

The Korean endemic land snail Koreanohadra kurodana (Gastropoda: Bradybaenidae) found in humid areas of broadleaf forests and shrubs have been considered vulnerable as the number of individuals are declining in recent years. The species is poorly characterized at the genomic level that limits the understanding of functions at the molecular and genetics level. In the present study, we performed de novo transcriptome sequencing to produce a comprehensive transcript dataset of visceral mass tissue of K. kurodana by the Illumina paired-end sequencing technology. Over 234 million quality reads were assembled to a total of 315,924 contigs and 191,071 unigenes, with an average and N50 length of 585.6 and 715 bp and 678 and 927 bp, respectively. Overall, 36.32 % of the unigenes found matches to known protein/nucleotide sequences in the public databases. The direction of the unigenes to functional categories was determined using COG, GO, KEGG, and InterProScan protein domain search. The GO analysis search resulted in 22,967 unigenes (12.02 %) being categorized into 40 functional groups. The KEGG annotation revealed that metabolism pathway genes were enriched. The most prominent protein motifs include the zinc finger, ribonuclease H, reverse transcriptase, and ankyrin repeat domains. The simple sequence repeats (SSRs) identified from >1 kb length of unigenes show a dominancy of dinucleotide repeat motifs followed with tri- and tetranucleotide motifs. A number of unigenes were putatively assessed to belong to adaptation and defense mechanisms including heat shock proteins 70, Toll-like receptor 4, AMP-activated protein kinase, aquaporin-2, etc. Our data provide a rich source for the identification and functional characterization of new genes and candidate polymorphic SSR markers in K. kurodana. The availability of transcriptome information ( http://bioinfo.sch.ac.kr/submission/ ) would promote the utilization of the resources for phylogenetics study and genetic diversity assessment.

Isolation of microsatellite loci and reliable genotyping using noninvasive samples of a critically endangered primate, Trachypithecus leucocephalus

Genetic information can be critical in identifying conservation priorities and developing conservation strategies. There is an urgent need for noninvasive genetic tools to study the wild populations of Asian colobine monkeys. The majority of these species are threatened with habitat destruction, population reduction and even extinction, but generally lack information on their genetic diversity and population structure. Genetic sampling and tissue collection have been scarce in these species owing to strict regulations on manipulation of endangered species, and the difficulties and risks associated with capturing these arboreal and fast-moving monkeys in the challenging environments that they inhabit. These difficulties have hindered the development of molecular genetic markers, which are usually derived from tissues or blood. In this study, we present a method for de novo microsatellite isolation and genotyping using DNA from noninvasive origins of a critically endangered Asian colobine, the white-headed langur (Trachypithecus leucocephalus). Genomic DNA isolated from hair was shown to be sufficient for microsatellite enrichment and isolation, with similar isolation efficiencies as from tissue DNA. We identified and characterized 20 polymorphic microsatellite loci, and evaluated their amplification success and genotyping reliability with 86 field-collected fecal samples. These results show that this panel of loci can produce reliable genotypes from fecal samples, and represent a useful tool for noninvasive investigation of genetic structure, individual identification and kinship assessment in this highly endangered species. Our approach can be applied to conservation genetic studies of other wild species that lack sequence information and tissue samples.

Transcriptome Profile of the Asian Giant Hornet (Vespa mandarinia) Using Illumina HiSeq 4000 Sequencing: De Novo Assembly, Functional Annotation, and Discovery of SSR Markers

Vespa mandarinia found in the forests of East Asia, including Korea, occupies the highest rank in the arthropod food web within its geographical range. It serves as a source of nutrition in the form of Vespa amino acid mixture and is listed as a threatened species, although no conservation measures have been implemented. Here, we performed de novo assembly of the V. mandarinia transcriptome by Illumina HiSeq 4000 sequencing. Over 60 million raw reads and 59,184,811 clean reads were obtained. After assembly, a total of 66,837 unigenes were clustered, 40,887, 44,455, and 22,390 of which showed homologous matches against the PANM, Unigene, and KOG databases, respectively. A total of 15,675 unigenes were assigned to Gene Ontology terms, and 5,132 unigenes were mapped to 115 KEGG pathways. The zinc finger domain (C2H2-like), serine/threonine/dual specificity protein kinase domain, and RNA recognition motif domain were among the top InterProScan domains predicted for V. mandarinia sequences. Among the unigenes, we identified 534,922 cDNA simple sequence repeats as potential markers. This is the first transcriptomic analysis of the wasp V. mandarinia using Illumina HiSeq 4000. The obtained datasets should promote the search for new genes to understand the physiological attributes of this wasp.

Microsatellites in Pursuit of Microbial Genome Evolution

Microsatellites or short sequence repeats are widespread genetic markers which are hypermutable 1-6 bp long short nucleotide motifs. Significantly, their applications in genetics are extensive due to their ceaseless mutational degree, widespread length variations and hypermutability skills. These features make them useful in determining the driving forces of evolution by using powerful molecular techniques. Consequently, revealing important questions, for example, what is the significance of these abundant sequences in DNA, what are their roles in genomic evolution? The answers of these important questions are hidden in the ways these short motifs contributed in altering the microbial genomes since the origin of life. Even though their size ranges from 1 -to- 6 bases, these repeats are becoming one of the most popular genetic probes in determining their associations and phylogenetic relationships in closely related genomes. Currently, they have been widely used in molecular genetics, biotechnology and evolutionary biology. However, due to limited knowledge; there is a significant gap in research and lack of information concerning hypermutational mechanisms. These mechanisms play a key role in microsatellite loci point mutations and phase variations. This review will extend the understandings of impacts and contributions of microsatellite in genomic evolution and their universal applications in microbiology.

Transcriptome Characterization for Non-Model Endangered Lycaenids, Protantigius superans and Spindasis takanosis, Using Illumina HiSeq 2500 Sequencing

The Lycaenidae butterflies, Protantigius superans and Spindasis takanosis, are endangered insects in Korea known for their symbiotic association with ants. However, necessary genomic and transcriptomics data are lacking in these species, limiting conservation efforts. In this study, the P. superans and S. takanosis transcriptomes were deciphered using Illumina HiSeq 2500 sequencing. The P. superans and S. takanosis transcriptome data included a total of 254,340,693 and 245,110,582 clean reads assembled into 159,074 and 170,449 contigs and 107,950 and 121,140 unigenes, respectively. BLASTX hits (E-value of 1.0 × 10⁻⁵) against the known protein databases annotated a total of 46,754 and 51,908 transcripts for P. superans and S. takanosis. Approximately 41.25% and 38.68% of the unigenes for P. superans and S. takanosis found homologous sequences in Protostome DB (PANM-DB). BLAST2GO analysis confirmed 18,611 unigenes representing Gene Ontology (GO) terms and a total of 5259 unigenes assigned to 116 pathways for P. superans. For S. takanosis, a total of 6697 unigenes were assigned to 119 pathways using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. Additionally, 382,164 and 390,516 Simple Sequence Repeats (SSRs) were compiled from the unigenes of P. superans and S. takanosis, respectively. This is the first report to record new genes and their utilization for conservation of lycaenid species population and as a reference information for closely related species.

De novo sequencing and characterization of the Bradysia odoriphaga (Diptera: Sciaridae) larval transcriptome

The most serious pestilent threat to the Chinese chive, Allium tuberosum Rottle ex Spreng (Liliaceae) is the Bradysia odoriphaga Yang and Zhang. There is limited genetic research focused on B. odoriphaga, partially due to the lack of genomic resources. The advent of high-throughput sequencing technologies has enabled generation of genomic resources in a short time frame and at minimal costs. In this study, we performed, for the first time, de novo transcriptome sequencing of the B. odoriphaga. Here, 16,829 unigenes were assembled from the total reads, 12,024 of these unigenes were annotated in the NCBI NR protein database, and 9784 were annotated in the Swiss-Prot database. Of these annotated unigenes, 7903 and 5060 unigenes have been assigned to gene ontology categories and clusters of orthologous groups, respectively. Furthermore, 8647 unigenes were mapped to 257 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database. We found that 408 unigenes were related to insecticide resistance and metabolism. In addition, 23,122 simple sequence repeats (SSRs) were identified in 11,009 unigenes, and 100 PCR primers of SSR loci were used to validate the assembly quality and polymorphisms. These results provide a good platform for further investigations into the insecticide resistance of B. odoriphaga. Finally, the SSRs identified in B. odoriphaga may be a useful genomic resource.

De novo sequencing, assembly and characterization of antennal transcriptome of Anomala corpulenta Motschulsky (Coleoptera: Rutelidae)

Background

Anomala corpulenta is an important insect pest and can cause enormous economic losses in agriculture, horticulture and forestry. It is widely distributed in China, and both larvae and adults can cause serious damage. It is difficult to control this pest because the larvae live underground. Any new control strategy should exploit alternatives to heavily and frequently used chemical insecticides. However, little genetic research has been carried out on A. corpulenta due to the lack of genomic resources. Genomic resources could be produced by next generation sequencing technologies with low cost and in a short time. In this study, we performed de novo sequencing, assembly and characterization of the antennal transcriptome of A. corpulenta.

Results

Illumina sequencing technology was used to sequence the antennal transcriptome of A. corpulenta. Approximately 76.7 million total raw reads and about 68.9 million total clean reads were obtained, and then 35,656 unigenes were assembled. Of these unigenes, 21,463 of them could be annotated in the NCBI nr database, and, among the annotated unigenes, 11,154 and 6,625 unigenes could be assigned to GO and COG, respectively. Additionally, 16,350 unigenes could be annotated in the Swiss-Prot database, and 14,499 unigenes could map onto 258 pathways in the KEGG Pathway database. We also found 24 unigenes related to OBPs, 6 to CSPs, and in total 167 unigenes related to chemodetection. We analyzed 4 OBPs and 3CSPs sequences and their RT-qPCR results agreed well with their FPKM values.

Conclusion

We produced the first large-scale antennal transcriptome of A. corpulenta, which is a species that has little genomic information in public databases. The identified chemodetection unigenes can promote the molecular mechanistic study of behavior in A. corpulenta. These findings provide a general sequence resource for molecular genetics research on A. corpulenta.

Analysis of transcriptomes of three orb-web spider species reveals gene profiles involved in silk and toxin

As an ancient arthropod with a history of 390 million years, spiders evolved numerous morphological forms resulting from adaptation to different environments. The venom and silk of spiders, which have promising commercial applications in agriculture, medicine and engineering fields, are of special interests to researchers. However, little is known about their genomic components, which hinders not only understanding spider biology but also utilizing their valuable genes. Here we report on deep sequenced and de novo assembled transcriptomes of three orb-web spider species, Gasteracantha arcuata, Nasoonaria sinensis and Gasteracantha hasselti which are distributed in tropical forests of south China. With Illumina paired-end RNA-seq technology, 54 871, 101 855 and 75 455 unigenes for the three spider species were obtained, respectively, among which 9 300, 10 001 and 10 494 unique genes are annotated, respectively. From these annotated unigenes, we comprehensively analyzed silk and toxin gene components and structures for the three spider species. Our study provides valuable transcriptome data for three spider species which previously lacked any genetic/genomic data. The results have laid the first fundamental genomic basis for exploiting gene resources from these spiders.

De novo characterization of the Dialeurodes citri transcriptome: mining genes involved in stress resistance and simple sequence repeats (SSRs) discovery

The citrus whitefly, Dialeurodes citri (Ashmead), is one of the three economically important whitefly species that infest citrus plants around the world; however, limited genetic research has been focused on D. citri, partly because of lack of genomic resources. In this study, we performed de novo assembly of a transcriptome using Illumina paired-end sequencing technology (Illumina Inc., San Diego, CA, USA). In total, 36,766 unigenes with a mean length of 497 bp were identified. Of these unigenes, we identified 17,788 matched known proteins in the National Center for Biotechnology Information database, as determined by Blast search, with 5731, 4850 and 14,441 unigenes assigned to clusters of orthologous groups (COG), gene ontology (GO), and SwissProt, respectively. In total, 7507 unigenes were assigned to 308 known pathways. In-depth analysis of the data showed that 117 unigenes were identified as potentially involved in the detoxification of xenobiotics and 67 heat shock protein (Hsp) genes were associated with environmental stress. In addition, these enzymes were searched against the GO and COG database, and the results showed that the three major detoxification enzymes and Hsps were classified into 18 and 3, 6, and 8 annotations, respectively. In addition, 149 simple sequence repeats were detected. The results facilitate the investigation of molecular resistance mechanisms to insecticides and environmental stress, and contribute to molecular marker development. The findings greatly improve our genetic understanding of D. citri, and lay the foundation for future functional genomics studies on this species.

De novo assembly, gene annotation, and marker discovery in stored-product pest Liposcelis entomophila (Enderlein) using transcriptome sequences

Background

As a major stored-product pest insect, Liposcelis entomophila has developed high levels of resistance to various insecticides in grain storage systems. However, the molecular mechanisms underlying resistance and environmental stress have not been characterized. To date, there is a lack of genomic information for this species. Therefore, studies aimed at profiling the L. entomophila transcriptome would provide a better understanding of the biological functions at the molecular levels.

Methodology/Principal Findings

We applied Illumina sequencing technology to sequence the transcriptome of L. entomophila. A total of 54,406,328 clean reads were obtained and that de novo assembled into 54,220 unigenes, with an average length of 571 bp. Through a similarity search, 33,404 (61.61%) unigenes were matched to known proteins in the NCBI non-redundant (Nr) protein database. These unigenes were further functionally annotated with gene ontology (GO), cluster of orthologous groups of proteins (COG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. A large number of genes potentially involved in insecticide resistance were manually curated, including 68 putative cytochrome P450 genes, 37 putative glutathione S-transferase (GST) genes, 19 putative carboxyl/cholinesterase (CCE) genes, and other 126 transcripts to contain target site sequences or encoding detoxification genes representing eight types of resistance enzymes. Furthermore, to gain insight into the molecular basis of the L. entomophila toward thermal stresses, 25 heat shock protein (Hsp) genes were identified. In addition, 1,100 SSRs and 57,757 SNPs were detected and 231 pairs of SSR primes were designed for investigating the genetic diversity in future.

Conclusions/Significance

We developed a comprehensive transcriptomic database for L. entomophila. These sequences and putative molecular markers would further promote our understanding of the molecular mechanisms underlying insecticide resistance or environmental stress, and will facilitate studies on population genetics for psocids, as well as providing useful information for functional genomic research in the future.

Characterization of the heart transcriptome of the white shark (Carcharodon carcharias)

BACKGROUND

The white shark (Carcharodon carcharias) is a globally distributed, apex predator possessing physical, physiological, and behavioral traits that have garnered it significant public attention. In addition to interest in the genetic basis of its form and function, as a representative of the oldest extant jawed vertebrate lineage, white sharks are also of conservation concern due to their small population size and threat from overfishing. Despite this, surprisingly little is known about the biology of white sharks, and genomic resources are unavailable. To address this deficit, we combined Roche-454 and Illumina sequencing technologies to characterize the first transciptome of any tissue for this species.

RESULTS

From white shark heart cDNA we generated 665,399 Roche 454 reads (median length 387-bp) that were assembled into 141,626 contigs (mean length 503-bp). We also generated 78,566,588 Illumina reads, which we aligned to the 454 contigs producing 105,014 454/Illumina consensus sequences. To these, we added 3,432 non-singleton 454 contigs. By comparing these sequences to the UniProtKB/Swiss-Prot database we were able to annotate 21,019 translated open reading frames (ORFs) of ≥ 20 amino acids. Of these, 19,277 were additionally assigned Gene Ontology (GO) functional annotations. While acknowledging the limitations of our single tissue transcriptome, Fisher tests showed the white shark transcriptome to be significantly enriched for numerous metabolic GO terms compared to the zebra fish and human transcriptomes, with white shark showing more similarity to human than to zebra fish (i.e. fewer terms were significantly different). We also compared the transcriptome to other available elasmobranch sequences, for signatures of positive selection and identified several genes of putative adaptive significance on the white shark lineage. The white shark transcriptome also contained 8,404 microsatellites (dinucleotide, trinucleotide, or tetranucleotide motifs ≥ five perfect repeats). Detailed characterization of these microsatellites showed that ORFs with trinucleotide repeats, were significantly enriched for transcription regulatory roles and that trinucleotide frequency within ORFs was lower than for a wide range of taxonomic groups including other vertebrates.

CONCLUSION

The white shark heart transcriptome represents a valuable resource for future elasmobranch functional and comparative genomic studies, as well as for population and other biological studies vital for effective conservation of this globally vulnerable species.

SSR_pipeline: a bioinformatic infrastructure for identifying microsatellites from paired-end Illumina high-throughput DNA sequencing data

SSR_pipeline is a flexible set of programs designed to efficiently identify simple sequence repeats (e.g., microsatellites) from paired-end high-throughput Illumina DNA sequencing data. The program suite contains 3 analysis modules along with a fourth control module that can automate analyses of large volumes of data. The modules are used to 1) identify the subset of paired-end sequences that pass Illumina quality standards, 2) align paired-end reads into a single composite DNA sequence, and 3) identify sequences that possess microsatellites (both simple and compound) conforming to user-specified parameters. The microsatellite search algorithm is extremely efficient, and we have used it to identify repeats with motifs from 2 to 25 bp in length. Each of the 3 analysis modules can also be used independently to provide greater flexibility or to work with FASTQ or FASTA files generated from other sequencing platforms (Roche 454, Ion Torrent, etc.). We demonstrate use of the program with data from the brine fly Ephydra packardi (Diptera: Ephydridae) and provide empirical timing benchmarks to illustrate program performance on a common desktop computer environment. We further show that the Illumina platform is capable of identifying large numbers of microsatellites, even when using unenriched sample libraries and a very small percentage of the sequencing capacity from a single DNA sequencing run. All modules from SSR_pipeline are implemented in the Python programming language and can therefore be used from nearly any computer operating system (Linux, Macintosh, and Windows).

Microsatellites for the marsh fritillary butterfly: de novo transcriptome sequencing, and a comparison with amplified fragment length polymorphism (AFLP) markers

BACKGROUND

Until recently the isolation of microsatellite markers from Lepidoptera has proved troublesome, expensive and time-consuming. Following on from a previous study of Edith's checkerspot butterfly, Euphydryas editha, we developed novel microsatellite markers for the vulnerable marsh fritillary butterfly, E. aurinia. Our goal was to optimize the process in order to reduce both time and cost relative to prevailing techniques. This was accomplished by using a combination of previously developed techniques: in silico mining of a de novo assembled transcriptome sequence, and genotyping the microsatellites found there using an economic method of fluorescently labelling primers.

PRINCIPAL FINDINGS

In total, we screened nine polymorphic microsatellite markers, two of which were previously published, and seven that were isolated de novo. These markers were able to amplify across geographically isolated populations throughout Continental Europe and the UK. Significant deviations from Hardy-Weinberg equilibrium were evident in some populations, most likely due to the presence of null alleles. However, we used an F(st) outlier approach to show that these markers are likely selectively neutral. Furthermore, using a set of 128 individuals from 11 populations, we demonstrate consistency in population differentiation estimates with previously developed amplified fragment length polymorphism (AFLP) markers (r = 0.68, p<0.001).

SIGNIFICANCE

Rapid development of microsatellite markers for difficult taxa such as Lepidoptera, and concordant results with other putatively neutral molecular markers, demonstrate the potential of de novo transcriptional sequencing for future studies of population structure and gene flow that are desperately needed for declining species across fragmented landscapes.

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).

Rapid development of microsatellite markers for Pandanus boninensis (Pandanaceae) by pyrosequencing technology

PREMISE OF THE STUDY

To facilitate rapid development of microsatellite or simple sequence repeat (SSR) markers, an expressed sequence tags (EST) database was constructed for Pandanus boninensis, an evergreen tree endemic to the Bonin Islands, using pyrosequencing technology.

METHODS AND RESULTS

We designed primers for 340 EST-SSRs identified from 109620 pyrosequencing reads, 48 of which were tested for PCR amplification. Thirty-four primers provided clear amplification, and 26 of those 34 displayed clear polymorphic patterns in sampled populations, with mean expected heterozygosity at the amplified loci ranging from 0.022 to 0.742 (average 0.262).

CONCLUSIONS

The developed markers are promising tools for future genetic studies of P. boninensis and related species.

Isolation and characterization of EST-SSR markers in Schima mertensiana (Theaceae) using pyrosequencing technology

PREMISE OF THE STUDY

Microsatellite or simple sequence repeat (SSR) markers from expressed sequence tags (ESTs) in Schima mertensiana, an evergreen tree that is endemic to the Bonin Islands, were developed to facilitate evaluation of the species's population genetic structure.

METHODS AND RESULTS

Using 149831 ESTs derived by pyrosequencing S. mertensiana complementary DNA (cDNA), 683 EST-SSRs were developed, 48 of which were tested for amplification. Thirty of the 48 showed clear amplification and detected polymorphism in sampled populations, with mean expected heterozygosity at the amplified loci ranging from 0.045 to 0.874 (average 0.545).

CONCLUSIONS

The markers developed in this study can be used for future studies of the genetic structure of S. mertensiana and related species.

Searching microsatellites in DNA sequences: approaches used and tools developed

Microsatellite instability associated genomic activities and evolutionary changes have led to a renewed focus on microsatellite research. In last decade, a number of microsatellite mining tools have been introduced based on different computational approaches. The choice is generally made between slow but exhaustive dynamic programming based approaches, or fast and incomplete heuristic methods. Tools based on stochastic approaches are more popular due to their simplicity and added ornamental features. We have performed a comparative evaluation of the relative efficiency of some microsatellite search tools with their default settings. The graphical user interface, the statistical analysis of the output and ability to mine imperfect repeats are the most important criteria in selecting a tool for a particular investigation. However, none of the available tools alone provides complete and accurate information about microsatellites, and a lot depends on the discretion of the user.

Microsatellite markers derived from Calophyllum inophyllum (Clusiaceae) expressed sequence tags

PREMISE OF THE STUDY

Robust markers are required (inter alia) for assessing origins of Calophyllum inophyllum populations on the Bonin Islands, Japan. Therefore, informative expressed sequence tag (EST)-based microsatellite or simple sequence repeat (SSRs) markers in the species were sought.

METHODS AND RESULTS

Using 135378 ESTs derived from de novo pyrosequencing, primers for 475 EST-SSRs were developed, 48 of which were tested for PCR amplification. Thirty-six of the 48 primers showed clear amplification, with 23 displaying polymorphism in sampled populations. Expected heterozygosity in the samples from the Bonin Islands and Ryukyu Islands populations ranged from 0.041 to 0.697 and from 0.041 to 0.773, respectively.

CONCLUSIONS

As EST-SSRs are potentially tightly linked with functional genes, and reportedly more transferable to related species than anonymous genomic SSRs, the developed primers have utility for future studies of the origins, genetic structure, and conservation of C. inophyllum and related species.

Multiplexed microsatellite recovery using massively parallel sequencing

Conservation and management of natural populations requires accurate and inexpensive genotyping methods. Traditional microsatellite, or simple sequence repeat (SSR), marker analysis remains a popular genotyping method because of the comparatively low cost of marker development, ease of analysis and high power of genotype discrimination. With the availability of massively parallel sequencing (MPS), it is now possible to sequence microsatellite-enriched genomic libraries in multiplex pools. To test this approach, we prepared seven microsatellite-enriched, barcoded genomic libraries from diverse taxa (two conifer trees, five birds) and sequenced these on one lane of the Illumina Genome Analyzer using paired-end 80-bp reads. In this experiment, we screened 6.1 million sequences and identified 356,958 unique microreads that contained di- or trinucleotide microsatellites. Examination of four species shows that our conversion rate from raw sequences to polymorphic markers compares favourably to Sanger- and 454-based methods. The advantage of multiplexed MPS is that the staggering capacity of modern microread sequencing is spread across many libraries; this reduces sample preparation and sequencing costs to less than $400 (USD) per species. This price is sufficiently low that microsatellite libraries could be prepared and sequenced for all 1373 organisms listed as 'threatened' and 'endangered' in the United States for under $0.5 M (USD).