A soft selective sweep during rapid evolution of gentle behaviour in an Africanized honeybee

Genomic and transcriptomic insights into vitamin A-induced thermogenesis and gene reuse as a cold adaptation strategy in wild boars

Wild boars inhabit diverse climates, including frigid regions like Siberia, but their migration history and cold adaptation mechanisms into high latitudes remain poorly understood. We constructed the most comprehensive wild boar whole-genome variant dataset to date, comprising 124 samples from tropical to frigid zones, among which 47 Russian, 8 South Chinese and 3 Vietnamese wild boars were newly supplemented. We also gathered 75 high-quality RNA-seq datasets from 10 tissues of 6 wild boars from Russia and 6 from southern China. Demographic analysis revealed the appearance of Russian wild boars in Far East of Asia (RUA) and Europe (RUE) after the last glacial maximum till ~ 10 thousand years ago. Recent gene flow (<100 years) from RUA to RUE reflects human-mediated introductions. Cold-region wild boars exhibit strong selection signatures indicative of genetic adaptation to cold climates. Further pathway and transcriptomic analyses reveal a novel cold resistance mechanism centered on enhanced vitamin A metabolism and catalysis, involving the reuse of UGT2B31 and rhythm regulation by ANGPTL8, RLN3 and ZBTB20. This may compensate for the pig's lack of brown fat/UCP1 thermogenesis. These findings provide new insights into the molecular basis of cold adaptation and improve our understanding of Eurasian wild boar migration history.

Evolution of parasites in the Anthropocene: new pressures, new adaptive directions

The Anthropocene is seeing the human footprint rapidly spreading to all of Earth's ecosystems. The fast-changing biotic and abiotic conditions experienced by all organisms are exerting new and strong selective pressures, and there is a growing list of examples of human-induced evolution in response to anthropogenic impacts. No organism is exempt from these novel selective pressures. Here, we synthesise current knowledge on human-induced evolution in eukaryotic parasites of animals, and present a multidisciplinary framework for its study and monitoring. Parasites generally have short generation times and huge fecundity, features that predispose them for rapid evolution. We begin by reviewing evidence that parasites often have substantial standing genetic variation, and examples of their rapid evolution both under conditions of livestock production and in serial passage experiments. We then present a two-step conceptual overview of the causal chain linking anthropogenic impacts to parasite evolution. First, we review the major anthropogenic factors impacting parasites, and identify the selective pressures they exert on parasites through increased mortality of either infective stages or adult parasites, or through changes in host density, quality or immunity. Second, we discuss what new phenotypic traits are likely to be favoured by the new selective pressures resulting from altered parasite mortality or host changes; we focus mostly on parasite virulence and basic life-history traits, as these most directly influence the transmission success of parasites and the pathology they induce. To illustrate the kinds of evolutionary changes in parasites anticipated in the Anthropocene, we present a few scenarios, either already documented or hypothetical but plausible, involving parasite taxa in livestock, aquaculture and natural systems. Finally, we offer several approaches for investigations and real-time monitoring of rapid, human-induced evolution in parasites, ranging from controlled experiments to the use of state-of-the-art genomic tools. The implications of fast-evolving parasites in the Anthropocene for disease emergence and the dynamics of infections in domestic animals and wildlife are concerning. Broader recognition that it is not only the conditions for parasite transmission that are changing, but the parasites themselves, is needed to meet better the challenges ahead.

Limited population structure but signals of recent selection in introduced African Fig Fly (Zaprionus indianus) in North America

Invasive species have devastating consequences for human health, food security, and the environment. Many invasive species adapt to new ecological niches following invasion, but little is known about the early steps of adaptation. Here we examine population genomics of a recently introduced drosophilid in North America, the African Fig Fly, Zaprionus indianus. This species is likely intolerant of subfreezing temperatures and recolonizes temperate environments yearly. We generated a new chromosome-level genome assembly for Z. indianus. Using resequencing of over 200 North American individuals collected over four years in temperate Virginia, plus a single collection from subtropical Florida, we tested for signatures of recolonization, population structure, and adaptation within invasive populations. We show founding populations are sometimes small and contain close genetic relatives, yet temporal population structure and differentiation of populations is mostly absent across recurrent recolonization events. Although we find limited signals of genome-wide spatial or temporal population structure, we identify haplotypes on the X chromosome that are repeatedly differentiated between Virginia and Florida populations. These haplotypes show signatures of natural selection and are not found in African populations. We also find evidence for several large structural polymorphisms segregating within North America populations and show X chromosome evolution in invasive populations is strikingly different from the autosomes. These results show that despite limited population structure, populations may rapidly evolve genetic differences early in an invasion. Further uncovering how these genomic regions influence invasive potential and success in new environments will advance our understanding of how organisms evolve in changing environments.

HBeeID: a molecular tool that identifies honey bee subspecies from different geographic populations

BACKGROUND

Honey bees are the principal commercial pollinators. Along with other arthropods, they are increasingly under threat from anthropogenic factors such as the incursion of invasive honey bee subspecies, pathogens and parasites. Better tools are needed to identify bee subspecies. Genomic data for economic and ecologically important organisms is increasing, but in its basic form its practical application to address ecological problems is limited.

RESULTS

We introduce HBeeID a means to identify honey bees. The tool utilizes a knowledge-based network and diagnostic SNPs identified by discriminant analysis of principle components and hierarchical agglomerative clustering. Tests of HBeeID showed that it identifies African, Americas-Africanized, Asian, and European honey bees with a high degree of certainty even when samples lack the full 272 SNPs of HBeeID. Its prediction capacity decreases with highly admixed samples.

CONCLUSION

HBeeID is a high-resolution genomic, SNP based tool, that can be used to identify honey bees and screen species that are invasive. Its flexible design allows for future improvements via sample data additions from other localities.

Post-invasion selection acts on standing genetic variation despite a severe founding bottleneck

Invasive populations often have lower genetic diversity relative to the native-range populations from which they derive.1,2 Despite this, many biological invaders succeed in their new environments, in part due to rapid adaptation.3,4,5,6 Therefore, the role of genetic bottlenecks in constraining the adaptation of invaders is debated.7,8,9,10 Here, we use whole-genome resequencing of samples from a 10-year time-series dataset, representing the natural invasion of the Asian honey bee (Apis cerana) in Australia, to investigate natural selection occurring in the aftermath of a founding event. We find that Australia's A. cerana population was founded by as few as one colony, whose arrival was followed by a period of rapid population expansion associated with an increase of rare variants.11 The bottleneck resulted in a steep loss of overall genetic diversity, yet we nevertheless detected loci with signatures of positive selection during the first years post-invasion. When we investigated the origin of alleles under selection, we found that selection acted primarily on the variation introduced by founders and not on the variants that arose post-invasion by mutation. In all, our data highlight that selection on standing genetic variation can occur in the early years post-invasion, even where founding bottlenecks are severe.

Mutation bias and adaptation in bacteria

Genetic mutation, which provides the raw material for evolutionary adaptation, is largely a stochastic force. However, there is ample evidence showing that mutations can also exhibit strong biases, with some mutation types and certain genomic positions mutating more often than others. It is becoming increasingly clear that mutational bias can play a role in determining adaptive outcomes in bacteria in both the laboratory and the clinic. As such, understanding the causes and consequences of mutation bias can help microbiologists to anticipate and predict adaptive outcomes. In this review, we provide an overview of the mechanisms and features of the bacterial genome that cause mutational biases to occur. We then describe the environmental triggers that drive these mechanisms to be more potent and outline the adaptive scenarios where mutation bias can synergize with natural selection to define evolutionary outcomes. We conclude by describing how understanding mutagenic genomic features can help microbiologists predict areas sensitive to mutational bias, and finish by outlining future work that will help us achieve more accurate evolutionary forecasts.

Single-cell dissection of aggression in honeybee colonies

Understanding how genotypic variation results in phenotypic variation is especially difficult for collective behaviour because group phenotypes arise from complex interactions among group members. A genome-wide association study identified hundreds of genes associated with colony-level variation in honeybee aggression, many of which also showed strong signals of positive selection, but the influence of these 'colony aggression genes' on brain function was unknown. Here we use single-cell (sc) transcriptomics and gene regulatory network (GRN) analyses to test the hypothesis that genetic variation for colony aggression influences individual differences in brain gene expression and/or gene regulation. We compared soldiers, which respond to territorial intrusion with stinging attacks, and foragers, which do not. Colony environment showed stronger influences on soldier-forager differences in brain gene regulation compared with brain gene expression. GRN plasticity was strongly associated with colony aggression, with larger differences in GRN dynamics detected between soldiers and foragers from more aggressive relative to less aggressive colonies. The regulatory dynamics of subnetworks composed of genes associated with colony aggression genes were more strongly correlated with each other across different cell types and brain regions relative to other genes, especially in brain regions involved with olfaction and vision and multimodal sensory integration, which are known to mediate bee aggression. These results show how group genetics can shape a collective phenotype by modulating individual brain gene regulatory network architecture.

Population structure and genome-wide evolutionary signatures reveal putative climate-driven habitat change and local adaptation in the large yellow croaker

The large yellow croaker (Larimichthys crocea) is one of the most economically valuable marine fish in China and is a notable species in ecological studies owing to a serious collapse of wild germplasm in the past few decades. The stock division and species distribution, which have important implications for ecological protection, germplasm recovery, and fishery resource management, have been debated since the 1960s. However, it is still uncertain even how many stocks exist in this species. To address this, we evaluated the fine-scale genetic structure of large yellow croaker populations distributed along the eastern and southern Chinese coastline based on 7.64 million SNP markers. Compared with the widely accepted stock boundaries proposed in the 1960s, our results revealed that a climate-driven habitat change probably occurred between the Naozhou (Nanhai) Stock and the Ming-Yuedong (Mindong) Stock. The boundary between these two stocks might have shifted northwards from the Pearl River Estuary to the northern area of the Taiwan Strait, accompanied by highly asymmetric introgression. In addition, we found divergent landscapes of natural selection between the stocks inhabiting northern and southern areas. The northern population exhibited highly agminated signatures of strong natural selection in genes related to developmental processes, whereas moderate and interspersed selective signatures were detected in many immune-related genes in the southern populations. These findings establish the stock status and genome-wide evolutionary landscapes of large yellow croaker, providing a basis for conservation, fisheries management and further evolutionary biology studies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00165-2.

Genome-wide evolutionary signatures of climate adaptation in spotted sea bass inhabiting different latitudinal regions

Consideration of the thermal adaptation of species is essential in both evolutionary biology and climate-change biology because it frequently leads to latitudinal gradients of various phenotypes among populations. The spotted sea bass (Lateolabrax maculatus) has a broad latitudinal distribution range along the marginal seas of the Northwest Pacific and thus provides an excellent teleost model for population genetic and climate adaptation studies. We generated over 8.57 million SNP loci using whole-genome resequencing from 100 samples collected at 14 geographic sites (five or ten samples per site). We estimated the genetic structure of the sampled fish and clustered them into three highly differentiated populations. The genetic differentiation pattern estimated by multivariable models combining geographic distance and sea surface temperature differences suggests that isolation by distance and isolation by environment both have significant effects on this species. Further investigation of genome-wide evolutionary signatures of climate adaptation identified many genes related to growth, muscle contraction, and vision that are under positive natural selection. Moreover, the contrasting patterns of natural selection in high-latitude and low-latitude populations prompted different strategies of trade-offs between growth rate and other traits that may play an essential role in adaptation to different local climates. Our results offer an opportunity to better understand the genetic basis of the phenotypic variation in eurythermal fishes inhabiting different climatic regions.

Characterization of the Bacterial Communities Inhabiting Tropical Propolis of Puerto Rico

Propolis is a resinous material produced by honeybees from different plant sources and used in the hive as a building material and to protect the colony from parasites and pathogens. Despite its antimicrobial properties, recent studies showed that propolis hosts diverse microbial strains, some with great antimicrobial potential. In this study, the first description of the bacterial community of propolis produced by the gentle Africanized honeybee was reported. Propolis was sampled from hives of two different geographic areas of Puerto Rico (PR, USA), and the associated microbiota investigated by both cultivation and metataxonomic approaches. Metabarcoding analysis showed appreciable bacterial diversity in both areas and statistically significant dissimilarity in the taxa composition of the two areas, probably due to the different climatic conditions. Both metabarcoding and cultivation data revealed the presence of taxa already detected in other hive components and compatible with the bee's foraging environment. Isolated bacteria and propolis extracts showed antimicrobial activity against Gram-positive and Gram-negative bacterial tester strains. These results support the hypothesis that the propolis microbiota could contribute to propolis' antimicrobial properties.

Rapid adaptation in a fast-changing world: Emerging insights from insect genomics

Many researchers have questioned the ability of biota to adapt to rapid anthropogenic environmental shifts. Here, we synthesize emerging genomic evidence for rapid insect evolution in response to human pressure. These new data reveal diverse genomic mechanisms (single locus, polygenic, structural shifts; introgression) underpinning rapid adaptive responses to a variety of anthropogenic selective pressures. While the effects of some human impacts (e.g. pollution; pesticides) have been previously documented, here we highlight startling new evidence for rapid evolutionary responses to additional anthropogenic processes such as deforestation. These recent findings indicate that diverse insect assemblages can indeed respond dynamically to major anthropogenic evolutionary challenges. Our synthesis also emphasizes the critical roles of genomic architecture, standing variation and gene flow in maintaining future adaptive potential. Broadly, it is clear that genomic approaches are essential for predicting, monitoring and responding to ongoing anthropogenic biodiversity shifts in a fast-changing world.

Russian honey bee genotype identification through enhanced marker panel set

Russian honey bees (RHB) are a breeding population developed by USDA-ARS as an effort to provide Varroa-resistant honey bees to beekeepers. The selection strategy for this breeding population was the first in honey bees to incorporate genetic stock identification (GSI). The original GSI approach has been in use for over a decade, and though effective, novel technologies and analytical approaches recently developed provide an opportunity for improvement. Here we outline a novel genotyping assay that capitalizes on the markers used in the GSI as well as new loci recently identified in a whole genome pooled study of commercial honey bee stocks. Our approach utilizes a microfluidic platform and machine learning analyses to arrive at an accurate, high throughput assay. This novel approach provides an improved tool that can be readily incorporated into breeding decisions towards healthier more productive bees.

Hard versus soft selective sweeps during domestication and improvement in soybean

Genome scans for selection can provide an efficient way to dissect the genetic basis of domestication traits and understand mechanisms of adaptation during crop evolution. Selection involving soft sweeps (simultaneous selection for multiple alleles) is probably common in plant genomes but is under-studied, and few if any studies have systematically scanned for soft sweeps in the context of crop domestication. Using genome resequencing data from 302 wild and domesticated soybean accessions, we conducted selection scans using five widely employed statistics to identify selection candidates under classical (hard) and soft sweeps. Across the genome, inferred hard sweeps are predominant in domesticated soybean landraces and improved varieties, whereas soft sweeps are more prevalent in a representative subpopulation of the wild ancestor. Six domestication-related genes, representing both hard and soft sweeps and different stages of domestication, were used as positive controls to assess the detectability of domestication-associated sweeps. Performance of various test statistics suggests that differentiation-based (FST ) methods are robust for detecting complete hard sweeps, and that LD-based strategies perform well for identifying recent/ongoing sweeps; however, none of the test statistics detected a known soft sweep we previously documented at the domestication gene Dt1. Genome scans yielded a set of 66 candidate loci that were identified by both differentiation-based and LD-based (iHH) methods; notably, this shared set overlaps with many previously identified QTLs for soybean domestication/improvement traits. Collectively, our results will help to advance genetic characterizations of soybean domestication traits and shed light on selection modes involved in adaptation in domesticated plant species.

The Movement of Western Honey Bees (Apis mellifera L.) Among United States and Territories: History, Benefits, Risks, and Mitigation Strategies

Beekeeping is a cornerstone activity that has led to the human-mediated, global spread of western honey bees (Apis mellifera L.) outside their native range of Europe, western Asia, and Africa. The exportation/importation of honey bees (i.e., transfer of honey bees or germplasm between countries) is regulated at the national level in many countries. Honey bees were first imported into the United States in the early 1600’s. Today, honey bee movement (i.e., transport of honey bees among states and territories) is regulated within the United States at the state, territory, and federal levels. At the federal level, honey bees present in the country (in any state or territory) can be moved among states and territories without federal restriction, with the exception of movement to Hawaii. In contrast, regulations at the state and territory levels vary substantially, ranging from no additional regulations beyond those stipulated at the federal level, to strict regulations for the introduction of live colonies, packaged bees, or queens. This variability can lead to inconsistencies in the application of regulations regarding the movement of honey bees among states and territories. In November 2020, we convened a technical working group (TWG), composed of academic and USDA personnel, to review and summarize the (1) history of honey bee importation into/movement within the United States, (2) current regulations regarding honey bee movement and case studies on the application of those regulations, (3) benefits associated with moving honey bees within the United States, (4) risks associated with moving honey bees within the United States, and (5) risk mitigation strategies. This review will be helpful for developing standardized best practices for the safe movement of honey bees between the 48 contiguous states and other states/territories within the United States.

Genomic Signatures of Recent Adaptation in a Wild Bumblebee

Environmental changes threaten insect pollinators, creating risks for agriculture and ecosystem stability. Despite their importance, we know little about how wild insects respond to environmental pressures. To understand the genomic bases of adaptation in an ecologically important pollinator, we analyzed genomes of Bombus terrestris bumblebees collected across Great Britain. We reveal extensive genetic diversity within this population, and strong signatures of recent adaptation throughout the genome affecting key processes including neurobiology and wing development. We also discover unusual features of the genome, including a region containing 53 genes that lacks genetic diversity in many bee species, and a horizontal gene transfer from a Wolbachia bacteria. Overall, the genetic diversity we observe and how it is distributed throughout the genome and the population should support the resilience of this important pollinator species to ongoing and future selective pressures. Applying our approach to more species should help understand how they can differ in their adaptive potential, and to develop conservation strategies for those most at risk.

Defensive Behavior and Morphometric Variation in Apis mellifera Colonies From Two Different Agro-Ecological Zones of North-Western Argentina

European lineages of Apis mellifera were first introduced into America for beekeeping purposes. A subsequent introduction and accidental release of A. m. scutellata resulted in hybridization events that gave rise to Africanized populations that rapidly spread throughout the continent. In Argentina, Africanized honey bees (AHBs) have been mostly detected in northern regions of the territory, and represent a valuable genetic resource for the selection of stocks with advantageous characteristics for beekeeping. The objective of the present study was to profile honey bee colonies of wild origin with potential beneficial traits for apiculture using morphological, molecular and behavioral traits. Honey bee colonies chosen for evaluation were located in two different agro-ecological regions in north-western Argentina (Tucumán province): The Chaco Depressed Plain (Leales apiary) and the Piedmont (Famaillá apiary). Each apiary was surveyed three times during the 2017–2018 season (mid-season, wintertime, and early spring) for: brood population, phoretic Varroa level and defensive behavior (run, fly, sting, and hang). At the midpoint of the beekeeping season colonies were also characterized by morphometry (45 variables) and mitochondrial haplotypes (COI–COII intergenic region). Apiaries studied showed similar patterns throughout the beekeeping season, for most of the characteristics monitored. However, significant variation in defensive behavior parameters was found between apiaries at the different times of evaluation. Twelve of 45 morphometric variables also showed significant differences between apiaries. The mitochondrial haplotype analysis revealed a high representation of African A4 and A1 haplotypes (91%) in both apiaries. Haplotype variation was associated with morphometric and behavioral traits. Multivariate analyses [principal component analysis (PCA) and principal coordinate analysis (PCoA)] including morphometric and behavior variables explained 65.3% (PCA) and 48.1% (PCoA) of the variability observed between colonies in the first two components. Several morphometric parameters and “fly” behavior were mainly associated with the separation of the colonies. The results from this study point to a possible association between morphometric and behavioral variation and the adaptation of honey bee colonies to differential agro-ecological conditions. We discuss how the detected variation between apiaries can be used for the selection and preservation of honey bee ecotypes in regional breeding programs.

Environmental correlates of genetic variation in the invasive European starling in North America

Populations of invasive species that colonize and spread in novel environments may differentiate both through demographic processes and local selection. European starlings (Sturnus vulgaris) were introduced to New York in 1890 and subsequently spread throughout North America, becoming one of the most widespread and numerous bird species on the continent. Genome-wide comparisons across starling individuals and populations can identify demographic and/or selective factors that facilitated this rapid and successful expansion. We investigated patterns of genomic diversity and differentiation using reduced-representation genome sequencing of 17 winter-season sampling sites. Consistent with this species' high dispersal rate and rapid expansion history, we found low geographical differentiation and few F_ST outliers even at a continental scale. Despite starting from a founding population of ~180 individuals, North American starlings show only a moderate genetic bottleneck, and models suggest a dramatic increase in effective population size since introduction. In genotype-environment associations we found that ~200 single-nucleotide polymorphisms are correlated with temperature and/or precipitation against a background of negligible genome- and range-wide divergence. Given this evidence, we suggest that local adaptation in North American starlings may have evolved rapidly even in this wide-ranging and evolutionarily young system. This survey of genomic signatures of expansion in North American starlings is the most comprehensive to date and complements ongoing studies of world-wide local adaptation in these highly dispersive and invasive birds.

Effects of aversive conditioning on expression of physiological stress in honey bees (Apis mellifera)

Stress is defined as any deviation from an organism's baseline physiological levels. Therefore, introduction of new stimuli and information, such as in learning, can be defined as a stressor. A large body of research exists examining the role that stress plays in learning, but virtually none addresses whether or not learning itself is a measurable cause of stress. The current study used a wide variety of learning centric stress responses. Researchers examined changes in expression of ten stress and learning related genes in various physiological systems in domesticated honey bees (Apis mellifera) as a result of exposure to an aversive conditioning task. Gene expression was examined using quantitative real-time polymerase chain reaction following the learning task. Results indicate that learning affects expression of some stress related genes.

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.

Selection and hybridization shaped the rapid spread of African honey bee ancestry in the Americas

Recent biological invasions offer 'natural' laboratories to understand the genetics and ecology of adaptation, hybridization, and range limits. One of the most impressive and well-documented biological invasions of the 20th century began in 1957 when Apis mellifera scutellata honey bees swarmed out of managed experimental colonies in Brazil. This newly-imported subspecies, native to southern and eastern Africa, both hybridized with and out-competed previously-introduced European honey bee subspecies. Populations of scutellata-European hybrid honey bees rapidly expanded and spread across much of the Americas in less than 50 years. We use broad geographic sampling and whole genome sequencing of over 300 bees to map the distribution of scutellata ancestry where the northern and southern invasions have presently stalled, forming replicated hybrid zones with European bee populations in California and Argentina. California is much farther from Brazil, yet these hybrid zones occur at very similar latitudes, consistent with the invasion having reached a climate barrier. At these range limits, we observe genome-wide clines for scutellata ancestry, and parallel clines for wing length that span hundreds of kilometers, supporting a smooth transition from climates favoring scutellata-European hybrid bees to climates where they cannot survive winter. We find no large effect loci maintaining exceptionally steep ancestry transitions. Instead, we find most individual loci have concordant ancestry clines across South America, with a build-up of somewhat steeper clines in regions of the genome with low recombination rates, consistent with many loci of small effect contributing to climate-associated fitness trade-offs. Additionally, we find no substantial reductions in genetic diversity associated with rapid expansions nor complete dropout of scutellata ancestry at any individual loci on either continent, which suggests that the competitive fitness advantage of scutellata ancestry at lower latitudes has a polygenic basis and that scutellata-European hybrid bees maintained large population sizes during their invasion. To test for parallel selection across continents, we develop a null model that accounts for drift in ancestry frequencies during the rapid expansion. We identify several peaks within a larger genomic region where selection has pushed scutellata ancestry to high frequency hundreds of kilometers past the present cline centers in both North and South America and that may underlie high-fitness traits driving the invasion.

Defense Response in Brazilian Honey Bees (Apis mellifera scutellata × spp.) Is Underpinned by Complex Patterns of Admixture

In 1957, an invasive and highly defensive honey bee began to spread across Brazil. In the previous year, Brazilian researchers hoped to produce a subtropical-adapted honey bee by crossing local commercial honey bees (of European origin) with a South African honey bee subspecies (Apis mellifera scutellata; an A-lineage honey bee subspecies). The resulting cross-African hybrid honey bees (AHBs)-escaped from their enclosure and spread through the Americas. Today, AHB is the most common honey bee from Northern Argentina to the Southern United States. AHBs are much more likely to sting nest intruders than managed European-derived honey bee colonies. Previous studies have explored how genetic variation contributes to differences in defense response between European-derived honey bee and AHB. Although this work demonstrated very strong genetic effects on defense response, they have yet to pinpoint which genes influence variation in defense response within AHBs, specifically. We quantified defense response for 116 colonies in Brazil and performed pooled sequencing on the most phenotypically divergent samples. We identified 65 loci containing 322 genes that were significantly associated with defense response. Loci were strongly associated with metabolic function, consistent with previous functional genomic analyses of this phenotype. Additionally, defense-associated loci had nonrandom and unexpected patterns of admixture. Defense response was not simply the product of more A-lineage honey bee ancestry as previously assumed, but rather an interaction between A-lineage and European alleles. Our results suggest that a combination of A-lineage and European alleles play roles in defensive behavior in AHBs.

Genomic regions influencing aggressive behavior in honey bees are defined by colony allele frequencies

For social animals, the genotypes of group members affect the social environment, and thus individual behavior, often indirectly. We used genome-wide association studies (GWAS) to determine the influence of individual vs. group genotypes on aggression in honey bees. Aggression in honey bees arises from the coordinated actions of colony members, primarily nonreproductive "soldier" bees, and thus, experiences evolutionary selection at the colony level. Here, we show that individual behavior is influenced by colony environment, which in turn, is shaped by allele frequency within colonies. Using a population with a range of aggression, we sequenced individual whole genomes and looked for genotype-behavior associations within colonies in a common environment. There were no significant correlations between individual aggression and specific alleles. By contrast, we found strong correlations between colony aggression and the frequencies of specific alleles within colonies, despite a small number of colonies. Associations at the colony level were highly significant and were very similar among both soldiers and foragers, but they covaried with one another. One strongly significant association peak, containing an ortholog of the Drosophila sensory gene dpr4 on linkage group (chromosome) 7, showed strong signals of both selection and admixture during the evolution of gentleness in a honey bee population. We thus found links between colony genetics and group behavior and also, molecular evidence for group-level selection, acting at the colony level. We conclude that group genetics dominates individual genetics in determining the fatal decision of honey bees to sting.

Evolutionary dynamics of recent selection on cognitive abilities

Cognitive abilities can vary dramatically among species. The relative importance of social and ecological challenges in shaping cognitive evolution has been the subject of a long-running and recently renewed debate, but little work has sought to understand the selective dynamics underlying the evolution of cognitive abilities. Here, we investigate recent selection related to cognition in the paper wasp Polistes fuscatus-a wasp that has uniquely evolved visual individual recognition abilities. We generate high quality de novo genome assemblies and population genomic resources for multiple species of paper wasps and use a population genomic framework to interrogate the probable mode and tempo of cognitive evolution. Recent, strong, hard selective sweeps in P. fuscatus contain loci annotated with functions in long-term memory formation, mushroom body development, and visual processing, traits which have recently evolved in association with individual recognition. The homologous pathways are not under selection in closely related wasps that lack individual recognition. Indeed, the prevalence of candidate cognition loci within the strongest selective sweeps suggests that the evolution of cognitive abilities has been among the strongest selection pressures in P. fuscatus' recent evolutionary history. Detailed analyses of selective sweeps containing candidate cognition loci reveal multiple cases of hard selective sweeps within the last few thousand years on de novo mutations, mainly in noncoding regions. These data provide unprecedented insight into some of the processes by which cognition evolves.

Colonization history and population differentiation of the Honey Bees (Apis mellifera L.) in Puerto Rico

Honey bees (Apis mellifera L.) are the primary commercial pollinators across the world. The subspecies A. m. scutellata originated in Africa and was introduced to the Americas in 1956. For the last 60 years, it hybridized successfully with European subspecies, previous residents in the area. The result of this hybridization was called Africanized honey bee (AHB). AHB has spread since then, arriving to Puerto Rico (PR) in 1994. The honey bee population on the island acquired a mosaic of features from AHB or the European honey bee (EHB). AHB in Puerto Rico shows a major distinctive characteristic, docile behavior, and is called gentle Africanized honey bees (gAHB). We used 917 SNPs to examine the population structure, genetic differentiation, origin, and history of range expansion and colonization of gAHB in PR. We compared gAHB to populations that span the current distribution of A. mellifera worldwide. The gAHB population is shown to be a single population that differs genetically from the examined populations of AHB. Texas and PR groups are the closest genetically. Our results support the hypothesis that the Texas AHB population is the source of gAHB in Puerto Rico.

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.