BackgroundA 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.ResultsWe 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.ConclusionsThe observed effects suggest that sex specific mitonuclear selection contributes to the maintenance of mtDNA polymorphism and to mitonuclear linkage disequilibrium in this model system.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-015-0421-2) contains supplementary material, which is available to authorized users.

The genetic structure of Drosophila subobscura from the Balkan Peninsula was studied with respect to restriction site polymorphism of mitochondrial DNA in populations from the Derventa River Gorge and Sicevo Gorge (Serbia). To investigate the role of cytonuclear interactions in shaping mitochondrial DNA variability in natural populations of this species, the study was complemented with the analysis of linkage disequilibria between mitochondrial haplotypes and chromosomal inversion arrangements. Similar to other populations of D. subobscura, two main haplotypes (I and II) were found, as well as a series of less common ones. The frequencies of haplotypes I and II accounted for 25.8% and 71.0%, respectively, in the population from the Derventa River Gorge, and for 32.4% and 58.1%, respectively, in the population from Sicevo Gorge. One of the haplotypes harbored a large insertion (2.7 kb) in the A+T rich region. The frequency distribution of both haplotypes did not depart from neutrality. Contrary to prior studies, we did not detect any significant linkage disequilibrium between the two most frequent mtDNA haplotypes and any of the chromosomal arrangements in either of the populations. We conclude that linkage disequilibrium is not a general occurrence in natural populations of D. subobscura, and we discuss how transient coadaptations, ecologically specific selective pressures, and demographics could contribute to population-specific patterns of linkage disequilibrium.

Differences in heavy metal tolerance among separate populations of the same species have often been interpreted as local adaptation. Persistence of differences after removing the stressor indicates that mechanisms responsible for the increased tolerance were genetically determined. Drosophila subobscura Collin (Diptera: Drosophilidae) populations were sampled from two localities with different history of heavy metal pollution, and reared for eight generations in the laboratory on a standard medium and on media with different concentrations of lead (Pb). To determine whether flies from different natural populations exposed to the Pb-contaminated media in the laboratory show population specific variability in fitness components over generations, experimental groups with different concentrations of lead were assayed in three generations (F2 , F5 , and F8 ) for fecundity, developmental time, and egg-to-adult viability. On the contaminated medium, fecundity was reduced in later generations and viability was increased, irrespective of the environmental origin of populations. For both populations, developmental time showed a tendency of slowing down on media with lead. Faster development was observed in later generations. Preadaptation to contamination, meaning higher fecundity, higher viability, and faster egg to adult development in all studied generations, was found in D. subobscura originating from the locality with a higher level of heavy metal pollution.

We analyzed changes in the genetic structure and effective population size of two ecologically distinct populations of Drosophila subobscura over several years. Population sizes of D. subobscura in beech and oak wood habitats for a period of 6 years were estimated by the capture-mark-release-recapture method. Inversion polymorphism parameters were also assessed in the same populations for a period of 3 years. Significant differences in the numbers of individuals were observed between sexes. This affected the effective population sizes between particular years. The ratio of the effective size over the cenzus dropped significantly in beech wood in 2 years. Although overall heterozygosity remained unchanged during the years in both habitats, frequencies of gene arrangements on five chromosomes show variability. After the bottleneck, some complex chromosomal arrangements appeared for the first time in both populations. Standard gene arrangements of chromosome A increased in frequency over the years in each habitat, while the complex arrangements remain rather stable and specific for each population. The results obtained indicate that the population structure may significantly change if the effective size of D. subobscura population is reduced, which is mostly related to microclimatic changes in habitats. Based on the results to date, monitoring of microevolutionary changes by using D. subobscura and its relatives seems a promising way to study the effects of global climate changes.

Local adaptation to environmental stress at different levels of genetic polymorphism in various plants and animals has been documented through evolution of heavy metal tolerance. We used samples of Drosophila subobscura populations from two differently polluted environments to analyze the change of chromosomal inversion polymorphism as genetic marker during laboratory exposure to lead. Exposure to environmental contamination can affect the genetic content within a particular inversion and produce targets for selection in populations from different environments. The aims were to discover whether the inversion polymorphism is shaped by the local natural environments, and if lead as a selection pressure would cause adaptive divergence of two populations during the multigenerational laboratory experiment. The results showed that populations retain signatures from past contamination events, and that heavy metal pollution can cause adaptive changes in population. Differences in inversion polymorphism between the two populations increased over generations under lead contamination in the laboratory. The inversion polymorphism of population originating from the more polluted natural environment was more stable during the experiment, both under conditions with and without lead. Therefore, results showed that inversion polymorphism as a genetic marker reflects a strong signature of adaptation to the local environment, and that historical demographic events and selection are important for both prediction of evolutionary potential and long-term viability of natural populations.

Adaptation to different environmental conditions is a natural phenomenon that potentially leads to population subdivision. We surveyed genetic differentiation in inversion polymorphism within populations of Drosophila subobscura sampled in three ecologically different forest communities. The analysis of inversion polymorphism revealed significant differences between some pairs of samples in some gene arrangement frequencies of the A, U, and E chromosomes and some karyotype combination frequencies of the U chromosome, but significant differentiation within populations was not observed. It cannot be decided which evolutionary forces are responsible for the observed variability in inversion polymorphism

Mitochondrial DNA variability of Drosophila subobscura Collin from Southeastern Serbia was studied with respect to Restriction Site Analysis (RSA) of complete mitochondrial genome and the nucleotide sequence of Cytochrome b (Cyt b) gene. The aim was to shed more light on the evolutionary forces that shape mtDNA variation of this species. Samples were collected from two sites in the foothills of the Balkan Mountains. No genetic differentiation was found between groups and most of the variation was observed within them. Restriction analysis revealed two main haplotypes and several rare ones. The sequencing of Cyt b gene showed larger number of haplotypes, among which, one is being the most common. The majority of singletons differed from the most frequent haplotype by one nucleotide change. Although some of the observed nucleotide differences may affect their host?s fitness, the observed pattern of variation is consistent with the seasonal fluctuations in population size. [Project of the Serbian Ministry of Education, Science and Technological Development, Grant no. 173012]

Populations of Drosophila subobscura from the urban area of Belgrade and from the locality, Deliblato, which is not under strong anthropogenic influence, were studied with the aim to characterize and compare their genetic structure by examining chromosomal inversion polymorphism. Additional analysis and comparison of this type of polymorphism with several other populations from different habitats in the central Balkans, was done. The obtained results indicate higher heterozygosity in the population from Belgrade. Despite being ecologically marginal and under strong and complex influences, this population did not show a decline in the number of inversions and it is not highly differentiated compared to the referent populations.