Rapid Evolution of Assortative Fertilization between Recently Allopatric Species of Drosophila

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.

Evolutionary Quantitative Proteomics of Reproductive Protein Divergence in Drosophila

Reproductive traits often evolve rapidly between species. Understanding the causes and consequences of this rapid divergence requires characterization of female and male reproductive proteins and their effect on fertilization success. Species in the Drosophila virilis clade exhibit rampant interspecific reproductive incompatibilities, making them ideal for studies on diversification of reproductive proteins and their role in speciation. Importantly, the role of intraejaculate protein abundance and allocation in interspecific divergence is poorly understood. Here, we identify and quantify the transferred male ejaculate proteome using multiplexed isobaric labeling of the lower female reproductive tract before and immediately after mating using three species of the virilis group. We identified over 200 putative male ejaculate proteins, many of which show differential abundance between species, suggesting that males transfer a species-specific allocation of seminal fluid proteins during copulation. We also identified over 2000 female reproductive proteins, which contain female-specific serine-type endopeptidases that showed differential abundance between species and elevated rates of molecular evolution, similar to that of some male seminal fluid proteins. Our findings suggest that reproductive protein divergence can also manifest in terms of species-specific protein abundance patterns.

Differences in Postmating Transcriptional Responses between Conspecific and Heterospecific Matings in Drosophila

In many animal species, females undergo physiological and behavioral changes after mating. Some of these changes are driven by male-derived seminal fluid proteins and are critical for fertilization success. Unfortunately, our understanding of the molecular interplay between female and male reproductive proteins remains inadequate. Here, we analyze the postmating response in a Drosophila species that has evolved strong gametic incompatibility with its sister species; Drosophila novamexicana females produce only ∼1% fertilized eggs in crosses with Drosophila americana males, compared to ∼98% produced in within-species crosses. This incompatibility is likely caused by mismatched male and female reproductive molecules. In this study, we use short-read RNA sequencing to examine the evolutionary dynamics of female reproductive genes and the postmating transcriptome response in crosses within and between species. First, we found that most female reproductive tract genes are slow-evolving compared to the genome average. Second, postmating responses in con- and heterospecific matings are largely congruent, but heterospecific matings induce expression of additional stress-response genes. Some of those are immunity genes that are activated by the Imd pathway. We also identify several genes in the JAK/STAT signaling pathway that are induced in heterospecific, but not conspecific mating. While this immune response was most pronounced in the female reproductive tract, we also detect it in the female head and ovaries. These results show that the female's postmating transcriptome-level response is determined in part by the genotype of the male, and that divergence in male reproductive genes and/or traits can have immunogenic effects on females.

Genetic architecture of a body colour cline in Drosophila americana

Phenotypic variation within a species is often structured geographically in clines. In Drosophila americana, a longitudinal cline for body colour exists within North America that appears to be due to local adaptation. The tan and ebony genes have been hypothesized to contribute to this cline, with alleles of both genes that lighten body colour found in D. americana. These alleles are similar in sequence and function to the allele fixed in D. americana's more lightly pigmented sister species, Drosophila novamexicana. Here, we examine the frequency and geographic distribution of these D. novamexicana-like alleles in D. americana. Among alleles from over 100 strains of D. americana isolated from 21 geographic locations, we failed to identify additional alleles of tan or ebony with as much sequence similarity to D. novamexicana as the D. novamexicana-like alleles previously described. However, using genetic analysis of 51 D. americana strains derived from 20 geographic locations, we identified one new allele of ebony and one new allele of tan segregating in D. americana that are functionally equivalent to the D. novamexicana allele. An additional 5 alleles of tan also showed marginal evidence of functional similarity. Given the rarity of these alleles, however, we conclude that they are unlikely to be driving the pigmentation cline. Indeed, phenotypic distributions of the 51 backcross populations analysed indicate a more complex genetic architecture, with diversity in the number and effects of loci altering pigmentation observed both within and among populations of D. americana. This genetic heterogeneity poses a challenge to association studies and genomic scans for clinal variation, but might be common in natural populations.

ebony Affects Pigmentation Divergence and Cuticular Hydrocarbons in Drosophila americana and D. novamexicana

Drosophila pigmentation has been a fruitful model system for understanding the genetic and developmental mechanisms underlying phenotypic evolution. For example, prior work has shown that divergence of the tan gene contributes to pigmentation differences between two members of the virilis group: Drosophila novamexicana, which has a light yellow body color, and D. americana, which has a dark brown body color. Quantitative trait locus (QTL) mapping and expression analysis has suggested that divergence of the ebony gene might also contribute to pigmentation differences between these two species. Here, we directly test this hypothesis by using CRISPR/Cas9 genome editing to generate ebony null mutants in D. americana and D. novamexicana and then using reciprocal hemizygosity testing to compare the effects of each species' ebony allele on pigmentation. We find that divergence of ebony does indeed contribute to the pigmentation divergence between species, with effects on both the overall body color as well as a difference in pigmentation along the dorsal abdominal midline. Motivated by recent work in D. melanogaster, we also used the ebony null mutants to test for effects of ebony on cuticular hydrocarbon (CHC) profiles. We found that ebony affects CHC abundance in both species, but does not contribute to qualitative differences in the CHC profiles between these two species. Additional transgenic resources for working with D. americana and D. novamexicana, such as white mutants of both species and yellow mutants in D. novamexicana, were generated in the course of this work and are also described. Taken together, this study advances our understanding of loci contributing to phenotypic divergence and illustrates how the latest genome editing tools can be used for functional testing in non-model species.

Post-ejaculatory modifications to sperm (PEMS)

Mammalian sperm must spend a minimum period of time within a female reproductive tract to achieve the capacity to fertilize oocytes. This phenomenon, termed sperm 'capacitation', was discovered nearly seven decades ago and opened a window into the complexities of sperm-female interaction. Capacitation is most commonly used to refer to a specific combination of processes that are believed to be widespread in mammals and includes modifications to the sperm plasma membrane, elevation of intracellular cyclic AMP levels, induction of protein tyrosine phosphorylation, increased intracellular Ca2+ levels, hyperactivation of motility, and, eventually, the acrosome reaction. Capacitation is only one example of post-ejaculatory modifications to sperm (PEMS) that are widespread throughout the animal kingdom. Although PEMS are less well studied in non-mammalian taxa, they likely represent the rule rather than the exception in species with internal fertilization. These PEMS are diverse in form and collectively represent the outcome of selection fashioning complex maturational trajectories of sperm that include multiple, sequential phenotypes that are specialized for stage-specific functionality within the female. In many cases, PEMS are critical for sperm to migrate successfully through the female reproductive tract, survive a protracted period of storage, reach the site of fertilization and/or achieve the capacity to fertilize eggs. We predict that PEMS will exhibit widespread phenotypic plasticity mediated by sperm-female interactions. The successful execution of PEMS thus has important implications for variation in fitness and the operation of post-copulatory sexual selection. Furthermore, it may provide a widespread mechanism of reproductive isolation and the maintenance of species boundaries. Despite their possible ubiquity and importance, the investigation of PEMS has been largely descriptive, lacking any phylogenetic consideration with regard to divergence, and there have been no theoretical or empirical investigations of their evolutionary significance. Here, we (i) clarify PEMS-related nomenclature; (ii) address the evolutionary origin, maintenance and divergence in PEMS in the context of the protracted life history of sperm and the complex, selective environment of the female reproductive tract; (iii) describe taxonomically widespread types of PEMS: sperm activation, chemotaxis and the dissociation of sperm conjugates; (iv) review the occurence of PEMS throughout the animal kingdom; (v) consider alternative hypotheses for the adaptive value of PEMS; (vi) speculate on the evolutionary implications of PEMS for genomic architecture, sexual selection, and reproductive isolation; and (vii) suggest fruitful directions for future functional and evolutionary analyses of PEMS.

The effects of heterospecific mating frequency on the strength of cryptic reproductive barriers

Heterospecific mating frequency is critical to hybrid zone dynamics and can directly impact the strength of reproductive barriers and patterns of introgression. The effectiveness of post-mating prezygotic (PMPZ) reproductive barriers, which include reduced fecundity via heterospecific matings and conspecific sperm precedence, may depend on the number, identity and order of mates. Studies of PMPZ barriers suggest that they may be important in many systems, but whether these barriers are effective at realistic heterospecific mating frequencies has not been tested. Here, we evaluate the strength of cryptic reproductive isolation in two leaf beetles (Chrysochus auratus and C. cobaltinus) in the context of a range of heterospecific mating frequencies observed in natural populations. We found both species benefited from multiple matings, but the benefits were greater in C. cobaltinus and extended to heterospecific matings. We found that PMPZ barriers greatly limited hybrid production by C. auratus females with moderate heterospecific mating frequencies, but that their effectiveness diminished at higher heterospecific mating frequencies. In contrast, there was no evidence for PMPZ barriers in C. cobaltinus females at any heterospecific mating frequency. We show that integrating realistic estimates of cryptic isolation with information on relative abundance and heterospecific mating frequency in the field substantially improves our understanding of the strong directional bias in F1 production previously documented in the Chrysochus hybrid zone. Our results demonstrate that heterospecific mating frequency is critical to understanding the impact of cryptic post-copulatory barriers on hybrid zone structure and dynamics, and that future studies of such barriers should incorporate field-relevant heterospecific mating frequencies.

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.

Evolutionary Dynamics of Male Reproductive Genes in the Drosophila virilis Subgroup

Postcopulatory sexual selection (PCSS) is a potent evolutionary force that can drive rapid changes of reproductive genes within species, and thus has the potential to generate reproductive incompatibilities between species. Male seminal fluid proteins (SFPs) are major players in postmating interactions, and are important targets of PCSS in males. The virilis subgroup of Drosophila exhibits strong interspecific gametic incompatibilities, and can serve as a model to study the genetic basis of PCSS and gametic isolation. However, reproductive genes in this group have not been characterized. Here we utilize short-read RNA sequencing of male reproductive organs to examine the evolutionary dynamics of reproductive genes in members of the virilis subgroup: D. americana, D. lummei, D. novamexicana, and D. virilis. We find that the majority of male reproductive transcripts are testes-biased, accounting for ∼15% of all annotated genes. Ejaculatory bulb (EB)-biased transcripts largely code for lipid metabolic enzymes, and contain orthologs of the D. melanogaster EB protein, Peb-me, which is involved in mating-plug formation. In addition, we identify 71 candidate SFPs, and show that this gene set has the highest rate of nonsynonymous codon substitution relative to testes- and EB-biased genes. Furthermore, we identify orthologs of 35 D. melanogaster SFPs that have conserved accessory gland expression in the virilis group. Finally, we show that several of the SFPs that have the highest rate of nonsynonymous codon substitution reside on chromosomal regions, which contributes to paternal gametic incompatibility between species. Our results show that SFPs rapidly diversify in the virilis group, and suggest that they likely play a role in PCSS and/or gametic isolation.

Multiple Genes Cause Postmating Prezygotic Reproductive Isolation in the Drosophila virilis Group

Understanding the genetic basis of speciation is a central problem in evolutionary biology. Studies of reproductive isolation have provided several insights into the genetic causes of speciation, especially in taxa that lend themselves to detailed genetic scrutiny. Reproductive barriers have usually been divided into those that occur before zygote formation (prezygotic) and after (postzygotic), with the latter receiving a great deal of attention over several decades. Reproductive barriers that occur after mating but before zygote formation [postmating prezygotic (PMPZ)] are especially understudied at the genetic level. Here, I present a phenotypic and genetic analysis of a PMPZ reproductive barrier between two species of the Drosophila virilis group: D. americana and D. virilis This species pair shows strong PMPZ isolation, especially when D. americana males mate with D. virilis females: ∼99% of eggs laid after these heterospecific copulations are not fertilized. Previous work has shown that the paternal loci contributing to this incompatibility reside on two chromosomes, one of which (chromosome 5) likely carries multiple factors. The other (chromosome 2) is fixed for a paracentric inversion that encompasses nearly half the chromosome. Here, I present two results. First, I show that PMPZ in this species cross is largely due to defective sperm storage in heterospecific copulations. Second, using advanced intercross and backcross mapping approaches, I identify genomic regions that carry genes capable of rescuing heterospecific fertilization. I conclude that paternal incompatibility between D. americana males and D. virilis females is underlain by four or more genes on chromosomes 2 and 5.

Postmating-prezygotic isolation between two allopatric populations of Drosophila montana: fertilisation success differs under sperm competition

Postmating but prezygotic (PMPZ) interactions are increasingly recognized as a potentially important early‐stage barrier in the evolution of reproductive isolation. A recent study described a potential example between populations of the same species: single matings between Drosophila montana populations resulted in differential fertilisation success because of the inability of sperm from one population (Vancouver) to penetrate the eggs of the other population (Colorado). As the natural mating system of D. montana is polyandrous (females remate rapidly), we set up double matings of all possible crosses between the same populations to test whether competitive effects between ejaculates influence this PMPZ isolation. We measured premating isolation in no‐choice tests, female fecundity, fertility and egg‐to‐adult viability after single and double matings as well as second‐male paternity success (P₂). Surprisingly, we found no PMPZ reproductive isolation between the two populations under a competitive setting, indicating no difficulty of sperm from Vancouver males to fertilize Colorado eggs after double matings. While there were subtle differences in how P₂ changed over time, suggesting that Vancouver males’ sperm are somewhat less competitive in a first‐male role within Colorado females, these effects did not translate into differences in overall P₂. Fertilisation success can thus differ dramatically between competitive and noncompetitive conditions, perhaps because the males that mate second produce higher quality ejaculates in response to sperm competition. We suggest that unlike in more divergent species comparisons, where sperm competition typically increases reproductive isolation, ejaculate tailoring can reduce the potential for PMPZ isolation when recently diverged populations interbreed.

Reproductive isolation among allopatric Drosophila montana populations

An outstanding goal in speciation research is to trace the mode and tempo of the evolution of barriers to gene flow. Such research benefits from studying incipient speciation, in which speciation between populations has not yet occurred, but where multiple potential mechanisms of reproductive isolation (RI: i.e., premating, postmating-prezygotic (PMPZ), and postzygotic barriers) may act. We used such a system to investigate these barriers among allopatric populations of Drosophila montana. In all heteropopulation crosses we found premating (sexual) isolation, which was either symmetric or asymmetric depending on the population pair compared. Postmating isolation was particularly strong in crosses involving males from one of the study populations, and while sperm were successfully transferred, stored, and motile, we experimentally demonstrated that the majority of eggs produced were unfertilized. Thus, we identified the nature of a PMPZ incompatibility. There was no evidence of intrinsic postzygotic effects. Measures of absolute and relative strengths of pre- and postmating barriers showed that populations differed in the mode and magnitude of RI barriers. Our results indicate that incipient RI among populations can be driven by different contributions of both premating and PMPZ barriers occurring between different population pairs and without the evolution of postzygotic barriers.

A single gene causes an interspecific difference in pigmentation in Drosophila

The genetic basis of species differences remains understudied. Studies in insects have contributed significantly to our understanding of morphological evolution. Pigmentation traits in particular have received a great deal of attention and several genes in the insect pigmentation pathway have been implicated in inter- and intraspecific differences. Nonetheless, much remains unknown about many of the genes in this pathway and their potential role in understudied taxa. Here we genetically analyze the puparium color difference between members of the virilis group of Drosophila. The puparium of Drosophila virilis is black, while those of D. americana, D. novamexicana, and D. lummei are brown. We used a series of backcross hybrid populations between D. americana and D. virilis to map the genomic interval responsible for the difference between this species pair. First, we show that the pupal case color difference is caused by a single Mendelizing factor, which we ultimately map to an ∼11-kb region on chromosome 5. The mapped interval includes only the first exon and regulatory region(s) of the dopamine N-acetyltransferase gene (Dat). This gene encodes an enzyme that is known to play a part in the insect pigmentation pathway. Second, we show that this gene is highly expressed at the onset of pupation in light brown taxa (D. americana and D. novamexicana) relative to D. virilis, but not in the dark brown D. lummei. Finally, we examine the role of Dat in adult pigmentation between D. americana (heavily melanized) and D. novamexicana (lightly melanized) and find no discernible effect of this gene in adults. Our results demonstrate that a single gene is entirely or almost entirely responsible for a morphological difference between species.