Population History Shapes Responses to Different Temperature Regimes in Drosophila subobscura.

Drosophila subobscura is considered a good model species for investigation of a population's ability to adapt and cope with climate changes. Decade long research has shown that inversion frequencies change in response to environmental factors indicating their role in adaptation to novel environments. The mechanisms behind organisms' responses to temperature are complex, involving changes in physiology, behavior, gene expression and regulation. On the other hand, a population's ability to respond to suboptimal conditions depends on standing genetic variation and population history. In order to elucidate the role of local adaptation in population response to the changing temperature, we investigated the response to temperature in D. subobscura individuals originating from two different altitudes by combining traditional cytogenetic techniques with assessing the levels of Hsp70 protein expression. Inversion polymorphism was assessed in the flies sampled from natural populations and in flies reared in laboratory conditions at three different temperatures after five and sixteen generations and Hsp70 protein expression profile in 12th generation flies at the basal level and after heat shock induction. Our results indicate that local adaptation and population history influence population response to the changing temperature.

Sex ratio and relatedness in the Griffon vulture (Gyps fulvus) population of Serbia.

Background: Once a widespread 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 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 programmes. The available genetic data for this valuable population are scarce and as a protected species that belongs to the highly endangered vulture group, it needs to be well described so that it can be properly managed and used as a restocking population. Considering the serious recent bottleneck event that the Griffon vulture population from Serbia experienced we estimated the overall relatedness among the birds from this population. Sex ratio, another important parameter that shows the vitality and strength of the population was evaluated as well. Methods: During the annual monitoring that was performed in the period from 2013-2021, we collected blood samples from individual birds that were marked in the nests. In total, 169 samples were collected and each was used for molecular sexing while 58 presumably unrelated birds from different nests were used for inbreeding and relatedness analyses. The relatedness was estimated using both biparentally (10 microsatellite loci) and uniparentally (Cytb and D-loop I of mitochondrial DNA) inherited markers. Results: The level of inbreeding was relatively high and on average it was 8.3% while the mean number of relatives for each bird was close to three. The sex ratio was close to 1:1 and for the analysed period of 9 years, it didn't demonstrate a statistically significant deviation from the expected ratio of 1:1, suggesting that this is a stable and healthy population. Our data suggest that, even though a relatively high level of inbreeding can be detected among the individual birds, the Griffon vulture population from Serbia can be used as a source population for restocking and reintroduction programmes in the region. These data combined with previously observed genetic differentiation between the populations from the Iberian and Balkan Peninsulas suggest that the introduction of foreign birds should be avoided and that local birds should be used instead.

Urban ecosystem drives genetic diversity in feral honey bee.

Urbanization can change biodiversity in both directions, positive and negative, and despite the rising global trend of urban beekeeping, little is known about the impact of urbanization on the genetic diversity of honey bees. We investigate how urbanization affects the genetic variability of feral and managed honey bee colonies that are spread throughout the entire city, even in highly urban areas, through genetic analysis of 82 worker bees. We found convincing evidence of high genetic differentiation between these two groups. Additionally, by comparing city samples with 241 samples from 46 apiaries in rural parts of the country, variations in mitochondrial tRNAleu-cox2 intergenic region and microsatellite loci indicated that feral colonies have distinct patterns of genetic diversity. These results, with evidence that feral honey bees find niches within highly modified and human-dominated urban landscapes, lead us to conclude that urbanization is a driver of the genetic diversity of feral honey bees in the city.

Stress Resistance Traits under Different Thermal Conditions in Drosophila subobscura from Two Altitudes.

Global warming and climate change are affecting many insect species in numerous ways. These species can develop diverse mechanisms as a response to variable environmental conditions. The rise in mean and extreme temperatures due to global warming and the importance of the population's ability to adapt to temperature stress will further increase. In this study, we investigated thermal stress response, which is considered to be one of the crucial elements of population fitness and survival in fast-changing environments. The dynamics and variation of thermal stress resistance traits in D. subobscura flies originating from two natural populations sampled from different altitudes were analysed. Three different temperature regimes (25 °C, 19 °C, and 16 °C) were used for the F1 progeny from both localities to establish six experimental groups and investigate stress resistance traits: desiccation resistance, heat knock-down resistance, starvation resistance, and chill-coma recovery time. We detected that laboratory thermal conditions and population origin may have an effect on the analysed traits, and that sex also significantly influences stress resistance. Individuals from the lower altitude reared at higher temperatures show inferior resistance to thermal shock.

Temperature-Specific and Sex-Specific Fitness Effects of Sympatric Mitochondrial and Mito-Nuclear Variation in Drosophila obscura.

The adaptive significance of sympatric mitochondrial (mtDNA) variation and the role of selective mechanisms that maintain it are debated to this day. Isofemale lines of Drosophila obscura collected from four populations were backcrossed within populations to construct experimental lines, with all combinations of mtDNA Cyt b haplotypes and nuclear genetic backgrounds (nuDNA). Individuals of both sexes from these lines were then subjected to four fitness assays (desiccation resistance, developmental time, egg-to-adult viability and sex ratio) on two experimental temperatures to examine the role of temperature fluctuations and sex-specific selection, as well as the part that interactions between the two genomes play in shaping mtDNA variation. The results varied across populations and fitness components. In the majority of comparisons, they show that sympatric mitochondrial variants affect fitness. However, their effect should be examined in light of interactions with nuDNA, as mito-nuclear genotype was even more influential on fitness across all components. We found both sex-specific and temperature-specific differences in mitochondrial and mito-nuclear genotype ranks in all fitness components. The effect of temperature-specific selection was found to be more prominent, especially in desiccation resistance. From the results of different components tested, we can also infer that temperature-specific mito-nuclear interactions rather than sex-specific selection on mito-nuclear genotypes have a more substantial role in preserving mtDNA variation in this model species.

Further Evidence of Population Admixture in the Serbian Honey Bee Population.

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.

Genetic Diversity Analysis of Mitochondrial Cytb Gene, Phylogeny and Phylogeography of Protected Griffon Vulture (Gyps fulvus) from Serbia.

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.

Life History Traits in Two Drosophila Species Differently Affected by Microbiota Diversity under Lead Exposure.

Life history traits determine the persistence and reproduction of each species. Factors that can affect life history traits are numerous and can be of different origin. We investigated the influence of population origin and heavy metal exposure on microbiota diversity and two life history traits, egg-to-adult viability and developmental time, in Drosophila melanogaster and Drosophila subobscura, grown in the laboratory on a lead (II) acetate-saturated substrate. We used 24 samples, 8 larval and 16 adult samples (two species × two substrates × two populations × two sexes). The composition of microbiota was determined by sequencing (NGS) of the V3-V4 variable regions of the 16S rRNA gene. The population origin showed a significant influence on life history traits, though each trait in the two species was affected differentially. Reduced viability in D. melanogaster could be a cost of fast development, decrease in Lactobacillus abundance and the presence of Wolbachia. The heavy metal exposure in D. subobscura caused shifts in developmental time but maintained the egg-to-adult viability at a similar level. Microbiota diversity indicated that the Komagataeibacter could be a valuable member of D. subobscura microbiota in overcoming the environmental stress. Research on the impact of microbiota on the adaptive response to heavy metals and consequently the potential tradeoffs among different life history traits is of great importance in evolutionary research.

Unprecedented Density and Persistence of Feral Honey Bees in Urban Environments of a Large SE-European City (Belgrade, Serbia).

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.

MtDNA Analysis Indicates Human-Induced Temporal Changes of Serbian Honey Bees Diversity.

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.

Drosophila Evolution over Space and Time (DEST): A New Population Genomics Resource.

Drosophila melanogaster is a leading model in population genetics and genomics, and a growing number of whole-genome data sets from natural populations of this species have been published over the last years. A major challenge is the integration of disparate data sets, often generated using different sequencing technologies and bioinformatic pipelines, which hampers our ability to address questions about the evolution of this species. Here we address these issues by developing a bioinformatics pipeline that maps pooled sequencing (Pool-Seq) reads from D. melanogaster to a hologenome consisting of fly and symbiont genomes and estimates allele frequencies using either a heuristic (PoolSNP) or a probabilistic variant caller (SNAPE-pooled). We use this pipeline to generate the largest data repository of genomic data available for D. melanogaster to date, encompassing 271 previously published and unpublished population samples from over 100 locations in >20 countries on four continents. Several of these locations have been sampled at different seasons across multiple years. This data set, which we call Drosophila Evolution over Space and Time (DEST), is coupled with sampling and environmental metadata. A web-based genome browser and web portal provide easy access to the SNP data set. We further provide guidelines on how to use Pool-Seq data for model-based demographic inference. Our aim is to provide this scalable platform as a community resource which can be easily extended via future efforts for an even more extensive cosmopolitan data set. Our resource will enable population geneticists to analyze spatiotemporal genetic patterns and evolutionary dynamics of D. melanogaster populations in unprecedented detail.

The discovery, distribution, and diversity of DNA viruses associated with Drosophila melanogaster in Europe.

Drosophila melanogaster is an important model for antiviral immunity in arthropods, but very few DNA viruses have been described from the family Drosophilidae. This deficiency limits our opportunity to use natural host-pathogen combinations in experimental studies, and may bias our understanding of the Drosophila virome. Here, we report fourteen DNA viruses detected in a metagenomic analysis of 6668 pool-sequenced Drosophila, sampled from forty-seven European locations between 2014 and 2016. These include three new nudiviruses, a new and divergent entomopoxvirus, a virus related to Leptopilina boulardi filamentous virus, and a virus related to Musca domestica salivary gland hypertrophy virus. We also find an endogenous genomic copy of galbut virus, a double-stranded RNA partitivirus, segregating at very low frequency. Remarkably, we find that Drosophila Vesanto virus, a small DNA virus previously described as a bidnavirus, may be composed of up to twelve segments and thus represent a new lineage of segmented DNA viruses. Two of the DNA viruses, Drosophila Kallithea nudivirus and Drosophila Vesanto virus are relatively common, found in 2 per cent or more of wild flies. The others are rare, with many likely to be represented by a single infected fly. We find that virus prevalence in Europe reflects the prevalence seen in publicly available datasets, with Drosophila Kallithea nudivirus and Drosophila Vesanto virus the only ones commonly detectable in public data from wild-caught flies and large population cages, and the other viruses being rare or absent. These analyses suggest that DNA viruses are at lower prevalence than RNA viruses in D.melanogaster, and may be less likely to persist in laboratory cultures. Our findings go some way to redressing an earlier bias toward RNA virus studies in Drosophila, and lay the foundation needed to harness the power of Drosophila as a model system for the study of DNA viruses.

Highly contiguous assemblies of 101 drosophilid genomes.

Over 100 years of studies in Drosophila melanogaster and related species in the genus Drosophila have facilitated key discoveries in genetics, genomics, and evolution. While high-quality genome assemblies exist for several species in this group, they only encompass a small fraction of the genus. Recent advances in long-read sequencing allow high-quality genome assemblies for tens or even hundreds of species to be efficiently generated. Here, we utilize Oxford Nanopore sequencing to build an open community resource of genome assemblies for 101 lines of 93 drosophilid species encompassing 14 species groups and 35 sub-groups. The genomes are highly contiguous and complete, with an average contig N50 of 10.5 Mb and greater than 97% BUSCO completeness in 97/101 assemblies. We show that Nanopore-based assemblies are highly accurate in coding regions, particularly with respect to coding insertions and deletions. These assemblies, along with a detailed laboratory protocol and assembly pipelines, are released as a public resource and will serve as a starting point for addressing broad questions of genetics, ecology, and evolution at the scale of hundreds of species.

Effect of lead pollution on fitness and its dependence on heterozygosity in Drosophila subobscura.

Lead is one of the most present contaminants in the environment, and different species respond differently to this type of polution. If combined with genomic stress, lead may act synergistically, causing significant decrease of fitness components. We used two genetically diverse Drosophila subobscura populations (regarding both putatively adaptive inversion and microsatellite loci polymorphisms) originating from two ecologically distinct habitats. To establish different levels of genome heterozygosity, series of intraline, intrapopulation and interpopulation crosses were made. The progeny were reared on a standard medium and a medium with 200 µg/mL of lead acetate. Development time was significantly extended to all groups reared on lead. The progeny of intraline crosses showed significantly extended development time compared to all other groups. The obtained results suggest that genome heterozygosity reduces the effect of lead pollution.

Sex-specific effects of sympatric mitonuclear variation on fitness in Drosophila subobscura.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.