1
|
Yang B, Zhou X, Liu S. Tracing the genealogy origin of geographic populations based on genomic variation and deep learning. Mol Phylogenet Evol 2024; 198:108142. [PMID: 38964594 DOI: 10.1016/j.ympev.2024.108142] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2023] [Revised: 05/30/2024] [Accepted: 07/01/2024] [Indexed: 07/06/2024]
Abstract
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.
Collapse
Affiliation(s)
- Bing Yang
- Department of Entomology, China Agricultural University, Beijing 100193, China
| | - Xin Zhou
- Department of Entomology, China Agricultural University, Beijing 100193, China.
| | - Shanlin Liu
- Department of Entomology, China Agricultural University, Beijing 100193, China; Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China.
| |
Collapse
|
2
|
Lefebre R, De Smet L, Tehel A, Paxton RJ, Bossuyt E, Verbeke W, van Dooremalen C, Ulgezen ZN, van den Bosch T, Schaafsma F, Valkenburg DJ, Dall’Olio R, Alaux C, Dezmirean DS, Giurgiu AI, Capela N, Simões S, Sousa JP, Bencsik M, McVeigh A, Ramsey MT, Ahmad S, Kumar T, Schäfer MO, Beaurepaire AL, Moro A, Flener CJ, Matthijs S, de Graaf DC. Allele Frequencies of Genetic Variants Associated with Varroa Drone Brood Resistance (DBR) in Apis mellifera Subspecies across the European Continent. INSECTS 2024; 15:419. [PMID: 38921134 PMCID: PMC11203681 DOI: 10.3390/insects15060419] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/30/2024] [Revised: 05/29/2024] [Accepted: 06/01/2024] [Indexed: 06/27/2024]
Abstract
Implementation of marker-assisted selection (MAS) in modern beekeeping would improve sustainability, especially in breeding programs aiming for resilience against the parasitic mite Varroa destructor. Selecting honey bee colonies for natural resistance traits, such as brood-intrinsic suppression of varroa mite reproduction, reduces the use of chemical acaricides while respecting local adaptation. In 2019, eight genomic variants associated with varroa non-reproduction in drone brood were discovered in a single colony from the Amsterdam Water Dune population in the Netherlands. Recently, a new study tested the applicability of these eight genetic variants for the same phenotype on a population-wide scale in Flanders, Belgium. As the properties of some variants varied between the two studies, one hypothesized that the difference in genetic ancestry of the sampled colonies may underly these contribution shifts. In order to frame this, we determined the allele frequencies of the eight genetic variants in more than 360 Apis mellifera colonies across the European continent and found that variant type allele frequencies of these variants are primarily related to the A. mellifera subspecies or phylogenetic honey bee lineage. Our results confirm that population-specific genetic markers should always be evaluated in a new population prior to using them in MAS programs.
Collapse
Affiliation(s)
- Regis Lefebre
- Department of Biochemistry and Microbiology, Ghent University, 9000 Ghent, Belgium; (L.D.S.); (E.B.); (D.C.d.G.)
| | - Lina De Smet
- Department of Biochemistry and Microbiology, Ghent University, 9000 Ghent, Belgium; (L.D.S.); (E.B.); (D.C.d.G.)
| | - Anja Tehel
- Institute for Biology, Martin Luther University Halle-Wittenberg, D-06126 Halle (Saale), Germany; (A.T.); (R.J.P.)
| | - Robert J. Paxton
- Institute for Biology, Martin Luther University Halle-Wittenberg, D-06126 Halle (Saale), Germany; (A.T.); (R.J.P.)
| | - Emma Bossuyt
- Department of Biochemistry and Microbiology, Ghent University, 9000 Ghent, Belgium; (L.D.S.); (E.B.); (D.C.d.G.)
| | - Wim Verbeke
- Department of Agricultural Economics, Ghent University, 9000 Ghent, Belgium;
| | - Coby van Dooremalen
- Wageningen Plant Research, Wageningen University & Research, 6708 PB Wageningen, The Netherlands (Z.N.U.); (T.v.d.B.); (F.S.); (D.-J.V.)
| | - Zeynep N. Ulgezen
- Wageningen Plant Research, Wageningen University & Research, 6708 PB Wageningen, The Netherlands (Z.N.U.); (T.v.d.B.); (F.S.); (D.-J.V.)
| | - Trudy van den Bosch
- Wageningen Plant Research, Wageningen University & Research, 6708 PB Wageningen, The Netherlands (Z.N.U.); (T.v.d.B.); (F.S.); (D.-J.V.)
| | - Famke Schaafsma
- Wageningen Plant Research, Wageningen University & Research, 6708 PB Wageningen, The Netherlands (Z.N.U.); (T.v.d.B.); (F.S.); (D.-J.V.)
| | - Dirk-Jan Valkenburg
- Wageningen Plant Research, Wageningen University & Research, 6708 PB Wageningen, The Netherlands (Z.N.U.); (T.v.d.B.); (F.S.); (D.-J.V.)
| | | | - Cedric Alaux
- Unité Abeilles et Environnement, Institut National de la Recherche pour l’Agriculture, CEDEX 9, 84914 Avignon, France;
| | - Daniel S. Dezmirean
- Apiculture and Sericulture Unit, University of Agricultural Sciences and Veterinary Medicine of Cluj-Napoca, 400372 Cluj Napoca, Romania; (D.S.D.); (A.I.G.)
| | - Alexandru I. Giurgiu
- Apiculture and Sericulture Unit, University of Agricultural Sciences and Veterinary Medicine of Cluj-Napoca, 400372 Cluj Napoca, Romania; (D.S.D.); (A.I.G.)
| | - Nuno Capela
- Centre for Functional Ecology, Associate Laboratory TERRA, Department of Life Sciences, University of Coimbra, 3000-456 Coimbra, Portugal; (N.C.); (S.S.); (J.P.S.)
| | - Sandra Simões
- Centre for Functional Ecology, Associate Laboratory TERRA, Department of Life Sciences, University of Coimbra, 3000-456 Coimbra, Portugal; (N.C.); (S.S.); (J.P.S.)
| | - José Paulo Sousa
- Centre for Functional Ecology, Associate Laboratory TERRA, Department of Life Sciences, University of Coimbra, 3000-456 Coimbra, Portugal; (N.C.); (S.S.); (J.P.S.)
| | - Martin Bencsik
- Physics and Mathematics, Nottingham Trent University, Nottingham NG8 11NS, UK; (M.B.); (A.M.); (M.T.R.); (S.A.); (T.K.)
| | - Adam McVeigh
- Physics and Mathematics, Nottingham Trent University, Nottingham NG8 11NS, UK; (M.B.); (A.M.); (M.T.R.); (S.A.); (T.K.)
| | - Michael Thomas Ramsey
- Physics and Mathematics, Nottingham Trent University, Nottingham NG8 11NS, UK; (M.B.); (A.M.); (M.T.R.); (S.A.); (T.K.)
| | - Sausan Ahmad
- Physics and Mathematics, Nottingham Trent University, Nottingham NG8 11NS, UK; (M.B.); (A.M.); (M.T.R.); (S.A.); (T.K.)
| | - Tarun Kumar
- Physics and Mathematics, Nottingham Trent University, Nottingham NG8 11NS, UK; (M.B.); (A.M.); (M.T.R.); (S.A.); (T.K.)
| | - Marc O. Schäfer
- Institute of Infectology, Federal Research Institute for Animal Health, Friedrich-Loeffler-Institut, 17493 Greifswald-Insel Riems, Germany;
| | - Alexis L. Beaurepaire
- Institute of Bee Health, University of Bern, 3003 Bern, Switzerland; (A.L.B.); (A.M.)
- Center for Bee Research, 3097 Liebefeld, Switzerland
| | - Arrigo Moro
- Institute of Bee Health, University of Bern, 3003 Bern, Switzerland; (A.L.B.); (A.M.)
- Galway Honey Bee Research Centre, Department of Zoology, University of Galway, H91 TK33 Galway, Ireland
| | | | - Severine Matthijs
- Department of Viral Reemerging, Enzootic and Bee Diseases, Sciensano, 1050 Elsene, Belgium;
| | - Dirk C. de Graaf
- Department of Biochemistry and Microbiology, Ghent University, 9000 Ghent, Belgium; (L.D.S.); (E.B.); (D.C.d.G.)
| |
Collapse
|
3
|
Blanchet G, Bellinger MR, Kearns AM, Cortes-Rodriguez N, Masuda B, Campana MG, Rutz C, Fleischer RC, Sutton JT. Reduction of genetic diversity in 'Alalā (Hawaiian crow; Corvus hawaiiensis) between the late 1800s and the late 1900s. J Hered 2024; 115:32-44. [PMID: 37846510 DOI: 10.1093/jhered/esad063] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2023] [Revised: 09/26/2023] [Accepted: 10/12/2023] [Indexed: 10/18/2023] Open
Abstract
Genetic and genomic data are increasingly used to aid conservation management of endangered species by providing insights into evolutionary histories, factors associated with extinction risks, and potential for future adaptation. For the 'Alalā, or Hawaiian crow (Corvus hawaiiensis), genetic concerns include negative correlations between inbreeding and hatching success. However, it is unclear if low genetic diversity and inbreeding depression are consequences of a historical population bottleneck, or if 'Alalā had historically low genetic diversity that predated human influence, perhaps as a result of earlier declines or founding events. In this study, we applied a hybridization-based sequence capture to generate a genome-wide single nucleotide polymorphism (SNP) dataset for comparing historical specimens collected in the 1890s, when 'Alalā were more numerous, to samples taken between 1973 and 1998, when 'Alalā population densities were near the lowest documented levels in the wild, prior to all individuals being collected for captive rearing. We found low genome-wide diversity in both sample groups, however, the modern sample group (1973 to 1998 cohort) exhibited relatively fewer polymorphic alleles, a lower proportion of polymorphic loci, and lower observed heterozygosity, consistent with a population decline and potential bottleneck effects. These results combined with a current low population size highlight the importance of continued efforts by conservation managers to mitigate inbreeding and maintain founder representation to preserve what genetic diversity remains.
Collapse
Affiliation(s)
- Geneviève Blanchet
- Department of Biology, University of Hawai'i at Hilo, 200 W Kāwili St, Hilo, Hawai'i 96720, United States
| | - M Renee Bellinger
- Department of Biology, University of Hawai'i at Hilo, 200 W Kāwili St, Hilo, Hawai'i 96720, United States
- U.S. Geological Survey, Pacific Island Ecosystems Research Center, PO Box 44, Hawai'i National Park, Hawai'i 96718, United States
| | - Anna M Kearns
- Center for Conservation Genomics, National Zoo and Conservation Biology Institute, Smithsonian Institution, Washington DC 20008, United States
| | - Nandadevi Cortes-Rodriguez
- Center for Conservation Genomics, National Zoo and Conservation Biology Institute, Smithsonian Institution, Washington DC 20008, United States
| | - Bryce Masuda
- San Diego Zoo Wildlife Alliance, P.O. Box 39, Volcano, HI 96785, United States
| | - Michael G Campana
- Center for Conservation Genomics, National Zoo and Conservation Biology Institute, Smithsonian Institution, Washington DC 20008, United States
| | - Christian Rutz
- Centre for Biological Diversity, School of Biology, University of St Andrews, St Andrews KY16 9TH, United Kingdom
| | - Robert C Fleischer
- Center for Conservation Genomics, National Zoo and Conservation Biology Institute, Smithsonian Institution, Washington DC 20008, United States
| | - Jolene T Sutton
- Department of Biology, University of Hawai'i at Hilo, 200 W Kāwili St, Hilo, Hawai'i 96720, United States
| |
Collapse
|
4
|
Cornman RS. Data mining reveals tissue-specific expression and host lineage-associated forms of Apis mellifera filamentous virus. PeerJ 2023; 11:e16455. [PMID: 38025724 PMCID: PMC10655722 DOI: 10.7717/peerj.16455] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2023] [Accepted: 10/23/2023] [Indexed: 12/01/2023] Open
Abstract
Background Apis mellifera filamentous virus (AmFV) is a large double-stranded DNA virus of uncertain phylogenetic position that infects honey bees (Apis mellifera). Little is known about AmFV evolution or molecular aspects of infection. Accurate annotation of open-reading frames (ORFs) is challenged by weak homology to other known viruses. This study was undertaken to evaluate ORFs (including coding-frame conservation, codon bias, and purifying selection), quantify genetic variation within AmFV, identify host characteristics that covary with infection rate, and examine viral expression patterns in different tissues. Methods Short-read data were accessed from the Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI). Sequence reads were downloaded from accessions meeting search criteria and scanned for kmers representative of AmFV genomic sequence. Samples with kmer counts above specified thresholds were downloaded in full for mapping to reference sequences and de novo assembly. Results At least three distinct evolutionary lineages of AmFV exist. Clade 1 predominates in Europe but in the Americas and Africa it is replaced by the other clades as infection level increases in hosts. Only clade 3 was found at high relative abundance in hosts with African ancestry, whereas all clades achieved high relative abundance in bees of non-African ancestry. In Europe and Africa, clade 2 was generally detected only in low-level infections but was locally dominant in some North American samples. The geographic distribution of clade 3 was consistent with an introduction to the Americas with 'Africanized' honey bees in the 1950s. Localized genomic regions of very high nucleotide divergence in individual isolates suggest recombination with additional, as-yet unidentified AmFV lineages. A set of 155 high-confidence ORFs was annotated based on evolutionary conservation in six AmFV genome sequences representative of the three clades. Pairwise protein-level identity averaged 94.6% across ORFs (range 77.1-100%), which generally exhibited low evolutionary rates and moderate to strong codon bias. However, no robust example of positive diversifying selection on coding sequence was found in these alignments. Most of the genome was detected in RNA short-read alignments. Transcriptome assembly often yielded contigs in excess of 50 kb and containing ORFs in both orientations, and the termini of long transcripts were associated with tandem repeats. Lower levels of AmFV RNA were detected in brain tissue compared to abdominal tissue, and a distinct set of ORFs had minimal to no detectable expression in brain tissue. A scan of DNA accessions from the parasitic mite Varroa destructor was inconclusive with respect to replication in that species. Discussion Collectively, these results expand our understanding of this enigmatic virus, revealing transcriptional complexity and co-evolutionary associations with host lineage.
Collapse
|
5
|
Kaskinova MD, Gaifullina LR, Saltykova ES. Haplotypes of the tRNAleu-COII mtDNA Region in Russian Apis mellifera Populations. Animals (Basel) 2023; 13:2394. [PMID: 37508171 PMCID: PMC10376158 DOI: 10.3390/ani13142394] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2023] [Revised: 07/19/2023] [Accepted: 07/22/2023] [Indexed: 07/30/2023] Open
Abstract
Analysis of the mtDNA tRNAleu-COII locus is a widely used tool to establish belonging to a particular evolutionary lineage of Apis mellifera L. (lineages A, M, C, O, and Y). In Russia, most of the area was once inhabited by Apis mellifera mellifera from the M evolutionary lineage, but the introduction of bee subspecies from the southern regions of Russia (A. m. caucasica, A. m. carnica) and from abroad (A. m. carnica, A. m. ligustica) led to fragmentation of their native range. In this study, the results of assessing the haplotype number for the tRNAleu-COII locus of mtDNA in Russian Apis mellifera populations were presented. We analyzed 269 colonies from 19 regions of Russia. As a result, two evolutionary lineages were identified: the East European lineage C (26.4%) and the Northwestern European lineage M (73.6%). A total of 29 haplotypes were identified, 8 of them were already reported, and 21 were found to be novel. From the C lineage, haplotypes C1, C2, C2c, C2j, and C3 were predominant. All M lineage samples from Russia belong to the M17 and M4' haplogroups but have only minor variations in the form of nucleotide substitutions. An analysis of publications devoted to the tRNAleu-COII locus haplotypes, as well as an analysis of the available tRNAleu-COII sequences in GenBank, showed that there is still a problem with the haplotype nomenclature.
Collapse
Affiliation(s)
- Milyausha D Kaskinova
- Institute of Biochemistry and Genetics, Ufa Federal Research Center, Russian Academy of Sciences, Prospekt Oktyabrya 71, 450054 Ufa, Russia
| | - Luisa R Gaifullina
- Institute of Biochemistry and Genetics, Ufa Federal Research Center, Russian Academy of Sciences, Prospekt Oktyabrya 71, 450054 Ufa, Russia
| | - Elena S Saltykova
- Institute of Biochemistry and Genetics, Ufa Federal Research Center, Russian Academy of Sciences, Prospekt Oktyabrya 71, 450054 Ufa, Russia
| |
Collapse
|
6
|
Carr SM. Multiple mitogenomes indicate Things Fall Apart with Out of Africa or Asia hypotheses for the phylogeographic evolution of Honey Bees (Apis mellifera). Sci Rep 2023; 13:9386. [PMID: 37296293 PMCID: PMC10256785 DOI: 10.1038/s41598-023-35937-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2023] [Accepted: 05/26/2023] [Indexed: 06/12/2023] Open
Abstract
Previous morpho-molecular studies of evolutionary relationships within the economically important genus of honey bees (Apis), including the Western Honey Bee (A. mellifera L.), have suggested Out of Africa or Asia origins and subsequent spread to Europe. I test these hypotheses by a meta-analysis of complete mitochondrial DNA coding regions (11.0 kbp) from 22 nominal subspecies represented by 78 individual sequences in A. mellifera. Parsimony, distance, and likelihood analyses identify six nested clades: Things Fall Apart with Out of Africa or Asia hypotheses. Molecular clock-calibrated phylogeographic analysis shows instead a basal origin of A. m. mellifera in Europe ~ 780 Kya, and expansion to Southeast Europe and Asia Minor ~ 720 Kya. Eurasian bees spread southward via a Levantine/Nilotic/Arabian corridor into Africa ~ 540 Kya. An African clade re-established in Iberia ~ 100 Kya spread thereafter to westerly Mediterranean islands and back into North Africa. Nominal subspecies within the Asia Minor and Mediterranean clades are less differentiated than are individuals within other subspecies. Names matter: paraphyletic anomalies are artefacts of mis-referral in GenBank of sequences to the wrong subspecies, or use of faulty sequences, which are clarified by inclusion of multiple sequences from available subspecies.
Collapse
Affiliation(s)
- Steven M Carr
- Genetics, Evolution, and Molecular Systematics Laboratory, Department of Biology, Memorial University of Newfoundland, St John's, NL, A1C 5S9, Canada.
| |
Collapse
|
7
|
Parejo M, Talenti A, Richardson M, Vignal A, Barnett M, Wragg D. AmelHap: Leveraging drone whole-genome sequence data to create a honey bee HapMap. Sci Data 2023; 10:198. [PMID: 37037860 PMCID: PMC10086014 DOI: 10.1038/s41597-023-02097-z] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2022] [Accepted: 03/22/2023] [Indexed: 04/12/2023] Open
Abstract
Honey bee, Apis mellifera, drones are typically haploid, developing from an unfertilized egg, inheriting only their queen's alleles and none from the many drones she mated with. Thus the ordered combination or 'phase' of alleles is known, making drones a valuable haplotype resource. We collated whole-genome sequence data for 1,407 drones, including 45 newly sequenced Scottish drones, collectively representing 19 countries, 8 subspecies and various hybrids. Following alignment to Amel_HAv3.1, variant calling and quality filtering, we retained 17.4 M high quality variants across 1,328 samples with a genotyping rate of 98.7%. We demonstrate the utility of this haplotype resource, AmelHap, for genotype imputation, returning >95% concordance when up to 61% of data is missing in haploids and up to 12% of data is missing in diploids. AmelHap will serve as a useful resource for the community for imputation from low-depth sequencing or SNP chip data, accurate phasing of diploids for association studies, and as a comprehensive reference panel for population genetic and evolutionary analyses.
Collapse
Affiliation(s)
- M Parejo
- Applied Genomics and Bioinformatics, University of the Basque Country (UPV/EHU), Leioa, Spain
| | - A Talenti
- The Roslin Institute, University of Edinburgh, Easter Bush Campus, Midlothian, UK
| | - M Richardson
- University of Edinburgh, King's Buildings Campus, Edinburgh, UK
- Beebytes Analytics CIC, Roslin Innovation Centre, Easter Bush Campus, Midlothian, UK
| | - A Vignal
- GenPhySE, Université de Toulouse, INRAE, INPT, INP-ENVT, 31326, Castanet Tolosan, France
| | - M Barnett
- Beebytes Analytics CIC, Roslin Innovation Centre, Easter Bush Campus, Midlothian, UK
| | - D Wragg
- The Roslin Institute, University of Edinburgh, Easter Bush Campus, Midlothian, UK.
- Beebytes Analytics CIC, Roslin Innovation Centre, Easter Bush Campus, Midlothian, UK.
| |
Collapse
|
8
|
Oleksa A, Căuia E, Siceanu A, Puškadija Z, Kovačić M, Pinto MA, Rodrigues PJ, Hatjina F, Charistos L, Bouga M, Prešern J, Kandemir İ, Rašić S, Kusza S, Tofilski A. Honey bee (Apis mellifera) wing images: a tool for identification and conservation. Gigascience 2023; 12:giad019. [PMID: 36971293 PMCID: PMC10041535 DOI: 10.1093/gigascience/giad019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2022] [Revised: 02/02/2023] [Accepted: 03/06/2023] [Indexed: 03/29/2023] Open
Abstract
BACKGROUND The honey bee (Apis mellifera) is an ecologically and economically important species that provides pollination services to natural and agricultural systems. The biodiversity of the honey bee in parts of its native range is endangered by migratory beekeeping and commercial breeding. In consequence, some honey bee populations that are well adapted to the local environment are threatened with extinction. A crucial step for the protection of honey bee biodiversity is reliable differentiation between native and nonnative bees. One of the methods that can be used for this is the geometric morphometrics of wings. This method is fast, is low cost, and does not require expensive equipment. Therefore, it can be easily used by both scientists and beekeepers. However, wing geometric morphometrics is challenging due to the lack of reference data that can be reliably used for comparisons between different geographic regions. FINDINGS Here, we provide an unprecedented collection of 26,481 honey bee wing images representing 1,725 samples from 13 European countries. The wing images are accompanied by the coordinates of 19 landmarks and the geographic coordinates of the sampling locations. We present an R script that describes the workflow for analyzing the data and identifying an unknown sample. We compared the data with available reference samples for lineage and found general agreement with them. CONCLUSIONS The extensive collection of wing images available on the Zenodo website can be used to identify the geographic origin of unknown samples and therefore assist in the monitoring and conservation of honey bee biodiversity in Europe.
Collapse
Affiliation(s)
- Andrzej Oleksa
- Department of Genetics, Faculty of Biological Sciences, Kazimierz Wielki University, Bydgoszcz 85-090, Poland
| | - Eliza Căuia
- Honeybee Genetics and Breeding Laboratory, Institute for Beekeeping Research and Development, Bucharest 013975, Romania
| | - Adrian Siceanu
- Honeybee Genetics and Breeding Laboratory, Institute for Beekeeping Research and Development, Bucharest 013975, Romania
| | - Zlatko Puškadija
- Faculty of Agrobiotechnical Sciences, Josip Juraj Strossmayer University of Osijek, Osijek 31000, Croatia
| | - Marin Kovačić
- Faculty of Agrobiotechnical Sciences, Josip Juraj Strossmayer University of Osijek, Osijek 31000, Croatia
| | - M Alice Pinto
- Centro de Investigação de Montanha, Instituto Politécnico de Bragança, Campus de Santa Apolónia, Bragança 5300-253, Portugal
- Laboratório Associado para a Sustentabilidade e Tecnologia em Regiões de Montanha (SusTEC), Instituto Politécnico de Bragança, Campus de Santa Apolónia, Bragança 5300-253, Portugal
| | - Pedro João Rodrigues
- Centre in Digitalization and Intelligent Robotics, Instituto Politécnico de Bragança, Campus de Santa Apolónia, Bragança 5300-253, Portugal
- Laboratório Associado para a Sustentabilidade e Tecnologia em Regiões de Montanha (SusTEC), Instituto Politécnico de Bragança, Campus de Santa Apolónia, Bragança 5300-253, Portugal
| | - Fani Hatjina
- Department of Apiculture, Institute of Animal Science–Ellinikos Georgikos Organismos ‘DIMITRA’, Nea Moudania 63200, Greece
| | - Leonidas Charistos
- Department of Apiculture, Institute of Animal Science–Ellinikos Georgikos Organismos ‘DIMITRA’, Nea Moudania 63200, Greece
| | - Maria Bouga
- Lab of Agricultural Zoology and Entomology, Agricultural University of Athens, Athens 11855, Greece
| | - Janez Prešern
- Agricultural Institute of Slovenia, Ljubljana SI-1000, Slovenia
| | - İrfan Kandemir
- Ankara University, Department of Biology, Faculty of Science, Ankara University, Beşevler-Ankara 06100, Turkey
| | - Slađan Rašić
- Faculty of Ecological Agriculture, EDUCONS University, Sremska Kamenica 21208, Serbia
| | - Szilvia Kusza
- Centre for Agricultural Genomics and Biotechnology, University of Debrecen, Debrecen 4032, Hungary
| | - Adam Tofilski
- Department of Zoology and Animal Welfare, University of Agriculture in Krakow, Krakow 31-425, Poland
| |
Collapse
|
9
|
Gmel AI, Guichard M, Dainat B, Williams GR, Eynard S, Vignal A, Servin B, Neuditschko M. Identification of runs of homozygosity in Western honey bees ( Apis mellifera) using whole-genome sequencing data. Ecol Evol 2023; 13:e9723. [PMID: 36694553 PMCID: PMC9843643 DOI: 10.1002/ece3.9723] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2021] [Revised: 12/15/2022] [Accepted: 12/19/2022] [Indexed: 01/19/2023] Open
Abstract
Runs of homozygosity (ROH) are continuous homozygous segments that arise through the transmission of haplotypes that are identical by descent. The length and distribution of ROH segments provide insights into the genetic diversity of populations and can be associated with selection signatures. Here, we analyzed reconstructed whole-genome queen genotypes, from a pool-seq data experiment including 265 Western honeybee colonies from Apis mellifera mellifera and Apis mellifera carnica. Integrating individual ROH patterns and admixture levels in a dynamic population network visualization allowed us to ascertain major differences between the two subspecies. Within A. m. mellifera, we identified well-defined substructures according to the genetic origin of the queens. Despite the current applied conservation efforts, we pinpointed 79 admixed queens. Genomic inbreeding (F ROH) strongly varied within and between the identified subpopulations. Conserved A. m. mellifera from Switzerland had the highest mean F ROH (3.39%), while queens originating from a conservation area in France, which were also highly admixed, showed significantly lower F ROH (0.45%). The majority of A. m. carnica queens were also highly admixed, except 12 purebred queens with a mean F ROH of 2.33%. Within the breed-specific ROH islands, we identified 14 coding genes for A. m. mellifera and five for A. m. carnica, respectively. Local adaption of A. m. mellifera could be suggested by the identification of genes involved in the response to ultraviolet light (Crh-BP, Uvop) and body size (Hex70a, Hex70b), while the A. m. carnica specific genes Cpr3 and Cpr4 are most likely associated with the lighter striping pattern, a morphological phenotype expected in this subspecies. We demonstrated that queen genotypes derived from pooled workers are useful tool to unravel the population dynamics in A. mellifera and provide fundamental information to conserve native honey bees.
Collapse
Affiliation(s)
- Annik Imogen Gmel
- Animal GenoPhenomics, Animal Production Systems and Animal HealthAgroscopePosieuxSwitzerland
| | - Matthieu Guichard
- Animal GenoPhenomics, Animal Production Systems and Animal HealthAgroscopePosieuxSwitzerland
- Swiss Bee Research CentreAgroscopeLiebefeldSwitzerland
| | | | | | - Sonia Eynard
- GenPhySEINRAE, INPT, INPENVTUniversité de ToulouseCastanet‐TolosanFrance
- UMT PrADEProtection des Abeilles Dans L'EnvironnementAvignonFrance
| | - Alain Vignal
- GenPhySEINRAE, INPT, INPENVTUniversité de ToulouseCastanet‐TolosanFrance
- UMT PrADEProtection des Abeilles Dans L'EnvironnementAvignonFrance
| | - Bertrand Servin
- GenPhySEINRAE, INPT, INPENVTUniversité de ToulouseCastanet‐TolosanFrance
- UMT PrADEProtection des Abeilles Dans L'EnvironnementAvignonFrance
| | | | - Markus Neuditschko
- Animal GenoPhenomics, Animal Production Systems and Animal HealthAgroscopePosieuxSwitzerland
| |
Collapse
|
10
|
Kaskinova M, Gaifullina L, Ilyasov R, Lelej A, Kwon HW, Thai PH, Saltykova E. Genetic Structure of Apis cerana Populations from South Korea, Vietnam and the Russian Far East Based on Microsatellite and Mitochondrial DNA Polymorphism. INSECTS 2022; 13:1174. [PMID: 36555084 PMCID: PMC9784541 DOI: 10.3390/insects13121174] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/26/2022] [Revised: 11/29/2022] [Accepted: 12/16/2022] [Indexed: 06/17/2023]
Abstract
In this article, we present the results of the genetic analysis of Apis cerana samples from the Russian Far East, South Korea and Vietnam. An analysis of the polymorphism of seven microsatellite loci and an assessment of the haplotype diversity of the mtDNA tRNAleu-COII locus were performed. A fragment of about 431 bp in tRNAleu-COII was sequenced. The analysis showed the presence of 14 haplotypes, while the predominant haplotype was Japan1. Microsatellite data revealed two differentiated clusters. The first cluster contained tropical climate A. cerana samples from Vietnam, and the second cluster combined temperate climate A. cerana samples from the Russian Far East and South Korea.
Collapse
Affiliation(s)
- Milyausha Kaskinova
- Institute of Biochemistry and Genetics, Ufa Federal Research Center, Russian Academy of Sciences, Prospekt Oktyabrya 71, 450054 Ufa, Russia
| | - Luisa Gaifullina
- Institute of Biochemistry and Genetics, Ufa Federal Research Center, Russian Academy of Sciences, Prospekt Oktyabrya 71, 450054 Ufa, Russia
| | - Rustem Ilyasov
- Scientific and Educational Center, Bashkir State Agrarian University, 50-Letiya Oktyabrya Str. 34, 450001 Ufa, Russia
- Department of Genetics and Biotechnology, Vavilov Institute of General Genetics, Russian Academy of Sciences, Gubkina Str. 3, 119333 Moscow, Russia
| | - Arkady Lelej
- Federal Scientific Center of the East Asia Terrestrial Biodiversity, Prospekt 100-let Vladivostoka, 159, 690022 Vladivostok, Russia
| | - Hyung Wook Kwon
- Department of Life Sciences and Convergence Research Center for Insect Vectors, Incheon National University, 119 Academy-ro, Yeonsu-gu, Incheon 22012, Republic of Korea
| | - Pham Hong Thai
- Research Center for Tropical Bees and Beekeeping, Vietnam National University of Agriculture, Trau Quy, Gia Lam, Hanoi 10000, Vietnam
| | - Elena Saltykova
- Institute of Biochemistry and Genetics, Ufa Federal Research Center, Russian Academy of Sciences, Prospekt Oktyabrya 71, 450054 Ufa, Russia
| |
Collapse
|
11
|
García CAY, Rodrigues PJ, Tofilski A, Elen D, McCormak GP, Oleksa A, Henriques D, Ilyasov R, Kartashev A, Bargain C, Fried B, Pinto MA. Using the Software DeepWings© to Classify Honey Bees across Europe through Wing Geometric Morphometrics. INSECTS 2022; 13:1132. [PMID: 36555043 PMCID: PMC9784756 DOI: 10.3390/insects13121132] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/02/2022] [Revised: 11/22/2022] [Accepted: 11/30/2022] [Indexed: 06/17/2023]
Abstract
DeepWings© is a software that uses machine learning to automatically classify honey bee subspecies by wing geometric morphometrics. Here, we tested the five subspecies classifier (A. m. carnica, Apis mellifera caucasia, A. m. iberiensis, Apis mellifera ligustica, and A. m. mellifera) of DeepWings© on 14,816 wing images with variable quality and acquired by different beekeepers and researchers. These images represented 2601 colonies from the native ranges of the M-lineage A. m. iberiensis and A. m. mellifera, and the C-lineage A. m. carnica. In the A. m. iberiensis range, 92.6% of the colonies matched this subspecies, with a high median probability (0.919). In the Azores, where the Iberian subspecies was historically introduced, a lower proportion (85.7%) and probability (0.842) were observed. In the A. m mellifera range, only 41.1 % of the colonies matched this subspecies, which is compatible with a history of C-derived introgression. Yet, these colonies were classified with the highest probability (0.994) of the three subspecies. In the A. m. carnica range, 88.3% of the colonies matched this subspecies, with a probability of 0.984. The association between wing and molecular markers, assessed for 1214 colonies from the M-lineage range, was highly significant but not strong (r = 0.31, p < 0.0001). The agreement between the markers was influenced by C-derived introgression, with the best results obtained for colonies with high genetic integrity. This study indicates the good performance of DeepWings© on a realistic wing image dataset.
Collapse
Affiliation(s)
- Carlos Ariel Yadró García
- Centro de Investigação de Montanha, Instituto Politécnico de Bragança, Campus de Santa Apolónia, 5300-253 Bragança, Portugal
- Laboratório Associado para a Sustentabilidade e Tecnologia em Regiões de Montanha (SusTEC), Instituto Politécnico de Bragança, Campus de Santa Apolónia, 5300-253 Bragança, Portugal
| | - Pedro João Rodrigues
- Laboratório Associado para a Sustentabilidade e Tecnologia em Regiões de Montanha (SusTEC), Instituto Politécnico de Bragança, Campus de Santa Apolónia, 5300-253 Bragança, Portugal
- Research Center in Digitalization and Intelligent Robotics (CeDRI), Instituto Politécnico de Bragança, Campus de Santa Apolónia, 5300-253 Bragança, Portugal
| | - Adam Tofilski
- Department of Zoology and Animal Welfare, University of Agriculture in Krakow, 31-425 Krakow, Poland
| | - Dylan Elen
- Department of Molecular Ecology & Evolution, School of Natural Sciences, Bangor University, Bangor, Gwynedd LL57 2DG, UK
- Taskforce Research, ZwarteBij.org VZW, 9890 Gavere, Belgium
| | - Grace P. McCormak
- Zoology, Earth and Life Sciences, School of Natural Sciences, University of Galway, H91 TK33 Galway, Ireland
| | - Andrzej Oleksa
- Department of Genetics, Faculty of Biological Sciences, Kazimierz Wielki University, Powstańców Wielkopolskich, 85-090 Bydgoszcz, Poland
| | - Dora Henriques
- Centro de Investigação de Montanha, Instituto Politécnico de Bragança, Campus de Santa Apolónia, 5300-253 Bragança, Portugal
- Laboratório Associado para a Sustentabilidade e Tecnologia em Regiões de Montanha (SusTEC), Instituto Politécnico de Bragança, Campus de Santa Apolónia, 5300-253 Bragança, Portugal
| | - Rustem Ilyasov
- Transplantology and Genotyping, Department of the Progressive Technologies, Bashkir State Agrarian University, 450001 Ufa, Russia
| | | | - Christian Bargain
- Association pour la Sauvegarde de l’Abeillee Noire, 56069 Ile de Groix, France
| | - Balser Fried
- Swiss Association of Mellifera Bee Friends, mellifera.ch, Ahornstrasse 7, 9533 Kirchberg, Switzerland
| | - Maria Alice Pinto
- Centro de Investigação de Montanha, Instituto Politécnico de Bragança, Campus de Santa Apolónia, 5300-253 Bragança, Portugal
- Laboratório Associado para a Sustentabilidade e Tecnologia em Regiões de Montanha (SusTEC), Instituto Politécnico de Bragança, Campus de Santa Apolónia, 5300-253 Bragança, Portugal
| |
Collapse
|
12
|
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, Volume 121, March 2021, 107626”. Food Control 2022. [DOI: 10.1016/j.foodcont.2022.109571] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
13
|
Fontanesi L, Taurisano V, Ribani A, Utzeri VJ. 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, Volume 121, March 2021, 107626”. Food Control 2022. [DOI: 10.1016/j.foodcont.2022.109570] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
14
|
Bovo S, Utzeri VJ, Ribani A, Taurisano V, Schiavo G, Fontanesi L. A genotyping by sequencing approach can disclose Apis mellifera population genomic information contained in honey environmental DNA. Sci Rep 2022; 12:19541. [PMID: 36379985 PMCID: PMC9666642 DOI: 10.1038/s41598-022-24101-z] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2022] [Accepted: 11/10/2022] [Indexed: 11/16/2022] Open
Abstract
Awareness has been raised over the last years on the genetic integrity of autochthonous honey bee subspecies. Genomic tools available in Apis mellifera can make it possible to measure this information by targeting individual honey bee DNA. Honey contains DNA traces from all organisms that contributed or were involved in its production steps, including the honey bees of the colony. In this study, we designed and tested a genotyping by sequencing (GBS) assay to analyse single nucleotide polymorphisms (SNPs) of A. mellifera nuclear genome using environmental DNA extracted from honey. A total of 121 SNPs (97 SNPs informative for honey bee subspecies identification and 24 SNPs associated with relevant traits of the colonies) were used in the assay to genotype honey DNA, which derives from thousands of honey bees. Results were integrated with information derived from previous studies and whole genome resequencing datasets. This GBS method is highly reliable in estimating honey bee SNP allele frequencies of the whole colony from which the honey derived. This assay can be used to identify the honey bee subspecies of the colony that produced the honey and, in turn, to authenticate the entomological origin of the honey.
Collapse
Affiliation(s)
- Samuele Bovo
- grid.6292.f0000 0004 1757 1758Division of Animal Sciences, Department of Agricultural and Food Sciences, University of Bologna, Viale Giuseppe Fanin 46, 40127 Bologna, Italy
| | - Valerio Joe Utzeri
- grid.6292.f0000 0004 1757 1758Division of Animal Sciences, Department of Agricultural and Food Sciences, University of Bologna, Viale Giuseppe Fanin 46, 40127 Bologna, Italy
| | - Anisa Ribani
- grid.6292.f0000 0004 1757 1758Division of Animal Sciences, Department of Agricultural and Food Sciences, University of Bologna, Viale Giuseppe Fanin 46, 40127 Bologna, Italy
| | - Valeria Taurisano
- grid.6292.f0000 0004 1757 1758Division of Animal Sciences, Department of Agricultural and Food Sciences, University of Bologna, Viale Giuseppe Fanin 46, 40127 Bologna, Italy
| | - Giuseppina Schiavo
- grid.6292.f0000 0004 1757 1758Division of Animal Sciences, Department of Agricultural and Food Sciences, University of Bologna, Viale Giuseppe Fanin 46, 40127 Bologna, Italy
| | - Luca Fontanesi
- grid.6292.f0000 0004 1757 1758Division of Animal Sciences, Department of Agricultural and Food Sciences, University of Bologna, Viale Giuseppe Fanin 46, 40127 Bologna, Italy
| |
Collapse
|
15
|
Mehdizadeh SA, Abdolahzare Z, Karaji FK, Mouazen A. Design and manufacturing a microcontroller based measurement device for honey adulteration detection. J Food Compost Anal 2022. [DOI: 10.1016/j.jfca.2022.105049] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
16
|
Kaskinova MD, Gaifullina LR, Kinzikeev AK, Gladkih AN, Sokolyanskaya MP, Saltykova ES. Genetic Structure of the Apis mellifera L. Population from Altai Krai. RUSS J GENET+ 2022. [DOI: 10.1134/s1022795422070092] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
17
|
DeepWings©: Automatic Wing Geometric Morphometrics Classification of Honey Bee (Apis mellifera) Subspecies Using Deep Learning for Detecting Landmarks. BIG DATA AND COGNITIVE COMPUTING 2022. [DOI: 10.3390/bdcc6030070] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Honey bee classification by wing geometric morphometrics entails the first step of manual annotation of 19 landmarks in the forewing vein junctions. This is a time-consuming and error-prone endeavor, with implications for classification accuracy. Herein, we developed a software called DeepWings© that overcomes this constraint in wing geometric morphometrics classification by automatically detecting the 19 landmarks on digital images of the right forewing. We used a database containing 7634 forewing images, including 1864 analyzed by F. Ruttner in the original delineation of 26 honey bee subspecies, to tune a convolutional neural network as a wing detector, a deep learning U-Net as a landmarks segmenter, and a support vector machine as a subspecies classifier. The implemented MobileNet wing detector was able to achieve a mAP of 0.975 and the landmarks segmenter was able to detect the 19 landmarks with 91.8% accuracy, with an average positional precision of 0.943 resemblance to manually annotated landmarks. The subspecies classifier, in turn, presented an average accuracy of 86.6% for 26 subspecies and 95.8% for a subset of five important subspecies. The final implementation of the system showed good speed performance, requiring only 14 s to process 10 images. DeepWings© is very user-friendly and is the first fully automated software, offered as a free Web service, for honey bee classification from wing geometric morphometrics. DeepWings© can be used for honey bee breeding, conservation, and even scientific purposes as it provides the coordinates of the landmarks in excel format, facilitating the work of research teams using classical identification approaches and alternative analytical tools.
Collapse
|
18
|
Marcelino J, Braese C, Christmon K, Evans JD, Gilligan T, Giray T, Nearman A, Niño EL, Rose R, Sheppard WS, vanEngelsdorp D, Ellis JD. The Movement of Western Honey Bees (Apis mellifera L.) Among United States and Territories: History, Benefits, Risks, and Mitigation Strategies. Front Ecol Evol 2022. [DOI: 10.3389/fevo.2022.850600] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Beekeeping is a cornerstone activity that has led to the human-mediated, global spread of western honey bees (Apis mellifera L.) outside their native range of Europe, western Asia, and Africa. The exportation/importation of honey bees (i.e., transfer of honey bees or germplasm between countries) is regulated at the national level in many countries. Honey bees were first imported into the United States in the early 1600’s. Today, honey bee movement (i.e., transport of honey bees among states and territories) is regulated within the United States at the state, territory, and federal levels. At the federal level, honey bees present in the country (in any state or territory) can be moved among states and territories without federal restriction, with the exception of movement to Hawaii. In contrast, regulations at the state and territory levels vary substantially, ranging from no additional regulations beyond those stipulated at the federal level, to strict regulations for the introduction of live colonies, packaged bees, or queens. This variability can lead to inconsistencies in the application of regulations regarding the movement of honey bees among states and territories. In November 2020, we convened a technical working group (TWG), composed of academic and USDA personnel, to review and summarize the (1) history of honey bee importation into/movement within the United States, (2) current regulations regarding honey bee movement and case studies on the application of those regulations, (3) benefits associated with moving honey bees within the United States, (4) risks associated with moving honey bees within the United States, and (5) risk mitigation strategies. This review will be helpful for developing standardized best practices for the safe movement of honey bees between the 48 contiguous states and other states/territories within the United States.
Collapse
|
19
|
Du M, Bernstein R, Hoppe A, Bienefeld K. Influence of model selection and data structure on the estimation of genetic parameters in honeybee populations. G3 (BETHESDA, MD.) 2022; 12:6500294. [PMID: 35100384 PMCID: PMC8824827 DOI: 10.1093/g3journal/jkab450] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/16/2021] [Accepted: 12/22/2021] [Indexed: 11/12/2022]
Abstract
Estimating genetic parameters of quantitative traits is a prerequisite for animal breeding. In honeybees, the genetic variance separates into queen and worker effects. However, under data paucity, parameter estimations that account for this peculiarity often yield implausible results. Consequently, simplified models that attribute all genetic contributions to either the queen (queen model) or the workers (worker model) are often used to estimate variance components in honeybees. However, the causes for estimations with the complete model (colony model) to fail and the consequences of simplified models for variance estimates are little understood. We newly developed the necessary theory to compare parameter estimates that were achieved by the colony model with those of the queen and worker models. Furthermore, we performed computer simulations to quantify the influence of model choice, estimation algorithm, true genetic parameters, rates of controlled mating, apiary sizes, and phenotype data completeness on the success of genetic parameter estimations. We found that successful estimations with the colony model were only possible if at least some of the queens mated controlled on mating stations. In that case, estimates were largely unbiased if more than 20% of the colonies had phenotype records. The simplified queen and worker models proved more stable and yielded plausible parameter estimates for almost all settings. Results obtained from these models were unbiased when mating was uncontrolled, but with controlled mating, the simplified models consistently overestimated heritabilities. This study elucidates the requirements for variance component estimation in honeybees and provides the theoretical groundwork for simplified honeybee models.
Collapse
Affiliation(s)
- Manuel Du
- Breeding and Behavior, Institute for Bee Research Hohen Neuendorf, 16540 Hohen Neuendorf, Germany
| | - Richard Bernstein
- Breeding and Behavior, Institute for Bee Research Hohen Neuendorf, 16540 Hohen Neuendorf, Germany
| | - Andreas Hoppe
- Breeding and Behavior, Institute for Bee Research Hohen Neuendorf, 16540 Hohen Neuendorf, Germany
| | - Kaspar Bienefeld
- Breeding and Behavior, Institute for Bee Research Hohen Neuendorf, 16540 Hohen Neuendorf, Germany
| |
Collapse
|
20
|
Chen C, Parejo M, Momeni J, Langa J, Nielsen RO, Shi W, Vingborg R, Kryger P, Bouga M, Estonba A, Meixner M. Population Structure and Diversity in European Honey Bees (Apis mellifera L.)—An Empirical Comparison of Pool and Individual Whole-Genome Sequencing. Genes (Basel) 2022; 13:genes13020182. [PMID: 35205227 PMCID: PMC8872436 DOI: 10.3390/genes13020182] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2021] [Revised: 12/27/2021] [Accepted: 12/30/2021] [Indexed: 01/27/2023] Open
Abstract
Background: Whole-genome sequencing has become routine for population genetic studies. Sequencing of individuals provides maximal data but is rather expensive and fewer samples can be studied. In contrast, sequencing a pool of samples (pool-seq) can provide sufficient data, while presenting less of an economic challenge. Few studies have compared the two approaches to infer population genetic structure and diversity in real datasets. Here, we apply individual sequencing (ind-seq) and pool-seq to the study of Western honey bees (Apis mellifera). Methods: We collected honey bee workers that belonged to 14 populations, including 13 subspecies, totaling 1347 colonies, who were individually (139 individuals) and pool-sequenced (14 pools). We compared allele frequencies, genetic diversity estimates, and population structure as inferred by the two approaches. Results: Pool-seq and ind-seq revealed near identical population structure and genetic diversities, albeit at different costs. While pool-seq provides genome-wide polymorphism data at considerably lower costs, ind-seq can provide additional information, including the identification of population substructures, hybridization, or individual outliers. Conclusions: If costs are not the limiting factor, we recommend using ind-seq, as population genetic structure can be inferred similarly well, with the advantage gained from individual genetic information. Not least, it also significantly reduces the effort required for the collection of numerous samples and their further processing in the laboratory.
Collapse
Affiliation(s)
- Chao Chen
- Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing 100093, China;
- Key Laboratory of Pollinating Insect Biology, Ministry of Agriculture and Rural Affairs, Beijing 100093, China
- Correspondence: (C.C.); (M.P.)
| | - Melanie Parejo
- Department of Genetics, Physical Anthropology and Animal Physiology, Faculty of Science and Technology, University of the Basque Country (UPV/EHU), 48940 Leioa, Spain; (J.L.); (A.E.)
- Swiss Bee Research Center, Agroscope, 3003 Bern, Switzerland
- Correspondence: (C.C.); (M.P.)
| | - Jamal Momeni
- Eurofins Genomics, 8200 Aarhus, Denmark; (J.M.); (R.O.N.); (R.V.)
| | - Jorge Langa
- Department of Genetics, Physical Anthropology and Animal Physiology, Faculty of Science and Technology, University of the Basque Country (UPV/EHU), 48940 Leioa, Spain; (J.L.); (A.E.)
| | | | - Wei Shi
- Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing 100093, China;
- Key Laboratory of Pollinating Insect Biology, Ministry of Agriculture and Rural Affairs, Beijing 100093, China
| | | | - Rikke Vingborg
- Eurofins Genomics, 8200 Aarhus, Denmark; (J.M.); (R.O.N.); (R.V.)
| | - Per Kryger
- Department of Agroecology, Aarhus University, 4200 Slagelse, Denmark;
| | - Maria Bouga
- Lab of Agricultural Zoology and Entomology, Agricultural University of Athens, 11855 Athens, Greece;
| | - Andone Estonba
- Department of Genetics, Physical Anthropology and Animal Physiology, Faculty of Science and Technology, University of the Basque Country (UPV/EHU), 48940 Leioa, Spain; (J.L.); (A.E.)
| | | |
Collapse
|
21
|
Liu Y, Zong W, Diaby M, Lin Z, Wang S, Gao B, Ji T, Song C. Diversity and Evolution of pogo and Tc1/mariner Transposons in the Apoidea Genomes. BIOLOGY 2021; 10:940. [PMID: 34571816 PMCID: PMC8472432 DOI: 10.3390/biology10090940] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/10/2021] [Revised: 09/14/2021] [Accepted: 09/16/2021] [Indexed: 11/16/2022]
Abstract
Bees (Apoidea), the largest and most crucial radiation of pollinators, play a vital role in the ecosystem balance. Transposons are widely distributed in nature and are important drivers of species diversity. However, transposons are rarely reported in important pollinators such as bees. Here, we surveyed 37 bee genomesin Apoidea, annotated the pogo and Tc1/mariner transposons in the genome of each species, and performed a phylogenetic analysis and determined their overall distribution. The pogo and Tc1/mariner families showed high diversity and low abundance in the 37 species, and their proportion was significantly higher in solitary bees than in social bees. DD34D/mariner was found to be distributed in almost all species and was found in Apis mellifera, Apis mellifera carnica, Apis mellifera caucasia, and Apis mellifera mellifera, and Euglossa dilemma may still be active. Using horizontal transfer analysis, we found that DD29-30D/Tigger may have experienced horizontal transfer (HT) events. The current study displayed the evolution profiles (including diversity, activity, and abundance) of the pogo and Tc1/mariner transposons across 37 species of Apoidea. Our data revealed their contributions to the genomic variations across these species and facilitated in understanding of the genome evolution of this lineage.
Collapse
Affiliation(s)
| | | | | | | | | | | | | | - Chengyi Song
- College of Animal Science and Technology, Yangzhou University, Yangzhou 225009, China; (Y.L.); (W.Z.); (M.D.); (Z.L.); (S.W.); (B.G.); (T.J.)
| |
Collapse
|
22
|
Tanasković M, Erić P, Patenković A, Erić K, Mihajlović M, Tanasić V, Stanisavljević L, Davidović S. MtDNA Analysis Indicates Human-Induced Temporal Changes of Serbian Honey Bees Diversity. INSECTS 2021; 12:insects12090767. [PMID: 34564207 PMCID: PMC8472511 DOI: 10.3390/insects12090767] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/15/2021] [Revised: 08/21/2021] [Accepted: 08/23/2021] [Indexed: 11/25/2022]
Abstract
Simple Summary The western honey bee is one of the most economically and ecologically important species currently facing serious challenges in its whole area of distribution. The honey bee is a highly diverse species with about 30 subspecies that are adapted to regional climate factors, vegetation, pests and pathogens. The local populations of honey bees are rapidly changing and their diversity is constantly manipulated by beekeepers through the import of foreign queens, selection and migratory beekeeping. This manipulation may lead to such changes that honey bees lose their ability to thrive in the areas that were previously suitable for their wellbeing. To see how this human interference changed the genetic variability of native honey bee populations from Serbia, we sequenced part of the mitochondrial genome and compared them with published sequences. Our results suggest that human influence significantly changes the natural composition of honey bees in Serbia and that the presence of some previously reported subspecies could not be confirmed. Abstract Local populations of Apis mellifera are rapidly changing by modern beekeeping through the introduction of nonnative queens, selection and migratory beekeeping. To assess the genetic diversity of contemporary managed honey bees in Serbia, we sequenced mitochondrial tRNAleu-cox2 intergenic region of 241 worker bees from 46 apiaries at eight localities. Nine haplotypes were observed in our samples, with C2d being the most common and widespread. To evaluate genetic diversity patterns, we compared our data with 1696 sequences from the NCBI GenBank from neighbouring countries and Serbia. All 32 detected haplotypes belonged to the Southeast Europe lineage C, with two newly described haplotypes from our sample. The most frequent haplotype was C2d, followed by C2c and C1a. To distinguish A. m. carnica from A. m. macedonica, both previously reported in Serbia, PCR-RFLP analysis on the COI gene segment of mtDNA was used, and the result showed only the presence of A.m. carnica subspecies. An MDS plot constructed on pairwise FST values showed significant geographical stratification. Our samples are grouped together, but distant from the Serbian dataset from the GenBank. This, with the absence of A. m. macedonica subspecies from its historic range of distribution in southern Serbia, indicates that honey bee populations are changing rapidly due to the anthropogenic influence.
Collapse
Affiliation(s)
- Marija Tanasković
- Department of Genetics of Populations and Ecogenotoxicology, Institute for Biological Research “Siniša Stanković”—National Institute of the Republic of Serbia, University of Belgrade, Bulevar Despota Stefana 142, 11060 Belgrade, Serbia; (P.E.); (A.P.); (K.E.); (S.D.)
- Correspondence:
| | - Pavle Erić
- Department of Genetics of Populations and Ecogenotoxicology, Institute for Biological Research “Siniša Stanković”—National Institute of the Republic of Serbia, University of Belgrade, Bulevar Despota Stefana 142, 11060 Belgrade, Serbia; (P.E.); (A.P.); (K.E.); (S.D.)
| | - Aleksandra Patenković
- Department of Genetics of Populations and Ecogenotoxicology, Institute for Biological Research “Siniša Stanković”—National Institute of the Republic of Serbia, University of Belgrade, Bulevar Despota Stefana 142, 11060 Belgrade, Serbia; (P.E.); (A.P.); (K.E.); (S.D.)
| | - Katarina Erić
- Department of Genetics of Populations and Ecogenotoxicology, Institute for Biological Research “Siniša Stanković”—National Institute of the Republic of Serbia, University of Belgrade, Bulevar Despota Stefana 142, 11060 Belgrade, Serbia; (P.E.); (A.P.); (K.E.); (S.D.)
| | - Milica Mihajlović
- Center for Forensic and Applied Molecular Genetics, Faculty of Biology, University of Belgrade, Studentski trg 16, 11000 Belgrade, Serbia; (M.M.); (V.T.)
| | - Vanja Tanasić
- Center for Forensic and Applied Molecular Genetics, Faculty of Biology, University of Belgrade, Studentski trg 16, 11000 Belgrade, Serbia; (M.M.); (V.T.)
| | - Ljubiša Stanisavljević
- Center for Bee Research, Faculty of Biology, University of Belgrade, Studentski trg 16, 11000 Belgrade, Serbia;
| | - Slobodan Davidović
- Department of Genetics of Populations and Ecogenotoxicology, Institute for Biological Research “Siniša Stanković”—National Institute of the Republic of Serbia, University of Belgrade, Bulevar Despota Stefana 142, 11060 Belgrade, Serbia; (P.E.); (A.P.); (K.E.); (S.D.)
| |
Collapse
|
23
|
Henriques D, Lopes AR, Chejanovsky N, Dalmon A, Higes M, Jabal-Uriel C, Le Conte Y, Reyes-Carreño M, Soroker V, Martín-Hernández R, Pinto MA. A SNP assay for assessing diversity in immune genes in the honey bee (Apis mellifera L.). Sci Rep 2021; 11:15317. [PMID: 34321557 PMCID: PMC8319136 DOI: 10.1038/s41598-021-94833-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2021] [Accepted: 07/12/2021] [Indexed: 02/07/2023] Open
Abstract
With a growing number of parasites and pathogens experiencing large-scale range expansions, monitoring diversity in immune genes of host populations has never been so important because it can inform on the adaptive potential to resist the invaders. Population surveys of immune genes are becoming common in many organisms, yet they are missing in the honey bee (Apis mellifera L.), a key managed pollinator species that has been severely affected by biological invasions. To fill the gap, here we identified single nucleotide polymorphisms (SNPs) in a wide range of honey bee immune genes and developed a medium-density assay targeting a subset of these genes. Using a discovery panel of 123 whole-genomes, representing seven A. mellifera subspecies and three evolutionary lineages, 180 immune genes were scanned for SNPs in exons, introns (< 4 bp from exons), 3' and 5´UTR, and < 1 kb upstream of the transcription start site. After application of multiple filtering criteria and validation, the final medium-density assay combines 91 quality-proved functional SNPs marking 89 innate immune genes and these can be readily typed using the high-sample-throughput iPLEX MassARRAY system. This medium-density-SNP assay was applied to 156 samples from four countries and the admixture analysis clustered the samples according to their lineage and subspecies, suggesting that honey bee ancestry can be delineated from functional variation. In addition to allowing analysis of immunogenetic variation, this newly-developed SNP assay can be used for inferring genetic structure and admixture in the honey bee.
Collapse
Affiliation(s)
- Dora Henriques
- Centro de Investigação de Montanha, Instituto Politécnico de Bragança, Campus de Santa Apolónia, 5300-253, Bragança, Portugal
| | - Ana R Lopes
- Centro de Investigação de Montanha, Instituto Politécnico de Bragança, Campus de Santa Apolónia, 5300-253, Bragança, Portugal
| | - Nor Chejanovsky
- Agricultural Research Organization, The Volcani Center, Rishon LeTsiyon, Israel
| | - Anne Dalmon
- INRAE, Unité Abeilles et Environnement, Avignon, France
| | - Mariano Higes
- IRIAF, Instituto Regional de Investigación y Desarrollo Agroalimentario y Forestal, Laboratorio de Patología Apícola, Centro de Investigación Apícola y Agroambiental (CIAPA), Consejería de Agricultura de la Junta de Comunidades de Castilla-La Mancha, Marchamalo, Spain
| | - Clara Jabal-Uriel
- IRIAF, Instituto Regional de Investigación y Desarrollo Agroalimentario y Forestal, Laboratorio de Patología Apícola, Centro de Investigación Apícola y Agroambiental (CIAPA), Consejería de Agricultura de la Junta de Comunidades de Castilla-La Mancha, Marchamalo, Spain
| | - Yves Le Conte
- INRAE, Unité Abeilles et Environnement, Avignon, France
| | | | - Victoria Soroker
- Agricultural Research Organization, The Volcani Center, Rishon LeTsiyon, Israel
| | - Raquel Martín-Hernández
- IRIAF, Instituto Regional de Investigación y Desarrollo Agroalimentario y Forestal, Laboratorio de Patología Apícola, Centro de Investigación Apícola y Agroambiental (CIAPA), Consejería de Agricultura de la Junta de Comunidades de Castilla-La Mancha, Marchamalo, Spain
- Instituto de Recursos Humanos para la Ciencia y la Tecnología (INCRECYT-FEDER), Fundación Parque Científico y Tecnológico de Castilla-La Mancha, 02006, Albacete, Spain
| | - M Alice Pinto
- Centro de Investigação de Montanha, Instituto Politécnico de Bragança, Campus de Santa Apolónia, 5300-253, Bragança, Portugal.
| |
Collapse
|
24
|
Tofilski A, Căuia E, Siceanu A, Vișan GO, Căuia D. Historical Changes in Honey Bee Wing Venation in Romania. INSECTS 2021; 12:insects12060542. [PMID: 34200932 PMCID: PMC8230453 DOI: 10.3390/insects12060542] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/29/2021] [Revised: 05/31/2021] [Accepted: 06/08/2021] [Indexed: 01/30/2023]
Abstract
Simple Summary Honey bees, in addition to producing honey, are important pollinators of wild and cultivated plants. Unfortunately, in some places, the population of honey bees is declining. One of the factors that affect their survival is adaptation to the local environment. Bees native to a particular area survive better than those imported. Despite this fact, some beekeepers import non-native bees and use them in their apiaries. Imported bees produce hybrids with bees from surrounding colonies because beekeepers do not control their mating. In consequence, the whole population can change. In this study, we verified how the population of Romanian bees has changed over the last four decades. We found significant temporal changes in wing venation. Despite these changes, the two major subpopulations of bees separated by mountains remain distinct. We provide a tool for the easy identification of native bees from Romania, which can help to protect them. Abstract The honey bee (Apis mellifera) is an ecologically and economically important species that provides pollination services to natural and agricultural systems. The biodiversity of the honey bee is being endangered by the mass import of non-native queens. In many locations, it is not clear how the local populations have been affected by hybridisation between native and non-native bees. There is especially little information about temporal changes in hybridisation. In Romania, A. m. carpatica naturally occurs, and earlier studies show that there are two subpopulations separated by the Carpathian Mountains. In this study, we investigated how the arrangement of veins in bees’ wings (venation) has changed in Romanian honey bees in the last four decades. We found that in the contemporary population of Romanian bees, there are still clear differences between the intra- and extra-Carpathian subpopulations, which indicates that natural variation among honey bees is still being preserved. We also found significant differences between bees collected before and after 2000. The observed temporal changes in wing venation are most likely caused by hybridisation between native bees and non-native bees sporadically introduced by beekeepers. In order to facilitate conservation and the monitoring of native Romanian bees, we developed a method facilitating their identification.
Collapse
Affiliation(s)
- Adam Tofilski
- Department of Zoology and Animal Welfare, University of Agriculture in Krakow, Al. 29 Listopada 54, 31-425 Krakow, Poland
- Correspondence:
| | - Eliza Căuia
- Honeybee Genetics and Breeding Laboratory, Institute for Beekeeping Research and Development, Blv Ficusului, No. 42, Sector 1, 013975 Bucharest, Romania; (E.C.); (A.S.); (G.O.V.); (D.C.)
| | - Adrian Siceanu
- Honeybee Genetics and Breeding Laboratory, Institute for Beekeeping Research and Development, Blv Ficusului, No. 42, Sector 1, 013975 Bucharest, Romania; (E.C.); (A.S.); (G.O.V.); (D.C.)
| | - Gabriela Oana Vișan
- Honeybee Genetics and Breeding Laboratory, Institute for Beekeeping Research and Development, Blv Ficusului, No. 42, Sector 1, 013975 Bucharest, Romania; (E.C.); (A.S.); (G.O.V.); (D.C.)
| | - Dumitru Căuia
- Honeybee Genetics and Breeding Laboratory, Institute for Beekeeping Research and Development, Blv Ficusului, No. 42, Sector 1, 013975 Bucharest, Romania; (E.C.); (A.S.); (G.O.V.); (D.C.)
| |
Collapse
|
25
|
Gray level co-occurrence matrix and extreme learning machine for Covid-19 diagnosis. INTERNATIONAL JOURNAL OF COGNITIVE COMPUTING IN ENGINEERING 2021; 2:93-103. [PMCID: PMC8177375 DOI: 10.1016/j.ijcce.2021.05.001] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/24/2021] [Revised: 05/24/2021] [Accepted: 05/30/2021] [Indexed: 06/01/2023]
Abstract
Background Chest CT is considered to be a more accurate method for diagnosing suspected patients. However, with the spread of the epidemic, traditional diagnostic methods have been unable to meet the requirements of efficiency and speed. Therefore, it is necessary to use artificial intelligence to help people make efficient and accurate judgments. A number of studies have shown that it is feasible to use deep learning methods to help people diagnose COVID-19. However, most of the existing methods are single-layer neural network structures, and their accuracy and efficiency need to be improved. Method In this scheme, a hybrid model is adopted. Firstly, the gray co-occurrence matrix is used to extract the features of the images, and then the extreme learning machine is used for classification. Results The experimental results show that the model proposed in this paper is feasible and can help medical staff to accurately determine suspected patients for subsequent isolation and treatment.
Collapse
|
26
|
More S, Bampidis V, Benford D, Bragard C, Halldorsson T, Hernández‐Jerez A, Bennekou SH, Koutsoumanis K, Machera K, Naegeli H, Nielsen SS, Schlatter J, Schrenk D, Silano V, Turck D, Younes M, Arnold G, Dorne J, Maggiore A, Pagani S, Szentes C, Terry S, Tosi S, Vrbos D, Zamariola G, Rortais A. A systems-based approach to the environmental risk assessment of multiple stressors in honey bees. EFSA J 2021; 19:e06607. [PMID: 34025804 PMCID: PMC8135085 DOI: 10.2903/j.efsa.2021.6607] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
The European Parliament requested EFSA to develop a holistic risk assessment of multiple stressors in honey bees. To this end, a systems-based approach that is composed of two core components: a monitoring system and a modelling system are put forward with honey bees taken as a showcase. Key developments in the current scientific opinion (including systematic data collection from sentinel beehives and an agent-based simulation) have the potential to substantially contribute to future development of environmental risk assessments of multiple stressors at larger spatial and temporal scales. For the monitoring, sentinel hives would be placed across representative climatic zones and landscapes in the EU and connected to a platform for data storage and analysis. Data on bee health status, chemical residues and the immediate or broader landscape around the hives would be collected in a harmonised and standardised manner, and would be used to inform stakeholders, and the modelling system, ApisRAM, which simulates as accurately as possible a honey bee colony. ApisRAM would be calibrated and continuously updated with incoming monitoring data and emerging scientific knowledge from research. It will be a supportive tool for beekeeping, farming, research, risk assessment and risk management, and it will benefit the wider society. A societal outlook on the proposed approach is included and this was conducted with targeted social science research with 64 beekeepers from eight EU Member States and with members of the EU Bee Partnership. Gaps and opportunities are identified to further implement the approach. Conclusions and recommendations are made on a way forward, both for the application of the approach and its use in a broader context.
Collapse
|