Tool for genomic selection and breeding to evolutionary adaptation: Development of a 100K single nucleotide polymorphism array for the honey bee

High‐throughput high‐density genotyping arrays continue to be a fast, accurate, and cost‐effective method for genotyping thousands of polymorphisms in high numbers of individuals. Here, we have developed a new high‐density SNP genotyping array (103,270 SNPs) for honey bees, one of the most ecologically and economically important pollinators worldwide. SNPs were detected by conducting whole‐genome resequencing of 61 honey bee drones (haploid males) from throughout Europe. Selection of SNPs for the chip was done in multiple steps using several criteria. The majority of SNPs were selected based on their location within known candidate regions or genes underlying a range of honey bee traits, including hygienic behavior against pathogens, foraging, and subspecies. Additionally, markers from a GWAS of hygienic behavior against the major honey bee parasite Varroa destructor were brought over. The chip also includes SNPs associated with each of three major breeding objectives—honey yield, gentleness, and Varroa resistance. We validated the chip and make recommendations for its use by determining error rates in repeat genotypings, examining the genotyping performance of different tissues, and by testing how well different sample types represent the queen's genotype. The latter is a key test because it is highly beneficial to be able to determine the queen's genotype by nonlethal means. The array is now publicly available and we suggest it will be a useful tool in genomic selection and honey bee breeding, as well as for GWAS of different traits, and for population genomic, adaptation, and conservation questions.     

Distinct subspecies or phenotypic plasticity? Genetic and morphological differentiation of mountain honey bees in East Africa

Identifying the forces shaping intraspecific phenotypic and genotypic divergence are of key importance in evolutionary biology. Phenotypic divergence may result from local adaptation or, especially in species with strong gene flow, from pronounced phenotypic plasticity. Here, we examine morphological and genetic divergence among populations of the western honey bee Apis mellifera in the topographically heterogeneous East African region. The currently accepted “mountain refugia hypothesis” states that populations living in disjunct montane forests belong to a different lineage than those in savanna habitats surrounding these forests. We obtained microsatellite data, mitochondrial sequences, and morphometric data from worker honey bees collected from feral colonies in three montane forests and corresponding neighboring savanna regions in Kenya. Honey bee colonies from montane forests showed distinct worker morphology compared with colonies in savanna areas. Mitochondrial sequence data did not support the existence of the two currently accepted subspecies. Furthermore, analyses of the microsatellite data with a Bayesian clustering method did not support the existence of two source populations as it would be expected under the mountain refugia scenario. Our findings suggest that phenotypic plasticity rather than distinct ancestry is the leading cause behind the phenotypic divergence observed between montane forest and savanna honey bees. Our study thus corroborates the idea that high gene flow may select for increased plasticity.     

Genomic changes underlying host specialization in the bee gut symbiont Lactobacillus Firm5

Bacteria that engage in long‐standing associations with particular hosts are expected to evolve host‐specific adaptations that limit their capacity to thrive in other environments. Consistent with this, many gut symbionts seem to have a limited host range, based on community profiling and phylogenomics. However, few studies have experimentally investigated host specialization of gut symbionts and the underlying mechanisms have largely remained elusive. Here, we studied host specialization of a dominant gut symbiont of social bees, Lactobacillus Firm5. We show that Firm5 strains isolated from honey bees and bumble bees separate into deep‐branching host‐specific phylogenetic lineages. Despite their divergent evolution, colonization experiments show that bumble bee strains are capable of colonizing the honey bee gut. However, they were less successful than honey bee strains, and competition with honey bee strains completely abolished their colonization. In contrast, honey bee strains of divergent phylogenetic lineages were able to coexist within individual bees. This suggests that both host selection and interbacterial competition play important roles in host specialization. Using comparative genomics of 27 Firm5 isolates, we found that the genomes of honey bee strains harbour more carbohydrate‐related functions than bumble bee strains, possibly providing a competitive advantage in the honey bee gut. Remarkably, most of the genes encoding carbohydrate‐related functions were not conserved among the honey bee strains, which suggests that honey bees can support a metabolically more diverse community of Firm5 strains than bumble bees. These findings advance our understanding of the genomic changes underlying host specialization.     

The prevalence of pathogens in honey bee (Apis mellifera) colonies infested with the parasitic mite Varroa jacobsoni

The prevalence of nine honey bee viruses in samples of dead adult bees from Apis mellifera colonies in the Netherlands and Germany infested with the parasitic mite Varroa jacobsoni was compared with virus incidence in uninfested colonies in Britain. In colonies with low mite populations the viruses present and their incidence during the year were similar to the results obtained from British colonies. However, in marked contrast with findings in Britain, acute paralysis virus (APV) was the primary cause of adult bee mortality in German honey bee colonies severely infested with V. jacobsoni. Dead brood from unsealed and sealed infested cells from German colonies with high mite populations also contained much APV. The evidence suggests that V. jacobsoni activates APV replication in adult bees by its feeding behaviour and transmits virus from adult honey bees to pupae. In addition, adult bees, in which APV is multiplying, transmit the virus to unsealed brood in the larval food.     

Local adaptation across a complex bioclimatic landscape in two montane bumble bee species

Understanding evolutionary responses to variation in temperature and precipitation across species ranges is of fundamental interest given ongoing climate change. The importance of temperature and precipitation for multiple aspects of bumble bee (Bombus) biology, combined with large geographic ranges that expose populations to diverse environmental pressures, make these insects well‐suited for studying local adaptation. Here, we analyzed genome‐wide sequence data from two widespread bumble bees, Bombus vosnesenskii and Bombus vancouverensis, using multiple environmental association analysis methods to investigate climate adaptation across latitude and altitude. The strongest signatures of selection were observed in B. vancouverensis, but despite unique responses between species for most loci, we detected several shared responses. Genes relating to neural and neuromuscular function and ion transport were especially evident with respect to temperature variables, while genes relating to cuticle formation, tracheal and respiratory system development, and homeostasis were associated with precipitation variables. Our data thus suggest that adaptive responses for tolerating abiotic variation are likely to be complex, but that several parallels among species can emerge even for these complex traits and landscapes. Results provide the framework for future work into mechanisms of thermal and desiccation tolerance in bumble bees and a set of genomic targets that might be monitored for future conservation efforts.     

Population structure of honey bees in the Carpathian Basin (Hungary) confirms introgression from surrounding subspecies

Carniolan honey bees (Apis mellifera carnica) are considered as an indigenous subspecies in Hungary adapted to most of the ecological and climatic conditions in this area. However, during the last decades Hungarian beekeepers have recognized morphological signs of the Italian honey bee (Apis mellifera ligustica). As the natural distribution of the honey bee subspecies can be affected by the importation of honey bee queens or by natural gene flow, we aimed at determining the genetic structure and characteristics of the local honey bee population using molecular markers. All together, 48 Hungarian and 84 foreign (Italian, Polish, Spanish, Liberian) pupae and/or workers were used for mitochondrial DNA analysis. Additionally, 53 sequences corresponding to 10 subspecies and the Buckfast hybrid were downloaded from GenBank. For the nuclear analysis, 236 Hungarian and 106 foreign honey bees were genotyped using nine microsatellites. Heterozygosity values, population‐specific alleles, FST values, principal coordinate analysis, assignment tests, structure analysis, and dendrograms were calculated. Haplotype and nucleotide diversity values showed moderate values. We found that one haplotype (H9) was dominant in Hungary. The presence of the black honey bee (Apis mellifera mellifera) was negligible, but a few individuals resembling other subspecies were identified. We proved that the Hungarian honey bee population is nearly homogeneous but also demonstrated introgression from the foreign subspecies. Both mitochondrial DNA and microsatellite analyses corroborated the observations of the beekeepers. Molecular analyses suggested that Carniolan honey bee in Hungary is slightly affected by Italian and black honey bee introgression. Genetic differences were detected between Polish and Hungarian Carniolan honey bee populations, suggesting the existence of at least two different gene pools within A. m. carnica.     

The influence of nest site selection on the population dynamics of Africanized honey bees in an urban landscape

Urban landscapes provide habitat for many species, including domesticated and feral honey bees, Apis mellifera L. (Hymenoptera: Apidae). With recent losses of managed honey bee colonies, there is increasing interest in feral honey bee colonies and their potential contribution to pollination services in agricultural, natural, and urban settings. However, in some regions the feral honey bee population consists primarily of Africanized honey bees. Africanized honey bees (AHB) are hybrids between European honey bees and the African honey bee, Apis mellifera scutellataLepeletier, and have generated economic, ecological, and human health concerns because of their aggressive behavior. In this study, we used two long‐term datasets (7–10 years) detailing the spatial and temporal distribution of AHB colonies in Tucson, AZ, USA, where feral colonies occupy a variety of cavities including water meter boxes. A stage‐structured matrix model was used to elucidate the implications of nest site selection and the effects of colony terminations on the structure and dynamics of the AHB population. Our results suggest that Tucson's AHB population is driven by a relatively small number of ‘source’ colonies that escape termination (ca. 0.165 colonies per km² or 125 colonies in total), although immigrating swarms and absconding colonies from the surrounding area may have also contributed to the stability of the Tucson AHB population. Furthermore, the structure of the population has likely been impacted by the number and spatial distribution of water meter boxes across the city. The study provides an example of how urban wildlife populations are driven by interactions among landscape structure, human management, and behavioral traits conferred by an invasive genotype.     

Effects of Imidacloprid and Varroa destructor on survival and health of European honey bees,Apis mellifera

There has been growing concern over declines in populations of honey bees and other pollinators which are a vital part to our food security. It is imperative to identify factors responsible for accelerated declines in bee populations and develop solutions for reversing bee losses. While exact causes of colony losses remain elusive, risk factors thought to play key roles are ectoparasitic mites Varroa destructor and neonicotinoid pesticides. The present study aims to investigate effects of a neonicotinoid pesticide Imidacloprid and Varroa mites individually on survivorship, growth, physiology, virus dynamics and immunity of honey bee workers. Our study provides clear evidence that the exposure to sublethal doses of Imidacloprid could exert a significantly negative effect on health and survival of honey bees. We observed a significant reduction in the titer of vitellogenin (Vg), an egg yolk precursor that regulates the honey bees development and behavior and often are linked to energy homeostasis, in bees exposed to Imidacloprid. This result indicates that sublethal exposure to neonicotinoid could lead to increased energy usage in honey bees as detoxification is a energy‐consuming metabolic process and suggests that Vg could be a useful biomarker for measuring levels of energy stress and sublethal effects of pesticides on honey bees. Measurement of the quantitative effects of different levels of Varroa mite infestation on the replication dynamic of Deformed wing virus (DWV), an RNA virus associated with Varroa infestation, and expression level of immune genes yields unique insights into how honey bees respond to stressors under laboratory conditions.     

A gene for resistance to the Varroa mite (Acari) in honey bee (Apis mellifera) pupae

Social insect colonies possess a range of defences which protect them against highly virulent parasites and colony collapse. The host–parasite interaction between honey bees (Apis mellifera) and the mite Varroa destructor is unusual, as honey bee colonies are relatively poorly defended against this parasite. The interaction has existed since the mid‐20th Century, when Varroa switched host to parasitize A. mellifera. The combination of a virulent parasite and relatively naïve host means that, without acaricides, honey bee colonies typically die within 3 years of Varroa infestation. A consequence of acaricide use has been a reduced selective pressure for the evolution of Varroa resistance in honey bee colonies. However, in the past 20 years, several natural‐selection‐based breeding programmes have resulted in the evolution of Varroa‐resistant populations. In these populations, the inhibition of Varroa's reproduction is a common trait. Using a high‐density genome‐wide association analysis in a Varroa‐resistant honey bee population, we identify an ecdysone‐induced gene significantly linked to resistance. Ecdysone both initiates metamorphosis in insects and reproduction in Varroa. Previously, using a less dense genetic map and a quantitative trait loci analysis, we have identified Ecdysone‐related genes at resistance loci in an independently evolved resistant population. Varroa cannot biosynthesize ecdysone but can acquire it from its diet. Using qPCR, we are able to link the expression of ecdysone‐linked resistance genes to Varroa's meals and reproduction. If Varroa co‐opts pupal compounds to initiate and time its own reproduction, mutations in the host's ecdysone pathway may represent a key selection tool for honey bee resistance and breeding.     

Signatures of local adaptation in lowland and highland teosintes from whole‐genome sequencing of pooled samples

Spatially varying selection triggers differential adaptation of local populations. Here, we mined the determinants of local adaptation at the genomewide scale in the two closest maize wild relatives, the teosintes Zea mays ssp parviglumis and ssp. mexicana. We sequenced 120 individuals from six populations: two lowland, two intermediate and two highland populations sampled along two altitudinal gradients. We detected 8 479 581 single nucleotide polymorphisms (SNPs) covered in the six populations with an average sequencing depth per site per population ranging from 17.0× to 32.2×. Population diversity varied from 0.10 to 0.15, and linkage disequilibrium decayed very rapidly. We combined two differentiation‐based methods, and correlation of allele frequencies with environmental variables to detect outlier SNPs. Outlier SNPs displayed significant clustering. From clusters, we identified 47 candidate regions. We further modified a haplotype‐based method to incorporate genotype uncertainties in haplotype calling, and applied it to candidate regions. We retrieved evidence for selection at the haplotype level in 53% of our candidate regions, and in 70% of the cases the same haplotype was selected in the two lowland or the two highland populations. We recovered a candidate region located within a previously characterized inversion on chromosome 1. We found evidence of a soft sweep at a locus involved in leaf macrohair variation. Finally, our results revealed frequent colocalization between our candidate regions and loci involved in the variation of traits associated with plant–soil interactions such as root morphology, aluminium and low phosphorus tolerance. Soil therefore appears to be a major driver of local adaptation in teosintes.     

The Potential Influence of Bumble Bee Visitation on Foraging Behaviors and Assemblages of Honey Bees on Squash Flowers in Highland Agricultural Ecosystems

Bee species interactions can benefit plant pollination through synergistic effects and complementary effects, or can be of detriment to plant pollination through competition effects by reducing visitation by effective pollinators. Since specific bee interactions influence the foraging performance of bees on flowers, they also act as drivers to regulate the assemblage of flower visitors. We selected squash (Cucurbita pepo L.) and its pollinators as a model system to study the foraging response of honey bees to the occurrence of bumble bees at two types of sites surrounded by a high amount of natural habitats (≥ 58% of land cover) and a low amount of natural habitats (≤ 12% of land cover) in a highland agricultural ecosystem in China. At the individual level, we measured the elapsed time from the departure of prior pollinator(s) to the arrival of another pollinator, the selection of honey bees for flowers occupied by bumble bees, and the length of time used by honey bees to explore floral resources at the two types of sites. At the community level, we explored the effect of bumble bee visitation on the distribution patterns of honey bees on squash flowers. Conclusively, bumble bee visitation caused an increase in elapsed time before flowers were visited again by a honey bee, a behavioral avoidance by a newly-arriving honey bee to select flowers occupied by bumble bees, and a shortened length of time the honey bee takes to examine and collect floral resources. The number of overall bumble bees on squash flowers was the most important factor explaining the difference in the distribution patterns of honey bees at the community level. Furthermore, decline in the number of overall bumble bees on the squash flowers resulted in an increase in the number of overall honey bees. Therefore, our study suggests that bee interactions provide an opportunity to enhance the resilience of ecosystem pollination services against the decline in pollinator diversity.     

Ecological genomics of local adaptation in Cornus florida L. by genotyping by sequencing

Discovering local adaptation, its genetic underpinnings, and environmental drivers is important for conserving forest species. Ecological genomic approaches coupled with next‐generation sequencing are useful means to detect local adaptation and uncover its underlying genetic basis in nonmodel species. We report results from a study on flowering dogwood trees (Cornus florida L.) using genotyping by sequencing (GBS). This species is ecologically important to eastern US forests but is severely threatened by fungal diseases. We analyzed subpopulations in divergent ecological habitats within North Carolina to uncover loci under local selection and associated with environmental–functional traits or disease infection. At this scale, we tested the effect of incorporating additional sequencing before scaling for a broader examination of the entire range. To test for biases of GBS, we sequenced two similarly sampled libraries independently from six populations of three ecological habitats. We obtained environmental–functional traits for each subpopulation to identify associations with genotypes via latent factor mixed modeling (LFMM) and gradient forests analysis. To test whether heterogeneity of abiotic pressures resulted in genetic differentiation indicative of local adaptation, we evaluated F_st per locus while accounting for genetic differentiation between coastal subpopulations and Piedmont‐Mountain subpopulations. Of the 54 candidate loci with sufficient evidence of being under selection among both libraries, 28–39 were Arlequin–BayeScan F_st outliers. For LFMM, 45 candidates were associated with climate (of 54), 30 were associated with soil properties, and four were associated with plant health. Reanalysis of combined libraries showed that 42 candidate loci still showed evidence of being under selection. We conclude environment‐driven selection on specific loci has resulted in local adaptation in response to potassium deficiencies, temperature, precipitation, and (to a marginal extent) disease. High allele turnover along ecological gradients further supports the adaptive significance of loci speculated to be under selection.     

The effects of landscape on bumblebees to ensure crop pollination in the highland agricultural ecosystems in China

Honey bees and wild bees provide critical pollination services to agricultural ecosystems; however, the relative contributions of different bee taxa are not well understood. The natural habitats surrounding farmland support food and nesting resources for wild bees and therefore play an important role in the maintenance of crop pollination. In this study, we selected Cucurbita pepo L. (squash) as a model crop to investigate the relative importance of honey bees and bumblebees in pollinating the crop. Thirteen fields, which were surrounded by a gradient of natural habitat, were investigated on the Yunnan‐Guizhou Plateau in China. We measured the visit densities of honey bees and bumblebees, the number of pollen grains deposited in a single visit by the two bee taxa, as well as the overall pollen grains deposited on stigmas during a flowering day, and then used Bayesian inference to decouple the pollen grains deposited by either the honey bees or the bumblebees. Compared with honey bees, bumblebees deposited a higher number of pollen grains on stigmas in a single visit, but had a lower visit density than honey bees. Meanwhile, the bumblebee visit density increased along the proportion of natural habitat, while the honey bee visit density was not affected by the surrounding natural habitat. Data simulations using Bayesian inference showed that on a flowering day, the number of pollen grains deposited by bumblebees increased with the proportion of natural habitat in the surrounding landscape, but the number of pollen grains deposited by honey bees did not. Moreover, the total numbers of pollen grains deposited by honey bees or bumblebees alone were all below 2000 (the critical level to satisfy the pollination requirement of this crop). Pollen calculations demonstrated that the number of pollen grains deposited by the two bee taxa was greater than 2000 in fields surrounded by more than 13% natural habitat (grasslands and forests). The results revealed that bumblebees ensured C. pepo pollination in combination with honey bees in the highland agricultural ecosystems.     

Orchard pollination in Capitol Reef National Park,Utah, USA. Honey bees or native bees?

Capitol Reef National Park in central Utah, USA surrounds 22 managed fruit orchards started over a century ago by Mormon pioneers. Honey bees are imported for pollination, although the area in which the Park is embedded has over 700 species of native bees, many of which are potential orchard pollinators. We studied the visitation of native bees to apple, pear, apricot, and sweet cherry over 2 years. Thirty species of bees visited the flowers but, except for pear flowers, most were uncommon compared to honey bees. Evidence that honey bees prevented native bees from foraging on orchard crop flowers was equivocal: generally, honey bee and native bee visitation rates to the flowers were not negatively correlated, nor were native bee visitation rates positively correlated with distance of orchards from honey bee hives. Conversely, competition was tentatively suggested by much larger numbers of honey bees than natives on the flowers of apples, apricots and cherry; and by the large increase of native bees on pears, where honey bee numbers were low. At least one-third of the native bee species visiting the flowers are potential pollinators, including cavity-nesting species such as Osmia lignaria propinqua, currently managed for small orchard pollination in the US, plus several fossorial species, including one rosaceous flower specialist (Andrena milwaukiensis). We suggest that gradual withdrawal of honey bees from the Park would help conserve native bee populations without decreasing orchard crop productivity, and would serve as a demonstration of the commercial value of native pollinators.     

Neutral and functionally important genes shed light on phylogeography and the history of high‐altitude colonization in a widespread New World duck

Phylogeographic studies often infer historical demographic processes underlying species distributions based on patterns of neutral genetic variation, but spatial variation in functionally important genes can provide additional insights about biogeographic history allowing for inferences about the potential role of adaptation in geographic range evolution. Integrating data from neutral markers and genes involved in oxygen (O₂)‐transport physiology, we test historical hypotheses about colonization and gene flow across low‐ and high‐altitude regions in the Ruddy Duck (Oxyura jamaicensis), a widely distributed species in the New World. Using multilocus analyses that for the first time include populations from the Colombian Andes, we also examined the hypothesis that Ruddy Duck populations from northern South America are of hybrid origin. We found that neutral and functional genes appear to have moved into the Colombian Andes from both North America and southern South America, and that high‐altitude Colombian populations do not exhibit evidence of adaptation to hypoxia in hemoglobin genes. Therefore, the biogeographic history of Ruddy Ducks is likely more complex than previously inferred. Our new data raise questions about the hypothesis that adaptation via natural selection to high‐altitude conditions through amino acid replacements in the hemoglobin protein allowed Ruddy Ducks to disperse south along the high Andes into southern South America. The existence of shared genetic variation with populations from both North America and southern South America as well as private alleles suggests that the Colombian population of Ruddy Ducks may be of old hybrid origin. This study illustrates the breadth of inferences one can make by combining data from nuclear and functionally important loci in phylogeography, and underscores the importance of complete range‐wide sampling to study species history in complex landscapes.     

Bumble bee abundance and richness improves honey bee pollination behaviour in sweet cherry

Pollination service in agricultural crops increases significantly with pollinator diversity and wild pollinator abundance. Differences in the foraging behaviour of pollinating insects are one of the reasons why pollinator diversity and abundance enhances crop pollination. Here, we focused on the foraging behaviour of honey bees and bumble bees in sweet cherry orchards. In addition, we studied the influence of bee diversity and abundance on the foraging behaviour of honey bees and bumble bees. Honey bees were found to visit fewer flowers than bumble bees. Bumble bees also showed a higher probability of changing trees between rows than honey bees. Both visitation rate and probability of row changes of honey bees increased with bumble bee diversity and with bumble bee abundance. We also found that the probability of row changes of honey bees increased with increasing bumble bee abundance. These effects of bumble bee richness and abundance on the pollination behaviour of honey bees can improve the pollination performance of honey bees in crops that depend on cross pollination. Our results highlight the higher pollination performance of bumble bees and the facilitative effect of wild pollinators to crop pollination.     

Molecular and phylogenetic characterization of honey bee viruses,Nosema microsporidia,protozoan parasites,and parasitic mites in China

China has the largest number of managed honey bee colonies, which produce the highest quantity of honey and royal jelly in the world; however, the presence of honey bee pathogens and parasites has never been rigorously identified in Chinese apiaries. We thus conducted a molecular survey of honey bee RNA viruses, Nosema microsporidia, protozoan parasites, and tracheal mites associated with nonnative Apis mellifera ligustica and native Apis cerana cerana colonies in China. We found the presence of black queen cell virus (BQCV), chronic bee paralysis virus (CBPV), deformed wing virus (DWV), Israeli acute paralysis virus (IAPV), and sacbrood virus (SBV), but not that of acute bee paralysis virus (ABPV) or Kashmir bee virus (KBV). DWV was the most prevalent in the tested samples. Phylogenies of Chinese viral isolates demonstrated that genetically heterogeneous populations of BQCV, CBPV, DWV, and A. cerana‐infecting SBV, and relatively homogenous populations of IAPV and A. meliifera‐infecting new strain of SBV with single origins, are spread in Chinese apiaries. Similar to previous observations in many countries, Nosema ceranae, but not Nosema apis, was prevalent in the tested samples. Crithidia mellificae, but not Apicystis bombi was found in five samples, including one A. c. cerana colony, demonstrating that C. mellificae is capable of infecting multiple honey bee species. Based on kinetoplast‐encoded cytochrome b sequences, the C. mellificae isolate from A. c. cerana represents a novel haplotype with 19 nucleotide differences from the Chinese and Japanese isolates from A. m. ligustica. This suggests that A. c. cerana is the native host for this specific haplotype. The tracheal mite, Acarapis woodi, was detected in one A. m. ligustica colony. Our results demonstrate that honey bee RNA viruses, N. ceranae, C. mellificae, and tracheal mites are present in Chinese apiaries, and some might be originated from native Asian honey bees.     

Genomic and phenotypic differentiation of Arabidopsis thaliana along altitudinal gradients in the North Italian Alps

Altitudinal gradients in mountain regions are short‐range clines of different environmental parameters such as temperature or radiation. We investigated genomic and phenotypic signatures of adaptation to such gradients in five Arabidopsis thaliana populations from the North Italian Alps that originated from 580 to 2350 m altitude by resequencing pools of 19–29 individuals from each population. The sample includes two pairs of low‐ and high‐altitude populations from two different valleys. High‐altitude populations showed a lower nucleotide diversity and negative Tajima's D values and were more closely related to each other than to low‐altitude populations from the same valley. Despite their close geographic proximity, demographic analysis revealed that low‐ and high‐altitude populations split between 260 000 and 15 000 years before present. Single nucleotide polymorphisms whose allele frequencies were highly differentiated between low‐ and high‐altitude populations identified genomic regions of up to 50 kb length where patterns of genetic diversity are consistent with signatures of local selective sweeps. These regions harbour multiple genes involved in stress response. Variation among populations in two putative adaptive phenotypic traits, frost tolerance and response to light/UV stress was not correlated with altitude. Taken together, the spatial distribution of genetic diversity reflects a potentially adaptive differentiation between low‐ and high‐altitude populations, whereas the phenotypic differentiation in the two traits investigated does not. It may resemble an interaction between adaptation to the local microhabitat and demographic history influenced by historical glaciation cycles, recent seed dispersal and genetic drift in local populations.     

Genomic signatures of fine‐scale local selection in Atlantic salmon suggest involvement of sexual maturation,energy homeostasis and immune defence‐related genes

Elucidating the genetic basis of adaptation to the local environment can improve our understanding of how the diversity of life has evolved. In this study, we used a dense SNP array to identify candidate loci potentially underlying fine‐scale local adaptation within a large Atlantic salmon (Salmo salar) population. By combining outlier, gene–environment association and haplotype homozygosity analyses, we identified multiple regions of the genome with strong evidence for diversifying selection. Several of these candidate regions had previously been identified in other studies, demonstrating that the same loci could be adaptively important in Atlantic salmon at subdrainage, regional and continental scales. Notably, we identified signals consistent with local selection around genes associated with variation in sexual maturation, energy homeostasis and immune defence. These included the large‐effect age‐at‐maturity gene vgll3, the known obesity gene mc4r, and major histocompatibility complex II. Most strikingly, we confirmed a genomic region on Ssa09 that was extremely differentiated among subpopulations and that is also a candidate for local selection over the global range of Atlantic salmon. This region colocalized with a haplotype strongly associated with spawning ecotype in sockeye salmon (Oncorhynchus nerka), with circumstantial evidence that the same gene (six6) may be the selective target in both cases. The phenotypic effect of this region in Atlantic salmon remains cryptic, although allelic variation is related to upstream catchment area and covaries with timing of the return spawning migration. Our results further inform management of Atlantic salmon and open multiple avenues for future research.