Multilocus analysis of nucleotide variation in Drosophila madeirensis, an endemic species of the Laurisilva forest in Madeira

Drosophila madeirensis is an endemic species of Madeira that inhabits the island Laurisilva forest. Nucleotide variation in D. madeirensis is analysed in six genomic regions and compared to that previously reported for the same regions in Drosophila subobscura, an abundant species in the Palearctic region that is closely related to D. madeirensis. The gene regions analysed are distributed along the O₃ inversion. The O₃ arrangement is monomorphic in D. madeirensis, and it was present in ancestral populations of D. subobscura but went extinct in this species after the origin of the derived O_ST and O₃₊₄ arrangements. Levels of nucleotide polymorphism in D. madeirensis are similar to those present in the O_ST and O₃₊₄ arrangements of D. subobscura, and the frequency spectrum is skewed towards rare variants. Purifying selection against deleterious nonsynonymous mutations is less effective in D. madeirensis. Although D. madeirensis and D. subobscura coexist at present in Madeira, no clear evidence of introgression was detected in the studied regions.     

Tracking changes in chromosomal arrangements and their genetic content during adaptation

There is considerable evidence for an adaptive role of inversions, but how their genetic content evolves and affects the subsequent evolution of chromosomal polymorphism remains controversial. Here, we track how life‐history traits, chromosomal arrangements and 22 microsatellites, within and outside inversions, change in three replicated populations of Drosophila subobscura for 30 generations of laboratory evolution since founding from the wild. The dynamics of fitness‐related traits indicated adaptation to the new environment concomitant with directional evolution of chromosomal polymorphism. Evidence of selective changes in frequency of inversions was obtained for seven of 23 chromosomal arrangements, corroborating a role for inversions in adaptation. The evolution of linkage disequilibrium between some microsatellites and chromosomes suggested that adaptive changes in arrangements involved changes in their genetic content. Several microsatellite alleles increased in frequency more than expected by drift in targeted inversions in all replicate populations. In particular, there were signs of selection in the O₃₊₄ arrangement favouring a combination of alleles in two loci linked to the inversion and changing along with it, although the lack of linkage disequilibrium between these loci precludes epistatic selection. Seven other alleles increased in frequency within inversions more than expected by drift, but were not in linkage disequilibrium with them. Possibly these alleles were hitchhiking along with alleles under selection that were not specific to those inversions. Overall, the selection detected on the genetic content of inversions, despite limited coverage of the genome, suggests that genetic changes within inversions play an important role in adaptation.     

Geographic pattern of allozyme and inversion polymorphism on chromosome O of Drosophila subobscura and its evolutionary origin

The chromosome O of Drosophila subobscura was studied with respect to genetic variability at three enzyme loci (Odh, Me, and Lap-4) and with respect to inversion polymorphism. Population samples were taken from seven localities along a north-south gradient from Sweden and Scotland to Tunisia.The chromosomal analysis revealed clinal frequency changes for gene arrangements from north to south. With the enzyme loci Odh and Me allele frequencies are similar throughout the distribution range. Both loci are located outside the common inversion complex O₃₊₄. On the other hand, frequency changes parallel to those of the gene arrangements were observed for the alleles of the Lap-locus. Nonrandom associations between Lap-alleles and the superimposed gene arrangements O_ST, O₃₊₄, O₃₊₄₊₈, and O₃₊₄₊₂₃ were found. These gene arrangements differ from each other with respect to allele frequencies at the Lap-locus but for a given gene arrangement the relative frequencies of Lap-alleles remain relatively constant along the north-south gradient. Thus allele frequencies at the Lap-locus can be predicted from inversion frequencies.These observations can be interpreted in such a way that the pattern of allozyme variation within gene arrangements is due to founder effects caused by the unique origin of inversions. The gene blocks in the different inversions seem to represent more or less separated gene pools. In polymorphic populations the coexistence of genetically differentiated inversions presumably gives rise to heterotic interaction.     

How much can history constrain adaptive evolution? A real‐time evolutionary approach of inversion polymorphisms in Drosophila subobscura

Chromosomal inversions are present in a wide range of animals and plants, having an important role in adaptation and speciation. Although empirical evidence of their adaptive value is abundant, the role of different processes underlying evolution of chromosomal polymorphisms is not fully understood. History and selection are likely to shape inversion polymorphism variation to an extent yet largely unknown. Here, we perform a real‐time evolution study addressing the role of historical constraints and selection in the evolution of these polymorphisms. We founded laboratory populations of Drosophila subobscura derived from three locations along the European cline and followed the evolutionary dynamics of inversion polymorphisms throughout the first 40 generations. At the beginning, populations were highly differentiated and remained so throughout generations. We report evidence of positive selection for some inversions, variable between foundations. Signs of negative selection were more frequent, in particular for most cold‐climate standard inversions across the three foundations. We found that previously observed convergence at the phenotypic level in these populations was not associated with convergence in inversion frequencies. In conclusion, our study shows that selection has shaped the evolutionary dynamics of inversion frequencies, but doing so within the constraints imposed by previous history. Both history and selection are therefore fundamental to predict the evolutionary potential of different populations to respond to global environmental changes.     

The banding pattern of polytene chromosomes of Drosophila guanche compared with that of D. subobscura

A detailed map of the salivary gland chromosomes of Drosophila guanche is presented and compared to the standard gene arrangements of D. subobscura. Generally, the polytene chromosomc banding patterns of the two species show a high degrce of homology. Only Segment I of the sex chromosome (Chromosome A) shows marked differences. The banding pattern proposed for this segment in D. guanche could have originated from a cluster of overlapping inversions including A₁ arrangement.     

New gene arrangements in natural populations of Drosophila subobscura

Four new gene arrangements were detected in natural populations of Drosophila subobscura near Zürich, Switzerland. Three of them are the result of inversions not recorded before, and one is a new combination of known gene arrangements. Another arrangement, I₃₊₄, was recorded for the first time in Switzerland. It may have been imported or possibly may have arisen de novo. It is suggested that inversion mutations are not so infrequent as so far supposed, if the presence of mobile elements is taken into account in D. subobscura.     

Individual inversions or their combinations: which is the main selective target in a natural population of Drosophila subobscura?

It is generally accepted that chromosomal inversions have been key elements in adaptation and speciation processes. In this context, Drosophila subobscura has been, and still is, an excellent model species due to its rich chromosomal polymorphism. In this species, many analyses from natural populations have demonstrated the adaptive potential of individual inversions (and their overlapped combinations, the so‐called arrangements). However, little information is available on the evolutionary role of combinations generated by inversions located in homologous and nonhomologous chromosomes. The aim of this research was to ascertain whether these combinations are also a target for natural selection. For this objective, we have studied the inversion composition of homologous and nonhomologous chromosomes from a D. subobscura sample collected in a well‐studied population, Mount Avala (Serbia). No significant deviation from H‐W expectations was detected, and when comparing particular karyotypic combinations, likelihood ratios close to 1 were obtained. Thus, it seems that for each pair of homologous chromosomes inversions no deviation from randomness was detected. Finally, no linkage disequilibrium was observed between inversions located in different chromosomes of the karyotype. For all these reasons, it can be assumed that, at the cytological level, the individual inversions rather than their combinations in different chromosomes are the main target of selection.     

Seasonal changes in chromosomal inversion polymorphism in a <Emphasis Type="Italic">Drosophila subobscura</Emphasis> natural population from a southeastern European continental refugium of the last glaciation period

Chromosomal inversion polymorphism was studied in a natural population of Drosophila subobscura from the Djerdap gorge of the Danube River, which is one of southeastern European continental refugia during the last glaciation period. The existence of changes in chromosomal inversion polymorphism during different periods of time was tested. The observed high abundance of chromosomal arrangements and lower IFR value allowed this population to be classified with central European populations. It was characterized by the O₇ and O_{3 + 4 + 5} arrangements, not found in other populations from the region analyzed. Within the region studied, the E_{1 + 2 + 9 + 12} chromosomal arrangement was found, apart from this population from a warm habitat, only in one population from a Mediterranean refugium. Many similarities have been observed between the population from the Djerdap gorge and other populations from the region analyzed, which indicates that postglacial migration may have started from this southeastern European refugium. Seasonal changes in the frequencies of the E_st and O_{3 + 4 + 7} chromosomal arrangements, E_st and O_{3 + 4 + 7}, being more frequent in spring and in summer, respectively, were found. The changes in the J chromosome (J_st/J_st), E chromosome (E_st/E_{1 + 2 + 9} and E₈/E_{1 + 2 + 9}), and O chromosome (O_st/O_{3 + 4 + 22}; O_{3 + 4}/O_{3 + 4}; and O_{3 + 4 + 1}/O_{3 + 4 + 7}), whose karyotypic combination frequencies were associated with seasonal changes, were also found. These changes are most probably associated with seasonal climatic changes in the population analyzed, which confirms that the chromosomal inversions in D. subobscura are adaptive. The text was submitted by the author in English.     

Gene activity of polytene chromosomes in Drosophila species of the obscura group

The polytene chromosome puffing patterns of Drosophila guanche were established and compared with those of Drosophila subobscura. A total of 150 loci, active in some of the 17 developmental stages studied, were described and 23 of them were found to form the characteristic puffing pattern of D. guanche. Taking into account the number of puffs as well as the gene activity of each chromosome and the total gene activity, D. guanche seems to be less active than D. subobscura. Although both species show a degree of homology in their puffing patterns lower than that found for sibling species, the degree of homology is stronger than that between species belonging to the same group but to different subgroups. Thus, D. guanche and D. subobscura must be considered as phylogenetically closely related species, belonging to the same subgroup.     

Genetic coadaptation in the chromosomal polymorphism of Drosophila subobscura II. Changes of gametic disequilibrium in experimental populations

Experimental populations were examined for temporal changes of gametic disequilibria between allozyme loci (Lap and Pept-1) and gene arrangements of the O chromosome of Drosophila subobscura (O _st and O ₃₊₄₊₇) under several environmental conditions. In the foundation of the experimental populations a genetic perturbation was carried out in order to test the relevance of the current hypotheses used to explain the allozyme-inversion associations observed in natural populations. Differential changes of gametic disequilibria were detected over generations under the different environmental conditions. Mere mechanical or stochastic factors cannot explain the results and natural selection is probably the major agent generating the detected gametic associations. The observations are interpreted as a proof of coadaptation of D. subobscura inversions.     

Contrasting patterns of phenotypic variation linked to chromosomal inversions in native and colonizing populations of Drosophila subobscura

In fewer than two decades after invading the Americas, the fly Drosophila subobscura evolved latitudinal clines for chromosomal inversion frequencies and wing size that are parallel to the long‐standing ones in native Palearctic populations. By sharp contrast, wing shape clines also evolved in the New World, but the relationship with latitude was opposite to that in the Old World. Previous work has suggested that wing trait differences among individuals are partially due to the association between chromosomal inversions and particular alleles which influence the trait under consideration. Furthermore, it is well documented that a few number of effective individuals founded the New World populations, which might have modified the biometrical effect of inversions on quantitative traits. Here we evaluate the relative contribution of chromosomal inversion clines in shaping the parallel clines in wing size and contrasting clines in wing shape in native and colonizing populations of the species. Our results reveal that inversion‐size and inversion‐shape associations in native and colonizing (South America) populations are generally different, probably due to the bottleneck effect. Contingent, unpredictable evolution was suggested as an explanation for the different details involved in the otherwise parallel wing size clines between Old and New World populations of D. subobscura. We challenge this assertion and conclude that contrasting wing shape clines came out as a correlated response of inversion clines that might have been predicted considering the genetic background of colonizers.     

Analysis of inbreeding in a colonizing population of Drosophila subobscura

An analysis of the effects of inbreeding on the genetic structure of a colonizing population of Drosophila subobscura has been carried out. Species of Drosophila, particularly D. subobscura, may have lethal alleles associated with chromosomal inversions and our aim was to assess the extent to which the genome is balanced in this way. The frequencies of chromosomal inversions were compared between a large population and a set of 72 lines that were maintained by brother-sister mating for 10 generations. Fisher's matrix method was used to calculate the expected homozygosity in these inbred lines for 5 allozyme loci (Aph, Hk-1, Lap, Odh, and Pept-1) used as markers of large chromosomal segments. Furthermore, the expected rates of fixation corresponding to these allozyme loci were also calculated. The results show that the amount of homozygosis observed did not differ significantly from expectations (with the corresponding loss of lines as a consequence of the reduction in viability). However, two deviations from strict neutrality were observed: there was a heterozygote excess at the Lap locus, and the frequency of the O ₅ inversion (always associated with a lethal gene in colonizing populations) was higher than expected.     

The role of genic selection in the establishment of inversion polymorphism in Drosophila subobscura

A model for explaining the establishment of newly arisen inversions in natural populations, in which the inverted segment may be selected for if its load of deleterious mutations is smaller than the average load of the noninverted segment in the population, is tested for Drosophila subobscura. The results show that for new inversions, originally with no deleterious alleles, the expected cumulative distribution of inversion lengths fits fairly well with the observed one. Therefore, genic selection may be an important cause of the establishment of newly arisen inversions in natural populations of D. subobscura. The applicability of the model to the maintenance of the inversion polymorphism present in this species; is discussed.     

Evolution of Chilean colonizing populations of <Emphasis Type="Italic">Drosophila subobscura</Emphasis>: lethal genes and chromosomal arrangements

Knowledge of the frequency, distribution, and fate of lethal genes in chromosomal inversions helps to illuminate the evolution of recently founded populations. We analyze the relationship between lethal genes and inversions in two colonizing populations of D. subobscura in Chile. In the ancestral Palearctic populations of this species, lethal genes seem distributed at random on chromosomes. But in colonizing American populations, some lethal genes are associated with specific chromosomal arrangements. Some of these associated lethals were detected only during the first stages of the colonization (O ₃₊₄₊₂ ), and never thereafter, whereas others have persisted (O ₃₊₄₊₇ and O₅). However, most lethal genes in American populations have been observed only once: they have arisen by novel mutation and soon disappear. Finally, recombination between different inversions has been observed in America. However, the persistence of lethal genes associated with the heterotic inversions O ₃₊₄₊₇ and O₅ could indicate that recombination inside these inversions is rare. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

STRUCTURE AND POPULATION GENETICS OF THE BREAKPOINTS OF A POLYMORPHIC INVERSION IN DROSOPHILA SUBOBSCURA

Drosophila subobscura is a paleartic species of the obscura group with a rich chromosomal polymorphism. To further our understanding on the origin of inversions and on how they regain variation, we have identified and sequenced the two breakpoints of a polymorphic inversion of D. subobscura—inversion 3 of the O chromosome—in a population sample. The breakpoints could be identified as two rather short fragments (～300 bp and 60 bp long) with no similarity to any known transposable element family or repetitive sequence. The presence of the ～300‐bp fragment at the two breakpoints of inverted chromosomes implies its duplication, an indication of the inversion origin via staggered double‐strand breaks. Present results and previous findings support that the mode of origin of inversions is neither related to the inversion age nor species‐group specific. The breakpoint regions do not consistently exhibit the lower level of variation within and stronger genetic differentiation between arrangements than more internal regions that would be expected, even in moderately small inversions, if gene conversion were greatly restricted at inversion breakpoints. Comparison of the proximal breakpoint region in species of the obscura group shows that this breakpoint lies in a small high‐turnover fragment within a long collinear region (～300 kb).     

INTROGRESSION IN THE DROSOPHILA SUBOBSCURA—D. MADEIRENSIS SISTER SPECIES: EVIDENCE OF GENE FLOW IN NUCLEAR GENES DESPITE MITOCHONDRIAL DIFFERENTIATION

Species hybridization, and thus the potential for gene flow, was once viewed as reproductive mistake. However, recent analysis based on large datasets and newly developed models suggest that gene exchange is not as rare as originally suspected. To investigate the history and speciation of the closely related species Drosophila subobscura, D. madeirensis, and D. guanche, we obtained polymorphism and divergence data for 26 regions throughout the genome, including the Y chromosome and mitochondrial DNA. We found that the D. subobscura X/autosome ratio of silent nucleotide diversity is significantly smaller than the 0.75 expected under neutrality. This pattern, if held genomewide, may reflect a faster accumulation of beneficial mutations on the X chromosome than on autosomes. We also detected evidence of gene flow in autosomal regions, while sex chromosomes remain distinct. This is consistent with the large X effect on hybrid male sterility seen in this system and the presence of two X chromosome inversions fixed between species. Overall, our data conform to chromosomal speciation models in which rearrangements are proposed to serve as gene flow barriers. Contrary to other observations in Drosophila, the mitochondrial genome appears resilient to gene flow in the presence of nuclear exchange.     

A molecular perspective on a complex polymorphic inversion system with cytological evidence of multiply reused breakpoints

Genome sequence comparison across the Drosophila genus revealed that some fixed inversion breakpoints had been multiply reused at this long timescale. Cytological studies of Drosophila inversion polymorphism had previously shown that, also at this shorter timescale, some breakpoints had been multiply reused. The paucity of molecularly characterized polymorphic inversion breakpoints has so far precluded contrasting whether cytologically shared breakpoints of these relatively young inversions are actually reused at the molecular level. The E chromosome of Drosophila subobscura stands out because it presents several inversion complexes. This is the case of the E_1+2+9+3 arrangement that originated from the ancestral E_st arrangement through the sequential accumulation of four inversions (E₁, E₂, E₉ and E₃) sharing some breakpoints. We recently identified the breakpoints of inversions E₁ and E₂, which allowed establishing reuse at the molecular level of the cytologically shared breakpoint of these inversions. Here, we identified and sequenced the breakpoints of inversions E₉ and E₃, because they share breakpoints at sections 58D and 64C with those of inversions E₁ and E₂. This has allowed establishing that E₉ and E₃ originated through the staggered-break mechanism. Most importantly, sequence comparison has revealed the multiple reuse at the molecular level of the proximal breakpoint (section 58D), which would have been used at least by inversions E₂, E₉ and E₃. In contrast, the distal breakpoint (section 64C) might have been only reused once by inversions E₁ and E₂, because the distal E₃ breakpoint is displaced >70 kb from the other breakpoint limits.     

Evolution of mating behavior between two populations adapting to common environmental conditions

Populations from the same species may be differentiated across contrasting environments, potentially affecting reproductive isolation among them. When such populations meet in a novel common environment, this isolation may be modified by biotic or abiotic factors. Curiously, the latter have been overlooked. We filled this gap by performing experimental evolution of three replicates of two populations of Drosophila subobscura adapting to a common laboratorial environment, and simulated encounters at three time points during this process. Previous studies showed that these populations were highly differentiated for several life‐history traits and chromosomal inversions. First, we show initial differentiation for some mating traits, such as assortative mating and male mating rate, but not others (e.g., female mating latency). Mating frequency increased during experimental evolution in both sets of populations. The assortative mating found in one population remained constant throughout the adaptation process, while disassortative mating of the other population diminished across generations. Additionally, differences in male mating rate were sustained across generations. This study shows that mating behavior evolves rapidly in response to adaptation to a common abiotic environment, although with a complex pattern that does not correspond to the quick convergence seen for life‐history traits.     

Efficient induction of heritable inversions in plant genomes using the CRISPR/Cas system

During the evolution of plant genomes, sequence inversions occurred repeatedly making the respective regions inaccessible for meiotic recombination and thus for breeding. Therefore, it is important to develop technologies that allow the induction of inversions within chromosomes in a directed and efficient manner. Using the Cas9 nuclease from Staphylococcus aureus (SaCas9), we were able to obtain scarless heritable inversions with high efficiency in the model plant Arabidopsis thaliana. Via deep sequencing, we defined the patterns of junction formation in wild‐type and in the non‐homologous end‐joining (NHEJ) mutant ku70‐1. Surprisingly, in plants deficient of KU70, inversion induction is enhanced, indicating that KU70 is required for tethering the local broken ends together during repair. However, in contrast to wild‐type, most junctions are formed by microhomology‐mediated NHEJ and thus are imperfect with mainly deletions, making this approach unsuitable for practical applications. Using egg‐cell‐specific expression of Cas9, we were able to induce heritable inversions at different genomic loci and at intervals between 3 and 18 kb, in the percentage range, in the T1 generation. By screening individual lines, inversion frequencies of up to the 10% range were found in T2. Most of these inversions had scarless junctions and were without any sequence change within the inverted region, making the technology attractive for use in crop plants. Applying our approach, it should be possible to reverse natural inversions and induce artificial ones to break or fix linkages between traits at will.     

Chromosomal polymorphism in natural populations of Drosophila subobscura near Zürich,Switzerland: a contribution to long-term comparisons

The chromosomal inversion polymorphism in natural populations of D. subobscura was investigated near Zürich, Seitzerland, in 1986 and 1987. These data are compared with earlier data from 1963, 1964 and 1984, collected in the same region. In all five acrocentric chromosomes, significant differences in the proportions of the dominant gene arrangements occurred; in the chromosomes A, I, O and U the standard order decreased from the sixties to the eighties, whereas the arrangements I₁, O₃₊₄, O₃₊₄₊₈ and U₁₊₂ became more frequent.