Search

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.

Evolutionary genetics has long struggled with understanding how functional genes under selection remain polymorphic in natural populations. Taking as a starting point that natural selection is ultimately a manifestation of ecological processes, we spotlight an underemphasized and potentially ubiquitous ecological effect that may have fundamental effects on the maintenance of genetic variation. Negative frequency dependency is a well-established emergent property of density dependence in ecology, because the relative profitability of different modes of exploiting or utilizing limiting resources tends to be inversely proportional to their frequency in a population. We suggest that this may often generate negative frequency-dependent selection (NFDS) on major effect loci that affect rate-dependent physiological processes, such as metabolic rate, that are phenotypically manifested as polymorphism in pace-of-life syndromes. When such a locus under NFDS shows stable intermediate frequency polymorphism, this should generate epistatic selection potentially involving large numbers of loci with more minor effects on life-history (LH) traits. When alternative alleles at such loci show sign epistasis with a major effect locus, this associative NFDS will promote the maintenance of polygenic variation in LH genes. We provide examples of the kind of major effect loci that could be involved and suggest empirical avenues that may better inform us on the importance and reach of this process.

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.

The patterns of reproductive timing and senescence vary within and across species owing to differences in reproductive strategies, but our understanding of the molecular underpinnings of such variation is incomplete. This is perhaps particularly true for sex differences. We investigated the evolution of sex-specific gene expression associated with life history divergence in replicated populations of the seed beetle Acanthoscelides obtectus, experimentally evolving under (E)arly or (L)ate life reproduction for >200 generations which has resulted in strongly divergent life histories. We detected 1,646 genes that were differentially expressed in E and L lines, consistent with a highly polygenic basis of life history evolution. Only 30% of differentially expressed genes were similarly affected in males and females. The evolution of long life was associated with significantly reduced sex differences in expression, especially in non-reproductive tissues. The expression differences were overall more pronounced in females, in accordance with their greater phenotypic divergence in lifespan. Functional enrichment analysis revealed differences between E and L beetles in gene categories previously implicated in aging, such as mitochondrial function and defense response. The results show that divergent life history evolution can be associated with profound changes in gene expression that alter the transcriptome in a sex-specific way, highlighting the importance of understanding the mechanisms of aging in each sex.

Oxygen (O2) 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 O2 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.

Male-female coevolution has taken different paths among closely related species, but our understanding of the factors that govern its direction is limited. While it is clear that ecological factors, life history, and the economics of reproduction are connected, the divergent links are often obscure. We propose that a complete understanding requires the conceptual integration of metabolic phenotypes. Metabolic rate, a nexus of life history evolution, is constrained by ecological factors and may exert important direct and indirect effects on the evolution of sexual dimorphism. We performed standardized experiments in 12 seed beetle species to gain a rich set of sex-specific measures of metabolic phenotypes, life history traits, and the economics of mating and analyzed our multivariate data using phylogenetic comparative methods. Resting metabolic rate (RMR) showed extensive evolution and evolved more rapidly in males than in females. The evolution of RMR was tightly coupled with a suite of life history traits, describing a pace-of-life syndrome (POLS), with indirect effects on the economics of mating. As predicted, high resource competition was associated with a low RMR and a slow POLS. The cost of mating showed sexually antagonistic coevolution, a hallmark of sexual conflict. The sex-specific costs and benefits of mating were predictably related to ecology, primarily through the evolution of male ejaculate size. Overall, our results support the tenet that resource competition affects metabolic processes that, in turn, have predictable effects on both life history evolution and reproduction, such that ecology shows both direct and indirect effects on male-female coevolution.

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 exceeding

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 F(1)s as well as inbred parental "selfs," and we estimated the 4 x 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.

Our understanding of coevolution between male genitalia and female traits remains incomplete. This is perhaps especially true for genital traits that cause internal injuries in females, such as the spiny genitalia of seed beetles where males with relatively long spines enjoy a high relative fertilization success. We report on a new set of experiments, based on extant selection lines, aimed at assessing the effects of long male spines on females in Callosobruchus maculatus . We first draw on an earlier study using microscale laser surgery, and demonstrate that genital spines have a direct negative (sexually antagonistic) effect on female fecundity. We then ask whether artificial selection for long versus short spines resulted in direct or indirect effects on female lifetime offspring production. Reference females mating with males from long-spine lines had higher offspring production, presumably due to an elevated allocation in males to those ejaculate components that are beneficial to females. Remarkably, selection for long male genital spines also resulted in an evolutionary increase in female offspring production as a correlated response. Our findings thus suggest that female traits that affect their response to male spines are both under direct selection to minimize harm but are also under indirect selection (a good genes effect), consistent with the evolution of mating and fertilization biases being affected by several simultaneous processes.

Statistics matter greatly in biology, whether we like it or not. As a discipline with an empirical inclination, we are faced with data every day and we rely on inferential statistical models to make sense of it and to provide us with novel insights. Much of the time, the growing level of complexity and sophistication of the models we put to use in ecology and evolution have led to more appropriate analyses of our data. However, this is not always the case. Here, I draw attention to a classic flaw in inferential statistics that has resurfaced in a new flavor as a result of increased reliance on complex linear mixed models -the multifaceted and disturbingly persistent problem of pseudoreplication.

Background: Understanding the forces that maintain diversity across a range of scales is at the very heart of biology. Frequency-dependent processes are generally recognized as the most central process for the maintenance of ecological diversity. The same is, however, not generally true for genetic diversity. Negative frequency dependent selection, where rare genotypes have an advantage, is often regarded as a relatively weak force in maintaining genetic variation in life history traits because recombination disassociates alleles across many genes. Yet, many regions of the genome show low rates of recombination and genetic variation in such regions (i.e., supergenes) may in theory be upheld by frequency dependent selection. Results: We studied what is essentially a ubiquitous life history supergene (i.e., mitochondrial DNA) in the fruit fly Drosophila subobscura, showing sympatric polymorphism with two main mtDNA genotypes co-occurring in populations world-wide. Using an experimental evolution approach involving manipulations of genotype starting frequencies, we show that negative frequency dependent selection indeed acts to maintain genetic variation in this region. Moreover, the strength of selection was affected by food resource conditions. Conclusions: Our work provides novel experimental support for the view that balancing selection through negative frequency dependency acts to maintain genetic variation in life history genes. We suggest that the emergence of negative frequency dependent selection on mtDNA is symptomatic of the fundamental link between ecological processes related to resource use and the maintenance of genetic variation.

Mitochondrial DNA (mtDNA) consists of few but vital maternally inherited genes that interact closely with nuclear genes to produce cellular energy. How important mtDNA polymorphism is for adaptation is still unclear. The assumption in population genetic studies is often that segregating mtDNA variation is selectively neutral. This contrasts with empirical observations of mtDNA haplotypes affecting fitness-related traits and thermal sensitivity, and latitudinal clines in mtDNA haplotype frequencies. Here, we experimentally test whether ambient temperature affects selection on mtDNA variation, and whether such thermal effects are influenced by intergenomic epistasis due to interactions between mitochondrial and nuclear genes, using replicated experimental evolution inCallosobruchus maculatusseed beetle populations seeded with a mixture of different mtDNA haplotypes. We also test for sex-specific consequences of mtDNA evolution on reproductive success, given that mtDNA mutations can have sexually antagonistic fitness effects. Our results demonstrate natural selection on mtDNA haplotypes, with some support for thermal environment influencing mtDNA evolution through mitonuclear epistasis. The changes in male and female reproductive fitness were both aligned with changes in mtDNA haplotype frequencies, suggesting that natural selection on mtDNA is sexually concordant in stressful thermal environments. We discuss the implications of our findings for the evolution of mtDNA.