Adaptive and neutral genetic differentiation among Scottish and endangered Irish red grouse (<Emphasis Type="Italic">Lagopus lagopus scotica</Emphasis>)

Studying patterns of intra-specific genetic variation among populations allows for a better understanding of population structure and local adaptation. However, those patterns may differ according to the genetic markers applied, as neutral genetic markers reflect demographic processes and random genetic drift, whereas adaptive markers also carry the footprint of selection. In combination, neutral and adaptive genetic markers permit to assess the relative roles of drift and selection in shaping population structure. Among the best understood adaptive genetic loci are the genes of the major histocompatibility complex (MHC). We here study variation and differentiation at neutral SNP markers and MHC class II genes in red grouse (Lagopus lagopus scotica) from Ireland and Scotland. Irish red grouse populations are fragmented and drastically declining, but red grouse are abundant in Scotland. We find evidence for positive selection acting on the MHC genes and variation in MHC gene copy numbers among Irish individuals. Furthermore, there was significant population differentiation among red grouse from Ireland and Scotland at the neutral SNP markers (F_ST = 0.084) and the MHC-BLB genes (F_ST: BLB1 = 0.116, BLB2 = 0.090, BLB3 = 0.104). Differentiation at the MHC-BLB1 was significantly higher than at the neutral SNP markers, suggesting that selection plays an important role in shaping MHC variation, in addition to genetic drift. We speculate that the observed differentiation pattern might be due to local adaptation to different parasite regimes. These findings have strong conservation implications and we advise against the introduction of Scottish red grouse to supplement Irish populations.     

Insights into the neutral and adaptive processes shaping the spatial distribution of genomic variation in the economically important Moroccan locust (Dociostaurus maroccanus)

Understanding the processes that shape neutral and adaptive genomic variation is a fundamental step to determine the demographic and evolutionary dynamics of pest species. Here, we use genomic data obtained via restriction site‐associated DNA sequencing to investigate the genetic structure of Moroccan locust (Dociostaurus maroccanus) populations from the westernmost portion of the species distribution (Iberian Peninsula and Canary Islands), infer demographic trends, and determine the role of neutral versus selective processes in shaping spatial patterns of genomic variation in this pest species of great economic importance. Our analyses showed that Iberian populations are characterized by high gene flow, whereas the highly isolated Canarian populations have experienced strong genetic drift and loss of genetic diversity. Historical demographic reconstructions revealed that all populations have passed through a substantial genetic bottleneck around the last glacial maximum (~21 ka BP) followed by a sharp demographic expansion at the onset of the Holocene, indicating increased effective population sizes during warm periods as expected from the thermophilic nature of the species. Genome scans and environmental association analyses identified several loci putatively under selection, suggesting that local adaptation processes in certain populations might not be impeded by widespread gene flow. Finally, all analyses showed few differences between outbreak and nonoutbreak populations. Integrated pest management practices should consider high population connectivity and the potential importance of local adaptation processes on population persistence.     

Linking genetic and ecological differentiation in an ungulate with a circumpolar distribution

Genetic differentiation among populations may arise from the disruption of gene flow due to local adaptation to distinct environments and/or neutral accumulation of mutations and genetic drift resulted from geographical isolation. Quantifying the role of these processes in determining the genetic structure of natural populations remains challenging. Here, we analyze the relative contribution of isolation‐by‐resistance (IBR), isolation‐by‐environment (IBE), genetic drift and historical isolation in allopatry during Pleistocene glacial cycles on shaping patterns of genetic differentiation in caribou/reindeer populations Rangifer tarandus across the entire distribution range of the species. Our study integrates analyses at range‐wide and regional scales to partial out the effects of historical and contemporary isolation mechanisms. At the circumpolar scale, our results indicate that genetic differentiation is predominantly explained by IBR and historical isolation. At a regional scale, we found that IBR, IBE and population size significantly explained the spatial distribution of genetic variation among populations belonging to the Euro‐Beringian lineage within North America. In contrast, genetic differentiation among populations within the North American lineage was predominantly explained by IBR and population size, but not IBE. We also found discrepancies between genetic and ecotype designation across the Holarctic species distribution range. Overall, these results indicate that multiple isolating mechanisms have played roles in shaping the spatial distribution of genetic variation across the distribution range of a large mammal with high potential for gene flow. Considering multiple spatial scales and simultaneously testing a comprehensive suite of potential isolating mechanisms, our study contributes to understand the ecological and evolutionary processes underlying organism–landscape interactions.     

Sometimes hidden but always there: the assumptions underlying genetic inference of demographic histories

Demographic processes directly affect patterns of genetic variation within contemporary populations as well as future generations, allowing for demographic inference from patterns of both present-day and past genetic variation. Advances in laboratory procedures, sequencing and genotyping technologies in the past decades have resulted in massive increases in high-quality genome-wide genetic data from present-day populations and allowed retrieval of genetic data from archaeological material, also known as ancient DNA. This has resulted in an explosion of work exploring past changes in population size, structure, continuity and movement. However, as genetic processes are highly stochastic, patterns of genetic variation only indirectly reflect demographic histories. As a result, past demographic processes need to be reconstructed using an inferential approach. This usually involves comparing observed patterns of variation with model expectations from theoretical population genetics. A large number of approaches have been developed based on different population genetic models that each come with assumptions about the data and underlying demography. In this article I review some of the key models and assumptions underlying the most commonly used approaches for past demographic inference and their consequences for our ability to link the inferred demographic processes to the archaeological and climate records.This article is part of the theme issue ‘Cross-disciplinary approaches to prehistoric demography’.     

Estimating sex-specific processes in human populations: Are XY-homologous markers an effective tool?

Homologous markers on the sex-specific regions of the X- and Y-chromosomes are differentially inherited through males and females, and have similar molecular characteristics. They may therefore be useful as a complement to the comparison of mtDNA and Y-chromosomal haplotypes for estimating sex-specific processes shaping human population structure. To test this idea, we analyzed XY-homologous microsatellite diversity in 33 human populations from Africa, Asia and Europe. Interpopulation comparisons suggest that the generally discordant pattern of genetic variation observed for X- and Y-linked markers could be an outcome of sex-specific migration processes (m(females)/m(males) approximately 3) or sex-specific demographic processes (N(females)/N(males) approximately 11) or a combination of both. However, intrapopulation diversity estimated by the X/Y ratio Watterson estimator (theta(H(Y))/theta(H(X))) suggests that the scenarios required to explain the global genetic variation of XY-homologous markers are many and complex, and that the sex-specific processes (effective population size and migration rate) shaping human population structures are likely to be specific to each population under study. XY-homologous markers provide an insight into the genuine complexity of sex-specific processes, and their further exploitation in human population studies seems worthwhile.     

Divergence genetics analysis reveals historical population genetic processes leading to contrasting phylogeographic patterns in co-distributed species

Coalescent samplers are computational time machines for inferring the historical demographic genetic processes that have given rise to observable patterns of spatial genetic variation among contemporary populations. We have used traditional characterizations of population structure and coalescent‐based inferences about demographic processes to reconstruct the population histories of two co‐distributed marine species, the frilled dog whelk, Nucella lamellosa, and the bat star, Patiria miniata. Analyses of population structure were consistent with previous work in both species except that additional samples of N. lamellosa showed a larger regional genetic break on Vancouver Island (VI) rather than between the southern Alexander Archipelago as in P. miniata. Our understanding of the causes, rather than just the patterns, of spatial genetic variation was dramatically improved by coalescent analyses that emphasized variation in population divergence times. Overall, gene flow was greater in bat stars (planktonic development) than snails (benthic development) but spatially homogeneous within species. In both species, these large phylogeographic breaks corresponded to relatively ancient divergence times between populations rather than regionally restricted gene flow. Although only N. lamellosa shows a large break on VI, population separation times on VI are congruent between species, suggesting a similar response to late Pleistocene ice sheet expansion. The absence of a phylogeographic break in P. miniata on VI can be attributed to greater gene flow and larger effective population size in this species. Such insights put the relative significance of gene flow into a more comprehensive historical biogeographic context and have important implications for conservation and landscape genetic studies that emphasize the role of contemporary gene flow and connectivity in shaping patterns of population differentiation.     

Genetic diversity and differentiation in a wide ranging anadromous fish,American shad (Alosa sapidissima), is correlated with latitude

Studies that span entire species ranges can provide insight into the relative roles of historical contingency and contemporary factors that influence population structure and can reveal patterns of genetic variation that might otherwise go undetected. American shad is a wide ranging anadromous clupeid fish that exhibits variation in demographic histories and reproductive strategies (both semelparity and iteroparity) and provides a unique perspective on the evolutionary processes that govern the genetic architecture of anadromous fishes. Using 13 microsatellite loci, we examined the magnitude and spatial distribution of genetic variation among 33 populations across the species' range to (i) determine whether signals of historical demography persist among contemporary populations and (ii) assess the effect of different reproductive strategies on population structure. Patterns of genetic diversity and differentiation among populations varied widely and reflect the differential influences of historical demography, microevolutionary processes and anthropogenic factors across the species' range. Sequential reductions of diversity with latitude among formerly glaciated rivers are consistent with stepwise postglacial colonization and successive population founder events. Weak differentiation among U.S. iteroparous populations may be a consequence of human‐mediated gene flow, while weak differentiation among semelparous populations probably reflects natural gene flow. Evidence for an effect of reproductive strategy on population structure suggests an important role for environmental variation and suggests that the factors that are responsible for shaping American shad life history patterns may also influence population genetic structure.     

Rooting Out Genetic Structure of Invasive Wild Pigs in Texas

Invasive wild pigs (Sus scrofa), also called feral swine or wild hogs, are recognized as among the most destructive invasive species in the world. Throughout the United States, invasive wild pigs have expanded rapidly over the past 40 years with populations now established in 38 states. Of the estimated 6.9 million wild pigs distributed throughout the United States, Texas supports approximately 40% of the population and similarly bears disproportionate ecological and economic costs. Genetic analyses are an effective tool for understanding invasion pathways and tracking dispersal of invasive species such as wild pigs and have been used recently in California and Florida, USA, which have similarly long-established populations and high densities of wild pigs. Our goals were to use molecular approaches to elucidate invasion and migration processes shaping wild pig populations throughout Texas, compare our results with patterns of genetic structure observed in California and Florida, and provide insights for effective management of this invasive species. We used a high-density single nucleotide polymorphism (SNP) array to evaluate population genetic structure. Genetic clusters of wild pigs throughout Texas demonstrate 2 distinct patterns: weakly resolved, spatially dispersed clusters and well-resolved, spatially localized clusters. The disparity in patterns of genetic structure suggests disparate processes are differentially shaping wild pig populations in various localities throughout the state. Our results differed from the patterns of genetic structure observed in California and Florida, which were characterized by localized genetic clusters. These differences suggest distinct biological and perhaps anthropogenic processes are shaping genetic structure in Texas. Further, these disparities demonstrate the need for location-specific management strategies for controlling wild pig populations and mitigating associated ecological and economic costs. © 2021 The Wildlife Society. This article has been contributed to by US Government employees and their work is in the public domain in the USA.     

Postglacial expansion and not human influence best explains the population structure in the endangered kea (Nestor notabilis)

Inferring past demography is a central question in evolutionary and conservation biology. It is, however, sometimes challenging to infer the processes that shaped the current patterns of genetic variation in endangered species. Population substructuring can occur as a result of survival in several isolated refugia and subsequent recolonization processes or via genetic drift following a population decline. The kea (Nestor notabilis) is an endemic parrot widely distributed in the mountains of the South Island of New Zealand that has gone through a major human‐induced population decline during the 1860s–1970s. The aims of this study were to understand the glacial and postglacial history of kea and to determine whether the recent population decline played a role in the shaping of the current genetic variation. We examined the distribution of genetic variation, differentiation and admixture in kea using 17 microsatellites and the mitochondrial control region. Mitochondrial data showed a shallow phylogeny and a genetic distinction between the North and South of the range consistent with the three genetic clusters identified with microsatellite data. Both marker types indicated an increase in genetic isolation by geographic distance. Approximate Bayesian Computation supported a scenario of postglacial divergence from a single ancestral glacial refugium, suggesting that the contemporary genetic structure has resulted from recolonization processes rather than from a recent population decline. The recent evolutionary origin of this genetic structure suggests that each genetic cluster does not need to be considered as independent conservation units.     

The impact of founder effects, gene flow, and European admixture on native American genetic diversity

Recent studies have concluded that the global pattern of neutral genetic diversity in humans reflects a series of founder effects and population movements associated with our recent expansion out of Africa. In contrast, regional studies tend to emphasize the significance of more complex patterns of colonization, gene flow, and secondary population movements in shaping patterns of diversity. Our objective in this study is to examine how founder effects, gene flow, and European admixture have molded patterns of neutral genetic diversity in the Americas. Our strategy is to test the fit of a serial founder effects process to the pattern of neutral autosomal genetic variation and to examine the contribution of gene flow and European admixture to departures from fit. The genetic data consist of 678 autosomal microsatellite loci assayed by Wang and colleagues in 530 individuals in 29 widely distributed Native American populations. We find that previous evidence for serial founder effects in the Americas may be driven in part by high levels of European admixture in northern North America, intermediate levels in Central America, and low levels in eastern South America. Geographically patterned admixture may also account for previously reported genetic differences between Andean and Amazonian groups. Though admixture has obscured the precise details of precontact evolutionary processes, we find that genetic diversity is still largely hierarchically structured and that gene flow between neighboring groups has had surprisingly little impact on macrogeographic patterns of genetic diversity in the Americas.     

Testing the role of ancient and contemporary landscapes on structuring genetic variation in a specialist grasshopper

Understanding the processes underlying spatial patterns of genetic diversity and structure of natural populations is a central topic in evolutionary biogeography. In this study, we combine data on ancient and contemporary landscape composition to get a comprehensive view of the factors shaping genetic variation across the populations of the scrub‐legume grasshopper (Chorthippus binotatus binotatus) from the biogeographically complex region of southeast Iberia. First, we examined geographical patterns of genetic structure and employed an approximate Bayesian computation (ABC) approach to compare different plausible scenarios of population divergence. Second, we used a landscape genetic framework to test for the effects of (1) Late Miocene paleogeography, (2) Pleistocene climate fluctuations, and (3) contemporary topographic complexity on the spatial patterns of population genetic differentiation. Genetic structure and ABC analyses supported the presence of three genetic clusters and a sequential west‐to‐east splitting model that predated the last glacial maximum (LGM, c. 21 Kya). Landscape genetic analyses revealed that population genetic differentiation was primarily shaped by contemporary topographic complexity, but was not explained by any paleogeographic scenario or resistance distances based on climate suitability in the present or during the LGM. Overall, this study emphasizes the need of integrating information on ancient and contemporary landscape composition to get a comprehensive view of their relative importance to explain spatial patterns of genetic variation in organisms inhabiting regions with complex biogeographical histories.     

Complex pattern of genetic structuring in the Atlantic salmon (Salmo salar L.) of the River Foyle system in northwest Ireland: disentangling the evolutionary signal from population stochasticity

Little is known about the microevolutionary processes shaping within river population genetic structure of aquatic organisms characterized by high levels of homing and spawning site fidelity. Using a microsatellite panel, we observed complex and highly significant levels of intrariver population genetic substructure and Isolation-by-Distance, in the Atlantic salmon stock of a large river system. Two evolutionary models have been considered explaining mechanisms promoting genetic substructuring in Atlantic salmon, the member-vagrant and metapopulation models. We show that both models can be simultaneously used to explain patterns and levels of population structuring within the Foyle system. We show that anthropogenic factors have had a large influence on contemporary population structure observed. In an analytical development, we found that the frequently used estimator of genetic differentiation, F(ST), routinely underestimated genetic differentiation by a factor three to four compared to the equivalent statistic Jost's D(est) (Jost 2008). These statistics also showed a near-perfect correlation. Despite ongoing discussions regarding the usefulness of "adjusted"F(ST) statistics, we argue that these could be useful to identify and quantify qualitative differences between populations, which are important from management and conservation perspectives as an indicator of existence of biologically significant variation among tributary populations or a warning of critical environmental damage.     

Spatially explicit models of dynamic histories: examination of the genetic consequences of Pleistocene glaciation and recent climate change on the American Pika

A central goal of phylogeography is to identify and characterize the processes underlying divergence. One of the biggest impediments currently faced is how to capture the spatiotemporal dynamic under which a species evolved. Here, we described an approach that couples species distribution models (SDMs), demographic and genetic models in a spatiotemporally explicit manner. Analyses of American Pika (Ochotona princeps) from the sky islands of the central Rocky Mountains of North America are used to provide insights into key questions about integrative approaches in landscape genetics, population genetics and phylogeography. This includes (i) general issues surrounding the conversion of time-specific SDMs into simple continuous, dynamic landscapes from past to current, (ii) the utility of SDMs to inform demographic models with deme-specific carrying capacities and migration potentials as well as (iii) the contribution of the temporal dynamic of colonization history in shaping genetic patterns of contemporary populations. Our results support that the inclusion of a spatiotemporal dynamic is an important factor when studying the impact of distributional shifts on patterns of genetic data. Our results also demonstrate the utility of SDMs to generate species-specific predictions about patterns of genetic variation that account for varying degrees of habitat specialization and life history characteristics of taxa. Nevertheless, the results highlight some key issues when converting SDMs for use in demographic models. Because the transformations have direct effects on the genetic consequence of population expansion by prescribing how habitat heterogeneity and spatiotemporal variation is related to the species-specific demographic model, it is important to consider alternative transformations when studying the genetic consequences of distributional shifts.     

Colonization history in Fennoscandian rodents

Fennoscandia probably constitutes one of the best places on earth to study faunal history. During the height of the most recent glacial period Fennoscandia was completely covered with ice. Thus, the majority of extant species must originate from ancestors who survived the latest glaciation in non-glaciated areas outside Fennoscandia. Moreover, the geography and geological history of Fennoscandia suggests that post-glacial recolonization by land mammals must have been restricted to specific routes in time and space. Phylogeographic surveys of mitochondrial DNA (mtDNA) variation in Fennoscandian rodents have demonstrated that glacial history and patterns of post-glacial colonization have played a major role in shaping present day patterns of genetic differentiation within species. Thus, the genetic imprints of historical demographic conditions and vicariant geographic events have been retained within species and can be used to infer the history of populations. The field vole (Microtus agrestis) is used to illustrate these data and processes. Comparisons are made with phylogeographic surveys of the bank vole (Clethrionomys glareolus) , the eastern house mouse (Mus musculus) and the wood lemming (Myopus schisticolor) as well as a few other species for which less extensive studies have been performed. The main patterns of post-glacial colonization of Fennoscandia by rodents are described. The effects of timing and patterns of colonization on contemporary population genetic structure and levels of genetic variation are discussed. Specifically, the effects of hybridization and introgression as well as founder events and bottlenecks are explored.     

The influence of historical landscape change on genetic variation and population structure of a terrestrial salamander (Plethodon cinereus)

Forest loss and fragmentation is expected to shape the genetic structure of amphibian populations and reduce genetic variation. Another factor widely understood to have impacted these same parameters in North America is the range expansion that occurred following glacial retreat at the end of the Pleistocene. The Eastern Red-Backed Salamander (Plethodon cinereus) has been subjected to both processes. In this context, we investigated the historical events that are likely to have shaped genetic variation in this species using a panel of six microsatellite markers screened on individuals sampled across ten localities in northeastern Indiana, USA. We found low genetic diversity across forest patches and minimal differentiation. We expected population structure associated with forest fragmentation to result from genetic drift in isolation but instead found that a balance between gene flow and drift was ~50 times more likely. Ratios of allele number and range (M), and coalescent modeling of population demography suggested the occurrence of marked historic decline in effective population size across the region. Taken together, the data point to a loss of genetic variation which preceded deforestation over the past 200 years. This result indicates an important role for ancient demographic processes in shaping current genetic variation that may make it difficult to detect the effect of recent habitat fragmentation.     

On the relative roles of selection and genetic drift in shaping MHC variation

Genes of the major histocompatibility complex (MHC) have provided some of the clearest examples of how natural selection generates discordances between adaptive and neutral variation in natural populations. The type and intensity of selection as well as the strength of genetic drift are believed to be important in shaping the resulting pattern of MHC diversity. However, evaluating the relative contribution of multiple microevolutionary forces is challenging, and empirical studies have reported contrasting results. For instance, balancing selection has been invoked to explain high levels of MHC diversity and low population differentiation in comparison with other nuclear markers. Other studies have shown that genetic drift can sometimes overcome selection and then patterns of genetic variation at adaptive loci cannot be discerned from those occurring at neutral markers. Both empirical and simulated data also indicate that loss of genetic diversity at adaptive loci can occur faster than at neutral loci when selection and population bottlenecks act simultaneously. Diversifying selection, on the other hand, explains accelerated MHC divergence as the result of spatial variation in pathogen‐mediated selective regimes. Because of all these possible scenarios and outcomes, collecting information from as many study systems as possible, is crucial to enhance our understanding about the evolutionary forces driving MHC polymorphism. In this issue, Miller and co‐workers present an illuminating contribution by combining neutral markers (microsatellites) and adaptive MHC class I loci during the investigation of genetic differentiation across island populations of tuatara Sphenodon punctatus. Their study of geographical variation reveals a major role of genetic drift in shaping MHC variation, yet they also discuss some support for diversifying selection.     

Genetic variation and structure of house sparrow populations: is there an island effect?

Population genetic structure and intrapopulation levels of genetic variation have important implications for population dynamics and evolutionary processes. Habitat fragmentation is one of the major threats to biodiversity. It leads to smaller population sizes and reduced gene flow between populations and will thus also affect genetic structure. We use a natural system of island and mainland populations of house sparrows along the coast of Norway to characterize the different population genetic properties of fragmented populations. We genotyped 636 individuals distributed across 14 populations at 15 microsatellite loci. The level of genetic differentiation was estimated using F‐statistics and specially designed Mantel tests were conducted to study the influence of population type (i.e. mainland or island) and geographic distance on the genetic population structure. Furthermore, the effects of population type, population size and latitude on the level of genetic variation within populations were examined. Our results suggest that genetic processes on islands and mainland differed in two important ways. First, the intrapopulation level of genetic variation tended to be lower and the occurrence of population bottlenecks more frequent on islands than the mainland. Second, although the general level of genetic differentiation was low to moderate, it was higher between island populations than between mainland populations. However, differentiation increased in mainland populations somewhat faster with geographical distance. These results suggest that population bottleneck events and genetic drift have been more important in shaping the genetic composition of island populations compared with populations on the mainland. Such knowledge is relevant for a better understanding of evolutionary processes and conservation of threatened populations.     

Genetic diversity and differentiation at MHC genes in island populations of tuatara (Sphenodon spp.)

Neutral genetic markers are commonly used to understand the effects of fragmentation and population bottlenecks on genetic variation in threatened species. Although neutral markers are useful for inferring population history, the analysis of functional genes is required to determine the significance of any observed geographical differences in variation. The genes of the major histocompatibility complex (MHC) are well‐known examples of genes of adaptive significance and are particularly relevant to conservation because of their role in pathogen resistance. In this study, we survey diversity at MHC class I loci across a range of tuatara populations. We compare the levels of MHC variation with that observed at neutral microsatellite markers to determine the relative roles of balancing selection, diversifying selection and genetic drift in shaping patterns of MHC variation in isolated populations. In general, levels of MHC variation within tuatara populations are concordant with microsatellite variation. Tuatara populations are highly differentiated at MHC genes, particularly between the northern and Cook Strait regions, and a trend towards diversifying selection across populations was observed. However, overall our results indicate that population bottlenecks and isolation have a larger influence on patterns of MHC variation in tuatara populations than selection.     

Diversity and divergence patterns in regulatory genes suggest differential gene flow in recently derived species of the Hawaiian silversword alliance adaptive radiation (Asteraceae)

The impact of gene flow and population size fluctuations in shaping genetic variation during adaptive radiation, at both the genome-wide and gene-specific levels, is very poorly understood. To examine how historical population size and gene flow patterns within and between loci have influenced lineage divergence in the Hawaiian silversword alliance, we have investigated the nucleotide sequence diversity and divergence patterns of four floral regulatory genes (ASAP1-A, ASAP1-B, ASAP3-A, ASAP3-B) and a structural gene (ASCAB9). Levels and patterns of molecular divergence across these five nuclear loci were estimated between two recently derived species (Dubautia ciliolata and Dubautia arborea) which are presumed to be sibling species. This multilocus analysis of genetic variation, haplotype divergence and historical demography indicates that population expansion and differential gene flow occurred subsequent to the divergence of these two lineages. Moreover, contrasting patterns of allele- sharing for regulatory loci vs. a structural locus between these two sibling species indicate alternative histories of genetic variation and partitioning among loci where alleles of the floral regulatory loci are shared primarily from D. arborea to D. ciliolata and alleles of the structural locus are shared in both directions. Taken together, these results suggest that adaptively radiating species can exhibit contrasting allele migration rates among loci such that allele movement at specific loci may supersede genetic divergence caused by drift and that lineage divergence during adaptive radiation can be associated with population expansion.     

Revisiting the Iberian honey bee (Apis mellifera iberiensis) contact zone: maternal and genome‐wide nuclear variations provide support for secondary contact from historical refugia

Dissecting diversity patterns of organisms endemic to Iberia has been truly challenging for a variety of taxa, and the Iberian honey bee is no exception. Surveys of genetic variation in the Iberian honey bee are among the most extensive for any honey bee subspecies. From these, differential and complex patterns of diversity have emerged, which have yet to be fully resolved. Here, we used a genome‐wide data set of 309 neutrally tested single nucleotide polymorphisms (SNPs), scattered across the 16 honey bee chromosomes, which were genotyped in 711 haploid males. These SNPs were analysed along with an intergenic locus of the mtDNA, to reveal historical patterns of population structure across the entire range of the Iberian honey bee. Overall, patterns of population structure inferred from nuclear loci by multiple clustering approaches and geographic cline analysis were consistent with two major clusters forming a well‐defined cline that bisects Iberia along a northeastern–southwestern axis, a pattern that remarkably parallels that of the mtDNA. While a mechanism of primary intergradation or isolation by distance could explain the observed clinal variation, our results are more consistent with an alternative model of secondary contact between divergent populations previously isolated in glacial refugia, as proposed for a growing list of other Iberian taxa. Despite current intense honey bee management, human‐mediated processes have seemingly played a minor role in shaping Iberian honey bee genetic structure. This study highlights the complexity of the Iberian honey bee patterns and reinforces the importance of Iberia as a reservoir of Apis mellifera diversity.