Allozyme constitution of two standard strains ofDrosophila subobscura

Long-read based assembly and synteny analysis of a reference Drosophila subobscura genome reveals signatures of structural evolution driven by inversions recombination-suppression effects

BACKGROUND

Drosophila subobscura has long been a central model in evolutionary genetics. Presently, its use is hindered by the lack of a reference genome. To bridge this gap, here we used PacBio long-read technology, together with the available wealth of genetic marker information, to assemble and annotate a high-quality nuclear and complete mitochondrial genome for the species. With the obtained assembly, we performed the first synteny analysis of genome structure evolution in the subobscura subgroup.

RESULTS

We generated a highly-contiguous ~ 129 Mb-long nuclear genome, consisting of six pseudochromosomes corresponding to the six chromosomes of a female haploid set, and a complete 15,764 bp-long mitogenome, and provide an account of their numbers and distributions of codifying and repetitive content. All 12 identified paracentric inversion differences in the subobscura subgroup would have originated by chromosomal breakage and repair, with some associated duplications, but no evidence of direct gene disruptions by the breakpoints. Between lineages, inversion fixation rates were 10 times higher in continental D. subobscura than in the two small oceanic-island endemics D. guanche and D. madeirensis. Within D. subobscura, we found contrasting ratios of chromosomal divergence to polymorphism between the A sex chromosome and the autosomes.

CONCLUSIONS

We present the first high-quality, long-read sequencing of a D. subobscura genome. Our findings generally support genome structure evolution in this species being driven indirectly, through the inversions' recombination-suppression effects in maintaining sets of adaptive alleles together in the face of gene flow. The resources developed will serve to further establish the subobscura subgroup as model for comparative genomics and evolutionary indicator of global change.

Temperature-related genetic changes in laboratory populations of Drosophila subobscura: evidence against simple climatic-based explanations for latitudinal clines

Parallel latitudinal clines to the long-standing ones in the original Palearctic populations have independently evolved at different rates for chromosomal polymorphism and body size in South and North American populations of Drosophila subobscura since colonization around 25 years ago. This strongly suggests that (micro) evolutionary changes are largely predictable, but the underlying mechanisms are unknown. The putative role of temperature per se was investigated by using three sets of populations at each of three temperatures (13 degrees , 18 degrees , and 22 degrees C) spanning much of the tolerable range for this species. We found a lower chromosomal diversity at the warmest temperature; a quick and consistent shift in gene arrangement frequencies in response to temperature; an evolutionary decrease in wing size, mediated by both cell area and cell number, at 18 degrees C; no relationship between wing size and those inversions involved in latitudinal clines; and a shortening of the basal length of longitudinal vein IV relative to its total length with increasing standard dose. The trends for chromosomal polymorphism and body size were generally inconsistent from simple climatic-based explanations of worldwide latitudinal patterns. The findings are discussed in the light of available information on D. subobscura and results from earlier thermal selection experiments with various Drosophila species.

Swift laboratory thermal evolution of wing shape (but not size) in Drosophila subobscura and its relationship with chromosomal inversion polymorphism

Latitudinal clinal variation in wing size and shape has evolved in North American populations of Drosophila subobscura within about 20 years since colonization. While the size cline is consistent to that found in original European populations (and globally in other Drosophila species), different parts of the wing have evolved on the two continents. This clearly suggests that 'chance and necessity' are simultaneously playing their roles in the process of adaptation. We report here rapid and consistent thermal evolution of wing shape (but not size) that apparently is at odds with that suggestion. Three replicated populations of D. subobscura derived from an outbred stock at Puerto Montt (Chile) were kept at each of three temperatures (13, 18 and 22 degrees C) for 1 year and have diverged for 27 generations at most. We used the methods of geometric morphometrics to study wing shape variation in both females and males from the thermal stocks, and rates of genetic divergence for wing shape were found to be as fast or even faster than those previously estimated for wing size on a continental scale. These shape changes did not follow a neat linear trend with temperature, and are associated with localized shifts of particular landmarks with some differences between sexes. Wing shape variables were found to differ in response to male genetic constitution for polymorphic chromosomal inversions, which strongly suggests that changes in gene arrangement frequencies as a response to temperature underlie the correlated changes in wing shape because of gene-inversion linkage disequilibria. In fact, we also suggest that the shape cline in North America likely predated the size cline and is consistent with the quite different evolutionary rates between inversion and size clines. These findings cast strong doubts on the supposed 'unpredictability' of the geographical cline for wing traits in D. subobscura North American colonizing populations.

Quantitative-genetic analysis of wing form and bilateral asymmetry in isochromosomal lines ofDrosophila subobscura using Procrustes methods

Fluctuating asymmetry (FA) is often used as a measure of underlying developmental instability (DI), motivated by the idea that morphological variance is maladaptive. Whether or not DI has evolutionary potential is a highly disputed topic, marred by methodological problems and fuzzy prejudices. We report here some results from an ongoing study of the effects of karyotype, homozygosity and temperature on wing form and bilateral asymmetry using isochromosomal lines of Drosophila subobscura. Our approach uses the recently developed methodologies in geometric morphometrics to analyse shape configurations of landmarks within the standard statistical framework employed in studies of bilateral asymmetries, and we have extended these methods to partition the individual variation and the variation in asymmetries into genetic and environmental causal components. The analyses revealed temperature-dependent expression of genetic variation for wing size and wing shape, directional asymmetry (DA) of wing size, increased asymmetries at suboptimal temperature, and a transition from FA to DA in males as a result of increase in the rearing temperature. No genetic variation was generally detected for FA in our samples, but these are preliminary results because no crosses between lines were carried out and, therefore, the contribution of dominance was not taken into account. In addition, only a subset of the standing genetic variation was represented in the experiments.

Seasonal cycles of allozyme-by-chromosomal-inversion gametic disequilibrium in Drosophila subobscura

Allozyme loci are frequently found non randomly associated to the chromosomal inversions in which they are included in Drosophila. Two opposite views compete to explain strong allozyme-by-inversion gametic disequilibria: they result from natural selection or, conversely, merely represent remnants of associations accidentally established at the origin of inversions. Empirical efforts aimed at deciding between adaptive and historical scenarios have focused on the spatial distribution of disequilibria. Yet, the evolutionary significance of these associations remains uncertain. I report here the results of a time-series analysis of the seasonal variation of alleles at six allozyme loci (Acph, Lap, Pept-1, Ao, Mpi, and Xdh) in connection with the O chromosomal polymorphisms of D. subobscura. The findings were: (1) in the segment I of the O chromosome, Lap and Pept-1 allozymes changed seasonally in a cyclical fashion within the ST gene arrangement, but they changed erratically within the 3 + 4 gene configuration; (2) the frequencies of Lap1.11 and Pept-1(0.40) within ST dropped to their lowest values in early and late summer, respectively, when the seasonal level of the ST arrangement is lowest. Furthermore, Lap1.11 and Pept-1(0.40) covary with ST only within these seasons, yet in a fashion inconsistent with these alleles having a major influence on the dynamics of the inversion; (3) seasonal cycling of alleles within inversions were not detected at Acph, Ao, Mpi, and Xdh, yet these loci are nearly monomorphic at the study population, and/or their sampled series were shorter than those for Lap and Pept-1; and (4) simply monitoring allozyme frequencies separately for each inversion proved to be superior, for evidencing the seasonal cycles of the disequilibria, to the use of the D' coefficient of association. Observed seasonal cycles of allozymes within inversions likely reflect natural selection.