Thibaut Malausa

QDD: a user-friendly program to select microsatellite markers and design primers from large sequencing projects

SUMMARYnQDD is an open access program providing a user-friendly tool for microsatellite detection and primer design from large sets of DNA sequences. The program is designed to deal with all steps of treatment of raw sequences obtained from pyrosequencing of enriched DNA libraries, but it is also applicable to data obtained through other sequencing methods, using FASTA files as input. The following tasks are completed by QDD: tag sorting, adapter/vector removal, elimination of redundant sequences, detection of possible genomic multicopies (duplicated loci or transposable elements), stringent selection of target microsatellites and customizable primer design. It can treat up to one million sequences of a few hundred base pairs in the tag-sorting step, and up to 50,000 sequences in a single input file for the steps involving estimation of sequence similarity.nnnAVAILABILITYnQDD is freely available under the GPL licence for Windows and Linux from the following web site: http://www.univ-provence.fr/gsite/Local/egee/dir/meglecz/QDD.html.nnnSUPPLEMENTARY INFORMATIONnSupplementary data are available at Bioinformatics online.

Current trends in microsatellite genotyping.

Microsatellites have been popular molecular markers ever since their advent in the late eighties. Despite growing competition from new genotyping and sequencing techniques, the use of these versatile and cost‐effective markers continues to increase, boosted by successive technical advances. First, methods for multiplexing PCR have considerably improved over the last years, thereby decreasing genotyping costs and increasing throughput. Second, next‐generation sequencing technologies allow the identification of large numbers of microsatellite loci at reduced cost in non‐model species. As a consequence, more stringent selection of loci is possible, thereby further enhancing multiplex quality and efficiency. However, current practices are lagging behind. By surveying recently published population genetic studies relying on simple sequence repeats, we show that more than half of the studies lack appropriate quality controls and do not make use of multiplex PCR. To make the most of the latest technical developments, we outline the need for a well‐established strategy including standardized high‐throughput bench protocols and specific bioinformatic tools, from primer design to allele calling.

Accuracy and quality assessment of 454 GS-FLX Titanium pyrosequencing

BackgroundThe rapid evolution of 454 GS-FLX sequencing technology has not been accompanied by a reassessment of the quality and accuracy of the sequences obtained. Current strategies for decision-making and error-correction are based on an initial analysis by Huse et al. in 2007, for the older GS20 system based on experimental sequences. We analyze here the quality of 454 sequencing data and identify factors playing a role in sequencing error, through the use of an extensive dataset for Roche control DNA fragments.ResultsWe obtained a mean error rate for 454 sequences of 1.07%. More importantly, the error rate is not randomly distributed; it occasionally rose to more than 50% in certain positions, and its distribution was linked to several experimental variables. The main factors related to error are the presence of homopolymers, position in the sequence, size of the sequence and spatial localization in PT plates for insertion and deletion errors. These factors can be described by considering seven variables. No single variable can account for the error rate distribution, but most of the variation is explained by the combination of all seven variables.ConclusionsThe pattern identified here calls for the use of internal controls and error-correcting base callers, to correct for errors, when available (e.g. when sequencing amplicons). For shotgun libraries, the use of both sequencing primers and deep coverage, combined with the use of random sequencing primer sites should partly compensate for even high error rates, although it may prove more difficult than previous thought to distinguish between low-frequency alleles and errors.

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 11u2003497 to 34u2003483 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.

The biology of small, introduced populations, with special reference to biological control.

Populations are introduced into novel environments in different contexts, one being the biological control of pests. Despite intense efforts, less than half introduced biological control agents establish. Among the possible approaches to improve biological control, one is to better understand the processes that underpin introductions and contribute to ecological and evolutionary success. In this perspective, we first review the demographic and genetic processes at play in small populations, be they stochastic or deterministic. We discuss the theoretical outcomes of these different processes with respect to individual fitness, population growth rate, and establishment probability. Predicted outcomes differ subtly in some cases, but enough so that the evaluating results of introductions have the potential to reveal which processes play important roles in introduced populations. Second, we attempt to link the theory we have discussed with empirical data from biological control introductions. A main result is that there are few available data, but we nonetheless report on an increasing number of well‐designed, theory‐driven, experimental approaches. Combining demography and genetics from both theoretical and empirical perspectives highlights novel and exciting avenues for research on the biology of small, introduced populations, and great potential for improving both our understanding and practice of biological control.

DNA markers to disentangle complexes of cryptic taxa in mealybugs (Hemiptera: Pseudococcidae)

Mealybugs (Hemiptera: Pseudococcidae) are major pests of a wide range of crops and ornamental plants worldwide. Their high degree of morphological similarity makes them difficult to identify and limits their study and management. We aimed to identify a set of markers for the genetic characterization and identification of complexes of taxa in the Pseudococcidae. We surveyed and tested the genetic markers used in previous studies and then identified new markers for particularly relevant genomic regions for which no satisfactory markers were available. We tested all markers on a subset of four taxa distributed worldwide. Five markers were retained after this first screening: two regions of the mitochondrial cytochrome oxidase I gene, 28S‐D2, the entire internal transcriber space 2 locus and the rpS15‐16S region of the primary mealybug endosymbiont Tremblaya princeps. We then assessed the utility of these markers for the characterization and identification of 239 samples from 43 sites in France and Brazil. The five markers studied (i) successfully distinguished all species identified by morphological examination, (ii) disentangled complexes of species by revealing intraspecific genetic variation and identified a set of closely related taxa for which taxonomic status requires clarification through further studies, and (iii) facilitated the inference of phylogenetic relationships between the characterized taxa.

Genetic structure and gene flow in French populations of two Ostrinia taxa: host races or sibling species?

Most models of ecological speciation concern phytophagous insects in which speciation is thought to be driven by host shifts and subsequent adaptations of populations. Despite the ever‐increasing number of studies, the current evolutionary status of most models remains incompletely resolved, as estimates of gene flow between taxa remain extremely rare. We studied the population genetics of two taxa of the Ostrinia genus — one feeding mainly on maize and the other on mugwort and hop — occurring in sympatry throughout France. The actual level of divergence of these taxa was unknown because the genetic structure of populations had been investigated over a limited geographical area and the magnitude of gene flow between populations had not been estimated. We used 11 microsatellite markers to investigate the genetic structure of populations throughout France and the extent of gene flow between the two Ostrinia taxa at several sites at which they are sympatric. We observed clear genetic differentiation between most populations collected on the typical respective hosts of each taxon. However, populations displaying intermediate allelic frequencies were found on hop plants in southern France. Individual assignments revealed that this result could be accounted for by the presence of both taxa on the same host. Gene flow, estimated by determining the proportion of hybrids detected, was low: probably < 1% per generation, regardless of site. This indicates that the two Ostrinia taxa have reached a high level of genetic divergence and should be considered sibling species rather than host races.

Genetic structure of European and Mediterranean maize borer populations on several wild and cultivated host plants

Target pests may become resistant to Bacillus thuringiensis (Bt) toxins produced by trangenic maize (Zea mays L.). Untreated refuge areas are set aside to conserve high frequencies of susceptibility alleles: a delay in resistance evolution is expected if susceptible individuals from refuges mate randomly with resistant individuals from Bt fields. In principle, refuges can be toxin‐free maize or any other plant, provided it hosts sufficiently large pest populations mating randomly with populations from Bt‐maize fields. Our aim was to examine the suitability of several cultivated or weedy plants [pepper (Capsicum frutescens L.), sorghum (Sorghum spec.), sunflower (Helianthus annuus L.), cocklebur (Xanthium spec.), cantaloupe (Cucumis melo L.), and hop (Humulus lupulus L.)] as refuges for Ostrinia nubilalis (Hübner) (Lepidoptera: Crambidae) and Sesamia nonagrioides Lefebvre (Lepidoptera: Noctuidae), two major maize pests in southern Europe. Larvae of both species were collected on these plants. Their genetic population structure was examined at several allozyme loci. We found little or no evidence for an influence of geographic distance, but detected a significant host‐plant effect on the genetic differentiation for both species. Ostrinia nubilalis populations from sunflower, pepper, cocklebur, and sorghum appear to belong to the same genetic entity as populations collected on maize, but to differ from populations on hop. Accordingly, females from pepper and cocklebur produced exclusively the ‘Z’ type sexual pheromone, which, in France, characterizes populations developing on maize. Qualitatively, these plants (except hop) could thus serve as refuges for O. nubilalis; however, they may be of little use quantitatively as they were found much less infested than maize. Sesamia nonagrioides populations on maize and sorghum reached comparable densities, but a slight genetic differentiation was detected between both. The degree of assortative mating between populations feeding on both hosts must therefore be assessed before sorghum can be considered as a suitable refuge for this species.

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 June 2010 - 31 July 2010.

This article documents the addition of 205 microsatellite marker loci to the Molecular Ecology Resources Database. Loci were developed for the following species: Bagassa guianensis, Bulweria bulwerii, Camelus bactrianus, Chaenogobius annularis, Creontiades dilutus, Diachasmimorpha tryoni, Dioscorea alata, Euhrychiopsis lecontei, Gmelina arborea, Haliotis discus hannai, Hirtella physophora, Melanaphis sacchari, Munida isos, Thaumastocoris peregrinus and Tuberolachnus salignus. These loci were cross‐tested on the following species: Halobaena caerulea, Procellaria aequinoctialis, Oceanodroma monteiroi, Camelus ferus, Creontiades pacificus, Dioscorea rotundata, Dioscorea praehensilis, Dioscorea abyssinica, Dioscorea nummularia, Dioscorea transversa, Dioscorea esculenta, Dioscorea pentaphylla, Dioscorea trifida, Hirtella bicornis, Hirtella glandulosa, Licania alba, Licania canescens, Licania membranaceae, Couepia guianensis and 7 undescribed Thaumastocoris species.

Genome-wide survey and analysis of microsatellites in nematodes, with a focus on the plant-parasitic species Meloidogyne incognita

BackgroundMicrosatellites are the most popular source of molecular markers for studying population genetic variation in eukaryotes. However, few data are currently available about their genomic distribution and abundance across the phylum Nematoda. The recent completion of the genomes of several nematode species, including Meloidogyne incognita, a major agricultural pest worldwide, now opens the way for a comparative survey and analysis of microsatellites in these organisms.ResultsUsing MsatFinder, the total numbers of 1-6 bp perfect microsatellites detected in the complete genomes of five nematode species (Brugia malayi, Caenorhabditis elegans, M. hapla, M. incognita, Pristionchus pacificus) ranged from 2,842 to 61,547, and covered from 0.09 to 1.20% of the nematode genomes. Under our search criteria, the most common repeat motifs for each length class varied according to the different nematode species considered, with no obvious relation to the AT-richness of their genomes. Overall, (AT)n, (AG)nand (CT)nwere the three most frequent dinucleotide microsatellite motifs found in the five genomes considered. Except for two motifs in P. pacificus, all the most frequent trinucleotide motifs were AT-rich, with (AAT)nand (ATT)nbeing the only common to the five nematode species. A particular attention was paid to the microsatellite content of the plant-parasitic species M. incognita. In this species, a repertoire of 4,880 microsatellite loci was identified, from which 2,183 appeared suitable to design markers for population genetic studies. Interestingly, 1,094 microsatellites were identified in 801 predicted protein-coding regions, 99% of them being trinucleotides. When compared against the InterPro domain database, 497 of these CDS were successfully annotated, and further assigned to Gene Ontology terms.ConclusionsContrasted patterns of microsatellite abundance and diversity were characterized in five nematode genomes, even in the case of two closely related Meloidogyne species. 2,245 di- to hexanucleotide loci were identified in the genome of M. incognita, providing adequate material for the future development of a wide range of microsatellite markers in this major plant parasite.