Correlated evolution of male and female reproductive traits drive a cascading effect of reinforcement in Drosophila yakuba

Synthesis and Scope of the Role of Postmating Prezygotic Isolation in Speciation

How barriers to gene flow arise and are maintained are key questions in evolutionary biology. Speciation research has mainly focused on barriers that occur either before mating or after zygote formation. In comparison, postmating prezygotic (PMPZ) isolation-a barrier that acts after gamete release but before zygote formation-is less frequently investigated but may hold a unique role in generating biodiversity. Here we discuss the distinctive features of PMPZ isolation, including the primary drivers and molecular mechanisms underpinning PMPZ isolation. We then present the first comprehensive survey of PMPZ isolation research, revealing that it is a widespread form of prezygotic isolation across eukaryotes. The survey also exposes obstacles in studying PMPZ isolation, in part attributable to the challenges involved in directly measuring PMPZ isolation and uncovering its causal mechanisms. Finally, we identify outstanding knowledge gaps and provide recommendations for improving future research on PMPZ isolation. This will allow us to better understand the nature of this often-neglected reproductive barrier and its contribution to speciation.

Rapid host switching of Wolbachia and even more rapid turnover of their phages and incompatibility-causing loci

About half of all insect species carry maternally inherited Wolbachia alphaproteobacteria, making Wolbachia the most common endosymbionts known in nature. Often Wolbachia spread to high frequencies within populations due to cytoplasmic incompatibility (CI), a Wolbachia-induced sperm modification caused by prophage-associated genes (cifs) that kill embryos without Wolbachia. Several Wolbachia variants also block viruses, including wMel from Drosophila melanogaster when transinfected into the mosquito Aedes aegypti. CI enables the establishment and stable maintenance of pathogen-blocking wMel in natural Ae. aegypti populations. These transinfections are reducing dengue disease incidence on multiple continents. While it has long been known that closely related Wolbachia occupy distantly related hosts, the timing of Wolbachia host switching and molecular evolution has not been widely quantified. We provide a new, conservative calibration for Wolbachia chronograms based on examples of co-divergence of Wolbachia and their insect hosts. Synthesizing publicly available and new genomic data, we use our calibration to demonstrate that wMel-like variants separated by only about 370,000 years have naturally colonized holometabolous dipteran and hymenopteran insects that diverged approximately 350 million years ago. Data from Wolbachia variants closely related to those currently dominant in D. melanogaster and D. simulans illustrate that cifs are rapidly acquired and lost among Wolbachia genomes, on a time scale of 104-105 years. This turnover occurs with and without the Wovirus prophages that contain them, with closely related cifs found in distantly related phages and distantly related cifs found in closely related phages. We present evidence for purifying selection on CI rescue function and on particular Cif protein domains. Our results quantify the tempo and mode of rapid host switching and horizontal gene transfer that underlie the spread and diversity of Wolbachia sampled from diverse host species. The wMel variants we highlight from hosts in different climates may offer new options for broadening Wolbachia-based biocontrol of diseases and pests.

The population genetics of speciation by cascade reinforcement

Species interactions drive diverse evolutionary outcomes. Speciation by cascade reinforcement represents one example of how species interactions can contribute to the proliferation of species. This process occurs when the divergence of mating traits in response to selection against interspecific hybridization incidentally leads to reproductive isolation among populations of the same species. Here, we investigated the population genetic outcomes of cascade reinforcement in North American chorus frogs (Hylidae: Pseudacris). Specifically, we estimated the frequency of hybridization among three taxa, assessed genetic structure within the focal species, P. feriarum, and ascertained the directionality of gene flow within P. feriarum across replicated contact zones via coalescent modeling. Through field observations and preliminary experimental crosses, we assessed whether hybridization is possible under natural and laboratory conditions. We found that hybridization occurs among P. feriarum and two conspecifics at a low rate in multiple contact zones, and that gene flow within the former species is unidirectional from allopatry into sympatry with these other species in three of four contact zones studied. We found evidence of substantial genetic structuring within P. feriarum including a divergent western allopatric cluster, a behaviorally-distinct sympatric South Carolina cluster, and several genetically-overlapping clusters from the remainder of the distribution. Furthermore, we found sub-structuring between reinforced and nonreinforced populations in the two most intensely-sampled contact zones. Our literature review indicated that P. feriarum hybridizes with at least five heterospecifics at the periphery of its range providing a mechanism for further intraspecific diversification. This work strengthens the evidence for cascade reinforcement in this clade, revealing the geographic and genetic landscape upon which this process can contribute to the proliferation of species.

Non-parallel morphological divergence following colonization of a new host plant

Adaptation to new ecological niches is known to spur population diversification and may lead to speciation if gene flow is ceased. While adaptation to the same ecological niche is expected to be parallel, it is more difficult to predict whether selection against maladaptive hybridization in secondary sympatry results in parallel divergence also in traits that are not directly related to the ecological niches. Such parallelisms in response to selection for reproductive isolation can be identified through estimating parallelism in reproductive character displacement across different zones of secondary contact. Here, we use a host shift in the phytophagous peacock fly Tephritis conura, with both host races represented in two geographically separate areas East and West of the Baltic Sea to investigate convergence in morphological adaptations. We asked (i) if there are consistent morphological adaptations to a host plant shift and (ii) if the response to secondary sympatry with the alternate host race is parallel across contact zones. We found surprisingly low and variable, albeit significant, divergence between host races. Only one trait, the length of the female ovipositor, which serves an important function in the interaction with the hosts, was consistently different between host races. Instead, co-existence with the other host race significantly affected the degree of morphological divergence, but the divergence was largely driven by different traits in different contact zones. Thus, local stochastic fixation or reinforcement could generate trait divergence, and additional evidence is needed to conclude whether divergence is locally adaptive.

Pure species discriminate against hybrids in the Drosophila melanogaster species subgroup

Introgression, the exchange of alleles between species, is a common event in nature. This transfer of alleles between species must happen through fertile hybrids. Characterizing the traits that cause defects in hybrids illuminates how and when gene flow is expected to occur. Inviability and sterility are extreme examples of fitness reductions but are not the only type of defects in hybrids. Some traits specific to hybrids are more subtle but are important to determine their fitness. In this report, we study whether F1 hybrids between two species pairs of Drosophila are as attractive as the parental species. We find that in both species pairs, the sexual attractiveness of the F1 hybrids is reduced and that pure species discriminate strongly against them. We also find that the cuticular hydrocarbon (CHC) profile of the female hybrids is intermediate between the parental species. Perfuming experiments show that modifying the CHC profile of the female hybrids to resemble pure species improves their chances of mating. Our results show that behavioral discrimination against hybrids might be an important component of the persistence of species that can hybridize.

A morphological trait involved in reproductive isolation between Drosophila sister species is sensitive to temperature

Male genitalia are usually extremely divergent between closely related species, but relatively constant within one species. Here we examine the effect of temperature on the shape of the ventral branches, a male genital structure involved in reproductive isolation, in the sister species Drosophila santomea and Drosophila yakuba. We designed a semi-automatic measurement machine learning pipeline that can reliably identify curvatures and landmarks based on manually digitized contours of the ventral branches. With this method, we observed that temperature does not affect ventral branches in D. yakuba but that in D. santomea ventral branches tend to morph into a D. yakuba-like shape at lower temperature. We found that male genitalia structures involved in reproductive isolation can be relatively variable within one species and can resemble the shape of closely related species' genitalia through plasticity to temperature. Our results suggest that reproductive isolation mechanisms can be dependent on the environmental context.

Heterospecific mating interactions as an interface between ecology and evolution

Reproductive interference (costly interspecific sexual interactions) is well-understood to promote divergence in mating-relevant traits (i.e. reproductive character displacement: RCD), but it can also reduce population growth, eventually leading to local extinction of one of the species. The ecological and evolutionary processes driven by reproductive interference can interact with each other. These interactions are likely to influence whether the outcome is coexistence or extinction, but remain little studied. In this paper, we first develop an eco-evolutionary perspective on reproductive interference by integrating ecological and evolutionary processes in a common framework. We also present a simple model to demonstrate the eco-evolutionary dynamics of reproductive interference. We then identify a number of factors that are likely to influence the relative likelihoods of extinction or RCD. We discuss particularly relevant factors by classifying them into four categories: the nature of the traits responding to selection, the mechanisms determining the expression of these traits, mechanisms of reproductive interference and the ecological background. We highlight previously underappreciated ways in which these factors may influence the relative likelihoods of RCD and local extinction. By doing so, we also identify questions and future directions that will increase our holistic understanding of the outcomes of reproductive interference.

Environmental and Genetic Contributions to Imperfect wMel-Like Wolbachia Transmission and Frequency Variation

Maternally transmitted Wolbachia bacteria infect about half of all insect species. They usually show imperfect maternal transmission and often produce cytoplasmic incompatibility (CI). Irrespective of CI, Wolbachia frequencies tend to increase when rare only if they benefit host fitness. Several Wolbachia, including wMel that infects Drosophila melanogaster, cause weak or no CI and persist at intermediate frequencies. On the island of São Tomé off West Africa, the frequencies of wMel-like Wolbachia infecting Drosophila yakuba (wYak) and Drosophila santomea (wSan) fluctuate, and the contributions of imperfect maternal transmission, fitness effects, and CI to these fluctuations are unknown. We demonstrate spatial variation in wYak frequency and transmission on São Tomé. Concurrent field estimates of imperfect maternal transmission do not predict spatial variation in wYak frequencies, which are highest at high altitudes where maternal transmission is the most imperfect. Genomic and genetic analyses provide little support for D. yakuba effects on wYak transmission. Instead, rearing at cool temperatures reduces wYak titer and increases imperfect transmission to levels observed on São Tomé. Using mathematical models of Wolbachia frequency dynamics and equilibria, we infer that temporally variable imperfect transmission or spatially variable effects on host fitness and reproduction are required to explain wYak frequencies. In contrast, spatially stable wSan frequencies are plausibly explained by imperfect transmission, modest fitness effects, and weak CI. Our results provide insight into causes of wMel-like frequency variation in divergent hosts. Understanding this variation is crucial to explain Wolbachia spread and to improve wMel biocontrol of human disease in transinfected mosquito systems.

The importance of intrinsic postzygotic barriers throughout the speciation process

Intrinsic postzygotic barriers can play an important and multifaceted role in speciation, but their contribution is often thought to be reserved to the final stages of the speciation process. Here, we review how intrinsic postzygotic barriers can contribute to speciation, and how this role may change through time. We outline three major contributions of intrinsic postzygotic barriers to speciation. (i) reduction of gene flow: intrinsic postzygotic barriers can effectively reduce gene exchange between sympatric species pairs. We discuss the factors that influence how effective incompatibilities are in limiting gene flow. (ii) early onset of species boundaries via rapid evolution: intrinsic postzygotic barriers can evolve between recently diverged populations or incipient species, thereby influencing speciation relatively early in the process. We discuss why the early origination of incompatibilities is expected under some biological models, and detail how other (and often less obvious) incompatibilities may also serve as important barriers early on in speciation. (iii) reinforcement: intrinsic postzygotic barriers can promote the evolution of subsequent reproductive isolation through processes such as reinforcement, even between relatively recently diverged species pairs. We incorporate classic and recent empirical and theoretical work to explore these three facets of intrinsic postzygotic barriers, and provide our thoughts on recent challenges and areas in the field in which progress can be made. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

To accept or reject heterospecific mates: behavioural decisions underlying premating isolation

Premating isolation in animals involves decision-making processes that affect whether individuals accept or reject heterospecific mates. An integrative understanding of the behavioural processes underlying heterospecific acceptance can clarify the conditions under which premating isolation evolves. As an illustration, we review how Reeve's (Reeve 1989 Am. Nat. 133, 407-435. (doi:10.1086/284926)) acceptance threshold model can help make sense of patterns of premating isolation in nature. This model derives a threshold trait value for acceptance for rejection of recipients of an action (e.g. mating) based on the fitness consequences of these decisions. We show that the maintenance of partial reproductive isolation can be an outcome of optimal acceptance thresholds, even in the face of reinforcement. We also use this model to clarify how the composition of multispecies communities can shape premating isolation. The acceptance threshold model can also be viewed as the behavioural underpinning of reproductive character displacement and cascading reinforcement. Finally, we highlight potential limitations of the acceptance threshold model with respect to investigating the role of sexual selection in speciation, and we propose that integration of behavioural models in speciation research will help us gain a full picture of the mechanisms underlying premating isolation. This article is part of the theme issue 'Signal detection theory in recognition systems: from evolving models to experimental tests'.

Rapid and Predictable Evolution of Admixed Populations Between Two Drosophila Species Pairs

The consequences of hybridization are varied, ranging from the origin of new lineages, introgression of some genes between species, to the extinction of one of the hybridizing species. We generated replicate admixed populations between two pairs of sister species of Drosophila: D. simulans and D. mauritiana; and D. yakuba and D. santomea Each pair consisted of a continental species and an island endemic. The admixed populations were maintained by random mating in discrete generations for over 20 generations. We assessed morphological, behavioral, and fitness-related traits from each replicate population periodically, and sequenced genomic DNA from the populations at generation 20. For both pairs of species, species-specific traits and their genomes regressed to those of the continental species. A few alleles from the island species persisted, but they tended to be proportionally rare among all sites in the genome and were rarely fixed within the populations. This paucity of alleles from the island species was particularly pronounced on the X-chromosome. These results indicate that nearly all foreign genes were quickly eliminated after hybridization and that selection against the minor species genome might be similar across experimental replicates.

The Past and Future of Experimental Speciation

Speciation is the result of evolutionary processes that generate barriers to gene flow between populations, facilitating reproductive isolation. Speciation is typically studied via theoretical models and snapshot tests in natural populations. Experimental speciation enables real-time direct tests of speciation theory and has been long touted as a critical complement to other approaches. We argue that, despite its promise to elucidate the evolution of reproductive isolation, experimental speciation has been underutilised and lags behind other contributions to speciation research. We review recent experiments and outline a framework for how experimental speciation can be implemented to address current outstanding questions that are otherwise challenging to answer. Greater uptake of this approach is necessary to rapidly advance understanding of speciation.

Wolbachia Acquisition by Drosophila yakuba-Clade Hosts and Transfer of Incompatibility Loci Between Distantly Related Wolbachia

Maternally transmitted Wolbachia infect about half of insect species, yet the predominant mode(s) of Wolbachia acquisition remains uncertain. Species-specific associations could be old, with Wolbachia and hosts codiversifying (i.e., cladogenic acquisition), or relatively young and acquired by horizontal transfer or introgression. The three Drosophila yakuba-clade hosts [(D. santomea, D. yakuba) D. teissieri] diverged ∼3 MYA and currently hybridize on the West African islands Bioko and São Tomé. Each species is polymorphic for nearly identical Wolbachia that cause weak cytoplasmic incompatibility (CI)-reduced egg hatch when uninfected females mate with infected males. D. yakuba-clade Wolbachia are closely related to wMel, globally polymorphic in D. melanogaster We use draft Wolbachia and mitochondrial genomes to demonstrate that D. yakuba-clade phylogenies for Wolbachia and mitochondria tend to follow host nuclear phylogenies. However, roughly half of D. santomea individuals, sampled both inside and outside of the São Tomé hybrid zone, have introgressed D. yakuba mitochondria. Both mitochondria and Wolbachia possess far more recent common ancestors than the bulk of the host nuclear genomes, precluding cladogenic Wolbachia acquisition. General concordance of Wolbachia and mitochondrial phylogenies suggests that horizontal transmission is rare, but varying relative rates of molecular divergence complicate chronogram-based statistical tests. Loci that cause CI in wMel are disrupted in D. yakuba-clade Wolbachia; but a second set of loci predicted to cause CI are located in the same WO prophage region. These alternative CI loci seem to have been acquired horizontally from distantly related Wolbachia, with transfer mediated by flanking Wolbachia-specific ISWpi1 transposons.

sppIDer: A Species Identification Tool to Investigate Hybrid Genomes with High-Throughput Sequencing

The genomics era has expanded our knowledge about the diversity of the living world, yet harnessing high-throughput sequencing data to investigate alternative evolutionary trajectories, such as hybridization, is still challenging. Here we present sppIDer, a pipeline for the characterization of interspecies hybrids and pure species, that illuminates the complete composition of genomes. sppIDer maps short-read sequencing data to a combination genome built from reference genomes of several species of interest and assesses the genomic contribution and relative ploidy of each parental species, producing a series of colorful graphical outputs ready for publication. As a proof-of-concept, we use the genus Saccharomyces to detect and visualize both interspecies hybrids and pure strains, even with missing parental reference genomes. Through simulation, we show that sppIDer is robust to variable reference genome qualities and performs well with low-coverage data. We further demonstrate the power of this approach in plants, animals, and other fungi. sppIDer is robust to many different inputs and provides visually intuitive insight into genome composition that enables the rapid identification of species and their interspecies hybrids. sppIDer exists as a Docker image, which is a reusable, reproducible, transparent, and simple-to-run package that automates the pipeline and installation of the required dependencies (https://github.com/GLBRC/sppIDer; last accessed September 6, 2018).

Aggression and discrimination among closely versus distantly related species of Drosophila

Fighting between different species is widespread in the animal kingdom, yet this phenomenon has been relatively understudied in the field of aggression research. Particularly lacking are studies that test the effect of genetic distance, or relatedness, on aggressive behaviour between species. Here we characterized male-male aggression within and between species of fruit flies across the Drosophila phylogeny. We show that male Drosophila discriminate between conspecifics and heterospecifics and show a bias for the target of aggression that depends on the genetic relatedness of opponent males. Specifically, males of closely related species treated conspecifics and heterospecifics equally, whereas males of distantly related species were overwhelmingly aggressive towards conspecifics. To our knowledge, this is the first study to quantify aggression between Drosophila species and to establish a behavioural bias for aggression against conspecifics versus heterospecifics. Our results suggest that future study of heterospecific aggression behaviour in Drosophila is warranted to investigate the degree to which these trends in aggression among species extend to broader behavioural, ecological and evolutionary contexts.

Strength of sexual and postmating prezygotic barriers varies between sympatric populations with different histories and species abundances

The impact of different reproductive barriers on species or population isolation may vary in different stages of speciation depending on evolutionary forces acting within species and through species' interactions. Genetic incompatibilities between interacting species are expected to reinforce prezygotic barriers in sympatric populations and lead to cascade reinforcement between conspecific populations living within and outside the areas of sympatry. We tested these predictions and studied whether and how the strength and target of reinforcement between Drosophila montana and Drosophila flavomontana vary between sympatric populations with different histories and species abundances. All barriers between D. montana females and D. flavomontana males were nearly complete, while in the reciprocal cross strong postzygotic isolation was accompanied by prezygotic barriers whose strength varied according to population composition. Sexual isolation between D. flavomontana females and D. montana males was increased in long-established sympatric populations, where D. flavomontana is abundant, while postmating prezygotic (PMPZ) barriers were stronger in populations where this species is a new invader and still rare and where female discrimination against heterospecific males was lower. Strengthening of sexual and PMPZ barriers in this cross also induced cascade reinforcement of respective barriers between D. flavomontana populations, which is a classic signature of reinforcement process.

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.

Genomic Signatures of Reinforcement

Reinforcement is the process by which selection against hybridization increases reproductive isolation between taxa. Much research has focused on demonstrating the existence of reinforcement, yet relatively little is known about the genetic basis of reinforcement or the evolutionary conditions under which reinforcement can occur. Inspired by reinforcement's characteristic phenotypic pattern of reproductive trait divergence in sympatry but not in allopatry, we discuss whether reinforcement also leaves a distinct genomic pattern. First, we describe three patterns of genetic variation we expect as a consequence of reinforcement. Then, we discuss a set of alternative processes and complicating factors that may make the identification of reinforcement at the genomic level difficult. Finally, we consider how genomic analyses can be leveraged to inform if and to what extent reinforcement evolved in the face of gene flow between sympatric lineages and between allopatric and sympatric populations of the same lineage. Our major goals are to understand if genome scans for particular patterns of genetic variation could identify reinforcement, isolate the genetic basis of reinforcement, or infer the conditions under which reinforcement evolved.

Moving Speciation Genetics Forward: Modern Techniques Build on Foundational Studies in Drosophila

The question of how new species evolve has been examined at every level, from macroevolutionary patterns of diversification to molecular population genetic analyses of specific genomic regions between species pairs. Drosophila has been at the center of many of these research efforts. Though our understanding of the speciation process has grown considerably over the past few decades, very few genes have been identified that contribute to barriers to reproduction. The development of advanced molecular genetic and genomic methods provides promising avenues for the rapid discovery of more genes that contribute to speciation, particularly those involving prezygotic isolation. The continued expansion of tools and resources, especially for species other than Drosophila melanogaster, will be most effective when coupled with comparative approaches that reveal the genetic basis of reproductive isolation across a range of divergence times. Future research programs in Drosophila have high potential to answer long-standing questions in speciation. These include identifying the selective forces that contribute to divergence between populations and the genetic basis of traits that cause reproductive isolation. The latter can be expanded upon to understand how the genetic basis of reproductive isolation changes over time and whether certain pathways and genes are more commonly involved.

Male-driven reproductive and agonistic character displacement in darters and its implications for speciation in allopatry

Selection against hybridization can cause mating traits to diverge between species in sympatry via reproductive character displacement (RCD). Additionally, selection against interspecific fighting can cause aggressive traits to diverge between sympatric species via agonistic character displacement (ACD). By directly affecting conspecific recognition traits, RCD and ACD between species can also incidentally cause divergence in mating and fighting traits among populations within a species [termed cascade RCD (CRCD) and cascade ACD]. Here, we demonstrate patterns consistent with male-driven RCD and ACD in 2 groups of darters (orangethroat darter clade Ceasia and rainbow darter Etheostoma caeruleum). In both groups, males that occur in sympatry (between Ceasia and E. caeruleum) have higher levels of preference for mating and fighting with conspecifics over heterospecifics than do males from allopatry. This is consistent with RCD and ACD. We also found patterns consistent with CRCD and cascade ACD among species of Ceasia. Ceasia males that are sympatric to E. caeruleum (but allopatric to one another) also have heightened preferences for mating and fighting with conspecific versus heterospecific Ceasia. In contrast, Ceasia males that are allopatric to E. caeruleum readily mate and fight with heterospecific Ceasia. We suggest that RCD and ACD between Ceasia and E. caeruleum has incidentally led to divergence in mating and fighting traits among Ceasia species. This study is unique in that male preferences evolve via both RCD (male preference for conspecific females) and ACD (male preference to fight conspecific males) which leads to subsequent divergence among allopatric lineages.

The Rate of Evolution of Postmating-Prezygotic Reproductive Isolation in Drosophila

Reproductive isolation is an intrinsic aspect of species formation. For that reason, the identification of the precise isolating traits, and the rates at which they evolve, is crucial to understanding how species originate and persist. Previous work has measured the rates of evolution of prezygotic and postzygotic barriers to gene flow, yet no systematic analysis has studied the rates of evolution of postmating-prezygotic (PMPZ) barriers. We measured the magnitude of two barriers to gene flow that act after mating occurs but before fertilization. We also measured the magnitude of a premating barrier (female mating rate in nonchoice experiments) and two postzygotic barriers (hybrid inviability and hybrid sterility) for all pairwise crosses of all nine known extant species within the melanogaster subgroup. Our results indicate that PMPZ isolation evolves faster than hybrid inviability but slower than premating isolation. Next, we partition postzygotic isolation into different components and find that, as expected, hybrid sterility evolves faster than hybrid inviability. These results lend support for the hypothesis that, in Drosophila, reproductive isolation mechanisms (RIMs) that act early in reproduction (or in development) tend to evolve faster than those that act later in the reproductive cycle. Finally, we tested whether there was evidence for reinforcing selection at any RIM. We found no evidence for generalized evolution of reproductive isolation via reinforcement which indicates that there is no pervasive evidence of this evolutionary process. Our results indicate that PMPZ RIMs might have important evolutionary consequences in initiating speciation and in the persistence of new species.

Fine scale mapping of genomic introgressions within the Drosophila yakuba clade

The process of speciation involves populations diverging over time until they are genetically and reproductively isolated. Hybridization between nascent species was long thought to directly oppose speciation. However, the amount of interspecific genetic exchange (introgression) mediated by hybridization remains largely unknown, although recent progress in genome sequencing has made measuring introgression more tractable. A natural place to look for individuals with admixed ancestry (indicative of introgression) is in regions where species co-occur. In west Africa, D. santomea and D. yakuba hybridize on the island of São Tomé, while D. yakuba and D. teissieri hybridize on the nearby island of Bioko. In this report, we quantify the genomic extent of introgression between the three species of the Drosophila yakuba clade (D. yakuba, D. santomea), D. teissieri). We sequenced the genomes of 86 individuals from all three species. We also developed and applied a new statistical framework, using a hidden Markov approach, to identify introgression. We found that introgression has occurred between both species pairs but most introgressed segments are small (on the order of a few kilobases). After ruling out the retention of ancestral polymorphism as an explanation for these similar regions, we find that the sizes of introgressed haplotypes indicate that genetic exchange is not recent (>1,000 generations ago). We additionally show that in both cases, introgression was rarer on X chromosomes than on autosomes which is consistent with sex chromosomes playing a large role in reproductive isolation. Even though the two species pairs have stable contemporary hybrid zones, providing the opportunity for ongoing gene flow, our results indicate that genetic exchange between these species is currently rare.

The Effect of Temperature on Drosophila Hybrid Fitness

Mechanisms of reproductive isolation inhibit gene flow between species and can be broadly sorted into two categories: prezygotic and postzygotic. While comparative studies suggest that prezygotic barriers tend to evolve first, postzygotic barriers are crucial for maintaining species boundaries and impeding gene flow that might otherwise cause incipient species to merge. Most, but not all, postzygotic barriers result from genetic incompatibilities between two or more loci from different species, and occur due to divergent evolution in allopatry. Hybrid defects result from improper allelic interactions between these loci. While some postzygotic barriers are environmentally-independent, the magnitude of others has been shown to vary in penetrance depending on environmental factors. We crossed Drosophila melanogaster mutants to two other species, D. simulans and D. santomea, and collected fitness data of the hybrids at two different temperatures. Our goal was to examine the effect of temperature on recessive incompatibility alleles in their genomes. We found that temperature has a stronger effect on the penetrance of recessive incompatibility alleles in the D. simulans genome than on those in the D. santomea genome. These results suggest that the penetrance of hybrid incompatibilities can be strongly affected by environmental context, and that the magnitude of such gene-by-environment interactions can be contingent on the genotype of the hybrid.

Wolbachia in the Drosophila yakuba Complex: Pervasive Frequency Variation and Weak Cytoplasmic Incompatibility, but No Apparent Effect on Reproductive Isolation

Three hybridizing species-the clade [(Drosophila yakuba, D. santomea), D. teissieri]-comprise the yakuba complex in the D. melanogaster subgroup. Their ranges overlap on Bioko and São Tomé, islands off west Africa. All three species are infected with Wolbachia-maternally inherited, endosymbiotic bacteria, best known for manipulating host reproduction to favor infected females. Previous analyses reported no cytoplasmic incompatibility (CI) in these species. However, we discovered that Wolbachia from each species cause intraspecific and interspecific CI. In D teissieri, analyses of F1 and backcross genotypes show that both host genotype and Wolbachia variation modulate CI intensity. Wolbachia-infected females seem largely protected from intraspecific and interspecific CI, irrespective of Wolbachia and host genotypes. Wolbachia do not affect host mating behavior or female fecundity, within or between species. The latter suggests little apparent effect of Wolbachia on premating or gametic reproductive isolation (RI) between host species. In nature, Wolbachia frequencies varied spatially for D. yakuba in 2009, with 76% (N = 155) infected on São Tomé, and only 3% (N = 36) infected on Bioko; frequencies also varied temporally in D. yakuba and D. santomea on São Tomé between 2009 and 2015. These temporal frequency fluctuations could generate asymmetries in interspecific mating success, and contribute to postzygotic RI. However, the fluctuations in Wolbachia frequencies that we observe also suggest that asymmetries are unlikely to persist. Finally, we address theoretical questions that our empirical findings raise about Wolbachia persistence when conditions fluctuate, and about the stable coexistence of Wolbachia and host variants that modulate Wolbachia effects.