Inference of chromosomal inversion dynamics from Pool‐Seq data in natural and laboratory populations of Drosophila melanogaster

Sequencing of pools of individuals (Pool‐Seq) represents a reliable and cost‐effective approach for estimating genome‐wide SNP and transposable element insertion frequencies. However, Pool‐Seq does not provide direct information on haplotypes so that, for example, obtaining inversion frequencies has not been possible until now. Here, we have developed a new set of diagnostic marker SNPs for seven cosmopolitan inversions in Drosophila melanogaster that can be used to infer inversion frequencies from Pool‐Seq data. We applied our novel marker set to Pool‐Seq data from an experimental evolution study and from North American and Australian latitudinal clines. In the experimental evolution data, we find evidence that positive selection has driven the frequencies of In(3R)C and In(3R)Mo to increase over time. In the clinal data, we confirm the existence of frequency clines for In(2L)t, In(3L)P and In(3R)Payne in both North America and Australia and detect a previously unknown latitudinal cline for In(3R)Mo in North America. The inversion markers developed here provide a versatile and robust tool for characterizing inversion frequencies and their dynamics in Pool‐Seq data from diverse D. melanogaster populations.     

Monophyly of Wolbachia pipientis genomes within Drosophila melanogaster: geographic structuring,titre variation and host effects across five populations

Wolbachia pipientis is one of the most widely studied endosymbionts today, yet we know little about its short‐term adaptation and evolution. Here, using a set of 91 inbred Drosophila melanogaster lines from five populations, we explore patterns of diversity and recent evolution in the Wolbachia strain wMel. Within the D. melanogaster lines, we identify six major mitochondrial clades and four wMel clades. Concordant with past studies, the Wolbachia haplotypes contain an overall low level of nucleotide diversity, yet they still display geographic structuring. Using Bayesian analysis informed with demographic estimates of colonization times, we estimate that all extant D. melanogaster mitochondrial haplotypes coalesce to a Wolbachia‐infected ancestor approximately 2200 years ago. Finally, we measure wMel titre within the infected flies and find that titre varies across populations, an effect attributable to host genetic factors. This demonstration of local phenotypic divergence suggests that intraspecific host genetic variation plays a key role in shaping this model symbiotic system.     

Extensive paternal mtDNA leakage in natural populations of Drosophila melanogaster

Strict maternal inheritance is considered a hallmark of animal mtDNA. Although recent reports suggest that paternal leakage occurs in a broad range of species, it is still considered an exceptionally rare event. To evaluate the impact of paternal leakage on the evolution of mtDNA, it is essential to reliably estimate the frequency of paternal leakage in natural populations. Using allele‐specific real‐time quantitative PCR (RT‐qPCR), we show that heteroplasmy is common in natural populations with at least 14% of the individuals carrying multiple mitochondrial haplotypes. However, the average frequency of the minor mtDNA haplotype is low (0.8%), which suggests that this pervasive heteroplasmy has not been noticed before due to a lack of power in sequencing surveys. Based on the distribution of mtDNA haplotypes in the offspring of heteroplasmic mothers, we found no evidence for strong selection against one of the haplotypes. We estimated that the rate of paternal leakage is 6% and that at least 100 generations are required for complete sorting of mtDNA haplotypes. Despite the high proportion of heteroplasmic individuals in natural populations, we found no evidence for recombination between mtDNA molecules, suggesting that either recombination is rare or recombinant haplotypes are counter‐selected. Our results indicate that evolutionary studies using mtDNA as a marker might be biased by paternal leakage in this species.     

Evidence for complex selection on four‐fold degenerate sites in Drosophila melanogaster

We considered genome‐wide four‐fold degenerate sites from an African Drosophila melanogaster population and compared them to short introns. To include divergence and to polarize the data, we used its close relatives Drosophila simulans, Drosophila sechellia, Drosophila erecta and Drosophila yakuba as outgroups. In D. melanogaster, the GC content at four‐fold degenerate sites is higher than in short introns; compared to its relatives, more AT than GC is fixed. The former has been explained by codon usage bias (CUB) favouring GC; the latter by decreased intensity of directional selection or by increased mutation bias towards AT. With a biallelic equilibrium model, evidence for directional selection comes mostly from the GC‐rich ancestral base composition. Together with a slight mutation bias, it leads to an asymmetry of the unpolarized allele frequency spectrum, from which directional selection is inferred. Using a quasi‐equilibrium model and polarized spectra, however, only purifying and no directional selection is detected. Furthermore, polarized spectra are proportional to those of the presumably unselected short introns. As we have no evidence for a decrease in effective population size, relaxed CUB must be due to a reduction in the selection coefficient. Going beyond the biallelic model and considering all four bases, signs of directional selection are stronger. In contrast to short introns, complementary bases show strand specificity and allele frequency spectra depend on mutation directions. Hence, the traditional biallelic model to describe the evolution of four‐fold degenerate sites should be replaced by more complex models assuming only quasi‐equilibrium and accounting for all four bases.     

Phenotypic plasticity,but not adaptive tracking,underlies seasonal variation in post‐cold hardening freeze tolerance of Drosophila melanogaster

In temperate regions, an organism's ability to rapidly adapt to seasonally varying environments is essential for its survival. In response to seasonal changes in selection pressure caused by variation in temperature, humidity, and food availability, some organisms exhibit plastic changes in phenotype. In other cases, seasonal variation in selection pressure can rapidly increase the frequency of genotypes that offer survival or reproductive advantages under the current conditions. Little is known about the relative influences of plastic and genetic changes in short‐lived organisms experiencing seasonal environmental fluctuations. Cold hardening is a seasonally relevant plastic response in which exposure to cool, but nonlethal, temperatures significantly increases the organism's ability to later survive at freezing temperatures. In the present study, we demonstrate seasonal variation in cold hardening in Drosophila melanogaster and test the extent to which plasticity and adaptive tracking underlie that seasonal variation. We measured the post‐cold hardening freeze tolerance of flies from outdoor mesocosms over the summer, fall, and winter. We bred outdoor mesocosm‐caught flies for two generations in the laboratory and matched each outdoor cohort to an indoor control cohort of similar genetic background. We cold hardened all flies under controlled laboratory conditions and then measured their post‐cold hardening freeze tolerance. Comparing indoor and field‐caught flies and their laboratory‐reared G1 and G2 progeny allowed us to determine the roles of seasonal environmental plasticity, parental effects, and genetic changes on cold hardening. We also tested the relationship between cold hardening and other factors, including age, developmental density, food substrate, presence of antimicrobials, and supplementation with live yeast. We found strong plastic responses to a variety of field‐ and laboratory‐based environmental effects, but no evidence of seasonally varying parental or genetic effects on cold hardening. We therefore conclude that seasonal variation in post‐cold hardening freeze tolerance results from environmental influences and not genetic changes.     

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

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

Kin recognition and co‐operative foraging in Drosophila melanogaster larvae

A long‐standing goal for biologists and social scientists is to understand the factors that lead to the evolution and maintenance of co‐operative behaviour between conspecifics. To that end, the fruit fly, Drosophila melanogaster, is becoming an increasingly popular model species to study sociality; however, most of the research to date has focused on adult behaviours. In this study, we set out to examine group‐feeding behaviour by larvae and to determine whether the degree of relatedness between individuals mediates the expression co‐operation. In a series of assays, we manipulated the average degree of relatedness in groups of third‐instar larvae that were faced with resource scarcity, and measured the size, frequency and composition of feeding clusters, as well as the fitness benefits associated with co‐operation. Our results suggest that larval D. melanogaster are capable of kin recognition (something that has not been previously described in this species), as clusters were more numerous, larger and involved more larvae, when more closely related kin were present in the social environment. These findings are discussed in the context of the correlated fitness‐associated benefits of co‐operation, the potential mechanisms by which individuals may recognize kin, and how that kinship may play an important role in facilitating the manifestation of this co‐operative behaviour.     

Inbreeding reveals mode of past selection on male reproductive characters in Drosophila melanogaster

Directional dominance is a prerequisite of inbreeding depression. Directionality arises when selection drives alleles that increase fitness to fixation and eliminates dominant deleterious alleles, while deleterious recessives are hidden from it and maintained at low frequencies. Traits under directional selection (i.e., fitness traits) are expected to show directional dominance and therefore an increased susceptibility to inbreeding depression. In contrast, traits under stabilizing selection or weakly linked to fitness are predicted to exhibit little‐to‐no inbreeding depression. Here, we quantify the extent of inbreeding depression in a range of male reproductive characters and then infer the mode of past selection on them. The use of transgenic populations of Drosophila melanogaster with red or green fluorescent‐tagged sperm heads permitted in vivo discrimination of sperm from competing males and quantification of characteristics of ejaculate composition, performance, and fate. We found that male attractiveness (mating latency) and competitive fertilization success (P₂) both show some inbreeding depression, suggesting they may have been under directional selection, whereas sperm length showed no inbreeding depression suggesting a history of stabilizing selection. However, despite having measured several sperm quality and quantity traits, our data did not allow us to discern the mechanism underlying the lowered competitive fertilization success of inbred (f = 0.50) males.     

The multigenerational effects of water contamination and endocrine disrupting chemicals on the fitness of Drosophila melanogaster

Water pollution due to human activities produces sedimentation, excessive nutrients, and toxic chemicals, and this, in turn, has an effect on the normal endocrine functioning of living beings. Overall, water pollution may affect some components of the fitness of organisms (e.g., developmental time and fertility). Some toxic compounds found in polluted waters are known as endocrine disruptors (ED), and among these are nonhalogenated phenolic chemicals such as bisphenol A and nonylphenol. To evaluate the effect of nonhalogenated phenolic chemicals on the endocrine system, we subjected two generations (F0 and F1) of Drosophila melanogaster to different concentrations of ED. Specifically, treatments involved wastewater, which had the highest level of ED (bisphenol A and nonylphenol) and treated wastewater from a constructed Heliconia psittacorum wetland with horizontal subsurface water flow (He); the treated wastewater was the treatment with the lowest level of ED. We evaluated the development time from egg to pupa and from pupa to adult as well as fertility. The results show that for individuals exposed to treated wastewater, the developmental time from egg to pupae was shorter in individuals of the F1 generation than in the F0 generation. Additionally, the time from pupae to adult was longer for flies growing in the H. psittacorum treated wastewater. Furthermore, fertility was lower in the F1 generation than in the F0 generation. Although different concentrations of bisphenol A and nonylphenol had no significant effect on the components of fitness of D. melanogaster (developmental time and fertility), there was a trend across generations, likely as a result of selection imposed on the flies. It is possible that the flies developed different strategies to avoid the effects of the various environmental stressors.     

Experimental evolution of slowed cognitive aging in Drosophila melanogaster

Reproductive output and cognitive performance decline in parallel during aging, but it is unknown whether this reflects a shared genetic architecture or merely the declining force of natural selection acting independently on both traits. We used experimental evolution in Drosophila melanogaster to test for the presence of genetic variation for slowed cognitive aging, and assess its independence from that responsible for other traits’ decline with age. Replicate experimental populations experienced either joint selection on learning and reproduction at old age (Old + Learning), selection on late‐life reproduction alone (Old), or a standard two‐week culture regime (Young). Within 20 generations, the Old + Learning populations evolved a slower decline in learning with age than both the Old and Young populations, revealing genetic variation for cognitive aging. We found little evidence for a genetic correlation between cognitive and demographic aging: although the Old + Learning populations tended to show higher late‐life fecundity than Old populations, they did not live longer. Likewise, selection for late reproduction alone did not result in improved late‐life learning. Our results demonstrate that Drosophila harbor genetic variation for cognitive aging that is largely independent from genetic variation for demographic aging and suggest that these two aspects of aging may not necessarily follow the same trajectories.     

Reverse evolution of fitness in Drosophila melanogaster

Abstract The evolution of fitness is central to evolutionary theory, yet few experimental systems allow us to track its evolution in genetically and environmentally relevant contexts. Reverse evolution experiments allow the study of the evolutionary return to ancestral phenotypic states, including fitness. This in turn permits well‐defined tests for the dependence of adaptation on evolutionary history and environmental conditions. In the experiments described here, 20 populations of heterogeneous evolutionary histories were returned to their common ancestral environment for 50 generations, and were then compared with both their immediate differentiated ancestors and populations which had remained in the ancestral environment. One measure of fitness returned to ancestral levels to a greater extent than other characters did. The phenotypic effects of reverse evolution were also contingent on previous selective history. Moreover, convergence to the ancestral state was highly sensitive to environmental conditions. The phenotypic plasticity of fecundity, a character directly selected for, evolved during the experimental time frame. Reverse evolution appears to force multiple, diverged populations to converge on a common fitness state through different life‐history and genetic changes.     

A test for Y‐linked additive and epistatic effects on surviving bacterial infections in Drosophila melanogaster

Y‐ and W‐chromosomes offer a theoretically powerful way for sexual dimorphism to evolve. Consistent with this possibility, Drosophila melanogaster Y‐chromosomes can influence gene regulation throughout the genome; particularly immune‐related genes. In order for Y‐linked regulatory variation (YRV) to contribute to adaptive evolution it must be comprised of additive genetic variance, such that variable Ys induce consistent phenotypic effects within the local gene pool. We assessed the potential for Y‐chromosomes to adaptively shape gram‐negative and gram‐positive bacterial defence by introgressing Ys across multiple genetic haplotypes from the same population. We found no Y‐linked additive effects on immune phenotypes, suggesting a restricted role for the Y to facilitate dimorphic evolution. We did find, however, a large magnitude Y by background interaction that induced rank order reversals of Y‐effects across the backgrounds (i.e. sign epistasis). Thus, Y‐chromosome effects appeared consistent within backgrounds, but highly variable among backgrounds. This large sign epistatic effect could constrain monomorphic selection in both sexes, considering that autosomal alleles under selection must spend half of their time in a male background where relative fitness values are altered. If the pattern described here is consistent for other traits or within other XY (or ZW) systems, then YRV may represent a universal constraint to autosomal trait evolution.     

Geographical delimitation of a partial selective sweep in African Drosophila melanogaster

Positive selection leaves characteristic footprints on DNA variation but detecting such patterns is challenging as the age, the intensity and the mode of selection as well as demography and evolutionary parameters (mutation and recombination rates) all play roles and these are difficult to disentangle. We recorded nucleotide variation in a sample of isogenic chromosomes from a western African population of Drosophila melanogaster at a locus (Fbp2) for which a partial selective sweep had previously been reported. We compared this locus to four other genes from the same chromosomes and from a European and an East African population. Then, we assessed Fbp2 variation in a sample of 370 chromosomes covering a comprehensive geographic sampling of 16 African localities. The signature of selection was tested while accounting for the demographic history of the populations. We found a significant signal of selection in two West African localities including Ivory Coast. Variation at Fpb2 would thus represent a case of an ongoing selective sweep in the range of this species. A weaker, nonsignificant, signal of selection was, however, apparent in some other populations, thus leaving open several possibilities: (i) the selective sweep originated in Ivory Coast and has spread to the rest of the continent; (ii) several African populations report the signature of a selective event having occurred in an ancestral population; (iii) this genome region is subject to independent selective events in African populations; and (iv) A neutral scenario with population subdivision and local bottleneck cannot be fully excluded to explain the molecular patterns observed in some populations.     

UCP4C mediates uncoupled respiration in larvae of Drosophila melanogaster

Larvae of Drosophila melanogaster reared at 23°C and switched to 14°C for 1 h are 0.5°C warmer than the surrounding medium. In keeping with dissipation of energy, respiration of Drosophila melanogaster larvae cannot be decreased by the F‐ATPase inhibitor oligomycin or stimulated by protonophore. Silencing of Ucp4C conferred sensitivity of respiration to oligomycin and uncoupler, and prevented larva‐to‐adult progression at 15°C but not 23°C. Uncoupled respiration of larval mitochondria required palmitate, was dependent on Ucp4C and was inhibited by guanosine diphosphate. UCP4C is required for development through the prepupal stages at low temperatures and may be an uncoupling protein.     

Dietary choices are influenced by genotype,mating status,and sex in Drosophila melanogaster

Mating causes many changes in physiology, behavior, and gene expression in a wide range of organisms. These changes are predicted to be sex specific, influenced by the divergent reproductive roles of the sexes. In female insects, mating is associated with an increase in egg production which requires high levels of nutritional input with direct consequences for the physiological needs of individual females. Consequently, females alter their nutritional acquisition in line with the physiological demands imposed by mating. Although much is known about the female mating‐induced nutritional response, far less is known about changes in males. In addition, it is unknown whether variation between genotypes translates into variation in dietary behavioral responses. Here we examine mating‐induced shifts in male and female dietary preferences across genotypes of Drosophila melanogaster. We find sex‐ and genotype‐specific effects on both the quantity and quality of the chosen diet. These results contribute to our understanding of sex‐specific metabolism and reveal genotypic variation that influences responses to physiological demands.     

Nutritional geometry and fitness consequences in Drosophila suzukii,the Spotted‐Wing Drosophila

Since its arrival to North America less than a decade ago, the invasive Spotted‐Wing Drosophila (Drosophila suzukii) has inflicted substantial economic losses on soft fruit agriculture due to its ability to oviposit into ripening fruits. More effective management approaches for this species are needed, but little is known about the factors that influence behavioral choices made by D. suzukii when selecting hosts, or the consequences that their offspring experience when developing in different environments. Using a nutritional geometry methodology, we found that the ratio of proteins‐to‐carbohydrates (P:C) present in media greatly influenced adult D. suzukii behavior and subsequent offspring development. Whereas adult flies showed a strong bias in their oviposition and association behaviors toward carbohydrate‐rich foods, larval survival and eclosion rate were strongly dependent on protein availability. Here, we explore the preference–performance hypothesis (PPH), in which females are predicted to oviposit on medias that provide the greatest offspring benefits, in regard to its relevance in D. suzukii behavior and consequences for management. Our results provide valuable insight into the ecology and evolution of this species that may hopefully lead to more effective management strategies.     

Reproductive and post‐reproductive life history of wild‐caught Drosophila melanogaster under laboratory conditions

The life history of the fruit fly (Drosophila melanogaster) is well understood, but fitness components are rarely measured by following single individuals over their lifetime, thereby limiting insights into lifetime reproductive success, reproductive senescence and post‐reproductive lifespan. Moreover, most studies have examined long‐established laboratory strains rather than freshly caught individuals and may thus be confounded by adaptation to laboratory culture, inbreeding or mutation accumulation. Here, we have followed the life histories of individual females from three recently caught, non‐laboratory‐adapted wild populations of D. melanogaster. Populations varied in a number of life‐history traits, including ovariole number, fecundity, hatchability and lifespan. To describe individual patterns of age‐specific fecundity, we developed a new model that allowed us to distinguish four phases during a female's life: a phase of reproductive maturation, followed by a period of linear and then exponential decline in fecundity and, finally, a post‐ovipository period. Individual females exhibited clear‐cut fecundity peaks, which contrasts with previous analyses, and post‐peak levels of fecundity declined independently of how long females lived. Notably, females had a pronounced post‐reproductive lifespan, which on average made up 40% of total lifespan. Post‐reproductive lifespan did not differ among populations and was not correlated with reproductive fitness components, supporting the hypothesis that this period is a highly variable, random ‘add‐on’ at the end of reproductive life rather than a correlate of selection on reproductive fitness. Most life‐history traits were positively correlated, a pattern that might be due to genotype by environment interactions when wild flies are brought into a novel laboratory environment but that is unlikely explained by inbreeding or positive mutational covariance caused by mutation accumulation.     

The adaptive significance of chromosomal inversion polymorphisms in Drosophila melanogaster

Chromosomal inversions, structural mutations that reverse a segment of a chromosome, cause suppression of recombination in the heterozygous state. Several studies have shown that inversion polymorphisms can form clines or fluctuate predictably in frequency over seasonal time spans. These observations prompted the hypothesis that chromosomal rearrangements might be subject to spatially and/or temporally varying selection. Here, we review what has been learned about the adaptive significance of inversion polymorphisms in the vinegar fly Drosophila melanogaster, the species in which they were first discovered by Sturtevant in 1917. A large body of work provides compelling evidence that several inversions in this system are adaptive; however, the precise selective mechanisms that maintain them polymorphic in natural populations remain poorly understood. Recent advances in population genomics, modelling and functional genetics promise to greatly improve our understanding of this long‐standing and fundamental problem in the near future.     

Population genomic analysis uncovers African and European admixture in Drosophila melanogaster populations from the south‐eastern United States and Caribbean Islands

Drosophila melanogaster is postulated to have colonized North America in the past several 100 years in two waves. Flies from Europe colonized the east coast United States while flies from Africa inhabited the Caribbean, which if true, make the south‐east US and Caribbean Islands a secondary contact zone for African and European D. melanogaster. This scenario has been proposed based on phenotypes and limited genetic data. In our study, we have sequenced individual whole genomes of flies from populations in the south‐east US and Caribbean Islands and examined these populations in conjunction with population sequences from the west coast US, Africa, and Europe. We find that west coast US populations are closely related to the European population, likely reflecting a rapid westward expansion upon first settlements into North America. We also find genomic evidence of African and European admixture in south‐east US and Caribbean populations, with a clinal pattern of decreasing proportions of African ancestry with higher latitude. Our genomic analysis of D. melanogaster populations from the south‐east US and Caribbean Islands provides more evidence for the Caribbean Islands as the source of previously reported novel African alleles found in other east coast US populations. We also find the border between the south‐east US and the Caribbean island to be the admixture hot zone where distinctly African‐like Caribbean flies become genomically more similar to European‐like south‐east US flies. Our findings have important implications for previous studies examining the generation of east coast US clines via selection.     

Local adaptation of a Drosophila parasitoid: habitat‐specific differences in thermal reaction norms

Local climate is an important source of selection on thermal reaction norms that has been well investigated in cline studies, where populations sampled along altitudinal or latitudinal gradients are compared. Several biotic factors vary with climate, but are rarely integrated as alternative agents of selection to climatic factors. We tested the hypothesis that habitat may select for thermal reaction norms and magnitude of phenotypic plasticity in a drosophila parasitoid, independently of the climate of origin. We sampled populations of Leptopilina boulardi, a Drosophila parasitoid in two different habitats, orchards and forests. Orchards offer laying opportunities over small distances for parasitoids, with a low variability in the number of hosts per patch, while forests offer more dispersed and more variable patches. The sampling was realized in a temperate and a Mediterranean climate. We measured egg load, volume of eggs, longevity and lipid content for parasitoids reared at two temperatures. Reaction norms were opposite for populations from forests and orchards for investment in reproduction, independently of the climate of origin. The maximal investment of resources in reproduction occurred at the lower temperature in orchards and the higher temperature in forests. Host distribution differences between habitats may explain these opposite reaction norms. We also observed a flatter reaction norm for egg load in forests than in orchards. This relative canalization may have been selected in response to the higher variability in laying opportunities observed in forests. Our results demonstrate the potential role of resource distribution in evolution of thermal plasticity.