Honeybees' novel complementary sex-determining system: function and origin.

Complementary sex determination regulates female and male development in honeybees (Apis mellifera) via heterozygous versus homo-/hemizygous genotypes of the csd (complementary sex determiner) gene involving numerous naturally occurring alleles. This lineage-specific function offers a rare opportunity to understand an undescribed regulatory mechanism and the molecular evolutionary path leading to this mechanism. We reviewed recent advances in understanding how Csd recognizes different versus identical protein variants, how these variants regulate downstream pathways and sexual differentiation, and how this mechanism has evolved and been shaped by evolutionary forces. Finally, we highlighted the shared regulatory principles of sex determination despite the diversity of primary signals and demonstrated that lineage-specific mutations are very informative for characterizing newly evolved functions.

A Secreted RNA Binding Protein Forms RNA-Stabilizing Granules in the Honeybee Royal Jelly.

RNA flow between organisms has been documented within and among different kingdoms of life. Recently, we demonstrated horizontal RNA transfer between honeybees involving secretion and ingestion of worker and royal jellies. However, how the jelly facilitates transfer of RNA is still unknown. Here, we show that worker and royal jellies harbor robust RNA-binding activity. We report that a highly abundant jelly component, major royal jelly protein 3 (MRJP-3), acts as an extracellular non-sequence-specific RNA-aggregating factor. Multivalent RNA binding stimulates higher-order assembly of MRJP-3 into extracellular ribonucleoprotein granules that protect RNA from degradation and enhance RNA bioavailability. These findings reveal that honeybees have evolved a secreted dietary RNA-binding factor to concentrate, stabilize, and share RNA among individuals. Our work identifies high-order ribonucleoprotein assemblies with functions outside cells and organisms.

Protein subcellular relocalization and function of duplicated flagellar calcium binding protein genes in honey bee trypanosomatid parasite.

The honey bee trypanosomatid parasite, Lotmaria passim, contains two genes that encode the flagellar calcium binding protein (FCaBP) through tandem duplication in its genome. FCaBPs localize in the flagellum and entire body membrane of L. passim through specific N-terminal sorting sequences. This finding suggests that this is an example of protein subcellular relocalization resulting from gene duplication, altering the intracellular localization of FCaBP. However, this phenomenon may not have occurred in Leishmania, as one or both of the duplicated genes have become pseudogenes. Multiple copies of the FCaBP gene are present in several Trypanosoma species and Leptomonas pyrrhocoris, indicating rapid evolution of this gene in trypanosomatid parasites. The N-terminal flagellar sorting sequence of L. passim FCaBP1 is in close proximity to the BBSome complex, while that of Trypanosoma brucei FCaBP does not direct GFP to the flagellum in L. passim. Deletion of the two FCaBP genes in L. passim affected growth and impaired flagellar morphogenesis and motility, but it did not impact host infection. Therefore, FCaBP represents a duplicated gene with a rapid evolutionary history that is essential for flagellar structure and function in a trypanosomatid parasite.

Multisite imaging of neural activity using a genetically encoded calcium sensor in the honey bee.

Understanding of the neural bases for complex behaviors in Hymenoptera insect species has been limited by a lack of tools that allow measuring neuronal activity simultaneously in different brain regions. Here, we developed the first pan-neuronal genetic driver in a Hymenopteran model organism, the honey bee, and expressed the calcium indicator GCaMP6f under the control of the honey bee synapsin promoter. We show that GCaMP6f is widely expressed in the honey bee brain, allowing to record neural activity from multiple brain regions. To assess the power of this tool, we focused on the olfactory system, recording simultaneous responses from the antennal lobe, and from the more poorly investigated lateral horn (LH) and mushroom body (MB) calyces. Neural responses to 16 distinct odorants demonstrate that odorant quality (chemical structure) and quantity are faithfully encoded in the honey bee antennal lobe. In contrast, odor coding in the LH departs from this simple physico-chemical coding, supporting the role of this structure in coding the biological value of odorants. We further demonstrate robust neural responses to several bee pheromone odorants, key drivers of social behavior, in the LH. Combined, these brain recordings represent the first use of a neurogenetic tool for recording large-scale neural activity in a eusocial insect and will be of utility in assessing the neural underpinnings of olfactory and other sensory modalities and of social behaviors and cognitive abilities.

A few essential genetic loci distinguish Penstemon species with flowers adapted to pollination by bees or hummingbirds.

In the formation of species, adaptation by natural selection generates distinct combinations of traits that function well together. The maintenance of adaptive trait combinations in the face of gene flow depends on the strength and nature of selection acting on the underlying genetic loci. Floral pollination syndromes exemplify the evolution of trait combinations adaptive for particular pollinators. The North American wildflower genus Penstemon displays remarkable floral syndrome convergence, with at least 20 separate lineages that have evolved from ancestral bee pollination syndrome (wide blue-purple flowers that present a landing platform for bees and small amounts of nectar) to hummingbird pollination syndrome (bright red narrowly tubular flowers offering copious nectar). Related taxa that differ in floral syndrome offer an attractive opportunity to examine the genomic basis of complex trait divergence. In this study, we characterized genomic divergence among 229 individuals from a Penstemon species complex that includes both bee and hummingbird floral syndromes. Field plants are easily classified into species based on phenotypic differences and hybrids displaying intermediate floral syndromes are rare. Despite unambiguous phenotypic differences, genome-wide differentiation between species is minimal. Hummingbird-adapted populations are more genetically similar to nearby bee-adapted populations than to geographically distant hummingbird-adapted populations, in terms of genome-wide dXY. However, a small number of genetic loci are strongly differentiated between species. These approximately 20 "species-diagnostic loci," which appear to have nearly fixed differences between pollination syndromes, are sprinkled throughout the genome in high recombination regions. Several map closely to previously established floral trait quantitative trait loci (QTLs). The striking difference between the diagnostic loci and the genome as whole suggests strong selection to maintain distinct combinations of traits, but with sufficient gene flow to homogenize the genomic background. A surprisingly small number of alleles confer phenotypic differences that form the basis of species identity in this species complex.

Improving bee health through genomics.

Declines in bee populations across the world threaten food security and ecosystem function. It is currently not possible to routinely predict which specific stressors lead to declines in different populations or contexts, hindering efforts to improve bee health. Genomics has the potential to dramatically improve our ability to identify, monitor and predict the effects of stressors, as well as to mitigate their impacts through the use of marker-assisted selection, RNA interference and potentially gene editing. Here we discuss the most compelling recent applications of genomics to investigate the mechanisms underpinning bee population declines and to improve the health of both wild and managed bee populations.

A postmeiotically bifurcated roadmap of honeybee spermatogenesis marked by phylogenetically restricted genes.

Haploid males of hymenopteran species produce gametes through an abortive meiosis I followed by meiosis II that can either be symmetric or asymmetric in different species. Thus, one spermatocyte could give rise to two spermatids with either equal or unequal amounts of cytoplasm. It is currently unknown what molecular features accompany these postmeiotic sperm cells especially in species with asymmetric meiosis II such as bees. Here we present testis single-cell RNA sequencing datasets from the honeybee (Apis mellifera) drones of 3 and 14 days after emergence (3d and 14d). We show that, while 3d testes exhibit active, ongoing spermatogenesis, 14d testes only have late-stage spermatids. We identify a postmeiotic bifurcation in the transcriptional roadmap during spermatogenesis, with cells progressing toward the annotated spermatids (SPT) and small spermatids (sSPT), respectively. Despite an overall similarity in their transcriptomic profiles, sSPTs express the fewest genes and the least RNA content among all the sperm cell types. Intriguingly, sSPTs exhibit a relatively high expression level for Hymenoptera-restricted genes and a high mutation load, suggesting that the special meiosis II during spermatogenesis in the honeybee is accompanied by phylogenetically young gene activities.

The Nasonia pair-rule gene regulatory network retains its function over 300 million years of evolution.

Insect segmentation is a well-studied and tractable system with which to investigate the genetic regulation of development. Though insects segment their germband using a variety of methods, modelling work implies that a single gene regulatory network can underpin the two main types of insect segmentation. This means limited genetic changes are required to explain significant differences in segmentation mode between different insects. This idea needs to be tested in a wider variety of species, and the nature of the gene regulatory network (GRN) underlying this model has not been tested. Some insects, e.g. Nasonia vitripennis and Apis mellifera segment progressively, a pattern not examined in previous studies of this segmentation model, producing stripes at different times progressively through the embryo, but not from a segment addition zone. Here, we aim to understand the GRNs patterning Nasonia using a simulation-based approach. We found that an existing model of Drosophila segmentation ( Clark, 2017) can be used to recapitulate the progressive segmentation of Nasonia, if provided with altered inputs in the form of expression of the timer genes Nv-caudal and Nv-odd paired. We predict limited topological changes to the pair-rule network and show, by RNAi knockdown, that Nv-odd paired is required for morphological segmentation. Together this implies that very limited changes to the Drosophila network are required to simulate Nasonia segmentation, despite significant differences in segmentation modes, implying that Nasonia use a very similar version of an ancestral GRN used by Drosophila, which must therefore have been conserved for at least 300 million years.

Behavioral genetics and genomics: Mendel's peas, mice, and bees.

The question of the heritability of behavior has been of long fascination to scientists and the broader public. It is now widely accepted that most behavioral variation has a genetic component, although the degree of genetic influence differs widely across behaviors. Starting with Mendel's remarkable discovery of "inheritance factors," it has become increasingly clear that specific genetic variants that influence behavior can be identified. This goal is not without its challenges: Unlike pea morphology, most natural behavioral variation has a complex genetic architecture. However, we can now apply powerful genome-wide approaches to connect variation in DNA to variation in behavior as well as analyses of behaviorally related variation in brain gene expression, which together have provided insights into both the genetic mechanisms underlying behavior and the dynamic relationship between genes and behavior, respectively, in a wide range of species and for a diversity of behaviors. Here, we focus on two systems to illustrate both of these approaches: the genetic basis of burrowing in deer mice and transcriptomic analyses of division of labor in honey bees. Finally, we discuss the troubled relationship between the field of behavioral genetics and eugenics, which reminds us that we must be cautious about how we discuss and contextualize the connections between genes and behavior, especially in humans.

Phylogenomic and functional characterization of an evolutionary conserved cytochrome P450-based insecticide detoxification mechanism in bees.

The regulatory process for assessing the risks of pesticides to bees relies heavily on the use of the honeybee, Apis mellifera, as a model for other bee species. However, the validity of using A. mellifera as a surrogate for other Apis and non-Apis bees in pesticide risk assessment has been questioned. Related to this line of research, recent work on A. mellifera has shown that specific P450 enzymes belonging to the CYP9Q subfamily act as critically important determinants of insecticide sensitivity in this species by efficiently detoxifying certain insecticide chemotypes. However, the extent to which the presence of functional orthologs of these enzymes is conserved across the diversity of bees is unclear. Here we used a phylogenomic approach to identify > 100 putative CYP9Q functional orthologs across 75 bee species encompassing all major bee families. Functional analysis of 26 P450s from 20 representative bee species revealed that P450-mediated detoxification of certain systemic insecticides, including the neonicotinoid thiacloprid and the butenolide flupyradifurone, is conserved across all major bee pollinator families. However, our analyses also reveal that CYP9Q-related genes are not universal to all bee species, with some Megachilidae species lacking such genes. Thus, our results reveal an evolutionary conserved capacity to metabolize certain insecticides across all major bee families while identifying a small number of bee species where this function may have been lost. Furthermore, they illustrate the potential of a toxicogenomic approach to inform pesticide risk assessment for nonmanaged bee species by predicting the capability of bee pollinator species to break down synthetic insecticides.

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.

Insights from Melipona bicolor hybrid genome assembly: a stingless bee genome with chromosome-level scaffold.

BACKGROUND: The highly eusocial stingless bees are crucial pollinators of native and agricultural ecosystems. Nevertheless, genomic studies within this bee tribe remain scarce. We present the genome assembly of the stingless bee Melipona bicolor. This bee is a remarkable exception to the typical single-queen colony structure, since in this species, multiple queens may coexist and share reproductive duties, resulting in genetically diverse colonies with weak kinship connections. As the only known genuinely polygynous bee, M. bicolor's genome provides a valuable resource for investigating sociality beyond kin selection. RESULTS: The genome was assembled employing a hybrid approach combining short and long reads, resulting in 241 contigs spanning 259 Mb (N50 of 6.2 Mb and 97.5% complete BUSCOs). Comparative analyses shed light on some evolutionary aspects of stingless bee genomics, including multiple chromosomal rearrangements in Melipona. Additionally, we explored the evolution of venom genes in M. bicolor and other stingless bees, revealing that, apart from two genes, the conserved repertoire of venom components remains under purifying selection in this clade. CONCLUSION: This study advances our understanding of stingless bee genomics, contributing to the conservation efforts of these vital pollinators and offering insights into the evolutionary mechanisms driving their unique adaptations.

Genomic diversity and population structure of Carniolan honey bee in its native habitat.

BACKGROUND: Research into the genetic diversity of honey bee (Apis mellifera L.) populations has become increasingly significant in recent decades, primarily due to population declines attributed to human activities and climate change. As a species of great importance, breeding programs that leverage understanding of genomic diversity could offer solutions to mitigate these challenges. The objective of this study was to examine the genomic diversity and population structure of Carniolan honey bees (Apis mellifera carnica) using the Illumina SNP chip on a large honey bee sample collected from Central and South-Eastern European countries. The study also aims to offer recommendations for future breeding programs. RESULTS: Our analysis involved Discriminant Analysis of Principal Components (DAPC), heterozygosity, admixture analysis, fixation indices (FST), Neighbour-Joining tree, gene flow and Isolation-by-distance analysis. DAPC indicated distinct separation between the Carniolan and Italian honey bee (Apis mellifera ligustica) populations, whereas the admixture analysis revealed varying levels of gene flow and genetic admixture within the Carniolan honey bee populations, demonstrating closer relationships between specific geographic regions (confirmed by Isolation-by-distance analysis). Furthermore, the research of heterozygosity, genomic inbreeding, pairwise FST values, and Neighbour-Joining tree provided insights into the patterns of genetic differentiation and similarity among the populations of Carniolan honey bee within its natural habitat. We have observed genetic homogeneity of the Carniolan honey bee population when considered in a broader genetic/geographical context. However, the Carniolan honey bee has sufficient genetic diversity in its geographical home range that needs to be carefully monitored and maintained. CONCLUSIONS: This study provides important insights into the genetic composition, differentiation, and relationships among Carniolan honey bee populations in Central and South-Eastern European countries. The findings are crucial for conservation efforts, breeding programs, and sustainable beekeeping practices. They emphasise the importance of considering genetic factors and population structure in the breeding and management of honey bees. By understanding these genetic relationships, we can develop strategies to preserve genetic diversity, improve breeding outcomes, and ensure the resilience of honey bee populations in the face of environmental changes and challenges. This knowledge can also inform policy makers and stakeholders on best practices to maintain healthy bee populations, which are vital for ecosystem services and agricultural productivity.

The nuclear and mitochondrial genome assemblies of Tetragonisca angustula (Apidae: Meliponini), a tiny yet remarkable pollinator in the Neotropics.

BACKGROUND: The field of bee genomics has considerably advanced in recent years, however, the most diverse group of honey producers on the planet, the stingless bees, are still largely neglected. In fact, only eleven of the ~ 600 described stingless bee species have been sequenced, and only three using a long-read (LR) sequencing technology. Here, we sequenced the nuclear and mitochondrial genomes of the most common, widespread and broadly reared stingless bee in Brazil and other neotropical countries-Tetragonisca angustula (popularly known in Brazil as jataí). RESULTS: A total of 48.01 Gb of DNA data were generated, including 2.31 Gb of Pacific Bioscience HiFi reads and 45.70 Gb of Illumina short reads (SRs). Our preferred assembly comprised 683 contigs encompassing 284.49 Mb, 62.84 Mb of which (22.09%) corresponded to 445,793 repetitive elements. N50, L50 and complete BUSCOs reached 1.02 Mb, 91 contigs and 97.1%, respectively. We predicted that the genome of T. angustula comprises 17,459 protein-coding genes and 4,108 non-coding RNAs. The mitogenome consisted of 17,410 bp, and all 37 genes were found to be on the positive strand, an unusual feature among bees. A phylogenomic analysis of 26 hymenopteran species revealed that six odorant receptor orthogroups of T. angustula were found to be experiencing rapid evolution, four of them undergoing significant contractions. CONCLUSIONS: Here, we provided the first nuclear and mitochondrial genome assemblies for the ecologically and economically important T. angustula, the fourth stingless bee species to be sequenced with LR technology thus far. We demonstrated that even relatively small amounts of LR data in combination with sufficient SR data can yield high-quality genome assemblies for bees.

Differences between uncapping and removal behaviors in Apis cerana from the perspective of long non-coding RNAs.

BACKGROUND: Hygienic behavior, a specialized form of immune response evolved in social insects, plays a crucial role in safeguarding colonies from disease spread. In honeybee colonies, such behavior typically entails the dual steps of uncapping and removal of unhealthy and deceased brood. Although in recent years, numerous studies have examined the development of hygienic behavior, the mechanisms underlying the division in the performance of uncapping and removal have yet to be sufficiently elucidated. In this regard, long non-coding RNAs (lncRNAs) have been evidenced to be engaged in regulating the physiological activities of honeybees; however, whether lncRNAs are likewise involved in the uncapping and removal tasks has not been clarified. RESULTS: In this study, the strong hygienic Apis cerana worker bees were used and the processes of uncapping and removal behaviors in three colonies were assayed with freeze-killed brood in the field. We then sequenced the antennal RNAs of honeybees to identify differentially expressed lncRNAs and performed lncRNA-mRNA association analysis to establish the differences between uncapping and removal. We detected 1,323 differentially expressed lncRNAs in the antennae, and the findings of lncRNA-mRNA association analyses revealed that the target genes of differentially expressed lncRNAs between uncapping and removal worker bees were predominantly linked to response to stimulus, receptor activity, and synapse. Notably, among the lncRNAs enriched in cellular response to stimulus, XR_001766094.2 was exclusively expressed in the uncapping worker bees. Based on these findings, we hypothesize that XR_001766094.2 plays a key role in distinguishing uncapping from removal behaviors by responding to external stimulus, thereby suggesting that the division of hygienic behaviors is governed by differential thresholds of responsiveness to environmental cues. CONCLUSION: We characterized differences in the uncapping and removal behaviors of worker bees from a perspective of lncRNAs. Uncapping bees may be equipped with a more rapid stimulatory response and more acute olfactory sensitivity, contributing to the rapid hygienic behavior in honeybee colonies. Our results thus establish a foundation for potential lncRNA-mediated gene expression regulation in hygienic behavior.

Regulatory roles of long non-coding RNAs in short-term heat stress in adult worker bees.

Long non-coding RNAs (lncRNAs) are crucial modulators of post-transcriptional gene expression regulation, cell fate determination, and disease development. However, lncRNA functions during short-term heat stress in adult worker bees are poorly understood. Here, we performed deep sequencing and bioinformatic analyses of honeybee lncRNAs. RNA interference was performed by using siRNA targeting the most highly expressed lncRNA. The silencing effect on lncRNA and the relative expression levels of seven heat shock protein (HSP) genes, were subsequently examined. Overall, 7,842 lncRNAs and 115 differentially expressed lncRNAs (DELs) were identified in adult worker bees following heat stress exposure. Structural analysis revealed that the overall expression abundance, length of transcripts, exon number, and open reading frames of lncRNAs were lower than those of mRNAs. GO analysis revealed that the target genes were mainly involved in "metabolism," "protein folding," "response to stress," and "signal transduction" pathways. KEGG analysis indicated that the "protein processing in endoplasmic reticulum" and "longevity regulating pathway-multiple species" pathways were most enriched. Quantitative real-time polymerase chain reaction (qRT-PCR) detection of the selected DELs confirmed the reliability of the sequencing data. Moreover, the siRNA experiment indicated that feeding siRNA yielded a silencing efficiency of 77.51% for lncRNA MSTRG.9645.5. Upon silencing this lncRNA, the expression levels of three HSP genes were significantly downregulated (p < 0.05), whereas those of three other HSP genes were significantly upregulated (p < 0.05). Our results provide a new perspective for understanding the regulatory mechanisms of lncRNAs in adult worker bees under short-term heat stress.

Post-transcriptional regulation of Dufour's gland reproductive signals in bumble bees.

Pheromone communication is a key mechanism by which the reproductive division of labor is maintained within insect communities. Understanding how pheromones evolved to regulate social behavior requires knowledge of the molecular regulation of their production. However, even in cases where pheromones were identified, our understanding of their biosynthesis and molecular regulation remains limited. Bumble bees provide a unique system to explore pheromone biosynthesis since workers produce ester sterility signals in their Dufour's gland that differ from gyne-specific esters and are not produced by queens. These esters are hypothesized to be produced in the exocrine gland where they are stored, and indeed queens, gynes and workers differ significantly in the expression of Dufour's gland genes coding to enzymes involved in the biosynthesis of esters. However, a previous transcriptome analysis revealed no gene expression differences in the Dufour's gland of workers despite differences in both ester production and ovarian activation, suggesting that ester production may be regulated lower down. Proteomics of the Dufour's gland of queens, gynes, and workers recovered over 2400 proteins and broadly matched the previous RNAseq data. However, more than 100 differentially expressed proteins were found between the worker groups, including key enzymes in fatty acid biosynthesis, indicating that the regulation of reproductive signal biosynthesis in workers is done post-transcription. Overall, our data provide evidence that pheromone biosynthesis in the Dufour's gland is caste specific, that gynes and workers are likely using different enzymes to make their respective wax esters, and that the regulation on pheromone production in queens, gynes and workers is likely done at different regulatory levels, with workers signals being subjected to regulation at the protein level.

Genetic Modification of a Hox Locus Drives Mimetic Color Pattern Variation in a Highly Polymorphic Bumble Bee.

Müllerian mimicry provides natural replicates ideal for exploring mechanisms underlying adaptive phenotypic divergence and convergence, yet the genetic mechanisms underlying mimetic variation remain largely unknown. The current study investigates the genetic basis of mimetic color pattern variation in a highly polymorphic bumble bee, Bombus breviceps (Hymenoptera, Apidae). In South Asia, this species and multiple comimetic species converge onto local Müllerian mimicry patterns by shifting the abdominal setal color from orange to black. Genetic crossing between the orange and black phenotypes suggested the color dimorphism being controlled by a single Mendelian locus, with the orange allele being dominant over black. Genome-wide association suggests that a locus at the intergenic region between 2 abdominal fate-determining Hox genes, abd-A and Abd-B, is associated with the color change. This locus is therefore in the same intergenic region but not the same exact locus as found to drive red black midabdominal variation in a distantly related bumble bee species, Bombus melanopygus. Gene expression analysis and RNA interferences suggest that differential expression of an intergenic long noncoding RNA between abd-A and Abd-B at the onset setal color differentiation may drive the orange black color variation by causing a homeotic shift late in development. Analysis of this same color locus in comimetic species reveals no sequence association with the same color shift, suggesting that mimetic convergence is achieved through distinct genetic routes. Our study establishes Hox regions as genomic hotspots for color pattern evolution in bumble bees and demonstrates how pleiotropic developmental loci can drive adaptive radiations in nature.

Putative Signals of Generalist Plant Species Adaptation to Local Pollinator Communities and Abiotic Factors.

The reproductive success of flowering plants with generalized pollination systems is influenced by interactions with a diverse pollinator community and abiotic factors. However, knowledge about the adaptative potential of plants to complex ecological networks and the underlying genetic mechanisms is still limited. Based on a pool-sequencing approach of 21 natural populations of Brassica incana in Southern Italy, we combined a genome-environmental association analysis with a genome scan for signals of population genomic differentiation to discover genetic variants associated with the ecological variation. We identified genomic regions putatively involved in the adaptation of B. incana to the identity of local pollinator functional categories and pollinator community composition. Interestingly, we observed several shared candidate genes associated with long-tongue bees, soil texture, and temperature variation. We established a genomic map of potential generalist flowering plant local adaptation to complex biotic interactions, and the importance of considering multiple environmental factors to describe the adaptive landscape of plant populations.

Enhanced sensitivity of chimeric insect olfactory co-receptors for detecting odorant molecules.

Insect olfactory receptors (ORs) are seven-transmembrane domain ion channels that function by forming heteromeric complexes with olfactory receptor co-receptors (Orcos). In this study, we investigated the potential for enhancing sensitivity of odor detection and responsivity through genetic modification of Orcos, considering its wider application in odor sensing. First, we measured the intensity of response to 1-octen-3-ol for the mosquito Aedes aegypti OR (AaOR8) when complexed individually with an Orco from the same mosquito (AaOrco), the honeybee Apis mellifera (AmOrco), the silkworm Bombyx mori (BmOrco), or the fruit fly Drosophila melanogaster (DmOrco). Relative to the other Orcos, AmOrco demonstrated higher sensitivity and responsivity, with a 1.8 to 21-fold decrease in the half-maximal effective concentration (EC50) and a 1.6-8.8-fold increase in the maximal effect (Emax), respectively. Furthermore, AmOrco co-expressed with AaOR10, BmOR56, or DmOR47a showed higher sensitivity and responsivity than AaOrco, BmOrco, or DmOrco co-expressed with their respective ORs. To further increase sensitivity and responsivity, we engineered chimeric Orcos by fusing AmOrco with DmOrco, considering the domain characteristics of Orcos. The response to 1-octen-3-ol was evaluated for AaOR8 when complexed individually with AmOrco, as well as for a mutant that combines DmOrco from the N-terminal (NT) to the C-terminal region of the fourth transmembrane domain (TM4) with the region of AmOrco following TM4 (Dm[NT-TM4]AmOrco). When compared to AmOrco, Dm(NT-TM4)AmOrco showed higher sensitivity and responsivity, with a 1.4-fold decrease in the EC50 and a 1.4-fold increase in the Emax, respectively. In addition, Dm(NT-TM4)AmOrco co-expressed with either DmOR47a or BmOR56 demonstrated higher sensitivity and responsivity than AmOrco co-expressed with their respective ORs. These results suggest that AmOrco could be a relatively more sensitive Orco, and further enhancement of sensitivity and responsivity could be achieved through recombination with heterologous Orcos near the TM4 of AmOrco.