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.

Identification and genetic analysis of a pervasive 'needle-eye' sperm phenotype in Drosophila sterile hybrid males.

Interspecies hybrid sterility has been extensively studied, especially in the genus Drosophila. Hybrid sterility is more often found in the heterogametic (XY or ZW) sex, a trend called Haldane's rule. Although this phenomenon is pervasive, identification of a common genetic mechanism remains elusive, with modest support found for a range of potential theories. Here, we identify a single precise morphological phenotype, which we call 'needle-eye sperm', that is associated with hybrid sterility in three separate species pairs that span the Drosophila genus. The nature of the phenotype indicates a common point of meiotic failure in sterile hybrid males. We used 10 generations of backcross selection paired with whole-genome pooled sequencing to genetically map the regions underlying the needle-eye (NE) sperm phenotype. Surprisingly, the sterility phenotype was present in ~50% of males even after 10 generations of backcrossing, and only a single region of the X chromosome was associated with sterility in one direction of backcross. Owing to the common phenotype among sterile male hybrids, and the strong effect of individual loci, further exploration of these findings may identify a universal mechanism for the evolution of hybrid sterility.

Population Density Affects Drosophila Male Pheromones in Laboratory-Acclimated and Natural Lines.

In large groups of vertebrates and invertebrates, aggregation can affect biological characters such as gene expression, physiological, immunological and behavioral responses. The insect cuticle is covered with hydrocarbons (cuticular hydrocarbons; CHCs) which reduce dehydration and increase protection against xenobiotics. Drosophila melanogaster and D. simulans flies also use some of their CHCs as contact pheromones. In these two sibling species, males also produce the volatile pheromone 11-cis-Vaccenyl acetate (cVa). To investigate the effect of insect density on the production of CHCs and cVa we compared the level of these male pheromones in groups of different sizes. These compounds were measured in six lines acclimated for many generations in our laboratory - four wild-type and one CHC mutant D. melanogaster lines plus one D. simulans line. Increasing the group size substantially changed pheromone amounts only in the four D. melanogaster wild-type lines. To evaluate the role of laboratory acclimation in this effect, we measured density-dependent pheromonal production in 21 lines caught in nature after 1, 12 and 25 generations in the laboratory. These lines showed varied effects which rarely persisted across generations. Although increasing group size often affected pheromone production in laboratory-established and freshly-caught D. melanogaster lines, this effect was not linear, suggesting complex determinants.

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 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.

Evolution of phenotypic variance in response to a novel hot environment.

Shifts in trait means are widely considered as evidence for adaptive responses, but the impact on phenotypic variance remains largely unexplored. Classic quantitative genetics provides a theoretical framework to predict how selection on phenotypic mean affects the variance. In addition to this indirect effect, it is also possible that the variance of the trait is the direct target of selection, but experimentally characterized cases are rare. Here, we studied gene expression variance of Drosophila simulans males before and after 100 generations of adaptation to a novel hot laboratory environment. In each of the two independently evolved populations, the variance of 125 and 97 genes was significantly reduced. We propose that the drastic loss in environmental complexity from nature to the laboratory may have triggered selection for reduced variance. Our observation that selection could drive changes in the variance of gene expression could have important implications for studies of adaptation processes in natural and experimental populations.

Comparing single- and mixed-species groups in fruit flies: differences in group dynamics, but not group formation.

Mixed-species groups describe active associations among individuals of 2 or more species at the same trophic level. Mixed-species groups are important to key ecological and evolutionary processes such as competition and predation, and research that ignores the presence of other species risks ignoring a key aspect of the environment in which social behavior is expressed and selected. Despite the defining emphasis of active formation for mixed-species groups, surprisingly little is known about the mechanisms by which mixed-species groups form. Furthermore, insects have been almost completely ignored in the study of mixed-species groups, despite their taxonomic importance and relative prominence in the study of single-species groups. Here, we measured group formation processes in Drosophila melanogaster and its sister species, Drosophila simulans. Each species was studied alone, and together, and one population of D. melanogaster was also studied both alone and with another, phenotypically distinct D. melanogaster population, in a nested-factorial design. This approach differs from typical methods of studying mixed-species groups in that we could quantitatively compare group formation between single-population, mixed-population, and mixed-species treatments. Surprisingly, we found no differences between treatments in the number, size, or composition of groups that formed, suggesting that single- and mixed-species groups form through similar mechanisms of active attraction. However, we found that mixed-species groups showed elevated interspecies male-male interactions, relative to interpopulation or intergenotype interactions in single-species groups. Our findings expand the conceptual and taxonomic study of mixed-species groups while raising new questions about the mechanisms of group formation broadly.

RNA virome diversity and Wolbachia infection in individual Drosophila simulans flies.

The endosymbiont bacteria of the genus Wolbachia are associated with multiple mutualistic effects on insect biology, including nutritional and antiviral properties. Members of the genus Wolbachia naturally occur in fly species of the genus Drosophila, providing an operational model host for studying how virome composition may be affected by its presence. Drosophila simulans populations can carry a variety of strains of members of the genus Wolbachia, with the wAu strain associated with strong antiviral protection under experimental conditions. We used D. simulans sampled from the Perth Hills, Western Australia, to investigate the potential virus protective effect of the wAu strain of Wolbachia on individual wild-caught flies. Our data revealed no appreciable variation in virus composition and abundance between individuals infected or uninfected with Wolbachia associated with the presence or absence of wAu. However, it remains unclear whether wAu might affect viral infection and host survival by increasing tolerance rather than inducing complete resistance. These data also provide new insights into the natural virome diversity of D. simulans. Despite the small number of individuals sampled, we identified a repertoire of RNA viruses, including nora virus, galbut virus, thika virus and La Jolla virus, that have been identified in other species of the genus Drosophila. Chaq virus-like sequences associated with galbut virus were also detected. In addition, we identified five novel viruses from the families Reoviridae, Tombusviridae, Mitoviridae and Bunyaviridae. Overall, this study highlights the complex interaction between Wolbachia and RNA virus infections and provides a baseline description of the natural virome of D. simulans.

Genes underlying species differences in cuticular hydrocarbon production between Drosophila melanogaster and D. simulans.

Surface chemical compounds are key components of survival and reproduction in many species. Cuticular hydrocarbons (CHCs) are chemical compounds produced by all insects that are used for both desiccation resistance and chemical communication, including communication related to mating. In the species pair of Drosophila melanogaster and D. simulans, female CHCs stimulate conspecific males to mate and repel heterospecific males. While CHCs are a critical contributor to both reproductive success within a species and isolation between species, few genes underlying species variation in CHC profiles are known. Here, we use genetic mapping of the 3rd chromosome to test a suite of candidate genes for interspecies variation in CHCs. Candidate gene CG5946 was found to be involved in species differences in the production of 7,11-heptacosadiene and 7-tricosene between D. melanogaster and D. simulans. This is therefore a new candidate locus contributing to species-specific variation in the CHC profile. In the process of mapping genes for CHCs, we also identified 29 candidate genes for the reduced survival or inviability of interspecies hybrids.

Sex-Ratio Meiotic Drive Shapes the Evolution of the Y Chromosome in Drosophila simulans.

The recent emergence and spread of X-linked segregation distorters-called "Paris" system-in the worldwide species Drosophila simulans has elicited the selection of drive-resistant Y chromosomes. Here, we investigate the evolutionary history of 386 Y chromosomes originating from 29 population samples collected over a period of 20 years, showing a wide continuum of phenotypes when tested against the Paris distorters, from high sensitivity to complete resistance (males sire ∼95% to ∼40% female progeny). Analyzing around 13 kb of Y-linked gene sequences in a representative subset of nine Y chromosomes, we identified only three polymorphic sites resulting in three haplotypes. Remarkably, one of the haplotypes is associated with resistance. This haplotype is fixed in all samples from Sub-Saharan Africa, the region of origin of the drivers. Exceptionally, with the spread of the drivers in Egypt and Morocco, we were able to record the replacement of the sensitive lineage by the resistant haplotype in real time, within only a few years. In addition, we performed in situ hybridization, using satellite DNA probes, on a subset of 21 Y chromosomes from six locations. In contrast to the low molecular polymorphism, this revealed extensive structural variation suggestive of rapid evolution, either neutral or adaptive. Moreover, our results show that intragenomic conflicts can drive astonishingly rapid replacement of Y chromosomes and suggest that the emergence of Paris segregation distorters in East Africa occurred less than half a century ago.

The fruitless gene affects female receptivity and species isolation.

Female mate rejection acts as a major selective force within species, and can serve as a reproductive barrier between species. In spite of its critical role in fitness and reproduction, surprisingly little is known about the genetic or neural basis of variation in female mate choice. Here, we identify fruitless as a gene affecting female receptivity within Drosophila melanogaster, as well as female Drosophila simulans rejection of male D. melanogaster. Of the multiple transcripts this gene produces, by far the most widely studied is the sex-specifically spliced transcript involved in the sex determination pathway. However, we find that female rejection behaviour is affected by a non-sex-specifically spliced fruitless transcript. This is the first implication of fruitless in female behaviour, and the first behavioural role identified for a fruitless non-sex-specifically spliced transcript. We found that this locus does not influence preferences via a single sensory modality, examining courtship song, antennal pheromone perception, or perception of substrate vibrations, and we conclude that fruitless influences mate choice via the integration of multiple signals or through another sensory modality.

A common suite of cellular abnormalities and spermatogenetic errors in sterile hybrid males in Drosophila.

When two species interbreed, the resulting hybrid offspring are often sterile, with the heterogametic (e.g. XY) hybrid usually being more severely affected. The prevailing theory for this pattern of sterility evokes divergent changes in separate lineages having maladaptive interactions when placed together in a hybrid individual, with recessive factors on the sex chromosome interacting with dominant factors on the autosomes. The effect of these interactions on gametogenesis should not be uniform across species pairs unless genetic divergence follows the same paths in different lineages or if a specific stage of gametogenesis is more susceptible to detrimental genetic interactions. Here, we perform a detailed cellular characterization of hybrid male sterility across three recently diverged species pairs of Drosophila. Across all three pairs, sterile hybrid sperm are alive but exhibit rapid nuclear de-condensation with age, with active, but non-differentiated, mitochondria. Surprisingly, all three sets of interspecies hybrids produce half of the number of sperm per round of spermatogenesis, with each sperm cell containing two tails. We identify non-disjunction failures during meiosis I as the likely cause. Thus, errors during meiosis I may be a general phenomenon underlying Drosophila male sterility, indicating either a heightened sensitivity of this spermatogenic stage to failure, or a basis to sterility other than the prevailing model.

Low concordance of short-term and long-term selection responses in experimental Drosophila populations.

Experimental evolution is becoming a popular approach to study the genomic selection response of evolving populations. Computer simulation studies suggest that the accuracy of the signature increases with the duration of the experiment. Since some assumptions of the computer simulations may be violated, it is important to scrutinize the influence of the experimental duration with real data. Here, we use a highly replicated Evolve and Resequence study in Drosophila simulans to compare the selection targets inferred at different time points. At each time point, approximately the same number of SNPs deviates from neutral expectations, but only 10% of the selected haplotype blocks identified from the full data set can be detected after 20 generations. Those haplotype blocks that emerge already after 20 generations differ from the others by being strongly selected at the beginning of the experiment and display a more parallel selection response. Consistent with previous computer simulations, our results demonstrate that only Evolve and Resequence experiments with a sufficient number of generations can characterize complex adaptive architectures.

Offspring sex ratios are stable across the life course in Drosophila simulans.

Within populations, adult sex ratios influence population growth and extinction risk, mating behaviours and parental care. Sex ratio adjustment can also have pronounced effects on individual fitness. Accordingly, it is important that we understand how often, and why, offspring sex ratios deviate from parity. In Drosophila melanogaster, females appear to improve their fitness by producing fewer sons when paired with older males. However, facultative sex ratio adjustment in D. melanogaster is controversial, and our understanding of how sex ratio skew affects fitness is hampered by pronounced sexual conflict in this species. Additionally, it is unclear whether maternal age or quality interacts with paternal age to influence offspring sex ratios. Here, we test whether offspring sex ratios vary as a function of maternal quality, and maternal and paternal age in Drosophila simulans, a sister species of D. melanogaster that lacks overt sexual conflict. We find that offspring sex ratios are slightly male-biased overall, but constant across the female life course, and independent of female quality, or paternal age. To really understand if, how and when females skew offspring sex ratios, we need studies linking offspring sex ratios to paternal and maternal phenotypes that are predicted to shift optimal investment in sons and daughters.

The recent invasion of natural Drosophila simulans populations by the P-element.

The P-element is one of the best understood eukaryotic transposable elements. It invaded Drosophila melanogaster populations within a few decades but was thought to be absent from close relatives, including Drosophila simulans. Five decades after the spread in D. melanogaster, we provide evidence that the P-element has also invaded D. simulans. P-elements in D. simulans appear to have been acquired recently from D. melanogaster probably via a single horizontal transfer event. Expression data indicate that the P-element is processed in the germ line of D. simulans, and genomic data show an enrichment of P-element insertions in putative origins of replication, similar to that seen in D. melanogaster. This ongoing spread of the P-element in natural populations provides a unique opportunity to understand the dynamics of transposable element spread and the associated piwi-interacting RNAs defense mechanisms.

Transcriptional polymorphism of piRNA regulatory genes underlies the mariner activity in Drosophila simulans testes.

During colonization of new areas, natural populations have to deal with changing environments, and transposable elements (TEs) can be useful "tools" in the adaptation process as they are major contributor to the structural and functional evolution of genomes. In this general context, the activity (copy number, transcriptional and excision rate) of the mariner mos1 element was estimated in 19 natural populations of D. simulans. It is shown (i) that mos1 expression is always higher and more variable in testes than in ovaries; (ii) that mos1 activity is higher in colonizing populations compared to the sub-Saharan African ones (ancestral populations); (iii) that mos1 variations in transcript levels and copy number are negatively correlated with transcriptional variations of piRNA genes, aubergine and argonaute3. Furthermore, mos1 levels of expression in testes highly contrast with the low expression patterns of ago3. These results strongly suggest that the expression polymorphism of piRNA genes could be responsible for the mos1 variations, first between male and female germlines and second, according to the status of natural populations (colonizing or not). These results provide new perspectives about TEs and piRNA genes co-evolution in Drosophila germlines.

The Drosophila simulans Y chromosome interacts with the autosomes to influence male fitness.

The Y chromosome should degenerate because it cannot recombine. However, male-limited transmission increases selection efficiency for male-benefit alleles on the Y, and therefore, Y chromosomes should contribute significantly to variation in male fitness. This means that although the Drosophila Y chromosome is small and gene-poor, Y-linked genes are vital for male fertility in Drosophila melanogaster and the Y chromosome has large male fitness effects. It is unclear whether the same pattern is seen in the closely related Drosophila simulans. We backcrossed Y chromosomes from three geographic locations into five genetic backgrounds and found strong Y and genetic background effects on male fertility. There was a significant Y-background interaction, indicating substantial epistasis between the Y and autosomal genes affecting male fertility. This supports accumulating evidence that interactions between the Y chromosome and the autosomes are key determinants of male fitness.

Landscape of standing variation for tandem duplications in Drosophila yakuba and Drosophila simulans.

We have used whole genome paired-end Illumina sequence data to identify tandem duplications in 20 isofemale lines of Drosophila yakuba and 20 isofemale lines of D. simulans and performed genome wide validation with PacBio long molecule sequencing. We identify 1,415 tandem duplications that are segregating in D. yakuba as well as 975 duplications in D. simulans, indicating greater variation in D. yakuba. Additionally, we observe high rates of secondary deletions at duplicated sites, with 8% of duplicated sites in D. simulans and 17% of sites in D. yakuba modified with deletions. These secondary deletions are consistent with the action of the large loop mismatch repair system acting to remove polymorphic tandem duplication, resulting in rapid dynamics of gain and loss in duplicated alleles and a richer substrate of genetic novelty than has been previously reported. Most duplications are present in only single strains, suggesting that deleterious impacts are common. Drosophila simulans shows larger numbers of whole gene duplications in comparison to larger proportions of gene fragments in D. yakuba. Drosophila simulans displays an excess of high-frequency variants on the X chromosome, consistent with adaptive evolution through duplications on the D. simulans X or demographic forces driving duplicates to high frequency. We identify 78 chimeric genes in D. yakuba and 38 chimeric genes in D. simulans, as well as 143 cases of recruited noncoding sequence in D. yakuba and 96 in D. simulans, in agreement with rates of chimeric gene origination in D. melanogaster. Together, these results suggest that tandem duplications often result in complex variation beyond whole gene duplications that offers a rich substrate of standing variation that is likely to contribute both to detrimental phenotypes and disease, as well as to adaptive evolutionary change.

Phenotypic plasticity is not affected by experimental evolution in constant, predictable or unpredictable fluctuating thermal environments.

The selective past of populations is presumed to affect the levels of phenotypic plasticity. Experimental evolution at constant temperatures is generally expected to lead to a decreased level of plasticity due to presumed costs associated with phenotypic plasticity when not needed. In this study, we investigated the effect of experimental evolution in constant, predictable and unpredictable daily fluctuating temperature regimes on the levels of phenotype plasticity in several life history and stress resistance traits in Drosophila simulans. Contrary to the expectation, evolution in the different regimes did not affect the levels of plasticity in any of the traits investigated even though the populations from the different thermal regimes had evolved different stress resistance and fitness trait means. Although costs associated with phenotypic plasticity are known, our results suggest that the maintenance of phenotypic plasticity might come at low and negligible costs, and thus, the potential of phenotypic plasticity to evolve in populations exposed to different environmental conditions might be limited.

Niche construction initiates the evolution of mutualistic interactions.

Niche construction theory explains how organisms' niche modifications may feed back to affect their evolutionary trajectories. In theory, the evolution of other species accessing the same modified niche may also be affected. We propose that this niche construction may be a general mechanism driving the evolution of mutualisms. Drosophilid flies benefit from accessing yeast-infested fruits, but the consequences of this interaction for yeasts are unknown. We reveal high levels of variation among strains of Saccharomyces cerevisiae in their ability to modify fruits and attract Drosophila simulans. More attractive yeasts are dispersed more frequently, both in the lab and in the field, and flies associated with more attractive yeasts have higher fecundity. Although there may be multiple natural yeast and fly species interactions, our controlled assays in the lab and field provide evidence of a mutualistic interaction, facilitated by the yeast's niche modification.