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.

Signatures of Positive Selection in the Genome of Apis Mellifera Carnica: A Subspecies of European Honeybees.

The Progress on Honeybee Genome Research

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

Pan‐genome analysis highlights the role of structural variation in the evolution and environmental adaptation of Asian honeybees

Geometric morphology and population genomics provide insights into the adaptive evolution of Apis cerana in Changbai Mountain

Positive selection in cytochrome P450 genes is associated with gonad phenotype and mating strategy in social bees

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

Population genomics and morphological features underlying the adaptive evolution of the eastern honey bee (Apis cerana)

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

A key gene for the climatic adaptation of Apis cerana populations in China according to selective sweep analysis

Genomic and transcriptomic analysis of the Asian honeybee Apis cerana provides novel insights into honeybee biology

Uncovering the novel characteristics of Asian honey bee, Apis cerana, by whole genome sequencing

Genome Sequencing and Comparative Genomics of Honey Bee Microsporidia, Nosema Apis Reveal Novel Insights into Host-Parasite Interactions

Dynamic Evolution of Repetitive Elements and Chromatin States in Apis mellifera Subspecies

Complete Mitochondrial Genome of Apis Cerana (hymenoptera: Apidae) from Two Geographical Regions: Insights into Structure and Genetic Differentiation

Genomic and morphological analyses of Asian honeybee subspecies, Apis cerana sanshasis, suggest it's evolutional origin and adaption to island

Genomic Analyses of Asian Honeybee from the Sansha Island in the South China Sea, Suggest It’s Evolutional Origin and Environmental Adaption

Transcriptomic analysis to uncover genes affecting cold resistance in the Chinese honey bee (Apis cerana cerana)

Defining Honeybee Subspecies in an Evolutionary Context Warrants Strategized Conservation

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

De Novo Genome Assembly of Chinese Plateau Honeybee Unravels Intraspecies Genetic Diversity in the Eastern Honeybee, Apis cerana