Vincent Dubut

High‐throughput microsatellite isolation through 454 GS‐FLX Titanium pyrosequencing of enriched DNA libraries

Microsatellites (or SSRs: simple sequence repeats) are among the most frequently used DNA markers in many areas of research. The use of microsatellite markers is limited by the difficulties involved in their de novo isolation from species for which no genomic resources are available. We describe here a high‐throughput method for isolating microsatellite markers based on coupling multiplex microsatellite enrichment and next‐generation sequencing on 454 GS‐FLX Titanium platforms. The procedure was calibrated on a model species (Apis mellifera) and validated on 13 other species from various taxonomic groups (animals, plants and fungi), including taxa for which severe difficulties were previously encountered using traditional methods. We obtained from 11 497 to 34 483 sequences depending on the species and the number of detected microsatellite loci ranged from 199 to 5791. We thus demonstrated that this procedure can be readily and successfully applied to a large variety of taxonomic groups, at much lower cost than would have been possible with traditional protocols. This method is expected to speed up the acquisition of high‐quality genetic markers for nonmodel organisms.

QDD version 3.1: a user-friendly computer program for microsatellite selection and primer design revisited: experimental validation of variables determining genotyping success rate.

Microsatellite marker development has been greatly simplified by the use of high‐throughput sequencing followed by in silico microsatellite detection and primer design. However, the selection of markers designed by the existing pipelines depends either on arbitrary criteria, or older studies on PCR success. Based on wet laboratory experiments, we have identified the following factors that are most likely to influence genotyping success rate: alignment score between the primers and the amplicon; the distance between primers and microsatellites; the length of the PCR product; target region complexity and the number of reads underlying the sequence. The QDD pipeline has been modified to include these most pertinent factors in the output to help the selection of markers. Furthermore, new features are also included in the present version: (i) not only raw sequencing reads are accepted as input, but also contigs, allowing the analysis of assembled high‐coverage data; (ii) input data can be both in fasta and fastq format to facilitate the use of Illumina and IonTorrent reads; (iii) A comparison to known transposable elements allows their detection; (iv) A contamination check can be carried out by BLASTing potential markers against the nucleotide (nt) database of NCBI; (v) QDD3 is now also available imbedded into a virtual machine making installation easier and operating system independent. It can be used both on command‐line version as well as integrated into a Galaxy server, providing a user‐friendly interface, as well as the possibility to utilize a large variety of NGS tools.

Representativeness of microsatellite distributions in genomes, as revealed by 454 GS-FLX Titanium pyrosequencing

BackgroundMicrosatellites are markers of choice in population genetics and genomics, as they provide useful insight into patterns and processes as diverse as genome evolutionary dynamics and demographic processes. The acquisition of microsatellites through multiplex-enriched libraries and 454 GS-FLX Titanium pyrosequencing is a promising new tool for the isolation of new markers in unknown genomes. This approach can also be used to evaluate the extent to which microsatellite-enriched libraries are representative of the genome from which they were isolated. In this study, we deciphered potential discrepancies in microsatellite content recovery for two reference genomes (Apis mellifera and Danio rerio), selected on the basis of their extreme heterogeneity in terms of the proportions and distributions of microsatellites on chromosomes.ResultsThe A. mellifera genome, in particular, was found to be highly heterogeneous, due to extremely high rates of recombination, with hotspots, but the only bias consistently introduced into pyrosequenced multiplex-enriched libraries concerned sequence length, with the overrepresentation of sequences 160 to 320 bp in length. Other deviations from expected proportions or distributions of motifs on chromosomes were observed, but the significance and intensity of these deviations was mostly limited. Furthermore, no consistent adverse competition between multiplexed probes was observed during the motif enrichment phase.ConclusionsThis approach therefore appears to be a promising strategy for improving the development of microsatellites, as it introduces no major bias in terms of the proportions and distribution of microsatellites.

Cross-species amplification of 41 microsatellites in European cyprinids: A tool for evolutionary, population genetics and hybridization studies

BackgroundCyprinids display the most abundant and widespread species among the European freshwater Teleostei and are known to hybridize quite commonly. Nevertheless, a limited number of markers for conducting comparative differentiation, evolutionary and hybridization dynamics studies are available to date.FindingsFive multiplex PCR sets were optimized in order to assay 41 cyprinid-specific polymorphic microsatellite loci (including 10 novel loci isolated from Chondrostoma nasus nasus, Chondrostoma toxostoma toxostoma and Leuciscus leuciscus) for 503 individuals (440 purebred specimens and 63 hybrids) from 15 European cyprinid species. The level of genetic diversity was assessed in Alburnus alburnus, Alburnoides bipunctatus, C. genei, C. n. nasus, C. soetta, C. t. toxostoma, L. idus, L. leuciscus, Pachychilon pictum, Rutilus rutilus, Squalius cephalus and Telestes souffia. The applicability of the markers was also tested on Abramis brama, Blicca bjoerkna and Scardinius erythrophtalmus specimens. Overall, between 24 and 37 of these markers revealed polymorphic for the investigated species and 23 markers amplified for all the 15 European cyprinid species.ConclusionsThe developed set of markers demonstrated its performance in discriminating European cyprinid species. Furthermore, it allowed detecting and characterizing hybrid individuals. These microsatellites will therefore be useful to perform comparative evolutionary and population genetics studies dealing with European cyprinids, what is of particular interest in conservation issues and constitutes a tool of choice to conduct hybridization studies.

Development of 55 novel polymorphic microsatellite loci for the critically endangered Zingel asper L. (Actinopterygii: Perciformes: Percidae) and cross-species amplification in five other percids

By combining biotin-enrichment protocol and next generation pyrosequencing, through 454 GS-FLX Titanium technology, 55 polymorphic microsatellites loci with perfect motif were isolated from the Rhone streber (Zingel asper), a critically endangered European fish species. Eight multiplex PCR kits were optimised in order to genotype a total of 58 polymorphic loci, including three previously published loci. The level of genetic diversity was assessed for 68 Z. asper, 30 Sander lucioperca, 33 Perca fluviatilis and four Gymnocephalus schraetzer individuals. Amplification success was also assessed on Romanichthys valsanicola and Zingel streber using single individuals. These markers will be useful to investigate the population structure of the highly fragmented Rhone streber. They represent a powerful tool for conservation issues and evolutionary approaches of this endemic species. Moreover, part of our markers demonstrated applicability to other percid species, allowing for potential applications to fisheries and aquaculture management.

Inter‐ and Extra‐Indian Admixture and Genetic Diversity in Reunion Island Revealed by Analysis of Mitochondrial DNA

Reunion Island is a French territory located in the western Indian Ocean. The genetic pattern of the Reunionese population has been shaped by contributions from highly contrasting regions of the world. Over the last 350 years, several migration waves and cultural and socio‐economic factors have led to the emergence of six main ethnic groups in Reunion. India is one of the principal regions that contributed to the setting up of the Reunionese population. Diversity, demographic and admixture analyses were performed on mtDNA variation of the Reunionese of Indian ancestry, including the Malbar and Zarab ethnic groups, in order to question their history. Using a phylogeographical approach, we generated and analysed quantitative data on the contribution of the Indian parental populations. Furthermore, we showed that the settlement of Reunion Island by Indians did not involve a founder effect, except in the very beginning of the Reunionese settlement (at the end of the 17th century). The accuracy of our results was optimised by a re‐evaluation of the classification of the Southern Asian mtDNA haplogroups. Finally, by comparing our results to a previous study dealing with the Reunionese population, we highlighted how ethno‐historical data are critical for reconstructing the complex history of multiethnic populations.

From Late Miocene to Holocene: Processes of Differentiation within the Telestes Genus (Actinopterygii: Cyprinidae)

Investigating processes and timing of differentiation of organisms is critical in the understanding of the evolutionary mechanisms involved in microevolution, speciation, and macroevolution that generated the extant biodiversity. From this perspective, the Telestes genus is of special interest: the Telestes species have a wide distribution range across Europe (from the Danubian district to Mediterranean districts) and have not been prone to translocation. Molecular data (mtDNA: 1,232 bp including the entire Cyt b gene; nuclear genome: 11 microsatellites) were gathered from 34 populations of the Telestes genus, almost encompassing the entire geographic range. Using several phylogenetic and molecular dating methods interpreted in conjunction with paleoclimatic and geomorphologic evidence, we investigated the processes and timing of differentiation of the Telestes lineages. The observed genetic structure and diversity were largely congruent between mtDNA and microsatellites. The Messinian Salinity Crisis (Late Miocene) seems to have played a major role in the speciation processes of the genus. Focusing on T. souffia, a species occurring in the Danube and Rhone drainages, we were able to point out several specific events from the Pleistocene to the Holocene that have likely driven the differentiation and the historical demography of this taxon. This study provides support for an evolutionary history of dispersal and vicariance with unprecedented resolution for any freshwater fish in this region.

Isolation and characterization of polymorphic microsatellite loci for the dace complex: Leuciscus leuciscus (Teleostei: Cyprinidae)

Ten novel polymorphic microsatellites were isolated from the dace complex (Leuciscus leuciscus), which is a European cyprinid species. Four multiplex polymerase chain reaction sets were optimized in order to genotype 26 polymorphic loci in all, including 16 previously published cyprinid‐specific loci. The level of genetic diversity was assessed in 142 dace individuals. We also successfully applied 26 of the microsatellites to 10 related species. These primers thus will be useful to assess population structure of the dace and other cyprinid species, with application for conservation issues and phylogeographical approaches.

Isolation and characterization of polymorphic microsatellite loci in the freshwater fishes Telestes souffia and Telestes muticellus (Teleostei: Cyprinidae)

Ten novel polymorphic microsatellites (seven with perfect motifs) were isolated from vairone species (Telestes souffia and Telestes muticellus), which are endangered European cyprinid species. Together with 11 previously published cyprinid‐specific loci, five multiplex sets were optimized, allowing the genotyping of 21 polymorphic loci. The level of genetic diversity was assessed in 97 individuals from the two species T. souffia and T. muticellus. We also successfully applied the 21 microsatellites to nine related species. These primers will thus be useful in assessing population structure of the vairone and other cyprinid species, with application for conservation issues and phylogeographical approaches.

Combining genetic and demographic data for prioritizing conservation actions: insights from a threatened fish species.

Abstract Prioritizing and making efficient conservation plans for threatened populations requires information at both evolutionary and ecological timescales. Nevertheless, few studies integrate multidisciplinary approaches, mainly because of the difficulty for conservationists to assess simultaneously the evolutionary and ecological status of populations. Here, we sought to demonstrate how combining genetic and demographic analyses allows prioritizing and initiating conservation plans. To do so, we combined snapshot microsatellite data and a 30-year-long demographic survey on a threatened freshwater fish species (Parachondrostoma toxostoma) at the river basin scale. Our results revealed low levels of genetic diversity and weak effective population sizes (<63 individuals) in all populations. We further detected severe bottlenecks dating back to the last centuries (200–800 years ago), which may explain the differentiation of certain populations. The demographic survey revealed a general decrease in the spatial distribution and abundance of P. toxostoma over the last three decades. We conclude that demo-genetic approaches are essential for (1) identifying populations for which both evolutionary and ecological extinction risks are high; and (2) proposing conservation plans targeted toward these at risk populations, and accounting for the evolutionary history of populations. We suggest that demo-genetic approaches should be the norm in conservation practices. We combined genetic and demographic data from a threatened freshwater fish species (Parachondrostoma toxostoma) at the river basin scale for conservation purposes. Genetic diversity and effective population sizes are very low, probably due to the strong genetic bottlenecks detected in this study. The species spatial distribution and abundance also decreased during the last decades.