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.

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.

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.

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.

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.

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.

The Drosophila simulans Genome Lacks the crystal-Stellate System.

The cry-Ste system is a genetic interaction system between heterochromatin and euchromatin in Drosophila melanogaster, regulated via the piRNA pathway. Deregulation of this system leads to meiotic defects and male sterility. Although the cry-Ste system is peculiar to D. melanogaster, ancestors of Ste and Su(Ste) elements are present in the three closely related species, D. simulans, D. sechellia, and D. mauritiana. The birth, evolution, and maintenance of this genetic system in Drosophila melanogaster are of interest. We investigate the presence of sequences homologous to cry and Ste elements in the simulans complex and describe their chromosomal distribution. The organization and expression of cry- and Ste-like sequences were further characterized in the D. simulans genome. Our results allow us to conclude that the cry-Ste genetic interaction system is absent in the D. simulans genome.

Testing the Drosophila maternal haploid gene for functional divergence and a role in hybrid incompatibility.

Crosses between Drosophila simulans females and Drosophila melanogaster males produce viable F1 sons and poorly viable F1 daughters. Unlike most hybrid incompatibilities, this hybrid incompatibility violates Haldane's rule, the observation that incompatibilities preferentially affect the heterogametic sex. Furthermore, it has a different genetic basis than hybrid lethality in the reciprocal cross, with the causal allele in Drosophila melanogaster being a large species-specific block of complex satellite DNA on its X chromosome known as the 359-bp satellite, rather than a protein-coding locus. The causal allele(s) in Drosophila simulans are unknown but likely involve maternally expressed genes or factors since the F1 females die during early embryogenesis. The maternal haploid (mh) gene is an intriguing candidate because it is expressed maternally and its protein product localizes to the 359-bp repeat. We found that this gene has diverged extensively between Drosophila melanogaster and Drosophila simulans. This observation led to the hypothesis that Drosophila melanogaster mh may have coevolved with the 359-bp repeat and that hybrid incompatibility thus results from the absence of a coevolved mh allele in Drosophila simulans. We tested for the functional divergence of mh by creating matched transformants of Drosophila melanogaster and Drosophila simulans orthologs in both Drosophila melanogaster and Drosophila simulans strains. Surprisingly, we find that Drosophila simulans mh fully complements the female sterile phenotype of Drosophila melanogaster mh mutations. Contrary to our hypothesis, we find no evidence that adding a Drosophila melanogaster mh gene to Drosophila simulans increases hybrid viability.

The Transposition Rate Has Little Influence on the Plateauing Level of the P-element.

The popular trap model assumes that the invasions of transposable elements (TEs) in mammals and invertebrates are stopped by piRNAs that emerge after insertion of the TE into a piRNA cluster. It remains, however, still unclear which factors influence the dynamics of TE invasions. The activity of the TE (i.e., transposition rate) is one frequently discussed key factor. Here we take advantage of the temperature-dependent activity of the P-element, a widely studied eukaryotic TE, to test how TE activity affects the dynamics of a TE invasion. We monitored P-element invasion dynamics in experimental Drosophila simulans populations at hot and cold culture conditions. Despite marked differences in transposition rates, the P-element reached very similar copy numbers at both temperatures. The reduction of the insertion rate upon approaching the copy number plateau was accompanied by similar amounts of piRNAs against the P-element at both temperatures. Nevertheless, we also observed fewer P-element insertions in piRNA clusters than expected, which is not compatible with a simple trap model. The ping-pong cycle, which degrades TE transcripts, becomes typically active after the copy number plateaued. We generated a model, with few parameters, that largely captures the observed invasion dynamics. We conclude that the transposition rate has at the most only a minor influence on TE abundance, but other factors, such as paramutations or selection against TE insertions are shaping the TE composition.

The fate of a suppressed X-linked meiotic driver: experimental evolution in Drosophila simulans.

Sex-ratio (SR) meiotic drivers are X-linked selfish genetic elements that promote their own transmission by preventing the production of Y-bearing sperm, which usually lowers male fertility. The spread of SR drivers in populations is expected to trigger the evolution of unlinked drive suppressors, a theoretically predicted co-evolution that has been observed in nature. Once completely suppressed, the drivers are expected either to decline if they still affect the fitness of their carriers, or to evolve randomly and possibly get fixed if the suppressors eliminate their deleterious effects. To explore this issue, we used the Paris sex-ratio system of Drosophila simulans in which drive results from the joint effect of two elements on the X chromosome: a segmental duplication and a deficient allele of the HP1D2 gene. We set up six experimental populations starting with 2/3 of X chromosomes carrying both elements (XSR) in a fully suppressing background. We let them evolve independently during almost a hundred generations under strong sexual competition, a condition known to cause the rapid disappearance of unsuppressed Paris XSR in previous experimental populations. In our study, the fate of XSR chromosomes varied among populations, from extinction to their maintenance at a frequency close to the starting one. While the reasons for these variable outcomes are still to be explored, our results show that complete suppression can prevent the demise of an otherwise deleterious XSR chromosome, turning a genetic conflict into cooperation between unlinked loci. Observations in natural populations suggest a contrasting fate of the two elements: disappearance of the duplication and maintenance of deficient HP1D2 alleles.

Pool-GWAS on reproductive dormancy in Drosophila simulans suggests a polygenic architecture.

The genetic basis of adaptation to different environments has been of long-standing interest to evolutionary biologists. Dormancy is a well-studied adaptation to facilitate overwintering. In Drosophila melanogaster, a moderate number of genes with large effects have been described, which suggests a simple genetic basis of dormancy. On the other hand, genome-wide scans for dormancy suggest a polygenic architecture in insects. In D. melanogaster, the analysis of the genetic architecture of dormancy is complicated by the presence of cosmopolitan inversions. Here, we performed a genome-wide scan to characterize the genetic basis of this ecologically extremely important trait in the sibling species of D. melanogaster, D. simulans that lacks cosmopolitan inversions. We performed Pool-GWAS in a South African D. simulans population for dormancy incidence at 2 temperature regimes (10 and 12°C, LD 10:14). We identified several genes with SNPs that showed a significant association with dormancy (P-value < 1e-13), but the overall modest response suggests that dormancy is a polygenic trait with many loci of small effect. Our results shed light on controversies on reproductive dormancy in Drosophila and have important implications for the characterization of the genetic basis of this trait.

Unique structure and positive selection promote the rapid divergence of Drosophila Y chromosomes.

Y chromosomes across diverse species convergently evolve a gene-poor, heterochromatic organization enriched for duplicated genes, LTR retrotransposons, and satellite DNA. Sexual antagonism and a loss of recombination play major roles in the degeneration of young Y chromosomes. However, the processes shaping the evolution of mature, already degenerated Y chromosomes are less well-understood. Because Y chromosomes evolve rapidly, comparisons between closely related species are particularly useful. We generated de novo long-read assemblies complemented with cytological validation to reveal Y chromosome organization in three closely related species of the Drosophila simulans complex, which diverged only 250,000 years ago and share >98% sequence identity. We find these Y chromosomes are divergent in their organization and repetitive DNA composition and discover new Y-linked gene families whose evolution is driven by both positive selection and gene conversion. These Y chromosomes are also enriched for large deletions, suggesting that the repair of double-strand breaks on Y chromosomes may be biased toward microhomology-mediated end joining over canonical non-homologous end-joining. We propose that this repair mechanism contributes to the convergent evolution of Y chromosome organization across organisms.