Van Yüzüncü Yıl University, Faculty of Agriculture, Department of Agricultural Biotechnology, Van
Abstract
Wing shape provides an important source of morphological information for evaluating population differentiation in honey bees, but the extent of geographic variation relative to allometry and within-population variation remains unclear. This study examined forewing shape and size variation in 1.040 honey bee (Apis mellifera) specimens collected from seven locations in Türkiye. Nineteen homologous landmarks were digitized and analysed using Generalized Procrustes Analysis, Procrustes ANOVA, allometric models, and clustering, and reduced landmark-displacement features were used to assess individual-level geographic assignment. Overall forewing size, measured as centroid size, did not differ significantly among locations (ANOVA: F₆,₁₀₃₃=0.89, P=0.502; Kruskal–Wallis: P=0.483). In contrast, forewing shape differed significantly among locations (Procrustes ANOVA: F₆,₁₀₃₃=3.06, Z=6.84, P=0.001), although the geographic effect was small, explaining only 1.75% of total shape variation. Allometry explained a larger proportion of shape variation (R²=4.83%, P=0.001), while the size × location interaction was not significant (P=0.072). Ward hierarchical clustering suggested two broad shape groups, but these were not meaningfully associated with geographic origin (Cramér's V=0.058, P=0.739). Individual-level geographic assignment based on reduced landmark-displacement features was limited, with the highest test accuracy reaching 27.6%, only marginally above the majority-class baseline of approximately 26.6%. Together, these results indicate that the sampled honey bee populations from seven locations in Türkiye show statistically detectable but biologically modest location-associated variation in forewing shape, which is outweighed by allometric and within-population variation and is insufficient for reliable individual-level geographic assignment.
TOZKAR, C. Özge, & ÇAKMAKÇI, C. (2026). Allometry and Geographic Variation in Forewing Shape and Size Among Honey Bee (Apis mellifera L.) Populations from Türkiye. ISPEC Journal of Agricultural Sciences, 10(3), 779–791. https://doi.org/10.5281/zenodo.21619414
📄Adams, D.C., Collyer, M.L., Kaliontzopoulou, A., Baken, E.K., 2021. Geomorph: Software for geometric morphometric analyses. R package version 4.0. https://cran.rproject.org/package=geomorph (Accessed: 01.01.2021).
📄Baken, E.K., Collyer, M.L., Kaliontzopoulou, A., Adams, D.C., 2021. Geomorph v4.0 and gmShiny: enhanced analytics and a new graphical interface for a comprehensive morphometric experience. Methods in Ecology and Evolution, 12(12): 2355-2363.
📄Bodur, Ç., Kence, M., Kence, A., 2007. Genetic structure of honeybee, Apis mellifera L. (Hymenoptera: Apidae) populations of Turkey inferred from microsatellite analysis. Journal of Apicultural Research, 46(1): 50-56.
📄Breiman, L., 2001. Random forests. Machine Learning, 45(1): 5-32.
📄Collyer, M.L., Adams, D.C., Freckleton, R., 2018. RRPP: an R package for fitting linear models to high-dimensional data using residual randomization. Methods in Ecology and Evolution, 9(7): 1772-1779.
📄Çavuş, E.D., Tozkar, C.Ö., 2026. Geographic classification of social bees using automated wing pattern analysis. Entomologia Experimentalis et Applicata: 1–18.
📄Dadgostar, S., Mobini Dehkordi, P., Khodarahmi, S., Ghorbani, M., 2025. Automated honey bee subspecies identification using advanced wing venation analysis and adaptive hierarchical clustering. Ecology and Evolution, 15(9): e72101.
📄De Nart, D., Costa, C., Di Prisco, G., Carpana, E., 2022. Image recognition using convolutional neural networks for classification of honey bee subspecies. Apidologie, 53(1): 5.
📄Dryden, I.L., Mardia, K.V., 2016. Statistical Shape Analysis, with Applications in R (2nd ed.). John Wiley & Sons, Chichester.
📄Francoy, T.M., Wittmann, D., Drauschke, M., Müller, S., Steinhage, V., Bezerra-Laure, M.A.F., De Jong, D., Gonçalves, L.S., 2008. Identification of Africanized honey bees through wing morphometrics: two fast and efficient procedures. Apidologie, 39(5): 488-494.
📄Greenwell, B.M., Boehmke, B.C., 2020. Variable importance plots—an introduction to the vip package. The R Journal, 12(1): 343-366.
📄Güler, A., 2010. A morphometric model for determining the effect of commercial queen bee usage on the native honeybee (Apis mellifera L.) population in a Turkish province. Apidologie, 41(6): 622-635.
📄Hailu, T.G., D'Alvise, P., Hasselmann, M., 2021. Disentangling Ethiopian honey bee (Apis mellifera) populations based on standard morphometric and genetic analyses. Insects, 12(3): 193.
📄Henriques, D., Chávez-Galarza, J., Teixeira, J.S.G., Ferreira, H., Neves, C.J., Francoy, T.M., Pinto, M.A., 2020. Wing geometric morphometrics of workers and drones and single nucleotide polymorphisms provide similar genetic structure in the Iberian honey bee (Apis mellifera iberiensis). Insects, 11(2): 89.
📄Kandemir, İ., Kence, M., Kence, A., 2000. Genetic and morphometric variation in honeybee (Apis mellifera L.) populations of Turkey. Apidologie, 31(3): 343-356.
📄Kekeçoğlu, M., Kambur, M., Uçak, M., Çaprazlı, T., Bir, S., 2020. Biodiversity of honey bees (Apis mellifera L.) in Turkey by geometric morphometric analysis. Biological Diversity and Conservation, 13(3): 282-289.
📄Klingenberg, C.P., 2016. Size, shape, and form: concepts of allometry in geometric morphometrics. Development Genes and Evolution, 226(3): 113-137.
📄Kuhn, M., 2008. Building predictive models in R using the caret package. Journal of Statistical Software, 28(5): 1-26.
📄Kursa, M.B., Rudnicki, W.R., 2010. Feature selection with the Boruta package. Journal of Statistical Software, 36(11): 1-13.
📄Masaquiza, D., Ferrán, M.O., Guamán, S., Naranjo, E., Vaca, M., Curbelo, L.M., Arenal, A., 2023. Geometric morphometric analysis of wing shape to identify populations of Apis mellifera in Camagüey, Cuba. Insects, 14(3): 306.
📄Masaquiza, D., Curbelo Rodríguez, L., Zapata, J., Monar, J., Vaca, M., Porrini, L., Eguaras, M., Daniele, M., Romero, D., Arenal, A., 2024. Use of wing geometric morphometric analysis and mtDNA to identify Africanization of Apis mellifera in the central highlands of Ecuador. Insects, 15(8): 628.
📄Nawrocka, A., Kandemir, İ., Fuchs, S., Tofilski, A., 2018. Computer software for identification of honey bee subspecies and evolutionary lineages. Apidologie, 49(2): 172-184.
📄Oleksa, A., Tofilski, A., 2015. Wing geometric morphometrics and microsatellite analysis provide similar discrimination of honey bee subspecies. Apidologie, 46(1): 49-60.
📄Özdil, F., Aytekin, İ., İlhan, F., Boztepe, S., 2012. Genetic variation in Turkish honeybees Apis mellifera anatoliaca, A. m. caucasica, A. m. meda (Hymenoptera: Apidae) inferred from RFLP analysis of three mtDNA regions (16S rDNA–COI–ND5). European Journal of Entomology, 109(2): 161-167.
📄R Core Team, 2025. R:A Language and Environment for Statistical Computing (Version 4.5.0). R Foundation for Statistical Computing, Vienna. https://www.R-project.org/ (Accessed: 02.02.2025).
📄Rohlf, F.J., 2015. The tps series of software. Hystrix, the Italian Journal of Mammalogy, 26(1): 9-12.
📄Ruttner, F., 1988. Biogeography and Taxonomy of Honeybees. Springer-Verlag, Berlin.
📄Tofilski, A., 2008. Using geometric morphometrics and standard morphometry to discriminate three honeybee subspecies. Apidologie, 39(5): 558-563.
📄Zelditch, M.L., Swiderski, D.L., Sheets, H.D., 2004. Geometric Morphometrics for Biologists: A Primer. Elsevier, San Diego.
📄Zhang, G.P., 2000. Neural networks for classification: a survey. IEEE Transactions on Systems, Man, and Cybernetics, Part C (Applications and Reviews), 30(4): 451-462.