Microsatellite flanking region similarities among different loci within insect species

Population genetic structure of Aedes aegypti subspecies in selected geographical locations in Sudan

Although knowledge of the composition and genetic diversity of disease vectors is important for their management, this is limiting in many instances. In this study, the population structure and phylogenetic relationship of the two Aedes aegypti subspecies namely Aedes aegypti aegypti (Aaa) and Aedes aegypti formosus (Aaf) in eight geographical areas in Sudan were analyzed using seven microsatellite markers. Hardy-Weinberg Equilibrium (HWE) for the two subspecies revealed that Aaa deviated from HWE among the seven microsatellite loci, while Aaf exhibited departure in five loci and no departure in two loci (A10 and M201). The Factorial Correspondence Analysis (FCA) plots revealed that the Aaa populations from Port Sudan, Tokar, and Kassala clustered together (which is consistent with the unrooted phylogenetic tree), Aaf from Fasher and Nyala populations clustered together, and Gezira, Kadugli, and Junaynah populations also clustered together. The Bayesian cluster analysis structured the populations into two groups suggesting two genetically distinct groups (subspecies). Isolation by distance test revealed a moderate to strong significant correlation between geographical distance and genetic variations (p = 0.003, r = 0.391). The migration network created using divMigrate demonstrated that migration and gene exchange between subspecies populations appear to occur based on their geographical proximity. The genetic structure of the Ae. aegypti subspecies population and the gene flow among them, which may be interpreted as the mosquito vector's capacity for dispersal, were revealed in this study. These findings will help in the improvement of dengue epidemiology research including information on the identity of the target vector/subspecies and the arboviruses vector surveillance program.

Novel Microsatellite Loci, Cross-Species Validation of Multiplex Assays, and By-Catch Mitochondrial Genomes on Ochthebius Beetles from Supratidal Rockpools

Here we focus on designing, for the first time, microsatellite markers for evolutionary and ecological research on aquatic beetles from the genus Ochthebius (Coleoptera, Hydraenidae). Some of these non-model species, with high cryptic diversity, exclusively inhabit supratidal rockpools, extreme and highly dynamic habitats with important anthropogenic threats. We analysed 15 individuals of four species (O. lejolisii, O. subinteger, O. celatus, and O. quadricollis) across 10 localities from the Mediterranean coasts of Spain and Malta. Using next-generation sequencing technology, two libraries were constructed to interpret the species of the two subgenera present consistently (Ochthebius s. str., O. quadricollis; and Cobalius, the rest of the species). Finally, 20 markers (10 for each subgenus) were obtained and successfully tested by cross-validation in the four species under study. As a by-catch, we could retrieve the complete mitochondrial genomes of O. lejolisii, O. quadricollis, and O. subinteger. Interestingly, the mitochondrial genome of O. quadricollis exhibited high genetic variability compared to already published data. The novel SSR panels and mitochondrial genomes for Ochthebius will be valuable in future research on species identification, diversity, genetic structure, and population connectivity in highly dynamic and threatened habitats such as supratidal coastal rockpools.

De Novo Mining and Validating Novel Microsatellite Markers to Assess Genetic Diversity in Maruca vitrata (F.), a Legume Pod Borer

Maruca vitrata (Fabricius) is an invasive insect pest capable of causing enormous economic losses to a broad spectrum of leguminous crops. Microsatellites are valuable molecular markers for population genetic studies; however, an inadequate number of M. vitrata microsatellite loci are available to carry out population association studies. Thus, we utilized this insect's public domain databases for mining expressed sequence tags (EST)-derived microsatellite markers. In total, 234 microsatellite markers were identified from 10053 unigenes. We discovered that trinucleotide repeats were the most predominant microsatellite motifs (61.53%), followed by dinucleotide repeats (23.50%) and tetranucleotide repeats (14.95%). Based on the analysis, twenty-five markers were selected for validation in M. vitrata populations collected from various regions of India. The number of alleles (Na), observed heterozygosity (Ho), and expected heterozygosity (He) ranged from 2 to 5; 0.00 to 0.80; and 0.10 to 0.69, respectively. The polymorphic loci showed polymorphism information content (PIC), ranging from 0.09 to 0.72. Based on the genetic distance matrix, the unrooted neighbor-joining dendrogram differentiated the selected populations into two discrete groups. The SSR markers developed and validated in this study will be helpful in population-level investigations of M. vitrata to understand the gene flow, demography, dispersal patterns, biotype differentiation, and host dynamics.

Large-scale genetic admixture suggests high dispersal in an insect pest, the apple fruit moth

Knowledge about population genetic structure and dispersal capabilities is important for the development of targeted management strategies for agricultural pest species. The apple fruit moth, Argyresthia conjugella (Lepidoptera, Yponomeutidae), is a pre-dispersal seed predator. Larvae feed on rowanberries (Sorbus aucuparia), and when rowanberry seed production is low (i.e., inter-masting), the moth switches from laying eggs in rowanberries to apples (Malus domestica), resulting in devastating losses in apple crops. Using genetic methods, we investigated if this small moth expresses any local genetic structure, or alternatively if gene flow may be high within the Scandinavian Peninsula (~850.000 km², 55^o - 69^o N). Genetic diversity was found to be high (n = 669, mean He = 0.71). For three out of ten tetranucleotide STRs, we detected heterozygote deficiency caused by null alleles, but tests showed little impact on the overall results. Genetic differentiation between the 28 sampling locations was very low (average FST = 0.016, P < 0.000). Surprisingly, we found that all individuals could be assigned to one of two non-geographic genetic clusters, and that a third, geographic cluster was found to be associated with 30% of the sampling locations, with weak but significant signals of isolation-by-distance. Conclusively, our findings suggest wind-aided dispersal and spatial synchrony of both sexes of the apple fruit moth over large areas and across very different climatic zones. We speculate that the species may recently have had two separate genetic origins caused by a genetic bottleneck after inter-masting, followed by rapid dispersal and homogenization of the gene pool across the landscape. We suggest further investigations of spatial genetic similarities and differences of the apple fruit moth at larger geographical scales, through life-stages, across inter-masting, and during attacks by the parasitoid wasp (Microgaster politus).

Role of microsatellites in genetic analysis of Bombyx mori silkworm: a review

In the genome of Bombyx mori Linnaeus (1758), the microsatellites, or simple sequence repeats (SSR), feature among their particular characteristics a high adenine and thymine (A/T) content, low number of repeats, low frequency, and a grouping in "families" with similar flanking regions. Such characteristics may be the result of a complex interaction between factors that limit the size and dispersion of SSR loci—such as their high association with transposons—and mean that microsatellites within this taxon suitable as molecular markers are relatively rare. The determination of genetic profiles in populations and cell lines has not been affected owing to the high level of polymorphism, nor has the analysis of diversity, structure and genetic relationships. However, the scarcity of suitable microsatellites has restricted their application in genetic mapping, limiting them to preliminary identification of gene location of genes or quantitative trait loci (QTLs) related to thermotolerance, resistance to viruses, pigmentation patterns, body development and the weight of the cocoon, the cortex, the pupa and the filament. The review confirms that, as markers, microsatellites are versatile and perform well. They could thus be useful both to advance research in emerging countries with few resources seeking to promote sericulture in their territories, and to advance in the genetic and molecular knowledge of characteristics of productive and biological interest, given the latest technological developments in terms of the sequencing, identification, isolation and genotyping of SSR loci.

Development of novel microsatellites for population genetic analysis of Phenacoccus solenopsis Tinsley (Hemipeta: Pseudoccoccidae) based on genomic analysis

The cotton mealybug, Phenacoccus solenopsis Tinsley (Hemipeta: Pseudoccoccidae), is an aggressively invasive pest causing huge economic losses of crops around the world. In this study, we developed genome-wide microsatellites for population genetic analysis of P. solenopsis. We obtained a random genome of P. solenopsis with a size of 267.07 Mb and scaffold N50 of 14.12 Kb. In total 115,639 microsatellites were isolated from the genome, of which those with trinucleotide motifs were the most abundant. Forty-two polymorphic loci were selected for primer validation based on three populations. Allele numbers varied from 2 to 5 with an average value of 2.5 per locus, and allelic richness ranged from 1.00 to 4.48. The observed heterozygosity (H₀) and expected heterozygosity (H_E) ranged from 0.00 to 0.92 and 0.00 to 0.73, respectively. Population genetic structure analysis based on the developed markers revealed strong differentiation between three populations of P. solenopsis collected from its invasive range in China. The microsatellites developed in our study should provide efficient genetic markers for population level studies of P. solenopsis to reveal invasion history and patterns of dispersal.

Using NextRAD sequencing to infer movement of herbivores among host plants

Herbivores often move among spatially interspersed host plants, tracking high-quality resources through space and time. This dispersal is of particular interest for vectors of plant pathogens. Existing molecular tools to track such movement have yielded important insights, but often provide insufficient genetic resolution to infer spread at finer spatiotemporal scales. Here, we explore the use of Nextera-tagmented reductively-amplified DNA (NextRAD) sequencing to infer movement of a highly-mobile winged insect, the potato psyllid (Bactericera cockerelli), among host plants. The psyllid vectors the pathogen that causes zebra chip disease in potato (Solanum tuberosum), but understanding and managing the spread of this pathogen is limited by uncertainty about the insect's host plant(s) outside of the growing season. We identified 1,978 polymorphic loci among psyllids separated spatiotemporally on potato or in patches of bittersweet nightshade (S. dulcumara), a weedy plant proposed to be the source of potato-colonizing psyllids. A subset of the psyllids on potato exhibited genetic similarity to insects on nightshade, consistent with regular movement between these two host plants. However, a second subset of potato-collected psyllids was genetically distinct from those collected on bittersweet nightshade; this suggests that a currently unrecognized source, i.e., other nightshade patches or a third host-plant species, could be contributing to psyllid populations in potato. Oftentimes, dispersal of vectors of pathogens must be tracked at a fine scale in order to understand, predict, and manage disease spread. We demonstrate that emerging sequencing technologies that detect genome-wide SNPs of a vector can be used to infer such localized movement.

Genetic differentiation of the regional Plutella xylostella populations across the Taiwan Strait based on identification of microsatellite markers

Movement of individuals through events, such as storms or crop transportation, may affect survival and distribution of insect pests, as well as population genetic structure at a regional scale. Understanding what factors contribute to gene flow in pest populations remains very important for sustainable pest management. The diamondback moth (Plutella xylostella) is an insect pest well known for its capacity of moving over short to long distances. Here, we used newly isolated microsatellite markers to analyze the genetic structure of nine populations across the Taiwan Strait of China (Taiwan and Fujian). A total of 12,152 simple sequence repeats (SSRs) were initially identified from the P. xylostella transcriptome (~94 Mb), with an average of 129 SSRs per Mb. Nine SSRs were validated to be polymorphic markers, and eight were used for this population genetic study. Our results showed that the P. xylostella populations could be divided into distinct two clusters, which is likely due to the year-round airflows in this region. A pattern of isolation by distance among the local populations within Fujian was found, and may be related to vegetable transportation. Considering the complexity of the P. xylostella population genetic structure from local and regional to global levels, we propose that developing ecologically sound strategies for managing this pest will require knowledge of the link between behavioral and population ecology and its genetic structure.

Genetic analysis of Indian tasar silkmoth (Antheraea mylitta) populations

Indian tasar silkmoth, Antheraea mylitta is an economically important wild silkmoth species distributed across India. A number of morphologically and ethologically well-defined ecotypes are known for this species that differ in their primary food plant specificity. Most of these ecotypes do not interbreed in nature, but are able to produce offspring under captive conditions. Microsatellite markers were developed for A. mylitta, and out of these, ten well-behaved microsatellite loci were used to analyze the population structure of different ecoraces. A total of 154 individual moths belonging to eight different ecoraces, were screened at each locus. Hierarchical analysis of population structure using Analysis of MOlecular VAriance (AMOVA) revealed significant structuring (F_ST = 0.154) and considerable inbreeding (F_IS = 0.505). A significant isolation by distance was also observed. The number of possible population clusters was investigated using distance method, Bayesian algorithm and self organization maps (SOM). The first two methods revealed two distinct clusters, whereas the SOM showed the different ecoraces not to be clearly differentiated. These results suggest that although there is a large degree of phenotypic variation among the different ecoraces of A. mylitta, genetically they are not very different, and the phenotypic differences may largely be a result of their respective ecology.

A novel set of microsatellite markers for the European Grapevine Moth Lobesia botrana isolated using next-generation sequencing and their utility for genetic characterization of populations from Europe and the Middle East

Using a high-throughput 454 pyrosequencing approach a novel set of microsatellite markers was developed for one of the key grapevine insect pests, the European grapevine moth Lobesia botrana (Lepidoptera: Tortricidae). 20 primer pairs flanking a microsatellite motif were designed based on the sequences obtained and were subsequently evaluated in a sample of 14 L. botrana populations from Europe and the Middle East. 11 markers showed stable and reproducible amplification patterns; however, one of the 11 markers was monomorphic in all L. botrana populations analysed. Estimated frequencies of null alleles of more than 20% were evident for two of the markers tested, but varied substantially depending on the respective L. botrana population. In 12 of the 14 L. botrana populations observed heterozygosities were lower to those expected under Hardy-Weinberg equilibrium, indicating a deficiency of heterozygotes in the respective populations. The overall F ST value of 0.075 suggested a moderate but significant genetic differentiation between the L. botrana populations included in this study. In addition, a clear geographic structure was detected in the set of samples, evident through a significant isolation by distance and through results from structure analysis. In structure analysis, L. botrana populations were grouped in two clearly separated clusters according to their European (Spain, Italy, Germany) or Middle Eastern (Israel, Syria, Turkey) origin. This novel set of microsatellite markers can now be applied to study the evolutionary ecology of this species including host shifts and host adaptation as well as spread of individuals across worldwide viticulture.

Landscape structure and the genetic effects of a population collapse

Both landscape structure and population size fluctuations influence population genetics. While independent effects of these factors on genetic patterns and processes are well studied, a key challenge is to understand their interaction, as populations are simultaneously exposed to habitat fragmentation and climatic changes that increase variability in population size. In a population network of an alpine butterfly, abundance declined 60-100% in 2003 because of low over-winter survival. Across the network, mean microsatellite genetic diversity did not change. However, patch connectivity and local severity of the collapse interacted to determine allelic richness change within populations, indicating that patch connectivity can mediate genetic response to a demographic collapse. The collapse strongly affected spatial genetic structure, leading to a breakdown of isolation-by-distance and loss of landscape genetic pattern. Our study reveals important interactions between landscape structure and temporal demographic variability on the genetic diversity and genetic differentiation of populations. Projected future changes to both landscape and climate may lead to loss of genetic variability from the studied populations, and selection acting on adaptive variation will likely occur within the context of an increasing influence of genetic drift.

Characterization of 12 Novel Microsatellite Markers of Sogatella furcifera (Hemiptera: Delphacidae) Identified From Next-Generation Sequence Data

The white-backed planthopper, Sogatella furcifera (Horváth) (Hemiptera: Delphacidae), is a major pest of rice and has long-range migratory behavior in Asia. Microsatellite markers (simple sequence repeats) have been widely used to determine the origins and genetic diversity of insect pests. We identified novel microsatellite loci for S. furcifera samples collected from Laos, Vietnam, and three localities in Bangladesh from next-generation Roche 454 pyrosequencing data. Size polymorphism at 12 microsatellite loci was verified for 40 adult individuals collected from Shinan, South Korea. The average number of alleles per locus was 7.92. The mean values of observed (H(o)) and expected heterozygosities (H(E)) were 0.615 and 0.757, respectively. These new microsatellite markers will be a resource for future ecological genetic studies of S. furcifera samples across more broad geographic regions in Asia and may assist in estimations of genetic differentiation and gene flow among populations for implementation of more effective management strategies to control this serious rice pest.

Characterization of novel microsatellite markers for Hyphantria cunea and implications for other Lepidoptera

This is the first report of microsatellite markers (simple sequence repeats, SSR) for fall webworm, Hyphantria cunea (Drury) (Lepidoptera: Arctiidae), an important quarantine pest in some European and Asian countries. Here, we developed 48 microsatellite markers for H. cunea from SSR enrichment libraries. Sequences isolated from libraries were sorted into four categories and analyzed. Our results suggest that sequences classified as Grouped should not be used for microsatellite primer design. The genetic diversity of microsatellite loci was assessed in 72 individuals from three populations. The number of alleles per locus ranged from 2 to 5 with an average of 3. The observed and expected heterozygosities of loci ranged from 0 to 0.958 and 0 to 0.773, respectively. A total of 18 out of 153 locus/population combinations deviated significantly from Hardy-Weinberg equilibrium. Moreover, significant linkage disequilibrium was detected in one pair of loci (1275 pairs in total). In the neutral test, two loci were grouped into the candidate category for positive selection and the remainder into the neutral category. In addition, a complex mutation pattern was observed for these loci, and F ST performed better than did R ST for the estimation of population differentiation in different mutation patterns. The results of the present study can be used for population genetic studies of H. cunea.

MuTAnT: a family of Mutator-like transposable elements targeting TA microsatellites in Medicago truncatula

Transposable elements (TEs) are mobile DNA segments, abundant and dynamic in plant genomes. Because their mobility can be potentially deleterious to the host, a variety of mechanisms evolved limiting that negative impact, one of them being preference for a specific target insertion site. Here, we describe a family of Mutator-like DNA transposons in Medicago truncatula targeting TA microsatellites. We identified 218 copies of MuTAnTs and an element carrying a complete ORF encoding a mudrA-like transposase. Most insertion sites are flanked by a variable number of TA tandem repeats, indicating that MuTAnTs are specifically targeting TA microsatellites. Other TE families flanked by TA repeats (e.g. TAFT elements in maize) were described previously, however we identified the first putative autonomous element sharing that characteristics with a related group of short non-autonomous transposons.

Development of Microsatellite Markers in the Deep-Sea Cup Coral Desmophyllum dianthus by 454 Sequencing and Cross-Species Amplifications in Scleractinia Order

Microsatellite loci were isolated for the first time for the deep-sea coral Desmophyllum dianthus, using 454 GS-FLX Titanium pyrosequencing. We developed conditions for amplifying 24 markers in 10 multiplex reactions. Three to 16 alleles per locus were detected across 25 samples analyzed from Santa Maria di Leuca coral province (Mediterranean Sea). For the 24 polymorphic loci, observed and expected heterozygosities ranged from 0.211 to 0.880 and 0.383 to 0.910, respectively; 3 loci deviated from Hardy-Weinberg equilibrium, after null allele and sequential Holm-Bonferroni corrections. These newly isolated microsatellites are very useful genetic markers that provide data for future conservation strategies. Cross-amplification of these microsatellites, tested in 46 coral species, representing 40 genera, and 10 families of the phylum Cnidaria, produced informative allelic profiles for 1 to 24 loci. The utility of extending analyses to cross-species amplifications is also discussed.

Developing Exon-Primed Intron-Crossing (EPIC) markers for population genetic studies in three Aedes disease vectors

Aedes aegypti, Aedes notoscriptus, and Aedes albopictus are important vectors of many arboviruses implicated in human disease such as dengue fever. Genetic markers applied across vector species can provide important information on population structure, gene flow, insecticide resistance, and taxonomy, however, robust microsatellite markers have proven difficult to develop in these species and mosquitoes generally. Here we consider the utility and transferability of 15 Ribosome protein (Rp) Exon-Primed Intron-Crossing (EPIC) markers for population genetic studies in these 3 Aedes species. Rp EPIC markers designed for Ae. aegypti also successfully amplified populations of the sister species, Ae. albopictus, as well as the distantly related species, Ae. notoscriptus. High SNP and good indel diversity in sequenced alleles plus support for amplification of the same regions across populations and species were additional benefits of these markers. These findings point to the general value of EPIC markers in mosquito population studies.

The genetic structure of Asian corn borer, Ostrinia furnacalis, populations in China: haplotype variance in northern populations and potential impact on management of resistance to transgenic maize

Asian corn borer, Ostrinia furnacalis (Guenée), is a severe pest that infests cultivated maize in the major production regions of China. Populations show genotype-by-environment variation in voltinism, such that populations with a single generation (univoltine) are fixed in Northern China where growing seasons are short. Low genetic differentiation was found among samples from 33 collection sites across China and one site from North Korea (n=1673) using variation at 6 nuclear microsatellite loci (ENA corrected global FST=0.020; P value<0.05). Analysis of molecular variance indicated that geographic region, number of generations or voltinism accounted for <0.38% of the total genetic variation at nuclear loci and was corroborated by clustering of co-ancestries among genotypes using the program STRUCTURE. In contrast, a mitochondrial haplotype network identified 4 distinct clusters, where 70.5% of samples from univoltine populations were within a single group. Univoltine populations were also placed into a unique cluster using Population Graph and Principal component analyses, which showed significant differentiation with multivoltine populations (φST=0.400; P value<0.01). This study suggests that gene flow among O. furnacalis in China may be high among regions, with the exception of northeastern localities. Haplotype variation may be due to random genetic drift resulting from partial reproductive isolation between univoltine and multivoltine O. furnacalis populations. Such reproductive isolation might impact the potential spread of alleles that confer resistance to transgenic maize in China.

Characterization of nineteen microsatellite markers and development of multiplex PCRs for the wedge clam Donax trunculus (Mollusca: Bivalvia)

The wedge clam Donax trunculus is an Atlantic-Mediterranean warm-temperate species found from Senegal to the northern coast of France, including the Mediterranean and Black Sea. It is commercially exploited in several European countries and constitutes an important fishing resource due to its high economical value. To contribute to its conservation and management, nineteen microsatellite markers were isolated from two enriched genomic libraries. These loci were characterized in 30 clams from a single population from northwest Spain. The number of alleles per locus ranged from 2 to 17 and observed and expected heterozygosity varied from zero to 0.714 and from 0.078 to 0.950, respectively. Linkage disequilibrium was not detected and nine loci were in agreement with Hardy-Weinberg equilibrium. Fifteen polymorphic markers were arranged into three multiplex PCR sets to reduce both time and cost of microsatellite genotyping. This is the first time that polymorphic microsatellite markers have been reported for D. trunculus. These new markers provide a valuable resource for future population genetics studies and management and culture of this species.

Characterization of microsatellite DNA libraries from three mealybug species and development of microsatellite markers for Pseudococcus viburni (Hemiptera: Pseudococcidae)

Mealybugs (Hemiptera: Pseudococcidae) are important pests for crops worldwide. Different species, cryptic taxa under the same species name or even populations within a species can differ in biological characteristics, such as phenology, resistance to insecticides, virus transmission and susceptibility to natural enemies. Therefore, their management efficacy depends on their accurate identification. Microsatellite genetic markers are efficient in revealing the fine-scale taxonomic status of insects, both at inter- and intra-specific level. Despite their potential uses, microsatellites have been developed only for one mealybug species so far. Hence, it is unclear whether microsatellites may be useful to assess mealybug population differentiation and structuring. In this work, we tested the feasibility of developing microsatellite markers in mealybugs by: (i) producing and characterizing microsatellite DNA libraries for three species: Pseudococcus viburni, Pseudococcus comstocki and Heliococcus bohemicus, and (ii) by developing and testing markers for Ps. viburni. The obtained libraries contained balanced percentages of dinucleotide (ranging from 15 to 25%) and trinucleotide (from 5 to 17%) motifs. The marker setup for Ps. viburni was successful, although 70% of the primers initially tested were discarded for a lack of polymorphism. Finally, 25 markers were combined in two multiplex polymerase chain reactions with 21 displaying no evidence of deviation from Hardy-Weinberg equilibrium. Ps. viburni markers were tested on one population from France and one from Chile. The markers revealed a significant genetic differentiation between the two populations with an Fst estimate of 0.266.

Random DNA libraries from three species of the stick insect genus Bacillus (Insecta: Phasmida): repetitive DNA characterization and first observation of polyneopteran MITEs

The repetitive DNA content of the stick insect species Bacillus rossius (facultative parthenogenetic), Bacillus grandii (gonochoric), and Bacillus atticus (obligate parthenogenetic) was analyzed through the survey of random genomic libraries roughly corresponding to 0.006% of the genome. By repeat masking, 19 families of transposable elements were identified (two LTR and six non-LTR retrotransposons; 11 DNA transposons). Moreover, a de novo analysis revealed, among the three libraries, the first MITE family observed in polyneopteran genomes. On the whole, transposable element abundance represented 23.3% of the genome in B. rossius, 22.9% in B. atticus, and 18% in B. grandii. Tandem repeat content in the three libraries is much lower: 1.32%, 0.64%, and 1.86% in B. rossius, B. grandii, and B. atticus, respectively. Microsatellites are the most abundant in all species. Minisatellites were only found in B. rossius and B. atticus, and five monomers belonging to the Bag320 satellite family were detected in B. atticus. Assuming the survey provides adequate representation of the relative genome, the obligate parthenogenetic species (B. atticus), compared with the other two species analyzed, does not show a lower transposable element content, as expected from some theoretical and empirical studies.

Isolation and characterization of ten polymorphic microsatellite loci in Ixodes arboricola, and cross-amplification in three other Ixodes species

We characterized ten polymorphic microsatellite loci from the tree-hole tick, Ixodes arboricola. Loci were screened in 11-18 individuals from three Belgian populations and five to ten alleles were found at each locus. Seven loci did not show deviations from Hardy-Weinberg equilibrium conditions and there were no indications for null alleles at these loci. The three other loci showed significant heterozygote deficiencies in at least one population, and a high potential for the occurrence of null alleles. We observed no effect of potential host DNA on the scoring of the microsatellites. Cross-amplification of the microsatellites was tested in eight specimens of three congeneric species: I. ricinus, I. hexagonus and I. frontalis. Depending on the species, six or seven of the loci were amplified in ≥ 4 of the 8 specimens and were polymorphic in each of these species (except for Ixaf 11 in I. frontalis and I. ricinus). These loci thus provide a tool for population genetic analysis of I. arboricola. The suitability of these markers needs to be further investigated in its congeners.

Cross-amplification of microsatellites from the codling moth Cydia pomonella to three other species of the tribe Grapholitini (Lepidoptera: Tortricidae)

This study examined cross-species amplification of 33 microsatellite markers, previously developed for Cydia pomonella, in three related fruit moth species of the same tribe (Grapholitini), namely Grapholita molesta, Grapholita funebrana and Grapholita lobarzewskii. Eight microsatellite loci yielded polymorphic products for G. molesta, nine for G. funebrana and 11 for G. lobarzewskii. At all these loci, the number of alleles ranged between four and 11 in G. molesta, and between four and nine in G. funebrana and G. lobarzewskii each. The successful cross-amplified loci can be used for research on population genetics and gene flow of the three target species.

Patterns of evolutionary conservation of microsatellites (SSRs) suggest a faster rate of genome evolution in Hymenoptera than in Diptera

Microsatellites, or simple sequence repeats (SSRs), are common and widespread DNA elements in genomes of many organisms. However, their dynamics in genome evolution is unclear, whereby they are thought to evolve neutrally. More available genome sequences along with dated phylogenies allowed for studying the evolution of these repetitive DNA elements along evolutionary time scales. This could be used to compare rates of genome evolution. We show that SSRs in insects can be retained for several hundred million years. Different types of microsatellites seem to be retained longer than others. By comparing Dipteran with Hymenopteran species, we found very similar patterns of SSR loss during their evolution, but both taxa differ profoundly in the rate. Relative to divergence time, Diptera lost SSRs twice as fast as Hymenoptera. The loss of SSRs on the Drosophila melanogaster X-chromosome was higher than on the other chromosomes. However, accounting for generation time, the Diptera show an 8.5-fold slower rate of SSR loss than the Hymenoptera, which, in contrast to previous studies, suggests a faster genome evolution in the latter. This shows that generation time differences can have a profound effect. A faster genome evolution in these insects could be facilitated by several factors very different to Diptera, which is discussed in light of our results on the haplodiploid D. melanogaster X-chromosome. Furthermore, large numbers of SSRs can be found to be in synteny and thus could be exploited as a tool to investigate genome structure and evolution.

Breakdown of phylogenetic signal: a survey of microsatellite densities in 454 shotgun sequences from 154 non model eukaryote species

Microsatellites are ubiquitous in Eukaryotic genomes. A more complete understanding of their origin and spread can be gained from a comparison of their distribution within a phylogenetic context. Although information for model species is accumulating rapidly, it is insufficient due to a lack of species depth, thus intragroup variation is necessarily ignored. As such, apparent differences between groups may be overinflated and generalizations cannot be inferred until an analysis of the variation that exists within groups has been conducted. In this study, we examined microsatellite coverage and motif patterns from 454 shotgun sequences of 154 Eukaryote species from eight distantly related phyla (Cnidaria, Arthropoda, Onychophora, Bryozoa, Mollusca, Echinodermata, Chordata and Streptophyta) to test if a consistent phylogenetic pattern emerges from the microsatellite composition of these species. It is clear from our results that data from model species provide incomplete information regarding the existing microsatellite variability within the Eukaryotes. A very strong heterogeneity of microsatellite composition was found within most phyla, classes and even orders. Autocorrelation analyses indicated that while microsatellite contents of species within clades more recent than 200 Mya tend to be similar, the autocorrelation breaks down and becomes negative or non-significant with increasing divergence time. Therefore, the age of the taxon seems to be a primary factor in degrading the phylogenetic pattern present among related groups. The most recent classes or orders of Chordates still retain the pattern of their common ancestor. However, within older groups, such as classes of Arthropods, the phylogenetic pattern has been scrambled by the long independent evolution of the lineages.

Long microsatellites and unusually high levels of genetic diversity in the Orthoptera

Much remains to be learned about the mutational processes governing the evolution of microsatellite repeat regions and the associated levels of genetic diversity observed at microsatellite markers across populations or species. An extensive survey of microsatellite variation in 210 insect species from six major orders revealed that within Orthopterans, which are characterized by giant genomes, levels of genetic diversity were ~20% higher and microsatellite repeat arrays were longer than in any other group. Because of the mutation dependence on repeat length, this result suggests a higher microsatellite loci mutation rate in the Orthoptera. We deem it plausible that differences among insect orders, either in mismatch repair systems or in abundance of transposable element-derived microsatellites, can shape the size distribution of both genomes and microsatellite repeat regions. Our findings emphasise that observed levels of genetic diversity can greatly vary across species (orders at least) because of molecular differences in the mechanisms that determine microsatellite size, and are therefore critical to conservation and population genetics studies, where microsatellite repeat variability is primarily interpreted in terms of population demography and history.

Breakdown of phylogenetic signal: a survey of microsatellite densities in 454 shotgun sequences from 154 non model eukaryote species

Microsatellites are ubiquitous in Eukaryotic genomes. A more complete understanding of their origin and spread can be gained from a comparison of their distribution within a phylogenetic context. Although information for model species is accumulating rapidly, it is insufficient due to a lack of species depth, thus intragroup variation is necessarily ignored. As such, apparent differences between groups may be overinflated and generalizations cannot be inferred until an analysis of the variation that exists within groups has been conducted. In this study, we examined microsatellite coverage and motif patterns from 454 shotgun sequences of 154 Eukaryote species from eight distantly related phyla (Cnidaria, Arthropoda, Onychophora, Bryozoa, Mollusca, Echinodermata, Chordata and Streptophyta) to test if a consistent phylogenetic pattern emerges from the microsatellite composition of these species. It is clear from our results that data from model species provide incomplete information regarding the existing microsatellite variability within the Eukaryotes. A very strong heterogeneity of microsatellite composition was found within most phyla, classes and even orders. Autocorrelation analyses indicated that while microsatellite contents of species within clades more recent than 200 Mya tend to be similar, the autocorrelation breaks down and becomes negative or non-significant with increasing divergence time. Therefore, the age of the taxon seems to be a primary factor in degrading the phylogenetic pattern present among related groups. The most recent classes or orders of Chordates still retain the pattern of their common ancestor. However, within older groups, such as classes of Arthropods, the phylogenetic pattern has been scrambled by the long independent evolution of the lineages.

Bioinformatic survey of Haliotis midae microsatellites reveals a non-random distribution of repeat motifs

Recent studies have shown the non-random distribution of microsatellite motifs between genomic regions within a particular species. This study investigates such microsatellite distributions in the genome of the economically important abalone Haliotis midae, via a bioinformatic survey. In particular, the association of specific repeat motifs to coding regions and transposable elements is investigated. An understanding of microsatellite genomic distribution will facilitate more efficient use and development of this popular molecular marker. A bias toward di- and tetranucleotide repeats was found in the H. midae genome. CA microsatellite units were the most abundant repeat motif, but were notably underrepresented in genic regions where GAGT repeats predominate. Approximately 17.5% and 21% of the microsatellites showed gene and/or transposable element associations, respectively. This could explain the high genomic frequencies of particular motifs across the genome and may allude to a possible functional role. The data presented in this study are the first to demonstrate such non-random dispersal of microsatellites in abalone and support previous findings arguing in favor of non-random distribution of repeat motifs.

Rise of the machines--recommendations for ecologists when using next generation sequencing for microsatellite development

Next generation sequencing is revolutionizing molecular ecology by simplifying the development of molecular genetic markers, including microsatellites. Here, we summarize the results of the large-scale development of microsatellites for 54 nonmodel species using next generation sequencing and show that there are clear differences amongst plants, invertebrates and vertebrates for the number and proportion of motif types recovered that are able to be utilized as markers. We highlight that the heterogeneity within each group is very large. Despite this variation, we provide an indication of what number of sequences and consequent proportion of a 454 run are required for the development of 40 designable, unique microsatellite loci for a typical molecular ecological study. Finally, to address the challenges of choosing loci from the vast array of microsatellite loci typically available from partial genome runs (average for this study, 2341 loci), we provide a microsatellite development flowchart as a procedural guide for application once the results of a partial genome run are obtained.

Mining and characterization of sequence tagged microsatellites from the brown planthopper Nilaparvata lugens

The brown planthopper, Nilaparvata lugens (Stål) (Hemiptera: Delphacidae), is an important pest of rice. To better understand the migration pattern and population structure of the Chinese populations of N. lugens, we developed and characterized 12 polymorphic microsatellites from the expressed sequence tags database of N. lugens. The occurrence of these simple sequence repeats was assessed in three populations collected from three provinces of China. The number of alleles per locus ranged from 3 to 13 with an average of 6.5 alleles per locus. The mean observed heterozygosity of the three populations ranged from 0.051 to 0.772 and the expected heterozygosity ranged from 0.074 to 0.766. The sequences of the 12 markers were highly variable. The polymorphism information content of the 12 markers was high and ranged from 0.074 to 0.807 (mean = 0.503). Sequencing of microsatellite alleles revealed that the fragment length differences were mainly due to the variation of the repeat motif. Significant genetic differentiation was detected among the three N. lugens populations as the Fst ranged from 0.034 to 0.273. Principle coordinates analysis also revealed significant genetic differentiation between populations of different years. We conclude that these microsatellite markers will be a powerful tools to study the migration routine of the N. lugens.

Molecular characterization and evolution of an interspersed repetitive DNA family of oysters

When genomic DNA from the European flat oyster Ostrea edulis L. was digested by BclI enzyme, a band of about 150 bp was observed in agarose gel. After cloning and sequencing this band and analysing their molecular characteristics and genomic organization by means of Southern blot, in situ hybridisation, and polymerase chain reaction (PCR) protocols, we concluded that this band is an interspersed highly repeated DNA element, which is related in sequence to the flanking regions of (CT)-microsatellite loci of the species O. edulis and Crassostrea gigas. Furthermore, we determined that this element forms part of a longer repetitive unit of 268 bp in length that, at least in some loci, is present in more than one copy. By Southern blot hybridisation and PCR amplifications-using primers designed for conserved regions of the 150-bp BclI clones of O. edulis-we determined that this repetitive DNA family is conserved in five other oyster species (O. stentina, C. angulata, C. gigas, C. ariakensis, and C. sikamea) while it is apparently absent in C. gasar. Finally, based on the analysis of the repetitive units in these oyster species, we discuss the slow degree of concerted evolution in this interspersed repetitive DNA family and its use for phylogenetic analysis.

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

Background

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

Results

Using 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)_n and (CT)_n were 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)_n and (ATT)_n being 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.

Conclusions

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

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

Background

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

Results

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

Conclusions

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

Genome-based microsatellite development in the Culex pipiens complex and comparative microsatellite frequency with Aedes aegypti and Anopheles gambiae

Background

Mosquitoes in the Culex pipiens complex are among the most medically important vectors for human disease worldwide and include major vectors for lymphatic filariasis and West Nile virus transmission. However, detailed genetic studies in the complex are limited by the number of genetic markers available. Here, we describe methods for the rapid and efficient identification and development of single locus, highly polymorphic microsatellite markers for Cx. pipiens complex mosquitoes via in silico screening of the Cx. quinquefasciatus genome sequence.

Methodology/Principal Findings

Six lab colonies representing four Cx. pipiens and two Cx. quinquefasciatus populations were utilized for preliminary assessment of 38 putative loci identified within 16 Cx. quinquefasciatus supercontig assemblies (CpipJ1) containing previously mapped genetic marker sequences. We identified and validated 12 new microsatellite markers distributed across all three linkage groups that amplify consistently among strains representing the complex. We also developed four groups of 3–5 microsatellite loci each for multiplex-ready PCR. Field collections from three cities in Indiana were used to assess the multiplex groups for their application to natural populations. All were highly polymorphic (Mean  = 13.0 alleles) per locus and reflected high polymorphism information content (PIC) (Mean  = 0.701). Pairwise F_ST indicated population structuring between Terre Haute and Fort Wayne and between Terre Haute and Indianapolis, but not between Fort Wayne and Indianapolis. In addition, we performed whole genome comparisons of microsatellite motifs and abundance between Cx. quinquefasciatus and the primary vectors for dengue virus and malaria parasites, Aedes aegypti and Anopheles gambiae, respectively.

Conclusions/Significance

We demonstrate a systematic approach for isolation and validation of microsatellites for the Cx. pipiens complex by direct screen of the Cx. quinquefasciatus genome supercontig assemblies. The genome density of microsatellites is greater in Cx. quinquefasciatus (0.26%) than in Ae. aegypti (0.14%), but considerably lower than in An. gambiae (0.77%).

Generation of microsatellite repeat families by RTE retrotransposons in lepidopteran genomes

Background

Developing lepidopteran microsatellite DNA markers can be problematical, as markers often exhibit multiple banding patterns and high frequencies of non-amplifying "null" alleles. Previous studies identified sequences flanking simple sequence repeat (SSR) units that are shared among many lepidopteran species and can be grouped into microsatellite-associated DNA families. These families are thought to be associated with unequal crossing-over during DNA recombination or with transposable elements (TEs).

Results

We identified full-length lepidopteran non-LTR retrotransposable elements of the RTE clade in Heliconius melpomene and Bombyx mori. These retroelements possess a single open reading frame encoding the Exonuclease/Endonuclease/Phosphatase and the Reverse Transcriptase/nLTR domains, a 5' UTR (untranslated region), and an extremely short 3' UTR that regularly consists of SSR units. Phylogenetic analysis supported previous suggestions of horizontal transfer among unrelated groups of organisms, but the diversity of lepidopteran RTE elements appears due to ancient divergence of ancestral elements rather than introgression by horizontal transfer. Similarity searches of lepidopteran genomic sequences in GenBank identified partial RTE elements, usually consisting of the 3' terminal region, in 29 species. Furthermore, we identified the C-terminal end of the Reverse Transcriptase/nLTR domain and the associated 3' UTR in over 190 microsatellite markers from 22 lepidopteran species, accounting for 10% of the lepidopteran microsatellites in GenBank. Occasional retrotransposition of autonomous elements, frequent retrotransposition of 3' partial elements, and DNA replication slippage during retrotransposition offers a mechanistic explanation for the association of SSRs with RTE elements in lepidopteran genomes.

Conclusions

Non-LTR retrotransposable elements of the RTE clade therefore join a diverse group of TEs as progenitors of SSR units in various organisms. When microsatellites are isolated using standard SSR enrichment protocols and primers designed at complementary repeated regions, amplification from multiple genomic sites can cause scoring difficulties that compromise their utility as markers. Screening against RTE elements in the isolation procedure provides one strategy for minimizing this problem.

Comparative genomic analysis reveals species-dependent complexities that explain difficulties with microsatellite marker development in molluscs

Reliable population DNA molecular markers are difficult to develop for molluscs, the reasons for which are largely unknown. Identical protocols for microsatellite marker development were implemented in three gastropods. Success rates were lower for Gibbula cineraria compared to Littorina littorea and L. saxatilis. Comparative genomic analysis of 47.2 kb of microsatellite containing sequences (MCS) revealed a high incidence of cryptic repetitive DNA in their flanking regions. The majority of these were novel, and could be grouped into DNA families based upon sequence similarities. Significant inter-specific variation in abundance of cryptic repetitive DNA and DNA families was observed. Repbase scans show that a large proportion of cryptic repetitive DNA was identified as transposable elements (TEs). We argue that a large number of TEs and their transpositional activity may be linked to differential rates of DNA multiplication and recombination. This is likely to be an important factor explaining inter-specific variation in genome stability and hence microsatellite marker development success rates. Gastropods also differed significantly in the type of TEs classes (autonomous vs non-autonomous) observed. We propose that dissimilar transpositional mechanisms differentiate the TE classes in terms of their propensity for transposition, fixation and/or silencing. Consequently, the phylogenetic conservation of non-autonomous TEs, such as CvA, suggests that dispersal of these elements may have behaved as microsatellite-inducing elements. Results seem to indicate that, compared to autonomous, non-autonomous TEs maybe have a more active role in genome rearrangement processes. The implications of the findings for genomic rearrangement, stability and marker development are discussed.

A helitron-like transposon superfamily from lepidoptera disrupts (GAAA)(n) microsatellites and is responsible for flanking sequence similarity within a microsatellite family

Transposable elements (TEs) are mobile DNA regions that alter host genome structure and gene expression. A novel 588 bp non-autonomous high copy number TE in the Ostrinia nubilalis genome has features in common with miniature inverted-repeat transposable elements (MITEs): high A + T content (62.3%), lack of internal protein coding sequence, and secondary structure consisting of subterminal inverted repeats (SIRs). The O. nubilalis TE has inserted at (GAAA)(n) microsatellite loci, and was named the microsatellite-associated interspersed nuclear element (MINE-1). Non-autonomous MINE-1 superfamily members also were identified downstream of (GAAA)(n) microsatellites within Bombyx mori and Pectinophora gossypiella genomes. Of 316 (GAAA)(n) microsatellites from the B. mori whole genome sequence, 201 (63.6%) have associated autonomous or non-autonomous MINE-1 elements. Autonomous B. mori MINE-1s a encode a helicase and endonuclease domain RepHel-like protein (BMHELp1) indicating their classification as Helitron-like transposons and were renamed Helitron1_BM. Transposition of MINE-1 members in Lepidoptera has resulted in the disruption of (GAAA)(n) microsatellite loci, has impacted the application of microsatellite-based genetic markers, and suggests genome sequence that flanks TT/AA dinucleotides may be required for target site recognition by RepHel endonuclease domains.

The distribution of microsatellites in the Nasonia parasitoid wasp genome

Microsatellites are important molecular markers used in numerous genetic contexts. Despite this widespread use, the evolutionary processes governing microsatellite distribution and diversity remain controversial. Here, we present results on the distribution of microsatellites of three species in the parasitic wasp genus Nasonia generated by an in silico data-mining approach. Our results show that the overall microsatellite density in Nasonia is comparable to that of the honey bee, but much higher than in eight non-Hymenopteran arthropods. Across the Nasonia vitripennis genome, microsatellite density varied both within and amongst chromosomes. In contrast to other taxa, dinucleotides are the most abundant repeat type in all four species of Hymenoptera studied. Whether the differences between the Hymenoptera and other taxa are of functional significance remains to be determined.