There is growing concern that massive loss of honey bees can cause serious negative effects on biodiversity and ecosystems. Surveys of colony losses have been performed worldwide to monitor the dynamic changes and health status of honey bee colonies. Here, we present the results of surveys regarding winter colony losses from 21 provinces in China from 2009 to 2021, with a total of 1,744,324 colonies managed by 13,704 beekeepers. The total colony losses were low (9.84%; 95% Confidence Interval (CI): 9.60–10.08%) but varied among years, provinces, and scales of apiaries. As little is known about the overwintering mortality of Apis cerana, in this study, we surveyed and compared the loss rates between Apis mellifera and A. cerana in China. We found colonies of A. mellifera suffered significantly lower losses than A. cerana in China. Larger apiaries resulted in higher losses in A. mellifera, whereas the opposite was observed in A. cerana. Furthermore, we used generalized linear mixed-effects models (GLMMs) to evaluate the effects of potential risk factors on winter colony losses and found that the operation size, species, migration, migration×species interaction, and queen problems were significantly related to the loss rates. New queens can increase their colony overwintering survival. Migratory beekeepers and large operations reported lower loss rates.

The geographic and biological diversity of China has resulted in the differential adaptation of the eastern honeybee, Apis cerana, to these varied habitats. A. cerana were collected from 14 locations in China. Their genomes were sequenced, and nucleotide polymorphisms were identified at more than 9 million sites. Both STRUCTURE and principal component analysis placed the bees into seven groups. Phylogenomic analysis groups the honeybees into many of the same clusters with high bootstrap values (91%–100%). Populations from Tibet and South Yunnan are sister taxa and together represent the earliest diverging lineage included in this study. We propose that the evolutionary origin of A. cerana in China was in the southern region of Yunnan Province and expanded from there into the southeastern regions and into the northeastern mountain regions. The Cold‐Temperate West Sichuan Plateau and Tropical Diannan populations were compared to identify genes under adaptive selection in these two habitats. Pathway enrichment analysis showing genes under selection, including the Hippo signaling pathway, GABAergic pathway, and trehalose‐phosphate synthase, indicates that most genes under selection pressure are involved in the process of signal transduction and energy metabolism. qRT‐PCR analysis reveals that one gene under selection, the AcVIAAT gene, involved in the GABAergic pathway, is responding to cold temperature stress. Through homologous recombination, we show that the AcVIAAT gene is able to replace the CNAG_01904 gene in the fungus Cryptococcus neoformans and that it makes the fungus less sensitive to conditions of oxidative stress and variations in temperature. Our results contribute to our understanding of the evolutionary origin of A. cerana in China and the molecular basis of environmental adaptation.

The microsporidian parasite designated here as Nosema sp. Isolate YNPr was isolated from the cabbage butterfly Pieris rapae collected in Honghe Prefecture, Yunnan Province, China. The genome was sequenced by Illumina sequencing and compared to those of two related members of the Nosema/Vairimorpha clade, Nosema ceranae and Nosema apis. Based upon assembly statistics, the Nosema sp. YNPr genome is 3.36 x 10⁶bp with a G+C content of 23.18% and 2,075 protein coding sequences. An “ACCCTT” motif is present approximately 50-bp upstream of the start codon, as reported from other members of the clade and from Encephalitozoon cuniculi, a sister taxon. Comparative small subunit ribosomal DNA (SSU rDNA) analysis as well as genome-wide phylogenetic analysis confirms a closer relationship between N. ceranae and Nosema sp. YNPr than between the two honeybee parasites N. ceranae and N. apis. The more closely related N. ceranae and Nosema sp. YNPr show similarities in a number of structural characteristics such as gene synteny, gene length, gene number, transposon composition and gene reduction. Based on transposable element content of the assemblies, the transposon content of Nosema sp. YNPr is 4.8%, that of N. ceranae is 3.7%, and that of N. apis is 2.5%, with large differences in the types of transposons present among these 3 species. Gene function annotation indicates that the number of genes participating in most metabolic activities is similar in all three species. However, the number of genes in the transcription, general function, and cysteine protease categories is greater in N. apis than in the other two species. Our studies further characterize the evolution of the Nosema/Vairimorpha clade of microsporidia. These organisms maintain variable but very reduced genomes. We are interested in understanding the effects of genetic drift versus natural selection on genome size in the microsporidia and in developing a testable hypothesis for further studies on the genomic ecology of this group.

The Himalayan giant honeybee, Apis laboriosa, is the largest individual honeybee with major ecological and economic importance in high-latitude environments. However, our understanding of its environmental adaptations is circumscribed by the paucity of genomic data for this species. Here, we provide a draft genome of wild A. laboriosa, along with a comparison to its closely related species, Apis dorsata. The draft genome of A. laboriosa based on the de novo assembly is 226.1 Mbp in length with a scaffold N50 size of 3.34 Mbp, a GC content of 32.2%, a repeat content of 6.86%, and a gene family number of 8,404. Comparative genomics analysis revealed that the genes in A. laboriosa genome have undergone stronger positive selection (2.5 times more genes) and more recent duplication/loss events (6.1 times more events) than those in the A. dorsata genome. Our study implies the potential molecular mechanisms underlying the high-altitude adaptation of A. laboriosa and will catalyze future comparative studies to understand the environmental adaptation of modern honeybees.

Apis laboriosa is the largest honeybee that lives mainly on cliff faces, with strong migratory ability. In this study, we firstly sequenced and assembled two complete mitochondrial genomes of A. laboriosa isolated from two distant locations in China (Chongqing and Shangri‐La regions). Combined with the published mitochondrial genome of A. laboriosa from Nepal, comparative genomic analyses were conducted to gain insight into the genetic diversity of giant honeybees from different geographical distributions. The mitochondrial genomes of A. laboriosa from Chongqing and Shangri‐La regions were 15,579 and 15,683 bp in length, respectively, both larger than that from Nepal with the length of 15,510 bp. Three mitochondrial genomes all harbor 37 common genes and present the same AT bias and the frequency of codon usage. However, the fragments including COX1, SSUrRNA, LSUrRNA, and the AT‐rich region of the mitochondrial genome from Shangri‐La region demonstrate distinctive insertions and deletions compared to those from Chongqing and Nepal regions. Phylogenetic trees of mitochondrial genomes show that A. laboriosa from Chongqing is most closely related to that from Nepal, rather than to Shangri‐La. Genetic distance between Shangri‐La and Chongqing or Nepal was even larger than that between the various subspecies of Apis mellifera. Overall, these results unmark that A. laboriosa in different geographical distributions can exhibit high genetic diversity at the mitochondrial genomic level, and therein, A. laboriosa from Shangri‐La may be the subspecies. All these studies will contribute to our understanding of the geographical distribution and genetic differentiation of black giant honeybee in Asian region.