The genomic footprint of sexual conflict

Sex-Specific Dominance of Gene Expression in Seed Beetles

When different alleles are favored in different environments, dominance reversal where alternate alleles are dominant in the environment in which they are favored can generate net balancing selection. The sexes represent two distinct genetic environments and sexually antagonistic (SA) selection can maintain genetic variation, especially when the alleles involved show sex-specific dominance. Sexual dimorphism in gene expression is pervasive and has been suggested to result from SA selection. Yet, whether gene-regulatory variation shows sex-specific dominance is poorly understood. We tested for sex-specific dominance in gene expression using three crosses between homozygous lines derived from a population of a seed beetle, where a previous study documented a signal of dominance reversal for fitness between the sexes. Overall, we found that the dominance effects of variants affecting gene expression were positively correlated between the sexes (r = 0.33 to 0.44). Yet, 586 transcripts showed significant differences in dominance between the sexes. Sex-specific dominance was significantly more common in transcripts with more sex-biased expression, in two of three of our crosses. Among transcripts showing sex-specific dominance, lesser sexual dimorphism in gene expression among heterozygotes was somewhat more common than greater. Gene ontology enrichment analyses showed that functional categories associated with known SA phenotypes in Callosobruchus maculatus were overrepresented among transcripts with sex-specific dominance, including genes involved in metabolic processes and the target-of-rapamycin pathway. Our results support the suggestion that sex-specific dominance of regulatory variants contributes to the maintenance of genetic variation in fitness mediated by SA selection in this species.

CYP4BN4v7 regulates the population density dependent oocyte maturity rate in bean beetles

The bean beetle (Callosobruchus maculatus) clearly exhibits population density-dependent polymorphism. Cytochrome P450 (CYP) is involved in many aspects of the physiological activities of insects. However, the role of CYP in population density-dependent polymorphisms remains unknown. The terminal oocyte maturity rate of high-population-density individuals (H) was faster than that of low-population-density individuals (L). A total of 56 CYP-like genes were identified from transcriptomic and genomic data, including seven clan 2 CYP-like genes, seven mitochondrial CYP-like genes, 19 clan 3 CYP-like genes, and 23 clan 4 CYP-like genes. Gene duplication might occur in CYP9Z4-like, CYP345A1-like, CYP345A2-like, CYP349A1-like, CYP349A2-like, and CYP4BN4-like. Thirteen and two CYP-like genes were up-regulated and down-regulated, respectively, in H. Among these CYP-like genes, CYP4BN4v7-like was the most abundant CYP. CYP4BN4v7-like was more highly expressed in the head than in the thorax and abdomen. Its mRNA levels in the head, thorax, and abdomen were greater in H than in L. After RNA interference decreased its mRNA level, the terminal oocyte maturity rate decreased. Moreover, the expression level of insulin-like peptide 1 (ILP1), which plays a vital role in regulating terminal oocyte development, decreased in the head. In conclusion, CYP4BN4v7-like modulated the population density-dependent terminal oocyte maturity rate by regulating the expression of ILP1.

Nonlinear transcriptomic responses to compounded environmental changes across temperature and resources in a pest beetle, Callosobruchus maculatus (Coleoptera: Chrysomelidae)

Many species are experiencing drastic and multidimensional changes to their environment due to anthropogenic events. These multidimensional changes may act nonadditively on physiological and life history responses, and thus may not be predicted by responses to single dimensional environmental changes. Therefore, work is needed to understand species' responses to multiple aspects of change. We used whole-transcriptomic RNA-Sequencing and life history assays to uncover responses to singly-applied shifts in resource or temperature environmental dimensions, in comparison to combined, multidimensional change, in the crop pest seed beetle, Callosobruchus maculatus. We found that multidimensional change caused larger fecundity, developmental period and offspring viability life history changes than predicted by additive effects of 1-dimensional changes. In addition, there was little overlap between genes differentially expressed under multidimensional treatment versus under altered resource or temperature conditions alone. Moreover, 115 genes exhibited significant resource × temperature interaction effects on expression, including those involved in energy metabolism, detoxification, and enhanced formation of cuticle structural components. We conclude that single dimensional changes alone cannot determine life history and transcriptomic responses to multidimensional environmental change. These results highlight the importance of studying multidimensional environmental change for understanding the molecular and phenotypic responses that may allow organisms including insects to rapidly adapt simultaneously to multiple aspects of environmental change.

A role of epigenetic mechanisms in regulating female reproductive responses to temperature in a pest beetle

Many species are threatened by climate change and must rapidly respond to survive in changing environments. Epigenetic modifications, such as DNA methylation, can facilitate plastic responses by regulating gene expression in response to environmental cues. Understanding epigenetic responses is therefore essential for predicting species' ability to rapidly adapt in the context of global environmental change. Here, we investigated the functional significance of different methylation-associated cellular processes on temperature-dependent life history in seed beetles, Callosobruchus maculatus Fabricius 1775 (Coleoptera: Bruchidae). We assessed changes under thermal stress in (1) DNA methyltransferase (Dnmt1 and Dnmt2) expression levels, (2) genome-wide methylation and (3) reproductive performance, with (2) and (3) following treatment with 3-aminobenzamide (3AB) and zebularine (Zeb) over two generations. These drugs are well-documented to alter DNA methylation across the tree of life. We found that Dnmt1 and Dnmt2 were expressed throughout the body in males and females, but were highly expressed in females compared with males and exhibited temperature dependence. However, whole-genome methylation did not significantly vary with temperature, and only marginally or inconclusively with drug treatment. Both 3AB and Zeb led to profound temperature-dependent shifts in female reproductive life history trade-off allocation, often increasing fitness compared with control beetles. Mismatch between magnitude of treatment effects on DNA methylation versus life history effects suggest potential of 3AB and Zeb to alter reproductive trade-offs via changes in DNA repair and recycling processes, rather than or in addition to (subtle) changes in DNA methylation. Together, our results suggest that epigenetic mechanisms relating to Dnmt expression, DNA repair and recycling pathways, and possibly DNA methylation, are strongly implicated in modulating insect life history trade-offs in response to temperature change.

Pesticide exposure triggers sex-specific inter- and transgenerational effects conditioned by past sexual selection

Environmental variation often induces plastic responses in organisms that can trigger changes in subsequent generations through non-genetic inheritance mechanisms. Such transgenerational plasticity thus consists of environmentally induced non-random phenotypic modifications that are transmitted through generations. Transgenerational effects may vary according to the sex of the organism experiencing the environmental perturbation, the sex of their descendants or both, but whether they are affected by past sexual selection is unknown. Here, we use experimental evolution on an insect model system to conduct a first test of the involvement of sexual selection history in shaping transgenerational plasticity in the face of rapid environmental change (exposure to pesticide). We manipulated evolutionary history in terms of the intensity of sexual selection for over 80 generations before exposing individuals to the toxicant. We found that sexual selection history constrained adaptation under rapid environmental change. We also detected inter- and transgenerational effects of pesticide exposure in the form of increased fitness and longevity. These cross-generational influences of toxicants were sex dependent (they affected only male descendants), and intergenerational, but not transgenerational, plasticity was modulated by sexual selection history. Our results highlight the complexity of intra-, inter- and transgenerational influences of past selection and environmental stress on phenotypic expression.

High-quality reference genome of cowpea beetle Callosobruchus maculatus

Callosobruchus maculatus is one of the most competitive stored grain pests, which causes a great loss to agricultural economy. However, due to an inadequacy of high-quality reference genome, the molecular mechanisms for olfactory and hypoxic adaptations to stored environments are unknown and require to be revealed urgently, which will contribute to the detection and prevention of the invasive pests C. maculatus. Here, we presented a high-quality chromosome-level genome of C. maculatus based on Illumina, Nanopore and Hi-C sequencing data. The total size was 1.2 Gb, and 65.17% (797.47 Mb) of it was identified to be repeat sequences. Among assembled chromosomes, chromosome 10 was considered the X chromosome according to the evidence of reads coverage and homologous genes among species. The current version of high-quality genome provides preferable data resources for the adaptive evolution research of C. maculatus.

Comparative Genomic Analysis of the Pattern of Evolution of Male and Female Reproductive Proteins in Seed Beetles

Male seminal fluid proteins often show signs of positive selection and divergent evolution, believed to reflect male-female coevolution. Yet, our understanding of the predicted concerted evolution of seminal fluid proteins and female reproductive proteins is limited. We sequenced, assembled, and annotated the genome of two species of seed beetles allowing a comparative analysis of four closely related species of these herbivorous insects. We compare the general pattern of evolution in genes encoding seminal fluid proteins and female reproductive proteins with those in digestive protein genes and well-conserved reference genes. We found that female reproductive proteins showed an overall ratio of nonsynonymous to synonymous substitutions (ω) similar to that of conserved genes, while seminal fluid proteins and digestive proteins exhibited higher overall ω values. Further, seminal fluid proteins and digestive proteins showed a higher proportion of sites putatively under positive selection, and explicit tests showed no difference in relaxed selection between protein types. Evolutionary rate covariation analyses showed that evolutionary rates among seminal fluid proteins were on average more closely correlated with those in female reproductive proteins than with either digestive or conserved genes. Gene expression showed the expected negative covariation with ω values, except for male-biased genes where this negative relationship was reversed. In conclusion, seminal fluid proteins showed relatively rapid evolution and signs of positive selection. In contrast, female reproductive proteins evolved at a lower rate under selective constraints, on par with genes known to be well conserved. Although our findings provide support for concerted evolution of seminal fluid proteins and female reproductive proteins, they also suggest that these two classes of proteins evolve under partly distinct selective regimes.

On the resolution of sexual conflict over shared traits

Anisogamy, different-sized male and female gametes, sits at the heart of sexual selection and conflict between the sexes. Sperm producers (males) and egg producers (females) of the same species generally share most, if not all, of the same genome, but selection frequently favours different trait values in each sex for traits common to both. The extent to which this conflict might be resolved, and the potential mechanisms by which this can occur, have been widely debated. Here, we summarize recent findings and emphasize that once the sexes evolve, sexual selection is ongoing, and therefore new conflict is always possible. In addition, sexual conflict is largely a multivariate problem, involving trait combinations underpinned by networks of interconnected genes. Although these complexities can hinder conflict resolution, they also provide multiple possible routes to decouple male and female phenotypes and permit sex-specific evolution. Finally, we highlight difficulty in the study of sexual conflict over shared traits and promising directions for future research.

Multivariate genetic architecture reveals testosterone-driven sexual antagonism in contemporary humans

Sex difference (SD) is ubiquitous in humans despite shared genetic architecture (SGA) between the sexes. A univariate approach, i.e., studying SD in single traits by estimating genetic correlation, does not provide a complete biological overview, because traits are not independent and are genetically correlated. The multivariate genetic architecture between the sexes can be summarized by estimating the additive genetic (co)variance across shared traits, which, apart from the cross-trait and cross-sex covariances, also includes the cross-sex-cross-trait covariances, e.g., between height in males and weight in females. Using such a multivariate approach, we investigated SD in the genetic architecture of 12 anthropometric, fat depositional, and sex-hormonal phenotypes. We uncovered sexual antagonism (SA) in the cross-sex-cross-trait covariances in humans, most prominently between testosterone and the anthropometric traits - a trend similar to phenotypic correlations. 27% of such cross-sex-cross-trait covariances were of opposite sign, contributing to asymmetry in the SGA. Intriguingly, using multivariate evolutionary simulations, we observed that the SGA acts as a genetic constraint to the evolution of SD in humans only when selection is sexually antagonistic and not concordant. Remarkably, we found that the lifetime reproductive success in both the sexes shows a positive genetic correlation with anthropometric traits, but not with testosterone. Moreover, we demonstrated that genetic variance is depleted along multivariate trait combinations in both the sexes but in different directions, suggesting absolute genetic constraint to evolution. Our results indicate that testosterone drives SA in contemporary humans and emphasize the necessity and significance of using a multivariate framework in studying SD.

Sexual conflict drive in the rapid evolution of new gametogenesis genes

The evolutionary forces underlying the rapid evolution in sequences and functions of new genes remain a mystery. Adaptation by natural selection explains the evolution of some new genes. However, many new genes perform sex-biased functions that have rapidly evolved over short evolutionary time scales, suggesting that new gene evolution may often be driven by conflicting selective pressures on males and females. It is well established that such sexual conflict (SC) plays a central role in maintaining phenotypic and genetic variation within populations, but the role of SC in driving new gene evolution remains essentially unknown. This review explores the connections between SC and new gene evolution through discussions of the concept of SC, the phenotypic and genetic signatures of SC in evolving populations, and the molecular mechanisms by which SC could drive the evolution of new genes. We synthesize recent work in this area with a discussion of the case of Apollo and Artemis, two extremely young genes (<200,000 years) in Drosophila melanogaster, which offered the first empirical insights into the evolutionary process by which SC could drive the evolution of new genes. These new duplicate genes exhibit the hallmarks of sexually antagonistic selection: rapid DNA and protein sequence evolution, essential sex-specific functions in gametogenesis, and complementary sex-biased expression patterns. Importantly, Apollo is essential for male fitness but detrimental to female fitness, while Artemis is essential for female fitness but detrimental to male fitness. These sexually antagonistic fitness effects and complementary changes to expression, sequence, and function suggest that these duplicates were selected for mitigating SC, but that SC has not been fully resolved. Finally, we propose Sexual Conflict Drive as a self-driven model to interpret the rapid evolution of new genes, explain the potential for SC and sexually antagonistic selection to contribute to long-term evolution, and suggest its utility for understanding the rapid evolution of new genes in gametogenesis.

The roles of gene duplications in the dynamics of evolutionary conflicts

Evolutionary conflicts occur when there is antagonistic selection between different individuals of the same or different species, life stages or between levels of biological organization. Remarkably, conflicts can occur within species or within genomes. In the dynamics of evolutionary conflicts, gene duplications can play a major role because they can bring very specific changes to the genome: changes in protein dose, the generation of novel paralogues with different functions or expression patterns or the evolution of small antisense RNAs. As we describe here, by having those effects, gene duplication might spark evolutionary conflict or fuel arms race dynamics that takes place during conflicts. Interestingly, gene duplication can also contribute to the resolution of a within-locus evolutionary conflict by partitioning the functions of the gene that is under an evolutionary trade-off. In this review, we focus on intraspecific conflicts, including sexual conflict and illustrate the various roles of gene duplications with a compilation of examples. These examples reveal the level of complexity and the differences in the patterns of gene duplications within genomes under different conflicts. These examples also reveal the gene ontologies involved in conflict and the genomic location of the elements of the conflict. The examples provide a blueprint for the direct study of these conflicts or the exploration of the presence of similar conflicts in other lineages.

Gene redundancy and gene compensation of insulin-like peptides in the oocyte development of bean beetle

Bean beetle (Callosobruchus maculatus) exhibits clear phenotypic plasticity depending on population density; However, the underlying molecular mechanism remains unknown. Compared to low-density individuals, high-density individuals showed a faster terminal oocyte maturity rate. Four insulin-like peptide (ILP) genes were identified in the bean beetle, which had higher expression levels in the head than in the thorax and abdomen. The population density could regulate the expression levels of CmILP1-3, CmILP2-3, and CmILP1 as well as CmILP3 in the head, thorax, and abdomen, respectively. RNA interference results showed that each CmILP could regulate terminal oocyte maturity rate, indicating that there was functional redundancy among CmILPs. Silencing each CmILP could lead to down-regulation of some other CmILPs, however, CmILP3 was up-regulated in the abdomen after silencing CmILP1 or CmILP2. Compared to single gene silencing, silencing CmILP3 with CmILP1 or CmILP2 at the same time led to more serious retardation in oocyte development, suggesting CmILP3 could be up-regulated to functionally compensate for the down-regulation of CmILP1 and CmILP2. In conclusion, population density-dependent plasticity in terminal oocyte maturity rate of bean beetle was regulated by CmILPs, which exhibited gene redundancy and gene compensation.

A chromosome-level assembly of the seed beetle Callosobruchus maculatus genome with annotation of its repetitive elements

Callosobruchus maculatus is a major agricultural pest of legume crops worldwide and an established model system in ecology and evolution. Yet, current molecular biological resources for this species are limited. Here, we employ Hi-C sequencing to generate a greatly improved genome assembly and we annotate its repetitive elements in a dedicated in-depth effort where we manually curate and classify the most abundant unclassified repeat subfamilies. We present a scaffolded chromosome-level assembly, which is 1.01 Gb in total length with 86% being contained within the 9 autosomes and the X chromosome. Repetitive sequences accounted for 70% of the total assembly. DNA transposons covered 18% of the genome, with the most abundant superfamily being Tc1-Mariner (9.75% of the genome). This new chromosome-level genome assembly of C. maculatus will enable future genetic and evolutionary studies not only of this important species but of beetles more generally.

The genomics and evolution of inter-sexual mimicry and female-limited polymorphisms in damselflies

Sex-limited morphs can provide profound insights into the evolution and genomic architecture of complex phenotypes. Inter-sexual mimicry is one particular type of sex-limited polymorphism in which a novel morph resembles the opposite sex. While inter-sexual mimics are known in both sexes and a diverse range of animals, their evolutionary origin is poorly understood. Here, we investigated the genomic basis of female-limited morphs and male mimicry in the common bluetail damselfly. Differential gene expression between morphs has been documented in damselflies, but no causal locus has been previously identified. We found that male mimicry originated in an ancestrally sexually dimorphic lineage in association with multiple structural changes, probably driven by transposable element activity. These changes resulted in ~900 kb of novel genomic content that is partly shared by male mimics in a close relative, indicating that male mimicry is a trans-species polymorphism. More recently, a third morph originated following the translocation of part of the male-mimicry sequence into a genomic position ~3.5 mb apart. We provide evidence of balancing selection maintaining male mimicry, in line with previous field population studies. Our results underscore how structural variants affecting a handful of potentially regulatory genes and morph-specific genes can give rise to novel and complex phenotypic polymorphisms.

Sex-specific trait architecture in a spider with sexual size dimorphism

Sexual dimorphism, or sex-specific trait expression, may evolve when selection favours different optima for the same trait between sexes, that is, under antagonistic selection. Intra-locus sexual conflict exists when the sexually dimorphic trait under antagonistic selection is based on genes shared between sexes. A common assumption is that the presence of sexual-size dimorphism (SSD) indicates that sexual conflict has been, at least partly, resolved via decoupling of the trait architecture between sexes. However, whether and how decoupling of the trait architecture between sexes has been realized often remains unknown. We tested for differences in architecture of adult body size between sexes in a species with extreme SSD, the African hermit spider (Nephilingis cruentata), where adult female body size greatly exceeds that of males. Specifically, we estimated the sex-specific importance of genetic and maternal effects on adult body size among individuals that we laboratory-reared for up to eight generations. Quantitative genetic model estimates indicated that size variation in females is to a larger extent explained by direct genetic effects than by maternal effects, but in males to a larger extent by maternal than by genetic effects. We conclude that this sex-specific body-size architecture enables body-size evolution to proceed much more independently than under a common architecture to both sexes.

Y-Linked Copy Number Polymorphism of Target of Rapamycin Is Associated with Sexual Size Dimorphism in Seed Beetles

The Y chromosome is theorized to facilitate evolution of sexual dimorphism by accumulating sexually antagonistic loci, but empirical support is scarce. Due to the lack of recombination, Y chromosomes are prone to degenerative processes, which poses a constraint on their adaptive potential. Yet, in the seed beetle, Callosobruchus maculatus segregating Y linked variation affects male body size and thereby sexual size dimorphism (SSD). Here, we assemble C. maculatus sex chromosome sequences and identify molecular differences associated with Y-linked SSD variation. The assembled Y chromosome is largely euchromatic and contains over 400 genes, many of which are ampliconic with a mixed autosomal and X chromosome ancestry. Functional annotation suggests that the Y chromosome plays important roles in males beyond primary reproductive functions. Crucially, we find that, besides an autosomal copy of the gene target of rapamycin (TOR), males carry an additional TOR copy on the Y chromosome. TOR is a conserved regulator of growth across taxa, and our results suggest that a Y-linked TOR provides a male specific opportunity to alter body size. A comparison of Y haplotypes associated with male size difference uncovers a copy number variation for TOR, where the haplotype associated with decreased male size, and thereby increased sexual dimorphism, has two additional TOR copies. This suggests that sexual conflict over growth has been mitigated by autosome to Y translocation of TOR followed by gene duplications. Our results reveal that despite of suppressed recombination, the Y chromosome can harbor adaptive potential as a male-limited supergene.

Why is measuring and predicting fitness under genomic conflict so hard?

Genomic conflict between the sexes is caused by differences in the optimal male and female reproductive strategies, and is a major contributor to genetic, phenotypic, and life history variation. While early experimental work appeared to strongly support the sexual conflict paradigm, recent work has produced more ambiguous results. Recent advances in experimental evolution studies combined with theoretical arguments can shed light on why measuring fitness under a conflict is so challenging, including the incidental alteration of mating dynamics, demographic effects, and inherent complexity in what quantity selection maximizes. We stress that non-intuitive results do not necessarily mean the absence of conflict, and follow-up experiments to determine why a priori predictions failed can ultimately teach us more than if they had been confirmed.

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.

Recent progress in understanding the genomic architecture of sexual conflict

Genomic conflict between the sexes over shared traits is widely assumed to be resolved through the evolution of sex-biased expression and the subsequent emergence of sexually dimorphic phenotypes. However, while there is support for a broad relationship between genome-wide patterns of expression level and sexual conflict, recent studies suggest that sex differences in the nature and strength of interactions between loci are instead key to conflict resolution. Furthermore, the advent of new technologies for measuring and perturbing expression means we now have much more power to detect genomic signatures of sexual conflict. Here, we review our current understanding of the genomic architecture of sexual conflict in the light of these new studies and highlight the potential for novel approaches to address outstanding knowledge gaps.

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.

Increased male investment in sperm competition results in reduced maintenance of gametes

Male animals often show higher mutation rates than their female conspecifics. A hypothesis for this male bias is that competition over fertilization of female gametes leads to increased male investment into reproduction at the expense of maintenance and repair, resulting in a trade-off between male success in sperm competition and offspring quality. Here, we provide evidence for this hypothesis by harnessing the power of experimental evolution to study effects of sexual selection on the male germline in the seed beetle Callosobruchus maculatus. We first show that 50 generations of evolution under strong sexual selection, coupled with experimental removal of natural selection, resulted in males that are more successful in sperm competition. We then show that these males produce progeny of lower quality if engaging in sociosexual interactions prior to being challenged to surveil and repair experimentally induced damage in their germline and that the presence of male competitors alone can be enough to elicit this response. We identify 18 candidate genes that showed differential expression in response to the induced germline damage, with several of these previously implicated in processes associated with DNA repair and cellular maintenance. These genes also showed significant expression changes across sociosexual treatments of fathers and predicted the reduction in quality of their offspring, with expression of one gene also being strongly correlated to male sperm competition success. Sex differences in expression of the same 18 genes indicate a substantially higher female investment in germline maintenance. While more work is needed to detail the exact molecular underpinnings of our results, our findings provide rare experimental evidence for a trade-off between male success in sperm competition and germline maintenance. This suggests that sex differences in the relative strengths of sexual and natural selection are causally linked to male mutation bias. The tenet advocated here, that the allocation decisions of an individual can affect plasticity of its germline and the resulting genetic quality of subsequent generations, has several interesting implications for mate choice processes.

Evolutionary dynamics of sex-biased gene expression in a young XY system: insights from the brown alga genus Fucus

Sex-biased gene expression is considered to be an underlying cause of sexually dimorphic traits. Although the nature and degree of sex-biased expression have been well documented in several animal and plant systems, far less is known about the evolution of sex-biased genes in more distant eukaryotic groups. Here, we investigate sex-biased gene expression in two brown algal dioecious species, Fucus serratus and Fucus vesiculosus, where male heterogamety (XX/XY) has recently emerged. We find that in contrast to evolutionary distant plant and animal lineages, male-biased genes do not experience high turnover rates, but instead reveal remarkable conservation of bias and expression levels between the two species, suggesting their importance in sexual differentiation. Genes with consistent male bias were enriched in functions related to gamete production, along with sperm competition and include three flagellar proteins under positive selection. We present one of the first reports, outside of the animal kingdom, showing that male-biased genes display accelerated rates of coding sequence evolution compared with female-biased or unbiased genes. Our results imply that evolutionary forces affect male and female sex-biased genes differently on structural and regulatory levels, resulting in unique properties of differentially expressed transcripts during reproductive development in Fucus algae.

Sexually antagonistic selection maintains genetic variance when sexual dimorphism evolves

Genetic variance (VG) in fitness related traits is often unexpectedly high, evoking the question how VG can be maintained in the face of selection. Sexually antagonistic (SA) selection favouring alternative alleles in the sexes is common and predicted to maintain VG, while directional selection should erode it. Both SA and sex-limited directional selection can lead to sex-specific adaptations but how each affect VG when sexual dimorphism evolves remain experimentally untested. Using replicated artificial selection on the seed beetle Callosobruchus maculatus body size we recently demonstrated an increase in size dimorphism under SA and male-limited (ML) selection by 50% and 32%, respectively. Here we test their consequences on genetic variation. We show that SA selection maintained significantly more ancestral, autosomal additive genetic variance than ML selection, while both eroded sex-linked additive variation equally. Ancestral female-specific dominance variance was completely lost under ML, while SA selection consistently sustained it. Further, both forms of selection preserved a high genetic correlation between the sexes (rm,f). These results demonstrate the potential for sexual antagonism to maintain more genetic variance while fuelling sex-specific adaptation in a short evolutionary time scale, and are in line with predicted importance of sex-specific dominance reducing sexual conflict over alternative alleles.

Omics in the Red Palm Weevil Rhynchophorus ferrugineus (Olivier) (Coleoptera: Curculionidae): A Bridge to the Pest

The red palm weevil (RPW), Rhynchophorus ferrugineus (Coleoptera: Curculionidae), is the most devastating pest of palm trees worldwide. Mitigation of the economic and biodiversity impact it causes is an international priority that could be greatly aided by a better understanding of its biology and genetics. Despite its relevance, the biology of the RPW remains poorly understood, and research on management strategies often focuses on outdated empirical methods that produce sub-optimal results. With the development of omics approaches in genetic research, new avenues for pest control are becoming increasingly feasible. For example, genetic engineering approaches become available once a species's target genes are well characterized in terms of their sequence, but also population variability, epistatic interactions, and more. In the last few years alone, there have been major advances in omics studies of the RPW. Multiple draft genomes are currently available, along with short and long-read transcriptomes, and metagenomes, which have facilitated the identification of genes of interest to the RPW scientific community. This review describes omics approaches previously applied to RPW research, highlights findings that could be impactful for pest management, and emphasizes future opportunities and challenges in this area of research.

Larger genomes show improved buffering of adult fitness against environmental stress in seed beetles

Our general understanding of the evolution of genome size (GS) is incomplete, and it has long been clear that GS does not reflect organismal complexity. Here, we assess the hypothesis that larger genomes may allow organisms to better cope with environmental variation. It is, for example, possible that genome expansion due to proliferation of transposable elements or gene duplications may affect the ability to regulate and fine-tune transcriptional profiles. We used 18 populations of the seed beetle Callosobruchus maculatus, which differ in GS by up to 4.5%, and exposed adults and juveniles to environmental stress in a series of experiments where stage-specific fitness was assayed. We found that populations with larger genomes were indeed better buffered against environmental stress for adult, but not for juvenile, fitness. The genetic correlation across populations between GS and canalization of adult fitness is consistent with a role for natural selection in the evolution of GS.

Evolutionary consequences of pesticide exposure include transgenerational plasticity and potential terminal investment transgenerational effects

Transgenerational plasticity, the influence of the environment experienced by parents on the phenotype and fitness of subsequent generations, is being increasingly recognized. Human-altered environments, such as those resulting from the increasing use of pesticides, may be major drivers of such cross-generational influences, which in turn may have profound evolutionary and ecological repercussions. Most of these consequences are, however, unknown. Whether transgenerational plasticity elicited by pesticide exposure is common, and the consequences of its potential carryover effects on fitness and population dynamics, remains to be determined. Here, we investigate whether exposure of parents to a common pesticide elicits intra-, inter-, and transgenerational responses (in F0, F1, and F2 generations) in life history (fecundity, longevity, and lifetime reproductive success), in an insect model system, the seed beetle Callosobruchus maculatus. We also assessed sex specificity of the effects. We found sex-specific and hormetic intergenerational and transgenerational effects on longevity and lifetime reproductive success, manifested both in the form of maternal and paternal effects. In addition, the transgenerational effects via mothers detected in this study are consistent with a new concept: terminal investment transgenerational effects. Such effects could underlie cross-generational responses to environmental perturbation. Our results indicate that pesticide exposure leads to unanticipated effects on population dynamics and have far-reaching ecological and evolutionary implications.

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.

Correlation scan: identifying genomic regions that affect genetic correlations applied to fertility traits

Although the genetic correlations between complex traits have been estimated for more than a century, only recently we have started to map and understand the precise localization of the genomic region(s) that underpin these correlations. Reproductive traits are often genetically correlated. Yet, we don't fully understand the complexities, synergism, or trade-offs between male and female fertility. In this study, we used reproductive traits in two cattle populations (Brahman; BB, Tropical Composite; TC) to develop a novel framework termed correlation scan (CS). This framework was used to identify local regions associated with the genetic correlations between male and female fertility traits. Animals were genotyped with bovine high-density single nucleotide polymorphisms (SNPs) chip assay. The data used consisted of ~1000 individual records measured through frequent ovarian scanning for age at first corpus luteum (AGECL) and a laboratory assay for serum levels of insulin growth hormone (IGF1 measured in bulls, IGF1b, or cows, IGF1c). The methodology developed herein used correlations of 500-SNP effects in a 100-SNPs sliding window in each chromosome to identify local genomic regions that either drive or antagonize the genetic correlations between traits. We used Fisher's Z-statistics through a permutation method to confirm which regions of the genome harboured significant correlations. About 30% of the total genomic regions were identified as driving and antagonizing genetic correlations between male and female fertility traits in the two populations. These regions confirmed the polygenic nature of the traits being studied and pointed to genes of interest. For BB, the most important chromosome in terms of local regions is often located on bovine chromosome (BTA) 14. However, the important regions are spread across few different BTA's in TC. Quantitative trait loci (QTLs) and functional enrichment analysis revealed many significant windows co-localized with known QTLs related to milk production and fertility traits, especially puberty. In general, the enriched reproductive QTLs driving the genetic correlations between male and female fertility are the same for both cattle populations, while the antagonizing regions were population specific. Moreover, most of the antagonizing regions were mapped to chromosome X. These results suggest regions of chromosome X for further investigation into the trade-offs between male and female fertility. We compared the CS with two other recently proposed methods that map local genomic correlations. Some genomic regions were significant across methods. Yet, many significant regions identified with the CS were overlooked by other methods.

Chromosome-level genome assembly of Monochamus saltuarius reveals its adaptation and interaction mechanism with pine wood nematode

Monochamus saltuarius (Coleoptera: Cerambycidae) was reported as the vector beetle of the pine wood nematode (PWN, Bursaphelenchus xylophilus) in Japan and Europe. It was first reported to transmitted the PWN to native Pinus species in 2018 in Liaoning Province, China. However, the lack of genomic resources has limited the in-depth understanding of its interspecific relationship with PWN. Here, we obtained a chromosome-level reference genome of M. saltuarius combining Illumina, Nanopore and Hi-C sequencing technologies. We assembled the scaffolds into ten chromosomes (including an X chromosome) and obtained a 682.23 Mb chromosome-level genome with a N50 of 73.69 Mb. In total, 427.67 Mb (62.69 %) repeat sequences were identified and 14, 492 protein-coding genes were predicted, of which 93.06 % were annotated. We described the mth/mthl, P450, OBP and OR gene families associated with the vector beetle's development and resistance, as well as the host selection and adaptation, which serve as a valuable resource for understanding the host adaptation in insects during evolution. This high quality reference genome of M. saltuarius also provide new avenues for researching the mechanism of this synergistic damage between vector beetles and PWN.

Evolution of resource generalism via generalized stress response confers increased reproductive thermal tolerance in a pest beetle

Generalism should be favoured evolutionarily when there is no genetic constraint or loss of fitness across alternative environments. However, evolution of generalism can require substantial evolutionary change, which can confer a general stress response to other aspects of the environment. We created generalist lineages from an ancestral, resource-specialized laboratory population of seed beetles (Callosobruchus maculatus) by rearing lines over 60 generations on a mixture of both ancestral and novel host species to test for costs associated with the evolution of generalism involving evolutionary changes in gene expression and correlated phenotypic responses during a shift to generalism. Evolved lines had higher fitness on the novel resource, with no loss of fitness on the ancestral resource, indicating that they overcame initial fitness trade-offs. This involved upregulation of major stress response (heat shock protein) genes and genes coding for metabolic enzymes, suggesting an underpinning metabolic and physiological cost. Resource generalist populations also evolved greater thermal tolerance breadth, highlighting that the evolution of resource generalism might pre-adapt species to respond favourably to other environmental stressors, following selection for generalized stress response gene upregulation. The rapid gain of novel hosts during a pest invasion might also confer greater thermal resilience to ongoing climate change.

X chromosomes show relaxed selection and complete somatic dosage compensation across Timema stick insect species

Sex chromosomes have evolved repeatedly across the tree of life. As they are present in different copy numbers in males and females, they are expected to experience different selection pressures than the autosomes, with consequences including a faster rate of evolution, increased accumulation of sexually antagonistic alleles and the evolution of dosage compensation. Whether these consequences are general or linked to idiosyncrasies of specific taxa is not clear as relatively few taxa have been studied thus far. Here, we use whole-genome sequencing to identify and characterize the evolution of the X chromosome in five species of Timema stick insects with XX:X0 sex determination. The X chromosome had a similar size (approximately 12% of the genome) and gene content across all five species, suggesting that the X chromosome originated prior to the diversification of the genus. Genes on the X showed evidence of relaxed selection (elevated dN/dS) and a slower evolutionary rate (dN + dS) than genes on the autosomes, likely due to sex-biased mutation rates. Genes on the X also showed almost complete dosage compensation in somatic tissues (heads and legs), but dosage compensation was absent in the reproductive tracts. Contrary to prediction, sex-biased genes showed little enrichment on the X, suggesting that the advantage X-linkage provides to the accumulation of sexually antagonistic alleles is weak. Overall, we found the consequences of X-linkage on gene sequences and expression to be similar across Timema species, showing the characteristics of the X chromosome are surprisingly consistent over 30 million years of evolution.

Antennal transcriptome analysis of chemosensory genes in the cowpea beetle, Callosobruchus maculatus (F.)

Olfaction, one of the most important sensory systems governing insect behavior, is a possible target for pest management. Therefore, in this study, we analyzed the antennal transcriptome of the cowpea beetle, Callosobruchus maculatus (F.) (Coleoptera: Chrysomelidae: Bruchinae), which is a major pest of stored pulses and legumes. The de novo antennal RNA-seq assembly results identified 17 odorant, 2 gustatory, and 10 ionotropic receptors, 1 sensory neuron membrane protein, and 12 odorant-binding and 7 chemosensory proteins. Moreover, differential gene expression analysis of virgin male and female antennal samples followed by qRT-PCR revealed 1 upregulated and 4 downregulated odorant receptors in males. We also performed homology searches using the coding sequences built from previously proposed amino acid sequences derived from genomic data and identified additional chemosensory-related genes.

In silico promoter analysis and functional validation identify CmZFH, the co-regulator of hypoxia-responsive genes CmScylla and CmLPCAT

Oxygen (O₂) plays an essential role in aerobic organisms including terrestrial insects. Under hypoxic stress, the cowpea bruchid (Callosobruchus maculatus) ceases feeding and growth. However, larvae, particularly 4th instar larvae exhibit very high tolerance to hypoxia and can recover normal growth once brought to normoxia. To better understand the molecular mechanism that enables insects to cope with low O₂ stress, we performed RNA-seq to distinguish hypoxia-responsive genes in midguts and subsequently identified potential common cis-elements in promoters of hypoxia-induced and -repressed genes, respectively. Selected elements were subjected to gel-shift and transient transfection assays to confirm their cis-regulatory function. Of these putative common cis-elements, AREB6 appeared to regulate the expression of CmLPCAT and CmScylla, two hypoxia-induced genes. CmZFH, the putative AREB6-binding protein, was hypoxia-inducible. Transient expression of CmZFH in Drosophila S2 cells activated CmLPCAT and CmScylla, and their induction was likely through interaction of CmZFH with AREB6. Binding to AREB6 was further confirmed by bacterially expressed CmZFH recombinant protein. Deletion analyses indicated that the N-terminal zinc-finger cluster of CmZFH was the key AREB6-binding domain. Through in silico and experimental exploration, we discovered novel transcriptional regulatory components associated with gene expression dynamics under hypoxia that facilitated insect survival.

Rapid evolution of sexual size dimorphism facilitated by Y-linked genetic variance

Sexual dimorphism is ubiquitous in nature but its evolution is puzzling given that the mostly shared genome constrains independent evolution in the sexes. Sex differences should result from asymmetries between the sexes in selection or genetic variation but studies investigating both simultaneously are lacking. Here, we combine a quantitative genetic analysis of body size variation, partitioned into autosomal and sex chromosome contributions and ten generations of experimental evolution to dissect the evolution of sexual body size dimorphism in seed beetles (Callosobruchus maculatus) subjected to sexually antagonistic or sex-limited selection. Female additive genetic variance (V_A) was primarily linked to autosomes, exhibiting a strong intersexual genetic correlation with males ([Formula: see text] = 0.926), while X- and Y-linked genes further contributed to the male V_A and X-linked genes contributed to female dominance variance. Consistent with these estimates, sexual body size dimorphism did not evolve in response to female-limited selection but evolved by 30-50% under male-limited and sexually antagonistic selection. Remarkably, Y-linked variance alone could change dimorphism by 30%, despite the C. maculatus Y chromosome being small and heterochromatic. Our results demonstrate how the potential for sexual dimorphism to evolve depends on both its underlying genetic basis and the nature of sex-specific selection.

Sex-biased gene expression and recent sex chromosome turnover

Cichlids are well known for their propensity to radiate generating arrays of morphologically and ecologically diverse species in short evolutionary time. Following this rapid evolutionary pace, cichlids show high rates of sex chromosome turnover. We here studied the evolution of sex-biased gene (SBG) expression in 14 recently diverged taxa of the Lake Tanganyika Tropheini cichlids, which show different XY sex chromosomes. Across species, sex chromosome sequence divergence predates divergence in expression between the sexes. Only one sex chromosome, the oldest, showed signs of demasculinization in gene expression and potentially contribution to the resolution of sexual conflict. SBGs in general showed high rates of turnovers and evolved mostly under drift. Sexual selection did not shape the rapid evolutionary changes of SBGs. Male-biased genes evolved faster than female-biased genes, which seem to be under more phylogenetic constraint. We found a relationship between the degree of sex bias and sequence evolution driven by sequence differences among the sexes. Consistent with other species, strong sex bias towards sex-limited expression contributes to resolving sexual conflict in cichlids. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Selection in males purges the mutation load on female fitness

Theory predicts that the ability of selection and recombination to purge mutation load is enhanced if selection against deleterious genetic variants operates more strongly in males than females. However, direct empirical support for this tenet is limited, in part because traditional quantitative genetic approaches allow dominance and intermediate-frequency polymorphisms to obscure the effects of the many rare and partially recessive deleterious alleles that make up the main part of a population's mutation load. Here, we exposed the partially recessive genetic load of a population of Callosobruchus maculatus seed beetles via successive generations of inbreeding, and quantified its effects by measuring heterosis-the increase in fitness experienced when masking the effects of deleterious alleles by heterozygosity-in a fully factorial sex-specific diallel cross among 16 inbred strains. Competitive lifetime reproductive success (i.e., fitness) was measured in male and female outcrossed F1s as well as inbred parental "selfs," and we estimated the 4 × 4 male-female inbred-outbred genetic covariance matrix for fitness using Bayesian Markov chain Monte Carlo simulations of a custom-made general linear mixed effects model. We found that heterosis estimated independently in males and females was highly genetically correlated among strains, and that heterosis was strongly negatively genetically correlated to outbred male, but not female, fitness. This suggests that genetic variation for fitness in males, but not in females, reflects the amount of (partially) recessive deleterious alleles segregating at mutation-selection balance in this population. The population's mutation load therefore has greater potential to be purged via selection in males. These findings contribute to our understanding of the prevalence of sexual reproduction in nature and the maintenance of genetic variation in fitness-related traits.

Metapopulation structure modulates sexual antagonism

Despite the far-reaching evolutionary implications of sexual conflict, the effects of metapopulation structure, when populations are subdivided into several demes connected to some degree by migration, on sexual conflict dynamics are unknown. Here, we used experimental evolution in an insect model system, the seed beetle Callosobruchus maculatus, to assess the independent and interacting effects of selection histories associated with mating system (monogamy vs. polygamy) and population subdivision on sexual conflict evolution. We confirm traditional predictions from sexual conflict theory by revealing increased resistance to male harm in females from populations with a history of intense sexual selection (polygamous populations) compared to females from populations with a history of relaxed sexual selection (monogamous populations). However, selection arising from metapopulation structure reversed the classic pattern of sexually antagonistic coevolution and led to reduced resistance in females from polygamous populations. These results underscore that population spatial structure moderates sexual selection and sexual conflict, and more broadly, that the evolution of sexual conflict is contingent on ecological context. The findings also have implications for population dynamics, conservation biology, and biological control.

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.

Genome Assembly of the Ragweed Leaf Beetle: A Step Forward to Better Predict Rapid Evolution of a Weed Biocontrol Agent to Environmental Novelties

Rapid evolution of weed biological control agents (BCAs) to new biotic and abiotic conditions is poorly understood and so far only little considered both in pre-release and post-release studies, despite potential major negative or positive implications for risks of nontargeted attacks or for colonizing yet unsuitable habitats, respectively. Provision of genetic resources, such as assembled and annotated genomes, is essential to assess potential adaptive processes by identifying underlying genetic mechanisms. Here, we provide the first sequenced genome of a phytophagous insect used as a BCA, that is, the leaf beetle Ophraella communa, a promising BCA of common ragweed, recently and accidentally introduced into Europe. A total 33.98 Gb of raw DNA sequences, representing ∼43-fold coverage, were obtained using the PacBio SMRT-Cell sequencing approach. Among the five different assemblers tested, the SMARTdenovo assembly displaying the best scores was then corrected with Illumina short reads. A final genome of 774 Mb containing 7,003 scaffolds was obtained. The reliability of the final assembly was then assessed by benchmarking universal single-copy orthologous genes (>96.0% of the 1,658 expected insect genes) and by remapping tests of Illumina short reads (average of 98.6 ± 0.7% without filtering). The number of protein-coding genes of 75,642, representing 82% of the published antennal transcriptome, and the phylogenetic analyses based on 825 orthologous genes placing O. communa in the monophyletic group of Chrysomelidae, confirm the relevance of our genome assembly. Overall, the genome provides a valuable resource for studying potential risks and benefits of this BCA facing environmental novelties.

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.

Genome Assembly of the Cold-Tolerant Leaf Beetle Gonioctena quinquepunctata, an Important Resource for Studying Its Evolution and Reproductive Barriers between Species

Coleoptera is the most species-rich insect order, yet is currently underrepresented in genomic databases. An assembly was generated for ca. 1.7 Gb genome of the leaf beetle Gonioctena quinquepunctata by first assembling long-sequence reads (Oxford Nanopore; ± 27-fold coverage) and subsequently polishing the resulting assembly with short sequence reads (Illumina; ± 85-fold coverage). The unusually large size (most Coleoptera species are associated with a reported size below 1 Gb) was at least partially attributed to the presence of a large fraction of repeated elements (73.8%). The final assembly was characterized by an N50 length of 432 kb and a BUSCO score of 95.5%. The heterozygosity rate was ± 0.6%. Automated genome annotation informed by RNA-Seq resulted in 40,568 predicted proteins, which is much larger than the typical range 17,000–23,000 predicted for other Coleoptera. However, no evidence of a genome duplication was detected. This new reference genome will contribute to our understanding of genetic variation in the Coleoptera. Among others, it will also allow exploring reproductive barriers between species, investigating introgression in the nuclear genome, and identifying genes involved in resistance to extreme climate conditions.

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.

Symbiont Genomic Features and Localization in the Bean Beetle Callosobruchus maculatus

A pervasive pest of stored leguminous products, the bean beetle Callosobruchus maculatus (Coleoptera: Chrysomelidae) associates with a simple bacterial community during adulthood. Despite its economic importance, little is known about the compositional stability, heritability, localization, and metabolic potential of the bacterial symbionts of C. maculatus. In this study, we applied community profiling using 16S rRNA gene sequencing to reveal a highly conserved bacterial assembly shared between larvae and adults. Dominated by Firmicutes and Proteobacteria, this community is localized extracellularly along the epithelial lining of the bean beetle's digestive tract. Our analysis revealed that only one species, Staphylococcus gallinarum (phylum Firmicutes), is shared across all developmental stages. Isolation and whole-genome sequencing of S. gallinarum from the beetle gut yielded a circular chromosome (2.8 Mb) and one plasmid (45 kb). The strain encodes complete biosynthetic pathways for the production of B vitamins and amino acids, including tyrosine, which is increasingly recognized as an important symbiont-supplemented precursor for cuticle biosynthesis in beetles. A carbohydrate-active enzyme search revealed that the genome codes for a number of digestive enzymes, reflecting the nutritional ecology of C. maculatus. The ontogenic conservation of the gut microbiota in the bean beetle, featuring a "core" community composed of S. gallinarum, may be indicative of an adaptive role for the host. In clarifying symbiont localization and metabolic potential, we further our understanding and study of a costly pest of stored products. IMPORTANCE From supplementing essential nutrients to detoxifying plant secondary metabolites and insecticides, bacterial symbionts are a key source of adaptations for herbivorous insect pests. Despite the pervasiveness and geographical range of the bean beetle Callosobruchus maculatus, the role of microbial symbioses in its natural history remains understudied. Here, we demonstrate that the bean beetle harbors a simple gut bacterial community that is stable throughout development. This community localizes along the insect's digestive tract and is largely dominated by Staphylococcus gallinarum. In elucidating symbiont metabolic potential, we highlight its possible adaptive significance for a widespread agricultural pest.

Impact of male trait exaggeration on sex-biased gene expression and genome architecture in a water strider

BACKGROUND

Exaggerated secondary sexual traits are widespread in nature and often evolve under strong directional sexual selection. Although heavily studied from both theoretical and empirical viewpoints, we have little understanding of how sexual selection influences sex-biased gene regulation during the development of exaggerated secondary sexual phenotypes, and how these changes are reflected in genomic architecture. This is primarily due to the limited availability of representative genomes and associated tissue and sex transcriptomes to study the development of these traits. Here we present the genome and developmental transcriptomes, focused on the legs, of the water strider Microvelia longipes, a species where males exhibit strikingly long third legs compared to females, which they use as weapons.

RESULTS

We generated a high-quality genome assembly with 90% of the sequence captured in 13 scaffolds. The most exaggerated legs in males were particularly enriched in both sex-biased and leg-biased genes, indicating a specific signature of gene expression in association with trait exaggeration. We also found that male-biased genes showed patterns of fast evolution compared to non-biased and female-biased genes, indicative of directional or relaxed purifying selection. By contrast to male-biased genes, female-biased genes that are expressed in the third legs, but not the other legs, are over-represented in the X chromosome compared to the autosomes. An enrichment analysis for sex-biased genes along the chromosomes revealed also that they arrange in large genomic regions or in small clusters of two to four consecutive genes. The number and expression of these enriched regions were often associated with the exaggerated legs of males, suggesting a pattern of common regulation through genomic proximity in association with trait exaggeration.

CONCLUSION

Our findings indicate how directional sexual selection may drive sex-biased gene expression and genome architecture along the path to trait exaggeration and sexual dimorphism.

The draft genome of the specialist flea beetle Altica viridicyanea (Coleoptera: Chrysomelidae)

BACKGROUND

Altica (Coleoptera: Chrysomelidae) is a highly diverse and taxonomically challenging flea beetle genus that has been used to address questions related to host plant specialization, reproductive isolation, and ecological speciation. To further evolutionary studies in this interesting group, here we present a draft genome of a representative specialist, Altica viridicyanea, the first Alticinae genome reported thus far.

RESULTS

The genome is 864.8 Mb and consists of 4490 scaffolds with a N50 size of 557 kb, which covered 98.6% complete and 0.4% partial insect Benchmarking Universal Single-Copy Orthologs. Repetitive sequences accounted for 62.9% of the assembly, and a total of 17,730 protein-coding gene models and 2462 non-coding RNA models were predicted. To provide insight into host plant specialization of this monophagous species, we examined the key gene families involved in chemosensation, detoxification of plant secondary chemistry, and plant cell wall-degradation.

CONCLUSIONS

The genome assembled in this work provides an important resource for further studies on host plant adaptation and functionally affiliated genes. Moreover, this work also opens the way for comparative genomics studies among closely related Altica species, which may provide insight into the molecular evolutionary processes that occur during ecological speciation.

Advancing Undergraduate Laboratory Education Using Non-Model Insect Species

Over the past decade, laboratory courses have made a fundamental shift to inquiry-based modules and authentic research experiences. In many cases, these research experiences emphasize addressing novel research questions. Insects are ideal for inquiry-based undergraduate laboratory courses because research on insects is not limited by regulatory, economic, and logistical constraints to the same degree as research on vertebrates. While novel research questions could be pursued with model insect species (e.g., Drosophila, Tribolium), the opportunities presented by non-model insects are much greater, as less is known about non-model species. We review the literature on the use of non-model insect species in laboratory education to provide a resource for faculty interested in developing new authentic inquiry-based laboratory modules using insects. Broader use of insects in undergraduate laboratory education will support the pedagogical goals of increased inquiry and resesarch experiences while at the same time fostering increased interest and research in entomology.

A mobile sex-determining region, male-specific haplotypes and rearing environment influence age at maturity in Chinook salmon

Variation in age at maturity is an important contributor to life history and demographic variation within and among species. The optimal age at maturity can vary by sex, and the ability of each sex to evolve towards its fitness optimum depends on the genetic architecture of maturation. Using GWAS of RAD sequencing data, we show that age at maturity in Chinook salmon exhibits sex-specific genetic architecture, with age at maturity in males influenced by large (up to 20 Mb) male-specific haplotypes. These regions showed no such effect in females. We also provide evidence for translocation of the sex-determining gene between two different chromosomes. This has important implications for sexually antagonistic selection, particularly that sex linkage of adaptive genes may differ within and among populations based on chromosomal location of the sex-determining gene. Our findings will facilitate research into the genetic causes of shifting demography in Chinook salmon as well as a better understanding of sex determination in this species and Pacific salmon in general.