The giant honeybee, Apis dorsata, is an important pollinator native to tropical and subtropical Asia. Understanding its potential response to climate change is important for pollinator conservation, ecological risk assessment, and biosecurity planning. This study used an optimized MaxEnt ecological niche modeling framework to predict the current and future global habitat suitability of A. dorsata. The model was developed using 1060 occurrence records and seven non-collinear bioclimatic variables and projected under three global climate models (IPSL-CM6A-LR, BCC-CSM2-MR, and MPI-ESM1-2-HR) and three Shared Socioeconomic Pathways (SSP126, SSP245, and SSP585) for 2041–2060, centered on 2050. The optimized model used linear, quadratic, and hinge features (LQH) with a regularization multiplier of 0.5 and demonstrated good predictive performance under 5-fold spatial cross-validation (mean AUC = 0.905 ± 0.007; TSS = 0.746 ± 0.010). The results indicate that the potential distribution of A. dorsata is primarily associated with the combined effects of seasonal temperature and moisture availability. Current projections identified high climatic suitability across South and Southeast Asia, while also revealing potentially suitable environments in parts of Africa, the Americas, and northern Australia. Future projections suggest that suitable climatic conditions will largely persist through 2050, with habitat gains generally exceeding losses and increasing under higher climate-forcing scenarios. Continued monitoring and proactive biosecurity are essential to address habitat loss within the native range and prevent establishment in newly suitable regions. This study highlights the potential redistribution of A. dorsata under future climate change.

Cuticular proteins (CPs)—key components of the insect exoskeleton—not only regulate development but also serve as structural barriers that enhance resistance against environmental stressors. This study identified CP gene families in Apis mellifera and analyzed their expression patterns during the worker capped brood development stages from mature larva to pre-eclosion. Using a comprehensive genome-wide bioinformatic approach, we identified 85 CP genes in A. mellifera which comprise six families: CPR (n = 43), CPAPs (n = 27), CPF (n = 2), Tweedle (n = 2), CPLCP (n = 8) and Apidermin (n = 3). Analysis of CP gene evolutionary relationship revealed that each CP family forms a distinct, relatively independent clade. Domain and motif analyses confirmed that all CPR members harbor a conserved Chitin_Bind_4 domain, consistent with CPR family structures in other taxa. Additionally, CPAP members possess one or three Chitin-binding Peritrophin-A domain (CBM_14), CPF members possess a conserved Pupal cuticle protein C1 domain (Cuticle_3), and Tweedle members contain a conserved domain of unknown function (DUF243). In addition, the analysis found no conserved domain within the CPLCP and Apidermin families. RNA-seq data revealed dynamic expression patterns of AmCPs during pupal development, with each gene family displaying a relatively characteristic temporal profile. Quantitative PCR validation of eight highly expressed CPR genes at 9 days post-capping confirmed the RNA-seq results. This work provides a comprehensive bioinformatic characterization and transcriptional analysis of CP genes in A. mellifera, offering a foundation for future functional studies on cuticle formation and identifying candidate genes potentially involved in cuticle development in honeybees. This work relies on transcriptomic data and in silico analyses. All proposed biological roles are hypothetical and require experimental validation.