Genetic differentiation of continental and island populations of Bombus terrestris (Hymenoptera: Apidae) in Europe

Ten microsatellite loci and a partial sequence of the COII mitochondrial gene were used to investigate genetic differentiation in B. terrestris, a bumble bee of interest for its high-value crop pollination. The analysis included eight populations from the European continent, five from Mediterranean islands (six subspecies altogether) and one from Tenerife (initially described as a colour form of B. terrestris but recently considered as a separate species, B. canariensis). Eight of the 10 microsatellite loci displayed high levels of polymorphism in most populations. In B. terrestris populations, the total number of alleles detected per polymorphic locus ranged from 3 to 16, with observed allelic diversity from 3.8 ± 0.5 to 6.5 ± 1.4 and average calculated heterozygosities from 0.41 ± 0.09 to 0.65 ± 0.07. B. canariensis showed a significantly lower average calculated heterozygosity (0.12 ± 0.08) and observed allelic diversity (1.5 ± 0.04) as compared to both continental and island populations of B. terrestris. No significant differentiation was found among populations of B. terrestris from the European continent. In contrast, island populations were all significantly and most of them strongly differentiated from continental populations. B. terrestris mitochondrial DNA is characterized by a low nucleotide diversity: 0.18%± 0.07%, 0.20%± 0.04% and 0.27%± 0.04% for the continental populations, the island populations and all populations together, respectively. The only haplotype found in the Tenerife population differs by a single nucleotide substitution from the most common continental haplotype of B. terrestris. This situation, identical to that of Tyrrhenian islands populations and quite different from that of B. lucorum (15 substitutions between terrestris and lucorum mtDNA) casts doubts on the species status of B. canariensis. The large genetic distance between the Tenerife and B. terrestris populations estimated from microsatellite data result, most probably, from a severe bottleneck in the Canary island population. Microsatellite and mitochondrial DNA data call for the protection of the island populations of B. terrestris against importation of bumble bees of foreign origin which are used as crop pollinators.     

Isolation and characterization of microsatellite loci in Eurycorymbus cavaleriei,a dioecious endemic tree species in China

Eleven dinucleotide microsatellites were developed and characterized for Eurycorymbus cavaleriei, a dioecious, rare and endemic tree species in China. A genomic DNA enrichment protocol was used to isolate microsatellite loci and polymorphism was explored using 30 trees from two natural populations. The observed heterozygosity (H_O) ranged from 0.379 to 0.931, with an average of 0.625. These markers provide powerful tools for the ongoing population and conservation genetics studies of E. cavaleriei.     

Remnant populations of the regal fritillary (Speyeria idalia) in Pennsylvania: Local genetic structure in a high gene flow species

The Regal Fritillary butterfly, Speyeria idalia (Drury) (Lepidoptera: Nymphalidae), has been described as a high gene flow species. Supporting this assertion, previous studies in the Great Plains, where it is still relatively widespread, have found evidence of gene flow across hundreds of kilometers. Using mitochondrial and microsatellite loci, we examined the spatial genetic structure of a very isolated Pennsylvania population of these butterflies that occupies three separate meadows located within ten kilometers of each other. We found restricted gene flow and a distinct structure, with each meadow having a unique genetic signature. Our findings indicate that even a species that normally exhibits high gene flow may show fine-scale genetic subdivision in areas where populations have been largely extirpated.Authors contributed equally.     

Contrasting microsatellite variation between subalpine and western larch, two closely related species with different distribution patterns

Subalpine larch (Larix lyallii Parl.) and western larch (Larix occidentalis Nutt.) represent two closely related species with contrasting abundance and distribution patterns in Western North America. Genetic diversity at seven informative microsatellite loci was determined for 19 populations of subalpine larch and nine populations of western larch. Contrasting genetic diversity and patterns of population differentiation were observed between the two species. The overall within-population genetic diversity parameters were lower in subalpine larch (A = 3.2; A(P) = 3.6; H(E) = 0.418) than in western larch (A(P) = 5.51; H(E) = 0.580), a pattern that is likely related to historical or demographic factors. No evidence of interspecific hybridization was observed. Significantly more population differentiation (theta = 0.15; R(ST) = 0.07), consistent with more restricted gene flow, was observed for subalpine larch as compared to western larch (theta = 0.05; R(ST) = 0.04). Under the assumption of an infinite allele mutation model, 12 of the 19 subalpine larch populations showed signs of deviation from the mutation-drift equilibrium, which suggests Holocene population bottlenecks and fluctuations in effective population size for this species. None of the western larch populations deviated significantly from the mutation-drift equilibrium. For both species, Mantel's test revealed a significant positive relationship between geographical and genetic distances indicative of isolation by distance. A similar geographical structure was detected in both species, suggesting at least two genetically distinct glacial populations in each species. The various implications for gene conservation are discussed.     

Genetic divergence between cave and surface populations of Astyanax in Mexico (Characidae,Teleostei)

A study of genetic diversity at microsatellite loci and the mitochondrial DNA (mtDNA) cytochrome b gene was carried out to assess genetic relationships among four Mexican cave (Pachon, Sabinos, Tinaja, Chica) and four surface populations of Astyanax fasciatus (Characidae) from northeast Mexico and the Yucatan. With the exception of Chica, the cave populations were all characterized by extremely low microsatellite variability, which most likely resulted from bottleneck events. Population analyses of the microsatellite data indicated no measurable levels of gene flow between all cave and surface populations (F(ST) > 0.0707). Phylogenetic analyses of mtDNA data showed that only two cave populations - Sabinos and Tinaja - group together to the exclusion of surface populations. From the microsatellite data these cave populations cluster with the Pachon cave fish population. The mtDNA thus appears to have been replaced in Pachon because of introgressive hybridization. It is likely that these three cave populations have descended from a surface ancestor in common with current surface populations, rather than evolving recently from one of the extant surface populations. Like Pachon, the Chica population clustered with the surface populations according to mtDNA data, but was not clearly associated with either the surface or the other cave populations according to the microsatellite data. Our data indicate that the Chica population evolved recently from a surface population, and subsequently hybridized with a phylogenetically older cave population. In conclusion, both the microsatellite and mtDNA data suggest multiple origins of cave populations and the Chica and Sabinos/Tinaja/Pachon were founded after at least two independent invasions from surface populations.     

High inbreeding and low connectivity among Ambystoma texanum populations in fragmented Ohio forests

Habitat loss and fragmentation negatively impact the size and diversity of many natural populations. Woodland amphibians require connected aquatic and terrestrial habitats to complete their life cycle, and often rely on metapopulation structure for long‐term persistence. Wetland loss and deforestation fragment amphibian populations, which may result in population isolation and its negative effects. The aim of this research was to analyze the population genetic structure of small‐mouthed salamanders (Ambystoma texanum) in western Ohio, where agriculture is now the dominant land use. Salamander tail tissue was collected from eight breeding pools. Three pools occur in the same forest; the other five are in forest patches at distances ranging from 250 m to 20 km from one another. Eight microsatellite loci were amplified by PCR and genotyped for allele size. Observed heterozygosities were lower than expected in all sampled populations; the two most isolated sites (Ha1, Ha2) had the highest inbreeding coefficients. Ha2 also had the lowest mean number of alleles and was found to be genetically differentiated from populations to which our data analysis indicates it was historically connected by gene flow. The most distant site (Ha1) had the highest number of private alleles and showed genetic differentiation from other populations both historically and currently. Geographic distance between pools was strongly correlated with the number of private alleles in a population. The results suggest that population isolation results in decreased genetic diversity and that a breakdown of metapopulation structure due to landscape change may contribute to differentiation between once‐connected populations.     

Genetic isolation of insular populations of the Maghrebian bat,Myotis punicus,in the Mediterranean Basin

Aim We investigate the population genetic structure of the Maghrebian bat, Myotis punicus, between the mainland and islands to assess the island colonization pattern and current gene flow between nearby islands and within the mainland. Location North Africa and the Mediterranean islands of Corsica and Sardinia. Methods We sequenced part of the control region (HVII) of 79 bats across 11 colonies. The phylogeographical pattern was assessed by analysing molecular diversity indices, examining differentiation among populations and estimating divergence time. In addition, we genotyped 182 bats across 10 colonies at seven microsatellite loci. We used analysis of molecular variance and a Bayesian approach to infer nuclear population structure. Finally, we estimated sex‐specific dispersal between Corsica and Sardinia. Results Mitochondrial analyses indicated that colonies between Corsica, Sardinia and North Africa are highly differentiated. Within islands there was no difference between colonies, while at the continental level Moroccan and Tunisian populations were highly differentiated. Analyses with seven microsatellite loci showed a similar pattern. The sole difference was the lack of nuclear differentiation between populations in North Africa, suggesting a male‐biased dispersal over the continental area. The divergence time of Sardinian and Corsican populations was estimated to date back to the early and mid‐Pleistocene. Main conclusions Island colonization by the Maghrebian bats seems to have occurred in a stepping‐stone manner and certainly pre‐dated human colonization. Currently, open water seems to prevent exchange of bats between the two islands, despite their ability to fly and the narrowness of the strait of Bonifacio. Corsican and Sardinian populations are thus currently isolated from any continental gene pool and must therefore be considered as different evolutionarily significant units (ESU).     

Population assignment reveals low migratory connectivity in a weakly structured songbird

Understanding migratory connectivity is essential for determining the drivers behind population dynamics and for implementing effective conservation strategies for migratory species. Genetic markers provide a means to describe migratory connectivity; however, they can be uninformative for species with weak population genetic structure, which has limited their application. Here, we demonstrated a genomic approach to describing migratory connectivity in the prothonotary warbler, Protonotaria citrea, a Neotropical songbird of conservation concern. Using 26,189 single nucleotide polymorphisms (SNPs), we revealed regional genetic structure between the Mississippi River Valley and the Atlantic Seaboard with overall weak genetic differentiation among populations (F_ST = 0.0055; 95% CI: 0.0051–0.0059). Genetic variation had a stronger association with geographic rather than environmental factors, with each explaining 14.5% and 8.2% of genetic variation, respectively. By varying the numbers of genomic markers used in population assignment models with individuals of known provenance, we identified a maximum assignment accuracy (89.7% to site, 94.3% to region) using a subset of 600 highly differentiated SNPs. We then assigned samples from nonbreeding sites to breeding region and found low migratory connectivity. Our results highlight the importance of filtering markers for informative loci in models of population assignment. Quantifying migratory connectivity for weakly structured species will be useful for expanding studies to a wider range of migratory species across taxonomic groups and may contribute to a deeper understanding of the evolution of migratory strategies.     

A fair fight between molecular marker types in a seascape genetics setting

From its inception, population genetics has been nearly as concerned with the genetic data type—to which analyses are brought to bear—as it is with the analysis methods themselves. The field has traversed allozymes, microsatellites, segregating sites in multilocus alignments and, currently, single nucleotide polymorphisms (SNPs) generated by high‐throughput genomic sequencing methods, primarily whole genome sequencing and reduced representation library (RRL) sequencing. As each emerging data type has gained traction, it has been compared to existing methods, based on its relative ability to discern population structural complexity at increasing levels of resolution. However, this is usually done by comparing the gold standard in one data type to the gold standard in the new data type. These gold standards frequently differ in power and in sampling density, both across a genome and throughout a spatial range. In this issue of Molecular Ecology, D’Aloia et al. apply the high‐throughput approach as fully as possible to microsatellites, nuclear loci and SNPs genotyped through an RRL method; this is coupled with a spatially dense sampling scheme. Completing a battery of population genetics analyses across data types (including a series of down‐sampled data sets), the authors find that SNP data are slightly more sensitive to fine‐scale genetic structure, and the results are more resilient to down‐sampling than microsatellites and nonrepetitive nuclear loci. However, their results are far from an unqualified victory for RRL SNP data over all previous data types: the authors note that modest additions to the microsatellites and nuclear loci data sets may provide the necessary analytical power to delineate the fine‐scale genetic structuring identified by SNPs. As always, as the field begins to fully embrace the newest thing, good science reminds us that traditional data types are far from useless, especially when combined with a well‐designed sampling scheme.     

Genotypic characteristics in populations of Sclerotinia sclerotiorum from New York State,USA

White mould, caused by the fungus Sclerotinia sclerotiorum, is one of the most destructive diseases of beans globally. In New York State, USA, white mould causes substantial losses in soybean, snap, dry and succulent baby lima beans, which are grown successively in intensive crop rotations. Management strategies for white mould in these crops are reliant upon the prophylactic use of fungicides. No complementary information on the genetic structure of the populations of S. sclerotiorum in New York State, USA is available. Twenty isolates of S. sclerotiorum were collected from symptomatic bean plants within each of 10 fields across New York State, USA in 2014. Eight microsatellite (SSR) markers were used to characterise the genotypic diversity of the hyphal‐tipped isolates. Twenty‐four multilocus genotypes (MLGs) were detected within the population but one MLG was most prevalent. Although STRUCTURE analysis identified two subpopulations, these subpopulations were not associated with geographic location, suggesting no spatial structure to the population. In addition, the pathogen populations were predominantly clonal, with some evidence of infrequent outcrossing. These findings may assist in understanding the durability of management strategies for white mould and support the selection of representative isolates for host resistance screening for pathogen populations in the sampling area.     

Characterization of microsatellite loci for a barrenwort species (Epimedium diphyllum,Berberidaceae)

Microsatellites were identified and characterized from Epimedium diphyllum, a species of barrenworts, both attractive garden plants and valuable medicinal plants. Some Japanese species of Epimedium are threatened with extinction, and are listed in the Red Data Book of plants in Japan. Natural hybrid zones also have been reported among some taxa of Japanese Epimedium. We developed eight polymorphic microsatellite primers for population genetic analyses of E. diphyllum. The numbers of alleles per locus ranged from 10 to 20, with observed levels of heterozygosity between 0.85 and 1.00. These primer sets yielded amplification in the other three Japanese Epimedium. These markers will be valuable for conservation genetics, evolutionary biology, pharmacognostic study, and horticultural study of Epimedium.     

Isolation and characterization of microsatellites in Sparganium emersum and cross‐species amplification in the related species S. erectum

We developed seven novel polymorphic microsatellite loci for the aquatic macrophyte Sparganium emersum (Sparganiaceae). These were characterized on 62 individuals collected from nine different populations. In this set of individuals, seven to 20 alleles per locus were detected and observed heterozygosity ranged between 0.16 and 0.95. Cross‐species amplification was tested in the related species Sparganium erectum, and was successful for five of the seven microsatellite loci.     

Molecular genetic diversity in populations of the stingless bee Plebeia remota: A case study

Genetic diversity is a major component of the biological diversity of an ecosystem. The survival of a population may be seriously threatened if its genetic diversity values are low. In this work, we measured the genetic diversity of the stingless bee Plebeia remota based on molecular data obtained by analyzing 15 microsatellite loci and sequencing two mitochondrial genes. Population structure and genetic diversity differed depending on the molecular marker analyzed: microsatellites showed low population structure and moderate to high genetic diversity, while mitochondrial DNA (mtDNA) showed high population structure and low diversity in three populations. Queen philopatry and male dispersal behavior are discussed as the main reasons for these findings.     

Microsatellite differentiation between Phillip Island and mainland Australian populations of the red fox Vulpes vulpes

Predation by the red fox Vulpes vulpes is believed to be threatening the little penguin Eudyptula minor on Phillip Island in Victoria. Polymorphism at seven microsatellite loci was examined to estimate the extent of differentiation between Phillip Island and mainland populations of V. vulpes. Loss of alleles has occurred on Phillip Island where foxes first appeared = 88 years ago compared with mainland populations. Genetic differentiation between the Phillip Island and mainland populations was high. The relatively high differentiation found between the two populations could be due to either low migration rates, the effect of the composition of founder animals or both effects. Further ecological and historical information about the populations is needed to explore the likely significance of these effects.     

Isolation and characterization of highly polymorphic microsatellite markers in Hypochaeris radicata (Asteraceae)

We developed five highly polymorphic dinucleotide microsatellite loci for the grassland species Hypochaeris radicata (Asteraceae). Polymorphism of these markers was examined in six populations in the Netherlands. All loci were polymorphic in all populations. The number of alleles per locus varied between 18 and 43. Expected heterozygosity was between 0.86 and 0.91. Cross‐species amplification was tested in six Hypochaeris species and was successful for three different loci in four species. These microsatellites are a useful tool in population genetic, dispersal and metapopulation studies or in testing levels of inbreeding.     

Extremely low effective population sizes, genetic structuring and reduced genetic diversity in a threatened bumblebee species, Bombus sylvarum (Hymenoptera: Apidae)

Habitat fragmentation may severely affect survival of social insect populations as the number of nests per population, not the number of individuals, represents population size, hence they may be particularly prone to loss of genetic diversity. Erosion of genetic diversity may be particularly significant among social Hymenoptera such as bumblebees (Bombus spp.), as this group may be susceptible to diploid male production, a suggested direct cost of inbreeding. Here, for the first time, we assess genetic diversity and population structuring of a threatened bumblebee species (Bombus sylvarum) which exists in highly fragmented habitat (rather than oceanic) islands. Effective population sizes, estimated from identified sisterhoods, were very low (range 21-72) suggesting that isolated populations will be vulnerable to loss of genetic variation through drift. Evidence of significant genetic structuring between populations (theta = 0.084) was found, but evidence of a bottleneck was detected in only one population. Comparison across highly fragmented UK populations and a continental population (where this species is more widespread) revealed significant differences in allelic richness attributable to a high degree of genetic diversity in the continental population. While not directly related to population size, this is perhaps explained by the high degree of isolation between UK populations relative to continental populations. We suggest that populations now existing on isolated habitat islands were probably linked by stepping-stone populations prior to recent habitat loss.     

Fine-scaled geographical population structuring in a highly mobile marine species: the Atlantic cod

Compared with many terrestrial and freshwater environments, dispersal and interbreeding is generally much less restricted in the marine environment. We studied the tendency for a marine species, the Atlantic cod, to be sub-structured into genetically differentiated populations on a fine geographical scale. We selected a coastal area free of any obvious physical barriers and restricted sampling to a 300-km region, well within the dispersal ability of this species. Screening 10 polymorphic microsatellite loci in 6 samples we detected a weak, but consistent, differentiation at all 10 loci. The average FST over loci was small (0.0023) but highly significant statistically, demonstrating that genetically differentiated populations can arise and persist in the absence of physical barriers or great distance. We found no geographical pattern in the genetic differentiation and there was no apparent trend of isolation by distance along the coastline. These findings lend support to the notion that low levels of differentiation are due to passive transport of eggs or larvae by the ocean currents rather than to adult dispersal, the latter being strongly dependent on distance.     

Gene flow between insular, coastal and interior populations of brown bears in Alaska

The brown bears of coastal Alaska have been recently regarded as comprising from one to three distinct genetic groups. We sampled brown bears from each of the regions for which hypotheses of genetic uniqueness have been made, including the bears of the Kodiak Archipelago and the bears of Admiralty, Baranof and Chichagof (ABC) Islands in southeast Alaska. These samples were analysed with a suite of nuclear microsatellite markers. The 'big brown bears' of coastal Alaska were found to be part of the continuous continental distribution of brown bears, and not genetically isolated from the physically smaller 'grizzly bears' of the interior. By contrast, Kodiak brown bears appear to have experienced little or no genetic exchange with continental populations in recent generations. The bears of the ABC Islands, which have previously been shown to undergo little or no female-mediated gene flow with mainland populations, were found not to be genetically isolated from mainland bears. The data from the four insular populations indicate that female and male dispersal can be reduced or eliminated by water barriers of 2–4 km and 7km in width, respectively.     

Eighteen polymorphic microsatellite markers for the highly endangered Spanish imperial eagle (Aquila adalberti) and related species

Here we describe the development of 18 polymorphic microsatellite markers for the endangered Spanish imperial eagle (Aquila adalberti). Microsatellites were tested in five other raptor species. These markers were revealed as good molecular tools for genetic population studies, individual identification and parentage assessment in Spanish imperial eagle and closely related species.     

Reduced genetic diversity and low effective size in peripheral northern European catfish Silurus glanis populations

Using 10 polymorphic microsatellites and 1251 individual samples (some dating back to the early 1980s), genetic structure and effective population size in all native and introduced Swedish populations of the European wels catfish or Silurus glanis were studied. Levels of genetic variability and phylogeographic relationships were compared with data from a previous study of populations in other parts of Europe. The genetically distinct Swedish populations displayed comparably low levels of genetic variability and according to one-sample estimates based on linkage disequilibrium and sib ship-reconstruction, current local effective population sizes were lower than minimum levels recommended for short-term genetic conservation. In line with a previous suggestion of postglacial colonisation from a single refugium, all Swedish populations were assembled on a common branch in a star-shaped dendrogram together with other European populations. Two distinct subpopulations were detected in upper and lower habitats of River Emån, indicating that even minor dispersal barriers may restrict gene flow for wels in running waters. Genetic assignment of specimens encountered in the brackish Baltic Sea and in lakes where the species does not occur naturally indicated presence of long-distance sea dispersal and confirmed unauthorised translocations, respectively.