An experimental test of the evolutionary consequences of sympatry in Drosophila subquinaria.

Prezygotic isolation is often stronger between sympatric as opposed to allopatric taxa, but the underlying cause can be difficult to infer from comparative studies alone. Experimental evolution, where evolutionary responses to treatments manipulating the presence/absence of heterospecific individuals are tracked, can provide a powerful complementary approach. We used experimental evolution to investigate a naturally occurring pattern of reproductive character displacement in the mushroom-feeding fly, Drosophila subquinaria. In nature, female D. subquinaria from populations sympatric with the closely related Drosophila recens discriminate more strongly against heterospecific males than do females from allopatric populations. Starting with 16 replicate allopatric populations of D. subquinaria, we manipulated the presence/absence of D. recens during mating (experimental sympatry vs. control) and, when present, we allowed hybrids to live or kill them each generation. Across 12 generations, heterospecific offspring production from no-choice mating trials between D. subquinaria females and D. recens males declined in both experimental sympatry treatments relative to the control, suggesting increased sexual isolation. Male cuticular hydrocarbon profiles also evolved, but only in the hybrids killed treatment. Our results strongly imply that the existing reproductive character displacement in wild D. subquinaria populations was an evolutionary response to selection arising from secondary contact with D. recens.

Investigating the phylogenetic history of toxin tolerance in mushroom-feeding Drosophila.

Understanding how and when key novel adaptations evolved is a central goal of evolutionary biology. Within the immigrans-tripunctata radiation of Drosophila, many mushroom-feeding species are tolerant of host toxins, such as cyclopeptides, that are lethal to nearly all other eukaryotes. In this study, we used phylogenetic and functional approaches to investigate the evolution of cyclopeptide tolerance in the immigrans-tripunctata radiation of Drosophila. First, we inferred the evolutionary relationships among 48 species in this radiation using 978 single copy orthologs. Our results resolved previous incongruities within species groups across the phylogeny. Second, we expanded on previous studies of toxin tolerance by assaying 16 of these species for tolerance to α-amanitin and found that six of them could develop on diet with toxin. Finally, we asked how α-amanitin tolerance might have evolved across the immigrans-tripunctata radiation, and inferred that toxin tolerance was ancestral in mushroom-feeding Drosophila and subsequently lost multiple times. Our findings expand our understanding of toxin tolerance across the immigrans-tripunctata radiation and emphasize the uniqueness of toxin tolerance in this adaptive radiation and the complexity of biochemical adaptations.

Investigating the phylogenetic history of toxin tolerance in mushroom-feeding Drosophila.

Understanding how and when key novel adaptations evolved is a central goal of evolutionary biology. Within the immigrans-tripunctata radiation of Drosophila , many mushroom-feeding species are tolerant of host toxins, such as cyclopeptides, that are lethal to nearly all other eukaryotes. In this study, we used phylogenetic and functional approaches to investigate the evolution of cyclopeptide tolerance in the immigrans-tripunctata radiation of Drosophila . We first inferred the evolutionary relationships among 48 species in this radiation using 978 single copy orthologs. Our results resolved previous incongruities within species groups across the phylogeny. Second, we expanded on previous studies of toxin tolerance by assaying 16 of these species for tolerance to α-amanitin and found that six of these species could develop on diet with toxin. Third, we examined fly development on a diet containing a natural mix of toxins extracted from the Death Cap Amanita phalloides mushroom. Both tolerant and susceptible species developed on diet with this mix, though tolerant species survived at significantly higher concentrations. Finally, we asked how cyclopeptide tolerance might have evolved across the immigrans-tripunctata radiation and inferred that toxin tolerance was ancestral and subsequently lost multiple times. Our results suggest the evolutionary history of cyclopeptide tolerance is complex, and simply describing this trait as present or absent does not fully capture the occurrence or impact on this adaptive radiation. More broadly, the evolution of novelty can be more complex than previously thought, and that accurate descriptions of such novelties are critical in studies examining their evolution.

Drosophila Evolution over Space and Time (DEST): A New Population Genomics Resource.

Drosophila melanogaster is a leading model in population genetics and genomics, and a growing number of whole-genome data sets from natural populations of this species have been published over the last years. A major challenge is the integration of disparate data sets, often generated using different sequencing technologies and bioinformatic pipelines, which hampers our ability to address questions about the evolution of this species. Here we address these issues by developing a bioinformatics pipeline that maps pooled sequencing (Pool-Seq) reads from D. melanogaster to a hologenome consisting of fly and symbiont genomes and estimates allele frequencies using either a heuristic (PoolSNP) or a probabilistic variant caller (SNAPE-pooled). We use this pipeline to generate the largest data repository of genomic data available for D. melanogaster to date, encompassing 271 previously published and unpublished population samples from over 100 locations in >20 countries on four continents. Several of these locations have been sampled at different seasons across multiple years. This data set, which we call Drosophila Evolution over Space and Time (DEST), is coupled with sampling and environmental metadata. A web-based genome browser and web portal provide easy access to the SNP data set. We further provide guidelines on how to use Pool-Seq data for model-based demographic inference. Our aim is to provide this scalable platform as a community resource which can be easily extended via future efforts for an even more extensive cosmopolitan data set. Our resource will enable population geneticists to analyze spatiotemporal genetic patterns and evolutionary dynamics of D. melanogaster populations in unprecedented detail.

Fitness consequences of a non-recombining sex-ratio drive chromosome can explain its prevalence in the wild.

Understanding the pleiotropic consequences of gene drive systems on host fitness is essential to predict their spread through a host population. Here, we study sex-ratio (SR) X-chromosome drive in the fly Drosophila recens, where SR causes the death of Y-bearing sperm in male carriers. SR males only sire daughters, which all carry SR, thus giving the chromosome a transmission advantage. The prevalence of the SR chromosome appears stable, suggesting pleiotropic costs. It was previously shown that females homozygous for SR are sterile, and here, we test for additional fitness costs of SR. We found that females heterozygous for SR have reduced fecundity and that male SR carriers have reduced fertility in conditions of sperm competition. We then use our fitness estimates to parametrize theoretical models of SR drive and show that the decrease in fecundity and sperm competition performance can account for the observed prevalence of SR in natural populations. In addition, we found that the expected equilibrium frequency of the SR chromosome is particularly sensitive to the degree of multiple mating and performance in sperm competition. Together, our data suggest that the mating system of the organism should be carefully considered during the development of gene drive systems.

Ongoing hybridization obscures phylogenetic relationships in the Drosophila subquinaria species complex.

Inferring evolutionary relationships among recently diverged lineages is necessary to understand how isolating barriers produce independent lineages. Here, we investigate the phylogenetic relationships between three incompletely isolated and closely related mushroom-feeding Drosophila species. These species form the Drosophila subquinaria species complex and consist of one Eurasian species (D. transversa) and two widespread North American species (D. subquinaria and D. recens) that are sympatric in central Canada. Although patterns of pre- and post-mating isolation among these species are well characterized, previous work on their phylogenetic relationships is limited and conflicting. In this study, we generated a multi-locus data set of 29 loci from across the genome sequenced in a population sample from each species, and then, we inferred species relationships and patterns of introgression. We find strong statistical support that D. subquinaria is paraphyletic, showing that samples from the geographic region sympatric with D. recens are most closely related to D. recens, whereas samples from the geographic region allopatric with D. recens are most closely related to D. transversa. We present several lines of evidence that both incomplete lineage sorting and gene flow are causing phylogenetic discordance. We suggest that ongoing gene flow primarily from D. recens into D. subquinaria in the sympatric part of their ranges causes phylogenetic uncertainty in the evolutionary history of these species. Our results highlight how population genetic data can be used to disentangle the sources of phylogenetic discordance among closely related species.

A phylogenetic examination of host use evolution in the quinaria and testacea groups of Drosophila.

Adaptive radiations provide an opportunity to examine complex evolutionary processes such as ecological specialization and speciation. While a well-resolved phylogenetic hypothesis is critical to completing such studies, the rapid rates of evolution in these groups can impede phylogenetic studies. Here we study the quinaria and testacea species groups of the immigrans-tripunctata radiation of Drosophila, which represent a recent adaptive radiation and are a developing model system for ecological genetics. We were especially interested in understanding host use evolution in these species. In order to infer a phylogenetic hypothesis for this group we sampled loci from both the nuclear genome and the mitochondrial DNA to develop a dataset of 43 protein-coding loci for these two groups along with their close relatives in the immigrans-tripunctata radiation. We used this dataset to examine their evolutionary relationships along with the evolution of feeding behavior. Our analysis recovers strong support for the monophyly of the testacea but not the quinaria group. Results from our ancestral state reconstruction analysis suggests that the ancestor of the testacea and quinaria groups exhibited mushroom-feeding. Within the quinaria group, we infer that transition to vegetative feeding occurred twice, and that this transition did not coincide with a genome-wide change in the rate of protein evolution.

A fast-evolving X-linked duplicate of importin-α2 is overexpressed in sex-ratio drive in Drosophila neotestacea.

Selfish genetic elements that manipulate gametogenesis to achieve a transmission advantage are known as meiotic drivers. Sex-ratio X chromosomes (SR) are meiotic drivers that prevent the maturation of Y-bearing sperm in male carriers to result in the production of mainly female progeny. The spread of an SR chromosome can affect host genetic diversity and genome evolution, and can even cause host extinction if it reaches sufficiently high prevalence. Meiotic drivers have evolved independently many times, though only in a few cases is the underlying genetic mechanism known. In this study we use a combination of transcriptomics and population genetics to identify widespread expression differences between the standard (ST) and sex-ratio (SR) X chromosomes of the fly Drosophila neotestacea. We found the X chromosome is enriched for differentially expressed transcripts and that many of these X-linked differentially expressed transcripts had elevated Ka /Ks values between ST and SR, indicative of potential functional differences. We identified a set of candidate transcripts, including a testis-specific, X-linked duplicate of the nuclear transport gene importin-α2 that is overexpressed in SR. We find suggestions of positive selection in the lineage leading to the duplicate and that its molecular evolutionary patterns are consistent with relaxed purifying selection in ST. As these patterns are consistent with involvement in the mechanism of drive in this species, this duplicate is a strong candidate worthy of further functional investigation. Nuclear transport may be a common target for genetic conflict, as the mechanism of the autosomal Segregation Distorter drive system in D. melanogaster involves the same pathway.

Fine-scale geographic patterns of gene flow and reproductive character displacement in Drosophila subquinaria and Drosophila recens.

When two species are incompletely isolated, strengthening premating isolation barriers in response to the production of low fitness hybrids may complete the speciation process. Here, we use the sister species Drosophila subquinaria and Drosophila recens to study the conditions under which this reinforcement of species boundaries occurs in natural populations. We first extend the region of known sympatry between these species, and then we conduct a fine-scale geographic survey of mate discrimination coupled with estimates of gene flow within and admixture between species. Within D. subquinaria, reinforcement is extremely effective: we find variation in mate discrimination both against D. recens males and against conspecific allopatric males on the scale of a few kilometres and in the face of gene flow both from conspecific populations and introgression from D. recens. In D. recens, we do not find evidence for increased mate discrimination in sympatry, even where D. recens is rare, consistent with substantial gene flow throughout the species' range. Finally, we find that introgression between species is asymmetric, with more from D. recens into D. subquinaria than vice versa. Within each species, admixture is highest in the geographic region where it is rare relative to the other species, suggesting that when hybrids are produced they are of low fitness. In sum, reinforcement within D. subquinaria is effective at maintaining species boundaries, but even when reinforcing selection is strong it may not always result in a pattern of strong reproductive character displacement due to variation in the frequency of hybridization and gene flow from neighbouring populations.

The Ecology and Evolutionary Dynamics of Meiotic Drive.

Meiotic drivers are genetic variants that selfishly manipulate the production of gametes to increase their own rate of transmission, often to the detriment of the rest of the genome and the individual that carries them. This genomic conflict potentially occurs whenever a diploid organism produces a haploid stage, and can have profound evolutionary impacts on gametogenesis, fertility, individual behaviour, mating system, population survival, and reproductive isolation. Multiple research teams are developing artificial drive systems for pest control, utilising the transmission advantage of drive to alter or exterminate target species. Here, we review current knowledge of how natural drive systems function, how drivers spread through natural populations, and the factors that limit their invasion.

Patterns of reproductive isolation in the Drosophila subquinaria complex: can reinforced premating isolation cascade to other species?

The reinforcement of premating barriers due to reduced hybrid fitness in sympatry may cause secondary sexual isolation within a species as a by-product. Consistent with this, in the fly Drosophila subquinaria, females that are sympatric with D. recens mate at very low rates not only with D. recens, but also with conspecific D. subquinaria males from allopatry. Here, we ask if these effects of reinforcement cascade more broadly to affect sexual isolation with other closely related species. We assay reproductive isolation of these species with D. transversa and find that choosy D. subquinaria females from the region sympatric with D. recens discriminate strongly against male D. transversa, whereas D. subquinaria from the allopatric region do not. This increased sexual isolation cannot be explained by natural selection to avoid mating with this species, as they are allopatric in geographic range and we do not identify any intrinsic postzygotic isolation between D. subquinaria and D. transversa. Variation in epicuticular hydrocarbons, which are used as mating signals in D. subquinaria, follow patterns of premating isolation: D. transversa and allopatric D. subquinaria are most similar to each other and differ from sympatric D. subquinaria, and those of D. recens are distinct from the other two species. We suggest that the secondary effects of reinforcement may cascade to strengthen reproductive isolation with other species that were not a target of selection. These effects may enhance the divergence that occurs in allopatry to help explain why some species are already sexually isolated upon secondary contact.

Phylogenetic Patterns of Codon Evolution in the ACTIN-DEPOLYMERIZING FACTOR/COFILIN (ADF/CFL) Gene Family.

The actin-depolymerizing factor/cofilin (ADF/CFL) gene family encodes a diverse group of relatively small proteins. Once known strictly as modulators of actin filament dynamics, recent research has demonstrated that these proteins are involved in a variety of cellular processes, from signal transduction to the cytonuclear trafficking of actin. In both plant and animal lineages, expression patterns of paralogs in the ADF/CFL gene family vary among tissue types and developmental stages. In this study we use computational approaches to investigate the evolutionary forces responsible for the diversification of the ADF/CFL gene family. Estimating the rate of non-synonymous to synonymous mutations (dN/dS) across phylogenetic lineages revealed that the majority of ADF/CFL codon positions were under strong purifying selection, with rare episodic events of accelerated protein evolution. In both plants and animals these instances of accelerated evolution were ADF/CFL subclass specific, and all of the sites under selection were located in regions of the protein that could serve in new functional roles. We suggest these sites may have been important in the functional diversification of ADF/CFL proteins.

PLANT EVOLUTION. Convergent evolution of strigolactone perception enabled host detection in parasitic plants.

Obligate parasitic plants in the Orobanchaceae germinate after sensing plant hormones, strigolactones, exuded from host roots. In Arabidopsis thaliana, the α/ß-hydrolase D14 acts as a strigolactone receptor that controls shoot branching, whereas its ancestral paralog, KAI2, mediates karrikin-specific germination responses. We observed that KAI2, but not D14, is present at higher copy numbers in parasitic species than in nonparasitic relatives. KAI2 paralogs in parasites are distributed into three phylogenetic clades. The fastest-evolving clade, KAI2d, contains the majority of KAI2 paralogs. Homology models predict that the ligand-binding pockets of KAI2d resemble D14. KAI2d transgenes confer strigolactone-specific germination responses to Arabidopsis thaliana. Thus, the KAI2 paralogs D14 and KAI2d underwent convergent evolution of strigolactone recognition, respectively enabling developmental responses to strigolactones in angiosperms and host detection in parasites.

Reproductive character displacement of female mate preferences for male cuticular hydrocarbons in Drosophila subquinaria.

Several lines of evidence implicate sexual isolation in both initiating and completing the speciation process. Although its existence is straightforward to demonstrate, understanding the evolution of sexual isolation requires identifying the underlying phenotypes responsible so that we can determine how these have diverged. Here, we study geographic variation in female mate preferences for male sexual displays in the fly Drosophila subquinaria. Female D. subquinaria that are sympatric with its sister species D. recens discriminate strongly against both D. recens and allopatric conspecific males, whereas females from allopatric populations do not. Furthermore, female mate preferences target at least in part a suite of cuticular hydrocarbons (CHCs) in males and geographic variation in CHCs mirrors the pattern of mate discrimination. In this study, we quantify female mate preferences for male CHCs from populations that span the geographic range of D. subquinaria. We find that the direction of linear sexual selection varies significantly between populations that are sympatric versus allopatric with D. recens in a pattern of reproductive character displacement. Differences in preference partially align with existing differences in CHCs and patterns of sexual isolation, although discrepancies remain that suggest the involvement of additional traits and/or more complex, nonlinear preference functions.

Asymmetrical sexual isolation but no postmating isolation between the closely related species Drosophila suboccidentalis and Drosophila occidentalis.

BACKGROUND: During the speciation process several types of isolating barriers can arise that limit gene flow between diverging populations. Studying recently isolated species can inform our understanding of how and when these barriers arise, and which barriers may be most important to limiting gene flow. Here we focus on Drosophila suboccidentalis and D. occidentalis, which are closely related mushroom-feeding species that inhabit western North America and are not known to overlap in geographic range. We investigate patterns of reproductive isolation between these species, including premating, postmating prezygotic, and postzygotic barriers to gene flow. RESULTS: Using flies that originate from a single population of each species, we find that the strength of premating sexual isolation between these species is asymmetric: while D. occidentalis females mate with D. suboccidentalis males at a reduced but moderate rate, D. suboccidentalis females discriminate strongly against mating with D. occidentalis males. Female hybrids will mate at high rates with males of either species, indicating that this discrimination has a recessive genetic basis. Hybrid males are accepted by females of both species. We do not find evidence for postmating prezygotic or postzygotic isolating barriers, as females use the sperm of heterospecific males and both male and female hybrids are fully fertile. CONCLUSIONS: Premating isolation is substantial but incomplete, and appears to be the primary form of reproductive isolation between these species. If these species do hybridize, the lack of postzygotic barriers may allow for gene flow between them. Given that these species are recently diverged and are not known to be sympatric, the level of premating isolation is relatively strong given the lack of intrinsic postzygotic isolation. Further work is necessary to characterize the geographic and genetic variation in reproductive isolating barriers, as well as to determine the factors that drive reproductive isolation and the consequences that isolating barriers as well as geographic isolation have had on patterns of gene flow between these species.

Reinforcement shapes clines in female mate discrimination in Drosophila subquinaria.

Reinforcement of species boundaries may alter mate recognition in a way that also affects patterns of mate preference among conspecific populations. In the fly Drosophila subquinaria, females sympatric with the closely related species D. recens reject mating with heterospecific males as well as with conspecific males from allopatric populations. Here, we assess geographic variation in behavioral isolation within and among populations of D. subquinaria and use cline theory to understand patterns of selection on reinforced discrimination and its consequences for sexual isolation within species. We find that selection has fixed rejection of D. recens males in sympatry, while significant genetic variation in this behavior occurs within allopatric populations. In conspecific matings sexual isolation is also asymmetric and stronger in populations that are sympatric with D. recens. The clines in behavioral discrimination within and between species are similar in shape and are maintained by strong selection in the face of gene flow, and we show that some of their genetic basis may be either shared or linked. Thus, while reinforcement can drive extremely strong phenotypic divergence, the long-term consequences for incipient speciation depend on gene flow, genetic linkage of discrimination traits, and the cost of these behaviors in allopatry.

Reproductive character displacement of epicuticular compounds and their contribution to mate choice in Drosophila subquinaria and Drosophila recens.

Interactions between species can alter selection on sexual displays used in mate choice within species. Here we study the epicuticular pheromones of two Drosophila species that overlap partially in geographic range and are incompletely reproductively isolated. Drosophila subquinaria shows a pattern of reproductive character displacement against Drosophila recens, and partial behavioral isolation between conspecific sympatric versus allopatric populations, whereas D. recens shows no such variation in mate choice. First, using manipulative perfuming experiments, we show that females use pheromones as signals for mate discrimination both between species and among populations of D. subquinaria. Second, we show that patterns of variation in epicuticular compounds, both across populations and between species, are consistent with those previously shown for mating probabilities: pheromone compositions differ between populations of D. subquinaria that are allopatric versus sympatric with D. recens, but are similar across populations of D. recens regardless of overlap with D. subquinaria. We also identify differences in pheromone composition among allopatric regions of D. subquinaria. In sum, our results suggest that epicuticular compounds are key signals used by females during mate recognition, and that these traits have diverged among D. subquinaria populations in response to reinforcing selection generated by the presence of D. recens.

Association of polyandry and sex-ratio drive prevalence in natural populations of Drosophila neotestacea.

Selfish genetic elements bias their own transmission to the next generation, even at the expense of the fitness of their carrier. Sex-ratio (SR) meiotic drive occurs when an X-chromosome causes Y-bearing sperm to die during male spermatogenesis, so that it is passed on to all of the male's offspring, which are all daughters. How SR is maintained as a stable polymorphism in the absence of genetic suppressors of drive is unknown. Here, we investigate the potential for the female remating rate to affect SR dynamics in natural populations, using the fly Drosophila neotestacea. In controlled laboratory conditions, females from populations where SR is rare mate more often than females from populations where SR is common. Furthermore, only when males mate multiply does the average fertility of SR males relative to wild-type males decrease to a level that can prevent SR from spreading. Our results suggest that differences in the female mating rate among populations may contribute to SR dynamics in the wild, and thus also affect the outcome of this intragenomic conflict. In line with this, we also present evidence of a localized population crash due to SR that may have resulted from habitat fragmentation along with a reduced mating rate.

Epicuticular compounds of Drosophila subquinaria and D. recens: identification, quantification, and their role in female mate choice.

The epicuticle of various Drosophila species consists of long-chain cuticular hydrocarbons (CHCs) and their derivatives that play a role in waterproofing and a dynamic means of chemical communication. Here, via gas chromatography and mass spectrometry, we identified and quantified the epicuticular composition of D. recens and D. subquinaria, two closely related species that show a pattern of reproductive character displacement in nature. Twenty-four compounds were identified with the most abundant, 11-cis-Vaccenyl acetate, present only in males of each species. Also exclusive to males were five tri-acylglycerides. The 18 remaining compounds were CHCs, all shared between the sexes and species. These CHCs were composed of odd carbon numbers (C29, C31, C33, and C35), with an increase in structural isomers in the C33 and C35 groups. Saturated hydrocarbons comprise only methyl-branched alkanes and were found only in the C29 and C31 groups. Alkenes were the least prevalent, with alkadienes dominating the chromatographic landscape in the longer chain lengths. Sexual dimorphism was extensive with 6/8 of the logcontrast CHCs differing significantly in relative concentration between males and females in D. recens and D. subquinaria, respectively. Males of the two species also differed significantly in relative concentration of six CHCs, while females differed in none. Female-choice mating trials revealed directional sexual selection on male CHCs in a population of each species, consistent with female mate preferences for these traits. The sexual selection vectors differed significantly in multivariate trait space, suggesting that different pheromone blends determine male attractiveness in each species.

Divergence of premating behaviors in the closely related species Drosophila subquinaria and D. recens.

Most animal species use distinctive courship patterns to choose among potential mates. Over time, the sensory signaling and preferences used during courtship can diverge among groups that are reproductively isolated. This divergence of signal traits and preferences is thought to be an important cause of behavioral isolation during the speciation process. Here, we examine the sensory modalities used in courtship by two closely related species, Drosophila subquinaria and Drosophila recens, which overlap in geographic range and are incompletely reproductively isolated. We use observational studies of courtship patterns and manipulation of male and female sensory modalities to determine the relative roles of visual, olfactory, gustatory, and auditory signals during conspecific mate choice. We find that sex-specific, species-specific, and population-specific cues are used during mate acquisition within populations of D. subquinaria and D. recens. We identify shifts in both male and female sensory modalities between species, and also between populations of D. subquinaria. Our results indicate that divergence in mating signals and preferences have occurred on a relatively short timescale within and between these species. Finally, we suggest that because olfactory cues are essential for D. subquinaria females to mate within species, they may also underlie variation in behavioral discrimination across populations and species.