Climatic selection on genes and traits after a 100 year-old invasion: a critical look at the temperate-tropical clines in Drosophila melanogaster from eastern Australia

Across two continents: The genomic basis of environmental adaptation in house mice (Mus musculus domesticus) from the Americas

Replicated clines across environmental gradients can be strong evidence of adaptation. House mice (Mus musculus domesticus) were introduced to the Americas by European colonizers and are now widely distributed from Tierra del Fuego to Alaska. Multiple aspects of climate, such as temperature, vary predictably across latitude in the Americas. Past studies of North American populations across latitudinal gradients provided evidence of environmental adaptation in traits related to body size, metabolism, and behavior and identified candidate genes using selection scans. Here, we investigate genomic signals of environmental adaptation on a second continent, South America, and ask whether there is evidence of parallel adaptation across multiple latitudinal transects in the Americas. We first identified loci across the genome showing signatures of selection related to climatic variation in mice sampled across a latitudinal transect in South America, accounting for neutral population structure. Consistent with previous results, most candidate SNPs were in putatively regulatory regions. Genes that contained the most extreme outliers relate to traits such as body weight or size, metabolism, immunity, fat, eye function, and the cardiovascular system. We then compared these results with the results of analyses of published data from two transects in North America. While most candidate genes were unique to individual transects, we found significant overlap among candidate genes identified independently in the three transects. These genes are diverse, with functions relating to metabolism, immunity, cardiac function, and circadian rhythm, among others. We also found parallel shifts in allele frequency in candidate genes across latitudinal gradients. Finally, combining data from all three transects, we identified several genes associated with variation in body weight. Overall, our results provide strong evidence of shared responses to selection and identify genes that likely underlie recent environmental adaptation in house mice across North and South America.

From sub-Saharan Africa to China: Evolutionary history and adaptation of Drosophila melanogaster revealed by population genomics

Drosophila melanogaster is a widely used model organism for studying environmental adaptation. However, the genetic diversity of populations in Asia is poorly understood, leaving a notable gap in our knowledge of the global evolution and adaptation of this species. We sequenced genomes of 292 D. melanogaster strains from various ecological settings in China and analyzed them along with previously published genome sequences. We have identified six global genetic ancestry groups, despite the presence of widespread genetic admixture. The strains from China represent a unique ancestry group, although detectable differentiation exists among populations within China. We deciphered the global migration and demography of D. melanogaster, and identified widespread signals of adaptation, including genetic changes in response to insecticides. We validated the effects of insecticide resistance variants using population cage trials and deep sequencing. This work highlights the importance of population genomics in understanding the genetic underpinnings of adaptation, an effort that is particularly relevant given the deterioration of ecosystems.

Geographic and seasonal variation of the for gene reveal signatures of local adaptation in Drosophila melanogaster

In the early 1980s, the observation that Drosophila melanogaster larvae differed in their foraging behaviour laid the foundation for the work that would later lead to the discovery of the foraging gene (for) and its associated foraging phenotypes, rover and sitter. Since then, the molecular characterization of the for gene and our understanding of the mechanisms that maintain its phenotypic variants in the laboratory have progressed enormously. However, the significance and dynamics of such variation are yet to be investigated in nature. With the advent of next-generation sequencing, it is now possible to identify loci underlying the adaptation of populations in response to environmental variation. Here, I present the results of a genotype-environment association analysis that quantifies variation at the for gene among samples of D. melanogaster structured across space and time. These samples consist of published genomes of adult flies collected worldwide, and at least twice per site of collection (during spring and fall). Both an analysis of genetic differentiation based on Fs⁢t values and an analysis of population structure revealed an east-west gradient in allele frequency. This gradient may be the result of spatially varying selection driven by the seasonality of precipitation. These results support the hypothesis that different patterns of gene flow as expected under models of isolation by distance and potentially isolation by environment are driving genetic differentiation among populations. Overall, this study is essential for understanding the mechanisms underlying the evolution of foraging behaviour in D. melanogaster.

A lack of genetic diversity and minimal adaptive evolutionary divergence in introduced Mysis shrimp after 50 years

The successes of introduced populations in novel habitats often provide powerful examples of evolution and adaptation. In the 1950s, opossum shrimp (Mysis diluviana) individuals from Clearwater Lake in Minnesota, USA were transported and introduced to Twin Lakes in Colorado, USA by fisheries managers to supplement food sources for trout. Mysis were subsequently introduced from Twin Lakes into numerous lakes throughout Colorado. Because managers kept detailed records of the timing of the introductions, we had the opportunity to test for evolutionary divergence within a known time interval. Here, we used reduced representation genomic data to investigate patterns of genetic diversity, test for genetic divergence between populations, and for evidence of adaptive evolution within the introduced populations in Colorado. We found very low levels of genetic diversity across all populations, with evidence for some genetic divergence between the Minnesota source population and the introduced populations in Colorado. There was little differentiation among the Colorado populations, consistent with the known provenance of a single founding population, with the exception of the population from Gross Reservoir, Colorado. Demographic modeling suggests that at least one undocumented introduction from an unknown source population hybridized with the population in Gross Reservoir. Despite the overall low genetic diversity we observed, F ST outlier and environmental association analyses identified multiple loci exhibiting signatures of selection and adaptive variation related to elevation and lake depth. The success of introduced species is thought to be limited by genetic variation, but our results imply that populations with limited genetic variation can become established in a wide range of novel environments. From an applied perspective, the observed patterns of divergence between populations suggest that genetic analysis can be a useful forensic tool to determine likely sources of invasive species.

An Ancestral Balanced Inversion Polymorphism Confers Global Adaptation

Since the pioneering work of Dobzhansky in the 1930s and 1940s, many chromosomal inversions have been identified, but how they contribute to adaptation remains poorly understood. In Drosophila melanogaster, the widespread inversion polymorphism In(3R)Payne underpins latitudinal clines in fitness traits on multiple continents. Here, we use single-individual whole-genome sequencing, transcriptomics, and published sequencing data to study the population genomics of this inversion on four continents: in its ancestral African range and in derived populations in Europe, North America, and Australia. Our results confirm that this inversion originated in sub-Saharan Africa and subsequently became cosmopolitan; we observe marked monophyletic divergence of inverted and noninverted karyotypes, with some substructure among inverted chromosomes between continents. Despite divergent evolution of this inversion since its out-of-Africa migration, derived non-African populations exhibit similar patterns of long-range linkage disequilibrium between the inversion breakpoints and major peaks of divergence in its center, consistent with balancing selection and suggesting that the inversion harbors alleles that are maintained by selection on several continents. Using RNA-sequencing, we identify overlap between inversion-linked single-nucleotide polymorphisms and loci that are differentially expressed between inverted and noninverted chromosomes. Expression levels are higher for inverted chromosomes at low temperature, suggesting loss of buffering or compensatory plasticity and consistent with higher inversion frequency in warm climates. Our results suggest that this ancestrally tropical balanced polymorphism spread around the world and became latitudinally assorted along similar but independent climatic gradients, always being frequent in subtropical/tropical areas but rare or absent in temperate climates.

Selection and geography shape male reproductive tract transcriptomes in Drosophila melanogaster

Transcriptome analysis of several animal clades suggests that male reproductive tract gene expression evolves quickly. However, the factors influencing the abundance and distribution of within-species variation, the ultimate source of interspecific divergence, are poorly known. Drosophila melanogaster, an ancestrally African species that has recently spread throughout the world and colonized the Americas in the last roughly 100 years, exhibits phenotypic and genetic latitudinal clines on multiple continents, consistent with a role for spatially varying selection in shaping its biology. Nevertheless, geographic expression variation in the Americas is poorly described, as is its relationship to African expression variation. Here, we investigate these issues through the analysis of two male reproductive tissue transcriptomes [testis and accessory gland (AG)] in samples from Maine (USA), Panama, and Zambia. We find dramatic differences between these tissues in differential expression between Maine and Panama, with the accessory glands exhibiting abundant expression differentiation and the testis exhibiting very little. Latitudinal expression differentiation appears to be influenced by the selection of Panama expression phenotypes. While the testis shows little latitudinal expression differentiation, it exhibits much greater differentiation than the accessory gland in Zambia vs American population comparisons. Expression differentiation for both tissues is non-randomly distributed across the genome on a chromosome arm scale. Interspecific expression divergence between D. melanogaster and D. simulans is discordant with rates of differentiation between D. melanogaster populations. Strongly heterogeneous expression differentiation across tissues and timescales suggests a complex evolutionary process involving major temporal changes in the way selection influences expression evolution in these organs.

Rapid adaptation in a fast-changing world: Emerging insights from insect genomics

Many researchers have questioned the ability of biota to adapt to rapid anthropogenic environmental shifts. Here, we synthesize emerging genomic evidence for rapid insect evolution in response to human pressure. These new data reveal diverse genomic mechanisms (single locus, polygenic, structural shifts; introgression) underpinning rapid adaptive responses to a variety of anthropogenic selective pressures. While the effects of some human impacts (e.g. pollution; pesticides) have been previously documented, here we highlight startling new evidence for rapid evolutionary responses to additional anthropogenic processes such as deforestation. These recent findings indicate that diverse insect assemblages can indeed respond dynamically to major anthropogenic evolutionary challenges. Our synthesis also emphasizes the critical roles of genomic architecture, standing variation and gene flow in maintaining future adaptive potential. Broadly, it is clear that genomic approaches are essential for predicting, monitoring and responding to ongoing anthropogenic biodiversity shifts in a fast-changing world.

Within-population variation in body size plasticity in response to combined nutritional and thermal stress is partially independent from variation in development time

Ongoing climate change has forced animals to face changing thermal and nutritional environments. Animals can adjust to such combinations of stressors via plasticity. Body size is a key trait influencing organismal fitness, and plasticity in this trait in response to nutritional and thermal conditions varies among genetically diverse, locally adapted populations. The standing genetic variation within a population can also influence the extent of body size plasticity. We generated near-isogenic lines from a newly collected population of Drosophila melanogaster at the mid-point of east coast Australia and assayed body size for all lines in combinations of thermal and nutritional stress. We found that isogenic lines showed distinct underlying patterns of body size plasticity in response to temperature and nutrition that were often different from the overall population response. We then tested whether plasticity in development time could explain, and therefore regulate, variation in body size to these combinations of environmental conditions. We selected five genotypes that showed the greatest variation in response to combined thermal and nutritional stress and assessed the correlation between response of developmental time and body size. While we found significant genetic variation in development time plasticity, it was a poor predictor of body size among genotypes. Our results therefore suggest that multiple developmental pathways could generate genetic variation in body size plasticity. Our study emphasizes the need to better understand genetic variation in plasticity within a population, which will help determine the potential for populations to adapt to ongoing environmental change.

The genomic basis of copper tolerance in Drosophila is shaped by a complex interplay of regulatory and environmental factors

BACKGROUND

Escalation in industrialization and anthropogenic activity have resulted in an increase of pollutants released into the environment. Of these pollutants, heavy metals such as copper are particularly concerning due to their bio-accumulative nature. Due to its highly heterogeneous distribution and its dual nature as an essential micronutrient and toxic element, the genetic basis of copper tolerance is likely shaped by a complex interplay of genetic and environmental factors.

RESULTS

In this study, we utilized the natural variation present in multiple populations of Drosophila melanogaster collected across Europe to screen for variation in copper tolerance. We found that latitude and the degree of urbanization at the collection sites, rather than any other combination of environmental factors, were linked to copper tolerance. While previously identified copper-related genes were not differentially expressed in tolerant vs. sensitive strains, genes involved in metabolism, reproduction, and protease induction contributed to the differential stress response. Additionally, the greatest transcriptomic and physiological responses to copper toxicity were seen in the midgut, where we found that preservation of gut acidity is strongly linked to greater tolerance. Finally, we identified transposable element insertions likely to play a role in copper stress response.

CONCLUSIONS

Overall, by combining genome-wide approaches with environmental association analysis, and functional analysis of candidate genes, our study provides a unique perspective on the genetic and environmental factors that shape copper tolerance in natural D. melanogaster populations and identifies new genes, transposable elements, and physiological traits involved in this complex phenotype.

Desiccation resistance differences in Drosophila species can be largely explained by variations in cuticular hydrocarbons

Maintaining water balance is a universal challenge for organisms living in terrestrial environments, especially for insects, which have essential roles in our ecosystem. Although the high surface area to volume ratio in insects makes them vulnerable to water loss, insects have evolved different levels of desiccation resistance to adapt to diverse environments. To withstand desiccation, insects use a lipid layer called cuticular hydrocarbons (CHCs) to reduce water evaporation from the body surface. It has long been hypothesized that the water-proofing capability of this CHC layer, which can confer different levels of desiccation resistance, depends on its chemical composition. However, it is unknown which CHC components are important contributors to desiccation resistance and how these components can determine differences in desiccation resistance. In this study, we used machine-learning algorithms, correlation analyses, and synthetic CHCs to investigate how different CHC components affect desiccation resistance in 50 Drosophila and related species. We showed that desiccation resistance differences across these species can be largely explained by variation in CHC composition. In particular, length variation in a subset of CHCs, the methyl-branched CHCs (mbCHCs), is a key determinant of desiccation resistance. There is also a significant correlation between the evolution of longer mbCHCs and higher desiccation resistance in these species. Given that CHCs are almost ubiquitous in insects, we suggest that evolutionary changes in insect CHC components can be a general mechanism for the evolution of desiccation resistance and adaptation to diverse and changing environments.

Clinal variation as a tool to understand climate change

Clines are observable gradients that reflect continuous change in biological traits of species across geographical ranges. Clinal gradients could vary at geographic scales (latitude and altitude). Since clinal variations represent active genomic responses at the population level they (clines) provide an immense power to address questions related to climatic change. With the fast pace of climate change i.e. warming, populations are also likely to exhibit rapid responses; at both the phenotypic and genotypic levels. We seek to understand how clinal variation could be used to anticipate climatic responses using Drosophila, a pervasively used inter-disciplinary model system owing to its molecular repertoire. The genomic information coupled with the phenotypic variation greatly facilitates our understanding of the Drosophilidae response to climate change. We discuss traits associated with clinal variation at the phenotypic level as well as their underlying genetic regulators. Given prevailing climatic conditions and future projections for climate change, clines could emerge as monitoring tools to track the cross-talk between climatic variables and organisms.

Multivariate selection mediated by aridity predicts divergence of drought-resistant traits along natural aridity gradients of an invasive weed

Geographical variation in the environment underpins selection for local adaptation and evolutionary divergence among populations. Because many environmental conditions vary across species' ranges, identifying the specific environmental variables underlying local adaptation is profoundly challenging. We tested whether natural selection mediated by aridity predicts clinal divergence among invasive populations of capeweed (Arctotheca calendula) that established and spread across southern Australia during the last two centuries. Using common garden experiments with two environmental treatments (wet and dry) that mimic aridity conditions across capeweed's invasive range, we estimated clinal divergence and effects of aridity on fitness and multivariate phenotypic selection in populations sampled along aridity gradients in Australia. We show that: (1) capeweed populations have relatively high fitness in aridity environments similar to their sampling locations; (2) the magnitude and direction of selection strongly differs between wet and dry treatments, with drought stress increasing the strength of selection; and (3) differences in directional selection between wet and dry treatments predict patterns of clinal divergence across the aridity gradient, particularly for traits affecting biomass, flowering phenology and putative antioxidant expression. Our results suggest that aridity-mediated selection contributes to trait diversification among invasive capeweed populations, possibly facilitating the expansion of capeweed across southern Australia.

A wMel Wolbachia variant in Aedes aegypti from field-collected Drosophila melanogaster with increased phenotypic stability under heat stress

Mosquito‐borne diseases remain a major cause of morbidity and mortality. Population replacement strategies involving the wMel strain of Wolbachia are being used widely to control mosquito‐borne diseases. However, these strategies may be influenced by temperature because wMel is vulnerable to heat. wMel infections in Drosophila melanogaster are genetically diverse, but few transinfections of wMel variants have been generated in Aedes aegypti. Here, we successfully transferred a wMel variant (termed wMelM) originating from a field‐collected D. melanogaster into Ae. aegypti. The new wMelM variant (clade I) is genetically distinct from the original wMel transinfection (clade III), and there are no genomic differences between wMelM in its original and transinfected host. We compared wMelM with wMel in its effects on host fitness, temperature tolerance, Wolbachia density, vector competence, cytoplasmic incompatibility and maternal transmission under heat stress in a controlled background. wMelM showed a higher heat tolerance than wMel, likely due to higher overall densities within the mosquito. Both wMel variants had minimal host fitness costs, complete cytoplasmic incompatibility and maternal transmission, and dengue virus blocking under laboratory conditions. Our results highlight phenotypic differences between Wolbachia variants and wMelM shows potential as an alternative strain in areas with strong seasonal temperature fluctuations.

The rapid tempo of adaptation

Selection in fruit flies leads to fast adaption to seasonal changes.

On the Trail of Spatial Patterns of Genetic Variation

The accurate determination of the spatial trends on the variability of a species’ gene pool is essential to elucidate the underlying demographic-evolutionary events, thus helping to unravel the microevolutionary history of the population under study. Herein we present a new software called GenoCline, mainly addressed to detect genetic clines from allele, haplotype, and genome-wide data. This program package allows identifying the geographic orientation of clinal genetic variation through a system of iterative rotation of a virtual coordinate axis. Besides, GenoCline can perform complementary analyses to explore the potential origin of the genetic clines observed, including spatial autocorrelation, isolation by distance, centroid method, multidimensional scaling and Sammon projection. Among the advantages of this software is the ease in data entry and potential interconnection with other programs. Genetic and geographic data can be entered in spreadsheet table formatting (.xls), whereas genome-wide data can be imported in Eigensoft format. Genetic frequencies can also be exported in a format compatible with other programs dealing with population genetic and evolutionary biology analyses. All illustrations of results are saved in.svg format so that there will be high quality and easily editable vectorial graphs available for the researcher. Being implemented in Java, GenoCline is highly portable, thus working in different operating systems.

Clinal and seasonal changes are correlated in Drosophila melanogaster natural populations

Spatial and seasonal variations in the environment are ubiquitous. Environmental heterogeneity can affect natural populations and lead to covariation between environment and allele frequencies. Drosophila melanogaster is known to harbor polymorphisms that change both with latitude and seasons. Identifying the role of selection in driving these changes is not trivial, because nonadaptive processes can cause similar patterns. Given the environment changes in similar ways across seasons and along the latitudinal gradient, one promising approach may be to look for parallelism between clinal and seasonal changes. Here, we test whether there is a genome-wide correlation between clinal and seasonal changes, and whether the pattern is consistent with selection. Allele frequency estimates were obtained from pooled samples from seven different locations along the east coast of the United States, and across seasons within Pennsylvania. We show that there is a genome-wide correlation between clinal and seasonal variations, which cannot be explained by linked selection alone. This pattern is stronger in genomic regions with higher functional content, consistent with natural selection. We derive a way to biologically interpret these correlations and show that around 3.7% of the common, autosomal variants could be under parallel seasonal and spatial selection. Our results highlight the contribution of natural selection in driving fluctuations in allele frequencies in natural fly populations and point to a shared genomic basis to climate adaptation that happens over space and time in D. melanogaster.

Developmental and adult acclimation impact cold and drought survival of invasive tropical Drosophila kikkawai

Narrow distribution patterns of tropical Drosophila species are limited by lower resistance to cold or drought. In the invasive tropical Drosophila kikkawai, we tested whether developmental and adult acclimations at cooler temperatures could enhance its stress resistance level. Adult acclimation of winter collected body color morphs revealed a significant increase in the level of cold resistance. For light morph, its abundance during winter is not consistent with thermal-melanism hypothesis. However, higher cold acclimation capacity, as well as storage of energy metabolites could support its winter survival. In the wild-caught light and intermediate morphs, there is a lack of trade-off between cold and heat resistance but not in the case of dark morph. Developmental plasticity (15°C) resulted in the fivefold increase of cold survival at 0°C; and a twofold increase in desiccation resistance but a modest reduction (∼28–35%) in heat resistance as compared to morph strains reared at 25°C. Drought acclimation changes were significantly higher as compared with cold or heat pretreatment. We observed a trade-off between basal resistance and acclimation capacity for cold, heat, or drought resistance. For homeostatic energy balance, adult acclimation responses (cold versus drought; heat versus drought) caused compensatory plastic changes in the levels of proline or trehalose (shared patterns) but different patterns for total body lipids. In contrast, rapid cold or heat hardening-induced changes in energy metabolites were different as compared to acclimation. The ability of D. kikkawai to significantly increase stress tolerance through plasticity is likely to support its invasion potential.

Summary: In body color morphs of tropical Drosophila kikkawai, plasticity induced a higher level of resistance to cold and drought as well as three energy metabolites, which are likely to support its invasive potential.

Low levels of genetic differentiation with isolation by geography and environment in populations of Drosophila melanogaster from across China

The fruit fly, Drosophila melanogaster, is a model species in evolutionary studies. However, population processes of this species in East Asia are poorly studied. Here we examined the population genetic structure of D. melanogaster across China. There were 14 mitochondrial haplotypes with 10 unique ones out of 23 known from around the globe. Pairwise F_ST values estimated from 15 novel microsatellites ranged from 0 to 0.11, with geographically isolated populations showing the highest level of genetic uniqueness. STRUCTURE analysis identified high levels of admixture at both the individual and population levels. Mantel tests indicated a strong association between genetic distance and geographical distance as well as environmental distance. Full redundancy analysis (RDA) showed that independent effects of environmental conditions and geography accounted for 62.10% and 31.58% of the total explained genetic variance, respectively. When geographic variables were constrained in a partial RDA analysis, the environmental variables bio2 (mean diurnal air temperature range), bio13 (precipitation of the wettest month), and bio15 (precipitation seasonality) were correlated with genetic distance. Our study suggests that demographic history, geographical isolation, and environmental factors have together shaped the population genetic structure of D. melanogaster after its introduction into China.

Parallel and Population-specific Gene Regulatory Evolution in Cold-Adapted Fly Populations

Changes in gene regulation at multiple levels may comprise an important share of the molecular changes underlying adaptive evolution in nature. However, few studies have assayed within- and between-population variation in gene regulatory traits at a transcriptomic scale, and therefore inferences about the characteristics of adaptive regulatory changes have been elusive. Here, we assess quantitative trait differentiation in gene expression levels and alternative splicing (intron usage) between three closely related pairs of natural populations of Drosophila melanogaster from contrasting thermal environments that reflect three separate instances of cold tolerance evolution. The cold-adapted populations were known to show population genetic evidence for parallel evolution at the SNP level, and here we find evidence for parallel expression evolution between them, with stronger parallelism at larval and adult stages than for pupae. We also implement a flexible method to estimate cis- vs trans-encoded contributions to expression or splicing differences at the adult stage. The apparent contributions of cis- vs trans-regulation to adaptive evolution vary substantially among population pairs. While two of three population pairs show a greater enrichment of cis-regulatory differences among adaptation candidates, trans-regulatory differences are more likely to be implicated in parallel expression changes between population pairs. Genes with significant cis-effects are enriched for signals of elevated genetic differentiation between cold- and warm-adapted populations, suggesting that they are potential targets of local adaptation. These findings expand our knowledge of adaptive gene regulatory evolution and our ability to make inferences about this important and widespread process.

Effect of body lipid content is linked to nutritional adaptation in the acclimation responses of mesic-adapted Paederus to seasonal variations in desiccation stress

Desiccation stress causes mesic-adapted arthropods to lose their body water content. However, mesic-adapted Paederus beetles can survive over prolonged periods under dry field conditions, suggesting that these beetles adopt an array of water conservation mechanisms. We investigated the water balance mechanisms of field-collected Paederus adults over a 14-month sampling period. We also assessed their nutritional adaptations by performing a stable isotope analysis to examine their diet. The water loss rate (WLR) of the beetles was significantly associated with the rice crop cycle and saturation deficit. The cuticular permeability (CP) of adult beetles was maintained at < 30 µg cm^-2h^-1 mmHg^-1; however, CP increased significantly with the WLR. This result indicates that CP might play a minor role in reducing excessive water loss in beetles. The beetles' body water content and percentage total body water content increased when the WLR was high. Trehalose, glucose, and glycogen did not appear to play a central role in enhancing the water reserves in the insects. The body lipid content ranged from 0.22 ± 0.06 to 0.87 ± 0.07 mg and was negatively associated with the WLR. This association indicates that the increase in internal metabolic water was mediated by lipid catabolism. Stable isotope analysis results revealed that the Paederus beetles shifted their diet to carbohydrate-rich plants when the saturation deficit increased and the associated WLR reached its peak; otherwise, they consumed a high amount of staple carbohydrate-poor herbivore prey. The accumulation of energy reserves in the form of lipids through seasonal dietary shifts may exert major effects on the survival and population success of mesic-adapted Paederus beetles.

Temperature, rainfall and wind variables underlie environmental adaptation in natural populations of Drosophila melanogaster

While several studies in a diverse set of species have shed light on the genes underlying adaptation, our knowledge on the selective pressures that explain the observed patterns lags behind. Drosophila melanogaster is a valuable organism to study environmental adaptation because this species originated in Southern Africa and has recently expanded worldwide, and also because it has a functionally well-annotated genome. In this study, we aimed to decipher which environmental variables are relevant for adaptation of D. melanogaster natural populations in Europe and North America. We analysed 36 whole-genome pool-seq samples of D. melanogaster natural populations collected in 20 European and 11 North American locations. We used the BayPass software to identify single nucleotide polymorphisms (SNPs) and transposable elements (TEs) showing signature of adaptive differentiation across populations, as well as significant associations with 59 environmental variables related to temperature, rainfall, evaporation, solar radiation, wind, daylight hours, and soil type. We found that in addition to temperature and rainfall, wind related variables are also relevant for D. melanogaster environmental adaptation. Interestingly, 23%-51% of the genes that showed significant associations with environmental variables were not found overly differentiated across populations. In addition to SNPs, we also identified 10 reference transposable element insertions associated with environmental variables. Our results showed that genome-environment association analysis can identify adaptive genetic variants that are undetected by population differentiation analysis while also allowing the identification of candidate environmental drivers of adaptation.

Climate stress resistance in male Queensland fruit fly varies among populations of diverse geographic origins and changes during domestication

Background

The highly polyphagous Queensland fruit fly (Bactrocera tryoni Froggatt) expanded its range substantially during the twentieth century and is now the most economically important insect pest of Australian horticulture, prompting intensive efforts to develop a Sterile Insect Technique (SIT) control program. Using a “common garden” approach, we have screened for natural genetic variation in key environmental fitness traits among populations from across the geographic range of this species and monitored changes in those traits induced during domestication.

Results

Significant variation was detected between the populations for heat, desiccation and starvation resistance and wing length (as a measure of body size). Desiccation resistance was correlated with both starvation resistance and wing length. Bioassay data for three resampled populations indicate that much of the variation in desiccation resistance reflects persistent, inherited differences among the populations. No latitudinal cline was detected for any of the traits and only weak correlations were found with climatic variables for heat resistance and wing length. All three stress resistance phenotypes and wing length changed significantly in certain populations with ongoing domestication but there was also a strong population by domestication interaction effect for each trait.

Conclusions

Ecotypic variation in heat, starvation and desiccation resistance was detected in Australian Qfly populations, and these stress resistances diminished rapidly during domestication. Our results indicate a need to select source populations for SIT strains which have relatively high climatic stress resistance and to minimise loss of that resistance during domestication.

The Evolution of Phenotypic Plasticity in Response to Temperature Stress

Phenotypic plasticity is the ability of a single genotype to produce different phenotypes in response to environmental variation. The importance of phenotypic plasticity in natural populations and its contribution to phenotypic evolution during rapid environmental change is widely debated. Here, we show that thermal plasticity of gene expression in natural populations is a key component of its adaptation: evolution to novel thermal environments increases ancestral plasticity rather than mean genetic expression. We determined the evolution of plasticity in gene expression by conducting laboratory natural selection on a Drosophila simulans population in hot and cold environments. After more than 60 generations in the hot environment, 325 genes evolved a change in plasticity relative to the natural ancestral population. Plasticity increased in 75% of these genes, which were strongly enriched for several well-defined functional categories (e.g., chitin metabolism, glycolysis, and oxidative phosphorylation). Furthermore, we show that plasticity in gene expression of populations exposed to different temperatures is rather similar across species. We conclude that most of the ancestral plasticity can evolve further in more extreme environments.

Pervasive Effects of Wolbachia on Host Temperature Preference

Heritable symbionts can modify a range of ecologically important host traits, including behavior. About half of all insect species are infected with maternally transmitted Wolbachia, a bacterial endosymbiont known to alter host reproduction, nutrient acquisition, and virus susceptibility. Here, we broadly test the hypothesis that Wolbachia modifies host behavior by assessing the effects of eight different Wolbachia strains on the temperature preference of six Drosophila melanogaster subgroup species. Four of the seven host genotypes infected with A-group Wolbachia strains (wRi in Drosophila simulans, wHa in D. simulans, wSh in Drosophila sechellia, and wTei in Drosophila teissieri) prefer significantly cooler temperatures relative to uninfected genotypes. Contrastingly, when infected with divergent B-group wMau, Drosophila mauritiana prefers a warmer temperature. For most strains, changes to host temperature preference do not alter Wolbachia titer. However, males infected with wSh and wTei tend to experience an increase in titer when shifted to a cooler temperature for 24 h, suggesting that Wolbachia-induced changes to host behavior may promote bacterial replication. Our results indicate that Wolbachia modifications to host temperature preference are likely widespread, which has important implications for insect thermoregulation and physiology. Understanding the fitness consequences of these Wolbachia effects is crucial for predicting evolutionary outcomes of host-symbiont interactions, including how Wolbachia spreads to become common.IMPORTANCE Microbes infect a diversity of species, influencing the performance and fitness of their hosts. Maternally transmitted Wolbachia bacteria infect most insects and other arthropods, making these bacteria some of the most common endosymbionts in nature. Despite their global prevalence, it remains mostly unknown how Wolbachia influence host physiology and behavior to proliferate. We demonstrate pervasive effects of Wolbachia on Drosophila temperature preference. Most hosts infected with A-group Wolbachia prefer cooler temperatures, whereas the one host species infected with divergent B-group Wolbachia prefers warmer temperatures, relative to uninfected genotypes. Changes to host temperature preference generally do not alter Wolbachia abundance in host tissues, but for some A-group strains, adult males have increased Wolbachia titer when shifted to a cooler temperature. This suggests that Wolbachia-induced changes to host behavior may promote bacterial replication. Our results help elucidate the impact of endosymbionts on their hosts amid the global Wolbachia pandemic.

Phenotypic coupling of sleep and starvation resistance evolves in D. melanogaster

BACKGROUND

One hypothesis for the function of sleep is that it serves as a mechanism to conserve energy. Recent studies have suggested that increased sleep can be an adaptive mechanism to improve survival under food deprivation in Drosophila melanogaster. To test the generality of this hypothesis, we compared sleep and its plastic response to starvation in a temperate and tropical population of Drosophila melanogaster.

RESULTS

We found that flies from the temperate population were more starvation resistant, and hypothesized that they would engage in behaviors that are considered to conserve energy, including increased sleep and reduced movement. Surprisingly, temperate flies slept less and moved more when they were awake compared to tropical flies, both under fed and starved conditions, therefore sleep did not correlate with population-level differences in starvation resistance. In contrast, total sleep and percent change in sleep when starved were strongly positively correlated with starvation resistance within the tropical population, but not within the temperate population. Thus, we observe unexpectedly complex relationships between starvation and sleep that vary both within and across populations. These observations falsify the simple hypothesis of a straightforward relationship between sleep and energy conservation. We also tested the hypothesis that starvation is correlated with metabolic phenotypes by investigating stored lipid and carbohydrate levels, and found that stored metabolites partially contributed towards variation starvation resistance.

CONCLUSIONS

Our findings demonstrate that the function of sleep under starvation can rapidly evolve on short timescales and raise new questions about the physiological correlates of sleep and the extent to which variation in sleep is shaped by natural selection.

Genomic Analysis of European Drosophila melanogaster Populations Reveals Longitudinal Structure, Continent-Wide Selection, and Previously Unknown DNA Viruses

Genetic variation is the fuel of evolution, with standing genetic variation especially important for short-term evolution and local adaptation. To date, studies of spatiotemporal patterns of genetic variation in natural populations have been challenging, as comprehensive sampling is logistically difficult, and sequencing of entire populations costly. Here, we address these issues using a collaborative approach, sequencing 48 pooled population samples from 32 locations, and perform the first continent-wide genomic analysis of genetic variation in European Drosophila melanogaster. Our analyses uncover longitudinal population structure, provide evidence for continent-wide selective sweeps, identify candidate genes for local climate adaptation, and document clines in chromosomal inversion and transposable element frequencies. We also characterize variation among populations in the composition of the fly microbiome, and identify five new DNA viruses in our samples.

Rejection of the beneficial acclimation hypothesis (BAH) for short term heat acclimation in Drosophila nepalensis

Beneficial acclimation hypothesis (BAH) is the phenotypic plasticity in response to changing environments which enables organisms to enhance their fitness. In recent years, however, BAH has received vigorous criticism and is still debatable. In this study, we tested thermal hardiness phenotypes (melanization, chill coma recovery, heat knockdown and percentage survival) on adult and pre-adult stages of Drosophila nepalensis, reared in different thermal environments (14, 17, 21 and 25 °C) to check whether increasing natural surrounding temperature and acclimation limit towards environmental change is detrimental or beneficial. Results showed that rearing D. nepalensis at higher temperatures (21 and 25 °C) reduces its melanization and cold hardiness but improves heat knockdown times. When temperature was raised to 26.2 °C (0.6 °C above the upper thermal maxima), to determine the short-term acclimation effects, survival and fitness of adults diminished approximately 1.5 to 2 folds. These results suggest that D. nepalensis has long-term developmental acclimation to both heat and cold which would be extremely beneficial as temperatures and climates alter in the region due to global warming. However, a lack of short-term heat acclimation suggests that rapid shifts in thermal extreme could be detrimental to D. nepalensis.

Life-History Evolution and the Genetics of Fitness Components in Drosophila melanogaster

Life-history traits or "fitness components"-such as age and size at maturity, fecundity and fertility, age-specific rates of survival, and life span-are the major phenotypic determinants of Darwinian fitness. Analyzing the evolution and genetics of these phenotypic targets of selection is central to our understanding of adaptation. Due to its simple and rapid life cycle, cosmopolitan distribution, ease of maintenance in the laboratory, well-understood evolutionary genetics, and its versatile genetic toolbox, the "vinegar fly" Drosophila melanogaster is one of the most powerful, experimentally tractable model systems for studying "life-history evolution." Here, I review what has been learned about the evolution and genetics of life-history variation in D. melanogaster by drawing on numerous sources spanning population and quantitative genetics, genomics, experimental evolution, evolutionary ecology, and physiology. This body of work has contributed greatly to our knowledge of several fundamental problems in evolutionary biology, including the amount and maintenance of genetic variation, the evolution of body size, clines and climate adaptation, the evolution of senescence, phenotypic plasticity, the nature of life-history trade-offs, and so forth. While major progress has been made, important facets of these and other questions remain open, and the D. melanogaster system will undoubtedly continue to deliver key insights into central issues of life-history evolution and the genetics of adaptation.

Evolutionary and functional genetics of insect diapause: a call for greater integration

Diapause in response to seasonality is an important model for rapid evolutionary adaptation that is highly genetically variable, and experiences strong natural selection. Forward genetic methods using various genomic and transcriptomic approaches have begun to characterize the genetic architecture and candidate genes underlying diapause evolution. Largely in parallel, reverse genetic studies have identified functional roles for candidate genes that may or may not be genetically variable. We illustrate the disconnect between the evolutionary and physiological literature using a suite of studies of the role of the circadian clock in diapause regulation. These extensive studies in two different disciplines provide excellent opportunities for integration, which should facilitate rapid progress in understanding both the regulation and evolution of diapause.

Recent climate-driven ecological change across a continent as perceived through local ecological knowledge

Documenting effects of climate change is an important step towards designing mitigation and adaptation responses. Impacts of climate change on terrestrial biodiversity and ecosystems have been well-documented in the Northern Hemisphere, but long-term data to detect change in the Southern Hemisphere are limited, and some types of change are generally difficult to measure. Here we present a novel approach using local ecological knowledge to facilitate a continent-scale view of climate change impacts on terrestrial biodiversity and ecosystems that people have perceived in Australia. We sought local knowledge using a national web-based survey, targeting respondents with close links to the environment (e.g. farmers, ecologists), and using a custom-built mapping tool to ask respondents to describe and attribute recent changes they had observed within an area they knew well. Results drawn from 326 respondents showed that people are already perceiving simple and complex climate change impacts on hundreds of species and ecosystems across Australia, significantly extending the detail previously reported for the continent. While most perceived trends and attributions remain unsubstantiated, >35 reported anecdotes concurred with examples in the literature, and >20 were reported more than once. More generally, anecdotes were compatible with expectations from global climate change impact frameworks, including examples across the spectrum from organisms (e.g. increased mortality in >75 species), populations (e.g. changes in recruitment or abundance in >100 species, phenological change in >50 species), and species (e.g. >80 species newly arriving or disappearing), to communities and landscapes (e.g. >50 examples of altered ecological interactions). The overarching pattern indicated by the anecdotes suggests that people are more often noticing climate change losers (typically native species) than winners in their local areas, but with observations of potential ‘adaptation in action’ via compositional and phenological change and through arrivals and range shifts (particularly for native birds and exotic plants). A high proportion of climate change-related anecdotes also involved cumulative or interactive effects of land use. We conclude that targeted elicitation of local ecological knowledge about climate change impacts can provide a valuable complement to data-derived knowledge, substantially extending the volume of explicit examples and offering a foundation for further investigation.

A clinal polymorphism in the insulin signaling transcription factor foxo contributes to life-history adaptation in Drosophila

A fundamental aim of adaptation genomics is to identify polymorphisms that underpin variation in fitness traits. In Drosophila melanogaster, latitudinal life-history clines exist on multiple continents and make an excellent system for dissecting the genetics of adaptation. We have previously identified numerous clinal single-nucleotide polymorphism in insulin/insulin-like growth factor signaling (IIS), a pathway known from mutant studies to affect life history. However, the effects of natural variants in this pathway remain poorly understood. Here we investigate how two clinal alternative alleles at foxo, a transcriptional effector of IIS, affect fitness components (viability, size, starvation resistance, fat content). We assessed this polymorphism from the North American cline by reconstituting outbred populations, fixed for either the low- or high-latitude allele, from inbred DGRP lines. Because diet and temperature modulate IIS, we phenotyped alleles across two temperatures (18°C, 25°C) and two diets differing in sugar source and content. Consistent with clinal expectations, the high-latitude allele conferred larger body size and reduced wing loading. Alleles also differed in starvation resistance and expression of insulin-like receptor, a transcriptional target of FOXO. Allelic reaction norms were mostly parallel, with few GxE interactions. Together, our results suggest that variation in IIS makes a major contribution to clinal life-history adaptation.

Linking local adaptation with the evolution of sex differences

Many conspicuous forms of evolutionary diversity occur within species. Two prominent examples include evolutionary divergence between populations differentially adapted to their local environments (local adaptation), and divergence between females and males in response to sex differences in selection (sexual dimorphism sensu lato). These two forms of diversity have inspired vibrant research programmes, yet these fields have largely developed in isolation from one another. Nevertheless, conceptual parallels between these research traditions are striking. Opportunities for local adaptation strike a balance between local selection, which promotes divergence, and gene flow-via dispersal and interbreeding between populations-which constrains it. Sex differences are similarly constrained by fundamental features of inheritance that mimic gene flow. Offspring of each sex inherit genes from same-sex and opposite-sex parents, leading to gene flow between each differentially selected half of the population, and raising the question of how sex differences arise and are maintained. This special issue synthesizes and extends emerging research at the interface between the research traditions of local adaptation and sex differences. Each field can promote understanding of the other, and interactions between local adaptation and sex differences can generate new empirical predictions about the evolutionary consequences of selection that varies across space, time, and between the sexes.This article is part of the theme issue 'Linking local adaptation with the evolution of sex differences'.

Changes in lipid classes of Drosophila melanogaster in response to selection for three stress traits

Laboratory selection on environmental stress traits is an evolutionary approach that is informative in the context of understanding stress adaptation. Here we characterize changes in a lipidome of Drosophila melanogaster in lines selected for increased heat (elevated heat knockdown refractoriness), cold (decreased time to recover from chill-coma) and desiccation survival. Selection for desiccation resistance resulted in changes in multiple lipid classes used to characterize a lipidome. This included a decrease in triacylglycerols (TAGs) which is relevant to interpretation of storage lipid levels in previous D. melanogaster desiccation survival selection experiments. Chill-coma recovery rate selection was expected to show extensive changes in lipid classes, but only phosphatidic acids exhibited significant change. Selection for increased heat knockdown resistance resulted in a substantial change in the abundance of a class of lipids (diacylglycerols) which could play a role in mediating the heat shock response or result in an increase in neutral lipid mobilization.

Implications of existing local (mal)adaptations for ecological forecasting under environmental change

Standing genetic variation represents a genetic load on population fitness but can also support a rapid response to short-term environmental change, and the greatest potential source of such standing genetic variation typically exists among locally adapted populations living along an environmental gradient. Here, we develop a spatially explicit simulation model to quantify the contribution of existing genetic variation arising from migration-mutation-selection-drift balance to time to extinction under environmental change. Simulations reveal that local adaptation across a species range associated with an underlying environmental gradient could extend time to extinction by nearly threefold irrespective of the rate of environmental change. The potential for preadapted alleles to increase the rate of adaptation changes the relative importance of established extinction risk factors; in particular, it reduced the importance of the breadth of environmental tolerance and it increased the relative importance of fecundity. Although migration of preadapted alleles generally increased persistence time, it decreased it at rates of environmental change close to the critical rate of change by creating a population bottleneck, which ultimately limited the rate at which de novo mutations could arise. An analysis of the extinction dynamics further revealed that one consequence of gene flow is the potential to maximize population growth rate in at least part of the species range, which is likely to have consequences for forecasting the consequences of ecological interactions. Our study shows that predictions of persistence time change fundamentally when existing local adaptations are explicitly taken into account, underscoring the need to preserve and manage genetic diversity.

Evolutionary insights from large scale resequencing datasets in Drosophila melanogaster

Drosophila melanogaster has long been used as an evolutionary model system. Its small genome size, well-annotated genome, and ease of sampling, also makes it a choice species for genome resequencing studies. Hundreds of genomic samples from populations worldwide are available and are currently being used to tackle a wide range of evolutionary questions. In this review, we focused on three insights that have increased our understanding of the evolutionary history of this species, and that have implications for the study of evolutionary processes in other species as well. Because of technical limitations, most of the studies so far have focused on SNP variants. However, long-read sequencing techniques should allow us in the near future to include other type of genomic variants that also influence genome evolution.

Functional Analysis of a Putative Target of Spatially Varying Selection in the Menin1 Gene of Drosophila melanogaster

While significant effort has been devoted to investigating the potential influence of spatially varying selection on genomic variation, relatively little effort has been devoted to experimental analysis of putative variants or genes experiencing such selection. Previous population genetic work identified an amino acid polymorphism in the Mnn1 gene as one of the most strongly latitudinally differentiated SNPs in the genome of Drosophila melanogaster in the United States and Australia. Here we report the results of our transgenic analysis of this amino acid polymorphism. Genotypes carrying alternative Mnn1 alleles differed in multiple phenotypes in a direction generally consistent with phenotypic differences previously observed along latitudinal clines. These results support inferences from earlier population genomic work that this variant influences fitness, and support the idea that the alleles exhibiting clines may be likely to have pleiotropic effects that are correlated along the axes favored by natural selection.

Quantifying the relative contributions of the X chromosome, autosomes, and mitochondrial genome to local adaptation

During local adaptation with gene flow, some regions of the genome are inherently more responsive to selection than others. Recent theory predicts that X-linked genes should disproportionately contribute to local adaptation relative to other genomic regions, yet this prediction remains to be tested. We carried out a multigeneration crossing scheme, using two cline-end populations of Drosophila melanogaster, to estimate the relative contributions of the X chromosome, autosomes, and mitochondrial genome to divergence in four traits involved in local adaptation (wing size, resistance to heat, desiccation, and starvation stresses). We found that the mitochondrial genome and autosomes contributed significantly to clinal divergence in three of the four traits. In contrast, the X made no significant contribution to divergence in these traits. Given the small size of the mitochondrial genome, our results indicate that it plays a surprisingly large role in clinal adaptation. In contrast, the X, which represents roughly 20% of the Drosophila genome, contributes negligibly-a pattern that conflicts with theoretical predictions. These patterns reinforce recent work implying a central role of mitochondria in climatic adaptation, and suggest that different genomic regions may play fundamentally different roles in processes of divergence with gene flow.

Genome-wide patterns of local adaptation in Western European Drosophila melanogaster natural populations

Signatures of spatially varying selection have been investigated both at the genomic and transcriptomic level in several organisms. In Drosophila melanogaster, the majority of these studies have analyzed North American and Australian populations, leading to the identification of several loci and traits under selection. However, several studies based mainly in North American populations showed evidence of admixture that likely contributed to the observed population differentiation patterns. Thus, disentangling demography from selection might be challenging when analyzing these populations. European populations could help identify loci under spatially varying selection provided that no recent admixture from African populations would have occurred. In this work, we individually sequence the genome of 42 European strains collected in populations from contrasting environments: Stockholm (Sweden) and Castellana Grotte (Southern Italy). We found low levels of population structure and no evidence of recent African admixture in these two populations. We thus look for patterns of spatially varying selection affecting individual genes and gene sets. Besides single nucleotide polymorphisms, we also investigated the role of transposable elements in local adaptation. We concluded that European populations are a good dataset to identify candidate loci under spatially varying selection. The analysis of the two populations sequenced in this work in the context of all the available D. melanogaster data allowed us to pinpoint genes and biological processes likely to be relevant for local adaptation. Identifying and analyzing populations with low levels of population structure and admixture should help to disentangle selective from non-selective forces underlying patterns of population differentiation in other species as well.

A Large Panel of Drosophila simulans Reveals an Abundance of Common Variants

The rapidly expanding availability of large NGS data sets provides an opportunity to investigate population genetics at an unprecedented scale. Drosophila simulans is the sister species of the model organism Drosophila melanogaster, and is often presumed to share similar demographic history. However, previous population genetic and ecological work suggests very different signatures of selection and demography. Here, we sequence a new panel of 170 inbred genotypes of a North American population of D. simulans, a valuable complement to the DGRP and other D. melanogaster panels. We find some unexpected signatures of demography, in the form of excess intermediate frequency polymorphisms. Simulations suggest that this is possibly due to a recent population contraction and selection. We examine the outliers in the D. simulans genome determined by a haplotype test to attempt to parse the contribution of demography and selection to the patterns observed in this population. Untangling the relative contribution of demography and selection to genomic patterns of variation is challenging, however, it is clear that although D. melanogaster was thought to share demographic history with D. simulans different forces are at work in shaping genomic variation in this population of D. simulans.

Parallel clinal variation in the mid-day siesta of Drosophila melanogaster implicates continent-specific targets of natural selection

Similar to many diurnal animals, Drosophila melanogaster exhibits a mid-day siesta that is more robust as ambient temperature rises, an adaptive response aimed at minimizing exposure to heat. Mid-day siesta levels are partly regulated by the thermosensitive splicing of a small intron (termed dmpi8) found in the 3’ untranslated region (UTR) of the circadian clock gene period (per). Using the well-studied D. melanogaster latitudinal cline along the eastern coast of Australia, we show that flies from temperate populations sleep less during the day compared to those from tropical regions. We identified combinations of four single nucleotide polymorphisms (SNPs) in the 3’ UTR of per that yield several different haplotypes. The two most abundant of these haplotypes exhibit a reciprocal tropical-temperate distribution in relative frequency. Intriguingly, transgenic flies with the major tropical isoform manifest increased daytime sleep and reduced dmpi8 splicing compared to those carrying the temperate variant. Our results strongly suggest that for a major portion of D. melanogaster in Australia, thermal adaptation of daily sleep behavior included spatially varying selection on ancestrally derived polymorphisms in the per 3’ UTR that differentially control dmpi8 splicing efficiency. Prior work showed that African flies from high altitudes manifest reduced mid-day siesta levels, indicative of parallel latitudinal and altitudinal adaptation across continents. However, geographical variation in per 3’ UTR haplotypes was not observed for African flies, providing a compelling case for inter-continental variation in factors targeted by natural selection in attaining a parallel adaptation. We propose that the ability to calibrate mid-day siesta levels to better match local temperature ranges is a key adaptation contributing to the successful colonization of D. melanogaster beyond its ancestral range in the lowlands of Sub-Saharan Africa.

In warm climates many animals, including humans, exhibit a mid-day siesta, almost certainly a behavior meant to minimize the harm from prolonged exposure to the hot mid-day sun. But what about animals that adapted to cooler more temperate climates, might they have a less pronounced siesta? Indeed, we show that in the common fruit fly, Drosophila melanogaster, those from temperate regions in Australia exhibit less mid-day siesta compared to their tropical counterparts. Prior work showed that mid-day sleep levels are partially regulated by a ‘clock’ gene called period (per), which controls the timing of wake-sleep cycles in addition to other daily rhythms. We identified several DNA differences in the per gene that show geographical variation and contribute to the daytime sleep differences in flies from tropical and temperate regions via a mechanism that involves how well a temperature-sensitive intron in per is removed. A similar reduction in mid-day sleep was previously observed in African flies that adapted to the cooler temperatures found at high altitudes. Together, our findings provide a rare example where latitude and altitude lead to a similar behavioral adaptation to temperature. Moreover, the results suggest inter-continental differences in the evolutionary solutions used to attain the same thermal adaptation to cooler climates.

No detectable effect of Wolbachia wMel on the prevalence and abundance of the RNA virome of Drosophila melanogaster

Wolbachia is an endosymbiotic bacterium that can block viral infections in arthropods, generating interest in its potential to control the spread of mosquito-borne disease. Drosophila melanogaster is model organism for Wolbachia infection, and the wMel strain of Wolbachia can improve host survival following viral infection. However, it is unclear whether wMel induces anti-viral blocking against the broader native virome of D. melanogaster, or whether the major effect of Wolbachia is a reduction in viral abundance rather than viral clearance. We examined the effect of Wolbachia on viral abundance by comparing the total transcriptome of wMel-positive and wMel-negative D. melanogaster populations sampled from six locations in Australia. In addition, we examined the impact of wMel on individual flies by obtaining transcriptome data from 20 wMel-positive and 20 wMel-negative D. melanogaster from the location (Melbourne) with highest density of wMel. These data revealed high viral abundance in both Wolbachia-positive and -negative populations and individuals. Notably, none of the viral species identified, representing RNA viruses from at least nine families/floating genera, showed evidence of protection by wMel. Although the viral loads of picorna-like viruses are reduced by wMel under experimental conditions, we observed no such effect here. These data show that D. melanogaster can harbour abundant RNA viruses regardless of its Wolbachia status and imply that the interaction between Wolbachia and viruses in nature is more complex than simple blocking.

Enhancing Ebony? Common Associations With a cis-Regulatory Haplotype for Drosophila melanogaster Thoracic Pigmentation in a Japanese Population and Australian Populations

The molecular underpinnings of pigmentation diversity in Drosophila have recently emerged as a model for understanding how the evolution of different cis-regulatory variants results in common adaptive phenotypes within species. We compared sequence variation in a 5' regulatory region harboring a modular enhancer containing a ∼0.7-kb core element contributing to abdominal melanisation in African, and a ∼0.5-kb core element contributing to thoracic pigmentation in D. melanogaster from Japan, to tropical and temperate populations from eastern Australia previously shown to be divergent in thoracic pigmentation and ebony expression. The Australian populations exhibited strong association with the core enhancer polymorphism cluster in complete association with Dark and Light phenotypes from Iriomote, Japan. Moreover, the Iriomote Light and Dark core enhancer haplotypes are common to the Australian populations in the direction predicted by pigmentation phenotype. We also confirmed the Japanese patterns of linkage disequilibrium and association of the tropical inversion In(3R)Payne with the Light enhancer haplotype in the Australian tropical light population. A worldwide survey of the ∼0.5-kb ebony control region SNPs and haplotypes in a subset of the Drosophila Genome Nexus (DGN) populations suggest origins in the sub-Saharan ancestral region surrounding Zambia and subsequent invasion following colonization out of Africa. A previous study demonstrated complex within and between population genetic architecture for abdominal pigmentation which is also correlated with thoracic pigmentation in melanized DGN sub-Saharan populations; however, the ∼0.5-kb ebony control region was not associated and both haplotypes are common even in the most intensely pigmented D. melanogaster from high altitude Ethiopia. In the Australian populations, the strong phenotypic association with the enhancer SNPs and haplotypes that at least partly regulates ebony expression in the Iriomote population, our previous work demonstrating opposing clines for thoracic pigmentation and ebony expression, where the expression cline parallels the In(3R)Payne cline, and the concerted evolution of pigmentation intensity and ebony expression under rapid experimental evolution, all point to a common adaptive evolutionary pathway in distinct populations.