Parallel Evolution of Cold Tolerance within Drosophila melanogaster

Courtship song differs between African and European populations of Drosophila melanogaster and involves a strong effect locus

The courtship song of Drosophila melanogaster has long served as an excellent model system for studies of animal communication and differences in courtship song have been demonstrated among populations and between species. Here, we report that flies of African and European origin, which diverged approximately 13,000 years ago, show significant genetic differentiation in the use of slow versus fast pulse song. Using a combination of quantitative trait mapping and population genetic analysis we detected a single strong QTL underlying this trait and we identified candidate genes that may contribute to the evolution of this trait. Song trait variation between parental strains of our recombinant inbred panel enabled detection of genomic intervals associated with six additional song traits, some of which include known courtship-related genes. These findings improve the prospects for further genetic insights into the evolution of reproductive behavior and the biology underlying courtship song.

Recombinant inbred line panels inform the genetic architecture and interactions of adaptive traits in Drosophila melanogaster

The distribution of allelic effects on traits, along with their gene-by-gene and gene-by-environment interactions, contributes to the phenotypes available for selection and the trajectories of adaptive variants. Nonetheless, uncertainty persists regarding the effect sizes underlying adaptations and the importance of genetic interactions. Herein, we aimed to investigate the genetic architecture and the epistatic and environmental interactions involving loci that contribute to multiple adaptive traits using two new panels of Drosophila melanogaster recombinant inbred lines (RILs). To better fit our data, we re-implemented functions from R/qtl (Broman et al. 2003) using additive genetic models. We found 14 quantitative trait loci (QTL) underlying melanism, wing size, song pattern, and ethanol resistance. By combining our mapping results with population genetic statistics, we identified potential new genes related to these traits. None of the detected QTLs showed clear evidence of epistasis, and our power analysis indicated that we should have seen at least one significant interaction if sign epistasis or strong positive epistasis played a pervasive role in trait evolution. In contrast, we did find roles for gene-by-environment interactions involving pigmentation traits. Overall, our data suggest that the genetic architecture of adaptive traits often involves alleles of detectable effect, that strong epistasis does not always play a role in adaptation, and that environmental interactions can modulate the effect size of adaptive alleles.

Inversions Can Accumulate Balanced Sexual Antagonism: Evidence from Simulations and Drosophila Experiments

Chromosomal inversion polymorphisms can be common, but the causes of their persistence are often unclear. We propose a model for the maintenance of inversion polymorphism, which requires that some variants contribute antagonistically to two phenotypes, one of which has negative frequency-dependent fitness. These conditions yield a form of frequency-dependent disruptive selection, favoring two predominant haplotypes segregating alleles that favor opposing antagonistic phenotypes. An inversion associated with one haplotype can reduce the fitness load incurred by generating recombinant offspring, reinforcing its linkage to the haplotype and enabling both haplotypes to accumulate more antagonistic variants than expected otherwise. We develop and apply a forward simulator to examine these dynamics under a tradeoff between survival and male display. These simulations indeed generate inversion-associated haplotypes with opposing sex-specific fitness effects. Antagonism strengthens with time, and can ultimately yield karyotypes at surprisingly predictable frequencies, with striking genotype frequency differences between sexes and between developmental stages. To test whether this model may contribute to well-studied yet enigmatic inversion polymorphisms in Drosophila melanogaster, we track inversion frequencies in laboratory crosses to test whether they influence male reproductive success or survival. We find that two of the four tested inversions show significant evidence for the tradeoff examined, with ln(3R)K favoring survival and ln(3L)Ok favoring male reproduction. In line with the apparent sex-specific fitness effects implied for both of those inversions, ln(3L)Ok was also found to be less costly to the viability and/or longevity of males than females, whereas ln(3R)K was more beneficial to female survival. Based on this work, we expect that balancing selection on antagonistically pleiotropic traits may provide a significant and underappreciated contribution to the maintenance of natural inversion polymorphism.

Genomic Diversity Illuminates the Environmental Adaptation of Drosophila suzukii

Biological invasions carry substantial practical and scientific importance and represent natural evolutionary experiments on contemporary timescales. Here, we investigated genomic diversity and environmental adaptation of the crop pest Drosophila suzukii using whole-genome sequencing data and environmental metadata for 29 population samples from its native and invasive range. Through a multifaceted analysis of this population genomic data, we increase our understanding of the D. suzukii genome, its diversity and its evolution, and we identify an appropriate genotype-environment association pipeline for our dataset. Using this approach, we detect genetic signals of local adaptation associated with nine distinct environmental factors related to altitude, wind speed, precipitation, temperature, and human land use. We uncover unique functional signatures for each environmental variable, such as the prevalence of cuticular genes associated with annual precipitation. We also infer biological commonalities in the adaptation to diverse selective pressures, particularly in terms of the apparent contribution of nervous system evolution to enriched processes (ranging from neuron development to circadian behavior) and to top genes associated with all nine environmental variables. Our findings therefore depict a finer-scale adaptive landscape underlying the rapid invasion success of this agronomically important species.

Heat tolerance, oxidative stress response tuning and robust gene activation in early-stage Drosophila melanogaster embryos

In organisms with complex life cycles, life stages that are most susceptible to environmental stress may determine species persistence in the face of climate change. Early embryos of Drosophila melanogaster are particularly sensitive to acute heat stress, yet tropical embryos have higher heat tolerance than temperate embryos, suggesting adaptive variation in embryonic heat tolerance. We compared transcriptomic responses to heat stress among tropical and temperate embryos to elucidate the gene regulatory basis of divergence in embryonic heat tolerance. The transcriptomes of tropical and temperate embryos differed in both constitutive and heat-stress-induced responses of the expression of relatively few genes, including genes involved in oxidative stress. Most of the transcriptomic response to heat stress was shared among all embryos. Embryos shifted the expression of thousands of genes, including increases in the expression of heat shock genes, suggesting robust zygotic gene activation and demonstrating that, contrary to previous reports, early embryos are not transcriptionally silent. The involvement of oxidative stress genes corroborates recent reports on the critical role of redox homeostasis in coordinating developmental transitions. By characterizing adaptive variation in the transcriptomic basis of embryonic heat tolerance, this study is a novel contribution to the literature on developmental physiology and developmental genetics.

Copy number variation introduced by a massive mobile element facilitates global thermal adaptation in a fungal wheat pathogen

Copy number variation (CNV) can drive rapid evolution in changing environments. In microbial pathogens, such adaptation is a key factor underpinning epidemics and colonization of new niches. However, the genomic determinants of such adaptation remain poorly understood. Here, we systematically investigate CNVs in a large genome sequencing dataset spanning a worldwide collection of 1104 genomes from the major wheat pathogen Zymoseptoria tritici. We found overall strong purifying selection acting on most CNVs. Genomic defense mechanisms likely accelerated gene loss over episodes of continental colonization. Local adaptation along climatic gradients was likely facilitated by CNVs affecting secondary metabolite production and gene loss in general. One of the strongest loci for climatic adaptation is a highly conserved gene of the NAD-dependent Sirtuin family. The Sirtuin CNV locus localizes to an ~68-kb Starship mobile element unique to the species carrying genes highly expressed during plant infection. The element has likely lost the ability to transpose, demonstrating how the ongoing domestication of cargo-carrying selfish elements can contribute to selectable variation within populations. Our work highlights how standing variation in gene copy numbers at the global scale can be a major factor driving climatic and metabolic adaptation in microbial species.

The Precision and Power of Population Branch Statistics in Identifying the Genomic Signatures of Local Adaptation

Population branch statistics, which estimate the branch lengths of focal populations with respect to two outgroups, have been used as an alternative to FST-based genome-wide scans for identifying loci associated with local selective sweeps. In addition to the original population branch statistic (PBS), there are subsequently proposed branch rescalings: normalized population branch statistic (PBSn1), which adjusts focal branch length with respect to outgroup branch lengths at the same locus, and population branch excess (PBE), which also incorporates median branch lengths at other loci. PBSn1 and PBE have been proposed to be less sensitive to allele frequency divergence generated by background selection or geographically ubiquitous positive selection rather than local selective sweeps. However, the accuracy and statistical power of branch statistics have not been systematically assessed. To do so, we simulate genomes in representative large and small populations with varying proportions of sites evolving under genetic drift or background selection (approximated using variable Ne), local selective sweeps, and geographically parallel selective sweeps. We then assess the probability that local selective sweep loci are correctly identified as outliers by FST and by each of the branch statistics. We find that branch statistics consistently outperform FST at identifying local sweeps. When background selection and/or parallel sweeps are introduced, PBSn1 and especially PBE correctly identify local sweeps among their top outliers at a higher frequency than PBS. These results validate the greater specificity of rescaled branch statistics such as PBE to detect population-specific positive selection, supporting their use in genomic studies focused on local adaptation.

Genomic Diversity Illuminates the Environmental Adaptation of Drosophila suzukii

Biological invasions carry substantial practical and scientific importance, and represent natural evolutionary experiments on contemporary timescales. Here, we investigated genomic diversity and environmental adaptation of the crop pest Drosophila suzukii using whole-genome sequencing data and environmental metadata for 29 population samples from its native and invasive range. Through a multifaceted analysis of this population genomic data, we increase our understanding of the D. suzukii genome, its diversity and its evolution, and we identify an appropriate genotype-environment association pipeline for our data set. Using this approach, we detect genetic signals of local adaptation associated with nine distinct environmental factors related to altitude, wind speed, precipitation, temperature, and human land use. We uncover unique functional signatures for each environmental variable, such as a prevalence of cuticular genes associated with annual precipitation. We also infer biological commonalities in the adaptation to diverse selective pressures, particularly in terms of the apparent contribution of nervous system evolution to enriched processes (ranging from neuron development to circadian behavior) and to top genes associated with all nine environmental variables. Our findings therefore depict a finer-scale adaptive landscape underlying the rapid invasion success of this agronomically important species.

Hybrid breakdown in male reproduction between recently diverged Drosophila melanogaster populations has a complex and variable genetic architecture

Secondary contact between formerly isolated populations may result in hybrid breakdown, in which untested allelic combinations in hybrids are maladaptive and limit genetic exchange. Studying early-stage reproductive isolation may yield key insights into the genetic architectures and evolutionary forces underlying the first steps toward speciation. Here, we leverage the recent worldwide expansion of Drosophila melanogaster to test for hybrid breakdown between populations that diverged within the last 13,000 years. We found clear evidence for hybrid breakdown in male reproduction, but not female reproduction or viability, supporting the prediction that hybrid breakdown affects the heterogametic sex first. The frequency of non-reproducing F2 males varied among different crosses involving the same southern African and European populations, as did the qualitative effect of cross direction, implying a genetically variable basis of hybrid breakdown and a role for uniparentally inherited factors. The levels of breakdown observed in F2 males were not recapitulated in backcrossed individuals, consistent with the existence of incompatibilities with at least three partners. Thus, some of the very first steps toward reproductive isolation could involve incompatibilities with complex and variable genetic architectures. Collectively, our findings emphasize this system's potential for future studies on the genetic and organismal basis of early-stage reproductive isolation.

Ethanol resistance in Drosophila melanogaster has increased in parallel cold-adapted populations and shows a variable genetic architecture within and between populations

Understanding the genetic properties of adaptive trait evolution is a fundamental crux of biological inquiry that links molecular processes to biological diversity. Important uncertainties persist regarding the genetic predictability of adaptive trait change, the role of standing variation, and whether adaptation tends to result in the fixation of favored variants. Here, we use the recurrent evolution of enhanced ethanol resistance in Drosophila melanogaster during this species' worldwide expansion as a promising system to add to our understanding of the genetics of adaptation. We find that elevated ethanol resistance has evolved at least three times in different cooler regions of the species' modern range-not only at high latitude but also in two African high-altitude regions. Applying a bulk segregant mapping framework, we find that the genetic architecture of ethanol resistance evolution differs substantially not only between our three resistant populations, but also between two crosses involving the same European population. We then apply population genetic scans for local adaptation within our quantitative trait locus regions, and we find potential contributions of genes with annotated roles in spindle localization, membrane composition, sterol and alcohol metabolism, and other processes. We also apply simulation-based analyses that confirm the variable genetic basis of ethanol resistance and hint at a moderately polygenic architecture. However, these simulations indicate that larger-scale studies will be needed to more clearly quantify the genetic architecture of adaptive evolution and to firmly connect trait evolution to specific causative loci.

Differences in temperature-sensitive expression of PcG-regulated genes among natural populations of Drosophila melanogaster

Environmental temperature can affect chromatin-based gene regulation, in particular in ectotherms such as insects. Genes regulated by the Polycomb group (PcG) vary in their transcriptional output in response to changes in temperature. Expression of PcG-regulated genes typically increases with decreasing temperatures. Here, we examined variations in temperature-sensitive expression of PcG target genes in natural populations from different climates of Drosophila melanogaster, and differences thereof across different fly stages and tissues. Temperature-induced expression plasticity was found to be stage- and sex-specific with differences in the specificity between the examined PcG target genes. Some tissues and stages, however, showed a higher number of PcG target genes with temperature-sensitive expression than others. Overall, we found higher levels of temperature-induced expression plasticity in African tropical flies from the ancestral species range than in flies from temperate Europe. We also observed differences between temperate flies, however, with more reduction of expression plasticity in warm-temperate than in cold-temperate populations. Although in general, temperature-sensitive expression appeared to be detrimental in temperate climates, there were also cases in which plasticity was increased in temperate flies, as well as no changes in expression plasticity between flies from different climates.

The effects of genetic drift and genomic selection on differentiation and local adaptation of the introduced populations of Aedes albopictus in southern Russia

Background

Asian tiger mosquito Aedes albopictus is an arbovirus vector that has spread from its native habitation areal in Southeast Asia throughout North and South Americas, Europe, and Africa. Ae. albopictus was first detected in the Southern Federal District of the Russian Federation in the subtropical town of Sochi in 2011. In subsequent years, this species has been described in the continental areas with more severe climate and lower winter temperatures.

Methods

Genomic analysis of pooled Ae. albopictus samples collected in the mosquito populations in the coastal and continental regions of the Krasnodar Krai was conducted to look for the genetic changes associated with the spread and potential cold adaptation in Ae. albopictus.

Results

The results of the phylogenetic analysis based on mitochondrial genomes corresponded well with the hypothesis that Ae. albopictus haplotype A1a2a1 was introduced into the region from a single source. Population analysis revealed the role of dispersal and genetic drift in the local adaptation of the Asian tiger mosquito. The absence of shared haplotypes between the samples and high fixation indices suggest that gene flow between samples was heavily restricted. Mitochondrial and genomic differentiation together with different distances between dispersal routes, natural and anthropogenic barriers and local effective population size reduction could lead to difficulties in local climatic adaptations due to reduced selection effectiveness. We have found genomic regions with selective sweep patterns which can be considered as having been affected by recent selection events. The genes located in these regions participate in neural protection, lipid conservation, and cuticle formation during diapause. These processes were shown to be important for cold adaptation in the previous transcriptomic and proteomic studies. However, the population history and relatively low coverage obtained in the present article could have negatively affect sweep detection.

Chill coma onset and recovery fail to reveal true variation in thermal performance among populations of Drosophila melanogaster

Species from colder climates tend to be more chill tolerant regardless of the chill tolerance trait measured, but for Drosophila melanogaster, population-level differences in chill tolerance among populations are not always found when a single trait is measured in the laboratory. We measured chill coma onset temperature, chill coma recovery time, and survival after chronic cold exposure in replicate lines derived from multiple paired African and European D. melanogaster populations. The populations in our study were previously found to differ in chronic cold survival ability, which is believed to have evolved independently in each population pair; however, they did not differ in chill coma onset temperature and chill coma recovery time in a manner that reflected their geographic origins, even though these traits are known to vary with origin latitude among Drosophila species and are among the most common metrics of thermal tolerance in insects. While it is common practice to measure only one chill tolerance trait when comparing chill tolerance among insect populations, our results emphasise the importance of measuring more than one thermal tolerance trait to minimize the risk of missing real adaptive variation in insect thermal tolerance.

False and true positives in arthropod thermal adaptation candidate gene lists

Genome-wide studies are prone to false positives due to inherently low priors and statistical power. One approach to ameliorate this problem is to seek validation of reported candidate genes across independent studies: genes with repeatedly discovered effects are less likely to be false positives. Inversely, genes reported only as many times as expected by chance alone, while possibly representing novel discoveries, are also more likely to be false positives. We show that, across over 30 genome-wide studies that reported Drosophila and Daphnia genes with possible roles in thermal adaptation, the combined lists of candidate genes and orthologous groups are rapidly approaching the total number of genes and orthologous groups in the respective genomes. This is consistent with the expectation of high frequency of false positives. The majority of these spurious candidates have been identified by one or a few studies, as expected by chance alone. In contrast, a noticeable minority of genes have been identified by numerous studies with the probabilities of such discoveries occurring by chance alone being exceedingly small. For this subset of genes, different studies are in agreement with each other despite differences in the ecological settings, genomic tools and methodology, and reporting thresholds. We provide a reference set of presumed true positives among Drosophila candidate genes and orthologous groups involved in response to changes in temperature, suitable for cross-validation purposes. Despite this approach being prone to false negatives, this list of presumed true positives includes several hundred genes, consistent with the "omnigenic" concept of genetic architecture of complex traits.

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.

Recurrent Collection of Drosophila melanogaster from Wild African Environments and Genomic Insights into Species History

A long-standing enigma concerns the geographic and ecological origins of the intensively studied vinegar fly, Drosophila melanogaster. This globally distributed human commensal is thought to originate from sub-Saharan Africa, yet until recently, it had never been reported from undisturbed wilderness environments that could reflect its precommensal niche. Here, we document the collection of 288 D. melanogaster individuals from multiple African wilderness areas in Zambia, Zimbabwe, and Namibia. The presence of D. melanogaster in these remote woodland environments is consistent with an ancestral range in southern-central Africa, as opposed to equatorial regions. After sequencing the genomes of 17 wilderness-collected flies collected from Kafue National Park in Zambia, we found reduced genetic diversity relative to town populations, elevated chromosomal inversion frequencies, and strong differences at specific genes including known insecticide targets. Combining these genomes with existing data, we probed the history of this species' geographic expansion. Demographic estimates indicated that expansion from southern-central Africa began ∼10,000 years ago, with a Saharan crossing soon after, but expansion from the Middle East into Europe did not begin until roughly 1,400 years ago. This improved model of demographic history will provide an important resource for future evolutionary and genomic studies of this key model organism. Our findings add context to the history of D. melanogaster, while opening the door for future studies on the biological basis of adaptation to human environments.

Decreased Temperature Sensitivity of Vestigial Gene Expression in Temperate Populations of Drosophila melanogaster

Drosophila melanogaster recently spread from its tropical origin in Africa and became a cosmopolitan species that has adapted to a wide range of different thermal environments, including temperate climates. An important limiting factor of temperate climates has probably been their low and varying temperatures. The transcriptional output of genes can vary across temperatures, which might have been detrimental while settling in temperate environments. The reduction of temperature-sensitive expression of functionally important genes to ensure consistent levels of gene expression might have been relevant while adapting to such environments. In this study, we focus on the gene vestigial (vg) whose product is a key factor in wing development. We provide evidence that temperature-sensitivity of vg has been buffered in populations from temperate climates. We investigated temperature-sensitivity of vg gene expression in six natural populations, including four temperate populations (three from Europe and one from high-altitude Africa), and two tropical populations from the ancestral species range. All temperate populations exhibited a lower degree of temperature-induced expression plasticity than the tropical populations.

How chromosomal rearrangements shape adaptation and speciation: Case studies in Drosophila pseudoobscura and its sibling species Drosophila persimilis

The gene arrangements of Drosophila have played a prominent role in the history of evolutionary biology from the original quantification of genetic diversity to current studies of the mechanisms for the origin and establishment of new inversion mutations within populations and their subsequent fixation between species supporting reproductive barriers. This review examines the genetic causes and consequences of inversions as recombination suppressors and the role that recombination suppression plays in establishing inversions in populations as they are involved in adaptation within heterogeneous environments. This often results in the formation of clines of gene arrangement frequencies among populations. Recombination suppression leads to the differentiation of the gene arrangements which may accelerate the accumulation of fixed genetic differences among populations. If these fixed mutations cause incompatibilities, then inversions pose important reproductive barriers between species. This review uses the evolution of inversions in Drosophila pseudoobscura and D. persimilis as a case study for how inversions originate, establish and contribute to the evolution of reproductive isolation.

Stress response, behavior, and development are shaped by transposable element-induced mutations in Drosophila

Most of the current knowledge on the genetic basis of adaptive evolution is based on the analysis of single nucleotide polymorphisms (SNPs). Despite increasing evidence for their causal role, the contribution of structural variants to adaptive evolution remains largely unexplored. In this work, we analyzed the population frequencies of 1,615 Transposable Element (TE) insertions annotated in the reference genome of Drosophila melanogaster, in 91 samples from 60 worldwide natural populations. We identified a set of 300 polymorphic TEs that are present at high population frequencies, and located in genomic regions with high recombination rate, where the efficiency of natural selection is high. The age and the length of these 300 TEs are consistent with relatively young and long insertions reaching high frequencies due to the action of positive selection. Besides, we identified a set of 21 fixed TEs also likely to be adaptive. Indeed, we, and others, found evidence of selection for 84 of these reference TE insertions. The analysis of the genes located nearby these 84 candidate adaptive insertions suggested that the functional response to selection is related with the GO categories of response to stimulus, behavior, and development. We further showed that a subset of the candidate adaptive TEs affects expression of nearby genes, and five of them have already been linked to an ecologically relevant phenotypic effect. Our results provide a more complete understanding of the genetic variation and the fitness-related traits relevant for adaptive evolution. Similar studies should help uncover the importance of TE-induced adaptive mutations in other species as well.

Muller "Elements" in Drosophila: How the Search for the Genetic Basis for Speciation Led to the Birth of Comparative Genomics

The concept of synteny, or conservation of genes on the same chromosome, traces its origins to the early days of Drosophila genetics. This discovery emerged from comparisons of linkage maps from different species of Drosophila with the goal of understanding the process of speciation. H. J. Muller published a landmark article entitled Bearings of the "Drosophila" work on systematics, where he synthesized genetic and physical map data and proposed a model of speciation and chromosomal gene content conservation. These models have withstood the test of time with the advent of molecular genetic analysis from protein to genome level variation. Muller's ideas provide a framework to begin to answer questions about the evolutionary forces that shape the structure of the genome.