Urban land use limits regional bumble bee gene flow

Potential declines in native pollinator communities and increased reliance on pollinator‐dependent crops have raised concerns about native pollinator conservation and dispersal across human‐altered landscapes. Bumble bees are one of the most effective native pollinators and are often the first to be extirpated in human‐altered habitats, yet little is known about how bumble bees move across fine spatial scales and what landscapes promote or limit their gene flow. In this study, we examine regional genetic differentiation and fine‐scale relatedness patterns of the yellow‐faced bumble bee, Bombus vosnesenskii, to investigate how current and historic habitat composition impact gene flow. We conducted our study across a landscape mosaic of natural, agricultural and urban/suburban habitats, and we show that B. vosnesenskii exhibits low but significant levels of differentiation across the study system (F_ST = 0.019, D_est = 0.049). Most importantly, we reveal significant relationships between pairwise F_ST and resistance models created from contemporary land use maps. Specifically, B. vosnesenskii gene flow is most limited by commercial, industrial and transportation‐related impervious cover. Finally, our fine‐scale analysis reveals significant but declining relatedness between individuals at the 1–9 km spatial scale, most likely due to local queen dispersal. Overall, our results indicate that B. vosnesenskii exhibits considerable local dispersal and that regional gene flow is significantly limited by impervious cover associated with urbanization.     

Fine‐scale spatial genetic structure of common and declining bumble bees across an agricultural landscape

Land‐use changes have threatened populations of many insect pollinators, including bumble bees. Patterns of dispersal and gene flow are key determinants of species' ability to respond to land‐use change, but have been little investigated at a fine scale (<10 km) in bumble bees. Using microsatellite markers, we determined the fine‐scale spatial genetic structure of populations of four common Bombus species (B. terrestris, B. lapidarius, B. pascuorum and B. hortorum) and one declining species (B. ruderatus) in an agricultural landscape in Southern England, UK. The study landscape contained sown flower patches representing agri‐environment options for pollinators. We found that, as expected, the B. ruderatus population was characterized by relatively low heterozygosity, number of alleles and colony density. Across all species, inbreeding was absent or present but weak (F_IS = 0.01–0.02). Using queen genotypes reconstructed from worker sibships and colony locations estimated from the positions of workers within these sibships, we found that significant isolation by distance was absent in B. lapidarius, B. hortorum and B. ruderatus. In B. terrestris and B. pascuorum, it was present but weak; for example, in these two species, expected relatedness of queens founding colonies 1 m apart was 0.02. These results show that bumble bee populations exhibit low levels of spatial genetic structure at fine spatial scales, most likely because of ongoing gene flow via widespread queen dispersal. In addition, the results demonstrate the potential for agri‐environment scheme conservation measures to facilitate fine‐scale gene flow by creating a more even distribution of suitable habitats across landscapes.     

Bumble bee pollen use and preference across spatial scales in human‐altered landscapes

相似文献   

Bumble bee colony dynamics: quantifying the importance of land use and floral resources for colony growth and queen production

Bumble bee (Bombus) species are ecologically and economically important pollinators, and many species are in decline. In this article, we develop a mechanistic model to analyse growth trajectories of Bombus vosnesenskii colonies in relation to floral resources and land use. Queen production increased with floral resources and was higher in semi‐natural areas than on conventional farms. However, the most important parameter for queen production was the colony growth rate per flower, as opposed to the average number of available flowers. This result indicates the importance of understanding mechanisms of colony growth, in order to predict queen production and enhance bumble bee population viability. Our work highlights the importance of interpreting bumble bee conservation efforts in the context of overall population dynamics and provides a framework for doing so.     

Distance,elevation and environment as drivers of diversity and divergence in bumble bees across latitude and altitude

Identifying drivers of dispersal limitation and genetic differentiation is a key goal in biogeography. We examine patterns of population connectivity and genetic diversity using restriction site‐associated DNA sequencing (RADseq) in two bumble bee species, Bombus vosnesenskii and Bombus bifarius, across latitude and altitude in mountain ranges from California, Oregon and Washington, U.S.A. Bombus vosnesenskii, which occurs across a broader elevational range at most latitudes, exhibits little population structure while B. bifarius, which occupies a relatively narrow higher elevation niche across most latitudes, exhibits much stronger population differentiation, although gene flow in both species is best explained by isolation with environmental niche resistance. A relationship between elevational habitat breadth and genetic diversity is also apparent, with B. vosnesenskii exhibiting relatively consistent levels of genetic diversity across its range, while B. bifarius has reduced genetic diversity at low latitudes, where it is restricted to high‐elevation habitat. The results of this study highlight the importance of the intersect between elevational range and habitat suitability in influencing population connectivity and suggest that future climate warming will have a fragmenting effect even on populations that are presently well connected, as they track their thermal niches upward in montane systems.     

Effect of Yangtze River on population genetic structure of the relict plant Parrotia subaequalis in eastern China

Parrotia subaequalis (Hamamelidaceae) is a Tertiary relic species endemic in eastern China. We used inter‐simple sequence repeat (ISSR) markers to access genetic diversity and population genetic structure in natural five populations of P. subaequalis. The levels of genetic diversity were higher at species level (H = 0.2031) but lower at population level (H = 0.1096). The higher genetic diversity at species levels might be attributed to the accumulation of distinctive genotypes which adapted to the different habitats after Quaternary glaciations. Meanwhile, founder effects on the early stage, and subsequent bottleneck of population regeneration due to its biological characteristics, environmental features, and human activities, seemed to explain the low population levels of genetic diversity. The hierarchical AMOVA revealed high levels (42.60%) of among‐population genetic differentiation, which was in congruence with the high levels of Nei's genetic differentiation index (G_ST = 0.4629) and limited gene flow (N_m = 0.5801) among the studied populations. Mantel test showed a significant isolation‐by‐distance, indicating that geographic isolation has a significant effect on genetic structure in this species. Unweighted pair‐group method with arithmetic average clustering, PCoA, and Bayesian analyses uniformly recovered groups that matched the geographical distribution of this species. In particular, our results suggest that Yangtze River has served as a natural barrier to gene flow between populations occurred on both riversides. Concerning the management of P. subaequalis, the high genetic differentiation among populations indicates that preserving all five natural populations in situ and collecting enough individuals from these populations for ex situ conservation are necessary.     

Reduced abundance and earlier collection of bumble bee workers under intensive cultivation of a mass‐flowering prairie crop

One of the most commonly seeded crops in Canada is canola, a cultivar of oilseed rape (Brassica napus). As a mass‐flowering crop grown intensively throughout the Canadian Prairies, canola has the potential to influence pollinator success across tens of thousands of square kilometers of cropland. Bumble bees (Bombus sp.) are efficient pollinators of many types of native and crop plants. We measured the influence of this mass‐flowering crop on the abundance and phenology of bumble bees, and on another species of social bee (a sweat bee; Halictus rubicundus), by continuously deploying traps at different levels of canola cultivation intensity, spanning the start and end of canola bloom. Queen bumble bees were more abundant in areas with more canola cover, indicating that this crop is attractive to queens. However, bumble bee workers were significantly fewer in these locations later in the season, suggesting reduced colony success. The median collection dates of workers of three bumble bee species were earlier near canola fields, suggesting a dynamic response of colonies to the increased floral resources. Different species experienced this shift to different extents. The sweat bee was not affected by canola cultivation intensity. Our findings suggest that mass‐flowering crops such as canola are attractive to bumble bee queens and therefore may lead to higher rates of colony establishment, but also that colonies established near this crop may be less successful. We propose that the effect on bumble bees can be mitigated by spacing the crop more evenly with respect to alternate floral resources.     

Annual flower strips support pollinators and potentially enhance red clover seed yield

Ecological intensification provides opportunity to increase agricultural productivity while minimizing negative environmental impacts, by supporting ecosystem services such as crop pollination and biological pest control. For this we need to develop targeted management solutions that provide critical resources to service‐providing organisms at the right time and place. We tested whether annual strips of early flowering phacelia Phacelia tanacetifolia support pollinators and natural enemies of seed weevils Protapion spp., by attracting and offering nectar and pollen before the crop flowers. This was expected to increase yield of red clover Trifolium pratense seed. We monitored insect pollinators, pests, natural enemies and seed yields in a total of 50 clover fields along a landscape heterogeneity gradient, over 2 years and across two regions in southern Sweden. About half of the fields were sown with flower strips of 125–2,000 m². The clover fields were pollinated by 60% bumble bees Bombus spp. and 40% honey bees Apis mellifera. The clover seed yield was negatively associated with weevil density, but was unrelated to bee species richness and density. Flower strips enhanced bumble bees species richness in the clover fields, with the strongest influence in heterogeneous landscapes. There were few detectable differences between crop fields with and without flower strips. However, long‐tongued bumble bees were redistributed toward field interiors and during phacelia bloom honey bees toward field edges. Clover seed yield also increased with increasing size of the flower strip. We conclude that annual flower strips of early flower resources can support bumble bee species richness and, if sufficiently large, possibly also increase crop yields. However, clover seed yield was mainly limited by weevil infestation, which was not influenced by the annual flower strips. A future goal should be to design targeted measures for pest control.     

High levels of fecundity in small and isolated populations of a self‐compatible Aloe pollinated by opportunistic birds and bees

• Plant species that are effective colonisers of transient habitats are expected to have a capacity for uniparental reproduction and show flexibility in pollination systems. Such traits may enable populations to be established from a small number of founding individuals without these populations succumbing to reductions in fecundity arising from pollinator limitation.
• We tested these predictions for Aloe thraskii (Xanthorrhoeaceae), a succulent treelet that colonises shifting coastal dunes and has both bird and bee pollinators. We performed hand‐pollination experiments, and selectively excluded bird visitors to determine differences in pollinator effectiveness. We measured pollinator visitation rates and fecundity in populations varying in their size, density and isolation distance.
• Controlled hand‐pollinations revealed that unlike most other Aloe species, A. thraskii is self‐compatible and thus capable of uniparental reproduction. The species does however depend on pollinators and is visited by various bird species as well as by bees. Fruit and seed set are not affected by selective exclusion of birds, thus indicating that bees are effective pollinators. Bird visitation rates increased with increasing plant height and population size, while bee visitation rates increased with increasing population size and density. We found that seed set per flower was lower in large populations than in small populations.
• These results suggest that establishment of populations of A. thraskii from a small number of individuals is unlikely to be limited by the fecundity of individual plants.

     

Evolution of resistance to single and combined floral phytochemicals by a bumble bee parasite

Repeated exposure to inhibitory compounds can drive the evolution of resistance, which weakens chemical defence against antagonists. Floral phytochemicals in nectar and pollen have antimicrobial properties that can ameliorate infection in pollinators, but evolved resistance among parasites could diminish the medicinal efficacy of phytochemicals. However, multicompound blends, which occur in nectar and pollen, present simultaneous chemical challenges that may slow resistance evolution. We assessed evolution of resistance by the common bumble bee gut parasite Crithidia bombi to two floral phytochemicals, singly and combined, over 6 weeks (~100 generations) of chronic exposure. Resistance of C. bombi increased under single and combined phytochemical exposure, without any associated costs of reduced growth under phytochemical‐free conditions. After 6 weeks’ exposure, phytochemical concentrations that initially inhibited growth by > 50%, and exceeded concentrations in floral nectar, had minimal effects on evolved parasite lines. Unexpectedly, the phytochemical combination did not impede resistance evolution compared to single compounds. These results demonstrate that repeated phytochemical exposure, which could occur in homogeneous floral landscapes or with therapeutic phytochemical treatment of managed hives, can cause rapid evolution of resistance in pollinator parasites. We discuss possible explanations for submaximal phytochemical resistance in natural populations. Evolved resistance could diminish the antiparasitic value of phytochemical ingestion, weakening an important natural defence against infection.     

Hybridization between Felis silvestris silvestris and Felis silvestris catus in two contrasted environments in France

European wildcat (Felis silvestris silvestris) populations are fragmented throughout most of the whole range of the subspecies and may be threatened by hybridization with the domestic cat F.s. catus. The underlying ecological processes promoting hybridization remain largely unknown. In France, wildcats are mainly present in the northeast and signs of their presence in the Pyrenees have been recently provided. However, no studies have been carried out in the French Pyrenees to assess their exposure to hybridization. We compared two local populations of wildcats, one living in a continuous forest habitat in the French Pyrenees, the other living in a highly fragmented forest‐agricultural landscape in northeastern France to get insights into the variability of hybridization rates. Strong evidence of hybridization was detected in northeastern France and not in the Pyrenees. Close kin in the Pyrenees were not found in the same geographic location contrary to what was previously reported for females in the northeastern wildcat population. The two wildcat populations were significantly differentiated (F_ST = 0.072) to an extent close to what has been reported (F_ST = 0.103) between the Iberian population, from which the Pyrenean population may originate, and the German population, which is connected to the northeastern population. The genetic diversity of the Pyrenean wildcats was lower than that of northeastern wildcat populations in France and in other parts of Europe. The lower hybridization in the Pyrenees may result from the continuity of natural forest habitats. Further investigations should focus on linking landscape features to hybridization rates working on local populations.     

Genetic connectivity among swarming sites in the wide ranging and recently declining little brown bat (Myotis lucifugus)

Characterizing movement dynamics and spatial aspects of gene flow within a species permits inference on population structuring. As patterns of structuring are products of historical and current demographics and gene flow, assessment of structure through time can yield an understanding of evolutionary dynamics acting on populations that are necessary to inform management. Recent dramatic population declines in hibernating bats in eastern North America from white‐nose syndrome have prompted the need for information on movement dynamics for multiple bat species. We characterized population genetic structure of the little brown bat, Myotis lucifugus, at swarming sites in southeastern Canada using 9 nuclear microsatellites and a 292‐bp region of the mitochondrial genome. Analyses of F_ST, Φ_ST, and Bayesian clustering (STRUCTURE) found weak levels of genetic structure among swarming sites for the nuclear and mitochondrial genome (Global F_ST = 0.001, P < 0.05, Global Φ_ST = 0.045, P < 0.01, STRUCTURE K = 1) suggesting high contemporary gene flow. Hierarchical AMOVA also suggests little structuring at a regional (provincial) level. Metrics of nuclear genetic structure were not found to differ between males and females suggesting weak asymmetries in gene flow between the sexes. However, a greater degree of mitochondrial structuring does support male‐biased dispersal long term. Demographic analyses were consistent with past population growth and suggest a population expansion occurred from approximately 1250 to 12,500 BP, following Pleistocene deglaciation in the region. Our study suggests high gene flow and thus a high degree of connectivity among bats that visit swarming sites whereby mainland areas of the region may be best considered as one large gene pool for management and conservation.     

High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.     

European wildcat populations are subdivided into five main biogeographic groups: consequences of Pleistocene climate changes or recent anthropogenic fragmentation?

Extant populations of the European wildcat are fragmented across the continent, the likely consequence of recent extirpations due to habitat loss and over‐hunting. However, their underlying phylogeographic history has never been reconstructed. For testing the hypothesis that the European wildcat survived the Ice Age fragmented in Mediterranean refuges, we assayed the genetic variation at 31 microsatellites in 668 presumptive European wildcats sampled in 15 European countries. Moreover, to evaluate the extent of subspecies/population divergence and identify eventual wild × domestic cat hybrids, we genotyped 26 African wildcats from Sardinia and North Africa and 294 random‐bred domestic cats. Results of multivariate analyses and Bayesian clustering confirmed that the European wild and the domestic cats (plus the African wildcats) belong to two well‐differentiated clusters (average Ф_ST = 0.159, R_st  = 0.392, P > 0.001; Analysis of molecular variance [AMOVA]). We identified from c. 5% to 10% cryptic hybrids in southern and central European populations. In contrast, wild‐living cats in Hungary and Scotland showed deep signatures of genetic admixture and introgression with domestic cats. The European wildcats are subdivided into five main genetic clusters (average Ф_ST = 0.103, R_st  = 0.143, P > 0.001; AMOVA) corresponding to five biogeographic groups, respectively, distributed in the Iberian Peninsula, central Europe, central Germany, Italian Peninsula and the island of Sicily, and in north‐eastern Italy and northern Balkan regions (Dinaric Alps). Approximate Bayesian Computation simulations supported late Pleistocene–early Holocene population splittings (from c. 60 k to 10 k years ago), contemporary to the last Ice Age climatic changes. These results provide evidences for wildcat Mediterranean refuges in southwestern Europe, but the evolution history of eastern wildcat populations remains to be clarified. Historical genetic subdivisions suggest conservation strategies aimed at enhancing gene flow through the restoration of ecological corridors within each biogeographic units. Concomitantly, the risk of hybridization with free‐ranging domestic cats along corridor edges should be carefully monitored.     

Single nucleotide polymorphism discovery in albacore and Atlantic bluefin tuna provides insights into worldwide population structure

The optimal management of the commercially important, but mostly over‐exploited, pelagic tunas, albacore (Thunnus alalunga Bonn., 1788) and Atlantic bluefin tuna (BFT; Thunnus thynnus L., 1758), requires a better understanding of population structure than has been provided by previous molecular methods. Despite numerous studies of both species, their population structures remain controversial. This study reports the development of single nucleotide polymorphisms (SNPs) in albacore and BFT and the application of these SNPs to survey genetic variability across the geographic ranges of these tunas. A total of 616 SNPs were discovered in 35 albacore tuna by comparing sequences of 54 nuclear DNA fragments. A panel of 53 SNPs yielded F_ST values ranging from 0.0 to 0.050 between samples after genotyping 460 albacore collected throughout the distribution of this species. No significant heterogeneity was detected within oceans, but between‐ocean comparisons (Atlantic, Pacific and Indian oceans along with Mediterranean Sea) were significant. Additionally, a 17‐SNP panel was developed in Atlantic BFT by cross‐species amplification in 107 fish. This limited number of SNPs discriminated between samples from the two major spawning areas of Atlantic BFT (F_ST = 0.116). The SNP markers developed in this study can be used to genotype large numbers of fish without the need for standardizing alleles among laboratories.     

Hitchhiking the high seas: Global genomics of rafting crabs

Population differentiation and diversification depend in large part on the ability and propensity of organisms to successfully disperse. However, our understanding of these processes in organisms with high dispersal ability is biased by the limited genetic resolution offered by traditional genotypic markers. Many neustonic animals disperse not only as pelagic larvae, but also as juveniles and adults while drifting or rafting at the surface of the open ocean. In theory, the heightened dispersal ability of these animals should limit opportunities for species diversification and population differentiation. To test these predictions, we used next‐generation sequencing of genomewide restriction‐site‐associated DNA tags (RADseq) and traditional mitochondrial DNA sequencing, to investigate the species‐level relationships and global population structure of Planes crabs collected from oceanic flotsam and sea turtles. Our results indicate that species diversity in this clade is low—likely three closely related species—with no evidence of cryptic or undescribed species. Moreover, our results indicate weak population differentiation among widely separated aggregations with genetic indices showing only subtle genetic discontinuities across all oceans of the world (RADseq F_ST = 0.08–0.16). The results of this study provide unprecedented resolution of the systematics and global biogeography of this group and contribute valuable information to our understanding of how theoretical dispersal potential relates to actual population differentiation and diversification among marine organisms. Moreover, these results demonstrate the limitations of single gene analyses and the value of genomic‐level resolution for estimating contemporary population structure in organisms with large, highly connected populations.     

Distribution of mating‐type alleles and genetic variability in field populations of Leptosphaeria maculans in western Canada

Leptosphaeria maculans is the most important fungal pathogen of canola (Brassica napus, oilseed rape) that causes the devastating stem canker in canola fields of western Canada. The population genetic structure of L. maculans, represented by nine subpopulations from a 6‐year period and three different provinces in western Canada, was determined using ten minisatellite markers. Isolates collected at different locations in six consecutive years had an even distribution of MAT1‐1 and MAT1‐2 across the nine subpopulations. All subpopulations of L. maculans exhibited a moderate gene diversity (H = 0.356–0.585). The majority of the genetic variation occurred within subpopulations. Approximately 8% and 4% of the variations were distributed between sampling year and location, respectively. Genetic distance (F_ST) results, using analysis of molecular variation (AMOVA), indicated that subpopulation pairing within isolates by year ranged from F_ST = 0.010 to 0.109, and the location subpopulation ranged from F_ST = 0.038 to 0.085. Bayesian clustering analyses of multiloci inferred two distinct clusters in all the subpopulations examined. This study indicates a relatively high degree of gene exchange between the different L. maculans isolates. Our results suggest that this can occur in the wide growing areas of canola fields in western Canada. This gene exchange produced different gene allele frequencies and divergence between populations.     

Spatial Patterns of Haplotype Variation in the Epiphytic Bromeliad Catopsis nutans

Identifying factors governing the origin, distribution, and maintenance of Neotropical plant diversity is an enduring challenge. To explore the complex and dynamic historical processes that shaped contemporary genetic patterns for a Central American plant species, we investigated the spatial distribution of chloroplast haplotypes of a geographically and environmentally widespread epiphytic bromeliad with wind‐dispersed seeds, Catopsis nutans, in Costa Rica. We hypothesized that genetic discontinuities occur between northwestern and southwestern Pacific slope populations, resembling patterns reported for other plant taxa in the region. Using non‐coding chloroplast DNA from 469 individuals and 23 populations, we assessed the influences of geographic and environmental distance as well as historical climatic variation on the genetic structure of populations spanning >1200 m in elevation. Catopsis nutans revealed seven haplotypes with low within‐population diversity (mean haplotype richness = 1.2) and moderate genetic structure (F_ST = 0.699). Pairwise F_ST was significantly correlated with both geographic and environmental distance. The frequency of dominant haplotypes was significantly correlated with elevation. A cluster of nine Pacific lowland populations exhibited a distinct haplotype profile and contained five of the seven haplotypes, suggesting historical isolation and limited seed‐mediated gene flow with other populations. Paleodistribution models indicated lowland and upland habitats in this region were contiguous through past climatic oscillations. Based on our paleodistribution analysis and comparable prior phylogeographic studies, the genetic signature of recent climatic oscillations are likely superimposed upon the distribution of anciently divergent lineages. Our study highlights the unique phylogeographic history of a Neotropical plant species spanning an elevation gradient.     

Do plant populations on distinct inselbergs talk to each other? A case study of genetic connectivity of a bromeliad species in an Ocbil landscape

Here, we explore the historical and contemporaneous patterns of connectivity among Encholirium horridum populations located on granitic inselbergs in an Ocbil landscape within the Brazilian Atlantic Forest, using both nuclear and chloroplast microsatellite markers. Beyond to assess the E. horridum population genetic structure, we built species distribution models across four periods (current conditions, mid‐Holocene, Last Glacial Maximum [LGM], and Last Interglacial) and inferred putative dispersal corridors using a least‐cost path analysis to elucidate biogeographic patterns. Overall, high and significant genetic divergence was estimated among populations for both nuclear and plastid DNA (Φ_ST(n) = 0.463 and Φ_ST(plastid) = 0.961, respectively, p < .001). For nuclear genome, almost total absence of genetic admixture among populations and very low migration rates were evident, corroborating with the very low estimates of immigration and emigration rates observed among E. horridum populations. Based on the cpDNA results, putative dispersal routes in Sugar Loaf Land across cycles of climatic fluctuations in the Quaternary period revealed that the populations’ connectivity changed little during those events. Genetic analyses highlighted the low genetic connectivity and long‐term persistence of populations, and the founder effect and genetic drift seemed to have been very important processes that shaped the current diversity and genetic structure observed in both genomes. The genetic singularity of each population clearly shows the need for in situ conservation of all of them.     

Differential introgression from a sister species explains high FST outlier loci within a mussel species

Scanning genomes for loci with high levels of population differentiation has become a standard of population genetics. F_ST outlier loci are most often interpreted as signatures of local selection, but outliers might arise for many other reasons too often left unexplored. Here, we tried to identify further the history and genetic basis underlying strong differentiation at F_ST outlier loci in a marine mussel. A genome scan of genetic differentiation has been conducted between Atlantic and Mediterranean populations of Mytilus galloprovincialis. The differentiation was low overall (F_ST = 0.03), but seven loci (2%) were strong F_ST outliers. We then analysed DNA sequence polymorphism at two outlier loci. The genetic structure proved to be the consequence of differential introgression of alleles from the sister‐hybridizing species Mytilus edulis. Surprisingly, the Mediterranean population was the most introgressed at these two loci, although the contact zone between the two species is nowadays localized along the Atlantic coasts of France and the British Isles. A historical contact between M. edulis and Mediterranean M. galloprovincialis should have happened during glacial periods. It proved difficult to disentangle two hypotheses: (i) introgression was adaptive, implying edulis alleles have been favoured in Mediterranean populations, or (ii) the genetic architecture of the barrier to edulis gene flow is different between the two M. galloprovincialis backgrounds. Five of the seven outliers between M. galloprovincialis populations were also outliers between M. edulis and Atlantic M. galloprovincialis, which would support the latter hypothesis. Differential introgression across semi‐permeable barriers to gene flow is a neglected scenario to interpret outlying loci that may prove more widespread than anticipated.