Isolation of microsatellite markers from the Adonis blue butterfly (Lysandra bellargus)

Population Genetic Structure and Demography of the Critically Endangered Chequered Blue Butterfly (Scolitantides orion) in a Highly Isolated Part of Its Distribution Range

Scolitantides orion is a butterfly species threatened in many European countries. In Poland, it survived in a single highly isolated area (Vistula River valley), which is an example of the dramatic decline in the population number. We studied the two largest remaining populations inhabiting opposite banks of the river. Mark-release-recapture studies showed that both populations were small, and they fluctuated in numbers, but adult individuals were twice as numerous on the western site. Genetic analyses were carried out using a mitochondrial (COI, ND5) and nuclear markers (Wgl, EF-1α, and microsatellite loci). We found out that genetic variation was low at both sites but higher in the smaller eastern population. This pattern is likely to be better explained by past distribution, when the butterfly, as a continental species used to be much more widespread in the east. However, the genetic differentiation between populations was low. This could suggest that the existing gene flow is facilitated by dominant regional wind direction, which may also contribute to a better genetic condition of the western population. Finally, a comparison of the obtained COI sequences with others available enabled us to reveal the phylogeographic pattern of the S. orion from different localities within its range.

32 species validation of a new Illumina paired-end approach for the development of microsatellites

Development and optimization of novel species-specific microsatellites, or simple sequence repeats (SSRs) remains an important step for studies in ecology, evolution, and behavior. Numerous approaches exist for identifying new SSRs that vary widely in terms of both time and cost investments. A recent approach of using paired-end Illumina sequence data in conjunction with the bioinformatics pipeline, PAL_FINDER, has the potential to substantially reduce the cost and labor investment while also improving efficiency. However, it does not appear that the approach has been widely adopted, perhaps due to concerns over its broad applicability across taxa. Therefore, to validate the utility of the approach we developed SSRs for 32 species representing 30 families, 25 orders, 11 classes, and six phyla and optimized SSRs for 13 of the species. Overall the IPE method worked extremely well and we identified 1000s of SSRs for all species (mean = 128,485), with 17% of loci being potentially amplifiable loci, and 25% of these met our most stringent criteria designed to that avoid SSRs associated with repetitive elements. Approximately 61% of screened primers yielded strong amplification of a single locus.

A meta-analysis of dispersal in butterflies

Dispersal has recently gained much attention because of its crucial role in the conservation and evolution of species facing major environmental changes such as habitat loss and fragmentation, climate change, and their interactions. Butterflies have long been recognized as ideal model systems for the study of dispersal and a huge amount of data on their ability to disperse has been collected under various conditions. However, no single 'best' method seems to exist leading to the co-occurrence of various approaches to study butterfly mobility, and therefore a high heterogeneity among data on dispersal across this group. Accordingly, we here reviewed the knowledge accumulated on dispersal and mobility in butterflies, to detect general patterns. This meta-analysis specifically addressed two questions. Firstly, do the various methods provide a congruent picture of how dispersal ability is distributed across species? Secondly, is dispersal species-specific? Five sources of data were analysed: multisite mark-recapture experiments, genetic studies, experimental assessments, expert opinions, and transect surveys. We accounted for potential biases due to variation in genetic markers, sample sizes, spatial scales or the level of habitat fragmentation. We showed that the various dispersal estimates generally converged, and that the relative dispersal ability of species could reliably be predicted from their relative vagrancy (records of butterflies outside their normal habitat). Expert opinions gave much less reliable estimates of realized dispersal but instead reflected migration propensity of butterflies. Within-species comparisons showed that genetic estimates were relatively invariable, while other dispersal estimates were highly variable. This latter point questions dispersal as a species-specific, invariant trait.

Genetic differentiation across a latitudinal gradient in two co-occurring butterfly species: revealing population differences in a context of climate change

Genetic differentiation within a species' range is determined by natural selection, genetic drift, and gene flow. Selection and drift enhance genetic differences if populations are sufficiently isolated, while gene flow precludes differentiation and local adaptation. Over large geographical areas, these processes can create a variety of scenarios, ranging from admixture to a high degree of population differentiation. Genetic differences among populations may signal functional differences within a species' range, potentially leading to population or ecotype-specific responses to global change. We investigated differentiation within the geographical range of two butterfly species along a broad latitudinal gradient. This gradient is the primary axis of climatic variation, and many ecologists expect populations at the poleward edge of this gradient to expand under climate change. Our study species inhabit a shared ecosystem and differ in body size and resource specialization; both also find their poleward range limit on an island. We find evidence for divergence of peripheral populations from the core in both taxa, suggesting the potential for genetic distinctiveness at the leading edge of climate change. We also find differences between the species in the extent of peripheral differentiation with the smaller and more specialized species showing greater population divergence (microsatellites and mtDNA) and reduced gene flow (mtDNA). Finally, gene flow estimates in both species differed strongly between two marker types. These findings suggest caution in assuming that populations are invariant across latitude and thus will respond as a single ecotype to climatic change.

Microsatellite flanking region similarities among different loci within insect species

Although microsatellites are ubiquitous in eukaryota, the number of available markers varies strongly among taxa. This meta-analysis was conducted on 32 insect species. Sequences were obtained from two assembled whole genomes, whole genome shotgun (WGS) sequences from 10 species and screening partial genomic libraries for microsatellites from 23 species. We have demonstrated: (1) strong differences in the abundance of microsatellites among species; (2) that microsatellites within species are often grouped into families based on similarities in their flanking sequences; (3) that the proportion of microsatellites grouped into families varies strongly among taxa; and (4) that microsatellite families were significantly more often associated with transposable elements - or their remnants - than unique microsatellite sequences.

Survey and analysis of microsatellites in the silkworm, Bombyx mori: frequency, distribution, mutations, marker potential and their conservation in heterologous species

We studied microsatellite frequency and distribution in 21.76-Mb random genomic sequences, 0.67-Mb BAC sequences from the Z chromosome, and 6.3-Mb EST sequences of Bombyx mori. We mined microsatellites of >/=15 bases of mononucleotide repeats and >/=5 repeat units of other classes of repeats. We estimated that microsatellites account for 0.31% of the genome of B. mori. Microsatellite tracts of A, AT, and ATT were the most abundant whereas their number drastically decreased as the length of the repeat motif increased. In general, tri- and hexanucleotide repeats were overrepresented in the transcribed sequences except TAA, GTA, and TGA, which were in excess in genomic sequences. The Z chromosome sequences contained shorter repeat types than the rest of the chromosomes in addition to a higher abundance of AT-rich repeats. Our results showed that base composition of the flanking sequence has an influence on the origin and evolution of microsatellites. Transitions/transversions were high in microsatellites of ESTs, whereas the genomic sequence had an equal number of substitutions and indels. The average heterozygosity value for 23 polymorphic microsatellite loci surveyed in 13 diverse silkmoth strains having 2-14 alleles was 0.54. Only 36 (18.2%) of 198 microsatellite loci were polymorphic between the two divergent silkworm populations and 10 (5%) loci revealed null alleles. The microsatellite map generated using these polymorphic markers resulted in 8 linkage groups. B. mori microsatellite loci were the most conserved in its immediate ancestor, B. mandarina, followed by the wild saturniid silkmoth, Antheraea assama.

High similarity between flanking regions of different microsatellites detected within each of two species of Lepidoptera: Parnassius apollo and Euphydryas aurinia

Microsatellite flanking regions have been compared in two butterfly species. Several microsatellite flanking regions showed high similarity to one another among different microsatellites within a same species, but very few similarities were found between species. This can be the consequence of either duplication/multiplication events involving large regions containing microsatellites or of microsatellites imbedded in minisatellite regions. The multiplication of microsatellites might also be linked to mobile elements. Furthermore, crossing over between nonhomologous microsatellites can lead to the exchange of the flanking regions between microsatellites. The same phenomenon was observed in both studied butterfly species but not in Aphis fabae (Hemiptera), which was screened at the same time using the same protocol. These findings might explain, at least partially, why microsatellite isolation in Lepidoptera has been relatively unsuccessful so far.

Microsatellite markers to assess the influence of population size, isolation and demographic change on the genetic structure of the UK butterflyPolyommatus bellargus

Five microsatellite DNA markers were isolated and used to quantify population genetic structure among a subset of UK populations of the Adonis blue (Polyommatus bellargus Rottemburg). Specifically, whether population size, degree of isolation or history of bottlenecking in 1976-1978 can explain current patterns of genetic variation. The butterfly is at its northern range limit in the UK, where it exists as a highly fragmented metapopulation on isolated pockets of calcareous grassland. Most populations were affected by a severe bottleneck in the late 1970s, when a drought caused the host plant (Hippocrepis comosa) to wilt. Mantel tests and spatial autocorrelation analysis indicated a significant effect of isolation by distance among the UK populations, a relationship that broke down at greater geographical scales (> 23.85 km), probably because of large areas of unsuitable habitat presenting barriers to gene flow. Similarly, amova revealed that variation among geographical regions was almost double that observed within regions. Larger populations were found to support significantly higher levels of genetic diversity, suggesting that small populations may lose genetic diversity through drift. If, as in other butterfly species, low genetic diversity increases the probability of population extinction, then these populations are likely to be under threat. Neither isolation nor a history of bottlenecks appeared to influence genetic diversity. The results indicate that adequate population size a crucial factor in the conservation of genetic diversity in P. bellargus in the UK.