Search

Abstract Background Research into the genetic diversity of honey bee (Apis mellifera L.) populations has become increasingly significant in recent decades, primarily due to population declines attributed to human activities and climate change. As a species of great importance, breeding programs that leverage understanding of genomic diversity could offer solutions to mitigate these challenges. The objective of this study was to examine the genomic diversity and population structure of Carniolan honey bees (Apis mellifera carnica) using the Illumina SNP chip on a large honey bee sample collected from Central and South-Eastern European countries. The study also aims to offer recommendations for future breeding programs. Results Our analysis involved Discriminant Analysis of Principal Components (DAPC), heterozygosity, admixture analysis, fixation indices (FST), Neighbour-Joining tree, gene flow and Isolation-by-distance analysis. DAPC indicated distinct separation between the Carniolan and Italian honey bee (Apis mellifera ligustica) populations, whereas the admixture analysis revealed varying levels of gene flow and genetic admixture within the Carniolan honey bee populations, demonstrating closer relationships between specific geographic regions (confirmed by Isolation-by-distance analysis). Furthermore, the research of heterozygosity, genomic inbreeding, pairwise FST values, and Neighbour-Joining tree provided insights into the patterns of genetic differentiation and similarity among the populations of Carniolan honey bee within its natural habitat. We have observed genetic homogeneity of the Carniolan honey bee population when considered in a broader genetic/geographical context. However, the Carniolan honey bee has sufficient genetic diversity in its geographical home range that needs to be carefully monitored and maintained. Conclusions This study provides important insights into the genetic composition, differentiation, and relationships among Carniolan honey bee populations in Central and South-Eastern European countries. The findings are crucial for conservation efforts, breeding programs, and sustainable beekeeping practices. They emphasise the importance of considering genetic factors and population structure in the breeding and management of honey bees. By understanding these genetic relationships, we can develop strategies to preserve genetic diversity, improve breeding outcomes, and ensure the resilience of honey bee populations in the face of environmental changes and challenges. This knowledge can also inform policy makers and stakeholders on best practices to maintain healthy bee populations, which are vital for ecosystem services and agricultural productivity.

A recent comprehensive genomic analysis based on 50K SNP profiles has shown that the regional Balkan sheep populations have considerable genetic overlap but are distinctly different from surrounding breeds. All eight Croatian sheep breeds were represented by a small number of individuals per breed. Here, we genotyped 220 individuals representing the native Croatian sheep breeds (Istrian Sheep, Krk Island Sheep, Cres Island Sheep, Rab Island Sheep, Lika Pramenka, Pag Island Sheep, Dalmatian Pramenka, Dubrovnik Sheep) and mouflon using the Ovine Infinium® HD SNP BeadChip (606,006 SNPs). In addition, we included publicly available Balkan Pramenka and other Mediterranean sheep breeds. Our analyses revealed the complex population structure of Croatian sheep breeds and their origin and geographic barriers (island versus mainland). Migration patterns confirmed the historical establishment of breeds and the pathways of gene flow. Inbreeding coefficients (FROH>2 Mb) between sheep populations ranged from 0.025 to 0.070, with lower inbreeding coefficients observed in Dalmatian Pramenka and Pag Island Sheep and higher inbreeding in Dubrovnik sheep. The estimated effective population size ranged from 61 to 1039 for Krk Island Sheep and Dalmatian Pramenka, respectively. Higher inbreeding levels and lower effective population size indicate the need for improved conservation management to maintain genetic diversity in some breeds. Our results will contribute to breeding and conservation strategies of native Croatian sheep breeds.

The majority of the nearly 400 existing local pig breeds are adapted to specific environments and human needs. The demand for large production quantities and the industrialized pig production have caused a rapid decline of many local pig breeds in recent decades. Black Slavonian pig and Turopolje pig, the latter highly threatened, are the two Croatian local indigenous breeds typically grown in extensive or semi-intensive systems. In order to guide a long-term breeding program to prevent the disappearance of these breeds, we analyzed their genetic diversity, inbreeding level and relationship with other local breeds across the world, as well as modern breeds and several wild populations, using high throughput genomic data obtained using the Illumina Infinium PorcineSNP60 v2 BeadChip. Multidimensional scaling analysis positioned Black Slavonian pigs close to the UK/North American breeds, while the Turopolje pig clustered within the Mediterranean breeds. Turopolje pig showed a very high inbreeding level (FROH>4Mb = 0.400 and FROH>8Mb = 0.332) that considerably exceeded the level of full-sib mating, while Black Slavonian pig showed much lower inbreeding (FROH>4Mb = 0.098 and FROH>8Mb = 0.074), indicating a planned mating strategy. In Croatian local breeds we identified several genome regions showing adaptive selection signals that were not present in commercial breeds. The results obtained in this study reflect the current genetic status and breeding management of the two Croatian indigenous local breeds. Given the small populations of both breeds, a controlled management activity has been implemented in Black Slavonian pigs since their commercial value has been recognized. In contrast, the extremely high inbreeding level observed in Turopolje pig argues for an urgent conservation plan with a long-term, diversity-oriented breeding program.