Mining non-model genomic libraries for microsatellites: BAC versus EST libraries and the generation of allelic richness

Streamlining of Simple Sequence Repeat Data Mining Methodologies and Pipelines for Crop Scanning

Genetic markers are powerful tools for understanding genetic diversity and the molecular basis of traits, ushering in a new era of molecular breeding in crops. Over the past 50 years, DNA markers have rapidly changed, moving from hybridization-based and second-generation-based to sequence-based markers. Simple sequence repeats (SSRs) are the ideal markers in plant breeding, and they have numerous desirable properties, including their repeatability, codominance, multi-allelic nature, and locus specificity. They can be generated from any species, which requires prior sequence knowledge. SSRs may serve as evolutionary tuning knobs, allowing for rapid identification and adaptation to new circumstances. The evaluations published thus far have mostly ignored SSR polymorphism and gene evolution due to a lack of data regarding the precise placements of SSRs on chromosomes. However, NGS technologies have made it possible to produce high-throughput SSRs for any species using massive volumes of genomic sequence data that can be generated fast and at a minimal cost. Though SNP markers are gradually replacing the erstwhile DNA marker systems, SSRs remain the markers of choice in orphan crops due to the lack of genomic resources at the reference level and their adaptability to resource-limited labor. Several bioinformatic approaches and tools have evolved to handle genomic sequences to identify SSRs and generate primers for genotyping applications in plant breeding projects. This paper includes the currently available methodologies for producing SSR markers, genomic resource databases, and computational tools/pipelines for SSR data mining and primer generation. This review aims to provide a 'one-stop shop' of information to help each new user carefully select tools for identifying and utilizing SSRs in genetic research and breeding programs.

Microsatellite evolutionary rate and pattern in Schistocerca gregaria inferred from direct observation of germline mutations

Unravelling variation among taxonomic orders regarding the rate of evolution in microsatellites is crucial for evolutionary biology and population genetics research. The mean mutation rate of microsatellites tends to be lower in arthropods than in vertebrates, but data are scarce and mostly concern accumulation of mutations in model species. Based on parent-offspring segregations and a hierarchical Bayesian model, the mean rate of mutation in the orthopteran insect Schistocerca gregaria was estimated at 2.1e(-4) per generation per untranscribed dinucleotide locus. This is close to vertebrate estimates and one order of magnitude higher than estimates from species of other arthropod orders, such as Drosophila melanogaster and Daphnia pulex. We also found evidence of a directional bias towards expansions even for long alleles and exceptionally large ranges of allele sizes. Finally, at transcribed microsatellites, the mean rate of mutation was half the rate found at untranscribed loci and the mutational model deviated from that usually considered, with most mutations involving multistep changes that avoid disrupting the reading frame. Our direct estimates of mutation rate were discussed in the light of peculiar biological and genomic features of S. gregaria, including specificities in mismatch repair and the dependence of its activity to allele length. Shedding new light on the mutational dynamics of grasshopper microsatellites is of critical importance for a number of research fields. As an illustration, we showed how our findings improve microsatellite application in population genetics, by obtaining a more precise estimation of S. gregaria effective population size from a published data set based on the same microsatellites.

Polyandry and postcopulatory sexual selection in a wild population

When females mate multiply, postcopulatory sexual selection can occur via sperm competition and cryptic female choice. Although postcopulatory selection has the potential to be a major force in driving evolution, few studies have estimated its strength in natural populations. Likewise, although polyandry is widespread across taxa and is the focus of a growing body of research, estimates of natural female mating rates are still limited in number. Microsatellites can be used to estimate the number of mates represented in females' sperm stores and the number of sires contributing to their offspring, enabling comparisons both of polyandry and of two components of postcopulatory selection: the proportion of males that mate but fail to sire offspring, and the degree of paternity skew among the males that do sire offspring. Here, we estimate the number of mates and sires among wild females in the Hawaiian swordtail cricket Laupala cerasina. We compare these estimates to the actual mating rates and paternity shares we observed in a semi-natural population. Our results show that postcopulatory sexual selection operates strongly in this species: wild females mated with an average minimum of 3.6 males but used the sperm from only 58% of them. Furthermore, among the males that did sire offspring, paternity was significantly skewed. These patterns were similar to those observed in the field enclosure, where females mated with an average of 5.7 males and used the sperm from 62% of their mates, with paternity significantly skewed among the sires.

High opportunity for postcopulatory sexual selection under field conditions

In polygamous systems, male fitness is determined not only by mating success but also by fertilization success. Despite the growing interest over the past several decades in postcopulatory sexual selection, its relative importance compared to precopulatory sexual selection remains a subject of debate. Here, we use extensive behavioral observations of a seminatural population of Hawaiian swordtail crickets, Laupala cerasina, and molecular paternity assignment to measure the opportunities for pre- and postcopulatory selection. Because postcopulatory selection has the potential to operate at multiple stages, we also separately attribute its effects to factors specific to mating events versus factors specific to males. We find that variance in postcopulatory success is over four times as great as variance in precopulatory success, with most of it unexplained by male mating order or the number of nuptial gifts given. Surprisingly, we also find that male singing effort is under postcopulatory selection, suggesting that males who sing more frequently also have more competitive ejaculates. Our results are consistent with the hypothesis that high polyandry levels promote greater relative postcopulatory selection. They also highlight the need for detailed behavioral observations under conditions as natural as possible when measuring mating and reproductive success.

Fast and cost-effective mining of microsatellite markers using NGS technology: an example of a Korean water deer Hydropotes inermis argyropus

Background

Microsatellites, a special class of repetitive DNA sequence, have become one of the most popular genetic markers for population/conservation genetic studies. However, its application to endangered species has been impeded by high development costs, a lack of available sequences, and technical difficulties. The water deer Hydropotes inermis is the sole existing endangered species of the subfamily Capreolinae. Although population genetics studies are urgently required for conservation management, no species-specific microsatellite marker has been reported.

Methods

We adopted next-generation sequencing (NGS) to elucidate the microsatellite markers of Korean water deer and overcome these impediments on marker developments. We performed genotyping to determine the efficiency of this method as applied to population genetics.

Results

We obtained 98 Mbp of nucleotide information from 260,467 sequence reads. A total of 20,101 di-/tri-nucleotide repeat motifs were identified; di-repeats were 5.9-fold more common than tri-repeats. [CA]_n and [AAC]_n/[AAT]_n repeats were the most frequent di- and tri-repeats, respectively. Of the 17,206 di-repeats, 12,471 microsatellite primer pairs were derived. PCR amplification of 400 primer pairs yielded 106 amplicons and 79 polymorphic markers from 20 individual Korean water deer. Polymorphic rates of the 79 new microsatellites varied from 2 to 11 alleles per locus (H_e: 0.050–0.880; H_o: 0.000–1.000), while those of known microsatellite markers transferred from cattle to Chinese water deer ranged from 4 to 6 alleles per locus (H_e: 0.279–0.714; H_o: 0.300–0.400).

Conclusions

Polymorphic microsatellite markers from Korean water deer were successfully identified using NGS without any prior sequence information and deposited into the public database. Thus, the methods described herein represent a rapid and low-cost way to investigate the population genetics of endangered/non-model species.

Enhancing genome investigations in the mosquito Culex quinquefasciatus via BAC library construction and characterization

BACKGROUND

Culex quinquefasciatus (Say) is a major species in the Culex pipiens complex and an important vector for several human pathogens including West Nile virus and parasitic filarial nematodes causing lymphatic filariasis. It is common throughout tropical and subtropical regions and is among the most geographically widespread mosquito species. Although the complete genome sequence is now available, additional genomic tools are needed to improve the sequence assembly.

FINDINGS

We constructed a bacterial artificial chromosome (BAC) library using the pIndigoBAC536 vector and HindIII partially digested DNA isolated from Cx. quinquefasciatus pupae, Johannesburg strain (NDJ). Insert size was estimated by NotI digestion and pulsed-field gel electrophoresis of 82 randomly selected clones. To estimate genome coverage, each 384-well plate was pooled for screening with 29 simple sequence repeat (SSR) and five gene markers. The NDJ library consists of 55,296 clones arrayed in 144 384-well microplates. Fragment insert size ranged from 50 to 190 kb in length (mean = 106 kb). Based on a mean insert size of 106 kb and a genome size of 579 Mbp, the BAC library provides ~10.1-fold coverage of the Cx. quinquefasciatus genome. PCR screening of BAC DNA plate pools for SSR loci from the genetic linkage map and for four genes associated with reproductive diapause in Culex pipiens resulted in a mean of 9.0 positive plate pools per locus.

CONCLUSION

The NDJ library represents an excellent resource for genome assembly enhancement and characterization in Culex pipiens complex mosquitoes.

The 19 genomes of Drosophila: a BAC library resource for genus-wide and genome-scale comparative evolutionary research

The genus Drosophila has been the subject of intense comparative phylogenomics characterization to provide insights into genome evolution under diverse biological and ecological contexts and to functionally annotate the Drosophila melanogaster genome, a model system for animal and insect genetics. Recent sequencing of 11 additional Drosophila species from various divergence points of the genus is a first step in this direction. However, to fully reap the benefits of this resource, the Drosophila community is faced with two critical needs: i.e., the expansion of genomic resources from a much broader range of phylogenetic diversity and the development of additional resources to aid in finishing the existing draft genomes. To address these needs, we report the first synthesis of a comprehensive set of bacterial artificial chromosome (BAC) resources for 19 Drosophila species from all three subgenera. Ten libraries were derived from the exact source used to generate 10 of the 12 draft genomes, while the rest were generated from a strategically selected set of species on the basis of salient ecological and life history features and their phylogenetic positions. The majority of the new species have at least one sequenced reference genome for immediate comparative benefit. This 19-BAC library set was rigorously characterized and shown to have large insert sizes (125-168 kb), low nonrecombinant clone content (0.3-5.3%), and deep coverage (9.1-42.9×). Further, we demonstrated the utility of this BAC resource for generating physical maps of targeted loci, refining draft sequence assemblies and identifying potential genomic rearrangements across the phylogeny.