Abstract:Background: The honey bee reference genome, HAv3.1, was produced from a commercial line sample that was thought to have a largely dominant Apis mellifera ligustica genetic background. Apis mellifera mellifera, often referred to as the black bee, has a separate evolutionary history and is the original type in western and northern Europe. Growing interest in this subspecies for conservation and non-professional apicultural practices, together with the necessity of deciphering genome backgrounds in hybrids, triggered the necessity for a specific genome assembly. Moreover, having several high-quality genomes is becoming key for taking structural variations into account in pangenome analyses. Results: Pacific Bioscience technology long reads were produced from a single haploid black bee drone. Scaffolding contigs into chromosomes was done using a high-density genetic map. This allowed for re-estimation of the recombination rate, which was over-estimated in some previous studies due to mis-assemblies, which resulted in spurious inversions in the older reference genomes. The sequence continuity obtained was very high and the only limit towards continuous chromosome-wide sequences seemed to be due to tandem repeat arrays that were usually longer than 10 kb and that belonged to two main families, the 371 and 91 bp repeats, causing problems in the assembly process due to high internal sequence similarity. Our assembly was used together with the reference genome to genotype two structural variants by a pangenome graph approach with Graphtyper2. Genotypes obtained were either correct or missing, when compared to an approach based on sequencing depth analysis, and genotyping rates were 89 and 76% for the two variants. Conclusions: Our new assembly for the Apis mellifera mellifera honey bee subspecies demonstrates the utility of multiple high-quality genomes for the genotyping of structural variants, with a test case on two insertions and deletions. It will therefore be an invaluable resource for future studies, for instance by including structural variants in GWAS. Having used a single haploid drone for sequencing allowed a refined analysis of very large tandem repeat arrays, raising the question of their function in the genome. High quality genome assemblies for multiple subspecies such as presented here, are crucial for emerging projects using pangenomes.

Inferring long-term and short-term determinants of genetic diversity in honey bees: Beekeeping impact and conservation strategies

The Influencing Factors and Primary Research Methods on Genetic Diversity of the Honeybee

Natural clines and human management impact the genetic structure of Algerian honey bee populations

Impacts of beekeeping on wild bee diversity and pollination networks in the Aegean Archipelago

Estimating genus-specific effects of non-native honey bees and urbanization on wild bee communities: A case study in Maryland, United States

A Melissopalynologic dynamics appraisal revealed the strategic foraging adaptation of the honeybees (Apis mellifera Linnaeus 1758) to the anthropogenic impacts in the Republic of Benin

Sequence-based genome-wide association studies reveal the polygenic architecture of resistance in Western honey bees

Genetic divergence and functional convergence of gut bacteria between the Eastern honey bee Apis cerana and the Western honey bee Apis mellifera

A century of wild bee sampling: historical data and neural network analysis reveal ecological traits associated with species loss

Defining Honeybee Subspecies in an Evolutionary Context Warrants Strategized Conservation

Genus-Wide Characterization of Bumblebee Genomes Provides Insights into Their Evolution and Variation in Ecological and Behavioral Traits

Whole-genome resequencing reveals genetic diversity and adaptive evolution in Chinese honeybee (Apis cerana cerana) in Guizhou, China

Population Structure, Demographic History, and Adaptation of Giant Honeybees in China Revealed by Population Genomic Data

Population Structure and Genetic Diversity of Chinese Honeybee (Apis Cerana Cerana) in Central China

Improving wild bee monitoring, sampling methods, and conservation

The black honey bee genome: insights on specific structural elements and a first step towards pangenomes

Genus-wide Characterization of Bumblebee Genomes Reveals Variation Associated with Key Ecological and Behavioral Traits of Pollinators

Nuclear and Mitochondrial DNA Suggest That Nature Reserve Maintains Novel Haplotypes and Genetic Diversity of Honeybees (Apis cerana)

Field Populations of Wild Apis cerana Honey Bees Exhibit Increased Genetic Diversity Under Pesticide Stress Along an Agricultural Intensification Gradient in Eastern India

Distinct genome architecture underlies fine-scale population differentiation in two common European bumblebees (Bombus pascuorum and Bombus lapidarius)

Crop pollination from native bees at risk from agricultural intensification