Relaxed Selection and the Rapid Evolution of Reproductive Genes

Genome-wide analysis reveals the contributors to fast molecular evolution of the Chinese hook snout carp (Opsariichthys bidens)

Variations in molecular evolutionary rate have been widely investigated among lineages and genes. However, it remains an open question whether fast rate of molecular evolution is driven by natural selection or random drift, and how the fast rate is linked to metabolic rate. Additionally, previous studies on fast molecular evolution have been largely restricted to concatenated matrix of genes or a few specifically selected genes, but less is known for individual genes at the genome-wide level. Here we addressed these questions using more than 5000 single-copy orthologous (SCO) genes through comparative genomic and phylogenetic analyses among fishes, with a special focus on a newly-sequenced clupeocephalan fish the Chinese hook snout carp Opsariichthys bidens. We showed O. bidens displays significantly higher mean substitution rate and more fast-evolving SCO genes (2172 genes) than most fishes studied here. The rapidly evolving genes are enriched in highly conserved and very basic functions such as translation and ribosome that are critical for biological fitness. We further revealed that ∼25 % of these fast-evolving genes exhibit a constant increase of substitution rate from the common ancestor down to the present, suggesting a neglected but important contribution from ancestral states. Model fitting showed that ∼85 % of fast-evolving genes exclusive to O. bidens and related species follow the adaptive evolutionary model rather than random-drift model, and 7.6 % of fast-evolving genes identified in O. bidens have experienced positive selection, both indicating the reflection of adaptive selection. Finally, metabolic rate was observed to be linked with substitution rate in a gene-specific manner. Overall, our findings reveal fast molecular evolution of SCO genes at genome-wide level in O. bidens, and uncover the evolutionary and ecological contributors to it.

Evolution and genetics of accessory gland transcriptome divergence between Drosophila melanogaster and D. simulans

Studies of allele-specific expression in interspecific hybrids have provided important insights into gene-regulatory divergence and hybrid incompatibilities. Many such investigations in Drosophila have used transcriptome data from complex mixtures of many tissues or from gonads, however, regulatory divergence may vary widely among species, sexes, and tissues. Thus, we lack sufficiently broad sampling to be confident about the general biological principles of regulatory divergence. Here, we seek to fill some of these gaps in the literature by characterizing regulatory evolution and hybrid misexpression in a somatic male sex organ, the accessory gland, in F1 hybrids between Drosophila melanogaster and D. simulans. The accessory gland produces seminal fluid proteins, which play an important role in male and female fertility and may be subject to adaptive divergence due to male-male or male-female interactions. We find that trans differences are relatively more abundant than cis, in contrast to most of the interspecific hybrid literature, though large effect-size trans differences are rare. Seminal fluid protein genes have significantly elevated levels of expression divergence and tend to be regulated through both cis and trans divergence. We find limited misexpression (over- or underexpression relative to both parents) in this organ compared to most other Drosophila studies. As in previous studies, male-biased genes are overrepresented among misexpressed genes and are much more likely to be underexpressed. ATAC-Seq data show that chromatin accessibility is correlated with expression differences among species and hybrid allele-specific expression. This work identifies unique regulatory evolution and hybrid misexpression properties of the accessory gland and suggests the importance of tissue-specific allele-specific expression studies.

Positive selection and relaxed purifying selection contribute to rapid evolution of male-biased genes in a dioecious flowering plant

Sex-biased genes offer insights into the evolution of sexual dimorphism. Sex-biased genes, especially those with male bias, show elevated evolutionary rates of protein sequences driven by positive selection and relaxed purifying selection in animals. Although rapid sequence evolution of sex-biased genes and evolutionary forces have been investigated in animals and brown algae, less is known about evolutionary forces in dioecious angiosperms. In this study, we separately compared the expression of sex-biased genes between female and male floral buds and between female and male flowers at anthesis in dioecious Trichosanthes pilosa (Cucurbitaceae). In floral buds, sex-biased gene expression was pervasive, and had significantly different roles in sexual dimorphism such as physiology. We observed higher rates of sequence evolution for male-biased genes in floral buds compared to female-biased and unbiased genes. Male-biased genes under positive selection were mainly associated with functions to abiotic stress and immune responses, suggesting that high evolutionary rates are driven by adaptive evolution. Additionally, relaxed purifying selection may contribute to accelerated evolution in male-biased genes generated by gene duplication. Our findings, for the first time in angiosperms, suggest evident rapid evolution of male-biased genes, advance our understanding of the patterns and forces driving the evolution of sexual dimorphism in dioecious plants.

Gametic selection favours polyandry and selfing

Competition among pollen or sperm (gametic selection) can cause evolution. Mating systems shape the intensity of gametic selection by determining the competitors involved, which can in turn cause the mating system itself to evolve. We model the bidirectional relationship between gametic selection and mating systems, focusing on variation in female mating frequency (monandry-polyandry) and self-fertilisation (selfing-outcrossing). First, we find that monandry and selfing both reduce the efficiency of gametic selection in removing deleterious alleles. This means that selfing can increase mutation load, in contrast to cases without gametic selection where selfing purges deleterious mutations and decreases mutation load. Second, we explore how mating systems evolve via their effect on gametic selection. By manipulating gametic selection, polyandry can evolve to increase the fitness of the offspring produced. However, this indirect advantage of post-copulatory sexual selection is weak and is likely to be overwhelmed by any direct fitness effects of mating systems. Nevertheless, gametic selection can be potentially decisive for selfing evolution because it significantly reduces inbreeding depression, which favours selfing. Thus, the presence of gametic selection could be a key factor driving selfing evolution.

Sex-biased gene content is associated with sex chromosome turnover in Danaini butterflies

Sex chromosomes play an outsized role in adaptation and speciation, and thus deserve particular attention in evolutionary genomics. In particular, fusions between sex chromosomes and autosomes can produce neo-sex chromosomes, which offer important insights into the evolutionary dynamics of sex chromosomes. Here, we investigate the evolutionary origin of the previously reported Danaus neo-sex chromosome within the tribe Danaini. We assembled and annotated genomes of Tirumala septentrionis (subtribe Danaina), Ideopsis similis (Amaurina), Idea leuconoe (Euploeina) and Lycorea halia (Itunina) and identified their Z-linked scaffolds. We found that the Danaus neo-sex chromosome resulting from the fusion between a Z chromosome and an autosome corresponding to the Melitaea cinxia chromosome (McChr) 21 arose in a common ancestor of Danaina, Amaurina and Euploina. We also identified two additional fusions as the W chromosome further fused with the synteny block McChr31 in I. similis and independent fusion occurred between ancestral Z chromosome and McChr12 in L. halia. We further tested a possible role of sexually antagonistic selection in sex chromosome turnover by analysing the genomic distribution of sex-biased genes in I. leuconoe and L. halia. The autosomes corresponding to McChr21 and McChr31 involved in the fusions are significantly enriched in female- and male-biased genes, respectively, which could have hypothetically facilitated fixation of the neo-sex chromosomes. This suggests a role of sexual antagonism in sex chromosome turnover in Lepidoptera. The neo-Z chromosomes of both I. leuconoe and L. halia appear fully compensated in somatic tissues, but the extent of dosage compensation for the ancestral Z varies across tissues and species.

Sexual conflict, heterochrony and tissue specificity as evolutionary problems of adaptive plasticity in development

Differential gene expression represents a fundamental cause and manifestation of phenotypic plasticity. Adaptive phenotypic plasticity in gene expression as a trait evolves when alleles that mediate gene regulation serve to increase organismal fitness by improving the alignment of variation in gene expression with variation in circumstances. Among the diverse circumstances that a gene encounters are distinct cell types, developmental stages and sexes, as well as an organism's extrinsic ecological environments. Consequently, adaptive phenotypic plasticity provides a common framework to consider diverse evolutionary problems by considering the shared implications of alleles that produce context-dependent gene expression. From this perspective, adaptive plasticity represents an evolutionary resolution to conflicts of interest that arise from any negatively pleiotropic effects of expression of a gene across ontogeny, among tissues, between the sexes, or across extrinsic environments. This view highlights shared properties within the general relation of fitness, trait expression and context that may nonetheless differ substantively in the grain of selection within and among generations to influence the likelihood of adaptive plasticity as an evolutionary response. Research programmes that historically have focused on these separate issues may use the insights from one another by recognizing their shared dependence on context-dependent gene regulatory evolution.

Disentangling Positive Selection from Relaxed Selection in Animal Mitochondrial Genomes

AbstractDisentangling different types of selection is a common goal in molecular evolution. Elevated dN/dS ratios (the ratio of nonsynonymous to synonymous substitution rates) in focal lineages are often interpreted as signs of positive selection. Paradoxically, relaxed purifying selection can also result in elevated dN/dS ratios, but tests to distinguish these two causes are seldomly implemented. Here, we reevaluated seven case studies describing elevated dN/dS ratios in animal mitochondrial DNA (mtDNA) and their accompanying hypotheses regarding selection. They included flightless lineages versus flighted lineages in birds, bats, and insects and physiological adaptations in snakes, two groups of electric fishes, and primates. We found that elevated dN/dS ratios were often not caused by the predicted mechanism, and we sometimes found strong support for the opposite mechanism. We discuss reasons why energetic hypotheses may be confounded by other selective forces acting on mtDNA and caution against overinterpreting singular molecular signals, including elevated dN/dS ratios.

A predominant role of genotypic variation in both expression of sperm competition genes and paternity success in Drosophila melanogaster

Sperm competition is a crucial aspect of male reproductive success in many species, including Drosophila melanogaster, and seminal fluid proteins (Sfps) can influence sperm competitiveness. However, the combined effect of environmental and genotypic variation on sperm competition gene expression remains poorly understood. Here, we used Drosophila Genetic Reference Panel (DGRP) inbred lines and manipulated developmental population density (i.e. larval density) to test the effects of genotype, environment and genotype-by-environment interactions (GEI) on the expression of the known sperm competition genes Sex Peptide, Acp36DE and CG9997. High larval density resulted in reduced adult body size, but expression of sperm competition genes remained unaffected. Furthermore, we found no significant GEI but genotypic effects in the expression of SP and Acp36DE. Our results also revealed GEI for relative competitive paternity success (second male paternity; P2), with genes' expression positively correlated with P2. Given the effect of genotype on the expression of genes, we conducted a genome-wide association study (GWAS) and identified polymorphisms in putative cis-regulatory elements as predominant factors regulating the expression of SP and Acp36DE. The association of genotypic variation with sperm competition outcomes, and the resilience of sperm competition genes' expression against environmental challenges, demonstrates the importance of genome variation background in reproductive fitness.

Phytochemical Cue for the Fitness Costs of Herbicide-Resistant Weeds

Despite increasing knowledge of the fitness costs of viability and fecundity involved in the herbicide-resistant weeds, relatively little is known about the linkage between herbicide resistance costs and phytochemical cues in weed species and biotypes. This study demonstrated relative fitness and phytochemical responses in six herbicide-resistant weeds and their susceptible counterparts. There were significant differences in the parameters of viability (growth and photosynthesis), fecundity fitness (flowering and seed biomass) and a ubiquitous phytochemical (-)-loliolide levels between herbicide-resistant weeds and their susceptible counterparts. Fitness costs occurred in herbicide-resistant Digitaria sanguinalis and Leptochloa chinensis but they were not observed in herbicide-resistant Alopecurus japonicas, Eleusine indica, Ammannia arenaria, and Echinochloa crus-galli. Correlation analysis indicated that the morphological characteristics of resistant and susceptible weeds were negatively correlated with (-)-loliolide concentration, but positively correlated with lipid peroxidation malondialdehyde and total phenol contents. Principal component analysis showed that the lower the (-)-loliolide concentration, the stronger the adaptability in E. crus-galli and E. indica. Therefore, not all herbicide-resistant weeds have fitness costs, but the findings showed several examples of resistance leading to improved fitness even in the absence of herbicides. In particular, (-)-loliolide may act as a phytochemical cue to explain the fitness cost of herbicide-resistant weeds by regulating vitality and fecundity.

Relative rates of evolution of male-beneficial nuclear compensatory mutations and male-harming Mother's Curse mitochondrial alleles

Mother's Curse alleles represent a significant source of potential male fitness defects. The maternal inheritance of mutations with the pattern of sex-specific fitness effects, s♀>0>s♂, allows Mother's Curse alleles to spread through a population even though they reduce male fitness. Although the mitochondrial genomes of animals contain only a handful of protein-coding genes, mutations in many of these genes have been shown to have a direct effect on male fertility. The evolutionary process of nuclear compensation is hypothesized to counteract the male-limited mitochondrial defects that spread via Mother's Curse. Here we use population genetic models to investigate the evolution of compensatory autosomal nuclear mutations that act to restore the loss of fitness caused by mitochondrial mutation pressures. We derive the rate of male fitness deterioration by Mother's Curse and the rate of restoration by nuclear compensatory evolution. We find that the rate of nuclear gene compensation is many times slower than that of its deterioration by cytoplasmic mutation pressure, resulting in a significant lag in the recovery of male fitness. Thus, the numbers of nuclear genes capable of restoring male mitochondrial fitness defects must be large in order to sustain male fitness in the face of mutation pressures.

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.

Speciation and development

Understanding general principles about the origin of species remains one of the foundational challenges in evolutionary biology. The genomic divergence between groups of individuals can spawn hybrid inviability and hybrid sterility, which presents a tantalizing developmental problem. Divergent developmental programs may yield either conserved or divergent phenotypes relative to ancestral traits, both of which can be responsible for reproductive isolation during the speciation process. The genetic mechanisms of developmental evolution involve cis- and trans-acting gene regulatory change, protein-protein interactions, genetic network structures, dosage, and epigenetic regulation, all of which also have roots in population genetic and molecular evolutionary processes. Toward the goal of demystifying Darwin's "mystery of mysteries," this review integrates microevolutionary concepts of genetic change with principles of organismal development, establishing explicit links between population genetic process and developmental mechanisms in the production of macroevolutionary pattern. This integration aims to establish a more unified view of speciation that binds process and mechanism.

The roles of sexual selection and sexual conflict in shaping patterns of genome and transcriptome variation

Sexual dimorphism is one of the most prevalent, and often the most extreme, examples of phenotypic variation within species, and arises primarily from genomic variation that is shared between females and males. Many sexual dimorphisms arise through sex differences in gene expression, and sex-biased expression is one way that a single, shared genome can generate multiple, distinct phenotypes. Although many sexual dimorphisms are expected to result from sexual selection, and many studies have invoked the possible role of sexual selection to explain sex-specific traits, the role of sexual selection in the evolution of sexually dimorphic gene expression remains difficult to differentiate from other forms of sex-specific selection. In this Review, we propose a holistic framework for the study of sex-specific selection and transcriptome evolution. We advocate for a comparative approach, across tissues, developmental stages and species, which incorporates an understanding of the molecular mechanisms, including genomic variation and structure, governing gene expression. Such an approach is expected to yield substantial insights into the evolution of genetic variation and have important applications in a variety of fields, including ecology, evolution and behaviour.

Sex-biased gene expression at single-cell resolution: cause and consequence of sexual dimorphism

Gene expression differences between males and females are thought to be key for the evolution of sexual dimorphism, and sex-biased genes are often used to study the molecular footprint of sex-specific selection. However, gene expression is often measured from complex aggregations of diverse cell types, making it difficult to distinguish between sex differences in expression that are due to regulatory rewiring within similar cell types and those that are simply a consequence of developmental differences in cell-type abundance. To determine the role of regulatory versus developmental differences underlying sex-biased gene expression, we use single-cell transcriptomic data from multiple somatic and reproductive tissues of male and female guppies, a species that exhibits extensive phenotypic sexual dimorphism. Our analysis of gene expression at single-cell resolution demonstrates that nonisometric scaling between the cell populations within each tissue and heterogeneity in cell-type abundance between the sexes can influence inferred patterns of sex-biased gene expression by increasing both the false-positive and false-negative rates. Moreover, we show that, at the bulk level, the subset of sex-biased genes that are the product of sex differences in cell-type abundance can significantly confound patterns of coding-sequence evolution. Taken together, our results offer a unique insight into the effects of allometry and cellular heterogeneity on perceived patterns of sex-biased gene expression and highlight the power of single-cell RNA-sequencing in distinguishing between sex-biased genes that are the result of regulatory change and those that stem from sex differences in cell-type abundance, and hence are a consequence rather than a cause of sexual dimorphism.

Domestication affects sex-biased gene expression evolution in the duck

Genes with sex-biased expression are thought to underlie sexually dimorphic phenotypes and are therefore subject to different selection pressures in males and females. Many authors have proposed that sexual conflict leads to the evolution of sex-biased expression, which allows males and females to reach separate phenotypic and fitness optima. The selection pressures associated with domestication may cause changes in population architectures and mating systems, which in turn can alter their direction and strength. We compared sex-biased expression and genetic signatures in wild and domestic ducks (Anas platyrhynchos), and observed changes of sexual selection and identified the genomic divergence affected by selection forces. The extent of sex-biased expression in both sexes is positively correlated with the level of both d _N /d _S and nucleotide diversity. This observed changing pattern may mainly be owing to relaxed genetic constraints. We also demonstrate a clear link between domestication and sex-biased evolutionary rate in a comparative framework. Decreased polymorphism and evolutionary rate in domesticated populations generally matched life-history phenotypes known to experience artificial selection. Taken together, our work suggests the important implications of domestication in sex-biased evolution and the roles of artificial selection and sexual selection for shaping the diversity and evolutionary rate of the genome.

Monsoon boosted radiation of the endemic East Asian carps

Major historical events often trigger the rapid flourishing of a few lineages, which in turn shape established biodiversity patterns. How did this process occur and develop? This study provides a window into this issue. The endemic East Asian carps (EEAC) dominated the ichthyofauna of East Asia and exhibited a high degree of adaptation to monsoonal river-lake ecosystems. A series of evidence, including ecogeography, phylogenetics, and macroevolution, suggests that the EEAC is a lineage that arose with the East Asian monsoon and thrived intimately with subsequent monsoon activities. We further deduce the evolution of the EEAC and find that a range of historical events in the monsoon setting (e.g., marine transgression and regression and glacial-interglacial cycle) have further reshaped the distribution patterns of EEAC's members. Comparative genomics analyses reveal that introgressions during the initial period of EEAC radiation and innovations in the regulation of the brain and nervous system may have aided their adaptation to river-lake ecosystems in a monsoon setting, which boosted radiation. Overall, this study strengthens knowledge of the evolutionary patterns of freshwater fishes in East Asia and provides a model case for understanding the impact of major historical events on the evolution of biota.

Discovering novel reproductive genes in a non-model fly using de novo GridION transcriptomics

Gene discovery has important implications for investigating phenotypic trait evolution, adaptation, and speciation. Male reproductive tissues, such as accessory glands (AGs), are hotspots for recruitment of novel genes that diverge rapidly even among closely related species/populations. These genes synthesize seminal fluid proteins that often affect post-copulatory sexual selection-they can mediate male-male sperm competition, ejaculate-female interactions that modify female remating and even influence reproductive incompatibilities among diverging species/populations. Although de novo transcriptomics has facilitated gene discovery in non-model organisms, reproductive gene discovery is still challenging without a reference database as they are often novel and bear no homology to known proteins. Here, we use reference-free GridION long-read transcriptomics, from Oxford Nanopore Technologies (ONT), to discover novel AG genes and characterize their expression in the widespread dung fly, Sepsis punctum. Despite stark population differences in male reproductive traits (e.g.: Body size, testes size, and sperm length) as well as female re-mating, the male AG genes and their secretions of S. punctum are still unknown. We implement a de novo ONT transcriptome pipeline incorporating quality-filtering and rigorous error-correction procedures, and we evaluate gene sequence and gene expression results against high-quality Illumina short-read data. We discover highly-expressed reproductive genes in AG transcriptomes of S. punctum consisting of 40 high-quality and high-confidence ONT genes that cross-verify against Illumina genes, among which 26 are novel and specific to S. punctum. Novel genes account for an average of 81% of total gene expression and may be functionally relevant in seminal fluid protein production. For instance, 80% of genes encoding secretory proteins account for 74% total gene expression. In addition, median sequence similarities of ONT nucleotide and protein sequences match within-Illumina sequence similarities. Read-count based expression quantification in ONT is congruent with Illumina's Transcript per Million (TPM), both in overall pattern and within functional categories. Rapid genomic innovation followed by recruitment of de novo genes for high expression in S. punctum AG tissue, a pattern observed in other insects, could be a likely mechanism of evolution of these genes. The study also demonstrates the feasibility of adapting ONT transcriptomics for gene discovery in non-model systems.

Mating strategy predicts gene presence/absence patterns in a genus of simultaneously hermaphroditic flatworms

Gene repertoire turnover is a characteristic of genome evolution. However, we lack well-replicated analyses of presence/absence patterns associated with different selection contexts. Here, we study ∼100 transcriptome assemblies across Macrostomum, a genus of simultaneously hermaphroditic flatworms exhibiting multiple convergent shifts in mating strategy and associated reproductive morphologies. Many species mate reciprocally, with partners donating and receiving sperm at the same time. Other species convergently evolved to mate by hypodermic injection of sperm into the partner. We find that for orthologous transcripts annotated as expressed in the body region containing the testes, sequences from hypodermically inseminating species diverge more rapidly from the model species, Macrostomum lignano, and have a lower probability of being observed in other species. For other annotation categories, simpler models with a constant rate of similarity decay with increasing genetic distance from M. lignano match the observed patterns well. Thus, faster rates of sequence evolution for hypodermically inseminating species in testis-region genes result in higher rates of homology detection failure, yielding a signal of rapid evolution in sequence presence/absence patterns. Our results highlight the utility of considering appropriate null models for unobserved genes, as well as associating patterns of gene presence/absence with replicated evolutionary events in a phylogenetic context.

Novel genome sequence of Chinese cavefish (Triplophysa rosa) reveals pervasive relaxation of natural selection in cavefish genomes

All cavefishes, living exclusively in caves across the globe, exhibit similar phenotypic traits, including the characteristic loss of eyes. To understand whether such phenotypic convergence shares similar genomic bases, here we investigated genome-wide evolutionary signatures of cavefish phenotypes by comparing whole-genome sequences of three pairs of cavefishes and their surface fish relatives. Notably, we newly sequenced and generated a whole-genome assembly of the Chinese cavefish Triplophysa rosa. Our comparative analyses revealed several shared features of cavefish genome evolution. Cavefishes had lower mutation rates than their surface fish relatives. In contrast, the ratio of nonsynonymous to synonymous substitutions (ω) was significantly elevated in cavefishes compared to in surface fishes, consistent with the relaxation of purifying selection. In addition, cavefish genomes had an increased mutational load, including mutations that alter protein hydrophobicity profiles, which were considered harmful. Interestingly, however, we found no overlap in positively selected genes among different cavefish lineages, indicating that the phenotypic convergence in cavefishes was not caused by positive selection of the same sets of genes. Analyses of previously identified candidate genes associated with cave phenotypes supported this conclusion. Genes belonging to the lipid metabolism functional ontology were under relaxed purifying selection in all cavefish genomes, which may be associated with the nutrient-poor habitat of cavefishes. Our work reveals previously uncharacterized patterns of cavefish genome evolution and provides comparative insights into the evolution of cave-associated phenotypic traits.

A possible genomic footprint of polygenic adaptation on population divergence in seed beetles?

Efforts to unravel the genomic basis of incipient speciation are hampered by a mismatch between our toolkit and our understanding of the ecology and genetics of adaptation. While the former is focused on detecting selective sweeps involving few independently acting or linked speciation genes, the latter states that divergence typically occurs in polygenic traits under stabilizing selection. Here, we ask whether a role of stabilizing selection on polygenic traits in population divergence may be unveiled by using a phenotypically informed integrative approach, based on genome‐wide variation segregating in divergent populations. We compare three divergent populations of seed beetles (Callosobruchus maculatus) where previous work has demonstrated a prominent role for stabilizing selection on, and population divergence in, key life history traits that reflect rate‐dependent metabolic processes. We derive and assess predictions regarding the expected pattern of covariation between genetic variation segregating within populations and genetic differentiation between populations. Population differentiation was considerable (mean F_ST = 0.23–0.26) and was primarily built by genes showing high selective constraints and an imbalance in inferred selection in different populations (positive Tajima's D_NS in one and negative in one), and this set of genes was enriched with genes with a metabolic function. Repeatability of relative population differentiation was low at the level of individual genes but higher at the level of broad functional classes, again spotlighting metabolic genes. Absolute differentiation (d_XY) showed a very different general pattern at this scale of divergence, more consistent with an important role for genetic drift. Although our exploration is consistent with stabilizing selection on polygenic metabolic phenotypes as an important engine of genome‐wide relative population divergence and incipient speciation in our study system, we note that it is exceedingly difficult to firmly exclude other scenarios.

Functional Diversity and Evolution of the Drosophila Sperm Proteome

Spermatozoa are central to fertilization and the evolutionary fitness of sexually reproducing organisms. As such, a deeper understanding of sperm proteomes (and associated reproductive tissues) has proven critical to the advancement of the fields of sexual selection and reproductive biology. Due to their extraordinary complexity, proteome depth-of-coverage is dependent on advancements in technology and related bioinformatics, both of which have made significant advancements in the decade since the last Drosophila sperm proteome was published. Here, we provide an updated version of the Drosophila melanogaster sperm proteome (DmSP3) using improved separation and detection methods and an updated genome annotation. Combined with previous versions of the sperm proteome, the DmSP3 contains a total of 3176 proteins, and we provide the first label-free quantitation of the sperm proteome for 2125 proteins. The top 20 most abundant proteins included the structural elements α- and β-tubulins and sperm leucyl-aminopeptidases. Both gene content and protein abundance were significantly reduced on the X chromosome, consistent with prior genomic studies of X chromosome evolution. We identified 9 of the 16 Y-linked proteins, including known testis-specific male fertility factors. We also identified almost one-half of known Drosophila ribosomal proteins in the DmSP3. The role of this subset of ribosomal proteins in sperm is unknown. Surprisingly, our expanded sperm proteome also identified 122 seminal fluid proteins (Sfps), proteins originally identified in the accessory glands. We show that a significant fraction of 'sperm-associated Sfps' are recalcitrant to concentrated salt and detergent treatments, suggesting this subclass of Sfps are expressed in testes and may have additional functions in sperm, per se. Overall, our results add to a growing landscape of both sperm and seminal fluid protein biology and in particular provides quantitative evidence at the protein level for prior findings supporting the meiotic sex-chromosome inactivation model for male-specific gene and X chromosome evolution.

Research gaps and new insights in the evolution of Drosophila seminal fluid proteins

While the striking effects of seminal fluid proteins (SFPs) on females are fairly conserved among Diptera, most SFPs lack detectable homologues among the SFP repertoires of phylogenetically distant species. How such a rapidly changing proteome conserves functions across taxa is a fascinating question. However, this and other pivotal aspects of SFPs' evolution remain elusive because discoveries on these proteins have been mainly restricted to the model Drosophila melanogaster. Here, we provide an overview of the current knowledge on the inter-specific divergence of the SFP repertoire in Drosophila and compile the increasing amount of relevant genomic information from multiple species. Capitalizing on the accumulated knowledge in D. melanogaster, we present novel sets of high-confidence SFP candidates and transcription factors presumptively involved in regulating the expression of SFPs. We also address open questions by performing comparative genomic analyses that failed to support the existence of many conserved SFPs shared by most dipterans and indicated that gene co-option is the most frequent mechanism accounting for the origin of Drosophila SFP-coding genes. We hope our update establishes a starting point to integrate further data and thus widen the understanding of the intricate evolution of these proteins.

When is Offspring Viability Fitness a Measure of Paternal Fitness and When is it not?

We use population genetics to model the evolution of a gene with an indirect effect owing to paternal care and with a second pleiotropic, direct effect on offspring viability. We use the model to illustrate how the common empirical practice of considering offspring viability as a component of parent fitness can confound a gene's direct and indirect fitness effects. We investigate when this confounding results in a distorted picture of overall evolution and when it does not. We find that the practice has no effect on mean fitness, W, but it does have an effect on the dynamics of gene frequency change, ∆q. We also find that, for some regions of parameter space associated with fitness trade-offs, the distortion is not only quantitative but also qualitative, obscuring the direction of gene frequency change. Because it affects the evolutionary dynamics, it also affects the expected amount of genetic variation at mutation-selection balance, an important consideration in molecular evolution. We discuss empirical techniques for separating direct from indirect effects and how field studies measuring the value of male paternal care might be improved by using them.

OUP accepted manuscript

Many animal species are haplodiploid: their fertilized eggs develop into diploid females and their unfertilized eggs develop into haploid males. The unique genetic features of haplodiploidy raise the prospect that these systems can be used to disentangle the population genetic consequences of haploid and diploid selection. To this end, sex-specific reproductive genes are of particular interest because, while they are shared within the same genome, they consistently experience selection in different ploidal environments. However, other features of these genes, including sex-specific expression and putative involvement in postcopulatory sexual selection, are potentially confounding factors because they may also impact the efficacy of selection asymmetrically between the sexes. Thus, to properly interpret evolutionary genomic patterns, it is necessary to generate a null expectation for the relative amount of polymorphism and divergence we expect to observe among sex-specific genes in haplodiploid species, given differences in ploidal environment, sex-limited expression, and their potential role in sexual selection. Here, we derive the theoretical expectation for the rate of evolution of sex-specific genes in haplodiploid species, under the assumption that they experience the same selective environment as genes expressed in both sexes. We find that the null expectation is that reproductive genes evolve more rapidly than constitutively expressed genes in haplodiploid genomes. However, despite the aforementioned differences, the null expectation does not differ between male- and female-specific reproductive genes, when assuming additivity. Our theoretical results provide an important baseline expectation that should be used in molecular evolution studies comparing rates of evolution among classes of genes in haplodiploid species.

OUP accepted manuscript

Many species have separate haploid and diploid phases. Theory predicts that each phase should experience the effects of evolutionary forces (like selection) differently. In the haploid phase, all fitness-affecting alleles are exposed to selection, whereas in the diploid phase, those same alleles can be masked by homologous alleles. This predicts that selection acting on genes expressed in haploids should be more effective than diploid-biased genes. Unfortunately, in arrhenotokous species, this prediction can be confounded with the effects of sex-specific expression, as haploids are usually reproductive males. Theory posits that, when accounting for ploidal- and sex-specific expression, selection should be equally efficient on haploid- and diploid-biased genes relative to constitutive genes. Here, we used a multiomic approach in honey bees to quantify the evolutionary rates of haploid-biased genes and test the relative effects of sexual- and haploid-expression on molecular evolution. We found that 16% of the honey bee’s protein-coding genome is highly expressed in haploid tissue. When accounting for ploidy and sex, haploid- and diploid-biased genes evolve at a lower rate than expected, indicating that they experience strong negative selection. However, the rate of molecular evolution of haploid-biased genes was higher than diploid-based genes. Genes associated with sperm storage are a clear exception to this trend with evidence of strong positive selection. Our results provide an important empirical test of theory outlining how selection acts on genes expressed in arrhenotokous species. We propose the haploid life history stage affects genome-wide patterns of diversity and divergence because of both sexual and haploid selection.

Faster Rates of Molecular Sequence Evolution in Reproduction-Related Genes and in Species with Hypodermic Sperm Morphologies

Sexual selection drives the evolution of many striking behaviors and morphologies and should leave signatures of selection at loci underlying these phenotypes. However, although loci thought to be under sexual selection often evolve rapidly, few studies have contrasted rates of molecular sequence evolution at such loci across lineages with different sexual selection contexts. Furthermore, work has focused on separate sexed animals, neglecting alternative sexual systems. We investigate rates of molecular sequence evolution in hermaphroditic flatworms of the genus Macrostomum. Specifically, we compare species that exhibit contrasting sperm morphologies, strongly associated with multiple convergent shifts in the mating strategy, reflecting different sexual selection contexts. Species donating and receiving sperm in every mating have sperm with bristles, likely to prevent sperm removal. Meanwhile, species that hypodermically inject sperm lack bristles, potentially as an adaptation to the environment experienced by hypodermic sperm. Combining functional annotations from the model, Macrostomum lignano, with transcriptomes from 93 congeners, we find genus-wide faster sequence evolution in reproduction-related versus ubiquitously expressed genes, consistent with stronger sexual selection on the former. Additionally, species with hypodermic sperm morphologies had elevated molecular sequence evolution, regardless of a gene's functional annotation. These genome-wide patterns suggest reduced selection efficiency following shifts to hypodermic mating, possibly due to higher selfing rates in these species. Moreover, we find little evidence for convergent amino acid changes across species. Our work not only shows that reproduction-related genes evolve rapidly also in hermaphroditic animals, but also that well-replicated contrasts of different sexual selection contexts can reveal underappreciated genome-wide effects.

Genomic Signatures of Sexual Selection on Pollen-Expressed Genes in Arabis alpina

Fertilization in angiosperms involves the germination of pollen on the stigma, followed by the extrusion of a pollen tube that elongates through the style and delivers two sperm cells to the embryo sac. Sexual selection could occur throughout this process when male gametophytes compete for fertilization. The strength of sexual selection during pollen competition should be affected by the number of genotypes deposited on the stigma. As increased self-fertilization reduces the number of mating partners, and the genetic diversity and heterozygosity of populations, it should thereby reduce the intensity of sexual selection during pollen competition. Despite the prevalence of mating system shifts, few studies have directly compared the molecular signatures of sexual selection during pollen competition in populations with different mating systems. Here we analyzed whole-genome sequences from natural populations of Arabis alpina, a species showing mating system variation across its distribution, to test whether shifts from cross- to self-fertilization result in molecular signatures consistent with sexual selection on genes involved in pollen competition. We found evidence for efficient purifying selection on genes expressed in vegetative pollen, and overall weaker selection on sperm-expressed genes. This pattern was robust when controlling for gene expression level and specificity. In agreement with the expectation that sexual selection intensifies under cross-fertilization, we found that the efficacy of purifying selection on male gametophyte-expressed genes was significantly stronger in genetically more diverse and outbred populations. Our results show that intra-sexual competition shapes the evolution of pollen-expressed genes, and that its strength fades with increasing self-fertilization rates.

Maternal-Offspring Interactions: Reciprocally Coevolving Social Environments

Maternal-zygotic co-evolution is one of the most common examples of indirect genetic effects. I investigate how maternal-zygotic gene interactions affect rates of evolution and adaptation. Using comparably parameterized population genetic models, I compare evolution to an abiotic environment with genotype-by-environment interaction (G × E) to evolution to a maternal environment with offspring genotype-by-maternal environment interaction (G × G_maternal). There are strong parallels between the 2 models in the components of fitness variance but they differ in their rates of evolution measured in terms of ∆p, gene frequency change, or of ∆W, change in mean fitness. The Price Equation is used to partition ∆W into 2 components, one owing to the genetic variance in fitness by natural selection and a second owing to change in environment. Adaptive evolution is faster in the 2-locus model with G × G_maternal with free recombination, than it is in the 1-locus model with G × E, because in the former the maternal genetic environment coevolves with the zygotic phenotype adapting to it. I discuss the relevance of these findings for the evolution of genes with indirect genetic effects.

Single-nucleus transcriptomes reveal evolutionary and functional properties of cell types in the Drosophila accessory gland

Many traits responsible for male reproduction evolve quickly, including gene expression phenotypes in germline and somatic male reproductive tissues. Rapid male evolution in polyandrous species is thought to be driven by competition among males for fertilizations and conflicts between male and female fitness interests that manifest in post-copulatory phenotypes. In Drosophila, seminal fluid proteins secreted by three major cell types of the male accessory gland and ejaculatory duct are required for female sperm storage and use, and influence female post-copulatory traits. Recent work has shown that these cell types have overlapping but distinct effects on female post-copulatory biology, yet relatively little is known about their evolutionary properties. Here we use single-nucleus RNA-Seq of the accessory gland and ejaculatory duct from Drosophila melanogaster and two closely related species to comprehensively describe the cell diversity of these tissues and their transcriptome evolution for the first time. We find that seminal fluid transcripts are strongly partitioned across the major cell types, and expression of many other genes additionally define each cell type. We also report previously undocumented diversity in main cells. Transcriptome divergence was found to be heterogeneous across cell types and lineages, revealing a complex evolutionary process. Furthermore, protein adaptation varied across cell types, with potential consequences for our understanding of selection on male post-copulatory traits.

Viviparity and habitat restrictions may influence the evolution of male reproductive genes in tsetse fly (Glossina) species

BACKGROUND

Glossina species (tsetse flies), the sole vectors of African trypanosomes, maintained along their long evolutionary history a unique reproductive strategy, adenotrophic viviparity. Viviparity reduces their reproductive rate and, as such, imposes strong selective pressures on males for reproductive success. These species live in sub-Saharan Africa, where the distributions of the main sub-genera Fusca, Morsitans, and Palpalis are restricted to forest, savannah, and riverine habitats, respectively. Here we aim at identifying the evolutionary patterns of the male reproductive genes of six species belonging to these three main sub-genera. We then interpreted the different patterns we found across the species in the light of viviparity and the specific habitat restrictions, which are known to shape reproductive behavior.

RESULTS

We used a comparative genomic approach to build consensus evolutionary trees that portray the selective pressure acting on the male reproductive genes in these lineages. Such trees reflect the long and divergent demographic history that led to an allopatric distribution of the Fusca, Morsitans, and Palpalis species groups. A dataset of over 1700 male reproductive genes remained conserved over the long evolutionary time scale (estimated at 26.7 million years) across the genomes of the six species. We suggest that this conservation may result from strong functional selective pressure on the male imposed by viviparity. It is noteworthy that more than half of these conserved genes are novel sequences that are unique to the Glossina genus and are candidates for selection in the different lineages.

CONCLUSIONS

Tsetse flies represent a model to interpret the evolution and differentiation of male reproductive biology under different, but complementary, perspectives. In the light of viviparity, we must take into account that these genes are constrained by a post-fertilization arena for genomic conflicts created by viviparity and absent in ovipositing species. This constraint implies a continuous antagonistic co-evolution between the parental genomes, thus accelerating inter-population post-zygotic isolation and, ultimately, favoring speciation. Ecological restrictions that affect reproductive behavior may further shape such antagonistic co-evolution.

Vulture Genomes Reveal Molecular Adaptations Underlying Obligate Scavenging and Low Levels of Genetic Diversity

Obligate scavenging on the dead and decaying animal matter is a rare dietary specialization that in extant vertebrates is restricted to vultures. These birds perform essential ecological services, yet many vulture species have undergone recent steep population declines and are now endangered. To test for molecular adaptations underlying obligate scavenging in vultures, and to assess whether genomic features might have contributed to their population declines, we generated high-quality genomes of the Himalayan and bearded vultures, representing both independent origins of scavenging within the Accipitridae, alongside a sister taxon, the upland buzzard. By comparing our data to published sequences from other birds, we show that the evolution of obligate scavenging in vultures has been accompanied by widespread positive selection acting on genes underlying gastric acid production, and immunity. Moreover, we find evidence of parallel molecular evolution, with amino acid replacements shared among divergent lineages of these scavengers. Our genome-wide screens also reveal that both the Himalayan and bearded vultures exhibit low levels of genetic diversity, equating to around a half of the mean genetic diversity of other bird genomes examined. However, demographic reconstructions indicate that population declines began at around the Last Glacial Maximum, predating the well-documented dramatic declines of the past three decades. Taken together, our genomic analyses imply that vultures harbor unique adaptations for processing carrion, but that modern populations are genetically depauperate and thus especially vulnerable to further genetic erosion through anthropogenic activities.

Variation in selective constraints along the Plasmodium life cycle

Plasmodium parasites, the cause of malaria, have a complex life cycle, infecting alternatively vertebrate hosts and female Anopheles mosquitoes and undergoing intra- and extra-cellular development in several organs of these hosts. Most of the ~5000 protein-coding genes present in Plasmodium genomes are only expressed at specific life stages, and different genes might therefore be subject to different selective pressures depending on the biological activity of the parasite and its microenvironment at this point in development. Here, we estimate the selective constraints on the protein-coding sequences of all annotated genes of rodent and primate Plasmodium parasites and, using data from scRNA-seq experiments spanning many developmental stages, analyze their variation with regard to when these genes are expressed in the parasite life cycle. Our study reveals extensive variation in selective constraints throughout the parasites' development and highlights stages that are evolving more rapidly than others. These findings provide novel insights into the biology of these parasites and could provide important information to develop better treatment strategies or vaccines against these medically-important organisms.

Nonadaptive molecular evolution of seminal fluid proteins in Drosophila

Seminal fluid proteins (SFPs) are a group of reproductive proteins that are among the most evolutionarily divergent known. As SFPs can impact male and female fitness, these proteins have been proposed to evolve under postcopulatory sexual selection (PCSS). However, the fast change of the SFPs can also result from nonadaptive evolution, and the extent to which selective constraints prevent SFPs rapid evolution remains unknown. Using intra‐ and interspecific sequence information, along with genomics and functional data, we examine the molecular evolution of approximately 300 SFPs in Drosophila. We found that 50–57% of the SFP genes, depending on the population examined, are evolving under relaxed selection. Only 7–12% showed evidence of positive selection, with no evidence supporting other forms of PCSS, and 35–37% of the SFP genes were selectively constrained. Further, despite associations of positive selection with gene location on the X chromosome and protease activity, the analysis of additional genomic and functional features revealed their lack of influence on SFPs evolving under positive selection. Our results highlight a lack of sufficient evidence to claim that most SFPs are driven to evolve rapidly by PCSS while identifying genomic and functional attributes that influence different modes of SFPs evolution.

Experimental evolution supports signatures of sexual selection in genomic divergence

Comparative genomics has contributed to the growing evidence that sexual selection is an important component of evolutionary divergence and speciation. Divergence by sexual selection is implicated in faster rates of divergence of the X chromosome and of genes thought to underlie sexually selected traits, including genes that are sex biased in expression. However, accurately inferring the relative importance of complex and interacting forms of natural selection, demography, and neutral processes that occurred in the evolutionary past is challenging. Experimental evolution provides an opportunity to apply controlled treatments for multiple generations and examine the consequent genomic divergence. Here, we altered sexual selection intensity, elevating sexual selection in polyandrous lines and eliminating it in monogamous lines, and examined patterns of allele frequency divergence in the genome of Drosophila pseudoobscura after more than 160 generations of experimental evolution. Divergence is not uniform across the genome but concentrated in "islands," many of which contain candidate genes implicated in mating behaviors and other sexually selected phenotypes. These are more often seen on the X chromosome, which also shows greater divergence in F ST than neutral expectations. There are characteristic signatures of selection seen in these regions, with lower diversity on the X chromosome than the autosomes, and differences in diversity on the autosomes between selection regimes. Reduced Tajima's D within some of the divergent regions may imply that selective sweeps have occurred, despite considerable recombination. These changes are associated with both differential gene expression between the lines and sex-biased gene expression within the lines. Our results are very similar to those thought to implicate sexual selection in divergence between species and natural populations, and hence provide experimental support for the likely role of sexual selection in driving such types of genetic divergence, but also illustrate how variable outcomes can be for different genomic regions.

Gene expression and alternative splicing dynamics are perturbed in female head transcriptomes following heterospecific copulation

BACKGROUND

Despite the growing interest in the female side of copulatory interactions, the roles played by differential expression and alternative splicing mechanisms of pre-RNA on tissues outside of the reproductive tract have remained largely unknown. Here we addressed these questions in the context of con- vs heterospecific matings between Drosophila mojavensis and its sister species, D. arizonae. We analyzed transcriptional responses in female heads using an integrated investigation of genome-wide patterns of gene expression, including differential expression (DE), alternative splicing (AS) and intron retention (IR).

RESULTS

Our results indicated that early transcriptional responses were largely congruent between con- and heterospecific matings but are substantially perturbed over time. Conspecific matings induced functional pathways related to amino acid balance previously associated with the brain's physiology and female postmating behavior. Heterospecific matings often failed to activate regulation of some of these genes and induced expression of additional genes when compared with those of conspecifically-mated females. These mechanisms showed functional specializations with DE genes mostly linked to pathways of proteolysis and nutrient homeostasis, while AS genes were more related to photoreception and muscle assembly pathways. IR seems to play a more general role in DE regulation during the female postmating response.

CONCLUSIONS

We provide evidence showing that AS genes substantially perturbed by heterospecific matings in female heads evolve at slower evolutionary rates than the genome background. However, DE genes evolve at evolutionary rates similar, or even higher, than those of male reproductive genes, which highlights their potential role in sexual selection and the evolution of reproductive barriers.

Rapid divergence of the male reproductive proteins in the Drosophila dunni group and implications for postmating incompatibilities between species

Proteins involved in post-copulatory interactions between males and females are among the fastest evolving genes in many species, usually attributed to their involvement in reproductive conflict. As a result, these proteins are thought to often be involved in the formation of postmating-prezygotic incompatibilities between species. The Drosophila dunni subgroup consists of a dozen recently diverged species found across the Caribbean islands with varying levels of hybrid incompatibility. We performed experimental crosses between species in the dunni group and see some evidence of hybrid incompatibilities. We also find evidence of reduced survival following hybrid mating, likely due to postmating-prezygotic incompatibilities. We assessed rates of evolution between these species genomes and find evidence of rapid evolution and divergence of some reproductive proteins, specifically the seminal fluid proteins. This work suggests the rapid evolution of seminal fluid proteins may be associated with postmating-prezygotic isolation, which acts as a barrier for gene flow between even the most closely related species.

Global climate change, diet, and the complex relationship between human host and microbiome: Towards an integrated picture

Dietary changes can alter the human microbiome with potential detrimental consequences for health. Given that environment, health, and evolution are interconnected, we ask: Could diet-driven microbiome perturbations have consequences that extend beyond their immediate impact on human health? We address this question in the context of the urgent health challenges posed by global climate change. Drawing on recent studies, we propose that not only can diet-driven microbiome changes lead to dysbiosis, they can also shape life-history traits and fuel human evolution. We posit that dietary shifts prompt mismatched microbiome-host genetics configurations that modulate human longevity and reproductive success. These mismatches can also induce a heritable intra-holobiont stress response, which encourages the holobiont to re-establish equilibrium within the changed nutritional environment. Thus, while mismatches between climate change-related genetic and epigenetic configurations within the holobiont increase the risk and severity of diseases, they may also affect life-history traits and facilitate adaptive responses. These propositions form a framework that can help systematize and address climate-related dietary challenges for policy and health interventions.

Male fertility in Arabidopsis requires active DNA demethylation of genes that control pollen tube function

Active DNA demethylation is required for sexual reproduction in plants but the molecular determinants underlying this epigenetic control are not known. Here, we show in Arabidopsis thaliana that the DNA glycosylases DEMETER (DME) and REPRESSOR OF SILENCING 1 (ROS1) act semi-redundantly in the vegetative cell of pollen to demethylate DNA and ensure proper pollen tube progression. Moreover, we identify six pollen-specific genes with increased DNA methylation as well as reduced expression in dme and dme;ros1. We further show that for four of these genes, reinstalling their expression individually in mutant pollen is sufficient to improve male fertility. Our findings demonstrate an essential role of active DNA demethylation in regulating genes involved in pollen function.

Reproductive Proteins Evolve Faster Than Non-reproductive Proteins Among Solanum Species

Elevated rates of evolution in reproductive proteins are commonly observed in animal species, and are thought to be driven by the action of sexual selection and sexual conflict acting specifically on reproductive traits. Whether similar patterns are broadly observed in other biological groups is equivocal. Here, we examine patterns of protein divergence among wild tomato species (Solanum section Lycopersicon), to understand forces shaping the evolution of reproductive genes in this diverse, rapidly evolving plant clade. By comparing rates of molecular evolution among loci expressed in reproductive and non-reproductive tissues, our aims were to test if: (a) reproductive-specific loci evolve more rapidly, on average, than non-reproductive loci; (b) 'male'-specific loci evolve at different rates than 'female'-specific loci; (c) genes expressed exclusively in gametophytic (haploid) tissue evolve differently from genes expressed in sporophytic (diploid) tissue or in both tissue types; and (d) mating system variation (a potential proxy for the expected strength of sexual selection and/or sexual conflict) affects patterns of protein evolution. We observed elevated evolutionary rates in reproductive proteins. However, this pattern was most evident for female- rather than male-specific loci, both broadly and for individual loci inferred to be positively selected. These elevated rates might be facilitated by greater tissue-specificity of reproductive proteins, as faster rates were also associated with more narrow expression domains. In contrast, we found little evidence that evolutionary rates are consistently different in loci experiencing haploid selection (gametophytic-exclusive loci), or in lineages with quantitatively different mating systems. Overall while reproductive protein evolution is generally elevated in this diverse plant group, some specific patterns of evolution are more complex than those reported in other (largely animal) systems, and include a more prominent role for female-specific loci among adaptively evolving genes.

Satyrization in Drosophila fruitflies

The satyr of Greek mythology was half‐man, half‐goat, with an animal persona signifying immoderate sexual appetites. In biology, satyrization is the disruption of reproduction in matings between closely related species. Interestingly, its effects are often reciprocally asymmetric, manifesting more strongly in one direction of heterospecific mating than the other. Heterospecific matings are well known to result in female fitness costs due to the production of sterile or inviable hybrid offspring and can also occur due to reduced female sexual receptivity, lowering the likelihood of any subsequent conspecific matings. Here we investigated the costs and mechanisms of satyrization in the Drosophila melanogaster species subgroup of fruitflies. The results showed that D. simulans females experienced higher fitness costs from a loss of remating opportunities due to significantly reduced post‐mating sexual receptivity than did D. melanogaster females, as a result of reciprocal heterospecific matings. Reciprocal tests of the effects of male reproductive accessory gland protein (Acp) injections on female receptivity in pairwise comparisons between D. melanogaster and five other species within the melanogaster species subgroup revealed significant post‐mating receptivity asymmetries. This was due to variation in the effects of heterospecific Acps within species with which D. melanogaster can mate, and significant but nonasymmetric Acp effects in species with which it cannot. We conclude that asymmetric satyrization due to post‐mating effects of Acps may be common among diverging and hybridising species. The findings are of interest in understanding the evolution of reproductive isolation and species divergence.