Comparative population genomics of latitudinal variation in Drosophila simulans and Drosophila melanogaster

Examples of clinal variation in phenotypes and genotypes across latitudinal transects have served as important models for understanding how spatially varying selection and demographic forces shape variation within species. Here, we examine the selective and demographic contributions to latitudinal variation through the largest comparative genomic study to date of Drosophila simulans and Drosophila melanogaster, with genomic sequence data from 382 individual fruit flies, collected across a spatial transect of 19 degrees latitude and at multiple time points over 2 years. Consistent with phenotypic studies, we find less clinal variation in D. simulans than D. melanogaster, particularly for the autosomes. Moreover, we find that clinally varying loci in D. simulans are less stable over multiple years than comparable clines in D. melanogaster. D. simulans shows a significantly weaker pattern of isolation by distance than D. melanogaster and we find evidence for a stronger contribution of migration to D. simulans population genetic structure. While population bottlenecks and migration can plausibly explain the differences in stability of clinal variation between the two species, we also observe a significant enrichment of shared clinal genes, suggesting that the selective forces associated with climate are acting on the same genes and phenotypes in D. simulans and D. melanogaster.     

Parallel effects of the inversion In(3R)Payne on body size across the North American and Australian clines in Drosophila melanogaster

Chromosomal inversions are thought to play a major role in climatic adaptation. In D. melanogaster, the cosmopolitan inversion In(3R)Payne exhibits latitudinal clines on multiple continents. As many fitness traits show similar clines, it is tempting to hypothesize that In(3R)P underlies observed clinal patterns for some of these traits. In support of this idea, previous work in Australian populations has demonstrated that In(3R)P affects body size but not development time or cold resistance. However, similar data from other clines of this inversion are largely lacking; finding parallel effects of In(3R)P across multiple clines would considerably strengthen the case for clinal selection. Here, we have analysed the phenotypic effects of In(3R)P in populations originating from the endpoints of the latitudinal cline along the North American east coast. We measured development time, egg‐to‐adult survival, several size‐related traits (femur and tibia length, wing area and shape), chill coma recovery, oxidative stress resistance and triglyceride content in homokaryon lines carrying In(3R)P or the standard arrangement. Our central finding is that the effects of In(3R)P along the North American cline match those observed in Australia: standard arrangement lines were larger than inverted lines, but the inversion did not influence development time or cold resistance. Similarly, In(3R)P did not affect egg‐to‐adult survival, oxidative stress resistance and lipid content. In(3R)P thus seems to specifically affect size traits in populations from both continents. This parallelism strongly suggests an adaptive pattern, whereby the inversion has captured alleles associated with growth regulation and clinal selection acts on size across both continents.     

Secondary contact and local adaptation contribute to genome‐wide patterns of clinal variation in Drosophila melanogaster

Populations arrayed along broad latitudinal gradients often show patterns of clinal variation in phenotype and genotype. Such population differentiation can be generated and maintained by both historical demographic events and local adaptation. These evolutionary forces are not mutually exclusive and can in some cases produce nearly identical patterns of genetic differentiation among populations. Here, we investigate the evolutionary forces that generated and maintain clinal variation genome‐wide among populations of Drosophila melanogaster sampled in North America and Australia. We contrast patterns of clinal variation in these continents with patterns of differentiation among ancestral European and African populations. Using established and novel methods we derive here, we show that recently derived North America and Australia populations were likely founded by both European and African lineages and that this hybridization event likely contributed to genome‐wide patterns of parallel clinal variation between continents. The pervasive effects of admixture mean that differentiation at only several hundred loci can be attributed to the operation of spatially varying selection using an F_ST outlier approach. Our results provide novel insight into the well‐studied system of clinal differentiation in D. melanogaster and provide a context for future studies seeking to identify loci contributing to local adaptation in a wide variety of organisms, including other invasive species as well as temperate endemics.     

Footprints of selection in wild populations of Bicyclus anynana along a latitudinal cline

One of the major questions in ecology and evolutionary biology is how variation in the genome enables species to adapt to divergent environments. Here, we study footprints of thermal selection in candidate genes in six wild populations of the afrotropical butterfly Bicyclus anynana sampled along a c. 3000 km latitudinal cline. We sequenced coding regions of 31 selected genes with known functions in metabolism, pigment production, development and heat shock responses. These include genes for which we expect a priori a role in thermal adaptation and, thus, varying selection pressures along a latitudinal cline, and genes we do not expect to vary clinally and can be used as controls. We identified amino acid substitution polymorphisms in 13 genes and tested these for clinal variation by correlation analysis of allele frequencies with latitude. In addition, we used two F_ST‐based outlier methods to identify loci with higher population differentiation than expected under neutral evolution, while accounting for potentially confounding effects of population structure and demographic history. Two metabolic enzymes of the glycolytic pathway, UGP and Treh, showed clinal variation. The same loci showed elevated population differentiation and were identified as significant outliers. We found no evidence of clines in the pigmentation genes, heat shock proteins and developmental genes. However, we identified outlier loci in more localized parts of the range in the pigmentation genes yellow and black. We discuss that the observed clinal variation and elevated population divergence in UGP and Treh may reflect adaptation to a geographic thermal gradient.     

Complex patterns of local adaptation in heat tolerance in Drosophila simulans from eastern Australia

Latitudinal clines are considered a powerful means of investigating evolutionary responses to climatic selection in nature. However, most clinal studies of climatic adaptation in Drosophila have involved species that contain cosmopolitan inversion polymorphisms that show clinal patterns themselves, making it difficult to determine whether the traits or inversions are under selection. Further, although climatic selection is unlikely to act on only one life stage in metamorphic organisms, a few studies have examined clinal patterns across life stages. Finally, clinal patterns of heat tolerance may also depend on the assay used. To unravel these potentially confounding effects on clinal patterns of thermal tolerance, we examined adult and larval heat tolerance traits in populations of Drosophila simulans from eastern Australia using static and dynamic (ramping 0.06 °C min^?1) assays. We also used microsatellites markers to clarify whether demographic factors or selection are responsible for population differentiation along clines. Significant cubic clinal patterns were observed for adult static basal, hardened and dynamic heat knockdown time and static basal heat survival in larvae. In contrast, static, hardened larval heat survival increased linearly with latitude whereas no clinal association was found for larval ramping survival. Significant associations between adult and larval traits and climatic variables, and low population differentiation at microsatellite loci, suggest a role for climatic selection, rather than demographic processes, in generating these clinal patterns. Our results suggest that adaptation to thermal stress may be species and life‐stage specific, complicating our efforts to understand the evolutionary responses to selection for increasing thermotolerance.     

Geographical differentiation of allozymic variability in natural Indian populations ofDrosophila melanogaster

Fourteen Indian populations ofD. melanogaster collected along a 22° latitudinal range were analyzed electrophoretically to compare the geographical patterns of allozymic variation at 13 loci. The data show higher genetic differentiation on the basis of moderate to higherF _ST values and significant statistical correlation of allelic frequencies at six polymorphic loci with latitude. The results add support to the hypothesis that the occurrence of parallel or complementary latitudinal clines across different continental populations provides evidence of natural selection maintaining such clinal variation.     

ALTITUDINAL CLINAL VARIATION IN WING SIZE AND SHAPE IN AFRICAN DROSOPHILA MELANOGASTER: ONE CLINE OR MANY?

Geographical patterns of morphological variation have been useful in addressing hypotheses about environmental adaptation. In particular, latitudinal clines in phenotypes have been studied in a number of Drosophila species. Some environmental conditions along latitudinal clines—for example, temperature—also vary along altitudinal clines, but these have been studied infrequently and it remains unclear whether these environmental factors are similar enough for convergence or parallel evolution. Most clinal studies in Drosophila have dealt exclusively with univariate phenotypes, allowing for the detection of clinal relationships, but not for estimating the directions of covariation between them. We measured variation in wing shape and size in D. melanogaster derived from populations at varying altitudes and latitudes across sub‐Saharan Africa. Geometric morphometrics allows us to compare shape changes associated with latitude and altitude, and manipulating rearing temperature allows us to quantify the extent to which thermal plasticity recapitulates clinal effects. Comparing effect vectors demonstrates that altitude, latitude, and temperature are only partly associated, and that the altitudinal shape effect may differ between Eastern and Western Africa. Our results suggest that selection responsible for these phenotypic clines may be more complex than just thermal adaptation.     

Seasonal variation in life history traits in two Drosophila species

Seasonal environmental heterogeneity is cyclic, persistent and geographically widespread. In species that reproduce multiple times annually, environmental changes across seasonal time may create different selection regimes that may shape the population ecology and life history adaptation in these species. Here, we investigate how two closely related species of Drosophila in a temperate orchard respond to environmental changes across seasonal time. Natural populations of Drosophila melanogaster and Drosophila simulans were sampled at four timepoints from June through November to assess seasonal change in fundamental aspects of population dynamics as well as life history traits. D. melanogaster exhibit pronounced change across seasonal time: early in the season, the population is inferred to be uniformly young and potentially represents the early generation following overwintering survivorship. D. melanogaster isofemale lines derived from the early population and reared in a common garden are characterized by high tolerance to a variety of stressors as well as a fast rate of development in the laboratory environment that declines across seasonal time. In contrast, wild D. simulans populations were inferred to be consistently heterogeneous in age distribution across seasonal collections; only starvation tolerance changed predictably over seasonal time in a parallel manner as in D. melanogaster. These results suggest fundamental differences in population and evolutionary dynamics between these two taxa associated with seasonal heterogeneity in environmental parameters and associated selection pressures.     

Patterns of transposable element variation and clinality in Drosophila

Natural populations often exist in spatially diverse environments and may experience variation in the strength and targets of natural selection across their ranges. Drosophila provides an excellent opportunity to study the effects of spatially varying selection in natural populations, as both Drosophila melanogaster and Drosophila simulans live across a wide range of environments in North America. Here, we characterize patterns of variation in transposable elements (TEs) from six populations of D. melanogaster and nine populations of D. simulans sampled from multiple latitudes across North America. We find a nearly twofold excess of TEs in D. melanogaster relative to D. simulans, with this difference largely driven by TEs segregating at the lowest and highest allele frequencies. We find no effect of latitude on either total TE abundance or average TE allele frequencies in either species. Moreover, we show that, as a class of mutations, the most common patterns of TE variation do not coincide with the sampled latitudinal gradient, nor are they consistent with local adaptation acting on environmental differences found in the most extreme latitudes. We also do not find a cline in ancestry for North American D. melanogaster—for either TEs or single nucleotide polymorphisms—suggesting a limited role for demography in shaping patterns of TE variation. Though we find little evidence for widespread clinality among TEs in Drosophila, this does not necessarily imply a limited role for TEs in adaptation. We discuss the need for improved models of adaptation to large‐scale environmental heterogeneity, and how these might be applied to TEs.     

An inversion supergene in Drosophila underpins latitudinal clines in survival traits

Chromosomal inversions often contribute to local adaptation across latitudinal clines, but the underlying selective mechanisms remain poorly understood. We and others have previously shown that a clinal inversion polymorphism in Drosophila melanogaster, In(3R)Payne, underpins body size clines along the North American and Australian east coasts. Here, we ask whether this polymorphism also contributes to clinal variation in other fitness‐related traits, namely survival traits (lifespan, survival upon starvation and survival upon cold shock). We generated homokaryon lines, either carrying the inverted or standard chromosomal arrangement, isolated from populations approximating the endpoints of the North American cline (Florida, Maine) and phenotyped the flies at two growth temperatures (18 °C, 25 °C). Across both temperatures, high‐latitude flies from Maine lived longer and were more stress resistant than low‐latitude flies from Florida, as previously observed. Interestingly, we find that this latitudinal pattern is partly explained by the clinal distribution of the In(3R)P polymorphism, which is at ~ 50% frequency in Florida but absent in Maine: inverted karyotypes tended to be shorter‐lived and less stress resistant than uninverted karyotypes. We also detected an interaction between karyotype and temperature on survival traits. As In(3R)P influences body size and multiple survival traits, it can be viewed as a ‘supergene’, a cluster of tightly linked loci affecting multiple complex phenotypes. We conjecture that the inversion cline is maintained by fitness trade‐offs and balancing selection across geography; elucidating the mechanisms whereby this inversion affects alternative, locally adapted phenotypes across the cline is an important task for future work.     

Complexity of the genetic basis of ageing in nature revealed by a clinal study of lifespan and methuselah,a gene for ageing,in Drosophila from eastern Australia

Clinal studies are a powerful tool for understanding the genetic basis of climatic adaptation. However, while clines in quantitative traits and genetic polymorphisms have been observed within and across continents, few studies have attempted to demonstrate direct links between them. The gene methuselah in Drosophila has been shown to have a major effect on stress response and longevity phenotypes based largely on laboratory studies of induced mutations in the mth gene. Clinal patterns in the most common mth haplotype and for lifespan (both increasing with latitude) have been observed in North American populations of D. melanogaster, implicating climatic selection. While these clinal patterns have led some to suggest that mth influences ageing in natural populations, limited evidence on the association between the two has so far been collected. Here, we describe a significant cline in the mth haplotype in eastern Australian D. melanogaster populations that parallel the cline in North America. We also describe a cline in mth gene expression. These findings further support the idea that mth is itself under selection. In contrast, we show that lifespan has a strong nonlinear clinal pattern, increasing southwards from the tropics, but then decreasing again from mid‐latitudes. Furthermore, in association studies, we find no evidence for a direct link between mth haplotype and lifespan. Thus, while our data support a role for mth variation being under natural selection, we found no link to naturally occurring variation in lifespan and ageing in Australian populations of D. melanogaster. Our results indicate that the mth locus likely has genetic background and environment‐specific effects.     

Conservation of social effects (Ψ) between two species of Drosophila despite reversal of sexual dimorphism

Indirect genetic effects (IGEs) describe the effect of the genes of social partners on the phenotype of a focal individual. Here, we measure indirect genetic effects using the “coefficient of interaction” (Ψ) to test whether Ψ evolved between Drosophila melanogaster and D. simulans. We compare Ψ for locomotion between ethanol and nonethanol environments in both species, but only D. melanogaster utilizes ethanol ecologically. We find that while sexual dimorphism for locomotion has been reversed in D. simulans, there has been no evolution of social effects between these two species. What did evolve was the interaction between genotype‐specific Ψ and the environment, as D. melanogaster varies unpredictably between environments and D. simulans does not. In this system, this suggests evolutionary lability of sexual dimorphism but a conservation of social effects, which brings forth interesting questions about the role of the social environment in sexual selection.     

Contrasting patterns of geographic variation in the cosmopolitan sibling species Drosophila melanogaster and Drosophila simulans

An electrophoretic study was carried out to compare the geographic pattern of genetic variation in Drosophila simulans with that of its sibling species, Drosophila melanogaster. An identical set of 32 gene-protein loci was studied in four geographically distant populations of D. simulans and two populations of D. melanogaster, all originating from Europe and Africa. The comparison yielded the following results: (1) tropical populations of D. simulans were, in terms of the number of unique alleles, average heterozygosity per locus, and percentage of loci polymorphic, more variable than conspecific-temperate populations; (2) some loci in both species showed interpopulation differences in allele frequencies that suggest latitudinal clines; and (3) temperate-tropical genetic differentiation between populations was much less in D. simulans than in D. melanogaster. Similar differences between these two species have previously been shown for chromosomal, quantitative, physiological, and middle-repetitive DNA variation. Estimates of N _m (number of migrants per generation) from the spatial distribution of rare alleles suggest that both species have similar levels of interpopulation gene flow. These observations lead us to propose two competing hypotheses: the low level of geographic differentiation in D. simulans is due to its evolutionarily recent worldwide colonization and, alternatively, D. simulans has a narrower niche than D. melanogaster. Geographic variation data on different genetic elements (e.g., mitochondrial DNA, two-dimensional proteins, etc.) are required before these hypotheses can be adequately tested.We thank the Natural Science and Engineering Research Council of Canada for financial support (Grant A0235 to R.S.S.).     

Sex‐biased transcriptome divergence along a latitudinal gradient

Sex‐dependent gene expression is likely an important genomic mechanism that allows sex‐specific adaptation to environmental changes. Among Drosophila species, sex‐biased genes display remarkably consistent evolutionary patterns; male‐biased genes evolve faster than unbiased genes in both coding sequence and expression level, suggesting sex differences in selection through time. However, comparatively little is known of the evolutionary process shaping sex‐biased expression within species. Latitudinal clines offer an opportunity to examine how changes in key ecological parameters also influence sex‐specific selection and the evolution of sex‐biased gene expression. We assayed male and female gene expression in Drosophila serrata along a latitudinal gradient in eastern Australia spanning most of its endemic distribution. Analysis of 11 631 genes across eight populations revealed strong sex differences in the frequency, mode and strength of divergence. Divergence was far stronger in males than females and while latitudinal clines were evident in both sexes, male divergence was often population specific, suggesting responses to localized selection pressures that do not covary predictably with latitude. While divergence was enriched for male‐biased genes, there was no overrepresentation of X‐linked genes in males. By contrast, X‐linked divergence was elevated in females, especially for female‐biased genes. Many genes that diverged in D. serrata have homologs also showing latitudinal divergence in Drosophila simulans and Drosophila melanogaster on other continents, likely indicating parallel adaptation in these distantly related species. Our results suggest that sex differences in selection play an important role in shaping the evolution of gene expression over macro‐ and micro‐ecological spatial scales.     

An experimental study of evolution in progress: Clines for quantitative traits in colonizing and Palearctic populations ofDrosophila

Summary Drosophila subobscura has recently colonized the American continent and is an excellent model for studying evolution in action. Previous analyses have shown that these colonizing populations have significant latitudinal clines for the frequencies of some chromosomal arrangements that parallel those clines found in the Old World. These results strongly suggest that this polymorphism is adaptive. In the present study, significant latitudinal clines for continuous morphometric variables (flies are larger in the north) have been detected in Old World populations ofD. subobscura. The adaptive nature of these clines is reinforced by the fact that parallel latitudinal clines for body size have also been detected inDrosophila obscura, a closely related sympatric species, as well as previously in otherDrosophila. On the other hand, no significant latitudinal clines for continuous morphometric traits, not even when using an overall size index, have been detected in colonizing populations ofD. subobscura. This is a rather surprising result given the number of generations that have elapsed since the species was detected in America and given that significant clines in chromosomal inversions are already established. Thus, the adaptive response of quantitative morphometric variables is not as rapid as that found for chromosomal inversions. Nevertheless, canonical correlation analysis suggests that significant latitudinal clines for body size might soon be detected in the American continent. The results obtained here are included in a projected time series with the aim of documenting size evolution in action.     

Clinal variation in developmental time and viability,and the response to thermal treatments in two species of Drosophila

The present study first addressed the question of whether developmental time (DT) and viability (VT) vary clinally along latitudinal and altitudinal gradients in Drosophila buzzatii, an autochthonous specialist and the generalist invasive Drosophila melanogaster. Coincident and positive altitudinal clines across species and, direct and inverse latitudinal clines were observed for DT in D. melanogaster and D. buzzatii, respectively. Opposing latitudinal and altitudinal clines were detected for VT only in D. melanogaster. The patterns observed along altitudinal gradients prompted us to investigate whether flies living at lowland and highland environments may respond differentially to thermal treatments consisting of regimes of constant and alternating temperatures. Flies reared at higher mean temperature developed faster than at lower mean temperature in both species. By contrast, the response in VT differed greatly between species. Highland D. melanogaster were more viable than lowland regardless the treatment, whereas, in D. buzzatii, highland flies were more viable than lowland in alternating thermal regimes and the reverse was true in treatments of constant temperature. The results obtained suggest that thermal amplitude may be an important factor that should be considered in investigations of thermal adaptation. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 233–245.     

Grow with the flow: a latitudinal cline in physiology is associated with more variable precipitation in Erythranthe cardinalis

Local adaptation is commonly observed in nature: organisms perform well in their natal environment, but poorly outside it. Correlations between traits and latitude, or latitudinal clines, are among the most common pieces of evidence for local adaptation, but identifying the traits under selection and the selective agents is challenging. Here, we investigated a latitudinal cline in growth and photosynthesis across 16 populations of the perennial herb Erythranthe cardinalis (Phrymaceae). Using machine learning methods, we identify interannual variation in precipitation as a likely selective agent: southern populations from more variable environments had higher photosynthetic rates and grew faster. We hypothesize that selection may favour a more annualized life history – grow now rather than save for next year – in environments where severe droughts occur more often. Thus, our study provides insight into how species may adapt if Mediterranean climates become more variable due to climate change.     

Are Drosophila preferences for yeasts stable or contextual?

Whether there are general mechanisms, driving interspecific chemical communication is uncertain. Saccharomycetaceae yeast and Drosophila fruit flies, both extensively studied research models, share the same fruit habitat, and it has been suggested their interaction comprises a facultative mutualism that is instigated and maintained by yeast volatiles. Using choice tests, experimental evolution, and volatile analyses, we investigate the maintenance of this relationship and reveal little consistency between behavioral responses of two isolates of sympatric Drosophila species. While D. melanogaster was attracted to a range of different Saccharomycetaceae yeasts and this was independent of fruit type, D. simulans preference appeared specific to a particular S. cerevisiae genotype isolated from a vineyard fly population. This response, however, was not consistent across fruit types and is therefore context‐dependent. In addition, D. simulans attraction to an individual S. cerevisiae isolate was pliable over ecological timescales. Volatile candidates were analyzed to identify a common signal for yeast attraction, and while D. melanogaster generally responded to fermentation profiles, D. simulans preference was more discerning and likely threshold‐dependent. Overall, there is no strong evidence to support the idea of bespoke interactions with specific yeasts for either of these Drosophila genotypes. Rather the data support the idea Drosophila are generally adapted to sense and locate fruits infested by a range of fungal microbes and/or that yeast–Drosophila interactions may evolve rapidly.     

Genome assembly and annotation of a Drosophila simulans strain from Madagascar

Drosophila simulans is a close relative of the genetic model D. melanogaster. Its worldwide distribution in combination with the absence of segregating chromosomal inversions makes this species an increasingly attractive model to study the molecular signatures of adaptation in natural and experimental populations. In an effort to improve the genomic resources for D. simulans, we assembled and annotated the genome of a strain originating from Madagascar (M252), the ancestral range of D. simulans. The comparison of the M252 genome to other available D. simulans assemblies confirmed its high quality, but also highlighted genomic regions that are difficult to assemble with NGS data. The annotation of M252 provides a clear improvement with alternative splicing for 52% of the multiple‐exon genes, UTRs for 70% of the genes, 225 novel genes and 781 pseudogenes being reported. We anticipate that the M252 genome will be a valuable resource for many research questions.     

Latitudinal Adh allozymic variation and ethanol tolerance in Indian populations of Drosophila melanogaster

Ten Indian geographical populations of D. melanogaster were assayed electrophoretically for Adh genic variation. The Indian geographical populations of D. melanogaster revealed significant clinal variation (3 % for 1 d? latitude) at Adh locus and Adh^F allelic frequency correlated significantly with increase in latitude. It was suggested that the abundance of secondary alcohols in the southern Indian tropical and humid environment might exert selective pressure favouring higher frequency of Adh^S allele. Patterns of ethanol utilization as well as ethanol tolerance were analyzed in larval and adult individuals of six geographical populations of D. melanogaster. Latitudinal variation in ethanol tolerance was observed in D. melanogaster populations from India. The parallel occurrence of latitudinal variation it Adh locus as well as ethanol tolerance in Indian geographical populations of D. melanogaster could be maintained by balancing natural selection varying spatially along the north-south axis of the Indian sub-continent.