Evolution of Phenotypic Variance Provides Insights into the Genetic Basis of Adaptation.

Most traits are polygenic, and the contributing loci can be identified by genome-wide association studies. The genetic basis of adaptation (adaptive architecture) is, however, difficult to characterize. Here, we propose to study the adaptive architecture of traits by monitoring the evolution of their phenotypic variance during adaptation to a new environment in well-defined laboratory conditions. Extensive computer simulations show that the evolution of phenotypic variance in a replicated experimental evolution setting can distinguish between oligogenic and polygenic adaptive architectures. We compared gene expression variance in male Drosophila simulans before and after 100 generations of adaptation to a novel hot environment. The variance change in gene expression was indistinguishable for genes with and without a significant change in mean expression after 100 generations of evolution. We suggest that the majority of adaptive gene expression evolution can be explained by a polygenic architecture. We propose that tracking the evolution of phenotypic variance across generations can provide an approach to characterize the adaptive architecture.

Wolbachia infection negatively impacts Drosophila simulans heat tolerance in a strain- and trait-specific manner.

The susceptibility of insects to rising temperatures has largely been measured by their ability to survive thermal extremes. However, the capacity for maternally inherited endosymbionts to influence insect heat tolerance has been overlooked. Further, while some studies have addressed the impact of heat on traits like fertility, which can decline at temperatures below lethal thermal limits, none have considered the impact of endosymbionts. Here, we assess the impact of three Wolbachia strains (wRi, wAu and wNo) on the survival and fertility of Drosophila simulans exposed to heat stress during development or as adults. The effect of Wolbachia infection on heat tolerance was generally small and trait/strain specific. Only the wNo infection significantly reduced the survival of adult males after a heat shock. When exposed to fluctuating heat stress during development, the wRi and wAu strains reduced egg-to-adult survival but only the wNo infection reduced male fertility. Wolbachia densities of all three strains decreased under developmental heat stress, but reductions occurred at temperatures above those that reduced host fertility. These findings emphasize the necessity to account for endosymbionts and their effect on both survival and fertility when investigating insect responses to heat stress.

Elevated ozone disrupts mating boundaries in drosophilid flies.

Animals employ different strategies to establish mating boundaries between closely related species, with sex pheromones often playing a crucial role in identifying conspecific mates. Many of these pheromones have carbon-carbon double bonds, making them vulnerable to oxidation by certain atmospheric oxidant pollutants, including ozone. Here, we investigate whether increased ozone compromises species boundaries in drosophilid flies. We show that short-term exposure to increased levels of ozone degrades pheromones of Drosophila melanogaster, D. simulans, D. mauritiana, as well as D. sechellia, and induces hybridization between some of these species. As many of the resulting hybrids are sterile, this could result in local population declines. However, hybridization between D. simulans and D. mauritiana as well as D. simulans and D. sechellia results in fertile hybrids, of which some female hybrids are even more attractive to the males of the parental species. Our experimental findings indicate that ozone pollution could potentially induce breakdown of species boundaries in insects.

The evolution of seminal fluid gene expression and postmating reproductive isolation in Drosophila.

Seminal fluid protein (Sfp) genes show, in general, a higher rate of sequence divergence than genes from other categories, which is often attributed to forms of postcopulatory sexual selection or sexual conflict. Recently, the relaxation of selective constraints has been proposed as an alternative explanation for the rapid sequence evolution of Sfps and other genes with sex-limited expression. The expression of Sfp genes is a likely target of selection, but the evolution of differences in their expression levels is less understood. Here, we explore both polymorphism and divergence in Sfp gene expression between Drosophila melanogaster and Drosophila simulans, how selection might have influenced their expression, and whether changes in expression might trigger the evolution of reproductive isolating barriers. In our analysis, Sfp genes showed higher divergence, but not higher polymorphism, in expression than the average male reproductive glands gene. Sfp genes with reproductive-tissue-specific expression were enriched for both directional and stabilizing selection, while relaxed selection was the predominant mode of evolution among Sfp genes with any other nonreproductive tissue-specific or nontissue-specific expression. The knockdown of single genes known to affect intraspecific sperm competition, and with patterns of expression divergence and polymorphism suggestive of directional selection, was not enough to break down postmating reproductive isolation barriers between species. Our results identify the expression of male-specific Sfp genes as an enriched target of selection and suggest a complex molecular relationship between postcopulatory sexual selection on a single gene's expression and its effect on the onset of speciation.

Wolbachia Induces Structural Defects Harmful to Drosophila simulans Riverside Spermiogenesis.

The relationship between cytoplasmic incompatibility and the obligate intracellular alphaproteobacteria Wolbachia has for a long time been reported. Although the molecular mechanisms responsible for this reproductive alteration are beginning to be understood, the effects of Wolbachia on germ cell structure and dynamics have not yet been fully investigated. We report here that the presence of Wolbachia in infected cysts of elongating spermatids is associated with major structural defects that become more evident in mature sperm. We find mitochondrial defects, an improper axoneme structure, reduced sperm numbers, and individualization failures. The large heterogeneous variety of the ultrastructural defects found in elongating spermatids and mature sperm provide the first cytological evidence for the reduced fertility associated with Wolbachia infection in Drosophila simulans males. The observed abnormalities could be the result of the mechanical stress induced by the high bacteria numbers during the process of spermatid elongation, rather than the result of the released factors affecting the proper morphogenesis of the germ cells. Moreover, high Wolbachia densities in male germ cells may not be appropriate for causing cytoplasmic incompatibility as the bacteria are harmful for spermatid differentiation, leading to abnormal sperm that is unlikely to be functional.

High levels of intra-strain structural variation in Drosophila simulans X pericentric heterochromatin.

Large genome structural variations can impact genome regulation and integrity. Repeat-rich regions like pericentric heterochromatin are vulnerable to structural rearrangements although we know little about how often these rearrangements occur over evolutionary time. Repetitive genome regions are particularly difficult to study with genomic approaches, as they are missing from most genome assemblies. However, cytogenetic approaches offer a direct way to detect large rearrangements involving pericentric heterochromatin. Here, we use a cytogenetic approach to reveal large structural rearrangements associated with the X pericentromeric region of Drosophila simulans. These rearrangements involve large blocks of satellite DNA-the 500-bp and Rsp-like satellites-which colocalize in the X pericentromeric heterochromatin. We find that this region is polymorphic not only among different strains, but between isolates of the same strain from different labs, and even within individual isolates. On the one hand, our observations raise questions regarding the potential impact of such variation at the phenotypic level and our ability to control for such genetic variability. On the other hand, this highlights the very rapid turnover of the pericentric heterochromatin most likely associated with genomic instability of the X pericentromere. It represents a unique opportunity to study the dynamics of pericentric heterochromatin, the evolution of associated satellites on a very short time scale, and to better understand how structural variation arises.

The genomic distribution of transposable elements is driven by spatially variable purifying selection.

It is widely accepted that the genomic distribution of transposable elements (TEs) mainly reflects the outcome of purifying selection and insertion bias (1). Nevertheless, the relative importance of these two evolutionary forces could not be tested thoroughly. Here, we introduce an experimental system, which allows separating purifying selection from TE insertion bias. We used experimental evolution to study the TE insertion patterns in Drosophila simulans founder populations harboring 1040 insertions of an active P-element. After 10 generations at a large population size, we detected strong selection against P-element insertions. The exception were P-element insertions in genomic regions for which a strong insertion bias has been proposed (2-4). Because recurrent P-element insertions cannot explain this pattern, we conclude that purifying selection, with variable strength along the chromosomes, is the major determinant of the genomic distribution of P-elements. Genomic regions with relaxed purifying selection against P-element insertions exhibit normal levels of purifying selection against base substitutions. This suggests that different types of purifying selection operate on base substitutions and P-element insertions. Our results highlight the power of experimental evolution to understand basic evolutionary processes, which are difficult to infer from patterns of natural variation alone.

A Quantitative, Genome-Wide Analysis in Drosophila Reveals Transposable Elements' Influence on Gene Expression Is Species-Specific.

Transposable elements (TEs) are parasite DNA sequences that are able to move and multiply along the chromosomes of all genomes. They can be controlled by the host through the targeting of silencing epigenetic marks, which may affect the chromatin structure of neighboring sequences, including genes. In this study, we used transcriptomic and epigenomic high-throughput data produced from ovarian samples of several Drosophila melanogaster and Drosophila simulans wild-type strains, in order to finely quantify the influence of TE insertions on gene RNA levels and histone marks (H3K9me3 and H3K4me3). Our results reveal a stronger epigenetic effect of TEs on ortholog genes in D. simulans compared with D. melanogaster. At the same time, we uncover a larger contribution of TEs to gene H3K9me3 variance within genomes in D. melanogaster, which is evidenced by a stronger correlation of TE numbers around genes with the levels of this chromatin mark in D. melanogaster. Overall, this work contributes to the understanding of species-specific influence of TEs within genomes. It provides a new light on the considerable natural variability provided by TEs, which may be associated with contrasted adaptive and evolutionary potentials.

Rapid evolutionary diversification of the flamenco locus across simulans clade Drosophila species.

Suppression of transposable elements (TEs) is paramount to maintain genomic integrity and organismal fitness. In D. melanogaster, the flamenco locus is a master suppressor of TEs, preventing the mobilization of certain endogenous retrovirus-like TEs from somatic ovarian support cells to the germline. It is transcribed by Pol II as a long (100s of kb), single-stranded, primary transcript, and metabolized into ~24-32 nt Piwi-interacting RNAs (piRNAs) that target active TEs via antisense complementarity. flamenco is thought to operate as a trap, owing to its high content of recent horizontally transferred TEs that are enriched in antisense orientation. Using newly-generated long read genome data, which is critical for accurate assembly of repetitive sequences, we find that flamenco has undergone radical transformations in sequence content and even copy number across simulans clade Drosophilid species. Drosophila simulans flamenco has duplicated and diverged, and neither copy exhibits synteny with D. melanogaster beyond the core promoter. Moreover, flamenco organization is highly variable across D. simulans individuals. Next, we find that D. simulans and D. mauritiana flamenco display signatures of a dual-stranded cluster, with ping-pong signals in the testis and/or embryo. This is accompanied by increased copy numbers of germline TEs, consistent with these regions operating as functional dual-stranded clusters. Overall, the physical and functional diversity of flamenco orthologs is testament to the extremely dynamic consequences of TE arms races on genome organization, not only amongst highly related species, but even amongst individuals.

The Evolutionary History of Drosophila simulans Y Chromosomes Reveals Molecular Signatures of Resistance to Sex Ratio Meiotic Drive.

The recent evolutionary history of the Y chromosome in Drosophila simulans, a worldwide species of Afrotropical origin, is closely linked to that of X-linked meiotic drivers (Paris system). The spread of the Paris drivers in natural populations has elicited the selection of drive-resistant Y chromosomes. To infer the evolutionary history of the Y chromosome in relation to the Paris drive, we sequenced 21 iso-Y lines, each carrying a Y chromosome from a different location. Among them, 13 lines carry a Y chromosome that is able to counteract the effect of the drivers. Despite their very different geographical origins, all sensitive Y's are highly similar, suggesting that they share a recent common ancestor. The resistant Y chromosomes are more divergent and segregate in four distinct clusters. The phylogeny of the Y chromosome confirms that the resistant lineage predates the emergence of Paris drive. The ancestry of the resistant lineage is further supported by the examination of Y-linked sequences in the sister species of D. simulans, Drosophila sechellia and Drosophila mauritiana. We also characterized the variation in repeat content among Y chromosomes and identified multiple simple satellites associated with resistance. Altogether, the molecular polymorphism allows us to infer the demographic and evolutionary history of the Y chromosome and provides new insights on the genetic basis of resistance.

Essential and recurrent roles for hairpin RNAs in silencing de novo sex chromosome conflict in Drosophila simulans.

Meiotic drive loci distort the normally equal segregation of alleles, which benefits their own transmission even in the face of severe fitness costs to their host organism. However, relatively little is known about the molecular identity of meiotic drivers, their strategies of action, and mechanisms that can suppress their activity. Here, we present data from the fruitfly Drosophila simulans that address these questions. We show that a family of de novo, protamine-derived X-linked selfish genes (the Dox gene family) is silenced by a pair of newly emerged hairpin RNA (hpRNA) small interfering RNA (siRNA)-class loci, Nmy and Tmy. In the w[XD1] genetic background, knockout of nmy derepresses Dox and MDox in testes and depletes male progeny, whereas knockout of tmy causes misexpression of PDox genes and renders males sterile. Importantly, genetic interactions between nmy and tmy mutant alleles reveal that Tmy also specifically maintains male progeny for normal sex ratio. We show the Dox loci are functionally polymorphic within D. simulans, such that both nmy-associated sex ratio bias and tmy-associated sterility can be rescued by wild-type X chromosomes bearing natural deletions in different Dox family genes. Finally, using tagged transgenes of Dox and PDox2, we provide the first experimental evidence Dox family genes encode proteins that are strongly derepressed in cognate hpRNA mutants. Altogether, these studies support a model in which protamine-derived drivers and hpRNA suppressors drive repeated cycles of sex chromosome conflict and resolution that shape genome evolution and the genetic control of male gametogenesis.

Regulatory logic of endogenous RNAi in silencing de novo genomic conflicts.

Although the biological utilities of endogenous RNAi (endo-RNAi) have been largely elusive, recent studies reveal its critical role in the non-model fruitfly Drosophila simulans to suppress selfish genes, whose unchecked activities can severely impair spermatogenesis. In particular, hairpin RNA (hpRNA) loci generate endo-siRNAs that suppress evolutionary novel, X-linked, meiotic drive loci. The consequences of deleting even a single hpRNA (Nmy) in males are profound, as such individuals are nearly incapable of siring male progeny. Here, comparative genomic analyses of D. simulans and D. melanogaster mutants of the core RNAi factor dcr-2 reveal a substantially expanded network of recently-emerged hpRNA-target interactions in the former species. The de novo hpRNA regulatory network in D. simulans provides insight into molecular strategies that underlie hpRNA emergence and their potential roles in sex chromosome conflict. In particular, our data support the existence of ongoing rapid evolution of Nmy/Dox-related networks, and recurrent targeting of testis HMG-box loci by hpRNAs. Importantly, the impact of the endo-RNAi network on gene expression flips the convention for regulatory networks, since we observe strong derepression of targets of the youngest hpRNAs, but only mild effects on the targets of the oldest hpRNAs. These data suggest that endo-RNAi are especially critical during incipient stages of intrinsic sex chromosome conflicts, and that continual cycles of distortion and resolution may contribute to speciation.

Temperature adaptation and its impact on the shape of performance curves in Drosophila populations.

Understanding how species adapt to different temperatures is crucial to predict their response to global warming, and thermal performance curves (TPCs) have been employed recurrently to study this topic. Nevertheless, fundamental questions regarding how thermodynamic constraints and evolution interact to shape TPCs in lineages inhabiting different environments remain unanswered. Here, we study Drosophila simulans along a latitudinal gradient spanning 3000 km to test opposing hypotheses based on thermodynamic constrains (hotter-is-better) versus biochemical adaptation (jack-of-all-temperatures) as primary determinants of TPCs variation across populations. We compare thermal responses in metabolic rate and the egg-to-adult survival as descriptors of organismal performance and fitness, respectively, and show that different descriptors of TPCs vary in tandem with mean environmental temperatures, providing strong support to hotter-is-better. Thermodynamic constraints also resulted in a strong negative association between maximum performance and thermal breadth. Lastly, we show that descriptors of TPCs for metabolism and egg-to-adult survival are highly correlated, providing evidence of co-adaptation, and that curves for egg-to-adult survival are systematically narrower and displaced toward lower temperatures. Taken together, our results support the pervasive role of thermodynamics constraining thermal responses in Drosophila populations along a latitudinal gradient, that are only partly compensated by evolutionary adaptation.

Evolution of Metabolome and Transcriptome Supports a Hierarchical Organization of Adaptive Traits.

Most organismal phenotypes have a polygenic basis, which enables adaptive phenotypic responses on ecological time scales. While adaptive phenotypic changes are highly parallel in replicate populations, this does not apply to the contributing loci. In particular for small populations, the same phenotypic shift can be fueled by different sets of alleles at alternative loci (genetic redundancy). Although this phenomenon is empirically well supported, the molecular basis of the genetic redundancy is not yet understood. To fill this gap, we compared the heterogeneity of the evolutionary transcriptomic and metabolomic response in ten Drosophila simulans populations which evolved parallel high-level phenotypic changes in a novel temperature environment but used different allelic combinations of alternative loci. We showed that the metabolome evolved more parallel than the transcriptome, confirming a hierarchical organization of molecular phenotypes. Different sets of genes responded in each evolved population but led to the enrichment of similar biological functions and a consistent metabolic profile. Since even the metabolomic response was still highly heterogeneous across evolved populations, we propose that selection may operate on pathways/networks.

Variation in mutation, recombination, and transposition rates in Drosophila melanogaster and Drosophila simulans.

The rates of mutation, recombination, and transposition are core parameters in models of evolution. They impact genetic diversity, responses to ongoing selection, and levels of genetic load. However, even for key evolutionary model species such as Drosophila melanogaster and Drosophila simulans, few estimates of these parameters are available, and we have little idea of how rates vary between individuals, sexes, or populations. Knowledge of this variation is fundamental for parameterizing models of genome evolution. Here, we provide direct estimates of mutation, recombination, and transposition rates and their variation in a West African and a European population of D. melanogaster and a European population of D. simulans Across 89 flies, we observe 58 single-nucleotide mutations, 286 crossovers, and 89 transposable element (TE) insertions. Compared to the European D. melanogaster, we find the West African population has a lower mutation rate (1.67 × 10-9 site-1 gen-1 vs. 4.86 × 10-9 site-1 gen-1) and a lower transposition rate (8.99 × 10-5 copy-1 gen-1 vs. 23.36 × 10-5 copy-1 gen-1), but a higher recombination rate (3.44 cM/Mb vs. 2.06 cM/Mb). The European D. simulans population has a similar mutation rate to European D. melanogaster, but a significantly higher recombination rate and a lower, but not significantly different, transposition rate. Overall, we find paternal-derived mutations are more frequent than maternal ones in both species. Our study quantifies the variation in rates of mutation, recombination, and transposition among different populations and sexes, and our direct estimates of these parameters in D. melanogaster and D. simulans will benefit future studies in population and evolutionary genetics.

Sex-Biased Expression Is Associated With Chromatin State in Drosophila melanogaster and Drosophila simulans.

In Drosophila melanogaster and D. simulans head tissue, 60% of orthologous genes show evidence of sex-biased expression in at least one species. Of these, ∼39% (2,192) are conserved in direction. We hypothesize enrichment of open chromatin in the sex where we see expression bias and closed chromatin in the opposite sex. Male-biased orthologs are significantly enriched for H3K4me3 marks in males of both species (∼89% of male-biased orthologs vs. ∼76% of unbiased orthologs). Similarly, female-biased orthologs are significantly enriched for H3K4me3 marks in females of both species (∼90% of female-biased orthologs vs. ∼73% of unbiased orthologs). The sex-bias ratio in female-biased orthologs was similar in magnitude between the two species, regardless of the closed chromatin (H3K27me2me3) marks in males. However, in male-biased orthologs, the presence of H3K27me2me3 in both species significantly reduced the correlation between D. melanogaster sex-bias ratio and the D. simulans sex-bias ratio. Male-biased orthologs are enriched for evidence of positive selection in the D. melanogaster group. There are more male-biased genes than female-biased genes in both species. For orthologs with gains/losses of sex-bias between the two species, there is an excess of male-bias compared to female-bias, but there is no consistent pattern in the relationship between H3K4me3 or H3K27me2me3 chromatin marks and expression. These data suggest chromatin state is a component of the maintenance of sex-biased expression and divergence of sex-bias between species is reflected in the complexity of the chromatin status.

[The P-Element Has Not Significant Effect on the Drosophila simulans Viability].

Cases of horizontal transfer of transposable elements (TEs) between species are known for the Drosophilidae family. In the middle of the last century, the case of horizontal transfer of the P-element from the Drosophila willistoni to the D. melanogaster was described. A novel P-element invasion into the D. simulans genome from D. melanogaster occurred approximately 10 years ago. Currently, the P-element has spread across all D. melanogaster population and 30% of D. simulans populations in Europe, Africa and America. In this paper, we investigated the presence of the P-element in D. simulans lines caught in different years in three Asian populations (Tashkent, Nalchik and Sakhalin Island). We also examined the physiological characteristics (cytotype, lifespan, fecundity and locomotor activity) of D. simulans lines with and without the P-element to determine the significance of this new mobile element in the genome. The P-element was found in lines isolated from nature after 2012. The number of P-element copies per genome (two-to-three dozen according to fluorescence in situ hybridization data) was greater than in the American and comparable to the African populations. There were signs of intraspecific hybrid dysgenesis for some pairs of lines. However, in general the presence of the P-element did not adversely affect the physiological characteristics. Either adaptation to the new TE occurs very quickly, or the rate of movement of the P-element is so insignificant that its appearance in the genome remains unnoticed.

Evolutionary conservation and diversification of auditory neural circuits that process courtship songs in Drosophila.

Acoustic communication signals diversify even on short evolutionary time scales. To understand how the auditory system underlying acoustic communication could evolve, we conducted a systematic comparison of the early stages of the auditory neural circuit involved in song information processing between closely-related fruit-fly species. Male Drosophila melanogaster and D. simulans produce different sound signals during mating rituals, known as courtship songs. Female flies from these species selectively increase their receptivity when they hear songs with conspecific temporal patterns. Here, we firstly confirmed interspecific differences in temporal pattern preferences; D. simulans preferred pulse songs with longer intervals than D. melanogaster. Primary and secondary song-relay neurons, JO neurons and AMMC-B1 neurons, shared similar morphology and neurotransmitters between species. The temporal pattern preferences of AMMC-B1 neurons were also relatively similar between species, with slight but significant differences in their band-pass properties. Although the shift direction of the response property matched that of the behavior, these differences are not large enough to explain behavioral differences in song preferences. This study enhances our understanding of the conservation and diversification of the architecture of the early-stage neural circuit which processes acoustic communication signals.

A comparative study of natural variation in hemolymph glucose levels under different dietary sugar conditions in Drosophila melanogaster and D. simulans.

Physiological responses to environmental changes play important roles in adaptive evolution. In particular, homeostatic regulatory systems that maintain constant circulating glucose levels are crucial in animals. However, variation in circulating glucose levels and the genetic effects on phenotypic variation in natural populations remain to be clarified. Here, we investigated the hemolymph glucose levels in natural populations of Drosophila melanogaster and its sibling species, D. simulans, in Japan. We quantified hemolymph glucose concentrations in third instar larvae of 27 lines for each species, which were reared on either glucose-free or glucose-rich food. In both species, genetic variation was not a major component of phenotypic variation on either glucose-free or glucose-rich food. The hemolymph glucose concentrations were much higher in D. simulans than in D. melanogaster. Genetic variance was larger in D. simulans than in D. melanogaster. The observed differences between the two species may be associated with the much more recent colonization history of D. simulans populations in Japan and/or the tolerance to environmental stresses. Our findings suggest that natural selection acting on hemolymph glucose levels in D. melanogaster is different from that in D. simulans.

Levels of P-element-induced hybrid dysgenesis in Drosophila simulans are uncorrelated with levels of P-element piRNAs.

Transposable elements (TEs) are genomic parasites that proliferate within host genomes, and which can also invade new species. The P-element, a DNA-based TE, recently invaded two Drosophila species: Drosophila melanogaster in the 20th century, and D. simulans in the 21st. In both species, lines collected before the invasion are susceptible to "hybrid dysgenesis", a syndrome of abnormal phenotypes apparently due to P-element-inflicted DNA damage. In D. melanogaster, lines collected after the invasion have evolved a maternally acting mechanism that suppresses hybrid dysgenesis, with extensive work showing that PIWI-interacting small RNAs (piRNAs) are a key factor in this suppression. Most of these studies use lines collected many generations after the initial P-element invasion. Here, we study D. simulans collected early, as well as late in the P-element invasion of this species. Like D. melanogaster, D. simulans from late in the invasion show strong resistance to hybrid dysgenesis and abundant P-element-derived piRNAs. Lines collected early in the invasion, however, show substantial variation in how much they suffer from hybrid dysgenesis, with some lines highly resistant. Surprisingly, although, these resistant lines do not show high levels of cognate maternal P-element piRNAs; in these lines, it may be that other mechanisms suppress hybrid dysgenesis.