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

Surviving high temperatures: The crucial role of vesicular inhibitory amino acid transporter in Asian honeybee, Apis cerana

Asian honeybees (Apis cerana) play a crucial role as pollinators to service for the ecological stability. However, their proliferation and growth are significantly impacted by environmental temperature stress. This study delves into the function of the Apis cerana vesicular inhibitory amino acid transporter gene (AcVIAAT) in safeguarding Asian honeybees against high-temperature stress. The AcVIAAT gene exhibits positive responsiveness in honeybees subjected to varying thermal conditions by triggering the genes associated with oxidative stress. Molecular docking, co-immunoprecipitation, and ELISA verify the capacity of the AcVIAAT protein to interact with γ-aminobutyric acid (GABA), a key inhibitory neurotransmitter. Administering GABA to honeybees significantly improves their survival rate under high-temperature stress and also simultaneously upregulating oxidative stress-related genes. Therefore, these findings reveal that the AcVIAAT gene enhances the thermoregulatory capacity of honeybees by modulating oxidative stress-related genes through facilitating GABA transport. The characterization of six non-synonymous SNPs in the AcVIAAT gene among A.cerana populations distributed across both the northern and southern regions indicates a potential association between gene variation and environmental adaptation. Our results contribute to elucidating the molecular mechanisms underlying high-temperature tolerance in Asian honeybees and provide a promising genetic marker for enhancing heat tolerance through genetic improvement.

Reminiscences on the honeybee genome project and the rise of epigenetic concepts in insect science

The sequencing of the honeybee genome in 2006 was an important technological and logistic achievement experience. But what benefits have flown from the honeybee genome project? What does the annotated genomic assembly mean for the study of behavioural complexity and organismal function in honeybees? Here, I discuss several lines of research that have arisen from this project and highlight the rapidly expanding studies on insect epigenomics, emergent properties of royal jelly, the mechanism of nutritional control of development and the contribution of epigenomic regulation to the evolution of sociality. I also argue that the term 'insect epigenetics' needs to be carefully redefined to reflect the diversity of epigenomic toolkits in insects and the impact of lineage-specific innovations on organismal outcomes. The honeybee genome project helped pioneer advances in social insect molecular biology, and fuelled breakthrough research into the role of flexible epigenomic control systems in linking genotype to phenotype.

Synergistic effects of Pleistocene geological and climatic events on complex phylogeographic history of widespread sympatric species of Megaloptera in East Asia

Unraveling the phylogeographic histories of species remains a key endeavor for comprehending the evolutionary processes contributing to the rich biodiversity and high endemism found in East Asia. In this study, we explored the phylogeographic patterns and demographic histories of three endemic fishfly and dobsonfly species ( Neochauliodes formosanus, Protohermes costalis, and Neoneuromus orientalis) belonging to the holometabolan order Megaloptera. These species, which share a broad and largely overlapping distribution, were analyzed using comprehensive mitogenomic data. Our findings revealed a consistent influence of vicariance on the population isolation of Neoc. formosanus and P. costalis between Hainan, Taiwan, and the East Asian mainland during the early Pleistocene, potentially hindering subsequent colonization of the later diverged Neon. orientalis to these islands. Additionally, we unveiled the dual function of the major mountain ranges in East Asia, serving both as barriers and conduits, in shaping the population structure of all three species. Notably, we demonstrated that these co-distributed species originated from Southwest, Southern, and eastern Central China, respectively, then subsequently migrated along multi-directional routes, leading to their sympatric distribution on the East Asian mainland. Furthermore, our results highlighted the significance of Pleistocene land bridges along the eastern coast of East Asia in facilitating the dispersal of mountain-dwelling insects with low dispersal ability. Overall, this study provides novel insight into the synergistic impact of Pleistocene geological and climatic events in shaping the diversity and distribution of aquatic insects in East Asia.

Chromosome-level genome assembly of Chinese water Scorpion Ranatra chinensis (Heteroptera: Nepidae)

Heteroptera (the true bugs), one of the most diverse lineages of insects, diversified in feeding strategies and living habitats, and thus become an ideal lineage for studies on adaptive evolution. Chinese water scorpion Ranatra chinensis (Heteroptera: Nepidae) is a predaceous bug living in lentic water systems, representing an ideal model for studying habitat transition and adaptation to water environment. However, genetic studies on this water bug remain limited. Here, we obtained a chromosome-level genome of R. chinensis using PacBio HiFi long reads and Hi-C sequencing reads. The total assembly size of genome is 867.89 Mb, with a scaffold N50 length of 26.48 Mb and the GC content of 39.50%. All contigs were assembled into 23 pseudo-chromosomes (N = 19 A + X1X2X3X4), and we predicted 18,424 protein-coding genes in this genome. This study will provide valuable genomic resources for future studies on the biology, water adaptation, and genome evolution of water bugs.

Lactobacillus Firm-5-derived succinate prevents honeybees from having diabetes-like symptoms

The gut microbiome plays an important role in honeybee hormonal regulation and growth, but the underlying mechanisms are poorly understood. Here, we showed that the depletion of gut bacteria resulted in reduced expression of insulin-like peptide gene (ilp) in the head, accompanied by metabolic syndromes resembling those of Type 1 diabetes in humans: hyperglycemia, impaired lipid storage, and decreased metabolism. These symptoms were alleviated by gut bacterial inoculation. Gut metabolite profiling revealed that succinate, produced by Lactobacillus Firm-5, played deterministic roles in activating ilp gene expression and in regulating metabolism in honeybees. Notably, we demonstrated that succinate modulates host ilp gene expression through stimulating gut gluconeogenesis, a mechanism resembling that of humans. This study presents evidence for the role of gut metabolite in modulating host metabolism and contributes to the understanding of the interactions between gut microbiome and bee hosts.

Tracing the genealogy origin of geographic populations based on genomic variation and deep learning

Assigning a query individual animal or plant to its derived population is a prime task in diverse applications related to organismal genealogy. Such endeavors have conventionally relied on short DNA sequences under a phylogenetic framework. These methods naturally show constraints when the inferred population sources are ambiguously phylogenetically structured, a scenario demanding substantially more informative genetic signals. Recent advances in cost-effective production of whole-genome sequences and artificial intelligence have created an unprecedented opportunity to trace the population origin for essentially any given individual, as long as the genome reference data are comprehensive and standardized. Here, we developed a convolutional neural network method to identify population origins using genomic SNPs. Three empirical datasets (an Asian honeybee, a red fire ant, and a chicken datasets) and two simulated populations are used for the proof of concepts. The performance tests indicate that our method can accurately identify the genealogy origin of query individuals, with success rates ranging from  93 % to 100 %. We further showed that the accuracy of the model can be significantly increased by refining the informative sites through FST filtering. Our method is robust to configurations related to batch sizes and epochs, whereas model learning benefits from the setting of a proper preset learning rate. Moreover, we explained the importance score of key sites for algorithm interpretability and credibility, which has been largely ignored. We anticipate that by coupling genomics and deep learning, our method will see broad potential in conservation and management applications that involve natural resources, invasive pests and weeds, and illegal trades of wildlife products.

Multiple stresses induced by chronic exposure to flupyradifurone affect honey bee physiological states

Flupyradifurone (FPF) has been reported to have a potential risk to terrestrial and aquatic ecosystems. In the present study, the effects of chronic FPF exposure on bees were systematically investigated at the individual behavioral, tissue, cell, enzyme activity, and the gene expression levels. Chronic exposure (14 d) to FPF led to reduced survival (12 mg/L), body weight gain (4 and 12 mg/L), and food utilization efficiency (4 and 12 mg/L). Additionally, FPF exposure (12 mg/L) impaired sucrose sensitivity and memory of bees. Morphological analysis revealed significant cellular and subcellular changes in brain neurons and midgut epithelial cells, including mitochondrial damage, nuclear disintegration, and apoptosis. FPF exposure (4 and 12 mg/L) led to oxidative stress, as evidenced by increased lipid peroxidation and alterations in antioxidant enzyme activity. Notably, gene expression analysis indicated significant dysregulation of apoptosis, immune, detoxification, sucrose responsiveness and memory-related genes, suggesting the involvement of different pathways in FPF-induced toxicity. The multiple stresses and potential mechanisms described here provide a basis for determining the intrinsic toxicity of FPF.

Drivers of genomic differentiation landscapes in populations of disparate ecological and geographical settings within mainland Apis cerana

Elucidating the evolutionary processes that drive population divergence can enhance our understanding of the early stages of speciation and inform conservation management decisions. The honeybee Apis cerana displays extensive population divergence, providing an informative natural system for exploring these processes. The mainland lineage A. cerana includes several peripheral subspecies with disparate ecological and geographical settings radiated from a central ancestor. Under this evolutionary framework, we can explore the patterns of genome differentiation and the evolutionary models that explain them. We can also elucidate the contribution of non-genomic spatiotemporal mechanisms (extrinsic features) and genomic mechanisms (intrinsic features) that influence these genomic differentiation landscapes. Based on 293 whole genomes, a small part of the genome is highly differentiated between central-peripheral subspecies pairs, while low and partial parallelism partly reflects idiosyncratic responses to environmental differences. Combined elements of recurrent selection and speciation-with-gene-flow models generate the heterogeneous genome landscapes. These elements weight differently between central-island and other central-peripheral subspecies pairs, influenced by glacial cycles superimposed on different geomorphologies. Although local recombination rates exert a significant influence on patterns of genomic differentiation, it is unlikely that low-recombination rates regions were generated by structural variation. In conclusion, complex factors including geographical isolation, divergent ecological selection and non-uniform genome features have acted concertedly in the evolution of reproductive barriers that could reduce gene flow in part of the genome and facilitate the persistence of distinct populations within mainland lineage of A. cerana.

Tropical and montane Apis cerana show distinct dance-distance calibration curves

Social bees have evolved sophisticated communication systems to recruit nestmates to newly found food sources. As foraging ranges can vary from a few hundred meters to several kilometers depending on the environment or season, populations of social bee species living in different climate zones likely show specific adaptations in their recruitment communication. Accordingly, studies in the western honey bee, Apis mellifera, demonstrated that temperate populations exhibit shallower dance-calibration curves compared with tropical populations. Here, we report the first comparison of calibration curves for three Indian Apis cerana lineages: the tropical Apis indica, and the two montane Himalayan populations Apis cerana cerana (Himachal Pradesh) and Apis cerana kashmirensis (Jammu and Kashmir). We found that the colonies of the two montane A. cerana populations show dance-distance calibration curves with significantly shallower slopes than those of the tropical A. indica. Next, we transferred A. c. cerana colonies to Bangalore (∼ 2600 km away) to obtain calibration curves in the same location as A. indica. The common garden experiment confirmed this difference in slopes, implying that the lineages exhibit genetically fixed differences in dance-distance coding. However, the slopes of the calibration curves of the transferred A. c. cerana colonies were also significantly higher than those of the colonies tested in their original habitat, indicating an important effect of the environment. The differences in dance-distance coding between temperate and tropical A. cerana lineages resemble those described for Apis mellifera, suggesting that populations of both species independently evolved similar adaptations.

Does coevolution in refugia drive mimicry in bumble bees? Insights from a South Asian mimicry group

Müllerian mimicry was proposed to be an example of a coevolved mutualism promoted by population isolation in glacial refugia. This, however, has not been well supported in butterfly models. Here, we use genomic data to test this theory while examining the population genetics behind mimetic diversification in a pair of co-mimetic bumble bees, Bombus breviceps Smith and Bombus trifasciatus Smith. In both lineages, populations were structured by geography but not as much by color pattern, suggesting sharing of color alleles across regions of restricted gene flow and formation of mimicry complexes in the absence of genetic differentiation. Demographic analyses showed mismatches between historical effective population size changes and glacial cycles, and niche modeling revealed only mild habitat retraction during glaciation. Moreover, mimetic subpopulations of the same color form in the two lineages only in some cases exhibit similar population history and genetic divergence. Therefore, the current study supports a more complex history in this comimicry than a simple refugium-coevolution model.

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.

The PWWP domain and the evolution of unique DNA methylation toolkits in Hymenoptera

DNMT3 in Hymenoptera has a unique duplication of the essential PWWP domain. Using GST-tagged PWWP fusion proteins and histone arrays we show that these domains have gained new properties and represent the first case of PWWP domains binding to H3K27 chromatin modifications, including H3K27me3, a key modification that is important during development. Phylogenetic analyses of 107 genomes indicate that the duplicated PWWP domains separated into two sister clades, and their distinct binding capacities are supported by 3D modeling. Other features of this unique DNA methylation system include variable copies, losses, and duplications of DNMT1 and DNMT3, and combinatorial generations of DNMT3 isoforms including variants missing the catalytic domain. Some of these losses and duplications of are found only in parasitic wasps. We discuss our findings in the context of the crosstalk between DNA methylation and histone methylation, and the expanded potential of epigenomic modifications in Hymenoptera to drive evolutionary novelties.

Understanding natural selection and similarity: Convergent, parallel and repeated evolution

Parallel and convergent evolution offer some of the most compelling evidence for the significance of natural selection in evolution, as the emergence of similar adaptive solutions is unlikely to occur by random chance alone. However, these terms are often employed inconsistently, leading to misinterpretation and confusion, and recently proposed definitions have unintentionally diminished the emphasis on the evolution of similar adaptive solutions. Here, I examine various conceptual frameworks and definitions related to parallel and convergent evolution and propose a consolidated framework that enhances our comprehension of these evolutionary patterns. The primary aim of this framework is to harmonize the concepts of parallel and convergent evolution together with natural selection and the idea of similarity. Both concepts involve the evolution of similar adaptive solutions as a result of environmental challenges. The distinction lies in ancestral phenotypes. Parallel evolution takes place when the ancestral phenotypes (before selection) of the lineages are similar. Convergent evolution happens when the lineages have distinct ancestral phenotypes (before selection). Because an ancestral-based distinction will inevitably lead to cases where uncertainty in the distinction may arise, the framework includes a general term, repeated evolution, which can be used as a term applying to the evolution of similar phenotypes and genotypes as well as similar responses to environmental pressures. Based on the argument that genetic similarity may frequently arise without selection, the framework posits that the similarity of genetic sequences is not of great interest unless linked to the actions of natural selection or to the origins (mutation, standing genetic variation, gene flow) and locations of the similar sequences.

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.

Adaptive evolution to the natural and anthropogenic environment in a global invasive crop pest, the cotton bollworm

The cotton bollworm, Helicoverpa armigera, is set to become the most economically devastating crop pest in the world, threatening food security and biosafety as its range expands across the globe. Key to understanding the eco-evolutionary dynamics of H. armigera, and thus its management, is an understanding of population connectivity and the adaptations that allow the pest to establish in unique environments. We assembled a chromosome-scale reference genome and re-sequenced 503 individuals spanning the species range to delineate global patterns of connectivity, uncovering a previously cryptic population structure. Using a genome-wide association study (GWAS) and cell line expression of major effect loci, we show that adaptive changes in a temperature- and light-sensitive developmental pathway enable facultative diapause and that adaptation of trehalose synthesis and transport underlies cold tolerance in extreme environments. Incorporating extensive pesticide resistance monitoring, we also characterize a suite of novel pesticide and Bt resistance alleles under selection in East China. These findings offer avenues for more effective management strategies and provide insight into how insects adapt to variable climatic conditions and newly colonized environments.

Defining honeybee subspecies in an evolutionary context warrants strategized conservation

Despite the urgent need for conservation consideration, strategic action plans for the preservation of the Asian honeybee, Apis cerana Fabricius, 1793, remain lacking. Both the convergent and divergent adaptations of this widespread insect have led to confusing phenotypical traits and inconsistent infraspecific taxonomy. Unclear subspecies boundaries pose a significant challenge to honeybee conservation efforts, as it is difficult to effectively prioritize conservation targets without a clear understanding of subspecies identities. Here, we investigated genome variations in 362 worker bees representing almost all populations of mainland A. cerana to understand how evolution has shaped its population structure. Whole-genome single nucleotide polymorphisms (SNPs) based on nuclear sequences revealed eight putative subspecies, with all seven peripheral subspecies exhibiting mutually exclusive monophyly and distinct genetic divergence from the widespread central subspecies. Our results demonstrated that most classic morphological traits, including body size, were related to the climatic variables of the local habitats and did not reflect the true evolutionary history of the organism. Thus, such morphological traits were not suitable for subspecific delineation. Conversely, wing vein characters showed relative independence to the environment and supported the subspecies boundaries inferred from nuclear genomes. Mitochondrial phylogeny further indicated that the present subspecies structure was a result of multiple waves of population divergence from a common ancestor. Based on our findings, we propose that criteria for subspecies delineation should be based on evolutionary independence, trait distinction, and geographic isolation. We formally defined and described eight subspecies of mainland A. cerana. Elucidation of the evolutionary history and subspecies boundaries enables a customized conservation strategy for both widespread and endemic honeybee conservation units, guiding colony introduction and breeding.

Effects of local domestication warrant attention in honey bee population genetics

Honey bees, Apis mellifera, have for millennia been managed and exploited by humans and introduced into most suitable regions worldwide. However, given the lack of records for many introduction events, treating A. mellifera populations as native would predictably bias genetic studies regarding origin and evolution. Here, we used the Dongbei bee, a well-documented population, introduced beyond the natural distribution range approximately 100 years ago, to elucidate the effects of local domestication on animal population genetic analyses. Strong domestication pressure was detected in this population, and the genetic divergence between Dongbei bee and its ancestral subspecies was found to have occurred at the lineage level. Results of phylogenetic and time divergence analyses could consequently be misinterpreted. Proposing new subspecies or lineages and performing analyses of origin should thus strive to eliminate anthropogenic effects. We highlight the need for definitions of landrace and breed in honey bee sciences and make preliminary suggestions.

Climatic oscillation promoted diversification of spinous assassin bugs during Pleistocene glaciation

Insect speciation is among the most fascinating topics in evolutionary biology; however, its underlying mechanisms remain unclear. Allopatric speciation represents one of the major types of speciation and is believed to have frequently occurred during glaciation periods, when climatic oscillation may have caused suitable habitats to be fragmented repeatedly, creating geographical isolation among populations. However, supporting evidence for allopatric speciation of insects in East Asia during the Pleistocene glaciation remains lacking. We aim to investigate the effect of climatic oscillation during the Pleistocene glaciation on the diversification pattern and evolutionary history of hemipteran insects and to test the hypothesis of Pleistocene species stability using spinous assassin bugs Sclomina (Hemiptera: Reduviidae), a small genus widely distributed in southern China but was later found to have cryptic species diversity. Here, using the whole mitochondrial genome (mitogenome) and nuclear ribosomal RNA genes, we investigated both interspecific and intraspecific diversification patterns of spinous assassin bugs. Approximate Bayesian computation, ecological niche modeling, and demographic history analyses were also applied to understand the diversification process and driven factors. Our data suggest that the five species of Sclomina are highly diverged, despite three of them currently being cryptic. Speciation occurred during the Pleistocene when suitable distribution areas were possibly fragmented. Six phylogeographic groups in the type species S. erinacea were identified, among which two groups underwent expansion during the early Last Glacial Period and after Last Glacier Maximum. Our analyses suggest that this genus may have experienced climate-driven habitat fragmentation and postglacial expansion in the Pleistocene, promoting allopatric speciation and intraspecific diversification. Our results reveal underestimated species diversity in a small insect group and illustrate a remarkable example of allopatric speciation of insects in East Asia promoted by Pleistocene climatic oscillations. These findings provide important insights into the speciation processes and aid the conservation of insect species diversity.

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.

Conservation Genomic Analysis of the Asian Honeybee in China Reveals Climate Factors Underlying Its Population Decline

The Asian honeybee, Apis cerana, is one of the most important native pollinators in Asia. Asian honeybees were believed to be under significant decline in China based on a report in 2005. On the contrary, a recent survey revealed that Asian honeybee populations in China are stable and even slightly increased in some regions. Therefore, the declining status of A. cerana populations in China is still unclear. Taking advantage of the abundant, publicly available genomic data for Asian honeybees in China, we employed conservation genomics methods to understand if Asian honeybee populations in China are declining and what the underlying climate factors are. We reconstructed the changes of effective population size (Ne) within the recent past for 6 population groups of Asian honeybees and found out that only one of them (population in Bomi, Tibet) showed a consistently declining Ne from the last 100 generations to 25 generations. Selective sweep analysis suggests that genes related to the tolerance of low temperatures and strong ultraviolet radiation are under selection in the declining population, indicating that these two climate factors most likely underlie the decline of BM populations during the recent past. Our study provides insights into the dynamic changes of Asian honeybee populations in China and identifies climate factors that underlie its population decline, which is valuable for the conservation of this important pollinator.

Identification of the Toxic Compounds in Camellia oleifera Honey and Pollen to Honey Bees (Apis mellifera)

Identifying the components of Camellia oleifera honey and pollen and conducting corresponding toxicological tests are essential to revealing the mechanism of Camellia oleifera toxicity to honey bees. In this research, we investigated the saccharides and alkaloids in honey, nectar, and pollen from Camellia oleifera, which were compared with honey, nectar, and pollen from Brassica napus, a widely planted flowering plant. The result showed that melibiose, manninotriose, raffinose, stachyose, and lower amounts of santonin and caffeine were found in Camellia oleifera nectar, pollen, and honey but not in B. napus nectar, pollen, and honey. Toxicological experiments indicated that manninotriose, raffinose, and stachyose in Camellia oleifera honey are toxic to bees, while alkaloids in Camellia oleifera pollen are not toxic to honey bees. The toxicity mechanism of oligosaccharides revealed by temporal metabolic profiling is that oligosaccharides cannot be further digested by honey bees and thus get accumulated in honey bees, disturbing the synthesis and metabolism of trehalose, ultimately causing honey bee mortality.

Genomic basis and phenotypic manifestation of (non-)parallel serpentine adaptation in Arabidopsis arenosa

Parallel evolution is common in nature and provides one of the most compelling examples of rapid environmental adaptation. In contrast to the recent burst of studies addressing genomic basis of parallel evolution, integrative studies linking genomic and phenotypic parallelism are scarce. Edaphic islands of toxic serpentine soils provide ideal systems for studying rapid parallel adaptation in plants, imposing strong, spatially replicated selection on recently diverged populations. We leveraged threefold independent serpentine adaptation of Arabidopsis arenosa and combined reciprocal transplants, ion uptake phenotyping, and available genome-wide polymorphisms to test if parallelism is manifested to a similar extent at both genomic and phenotypic levels. We found pervasive phenotypic parallelism in functional traits yet with varying magnitude of fitness differences that was congruent with neutral genetic differentiation between populations. Limited costs of serpentine adaptation suggest absence of soil-driven trade-offs. On the other hand, the genomic parallelism at the gene level was significant, although relatively minor. Therefore, the similarly modified phenotypes, for example, of ion uptake arose possibly by selection on different loci in similar functional pathways. In summary, we bring evidence for the important role of genetic redundancy in rapid adaptation involving traits with polygenic architecture.

Significant compositional and functional variation reveals the patterns of gut microbiota evolution among the widespread Asian honeybee populations

The gut microbiome is a crucial element that facilitates a host’s adaptation to a changing environment. Compared to the western honeybee Apis mellifera, the Asian honeybee, Apis cerana populations across its natural range remain mostly semi-feral and are less affected by bee management, which provides a good system to investigate how gut microbiota evolve under environmental heterogeneity on large geographic scales. We compared and analyzed the gut microbiomes of 99 Asian honeybees, from genetically diverged populations covering 13 provinces across China. Bacterial composition varied significantly across populations at phylotype, sequence-discrete population (SDP), and strain levels, but with extensive overlaps, indicating that the diversity of microbial community among A. cerana populations is driven by nestedness. Pollen diets were significantly correlated with both the composition and function of the gut microbiome. Core bacteria, Gilliamella and Lactobacillus Firm-5, showed antagonistic turnovers and contributed to the enrichment in carbohydrate transport and metabolism. By feeding and inoculation bioassays, we confirmed that the variations in pollen polysaccharide composition contributed to the trade-off of these core bacteria. Progressive change, i.e., nestedness, is the foundation of gut microbiome evolution among the Asian honeybee. Such a transition during the co-diversification of gut microbiomes is affected by environmental factors, diets in general, and pollen polysaccharides in particular.

Population Structure and Genetic Diversity of Chinese Honeybee (Apis Cerana Cerana) in Central China

Central China has a rich terrain with a temperate monsoon climate and varied natural environments for the Chinese honeybee (Apis cerana cerana). However, little comprehensive research on population genetic diversity has been done in this area. A population survey of the structure and genetic diversity of Apis cerana cerana in this area is deeply needed for understanding adaptation to variable environments and providing more references for the protection of honeybee biodiversity. In this study, we present a dataset of 72 populations of Chinese honeybees collected from nine sites by whole genome sequencing in Central China. We obtained 2,790,214,878 clean reads with an average covering a depth of 22×. A total of 27,361,052 single nucleotide polymorphisms (SNPs) were obtained by mapping to the reference genome with an average mapping rate of 93.03%. Genetic evolution analysis was presented via the population structure and genetic diversity based on the datasets of SNPs. It showed that Apis cerana cerana in plains exhibited higher genetic diversity than in mountain areas. The mantel test between Apis cerana cerana groups revealed that some physical obstacles, especially the overurbanization of the plains, contributed to the differentiation. This study is conducive to elucidating the evolution of Apis cerana in different environments and provides a theoretical basis for investigating and protecting the Chinese honeybee.

Geographical resistome profiling in the honeybee microbiome reveals resistance gene transfer conferred by mobilizable plasmids

BACKGROUND

The spread of antibiotic resistance genes (ARGs) has been of global concern as one of the greatest environmental threats. The gut microbiome of animals has been found to be a large reservoir of ARGs, which is also an indicator of the environmental antibiotic spectrum. The conserved microbiota makes the honeybee a tractable and confined ecosystem for studying the maintenance and transfer of ARGs across gut bacteria. Although it has been found that honeybee gut bacteria harbor diverse sets of ARGs, the influences of environmental variables and the mechanism driving their distribution remain unclear.

RESULTS

We characterized the gut resistome of two closely related honeybee species, Apis cerana and Apis mellifera, domesticated in 14 geographic locations across China. The composition of the ARGs was more associated with host species rather than with geographical distribution, and A. mellifera had a higher content of ARGs in the gut. There was a moderate geographic pattern of resistome distribution, and several core ARG groups were found to be prevalent among A. cerana samples. These shared genes were mainly carried by the honeybee-specific gut members Gilliamella and Snodgrassella. Transferrable ARGs were frequently detected in honeybee guts, and the load was much higher in A. mellifera samples. Genomic loci of the bee gut symbionts containing a streptomycin resistance gene cluster were nearly identical to those of the broad-host-range IncQ plasmid, a proficient DNA delivery system in the environment. By in vitro conjugation experiments, we confirmed that the mobilizable plasmids could be transferred between honeybee gut symbionts by conjugation. Moreover, "satellite plasmids" with fragmented genes were identified in the integrated regions of different symbionts from multiple areas.

CONCLUSIONS

Our study illustrates that the gut microbiota of different honeybee hosts varied in their antibiotic resistance structure, highlighting the role of the bee microbiome as a potential bioindicator and disseminator of antibiotic resistance. The difference in domestication history is highly influential in the structuring of the bee gut resistome. Notably, the evolution of plasmid-mediated antibiotic resistance is likely to promote the probability of its persistence and dissemination. Video Abstract.

Population genomics and phylogeography of Colletes gigas , a wild bee specialized on winter flowering plants

Diet specialization may affect the population genetic structure of pollinators by reducing gene flow and driving genetic differentiation, especially in pollen‐specialist bees. Colletes gigas is a pollen‐specialist pollinator of Camellia oleifera, one of the most important staple oil crops in China. Ca. oleifera blooms in cold climates and contains special compounds that make it an unusable pollen source to other pollinators. Thus, C. gigas undoubtedly plays a key role as the main pollinator of Ca. oleifera, with biological and economic significance. Here, we use a population genomic approach to analyze the roles of geography and climate on the genetic structure, genetic diversity, and demographic history of C. gigas. A total of 1,035,407 SNPs were identified from a 582.77 Gb dataset. Clustering and phylogenetic analyses revealed a marked genetic structure, with individuals grouped into nine local clusters. A significant isolation by distance was detected by both the Mantel test (R = .866, p = .008) and linear regression (R² = .616, p < .001). Precipitation and sunshine duration were positively and significantly (R ≥ .765, p ≤ .016) correlated with observed heterozygosity (H_o) and expected heterozygosity (H_e). These results showed that C. gigas populations had a distinct phylogeographic pattern determined by geographical distance and environmental factors (precipitation and sunshine duration). In addition, an analysis of paleogeographic dynamics indicated that C. gigas populations exhibited patterns of glacial expansion and interglacial contraction, likely resulting from post‐glacial habitat contraction and fragmentation. Our results indicated that the peculiar phylogeographic patterns in C. gigas populations may be related to their specialization under long‐term adaptation to host plants. This work improves our understanding of the population genetics in pollen‐specialist bees. The distinct genetic clusters identified in this study should be taken into consideration for the protection and utilization of this specialized crop pollinator.

The first steps toward a global pandemic: Reconstructing the demographic history of parasite host switches in its native range

Host switching allows parasites to expand their niches. However, successful switching may require suites of adaptations and also may decrease performance on the old host. As a result, reductions in gene flow accompany many host switches, driving speciation. Because host switches tend to be rapid, it is difficult to study them in real-time, and their demographic parameters remain poorly understood. As a result, fundamental factors that control subsequent parasite evolution, such as the size of the switching population or the extent of immigration from the original host, remain largely unknown. To shed light on the host switching process, we explored how host switches occur in independent host shifts by two ectoparasitic honey bee mites (Varroa destructor and V. jacobsoni). Both switched to the western honey bee (Apis mellifera) after being brought into contact with their ancestral host (Apis cerana), ~70 and ~12 years ago, respectively. Varroa destructor subsequently caused worldwide collapses of honey bee populations. Using whole-genome sequencing on 63 mites collected in their native ranges from both the ancestral and novel hosts, we were able to reconstruct the known temporal dynamics of the switch. We further found multiple previously undiscovered mitochondrial lineages on the novel host, along with the genetic equivalent of tens of individuals that were involved in the initial host switch. Despite being greatly reduced, some gene flow remains between mites adapted to different hosts. Our findings suggest that while reproductive isolation may facilitate the fixation of traits beneficial for exploiting the new host, ongoing genetic exchange may allow genetic amelioration of inbreeding effects.

LncRNAs at the heart of development and disease

Long noncoding RNAs (LncRNAs) have emerged as a diverse class of functional molecules that contribute to nearly every facet of mammalian cardiac development and disease. Recent examples show that lncRNAs can be important co-regulators of cardiac patterning and morphogenesis and modulators of the pathogenic signaling that drives heart disease. The flexibility and chemical nature of RNA allows lncRNAs to utilize diverse mechanisms, mediating their effects through their sequence, structure, and molecular interactions with DNA, protein, and other RNAs. In vivo, i.e., animal, studies of individual lncRNAs highlight their ability to balance conserved cardiac gene expression networks, serve as specific and early biomarkers, and indicate their promise as useful therapeutic targets to treat human heart disease. Here, we review recent functionally characterized lncRNAs in cardiac biology and pathology and provide a perspective on emerging approaches to decipher the role of lncRNAs in the heart.

Geometric morphology and population genomics provide insights into the adaptive evolution of Apis cerana in Changbai Mountain

Background

Exploration of adaptive evolution of organisms in response to environmental change can help to understand the evolutionary history of species and the underlying mechanisms of adaptation to local environments, thus guiding future conservation programmes. Before the introduction of Apis mellifera in China, eastern honey bees (Apis cerana) were the only species used for beekeeping in this region. In the mountains of Changbai, populations of A. cerana are considered a distinct ecotype of the species which formed through the distinct selective pressures in this area over time.

Result

We performed a measure of 300 wing specimens of eastern honey bees and obtained the geometric morphological variation in the wing of A. cerana in Changbai Mountain. A total of 3,859,573 high-quality SNP loci were yielded via the whole-genome resequencing of 130 individuals in 5 geographic regions.

Conclusion

Corresponding geometric morphology and population genomics confirmed the particularity of the A. cerana in Changbai Mountain. Genetic differentiation at the subspecies level exists between populations in Changbai Mountain and remaining geographic regions, and a significant reduction in the effective population size and an excessive degree of inbreeding may be responsible for a substantial loss of population genetic diversity. Candidate genes potentially associated with cold environmental adaptations in populations under natural selection were identified, which may represent local adaptations in populations. Our study provided insights into the evolutionary history and adaptation of A. cerana in Changbai Mountain, as well as its conservation.

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.

De Novo Genome Assembly of Chinese Plateau Honeybee Unravels Intraspecies Genetic Diversity in the Eastern Honeybee, Apis cerana

Simple Summary

In this study, we obtained a chromosome-scale assembly genome of Apis cerana abansis, which lives in the southeastern margin of the Titan Plateau, by using PacBio, Illumina and high-throughput chromatin conformation capture (Hi-C) sequencing technologies. With a more comprehensive annotation pipeline, we obtained an ampler and more accurate Apis cerana genome than previous studies. Comparative genomic analysis was performed to identify the divergence among different A. cerana genomes by studying two aspects: the differential content of repeat content and the gene loss/gain events occurred in chemosensory receptors and immune-related proteins. Our results show that the content of repetitive sequences differ in types and quantity among four A. cerana strains; the gene loss/gain events in chemoreceptor- and immune-related proteins occur in different A. cerana strains, especially in A. cerana abansis (Aba strain). Specifically, while compared with the other three published genomes, the Aba strain contains the largest number of repeat contents and loses the largest number of both chemosensory-receptor- and immune-related proteins, as well as subfamilies, whereas the Baisha strain contains the largest number of chemoreceptor- and immune-related proteins. We hypothesized that gene loss/gain may be evolutionary strategies used by the different A. cerana strains to adapt to their respective environments.

Abstract

Apis cerana abansis, widely distributed in the southeastern margin of the Qinghai-Tibet Plateau, is considered an excellent model to study the phenotype and genetic variation for highland adaptation of Asian honeybee. Herein, we assembled and annotated the chromosome-scale assembly genome of A. cerana abansis with the help of PacBio, Illumina and Hi-C sequencing technologies in order to identify the genome differences between the A. cerana abansis and the published genomes of different A. cerana strains. The sequencing methods, assembly and annotation strategies of A. cerana abansis were more comprehensive than previously published A. cerana genomes. Then, the intraspecific genetic diversity of A. cerana was revealed at the genomic level. We re-identified the repeat content in the genome of A. cerana abansis, as well as the other three A. cerana strains. The chemosensory and immune-related proteins in different A. cerana strains were carefully re-identified, so that 132 odorant receptor subfamilies, 12 gustatory receptor subfamilies and 22 immune-related pathways were found. We also discovered that, compared with other published genomes, the A. ceranaabansis lost the largest number of chemoreceptors compared to other strains, and hypothesized that gene loss/gain might help different A. cerana strains to adapt to their respective environments. Our work contains more complete and precise assembly and annotation results for the A. cerana genome, thus providing a resource for subsequent in-depth related studies.

Identification of Species by Combining Molecular and Morphological Data Using Convolutional Neural Networks

Integrative taxonomy is central to modern taxonomy and systematic biology, including behaviour, niche preference, distribution, morphological analysis and DNA barcoding. However, decades of use demonstrate that these methods can face challenges when used in isolation, for instance, potential misidentifications due to phenotypic plasticity for morphological methods, and incorrect identifications because of introgression, incomplete lineage sorting and horizontal gene transfer for DNA barcoding. Although researchers have advocated the use of integrative taxonomy, few detailed algorithms have been proposed. Here, we develop a convolutional neural network method (morphology-molecule network (MMNet)) that integrates morphological and molecular data for species identification. The newly proposed method (MMNet) worked better than four currently-available alternative methods when tested with 10 independent datasets representing varying genetic diversity from different taxa. High accuracies were achieved for all groups, including beetles (98.1% of 123 species), butterflies (98.8% of 24 species), fishes (96.3% of 214 species) and moths (96.4% of 150 total species). Further, MMNet demonstrated a high degree of accuracy (>98%) in four datasets including closely related species from the same genus. The average accuracy of two modest sub-genomic (single nucleotide polymorphism) datasets, comprising eight putative subspecies respectively, is 90%. Additional tests show that the success rate of species identification under this method most strongly depends on the amount of training data, and is robust to sequence length and image size. Analyses on the contribution of different data types (image versus gene) indicate that both morphological and genetic data are important to the model, and that genetic data contribute slightly more. The approaches developed here serve as a foundation for the future integration of multi-modal information for integrative taxonomy, such as image, audio, video, 3D scanning and biosensor data, to characterize organisms more comprehensively as a basis for improved investigation, monitoring and conservation of biodiversity.

Parallel adaptation in autopolyploid Arabidopsis arenosa is dominated by repeated recruitment of shared alleles

Relative contributions of pre-existing vs de novo genomic variation to adaptation are poorly understood, especially in polyploid organisms. We assess this in high resolution using autotetraploid Arabidopsis arenosa, which repeatedly adapted to toxic serpentine soils that exhibit skewed elemental profiles. Leveraging a fivefold replicated serpentine invasion, we assess selection on SNPs and structural variants (TEs) in 78 resequenced individuals and discover significant parallelism in candidate genes involved in ion homeostasis. We further model parallel selection and infer repeated sweeps on a shared pool of variants in nearly all these loci, supporting theoretical expectations. A single striking exception is represented by TWO PORE CHANNEL 1, which exhibits convergent evolution from independent de novo mutations at an identical, otherwise conserved site at the calcium channel selectivity gate. Taken together, this suggests that polyploid populations can rapidly adapt to environmental extremes, calling on both pre-existing variation and novel polymorphisms.

Relative contributions of pre-existing versus de novo genomic variation to adaptation remain unclear. Here, the authors address this problem by examining the adaptation of autotetraploid Arabidopsis arenosa to serpentine soils and find that both types of variations contribute to rapid adaptation.

Tachykinin signaling inhibits task-specific behavioral responsiveness in honeybee workers

Behavioral specialization is key to the success of social insects and leads to division of labor among colony members. Response thresholds to task-specific stimuli are thought to proximally regulate behavioral specialization, but their neurobiological regulation is complex and not well understood. Here, we show that response thresholds to task-relevant stimuli correspond to the specialization of three behavioral phenotypes of honeybee workers in the well-studied and important Apis mellifera and Apis cerana. Quantitative neuropeptidome comparisons suggest two tachykinin-related peptides (TRP2 and TRP3) as candidates for the modification of these response thresholds. Based on our characterization of their receptor binding and downstream signaling, we confirm a functional role of tachykinin signaling in regulating specific responsiveness of honeybee workers: TRP2 injection and RNAi-mediated downregulation cause consistent, opposite effects on responsiveness to task-specific stimuli of each behaviorally specialized phenotype but not to stimuli that are unrelated to their tasks. Thus, our study demonstrates that TRP signaling regulates the degree of task-specific responsiveness of specialized honeybee workers and may control the context specificity of behavior in animals more generally.

Leucokinin and Associated Neuropeptides Regulate Multiple Aspects of Physiology and Behavior in Drosophila

Leucokinins (LKs) constitute a family of neuropeptides identified in numerous insects and many other invertebrates. LKs act on G-protein-coupled receptors that display only distant relations to other known receptors. In adult Drosophila, 26 neurons/neurosecretory cells of three main types express LK. The four brain interneurons are of two types, and these are implicated in several important functions in the fly's behavior and physiology, including feeding, sleep-metabolism interactions, state-dependent memory formation, as well as modulation of gustatory sensitivity and nociception. The 22 neurosecretory cells (abdominal LK neurons, ABLKs) of the abdominal neuromeres co-express LK and a diuretic hormone (DH44), and together, these regulate water and ion homeostasis and associated stress as well as food intake. In Drosophila larvae, LK neurons modulate locomotion, escape responses and aspects of ecdysis behavior. A set of lateral neurosecretory cells, ALKs (anterior LK neurons), in the brain express LK in larvae, but inconsistently so in adults. These ALKs co-express three other neuropeptides and regulate water and ion homeostasis, feeding, and drinking, but the specific role of LK is not yet known. This review summarizes Drosophila data on embryonic lineages of LK neurons, functional roles of individual LK neuron types, interactions with other peptidergic systems, and orchestrating functions of LK.

Leucokinins: Multifunctional Neuropeptides and Hormones in Insects and Other Invertebrates

Leucokinins (LKs) constitute a neuropeptide family first discovered in a cockroach and later identified in numerous insects and several other invertebrates. The LK receptors are only distantly related to other known receptors. Among insects, there are many examples of species where genes encoding LKs and their receptors are absent. Furthermore, genomics has revealed that LK signaling is lacking in several of the invertebrate phyla and in vertebrates. In insects, the number and complexity of LK-expressing neurons vary, from the simple pattern in the Drosophila larva where the entire CNS has 20 neurons of 3 main types, to cockroaches with about 250 neurons of many different types. Common to all studied insects is the presence or 1-3 pairs of LK-expressing neurosecretory cells in each abdominal neuromere of the ventral nerve cord, that, at least in some insects, regulate secretion in Malpighian tubules. This review summarizes the diverse functional roles of LK signaling in insects, as well as other arthropods and mollusks. These functions include regulation of ion and water homeostasis, feeding, sleep-metabolism interactions, state-dependent memory formation, as well as modulation of gustatory sensitivity and nociception. Other functions are implied by the neuronal distribution of LK, but remain to be investigated.