It is assumed that wild honey bees have become largely extinct across Europe since the 1980s, following the introduction of exotic ectoparasitic mite (Varroa) and the associated spillover of various pathogens. However, several recent studies reported on unmanaged colonies that survived the Varroa mite infestation. Herewith, we present another case of unmanaged, free-living population of honey bees in SE Europe, a rare case of feral bees inhabiting a large and highly populated urban area: Belgrade, the capital of Serbia. We compiled a massive data-set derived from opportunistic citizen science (>1300 records) during the 2011–2017 period and investigated whether these honey bee colonies and the high incidence of swarms could be a result of a stable, self-sustaining feral population (i.e., not of regular inflow of swarms escaping from local managed apiaries), and discussed various explanations for its existence. We also present the possibilities and challenges associated with the detection and effective monitoring of feral/wild honey bees in urban settings, and the role of citizen science in such endeavors. Our results will underpin ongoing initiatives to better understand and support naturally selected resistance mechanisms against the Varroa mite, which should contribute to alleviating current threats and risks to global apiculture and food production security.

Large-scale genomic resources can place genetic variation into an ecologically informed context. To advance our understanding of the population genetics of the fruit fly Drosophila melanogaster, we present an expanded release of the community-generated population genomics resource Drosophila Evolution over Space and Time (DEST 2.0; https://dest.bio/). This release includes 530 high-quality pooled libraries from flies collected across six continents over more than a decade (2009 to 2021), most at multiple time points per year; 211 of these libraries are sequenced and shared here for the first time. We used this enhanced resource to elucidate several aspects of the species' demographic history and identify novel signs of adaptation across spatial and temporal dimensions. For example, we showed that the spatial genetic structure of populations is stable over time, but that drift due to seasonal contractions of population size causes populations to diverge over time. We identified signals of adaptation that vary between continents in genomic regions associated with xenobiotic resistance, consistent with independent adaptation to common pesticides. Moreover, by analyzing samples collected during spring and fall across Europe, we provide new evidence for seasonal adaptation related to loci associated with pathogen response. Furthermore, we have also released an updated version of the DEST genome browser. This is a useful tool for studying spatiotemporal patterns of genetic variation in this classic model system.

A number of recent studies have shown that the pattern of mitochondrial DNA variation and evolution is at odds with a neutral equilibrium model. Theory has suggested that selection on mitonuclear genotypes can act to maintain stable mitonuclear polymorphism within populations. However, this effect largely relies upon selection being either sex-specific or frequency dependent. Here, we use mitonuclear introgression lines to assess differences in a series of key life-history traits (egg-to-adult developmental time, viability, offspring sex-ratio, adult longevity and resistance to desiccation) in Drosophila subobscura fruit flies carrying one of three different sympatric mtDNA haplotypes. We found functional differences between these sympatric mtDNA haplotypes, but these effects were contingent upon the nuclear genome with which they were co-expressed. Further, we demonstrate a significant mitonuclear genetic effect on adult sex ratio, as well as a sex × mtDNA × nuDNA interaction for adult longevity. The observed effects suggest that sex specific mitonuclear selection contributes to the maintenance of mtDNA polymorphism and to mitonuclear linkage disequilibrium in this model system.

Once a widespread and common species across the region of southeast Europe, the Griffon vulture is now confined to small and isolated populations across the Balkan Peninsula. The population from Serbia with 290 couples represents its biggest and most viable population that can serve as an important reservoir of genetic diversity from which the birds can be used for the region’s reintroduction or recolonization programs. To estimate the level of genetic diversity, the mitochondrial Cytb gene from 58 unrelated birds sampled during the marking in the nests was sequenced and compared to the homologous Griffon vulture sequences available in publicly accessible online databases. Phylogeographic analysis based on Cytb sequences showed that the most frequent haplotype is found in all Griffon vulture populations and that each population possesses private haplotypes. Our data suggest that the Griffon vulture population from Serbia should be used as a source population for restocking and reintroduction programs in the region. The observed genetic differentiation between the populations from the Iberian and Balkan Peninsulas suggest that the introduction of foreign birds from remote populations should be avoided and that birds from indigenous or neighboring populations, if available, should be used instead.

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 tRNA^leu-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 F_ST 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.

Socioeconomic interests and beekeeper preferences have often taken precedence over the conservation of locally native honey bee subspecies, leading to the predominance of admixture populations in human-dominated areas. To assess the genetic diversity of contemporary managed Serbian honey bee colonies, we used 14 microsatellite loci and analyzed 237 worker bees from 46 apiaries in eight localities of northern and southern Serbia. Furthermore, we compared data for nine microsatellite loci with 338 individuals from Italy, Hungary, Poland, and Spain. The standard parameters of genetic diversity in Serbian honey bee populations were in line with other analyses, although somewhat smaller. STRUCTURE analysis showed the existence of two equally distributed genetic clusters and Analysis of molecular variances could not confirm the presence of a geographically discrete population but showed local differences. Discriminant analysis of principal components showed overlapping of worker bees from different parts of Serbia. Clear genetic differentiation can be observed when comparing all populations between geographical regions and their corresponding subspecies. The absence of the A. m. macedonica subspecies from its historical distribution range in southern Serbia as well as the lack of distinctive geographical groups suggest that selective breeding, queen import, and migratory beekeeping practices strongly influence the genetic structure and diversity of honey bees, leading to the genetic uniformization and creation of the admixture population.