Abstract
Honey carries environmental DNA (eDNA) from organisms encountered by honey bees and can therefore encode an “entomological fingerprint” of agricultural and forest landscapes. We compared two metabarcoding assays targeting mitochondrial regions (cytochrome c oxidase subunit 1, COI; and cytochrome b, CYTB) to profile plant-sucking insects (Hemiptera) in honey samples from Italy and Türkiye, including citrus blossom, polyfloral, and honeydew honey samples. We designed a new short CYTB primer pair optimised for Aphididae and evaluated it alongside our previously published Hemiptera-COI assay. Across 3.51 million reads, our pipeline assigned 1.25 million COI and 1.48 million CYTB reads to Neoptera. CYTB resolved more taxa per sample (mean 7.9 families, 36 species) than COI (4.5 families, 11.7 species) and yielded a higher proportion of species-level assignments (98.3% vs. 43.7%). COI captured Metcalfa pruinosa (Flatidae) where abundant—often dominating polyfloral honey samples—whereas CYTB revealed fine-scale diversity of aphids relevant to citrus orchards and conifer/oak forests. Short-amplicon haplotype screening recovered multiple COI haplotypes for M. pruinosa, Thelaxes suberi, Cinara cedri and Aphis gossypii, suggesting potential for population monitoring from honey eDNA. Our results show that combining complementary metabarcodes mitigates primer and database biases and enhances landscape-scale inference, honey authentication, and surveillance of invasive or pest hemipterans.
Data availability
The datasets generated during the current study showing alignments of CYTB and COI primers with Hemiptera sequences are available in Zenodo with the following DOI: https://doi.org/10.5281/zenodo.19183389 . The sequence datasets generated and analysed during the current study are available in the EMBL-EBI European Nucleotide Archive (ENA) repository http://www.ebi.ac.uk/ena , with the project accession number PRJEB102161.
References
Codex Alimentarius – International Food Standards. Standard for honey. Adopted in 1981. Revised in 1987, 2001. Amended in 2019, 2022. Food and Agriculture Organization of the United Nations-World Health Organization. CXS 12–1981, 1-8 (2022).
Pita-Calvo, C. & Vázquez, M. Differences between honeydew and blossom honeys: a review. Trends Food Sci. Technol. 57, 79–87 (2017).
da Silva, P. M., Gauche, C., Gonzaga, L. V., Costa, A. C. O. & Fett, R. Honey: Chemical composition, stability and authenticity. Food Chem. 196, 309–323 (2016).
De-Melo, M., Almeida-Muradian, A. A., Sancho, L. B. D., Pascual-Maté, A. & M. T., & Composition and properties of Apis mellifera honey: A review. J. Apic. Res. 57, 5–37 (2018).
Laube, I. et al. Development of primer and probe sets for the detection of plant species in honey. Food Chem. 118, 979–986 (2010).
Bruni, I. et al. A DNA barcoding approach to identify plant species in multiflower honey. Food Chem. 170, 308–315 (2015).
Hawkins, J. et al. Using DNA metabarcoding to identify the floral composition of honey: a new tool for investigating honey bee foraging preferences. PLoS One 10, e0134735 (2015).
Torricelli, M., Pierboni, E., Tovo, G. R., Curcio, L. & Rondini, C. In-house validation of a DNA extraction protocol from honey and bee pollen and analysis in fast Real-Time PCR of commercial honey samples using a knowledge-based approach. Food Anal. Methods. 9, 3439–3450 (2016).
Kast, C. & Roetschi, A. Evaluation of baker’s yeast in honey using a real-time PCR assay. Food Microbiol. 62, 282–288 (2017).
Bovo, S. et al. Shotgun metagenomics of honey DNA: Evaluation of a methodological approach to describe a multi-kingdom honey bee derived environmental DNA signature. PLoS One 13, e0205575 (2018).
Utzeri, V. J. et al. Application of next generation semiconductor based sequencing to detect the botanical composition of monofloral, polyfloral and honeydew honey. Food Control 86, 342–349 (2018).
Utzeri, V. J., Ribani, A. & Fontanesi, L. Authentication of honey based on a DNA method to differentiate Apis mellifera subspecies: Application to Sicilian honey bee (A. m. siciliana) and Iberian honey bee (A. m. iberiensis) honeys. Food Control 91, 294–301 (2018).
Ribani, A., Utzeri, V. J., Taurisano, V. & Fontanesi, L. Honey as a source of environmental DNA for the detection and monitoring of honey bee pathogens and parasites. Vet. Sci. 7, 113 (2020).
Utzeri, V. J. et al. A next generation sequencing approach for targeted Varroa destructor (Acari: Varroidae) mitochondrial DNA analysis based on honey derived environmental DNA. J. Invertebr. Pathol. 161, 47–53 (2019).
Bovo, S., Utzeri, V. J., Ribani, A., Cabbri, R. & Fontanesi, L. Shotgun sequencing of honey DNA can describe honey bee derived environmental signatures and the honey bee hologenome complexity. Sci. Rep. 10, 9279 (2020).
Ribani, A., Utzeri, V. J., Taurisano, V., Galuppi, R. & Fontanesi, L. Analysis of honey environmental DNA indicates that the honey bee (Apis mellifera L.) trypanosome parasite Lotmaria passim is widespread in the apiaries of the North of Italy. J. Invertebr. Pathol. 184, 107628 (2021).
Ribani, A., Taurisano, V., Utzeri, V. J. & Fontanesi, L. Honey environmental DNA can be used to detect and monitor honey bee pests: Development of methods useful to identify Aethina tumida and Galleria mellonella infestations. Vet. Sci. 9, 213 (2022).
Utzeri, V. J. et al. Entomological signatures in honey: An environmental DNA metabarcoding approach can disclose information on plant-sucking insects in agricultural and forest landscapes. Sci. Rep. 8, 9996 (2018).
Utzeri, V. J., Ribani, A., Taurisano, V., Banqué, C. H. I. & Fontanesi, L. Distribution of the main Apis mellifera mitochondrial DNA lineages in Italy assessed using an environmental DNA approach. Insects 12, 620 (2021).
Bovo, S. et al. A genotyping by sequencing approach can disclose Apis mellifera population genomic information contained in honey environmental DNA. Sci. Rep. 12, 19541 (2022).
Utzeri, V. J., Ribani, A., Taurisano, V. & Fontanesi, L. Entomological authentication of honey based on a DNA method that distinguishes Apis mellifera mitochondrial C mitotypes: Application to honey produced by A. m. ligustica and A. m. carnica. Food Control 134, 108713 (2022).
Kek, S. P., Chin, N. L., Tan, S. W., Yusof, Y. A. & Chua, L. S. Molecular identification of honey entomological origin based on bee mitochondrial 16S rRNA and COI gene sequences. Food Control 78, 150–159 (2017).
Prosser, S. W. & Hebert, P. D. Rapid identification of the botanical and entomological sources of honey using DNA metabarcoding. Food Chem. 214, 183–191 (2017).
Khansaritoreh, E. et al. Employing DNA metabarcoding to determine the geographical origin of honey. Heliyon 6, e05596 (2020).
Milla, L. et al. Pollen DNA metabarcoding identifies regional provenance and high plant diversity in Australian honey. Ecol. Evol. 11, 8683–8698 (2021).
Wirta, H., Abrego, N., Miller, K., Roslin, T. & Vesterinen, E. DNA traces the origin of honey by identifying plants, bacteria and fungi. Sci. Rep. 11, 4798 (2021).
Laha, R. C. et al. Meta-barcoding in combination with palynological inference is a potent diagnostic marker for honey floral composition. AMB Express 7, 132 (2017).
Espinosa Prieto, A., Hardion, L., Debortoli, N. & Beisel, J. N. Finding the perfect pairs: A matchmaking of plant markers and primers for multi-marker eDNA metabarcoding. Mol. Ecol. Resour. 24, e13937 (2024).
Keck, F., Couton, M. & Altermatt, F. Navigating the seven challenges of taxonomic reference databases in metabarcoding analyses. Mol. Ecol. Resour. 23, 742–755 (2023).
Fontanesi, L., Taurisano, V., Ribani, A. & Utzeri, V. J. A reply to the Letter to the Editor of Moškrič et al. entitled A comment on the paper from Utzeri et al. (2022) Entomological authentication of honey based on a DNA method that distinguishes Apis mellifera mitochondrial C mitotypes: Application to honey produced by A. m. ligustica and A. m. carnica, Food Control 121, 2021, 107626. Food Control 147, 109570 (2023).
Taurisano, V. et al. Distribution of honey bee mitochondrial DNA haplotypes in an Italian region where a legislative act is protecting the Apis mellifera ligustica subspecies. Sci. Rep. 14, 20583 (2024).
Taurisano, V. et al. Comparing morphometric and mitochondrial DNA data from honeybees and honey samples for identifying Apis mellifera ligustica subspecies at the colony level. Animals 15, 1743 (2025).
NCBI Resource Coordinators Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 46, D8–D13 (2018).
Camacho, C. et al. BLAST+: architecture and applications. BMC Bioinform. 10, 421 (2009).
Fu, L., Niu, B., Zhu, Z., Wu, S. & Li, W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics 28, 3150–3152 (2012).
Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 17, 10–12 (2011).
Schmieder, R. & Edwards, R. Quality control and preprocessing of metagenomic datasets. Bioinformatics 27, 863–864 (2011).
Tamura, K., Stecher, G. & Kumar, S. MEGA 11: Molecular Evolutionary Genetics Analysis Version 11. Mol. Biol. Evol. 38, 3022–3027 (2021).
Tamura, K., Nei, M. & Kumar, S. Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc. Natl. Acad. Sci. USA 101, 11030–11035 (2004).
Park, C. G. et al. Genetic variability of the invasive species Metcalfa pruinosa (Hemiptera: Flatidae) in the Republic of Korea. J. Econ. Entomol. 109, 1897–1906 (2016).
Lee, K. H. et al. Tracing the invasion and expansion characteristics of the flatid planthopper, Metcalfa pruinosa (Hemiptera: Flatidae), in Korea using mitochondrial DNA sequences. Insects 12, 4 (2020).
Coeur d’acier, A. et al. DNA barcoding and the associated PhylAphidB@ se website for the identification of European aphids (Insecta: Hemiptera: Aphididae). PLoS One 9, e97620 (2014).
Liu, M., Clarke, L. J., Baker, S. C., Jordan, G. J. & Burridge, C. P. A practical guide to DNA metabarcoding for entomological ecologists. Ecol. Entomol. 45, 373–385 (2020).
Yang, J. et al. Zooplankton community profiling in a eutrophic freshwater ecosystem-lake tai basin by DNA metabarcoding. Sci. Rep. 7, 1773 (2017).
Keller, A. et al. Evaluating multiplexed next-generation sequencing as a method in palynology for mixed pollen samples. Plant. Biol. 17, 558–566 (2015).
Zangheri, S. & Donadini, P. Comparsa nel Veneto di un Omottero neartico: Metcalfa pruinosa Say (Homoptera, Flatidae). Redia 63, 301–304 (1980).
Lucchi, A. Metcalfa pruinosa and honey production in Italy. Am. Bee J. 137, 532–535 (1997).
Bagnoli, B. & Lucchi, A. Dannosità e Misure di Controllo Integrato. In: A Lucchi (ed.) La Metcalfa Negli Ecosistemi Italiani, 65–88. ARSIA, Firenze (2000).
Della Giustina, W. & Navarro, E. Metcalfa pruinosa, un nouvel envahisseur ? Phytoma 451, 30–32 (1993).
Dean, H. A. & Bailey, J. C. A flatid planthopper, Metcalfa pruinosa. J. Econ. Entomol. 54, 1104–1106 (1961).
Li, R., Li, M., Yan, J., Bai, M. & Zhang, H. Five mitochondrial genomes of the genus Eysarcoris Hahn, 1834 with phylogenetic implications for the Pentatominae (Hemiptera: Pentatomidae). Insects 12, 597 (2021).
Bai, H., Zhang, J., Dietrich, C. H., Li, Y. & Yuan, X. Structural characteristics of mitochondrial genomes of eight treehoppers (Hemiptera: Membracidae: Centrotinae) and their phylogenetic implications. Genes 14, 1510 (2023).
Zhou, X., Lei, Y., Dietrich, C. H. & Huang, M. Investigating monophyly of Typhlocybini based on complete mitochondrial genomes with characterization and comparative analysis of 19 species (Hemiptera: Cicadellidae: Typhlocybinae). Insects 14, 842 (2023).
Lin, S., Zhang, H., Hou, Y., Zhuang, Y. & Miranda, L. High-level diversity of dinoflagellates in the natural environment, revealed by assessment of mitochondrial cox1 and cob genes for dinoflagellate DNA barcoding. Appl. Environ. Microbiol. 75, 1279–1290 (2009).
Matsumoto, Y., Yanase, T., Tsuda, T. & Noda, H. Species-specific mitochondrial gene rearrangements in biting midges and vector species identification. Med. Vet. Entomol. 23, 47–55 (2009).
Stuart, O. P. et al. Morphological and molecular analysis of Australian earwigs (Dermaptera) points to unique species and regional endemism in the Anisolabididae family. Insects 10, 72 (2019).
Görür, G. & Demirtaş, E. Phenotypic plasticity in morphological features of Thelaxes suberi (Del Guercio, 1911) (Hemiptera: Aphidoidea) feed on Quercus (L.) species. J. Insect Biodivers. 5, 1–9 (2017).
Lieutier, F. & Ghaioule, D. Entomological research in Mediterranean forest ecosystems (Editions Quae., 2005).
Akyıldırım Beğen, H., Görür, G. & Şenol, Ö. Cinara (Hemiptera: Aphidoidea) species distributed in Turkey and their host plants. Turk. J. Biod. 2, 24–33 (2019).
Oğuzoğlu, Ş., Avcı, M. & Şenol, Ö. New contributions to the Turkish aphid fauna and species composition (Hemiptera: Aphididomorpha) in Isparta forests. Turk. J. Entomol. 46, 453–467 (2023).
Obrycki, J. J. & Kring, T. J. Predaceous Coccinellidae in biological control. Annu. Rev. Entomol. 43, 295–321 (1998).
Wang, X. Y., Yang, X. M., Lu, B., Zhou, L. H. & Wu, K. M. Genetic variation and phylogeographic structure of the cotton aphid, Aphis gossypii, based on mitochondrial DNA and microsatellite markers. Sci. Rep. 7, 1920 (2017).
Komazaki, S., Shigehara, T. & Toda, S. Diversity of Japanese Aphis gossypii and comparison with other Aphis species based on the mitochondrial cytochrome oxidase I sequence. Ann. Entomol. Soc. Am. 103, 916–924 (2010).
Rebijith, K. B., Asokan, R., Krishna, V., Kumar, N. K. & Ramamurthy, V. V. Development of species-specific markers and molecular differences in mitochondrial and nuclear DNA sequences of Aphis gossypii and Myzus persicae (Hemiptera: Aphididae). Fla. Entomol. 95, 674–682 (2012).
Lagos-Kutz, D., Favret, C., Giordano, R. & Voegtlin, D. J Molecular and morphological differentiation between Aphis gossypii Glover (Hemiptera, Aphididae) and related species, with particular reference to the North American Midwest. ZooKeys 459, 49 (2014).
Dong, Z., Li, Y., Liu, S. & Zhang, Z. Genetic diversity of the melon aphid Aphis gossypii (Hemiptera: Aphididae) in a diversified vegetable growing area. Appl. Entomol. Zool. 53, 259–266 (2018).
Acknowledgements
We thank beekeepers who provided honey samples, in particular Nicola Mele and Zeid Nabulsi.
Funding
This research was funded by the University of Bologna 2023–2025 RFO programs, by Regione Emilia-Romagna, BEE-RER-3 projects—CUP E37G22000030007—del Regolamento (UE) No. 1308/2013 (OCM Apicoltura), and by the Italian Ministry of Agriculture, Food Sovereignty, and Forestry (MASAF)—project ItaHoneyDNA—CUP J39I24001640006. KENJ has been supported by a Fulbright Fellowship.
Author information
Authors and Affiliations
Contributions
LF designed the study, interpreted the results, drafted the manuscript, provided funds and supervised the project. AR organised laboratory analyses, compiled data, coordinated the sampling and drafted the manuscript. SB performed bioinformatic analyses, curated the internal databases, carried out data analyses and drafted the manuscript. VT performed laboratory analyses and drafted the manuscript. KENJ performed laboratory analyses and revised the manuscript. AÖK performed laboratory analyses, collected samples and metadata. GS, VJU and FB provided metadata, tested primer pairs, and elaborated data. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Below is the link to the electronic supplementary material.
Supplementary Material 1 (download PDF )
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and permissions
About this article
Cite this article
Ribani, A., Bovo, S., Taurisano, V. et al. Honey environmental DNA reveals entomological fingerprints through dual mitochondrial cytochrome c oxidase subunit 1 (COI) and cytochrome b (CYTB) metabarcoding.
Sci Rep (2026). https://doi.org/10.1038/s41598-026-46493-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-026-46493-y
Keywords
- Apiculture
- Beekeeping
- DNA fingerprint
- eDNA
- Hemiptera
- Honeydew
Source: Ecology - nature.com
