Intergenomic interactions between mitochondrial and Y-linked genes shape male mating patterns and fertility in Drosophila melanogaster

Support for Y-compensation of mother's curse affecting lifespan in Drosophila melanogaster

Mother's curse refers to male-biased deleterious mutations that may accumulate on mitochondria due to its strict maternal inheritance. If these mutations persist, males should ideally compensate through mutations on Y-chromosomes given its strict paternal inheritance. Previous work addressed this hypothesis by comparing coevolved and non-coevolved Y-mitochondria pairs placed alongside completely foreign autosomal backgrounds, expecting males with coevolved pairs to exhibit greater fitness due to Y-compensation. To date, no evidence for Y-compensation has been found. That experimental design assumes Y-chromosomes compensate via direct interaction with mitochondria and/or coevolved autosomes are unimportant in its function or elucidation. If Y-chromosomes instead compensate by modifying autosomal targets (or its elucidation requires coevolved autosomes), then this design could fail to detect Y-compensation. Here we address if Y-chromosomes ameliorate mitochondrial mutations affecting male lifespan in Drosophila melanogaster. Using three disparate populations we compared lifespan among males with coevolved and non-coevolved Y-mitochondria pairs placed alongside autosomal backgrounds coevolved with mitochondria. We found coevolved pairs exhibited lower mortality risk relative to non-coevolved pairs. In contrast, no such pattern was observed when coevolved and non-coevolved pairs were placed alongside non-coevolved autosomes, as with previous studies. These data are consistent with Y-compensation and highlight the importance of autosomes in this capacity. However, we cannot fully exclude the possibility that Y-autosomal coevolution independent of mitochondrial mutations contributed to our results. Regardless, modern practices in medicine, conservation, and agriculture that introduce foreign Y-chromosomes into non-coevolved backgrounds should be used with caution, as they may disrupt Y-autosome coadaptation and/or inadvertently unbridle mother's curse.

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.

Evolutionary genetics of the mitochondrial genome: insights from Drosophila

Mitochondria are key to energy conversion in virtually all eukaryotes. Intriguingly, despite billions of years of evolution inside the eukaryote, mitochondria have retained their own small set of genes involved in the regulation of oxidative phosphorylation (OXPHOS) and protein translation. Although there was a long-standing assumption that the genetic variation found within the mitochondria would be selectively neutral, research over the past 3 decades has challenged this assumption. This research has provided novel insight into the genetic and evolutionary forces that shape mitochondrial evolution and broader implications for evolutionary ecological processes. Many of the seminal studies in this field, from the inception of the research field to current studies, have been conducted using Drosophila flies, thus establishing the species as a model system for studies in mitochondrial evolutionary biology. In this review, we comprehensively review these studies, from those focusing on genetic processes shaping evolution within the mitochondrial genome, to those examining the evolutionary implications of interactions between genes spanning mitochondrial and nuclear genomes, and to those investigating the dynamics of mitochondrial heteroplasmy. We synthesize the contribution of these studies to shaping our understanding of the evolutionary and ecological implications of mitochondrial genetic variation.

Y chromosome-linked variation affects locomotor activity in male Drosophila melanogaster and is robust to differences in thermal environment

Although containing genes important for sex determination, genetic variation within the Y chromosome was traditionally predicted to contribute little to the expression of sexually dimorphic traits. This prediction was shaped by the assumption that the chromosome harbours few protein-coding genes, and that capacity for Y-linked variation to shape adaptation would be hindered by the chromosome's lack of recombination and holandric inheritance. Consequently, most studies exploring the genotypic contributions to sexually dimorphic traits have focused on the autosomes and X chromosome. Yet, several studies have now demonstrated that the Y chromosome harbours variation affecting male fitness, moderating the expression of hundreds of genes across the nuclear genome. Furthermore, emerging results have shown that expression of this Y-linked variation may be sensitive to environmental heterogeneity, leading to the prediction that Y-mediated gene-by-environment interactions will shape the expression of sexually dimorphic phenotypes. We tested this prediction, investigating whether genetic variation across six distinct Y chromosome haplotypes affects the expression of locomotor activity, at each of two temperatures (20 and 28 °C) in male fruit flies (Drosophila melanogaster). Locomotor activity is a sexually dimorphic trait in this species, previously demonstrated to be under intralocus sexual conflict. We demonstrate Y haplotype effects on male locomotor activity, but the rank order and magnitude of these effects were unaltered by differences in temperature. Our study contributes to a growing number of studies demonstrating Y-linked effects moderating expression of traits evolving under sexually antagonistic selection, suggesting a role for the Y chromosome in shaping outcomes of sexual conflict.

Migration restores hybrid incompatibility driven by mitochondrial-nuclear sexual conflict

In the mitochondrial genome, sexual asymmetry in transmission allows the accumulation of male-harming mutations since selection acts only on the effect of the mutation in females. Called the 'Mother's Curse', this phenomenon induces a selective pressure for nuclear variants that compensate for this reduction in male fitness. Previous work has demonstrated the existence of these interactions and their potential to act as Dobzhansky-Muller incompatibilities, contributing to reproductive isolation between populations. However, it is not clear how readily they would give rise to and sustain hybrid incompatibilities. Here, we use computer simulations in SLiM 3 to investigate the consequences of sexually antagonistic mitochondrial-nuclear interactions in a subdivided population. We consider distinct migration schemes and vary the chromosomal location, and consequently the transmission pattern, of nuclear restorers. Disrupting these co-evolved interactions results in less-fit males, skewing the sex ratio toward females. Restoration of male fitness depends on both the chromosomal location of nuclear restorer loci and the migration scheme. Our results show that these interactions may act as Dobzhansky-Muller incompatibilities, but their strength is not enough to drive population isolation. Overall, this model shows the varied ways in which populations can respond to migration's disruption of co-evolved mitochondrial-nuclear interactions.

The decreased expression of GW2 homologous genes contributed to the increased grain width and thousand‑grain weight in wheat-Dasypyrum villosum 6VS·6DL translocation lines

This study demonstrated that the aberrant transcription of DvGW2 contributed to the increased grain width and thousand-grain weight in wheat-Dasypyrum villosum T6VS·6DL translocation lines. Due to the high immunity to powdery mildew, Dasypyrum villosum 6VS has been one of the most successful applications of the wild relatives in modern wheat breeding. Along with the desired traits, side-effects could be brought when large alien chromosome fragments are introduced into wheat, but little is known about effects of 6VS on agronomic traits. Here, we found that T6VS·6DL translocation had significantly positive effects on grain weight, plant heightand spike length, and small negative effects on total spikelet number and spikelet compactness using recipient and wheat-D. villosum T6VS·6DL allohexaploid wheats, Wan7107 and Pm97033. Further analysis showed that the 6VS segment might exert direct genetic effect on grain width, then driving the increase of thousand-grain weight. Furthermore, comparative transcriptome analysis identified 2549 and 1282 differentially expressed genes (DEGs) and 2220 and 1496 specifically expressed genes (SEGs) at 6 days after pollination (DAP) grains and 15 DAP endosperms, respectively. Enrichment analysis indicated that the process of cell proliferation category was over-represented in the DEGs. Notably, two homologous genes, TaGW2-D1 and DvGW2, were identified as putative candidate genes associated with grain weight and yield. The expression analysis showed that DvGW2 had an aberrant expression in Pm97033, resulting in significantly lower total expression level of GW2 than Wan7107, which drives the increase of grain weight and width in Pm97033. Collectively, our data indicated that the compromised expression of DvGW2 is critical for increased grain width and weight in T6VS·6DL translocation lines.

Adaptive co-evolution of mitochondria and the Y-chromosome: A resolution to conflict between evolutionary opponents

In most species with motile sperm, male fertility depends upon genes located on the Y-chromosome and in the mitochondrial genome. Coordinated adaptive evolution for the function of male fertility between genes on the Y and the mitochondrion is hampered by their uniparental inheritance in opposing sexes: The Y-chromosome is inherited uniparentally, father to son, and the mitochondrion is inherited maternally, mother to offspring. Preserving male fertility is problematic, because maternal inheritance permits mitochondrial mutations advantageous to females, but deleterious to male fertility, to accumulate in a population. Although uniparental inheritance with sex-restricted adaptation also affects genes on the Y-chromosome, females lack a Y-chromosome and escape the potential maladaptive consequences of male-limited selection. Evolutionary models have shown that mitochondrial mutations deleterious to male fertility can be countered by compensatory evolution of Y-linked mutations that restore it. However, direct adaptive coevolution of Y- and mitochondrial gene combinations has not yet been mathematically characterized. We use population genetic models to show that adaptive coevolution of Y and mitochondrial genes are possible when Y-mt gene combinations have positive effects on male fertility and populations are inbred.

Mitochondrial-Y chromosome epistasis in Drosophila melanogaster

The coordination between mitochondrial and nuclear genes is crucial to eukaryotic organisms. Predicting the nature of these epistatic interactions can be difficult because of the transmission asymmetry of the genes involved. While autosomes and X-linked genes are transmitted through both sexes, genes on the Y chromosome and in the mitochondrial genome are uniparentally transmitted through males and females, respectively. Here, we generate 36 otherwise isogenic Drosophila melanogaster strains differing only in the geographical origin of their mitochondrial genome and Y chromosome, to experimentally examine the effects of the uniparentally inherited parts of the genome, as well as their interaction, in males. We assay longevity and gene expression through RNA-sequencing. We detect an important role for both mitochondrial and Y-linked genes, as well as extensive mitochondrial-Y chromosome epistasis. In particular, genes involved in male reproduction appear to be especially sensitive to such interactions, and variation on the Y chromosome is associated with differences in longevity. Despite these interactions, we find no evidence that the mitochondrial genome and Y chromosome are co-adapted within a geographical region. Overall, our study demonstrates a key role for the uniparentally inherited parts of the genome for male biology, but also that mito-nuclear interactions are complex and not easily predicted from simple transmission asymmetries.

An experimental test of temperature-dependent selection on mitochondrial haplotypes in Callosobruchus maculatus seed beetles

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 in Callosobruchus maculatus seed 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.

Mitonuclear gene X environment effects on lifespan and health: How common, how big?

Mitochondrial genetic variation can have profound effects on fitness, and the mitotype must interact with both the nuclear genes and the environment. We used Drosophila to investigate the extent to which mitotype effects on lifespan and activity are modulated by nucleotype and environmental variation. When nucleotype is varied, mitochondrial effects on lifespan persisted but were relatively small, and still male biased. Varying food as well, mitotype had substantial effects on male climbing speed, modifiable by nucleotype but less so by diet. Finally, mitotype affected fly lifespan much more in a cage environment compared with a vial, also modifiable by nucleotype and diet. The cage may represent a stressful environment. Mitochondrial genotype may affect fitness much more in conditions of stress, which may have implications for human health.

Sexual conflict through mother's curse and father's curse

In contrast with autosomes, lineages of sex chromosomes reside for different amounts of time in males and females, and this transmission asymmetry makes them hotspots for sexual conflict. Similarly, the maternal inheritance of the mitochondrial genome (mtDNA) means that mutations that are beneficial in females can spread in a population even if they are deleterious in males, a form of sexual conflict known as Mother's Curse. While both Mother's Curse and sex chromosome induced sexual conflict have been well studied on their own, the interaction between mitochondrial genes and genes on sex chromosomes is poorly understood. Here, we use analytical models and computer simulations to perform a comprehensive examination of how transmission asymmetries of nuclear, mitochondrial, and sex chromosome-linked genes may both cause and resolve sexual conflicts. For example, the accumulation of male-biased Mother's Curse mtDNA mutations will lead to selection in males for compensatory nuclear modifier loci that alleviate the effect. We show how the Y chromosome, being strictly paternally transmitted provides a particularly safe harbor for such modifiers. This analytical framework also allows us to discover a novel kind of sexual conflict, by which Y chromosome-autosome epistasis may result in the spread of male beneficial but female deleterious mutations in a population. We christen this phenomenon Father's Curse. Extending this analytical framework to ZW sex chromosome systems, where males are the heterogametic sex, we also show how W-autosome epistasis can lead to a novel kind of nuclear Mother's Curse. Overall, this study provides a comprehensive framework to understand how genetic transmission asymmetries may both cause and resolve sexual conflicts.

Mitochondrial genetic effects on reproductive success: signatures of positive intrasexual, but negative intersexual pleiotropy

Theory predicts that maternal inheritance of mitochondria will facilitate the accumulation of mtDNA mutations that are male biased, or even sexually antagonistic, in effect. While there are many reported cases of mtDNA mutations conferring cytoplasmic male sterility in plants, historically it was assumed such mutations would not persist in the streamlined mitochondrial genomes of bilaterian metazoans. Intriguingly, recent cases of mitochondrial variants exerting male biases in effect have come to light in bilaterians. These cases aside, it remains unknown whether the mitochondrial genetic variation affecting phenotypic expression, and in particular reproductive performance, in bilaterians is routinely composed of sex-biased or sex-specific variation. If selection consistently favours mtDNA variants that augment female fitness, but at cost to males, this could shape patterns of pleiotropy and lead to negative intersexual correlations across mtDNA haplotypes. Here, we show that genetic variation across naturally occurring mitochondrial haplotypes affects components of reproductive success in both sexes, in the fruit fly Drosophila melanogaster We find that intrasexual correlations across mitochondrial haplotypes, for components of reproductive success, are generally positive, while intersexual correlations are negative. These results accord with theoretical predictions, suggesting that maternal inheritance has led to the fixation of numerous mutations of sexually antagonistic effect.

Support for the Dominance Theory in Drosophila Transcriptomes

Interactions among divergent elements of transcriptional networks from different species can lead to misexpression in hybrids through regulatory incompatibilities, some with the potential to generate sterility. While the possible contribution of faster-male evolution to this misexpression has been explored, the role of the hemizygous X chromosome (i.e., the dominance theory for transcriptomes) remains yet to be determined. Here, we study genome-wide patterns of gene expression in females and males of Drosophila yakuba, Drosophila santomea and their hybrids. We used attached-X stocks to specifically test the dominance theory, and we uncovered a significant contribution of recessive alleles on the X chromosome to hybrid misexpression. Our analyses also suggest a contribution of weakly deleterious regulatory mutations to gene expression divergence in genes with sex-biased expression, but only in the sex toward which the expression is biased (e.g., genes with female-biased expression when analyzed in females). In the opposite sex, we found stronger selective constraints on gene expression divergence. Although genes with a high degree of male-biased expression show a clear signal of faster-X evolution of gene expression, we also detected slower-X evolution in other gene classes (e.g., female-biased genes). This slower-X effect is mediated by significant decreases in cis- and trans-regulatory divergence. The distinct behavior of X-linked genes with a high degree of male-biased expression is consistent with these genes experiencing a higher incidence of positively selected regulatory mutations than their autosomal counterparts.

Sex-specific adaptation and genomic responses to Y chromosome presence in female reproductive and neural tissues

Y chromosomes typically harbour a small number of genes and an abundance of repetitive sequences. In Drosophila, the Y chromosome comprises multimegabase long segments of repetitive DNA and a handful of protein-coding genes. In mammals, the Y chromosome also harbours a disproportionally high abundance of repeats. Here, we built on a Drosophila melanogaster model in which the Y chromosome is decoupled from sexual determination. Genotypes were genetically identical for the autosomes, X chromosome, and mitochondria, but differ by the presence or dose of the Y chromosome. Addition of an extra Y chromosome had limited impact in males. However, the presence of a Y chromosome in females induced a disproportionate response in genes expressed in the ovaries as well as genes encoded by the mitochondrial genome. Furthermore, the data revealed significant consequences of Y chromosome presence in larvae neuronal tissue. This included the repression of genes implicated in reproductive behaviour, courtship, mating and synaptic function. Our findings exhibit the Y chromosome as a hotspot for sex-specific adaptation. They suggest roles for natural selection on Y-linked genetic elements exerting impact on sex-specific tissues as well as somatic tissues shared by males and females.

Coadaptation of mitochondrial and nuclear genes, and the cost of mother's curse

Strict maternal inheritance renders the mitochondrial genome susceptible to accumulating mutations that harm males, but are otherwise benign or beneficial for females. This 'mother's curse' effect can degrade male survival and fertility if unopposed by counteracting evolutionary processes. Coadaptation between nuclear and mitochondrial genomes-with nuclear genes evolving to compensate for male-harming mitochondrial substitutions-may ultimately resolve mother's curse. However, males are still expected to incur a transient fitness cost during mito-nuclear coevolution, and it remains unclear how severe such costs should be. We present a population genetic analysis of mito-nuclear coadaptation to resolve mother's curse effects, and show that the magnitude of the 'male mitochondrial load'-the negative impact of mitochondrial substitutions on male fitness components-may be large, even when genetic variation for compensatory evolution is abundant. We also find that the male load is surprisingly sensitive to population size: male fitness costs of mito-nuclear coevolution are particularly pronounced in both small and large populations, and minimized in populations of intermediate size. Our results reveal complex interactions between demography and genetic constraints during the resolution of mother's curse, suggesting potentially widespread species differences in susceptibility to mother's curse effects.

Reducing cryptic relatedness in genomic data sets via a central node exclusion algorithm

Cryptic relatedness is a confounding factor in genetic diversity and genetic association studies. Development of strategies to reduce cryptic relatedness in a sample is a crucial step for downstream genetic analyses. This study uses a node selection algorithm, based on network degrees of centrality, to evaluate its applicability and impact on evaluation of genetic diversity and population stratification. 1,036 Guzerá (Bos indicus) females were genotyped using Illumina Bovine SNP50 v2 BeadChip. Four strategies were compared. The first and second strategies consist on a iterative exclusion of most related individuals based on PLINK kinship coefficient (φij) and VanRaden's φij, respectively. The third and fourth strategies were based on a node selection algorithm. The fourth strategy, Network G matrix, preserved the larger number of individuals with a better diversity and representation from the initial sample. Determining the most probable number of populations was directly affected by the kinship metric. Network G matrix was the better strategy for reducing relatedness due to producing a larger sample, with more distant individuals, a more similar distribution when compared with the full data set in the MDS plots and keeping a better representation of the population structure. Resampling strategies using VanRaden's φij as a relationship metric was better to infer the relationships among individuals. Moreover, the resampling strategies directly impact the genomic inflation values in genomewide association studies. The use of the node selection algorithm also implies better selection of the most central individuals to be removed, providing a more representative sample.

No effect of mitochondrial genotype on reproductive plasticity following exposure to a non-infectious pathogen challenge in female or male Drosophila

Mitochondrial genetic variation shapes the expression of life-history traits associated with reproduction, development and survival, and has also been associated with the prevalence and progression of infectious bacteria and viruses in humans. The breadth of these effects on multifaceted components of health, and their link to disease susceptibility, led us to test whether variation across mitochondrial haplotypes affected reproductive success following an immune challenge in the form of a non-infectious pathogen. We test this, by challenging male and female fruit flies (Drosophila melanogaster), harbouring each of three distinct mitochondrial haplotypes in an otherwise standardized genetic background, to either a mix of heat-killed bacteria, or a procedural control, prior to measuring their subsequent reproductive performance. The effect of the pathogen challenge on reproductive success did not differ across mitochondrial haplotypes; thus there was no evidence that patterns of reproductive plasticity were modified by the mitochondrial genotype following a non-infectious pathogen exposure. We discuss the implications of our data, and suggest future research avenues based on these results.

Interspecific Y chromosome variation is sufficient to rescue hybrid male sterility and is influenced by the grandparental origin of the chromosomes

Y chromosomes display population variation within and between species. Co-evolution within populations is expected to produce adaptive interactions between Y chromosomes and the rest of the genome. One consequence is that Y chromosomes from disparate populations could disrupt harmonious interactions between co-evolved genetic elements and result in reduced male fertility, sterility or inviability. Here we address the contribution of 'heterospecific Y chromosomes' to fertility in hybrid males carrying a homozygous region of Drosophila mauritiana introgressed in the Drosophila simulans background. In order to detect Y chromosome-autosome interactions, which may go unnoticed in a single-species background of autosomes, we constructed hybrid genotypes involving three sister species: Drosophila simulans, D. mauritiana, and D. sechellia. These engineered strains varied due to: (i) species origin of the Y chromosome (D. simulans or D. sechellia); (ii) location of the introgressed D. mauritiana segment on the D. simulans third chromosome, and (iii) grandparental genomic background (three genotypes of D. simulans). We find complex interactions between the species origin of the Y chromosome, the identity of the D. mauritiana segment and the grandparental genetic background donating the chromosomes. Unexpectedly, the interaction of the Y chromosome and one segment of D. mauritiana drastically reduced fertility in the presence of Ysim, whereas the fertility is partially rescued by the Y chromosome of D. sechellia when it descends from a specific grandparental genotype. The restoration of fertility occurs in spite of an autosomal and X-linked genome that is mostly of D. simulans origin. These results illustrate the multifactorial basis of genetic interactions involving the Y chromosome. Our study supports the hypothesis that the Y chromosome can contribute significantly to the evolution of reproductive isolation and highlights the conditional manifestation of infertility in specific genotypic combinations.

Direct and indirect genetic effects of sex-specific mitonuclear epistasis on reproductive ageing

Mitochondria are involved in ageing and their function requires coordinated action of both mitochondrial and nuclear genes. Epistasis between the two genomes can influence lifespan but whether this also holds for reproductive senescence is unclear. Maternal inheritance of mitochondria predicts sex differences in the efficacy of selection on mitonuclear genotypes that should result in differences between females and males in mitochondrial genetic effects. Mitonuclear genotype of a focal individual may also indirectly affect trait expression in the mating partner. We tested these predictions in the seed beetle Callosobruchus maculatus, using introgression lines harbouring distinct mitonuclear genotypes. Our results reveal both direct and indirect sex-specific effects of mitonuclear epistasis on reproductive ageing. Females harbouring coadapted mitonuclear genotypes showed higher lifetime fecundity due to slower senescence relative to novel mitonuclear combinations. We found no evidence for mitonuclear coadaptation in males. Mitonuclear epistasis not only affected age-specific ejaculate weight, but also influenced male age-dependent indirect effects on traits expressed by their female partners (fecundity, egg size, longevity). These results demonstrate important consequences of sex-specific mitonuclear epistasis for both mating partners, consistent with a role for mitonuclear genetic constraints upon sex-specific adaptive evolution.

Complex mitonuclear interactions and metabolic costs of mating in male seed beetles

The lack of evolutionary response to selection on mitochondrial genes through males predicts the evolution of nuclear genetic influence on male-specific mitochondrial function, for example by gene duplication and evolution of sex-specific expression of paralogs involved in metabolic pathways. Intergenomic epistasis may therefore be a prevalent feature of the genetic architecture of male-specific organismal function. Here, we assess the role of mitonuclear genetic variation for male metabolic phenotypes [metabolic rate and respiratory quotient (RQ)] associated with ejaculate renewal, in the seed beetle Callosobruchus maculatus, by assaying lines with crossed combinations of distinct mitochondrial haplotypes and nuclear lineages. We found a significant increase in metabolic rate following mating relative to virgin males. Moreover, processes associated with ejaculate renewal showed variation in metabolic rate that was affected by mitonuclear interactions. Mitochondrial haplotype influenced mating-related changes in RQ, but this pattern varied over time. Mitonuclear genotype and the energy spent during ejaculate production affected the weight of the ejaculate, but the strength of this effect varied across mitochondrial haplotypes showing that the genetic architecture of male-specific reproductive function is complex. Our findings unveil hitherto underappreciated metabolic costs of mating and ejaculate renewal, and provide the first empirical demonstration of mitonuclear epistasis on male reproductive metabolic processes.