Two extended haplotype blocks are associated with adaptation to high altitude habitats in East African honey bees

Quantifying Abdominal Coloration of Worker Honey Bees

The main drawback in using coloration to identify honey bee subspecies is the lack of knowledge regarding genetic background, subjectivity of coloration grading, and the effect of the environment. The aim of our study was to evaluate the effect of environmental temperature on the abdominal coloration of honey bee workers and to develop a tool for quantifying abdominal coloration. We obtained four frames of honey bee brood from two colonies and incubated them at two different temperatures (30 and 34 °C). One colony had workers exhibiting yellow marks on the abdomen, while the other did not. We collected hatched workers and photographed abdomens. Images were analyzed using custom-written R script to obtain vectors that summarize the coloration over the abdomen length in a single value-coloration index. We used UMAP to reduce the dimensions of the vectors and to develop a classification procedure with the support vector machine method. We tested the effect of brood origin and temperature on coloration index with ANOVA. UMAP did not distinguish individual abdomens according to experimental group. The trained classifier sufficiently separated abdomens incubated at different temperatures. We improved the performance by preprocessing data with UMAP. The differences among the mean coloration index values were not significant between the gray groups incubated at different temperatures nor between the yellow groups. However, the differences between the gray and yellow groups were significant, permitting options for application of our tool and the newly developed coloration index. Our results indicate that the environmental temperature in the selected range during development does not seem to impact honey bee coloration significantly. The developed color-recording protocol and statistical analysis provide useful tools for quantifying abdominal coloration in honey bees.

Chromosomal Inversions and the Demography of Speciation in Drosophila montana and Drosophila flavomontana

Chromosomal inversions may play a central role in speciation given their ability to locally reduce recombination and therefore genetic exchange between diverging populations. We analyzed long- and short-read whole-genome data from sympatric and allopatric populations of 2 Drosophila virilis group species, Drosophila montana and Drosophila flavomontana, to understand if inversions have contributed to their divergence. We identified 3 large alternatively fixed inversions on the X chromosome and one on each of the autosomes 4 and 5. A comparison of demographic models estimated for inverted and noninverted (colinear) chromosomal regions suggests that these inversions arose before the time of the species split. We detected a low rate of interspecific gene flow (introgression) from D. montana to D. flavomontana, which was further reduced inside inversions and was lower in allopatric than in sympatric populations. Together, these results suggest that the inversions were already present in the common ancestral population and that gene exchange between the sister taxa was reduced within inversions both before and after the onset of species divergence. Such ancestrally polymorphic inversions may foster speciation by allowing the accumulation of genetic divergence in loci involved in adaptation and reproductive isolation inside inversions early in the speciation process, while gene exchange at colinear regions continues until the evolving reproductive barriers complete speciation. The overlapping X inversions are particularly good candidates for driving the speciation process of D. montana and D. flavomontana, since they harbor strong genetic incompatibilities that were detected in a recent study of experimental introgression.

Pan-genome analysis highlights the role of structural variation in the evolution and environmental adaptation of Asian honeybees

The Asian honeybee, Apis cerana, is an ecologically and economically important pollinator. Mapping its genetic variation is key to understanding population-level health, histories and potential capacities to respond to environmental changes. However, most efforts to date were focused on single nucleotide polymorphisms (SNPs) based on a single reference genome, thereby ignoring larger scale genomic variation. We employed long-read sequencing technologies to generate a chromosome-scale reference genome for the ancestral group of A. cerana. Integrating this with 525 resequencing data sets, we constructed the first pan-genome of A. cerana, encompassing almost the entire gene content. We found that 31.32% of genes in the pan-genome were variably present across populations, providing a broad gene pool for environmental adaptation. We identified and characterized structural variations (SVs) and found that they were not closely linked with SNP distributions; however, the formation of SVs was closely associated with transposable elements. Furthermore, phylogenetic analysis using SVs revealed a novel A. cerana ecological group not recoverable from the SNP data. Performing environmental association analysis identified a total of 44 SVs likely to be associated with environmental adaptation. Verification and analysis of one of these, a 330 bp deletion in the Atpalpha gene, indicated that this SV may promote the cold adaptation of A. cerana by altering gene expression. Taken together, our study demonstrates the feasibility and utility of applying pan-genome approaches to map and explore genetic feature variations of honeybee populations, and in particular to examine the role of SVs in the evolution and environmental adaptation of A. cerana.

Ecology and Management of African Honey Bees (Apis mellifera L.)

In Africa, humans evolved as honey hunters of honey bee subspecies adapted to diverse geographical regions. Beekeeping today is practiced much as it was when Africans moved from honey hunting to beekeeping nearly 5,000 years ago, with beekeepers relying on seasonally available wild bees. Research suggests that populations are resilient, able to resist diseases and novel parasites. Distinct biomes, as well as environmental pressures, shaped the behavior and biology of these bees and in turn influenced how indigenous beekeeping developed. It appears that passive beekeeping practices that enabled free-living populations contributed to the overall resilience and health of the bee. There is clearly a need for research aimed at a deeper understanding of bee biology and the ecosystems from which they benefit and on which humans depend, as well as a growing realization that the management of these bees requires an indigenous approach that reflects a broader knowledge base and the economics of local communities and markets.

Chromosome study of the Hymenoptera (Insecta): from cytogenetics to cytogenomics

A brief overview of the current stage of the chromosome study of the insect order Hymenoptera is given. It is demonstrated that, in addition to routine staining and other traditional techniques of chromosome research, karyotypes of an increasing number of hymenopterans are being studied using molecular methods, e.g., staining with base-specific fluorochromes and fluorescence in situ hybridization (FISH), including microdissection and chromosome painting. Due to the advent of whole genome sequencing and other molecular techniques, together with the "big data" approach to the chromosomal data, the current stage of the chromosome research on Hymenoptera represents a transition from Hymenoptera cytogenetics to cytogenomics.

Genetic differentiation and local adaptation of the Japanese honeybee, Apis cerana japonica

We examine the population genetic structure and divergence among the regional populations of the Japanese honeybee, Apis cerana japonica, by re-sequencing the genomes of 105 individuals from the three main Japanese islands with diverse climates. The genetic structure results indicated that these individuals are distinct from the mainland Chinese A. cerana samples. Furthermore, population structure analyses have identified three genetically distinct geographic regions in Japan: Northern (Tohoku-Kanto-Chubu districts), Central (Chugoku district), and Southern (Kyushu district). In some districts, "possible non-native" individuals, likely introduced from other regions in recent years, were discovered. Then, genome-wide scans were conducted to detect candidate genes for adaptation by two different approaches. We performed a population branch statistics (PBS) analysis to identify candidate genes for population-specific divergence. A latent factor mixed model (LFMM) was used to identify genes associated with climatic variables along a geographic gradient. The PBSmax analysis identified 25 candidate genes for population-specific divergence whereas the LFMM analysis identified 73 candidate genes for adaptation to climatic variables along a geographic gradient. However, no common genes were identified by both methods.

The Genomic Basis of Adaptation to High Elevations in Africanized Honey Bees

A range of different genetic architectures underpin local adaptation in nature. Honey bees (Apis mellifera) in the Eastern African Mountains harbor high frequencies of two chromosomal inversions that likely govern adaptation to this high-elevation habitat. In the Americas, honey bees are hybrids of European and African ancestries and adaptation to latitudinal variation in climate correlates with the proportion of these ancestries across the genome. It is unknown which, if either, of these forms of genetic variation governs adaptation in honey bees living at high elevations in the Americas. Here, we performed whole-genome sequencing of 29 honey bees from both high- and low-elevation populations in Colombia. Analysis of genetic ancestry indicated that both populations were predominantly of African ancestry, but the East African inversions were not detected. However, individuals in the higher elevation population had significantly higher proportions of European ancestry, likely reflecting local adaptation. Several genomic regions exhibited particularly high differentiation between highland and lowland bees, containing candidate loci for local adaptation. Genes that were highly differentiated between highland and lowland populations were enriched for functions related to reproduction and sperm competition. Furthermore, variation in levels of European ancestry across the genome was correlated between populations of honey bees in the highland population and populations at higher latitudes in South America. The results are consistent with the hypothesis that adaptation to both latitude and elevation in these hybrid honey bees are mediated by variation in ancestry at many loci across the genome.

Divergent contributions of coding and noncoding sequences to initial high-altitude adaptation in passerine birds endemic to the Qinghai-Tibet Plateau

Early events in the evolution of an ancestral lineage can shape the adaptive patterns of descendant species, but the evolutionary mechanisms driving initial adaptation from an ancestor remain largely unexplored. High-altitude adaptations have been extensively explored from the viewpoint of protein-coding genes; however, the contribution of noncoding regions remains relatively neglected. Here, we integrate genomic and transcriptomic data to investigate adaptive evolution in the ancestor of three high-altitude snowfinch species endemic to the Qinghai-Tibet Plateau. Our genome-wide scan for adaptation in the snowfinch ancestor identifies strong adaptation signals in functions of development and metabolism for the coding genes, but in functions of the nervous system development for noncoding regions. This pattern is exclusive to the snowfinch ancestor compared to a control ancestral lineage subject to weak selection. Changes in noncoding regions in the snowfinch ancestor, especially those nearest to coding genes, may be disproportionately associated with the differential expression of genes in the brain tissue compared to other tissues. Extensive gene expression in the brain tissue can be further altered via genetic regulatory networks of transcription factors harbouring potential accelerated regulatory regions (e.g., the development-related transcription factor YEATS4). Altogether, our study provides new evidence concerning how coding and noncoding sequences work through decoupled pathways in initial adaptation to the selective pressure of high-altitude environments. The analysis highlights the idea that noncoding sequences may be promising elements in facilitating the rapid evolution and adaptation to high altitudes.

Impacts of oxygen deficiency on embryo life-history traits of migratory locust Locusta migratoria from low and high altitudes

Hypoxia challenges aerobic organisms in numerous environments, and hypoxic conditions may become more severe under future climate-change scenarios. The impact of hypoxia on the development of terrestrial insect embryos is not well understood. Here, to address this gap, embryonic life-history traits of migratory locust Locusta migratoria from low-altitude and high-altitude regions were compared under 2 oxygen levels: normoxia (i.e., 21 kPa oxygen partial pressure and mild hypoxia (i.e., 10 kPa oxygen partial pressure). Our results demonstrated that, whether reared under normoxia or mild hypoxia, L. migratoria from high-altitude populations had longer developmental times, reduced weight, and lower mean relative growth rate as compared with those from low-altitude populations. When transferred from normoxia to mild hypoxia, nearly all the tested life-history traits presented significant negative changes in the low-altitude populations, but not in the high-altitude populations. The factor 'strain' alone explained 18.26%-54.59% of the total variation for traits, suggesting that the phenotypic differences between L. migratoria populations from the 2 altitudes could be driven by genetic variation. Significant genetic correlations were found between life-history traits, and most of these showed differentiation between the 2 altitudinal gradients. G-matrix comparisons showed significant structural differences between L. migratoria from the 2 regions, as well as several negative covariances (i.e., trade-offs) between traits in the low-altitude populations. Overall, our study provides clear evidence that evolutionary divergence of embryonic traits between L. migratoria populations from different altitudes has occurred.

The impact of geography and climate on the population structure and local adaptation in a wild bee

Deciphering processes that contribute to genetic differentiation and divergent selection of natural populations is useful for evaluating the adaptive potential and resilience of organisms faced with various anthropogenic stressors. Insect pollinator species, including wild bees, provide critical ecosystem services but are highly susceptible to biodiversity declines. Here, we use population genomics to infer the genetic structure and test for evidence of local adaptation in an economically important native pollinator, the small carpenter bee (Ceratina calcarata). Using genome-wide SNP data (n = 8302), collected from specimens across the species' entire distribution, we evaluated population differentiation and genetic diversity and identified putative signatures of selection in the context of geographic and environmental variation. Results of the analyses of principal component and Bayesian clustering were concordant with the presence of two to three genetic clusters, associated with landscape features and inferred phylogeography of the species. All populations examined in our study demonstrated a heterozygote deficit, along with significant levels of inbreeding. We identified 250 robust outlier SNPs, corresponding to 85 annotated genes with known functional relevance to thermoregulation, photoperiod, and responses to various abiotic and biotic stressors. Taken together, these data provide evidence for local adaptation in a wild bee and highlight genetic responses of native pollinators to landscape and climate features.

AmelHap: Leveraging drone whole-genome sequence data to create a honey bee HapMap

Honey bee, Apis mellifera, drones are typically haploid, developing from an unfertilized egg, inheriting only their queen's alleles and none from the many drones she mated with. Thus the ordered combination or 'phase' of alleles is known, making drones a valuable haplotype resource. We collated whole-genome sequence data for 1,407 drones, including 45 newly sequenced Scottish drones, collectively representing 19 countries, 8 subspecies and various hybrids. Following alignment to Amel_HAv3.1, variant calling and quality filtering, we retained 17.4 M high quality variants across 1,328 samples with a genotyping rate of 98.7%. We demonstrate the utility of this haplotype resource, AmelHap, for genotype imputation, returning >95% concordance when up to 61% of data is missing in haploids and up to 12% of data is missing in diploids. AmelHap will serve as a useful resource for the community for imputation from low-depth sequencing or SNP chip data, accurate phasing of diploids for association studies, and as a comprehensive reference panel for population genetic and evolutionary analyses.

A key gene for the climatic adaptation of Apis cerana populations in China according to selective sweep analysis

BACKGROUND

Apis cerana is widely distributed in China and, prior to the introduction of western honeybees, was the only bee species kept in China. During the long-term natural evolutionary process, many unique phenotypic variations have occurred among A. cerana populations in different geographical regions under varied climates. Understanding the molecular genetic basis and the effects of climate change on the adaptive evolution of A. cerana can promote A. cerana conservation in face of climate change and allow for the effective utilization of its genetic resources.

RESULT

To investigate the genetic basis of phenotypic variations and the impact of climate change on adaptive evolution, A. cerana workers from 100 colonies located at similar geographical latitudes or longitudes were analyzed. Our results revealed an important relationship between climate types and the genetic variation of A. cerana in China, and a greater influence of latitude compared with longitude was observed. Upon selection and morphometry analyses combination for populations under different climate types, we identified a key gene RAPTOR, which was deeply involved in developmental processes and influenced the body size.

CONCLUSION

The selection of RAPTOR at the genomic level during adaptive evolution could allow A. cerana to actively regulate its metabolism, thereby fine-tuning body sizes in response to harsh conditions caused by climate change, such as food shortages and extreme temperatures, which may partially elucidate the size differences of A. cerana populations. This study provides crucial support for the molecular genetic basis of the expansion and evolution of naturally distributed honeybee populations.

On the origin and structure of haplotype blocks

The term "haplotype block" is commonly used in the developing field of haplotype-based inference methods. We argue that the term should be defined based on the structure of the Ancestral Recombination Graph (ARG), which contains complete information on the ancestry of a sample. We use simulated examples to demonstrate key features of the relationship between haplotype blocks and ancestral structure, emphasizing the stochasticity of the processes that generate them. Even the simplest cases of neutrality or of a "hard" selective sweep produce a rich structure, often missed by commonly used statistics. We highlight a number of novel methods for inferring haplotype structure, based on the full ARG, or on a sequence of trees, and illustrate how they can be used to define haplotype blocks using an empirical data set. While the advent of new, computationally efficient methods makes it possible to apply these concepts broadly, they (and additional new methods) could benefit from adding features to explore haplotype blocks, as we define them. Understanding and applying the concept of the haplotype block will be essential to fully exploit long and linked-read sequencing technologies.

Population Genomics for Insect Conservation

Insects constitute vital components of ecosystems. There is alarming evidence for global declines in insect species diversity, abundance, and biomass caused by anthropogenic drivers such as habitat degradation or loss, agricultural practices, climate change, and environmental pollution. This raises important concerns about human food security and ecosystem functionality and calls for more research to assess insect population trends and identify threatened species and the causes of declines to inform conservation strategies. Analysis of genetic diversity is a powerful tool to address these goals, but so far animal conservation genetics research has focused strongly on endangered vertebrates, devoting less attention to invertebrates, such as insects, that constitute most biodiversity. Insects' shorter generation times and larger population sizes likely necessitate different analytical methods and management strategies. The availability of high-quality reference genome assemblies enables population genomics to address several key issues. These include precise inference of past demographic fluctuations and recent declines, measurement of genetic load levels, delineation of evolutionarily significant units and cryptic species, and analysis of genetic adaptation to stressors. This enables identification of populations that are particularly vulnerable to future threats, considering their potential to adapt and evolve. We review the application of population genomics to insect conservation and the outlook for averting insect declines.

Reconstructing queen genotypes by pool sequencing colonies in eusocial insects: Statistical Methods and their application to honeybee

Eusocial insects are crucial to many ecosystems, and particularly the honeybee (Apis mellifera). One approach to facilitate their study in molecular genetics, is to consider whole-colony genotyping by combining DNA of multiple individuals in a single pool sequencing experiment. Cheap and fast, this technique comes with the drawback of producing data requiring dedicated methods to be fully exploited. Despite this limitation, pool sequencing data have been shown to be informative and cost-effective when working on random mating populations. Here, we present new statistical methods for exploiting pool sequencing of eusocial colonies in order to reconstruct the genotypes of the queen of such colony. This leverages the possibility to monitor genetic diversity, perform genomic-based studies or implement selective breeding. Using simulations and honeybee real data, we show that the new methods allow for a fast and accurate estimation of the queen's genetic ancestry, with correlations of about 0.9 to that obtained from individual genotyping. Also, it allows for an accurate reconstruction of the queen genotypes, with about 2% genotyping error. We further validate these inferences using experimental data on colonies with both pool sequencing and individual genotyping of drones. In brief, in this study we present statistical models to accurately estimate the genetic ancestry and reconstruct the genotypes of the queen from pool sequencing data from workers of an eusocial colony. Such information allows to exploit pool sequencing for traditional population genetics analyses, association studies and for selective breeding. While validated in Apis mellifera, these methods are applicable to other eusocial hymenopterans.

Inversions and parallel evolution

Local adaptation leads to differences between populations within a species. In many systems, similar environmental contrasts occur repeatedly, sometimes driving parallel phenotypic evolution. Understanding the genomic basis of local adaptation and parallel evolution is a major goal of evolutionary genomics. It is now known that by preventing the break-up of favourable combinations of alleles across multiple loci, genetic architectures that reduce recombination, like chromosomal inversions, can make an important contribution to local adaptation. However, little is known about whether inversions also contribute disproportionately to parallel evolution. Our aim here is to highlight this knowledge gap, to showcase existing studies, and to illustrate the differences between genomic architectures with and without inversions using simple models. We predict that by generating stronger effective selection, inversions can sometimes speed up the parallel adaptive process or enable parallel adaptation where it would be impossible otherwise, but this is highly dependent on the spatial setting. We highlight that further empirical work is needed, in particular to cover a broader taxonomic range and to understand the relative importance of inversions compared to genomic regions without inversions.

This article is part of the theme issue ‘Genomic architecture of supergenes: causes and evolutionary consequences’.

Octopamine drives honeybee thermogenesis

In times of environmental change species have two options to survive: they either relocate to a new habitat or they adapt to the altered environment. Adaptation requires physiological plasticity and provides a selection benefit. In this regard, the Western honeybee (Apis mellifera) protrudes with its thermoregulatory capabilities, which enables a nearly worldwide distribution. Especially in the cold, shivering thermogenesis enables foraging as well as proper brood development and thus survival. In this study, we present octopamine signaling as a neurochemical prerequisite for honeybee thermogenesis: we were able to induce hypothermia by depleting octopamine in the flight muscles. Additionally, we could restore the ability to increase body temperature by administering octopamine. Thus, we conclude that octopamine signaling in the flight muscles is necessary for thermogenesis. Moreover, we show that these effects are mediated by β octopamine receptors. The significance of our results is highlighted by the fact the respective receptor genes underlie enormous selective pressure due to adaptation to cold climates. Finally, octopamine signaling in the service of thermogenesis might be a key strategy to survive in a changing environment.

Further Evidence of Population Admixture in the Serbian Honey Bee Population

Simple Summary

The western honey bee is one of the most ecologically and economically important pollinator species. Due to human interference, it faces serious challenges, not only in number decline and habitat loss, but also in natural subspecies diversity and distribution. The conservation of genetic diversity and perseverance of locally adapted populations and subspecies becomes a crucial task in the face of rapid environmental changes. In order to further assess present genetic variability in Serbian honey bee populations, we analyzed 14 microsatellite loci and then compared nine of them with previously published data. Our results suggest that Serbia now harbors a distinct, relatively homogenous honey bee population, although some local differences are still preserved.

Abstract

Socioeconomic interests and beekeeper preferences have often taken precedence over the conservation of locally native honey bee subspecies, leading to the predominance of admixture populations in human-dominated areas. To assess the genetic diversity of contemporary managed Serbian honey bee colonies, we used 14 microsatellite loci and analyzed 237 worker bees from 46 apiaries in eight localities of northern and southern Serbia. Furthermore, we compared data for nine microsatellite loci with 338 individuals from Italy, Hungary, Poland, and Spain. The standard parameters of genetic diversity in Serbian honey bee populations were in line with other analyses, although somewhat smaller. STRUCTURE analysis showed the existence of two equally distributed genetic clusters and Analysis of molecular variances could not confirm the presence of a geographically discrete population but showed local differences. Discriminant analysis of principal components showed overlapping of worker bees from different parts of Serbia. Clear genetic differentiation can be observed when comparing all populations between geographical regions and their corresponding subspecies. The absence of the A. m. macedonica subspecies from its historical distribution range in southern Serbia as well as the lack of distinctive geographical groups suggest that selective breeding, queen import, and migratory beekeeping practices strongly influence the genetic structure and diversity of honey bees, leading to the genetic uniformization and creation of the admixture population.

Attack of the dark clones the genetics of reproductive and color traits of South African honey bees (Apis mellifera spp.)

The traits of two subspecies of western honey bees, Apis mellifera scutellata and A.m. capensis, endemic to the Republic of South Africa (RSA), are of biological and commercial relevance. Nevertheless, the genetic basis of important phenotypes found in these subspecies remains poorly understood. We performed a genome wide association study on three traits of biological relevance in 234 A.m. capensis, 73 A.m. scutellata and 158 hybrid individuals. Thirteen markers were significantly associated to at least one trait (P ≤ 4.28 × 10⁻⁶): one for ovariole number, four for scutellar plate and eight for tergite color. We discovered two possible causative variants associated to the respective phenotypes: a deletion in GB46429 or Ebony (NC_007070.3:g.14101325G>del) (R69Efs*85) and a nonsense on GB54634 (NC_007076.3:g.4492792A>G;p.Tyr128*) causing a premature stop, substantially shortening the predicted protein. The mutant genotypes are significantly associated to phenotypes in A.m. capensis. Loss-of-function of Ebony can cause accumulation of circulating dopamine, and increased dopamine levels correlate to ovary development in queenless workers and pheromone production. Allelic association (P = 1.824 x 10⁻⁵) of NC_007076.3:g.4492792A>G;p.Tyr128* to ovariole number warrants further investigation into function and expression of the GB54634 gene. Our results highlight genetic components of relevant production/conservation behavioral phenotypes in honey bees.

Stable inversion clines in a grasshopper species group despite complex geographical history

Chromosomal inversions are known to play roles in adaptation and differentiation in many species. They involve clusters of correlated genes (i.e., loci in linkage disequilibrium, LD) possibly associated with environmental variables. The grasshopper "species complex" Trimerotropis pallidipennis comprises several genetic lineages distributed from North to South America in arid and semi-arid high-altitude environments. The southernmost lineage, Trimerotropis sp., segregates for four to seven putative inversions that display clinal variation, possibly through adaptation to temperate environments. We analysed chromosomal, mitochondrial and genome-wide single nucleotide polymorphism data in 19 Trimerotropis sp. populations mainly distributed along two altitudinal gradients (MS and Ju). Populations across Argentina comprise two main chromosomally and genetically differentiated lineages: one distributed across the southernmost border of the "Andes Centrales," adding evidence for a differentiation hotspot in this area; and the other widely distributed in Argentina. Within the latter, network analytical approaches to LD found three clusters of correlated loci (LD-clusters), with inversion karyotypes explaining >79% of the genetic variation. Outlier loci associated with environmental variables mapped to two of these LD-clusters. Furthermore, despite the complex geographical history indicated by population genetic analyses, the clines in inversion karyotypes have remained stable for more than 20 generations, implicating their role in adaptation and differentiation within this lineage. We hypothesize that these clines could be the consequence of a coupling between extrinsic postzygotic barriers and spatially varying selection along environmental gradients resulting in a hybrid zone. These results provide a framework for future investigations about candidate genes implicated in rapid adaptation to new environments.

Thrice out of Asia and the adaptive radiation of the western honey bee

The origin of the western honey bee Apis mellifera has been intensely debated. Addressing this knowledge gap is essential for understanding the evolution and genetics of one of the world’s most important pollinators. By analyzing 251 genomes from 18 native subspecies, we found support for an Asian origin of honey bees with at least three expansions leading to African and European lineages. The adaptive radiation of honey bees involved selection on a few genomic “hotspots.” We found 145 genes with independent signatures of selection across all bee lineages, and these genes were highly associated with worker traits. Our results indicate that a core set of genes associated with worker and colony traits facilitated the adaptive radiation of honey bees across their vast distribution.

Conservation insights from wild bee genetic studies: Geographic differences, susceptibility to inbreeding, and signs of local adaptation

Conserving bees are critical both ecologically and economically. Genetic tools are valuable for monitoring these vital pollinators since tracking these small, fast-flying insects by traditional means is difficult. By surveying the current state of the literature, this review discusses how recent advances in landscape genetic and genomic research are elucidating how wild bees respond to anthropogenic threats. Current literature suggests that there may be geographic differences in the vulnerability of bee species to landscape changes. Populations of temperate bee species are becoming more isolated and more genetically depauperate as their landscape becomes more fragmented, but tropical bee species appear unaffected. These differences may be an artifact of historical differences in land-use, or it suggests that different management plans are needed for temperate and tropical bee species. Encouragingly, genetic studies on invasive bee species indicate that low levels of genetic diversity may not lead to rapid extinction in bees as once predicted. Additionally, next-generation sequencing has given researchers the power to identify potential genes under selection, which are likely critical to species' survival in their rapidly changing environment. While genetic studies provide insights into wild bee biology, more studies focusing on a greater phylogenetic and life-history breadth of species are needed. Therefore, caution should be taken when making broad conservation decisions based on the currently few species examined.

Disentangling Ethiopian Honey Bee (Apis mellifera) Populations Based on Standard Morphometric and Genetic Analyses

Simple Summary

We conducted this population study of Ethiopian honey bees, using morphometric and genetic methods, to decipher their controversial classification. These honey bees are highly diverse and showed differentiation based on size and genetic information according to prevailing agro-ecological conditions, demonstrating morphological and molecular signatures of local adaptation. The results of both morphometric and genetic analyses suggest that Ethiopian honey bees differ from populations in the neighboring geographic regions and are characterized by extensive gene flow within the country, enhanced by honey bee colony trade. Consequently, future research that includes studying traits of vitality, behavior and colony performance of honey bees in remaining pocket areas of highland agro-ecological zones could contribute to the development of appropriate conservation management.

Abstract

The diversity and local differentiation of honey bees are subjects of broad general interest. In particular, the classification of Ethiopian honey bees has been a subject of debate for decades. Here, we conducted an integrated analysis based on classical morphometrics and a putative nuclear marker (denoted r7-frag) for elevational adaptation to classify and characterize these honey bees. Therefore, 660 worker bees were collected out of 66 colonies from highland, midland and lowland agro-ecological zones (AEZs) and were analyzed in reference to populations from neighboring countries. Multivariate morphometric analyses show that our Ethiopian samples are separate from Apis mellifera scutellata, A. m. jemenitica, A. m. litorea and A. m. monticola, but are closely related to A. m. simensis reference. Linear discriminant analysis showed differentiation according to AEZs in the form of highland, midland and lowland ecotypes. Moreover, size was positively correlated with elevation. Similarly, our Ethiopian samples were differentiated from A. m. monticola and A. m. scutellata based on r7-frag. There was a low tendency towards genetic differentiation between the Ethiopian samples, likely impacted by increased gene flow. However, the differentiation slightly increased with increasing elevational differences, demonstrated by the highland bees that showed higher differentiation from the lowland bees (F_ST = 0.024) compared to the midland bees (F_ST = 0.015). An allelic length polymorphism was detected (denoted as d) within r7-frag, showing a patterned distribution strongly associated with AEZ (X² = 11.84, p < 0.01) and found predominantly in highland and midland bees of some pocket areas. In conclusion, the Ethiopian honey bees represented in this study are characterized by high gene flow that suppresses differentiation.

Genetic markers for the resistance of honey bee to Varroa destructor

In the mid-20th century, the first case of infection of European bees Apis mellifera L. with the ectoparasite mite Varroa destructor was recorded. The original host of this mite is the Asian bee Apis cerana. The mite V. destructor was widespread throughout Europe, North and South America, and Australia remained the only continent free from this parasite. Without acaricide treatment any honeybee colony dies within 1–4 years. The use of synthetic acaricides has not justified itself – they make beekeeping products unsuitable and mites develop resistance to them, which forces the use of even greater concentrations that can be toxic to the bees. Therefore, the only safe measure to combat the mite is the use of biological control methods. One of these methods is the selection of bee colonies with natural mite resistance. In this article we summarize publications devoted to the search for genetic markers associated with resistance to V. destructor. The first part discusses the basic mechanisms of bee resistance (Varroa sensitive hygienic behavior and grooming) and methods for their assessment. The second part focuses on research aimed at searching for loci and candidate genes associated with resistance to varroosis by mapping quantitative traits loci and genome-wide association studies. The third part summarizes studies of the transcriptome profile of Varroa resistant bees. The last part discusses the most likely candidate genes – potential markers for breeding Varroa resistant bees. Resistance to the mite is manifested in a variety of phenotypes and is under polygenic control. The establishing of gene pathways involved in resistance to Varroa will help create a methodological basis for the selection of Varroa resistant honeybee colonies.

Gene reuse facilitates rapid radiation and independent adaptation to diverse habitats in the Asian honeybee

Animals with recent shared ancestry frequently adapt in parallel to new but similar habitats, a process often underlined by repeated selection of the same genes. Yet, in contrast, few examples have demonstrated the significance of gene reuse in colonization of multiple disparate habitats. By analyzing 343 genomes of the widespread Asian honeybee, Apis cerana, we showed that multiple peripheral subspecies radiated from a central ancestral population and adapted independently to diverse habitats. We found strong evidence of gene reuse in the Leucokinin receptor (Lkr), which was repeatedly selected in almost all peripheral subspecies. Differential expression and RNA interference knockdown revealed the role of Lkr in influencing foraging labor division, suggesting that Lkr facilitates collective tendency for pollen/nectar collection as an adaptation to floral changes. Our results suggest that honeybees may accommodate diverse floral shifts during rapid radiation through fine-tuning individual foraging tendency, a seemingly complex process accomplished by gene reuse.

Extensive haplotypes are associated with population differentiation and environmental adaptability in Upland cotton (Gossypium hirsutum)

Three extensive eco-haplotypes associated with population differentiation and environmental adaptability in Upland cotton were identified, with A06_85658585, A08_43734499 and A06_113104285 considered the eco-loci for environmental adaptability. Population divergence is suggested to be the primary force driving the evolution of environmental adaptability in various species. Chromosome inversion increases reproductive isolation between subspecies and accelerates population divergence to adapt to new environments. Although modern cultivated Upland cotton (Gossypium hirsutum L.) has spread worldwide, the noticeable phenotypic differences still existed among cultivars grown in different areas. In recent years, the long-distance migration of cotton cultivation areas throughout China has demanded that breeders better understand the genetic basis of environmental adaptability in Upland cotton. Here, we integrated the genotypes of 419 diverse accessions, long-term environment-associated variables (EAVs) and environment-associated traits (EATs) to evaluate subgroup differentiation and identify adaptive loci in Upland cotton. Two highly divergent genomic regions were found on chromosomes A06 and A08, which likely caused by extensive chromosome inversions. The subgroups could be geographically classified based on distinct haplotypes in the divergent regions. A genome-wide association study (GWAS) also confirmed that loci located in these regions were significantly associated with environmental adaptability in Upland cotton. Our study first revealed the cause of population divergence in Upland cotton, as well as the consequences of variation in its environmental adaptability. These findings provide new insights into the genetic basis of environmental adaptability in Upland cotton, which could accelerate the development of molecular markers for adaptation to climate change in future cotton breeding.

Characterization of the Kenyan Honey Bee (Apis mellifera) Gut Microbiota: A First Look at Tropical and Sub-Saharan African Bee Associated Microbiomes

Gut microbiota plays important roles in many physiological processes of the host including digestion, protection, detoxification, and development of immune responses. The honey bee (Apis mellifera) has emerged as model for gut-microbiota host interaction studies due to its gut microbiota being highly conserved and having a simple composition. A key gap in this model is understanding how the microbiome differs regionally, including sampling from the tropics and in particular from Africa. The African region is important from the perspective of the native diversity of the bees, and differences in landscape and bee management. Here, we characterized the honey bee gut microbiota in sub-Saharan Africa using 16S rRNA amplicon sequencing. We confirm the presence of the core gut microbiota members and highlight different compositions of these communities across regions. We found that bees from the coastal regions harbor a higher relative abundance and diversity on core members. Additionally, we showed that Gilliamella, Snodgrassella, and Frischella dominate in all locations, and that altitude and humidity affect Gilliamella abundance. In contrast, we found that Lactobacillus was less common compared temperate regions of the world. This study is a first comprehensive characterization of the gut microbiota of honey bees from sub-Saharan Africa and underscores the need to study microbiome diversity in other indigenous bee species and regions.

Genome-wide patterns of differentiation within and among U.S. commercial honey bee stocks

BACKGROUND

The population genetics of U.S. honey bee stocks remain poorly characterized despite the agricultural importance of Apis mellifera as the major crop pollinator. Commercial and research-based breeding programs have made significant improvements of favorable genetic traits (e.g. production and disease resistance). The variety of bees produced by artificial selection provides an opportunity to characterize the genetic diversity and regions of the genome undergoing selection in commonly managed stocks.

RESULTS

Pooled sequencing of eight honey bee stocks found strong genetic similarity among six of the stocks. Two stocks, Pol-line and Hilo, showed significant differentiation likely due to their intense and largely closed breeding for resistance to the parasitic Varroa mite. Few variants were identified as being specific to any one stock, indicating potential admixture among the sequenced stocks. Juxtaposing the underlying genetic variation of stocks selected for disease- and parasite-resistance behavior, we identified genes and candidate regions putatively associated with resistance regulated by hygienic behavior.

CONCLUSION

This study provides important insights into the distinct genetic characteristics and population diversity of honey bee stocks used in the United States, and provides further evidence of high levels of admixture in commercially managed honey bee stocks. Furthermore, breeding efforts to enhance parasite resistance in honey bees may have created unique genetic profiles. Genomic regions of interest have been highlighted for potential future work related to developing genetic markers for selection of disease and parasite resistance traits. Due to the vast genomic similarities found among stocks in general, our findings suggest that additional data regarding gene expression, epigenetic and regulatory information are needed to more fully determine how stock phenotypic diversity is regulated.

Eusociality influences the strength of negative selection on insect genomes

While much of the focus of sociobiology concerns identifying genomic changes that influence social behaviour, we know little about the consequences of social behaviour on genome evolution. It has been hypothesized that social evolution can influence the strength of negative selection via two mechanisms. First, division of labour can influence the efficiency of negative selection in a caste-specific manner; indirect negative selection on worker traits is theoretically expected to be weaker than direct selection on queen traits. Second, increasing social complexity is expected to lead to relaxed negative selection because of its influence on effective population size. We tested these two hypotheses by estimating the strength of negative selection in honeybees, bumblebees, paper wasps, fire ants and six other insects that span the range of social complexity. We found no consistent evidence that negative selection was significantly stronger on queen-biased genes relative to worker-biased genes. However, we found strong evidence that increased social complexity reduced the efficiency of negative selection. Our study clearly illustrates how changes in behaviour can influence patterns of genome evolution by modulating the strength of natural selection.

Footprints of local adaptation span hundreds of linked genes in the Atlantic silverside genome

The study of local adaptation in the presence of ongoing gene flow is the study of natural selection in action, revealing the functional genetic diversity most relevant to contemporary pressures. In addition to individual genes, genome‐wide architecture can itself evolve to enable adaptation. Distributed across a steep thermal gradient along the east coast of North America, Atlantic silversides (Menidia menidia) exhibit an extraordinary degree of local adaptation in a suite of traits, and the capacity for rapid adaptation from standing genetic variation, but we know little about the patterns of genomic variation across the species range that enable this remarkable adaptability. Here, we use low‐coverage, whole‐transcriptome sequencing of Atlantic silversides sampled along an environmental cline to show marked signatures of divergent selection across a gradient of neutral differentiation. Atlantic silversides sampled across 1371 km of the southern section of its distribution have very low genome‐wide differentiation (median F_ST = 0.006 across 1.9 million variants), consistent with historical connectivity and observations of recent migrants. Yet almost 14,000 single nucleotide polymorphisms (SNPs) are nearly fixed (F_ST > 0.95) for alternate alleles. Highly differentiated SNPs cluster into four tight linkage disequilibrium (LD) blocks that span hundreds of genes and several megabases. Variants in these LD blocks are disproportionately nonsynonymous and concentrated in genes enriched for multiple functions related to known adaptations in silversides, including variation in lipid storage, metabolic rate, and spawning behavior. Elevated levels of absolute divergence and demographic modeling suggest selection maintaining divergence across these blocks under gene flow. These findings represent an extreme case of heterogeneity in levels of differentiation across the genome, and highlight how gene flow shapes genomic architecture in continuous populations. Locally adapted alleles may be common features of populations distributed along environmental gradients, and will likely be key to conserving variation to enable future responses to environmental change.

Consequences of Single-Locus and Tightly Linked Genomic Architectures for Evolutionary Responses to Environmental Change

Genetic and genomic architectures of traits under selection are key factors influencing evolutionary responses. Yet, knowledge of their impacts has been limited by a widespread assumption that most traits are controlled by unlinked polygenic architectures. Recent advances in genome sequencing and eco-evolutionary modeling are unlocking the potential for integrating genomic information into predictions of population responses to environmental change. Using eco-evolutionary simulations, we demonstrate that hypothetical single-locus control of a life history trait produces highly variable and unpredictable harvesting-induced evolution relative to the classically applied multilocus model. Single-locus control of complex traits is thought to be uncommon, yet blocks of linked genes, such as those associated with some types of structural genomic variation, have emerged as taxonomically widespread phenomena. Inheritance of linked architectures resembles that of single loci, thus enabling single-locus-like modeling of polygenic adaptation. Yet, the number of loci, their effect sizes, and the degree of linkage among them all occur along a continuum. We review how linked architectures are often associated, directly or indirectly, with traits expected to be under selection from anthropogenic stressors and are likely to play a large role in adaptation to environmental disturbance. We suggest using single-locus models to explore evolutionary extremes and uncertainties when the trait architecture is unknown, refining parameters as genomic information becomes available, and explicitly incorporating linkage among loci when possible. By overestimating the complexity (e.g., number of independent loci) of the genomic architecture of traits under selection, we risk underestimating the complexity (e.g., nonlinearity) of their evolutionary dynamics.

Temporal Variation in Genetic Composition of Migratory Helicoverpa Zea in Peripheral Populations

Migrant populations of Helicoverpa zea (Boddie) captured during 2002, 2005, 2016, and 2018 from Landisville and Rock Springs in Pennsylvania, USA were genotyped using 85 single nucleotide polymorphism (SNP) markers. Samples (n = 702) genotyped were divided into 16 putative populations based on collection time and site. Fixation indices (F-statistics), analysis of molecular variance, and discriminant analysis of principal components were used to examine within and among population genetic variation. The observed and expected heterozygosity in putative populations ranged from 0.317-0.418 and 0.320-0.359, respectively. Broad range of FST (0.0-0.2742) and FIS (0.0-0.2330) values indicated different genotype frequencies between and within the populations, respectively. High genetic diversity within and low genetic differentiation between populations was found in 2002 and 2005. Interestingly, high genetic differentiation between populations from two collection sites observed in 2018 populations was not evident in within-site comparisons of putative populations collected on different dates during the season. The shift of H. zea population genetic makeup in 2018 may be influenced by multiple biotic and abiotic factors including tropical storms. Continued assessment of these peripheral populations of H. zea will be needed to assess the impacts of genetic changes on pest control and resistance management tactics.

How Natural Genetic Variation Shapes Behavior

Nervous systems allow animals to acutely respond and behaviorally adapt to changes and recurring patterns in their environment at multiple timescales-from milliseconds to years. Behavior is further shaped at intergenerational timescales by genetic variation, drift, and selection. This sophistication and flexibility of behavior makes it challenging to measure behavior consistently in individual subjects and to compare it across individuals. In spite of these challenges, careful behavioral observations in nature and controlled measurements in the laboratory, combined with modern technologies and powerful genetic approaches, have led to important discoveries about the way genetic variation shapes behavior. A critical mass of genes whose variation is known to modulate behavior in nature is finally accumulating, allowing us to recognize emerging patterns. In this review, we first discuss genetic mapping approaches useful for studying behavior. We then survey how variation acts at different levels-in environmental sensation, in internal neuronal circuits, and outside the nervous system altogether-and then discuss the sources and types of molecular variation linked to behavior and the mechanisms that shape such variation. We end by discussing remaining questions in the field.

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.

Adaptation to the abiotic environment in insects: the influence of variability on ecophysiology and evolutionary genomics

Advances in tools to gather environmental, phenotypic, and molecular data have accelerated our ability to detect abiotic drivers of variation across the genome-to-phenome spectrum in model and non-model insects. However, differences in the spatial and temporal resolution of these data sets may create gaps in our understanding of linkages between environment, genotype, and phenotype that yield missed or misleading results about adaptive variation. In this review we highlight sources of variability that might impact studies of phenotypic and 'omic environmental adaptation, challenges to collecting data at relevant scales, and possible solutions that link intensive fine-scale reductionist studies of mechanisms to large-scale biogeographic patterns.

Four Novel SNPs of MYO1A Gene Associated with Heat-Tolerance in Chinese Cattle

Based on the previous GWAS research related to bovine heat tolerance trait, this study aimed to explore the effect of myosin-1a (MYO1A) gene on bovine heat tolerance trait, and find the molecular markers related to the heat tolerance of Chinese cattle. In our study, four novel candidate SNPs highly conserved in B. indicus breeds but barely existed in B. taurus were identified in MYO1A gene according to Bovine Genome Variation Database and Selective Signatures (BGVD). PCR and DNA sequencing were used to genotype 1072 individuals including 34 Chinese indigenous cattle breeds as well as Angus and Indian zebu. Two synonymous mutations (rs208210464 and rs110123931), one missense mutation (rs209999142; Phe172Ser), and one intron mutation (rs135771836) were detected. The frequencies of mutant alleles of the four SNPs gradually increased from northern groups to southern groups of Chinese cattle, which was consistent with the distribution of various climatic conditions of China. Additionally, four SNPs were significantly associated with four climatic conditions including annual mean temperature (T), relative humidity (H), temperature-humidity index (THI), and average annual sunshine hours (100-cloudiness) (SR). Among these, rs209999142 and Hap 1/1 had better performance than others. Our results suggested that rs209999142 was associated with heat-tolerance trait and rs208210464, rs110123931, and rs135771836 showed high phenotypic effect on heat-tolerance trait because of the strong linkage with rs209999142. These SNPs could be used as candidates for marker-assisted selection (MAS) in cattle breeding.

Extreme Differences in Recombination Rate between the Genomes of a Solitary and a Social Bee

Social insect genomes exhibit the highest rates of crossing over observed in plants and animals. The evolutionary causes of these extreme rates are unknown. Insight can be gained by comparing recombination rate variation across the genomes of related social and solitary insects. Here, we compare the genomic recombination landscape of the highly social honey bee, Apis mellifera, with the solitary alfalfa leafcutter bee, Megachile rotundata, by analyzing patterns of linkage disequilibrium in population-scale genome sequencing data. We infer that average recombination rates are extremely elevated in A. mellifera compared with M. rotundata. However, our results indicate that similar factors control the distribution of crossovers in the genomes of both species. Recombination rate is significantly reduced in coding regions in both species, with genes inferred to be germline methylated having particularly low rates. Genes with worker-biased patterns of expression in A. mellifera and their orthologs in M. rotundata have higher than average recombination rates in both species, suggesting that selection for higher diversity in genes involved in worker caste functions in social taxa is not the explanation for these elevated rates. Furthermore, we find no evidence that recombination has modulated the efficacy of selection among genes during bee evolution, which does not support the hypothesis that high recombination rates facilitated positive selection for new functions in social insects. Our results indicate that the evolution of sociality in insects likely entailed selection on modifiers that increased recombination rates genome wide, but that the genomic recombination landscape is determined by the same factors.

Coevolution of Genome Architecture and Social Behavior

Although social behavior can have a strong genetic component, it can also result in selection on genome structure and function, thereby influencing the evolution of the genome itself. Here we explore the bidirectional links between social behavior and genome architecture by considering variation in social and/or mating behavior among populations (social polymorphisms) and across closely related species. We propose that social behavior can influence genome architecture via associated demographic changes due to social living. We establish guidelines to exploit emerging whole-genome sequences using analytical approaches that examine genome structure and function at different levels (regulatory vs structural variation) from the perspective of both molecular biology and population genetics in an ecological context.

Social Parasitism in the Honeybee (Apis mellifera) Is Not Controlled by a Single SNP

The Cape bee (Apis mellifera capensis) is a subspecies of the honeybee, in which workers commonly lay diploid unfertilized eggs via a process known as thelytoky. A recent study aimed to map the genetic basis of this trait in the progeny of a single capensis queen where workers laid either diploid (thelytokous) or haploid (arrhenotokous) eggs. A nonsynonymous single nucleotide polymorphism (SNP) in a gene of unknown function was reported to be strongly associated with thelytoky in this colony. Here, we analyze genome sequences from a global sample of A. mellifera and identify populations where the proposed thelytoky allele at this SNP is common but thelytoky is absent. We also analyze genome sequences of three capensis queens produced by thelytoky and find that, contrary to predictions, they do not carry the proposed thelytoky allele. The proposed SNP is therefore neither sufficient nor required to produce thelytoky in A. mellifera.

Intercontinental karyotype-environment parallelism supports a role for a chromosomal inversion in local adaptation in a seaweed fly

Large chromosomal rearrangements are thought to facilitate adaptation to heterogeneous environments by limiting genomic recombination. Indeed, inversions have been implicated in adaptation along environmental clines and in ecotype specialization. Here, we combine classical ecological studies and population genetics to investigate an inversion polymorphism previously documented in Europe among natural populations of the seaweed fly Coelopa frigida along a latitudinal cline in North America. We test if the inversion is present in North America and polymorphic, assess which environmental conditions modulate the inversion karyotype frequencies, and document the relationship between inversion karyotype and adult size. We sampled nearly 2000 flies from 20 populations along several environmental gradients to quantify associations of inversion frequencies to heterogeneous environmental variables. Genotyping and phenotyping showed a widespread and conserved inversion polymorphism between Europe and America. Variation in inversion frequency was significantly associated with environmental factors, with parallel patterns between continents, indicating that the inversion may play a role in local adaptation. The three karyotypes of the inversion are differently favoured across micro-habitats and represent life-history strategies likely to be maintained by the collective action of several mechanisms of balancing selection. Our study adds to the mounting evidence that inversions are facilitators of adaptation and enhance within-species diversity.

A hybrid de novo genome assembly of the honeybee, Apis mellifera, with chromosome-length scaffolds

Background

The ability to generate long sequencing reads and access long-range linkage information is revolutionizing the quality and completeness of genome assemblies. Here we use a hybrid approach that combines data from four genome sequencing and mapping technologies to generate a new genome assembly of the honeybee Apis mellifera. We first generated contigs based on PacBio sequencing libraries, which were then merged with linked-read 10x Chromium data followed by scaffolding using a BioNano optical genome map and a Hi-C chromatin interaction map, complemented by a genetic linkage map.

Results

Each of the assembly steps reduced the number of gaps and incorporated a substantial amount of additional sequence into scaffolds. The new assembly (Amel_HAv3) is significantly more contiguous and complete than the previous one (Amel_4.5), based mainly on Sanger sequencing reads. N50 of contigs is 120-fold higher (5.381 Mbp compared to 0.053 Mbp) and we anchor > 98% of the sequence to chromosomes. All of the 16 chromosomes are represented as single scaffolds with an average of three sequence gaps per chromosome. The improvements are largely due to the inclusion of repetitive sequence that was unplaced in previous assemblies. In particular, our assembly is highly contiguous across centromeres and telomeres and includes hundreds of AvaI and AluI repeats associated with these features.

Conclusions

The improved assembly will be of utility for refining gene models, studying genome function, mapping functional genetic variation, identification of structural variants, and comparative genomics.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5642-0) contains supplementary material, which is available to authorized users.

How chromosomal rearrangements shape adaptation and speciation: Case studies in Drosophila pseudoobscura and its sibling species Drosophila persimilis

The gene arrangements of Drosophila have played a prominent role in the history of evolutionary biology from the original quantification of genetic diversity to current studies of the mechanisms for the origin and establishment of new inversion mutations within populations and their subsequent fixation between species supporting reproductive barriers. This review examines the genetic causes and consequences of inversions as recombination suppressors and the role that recombination suppression plays in establishing inversions in populations as they are involved in adaptation within heterogeneous environments. This often results in the formation of clines of gene arrangement frequencies among populations. Recombination suppression leads to the differentiation of the gene arrangements which may accelerate the accumulation of fixed genetic differences among populations. If these fixed mutations cause incompatibilities, then inversions pose important reproductive barriers between species. This review uses the evolution of inversions in Drosophila pseudoobscura and D. persimilis as a case study for how inversions originate, establish and contribute to the evolution of reproductive isolation.

A Single SNP Turns a Social Honey Bee (Apis mellifera) Worker into a Selfish Parasite

The evolution of altruism in complex insect societies is arguably one of the major transitions in evolution and inclusive fitness theory plausibly explains why this is an evolutionary stable strategy. Yet, workers of the South African Cape honey bee (Apis mellifera capensis) can reverse to selfish behavior by becoming social parasites and parthenogenetically producing female offspring (thelytoky). Using a joint mapping and population genomics approach, in combination with a time-course transcript abundance dynamics analysis, we show that a single nucleotide polymorphism at the mapped thelytoky locus (Th) is associated with the iconic thelytokous phenotype. Th forms a linkage group with the ecdysis-triggering hormone receptor (Ethr) within a nonrecombining region under strong selection in the genome. A balanced detrimental allele system plausibly explains why the trait is specific to A. m. capensis and cannot easily establish itself into genomes of other honey bee subspecies.

The genomic basis of adaptation to high-altitude habitats in the eastern honey bee (Apis cerana)

The eastern honey bee (Apis cerana) is of central importance for agriculture in Asia. It has adapted to a wide variety of environmental conditions across its native range in southern and eastern Asia, which includes high-altitude regions. eastern honey bees inhabiting mountains differ morphologically from neighbouring lowland populations and may also exhibit differences in physiology and behaviour. We compared the genomes of 60 eastern honey bees collected from high and low altitudes in Yunnan and Gansu provinces, China, to infer their evolutionary history and to identify candidate genes that may underlie adaptation to high altitude. Using a combination of F_ST -based statistics, long-range haplotype tests and population branch statistics, we identified several regions of the genome that appear to have been under positive selection. These candidate regions were strongly enriched for coding sequences and had high haplotype homozygosity and increased divergence specifically in highland bee populations, suggesting they have been subjected to recent selection in high-altitude habitats. Candidate loci in these genomic regions included genes related to reproduction and feeding behaviour in honey bees. Functional investigation of these candidate loci is necessary to fully understand the mechanisms of adaptation to high-altitude habitats in the eastern honey bee.

Recent advances in population and quantitative genomics of honey bees

The increase in the availability of individual Apis mellifera genomes has resulted in significant progress toward understanding the evolution and adaptation of the honey bee. These efforts have identified new subspecies, evolutionary lineages, and a significant number of genes involved with adaptations and colony-level quantitative traits. Many studies have also developed genetic assays that are being used to monitor the movement and admixture of honey bee populations. These resources are valuable for conservation and breeding programs that seek to improve the economic value of colonies or preserve locally adapted populations and subspecies. This review provides a brief discussion on how population and quantitative genomic studies has improved our understanding of the honey bee.