Accumulation of Deleterious Mutations Near Sexually Antagonistic Genes

Mutation generates a steady supply of genetic variation that, while occasionally useful for adaptation, is more often deleterious for fitness. Recent research has emphasized that the fitness effects of mutations often differ between the sexes, leading to important evolutionary consequences for the maintenance of genetic variation and long-term population viability. Some forms of sex-specific selection—i.e., stronger purifying selection in males than females—can help purge a population’s load of female-harming mutations and promote population growth. Other scenarios—e.g., sexually antagonistic selection, in which mutations that harm females are beneficial for males—inflate genetic loads and potentially dampen population viability. Evolutionary processes of sexual antagonism and purifying selection are likely to impact the evolutionary dynamics of different loci within a genome, yet theory has mostly ignored the potential for interactions between such loci to jointly shape the evolutionary genetic basis of female and male fitness variation. Here, we show that sexually antagonistic selection at a locus tends to elevate the frequencies of deleterious alleles at tightly linked loci that evolve under purifying selection. Moreover, haplotypes that segregate for different sexually antagonistic alleles accumulate different types of deleterious mutations. Haplotypes that carry female-benefit sexually antagonistic alleles preferentially accumulate mutations that are primarily male harming, whereas male-benefit haplotypes accumulate mutations that are primarily female harming. The theory predicts that sexually antagonistic selection should shape the genomic organization of genetic variation that differentially impacts female and male fitness, and contribute to sexual dimorphism in the genetic basis of fitness variation.

. This scenario arises when deleterious alleles are, on average, more harmful to male fitness than to female fitnessa pattern that has been repeatedly documented in lab populations of Drosophila Agrawal 2008, 2013;Hollis et al. 2009;, though less consistently in other animal systems (e.g., Almbro and Simmons 2014;Power and Holman 2015). Sex-limited deleterious mutations, which are harmful for one sex and neutral for the other, represent an extreme class of mutations with unequal deleterious effects between the sexes. Such mutations readily accumulate within populations due to being sheltered from purifying selection within one of the two sexes (Van Dyken and Wade 2010; ; the accumulation of both female-and male-limited mutations within a population is expected to elevate the genetic loads of both sexes (see Day and Bonduriansky 2004;.
Sexually antagonistic mutations, which individually improve fitness within one sex while reducing it in the other (see Rice 1992), represent an additional source of genetic variation for fitness in animal and plant populations (Chippindale et al. 2001;Delph et al. 2004;Foerster et al. 2007;Abbott 2011). Sexually antagonistic alleles are expected to segregate at relatively high population frequencies compared to unconditionally deleterious mutations (Connallon and Clark 2012). Some sexually antagonistic alleles may be stably maintained at intermediate frequencies by balancing selection (e.g., Kidwell et al. 1977;Jordan and Connallon 2014;Tazzyman and Abbott 2015), which occurs most readily when the strength of antagonistic selection is strong and roughly equal between the sexes (Kidwell et al. 1977;Patten and Haig 2009;Fry 2010). In such cases, sexually antagonistic alleles will similarly inflate the genetic loads of both sexes (Whitlock and Agrawal 2009;. Finally, sex differences in selection appear to indirectly favor the evolution of male-biased mutation rates, which potentially elevate the total mutation rate in species with separate sexes (Ellegren 2007;Hedrick 2007;Whitlock and Agrawal 2009). Male-biased mutation rates are common among vertebrates (Ellegren 2007;Hedrick 2007), and are occasionally reported within invertebrates (Bachtrog 2008). Many species exhibit a twofold or higher mutation rate in males than females, and some taxa exceed ratios of 8-to 10-fold (Hedrick 2007). While male-biased mutation can potentially elevate the total deleterious mutation rate of a species, and thereby increase the average load of deleterious mutations per genome, this theoretical prediction currently lacks direct empirical support, and remains open to investigation (for discussion, see Whitlock and Agrawal 2009).
Here, we demonstrate that sexually antagonistic genetic variation, sexually asymmetric selection against deleterious alleles, and male-biased mutation rates interact to shape the burden of deleterious mutations that are carried by individuals of a population. Inspired by recent theoretical work on multilocus balancing selection on sexually antagonistic alleles (Patten et al. 2010(Patten et al. , 2013Úbeda et al. 2011;Haig et al. 2014;Patten 2014), we demonstrate that sexually antagonistic polymorphism invariably promotes the accumulation of deleterious alleles in genes that are physically linked to sexually antagonistic loci. The interaction between sexual antagonism and linkage arises from the unequal transmission patterns of sexually antagonistic variation, where male-beneficial alleles are preferentially transmitted from fathers to offspring, and female-beneficial alleles are preferentially maternally inherited (Day and Bonduriansky 2004;Patten 2014). This transmission bias mediates the strength of selection against deleterious mutations that cosegregate with each sexually antagonistic allele. We show that femalelimited deleterious alleles readily accumulate in tight linkage with malebenefit sexually antagonistic alleles. Likewise, male-limited deleterious mutations accumulate near antagonistic alleles with female-beneficial effects. We also show that sexually antagonistic polymorphism distorts the mutation rates of linked genes, with elevated mutation rates upon haplotypes that carry male-benefit sexually antagonistic alleles. The theory predicts an elevated deleterious mutation load at loci that are linked to sexually antagonistic genes, with haplotypes that segregate for sexually antagonistic alleles each enriched in deleterious mutations that are particularly damaging for one sex and relatively benign for the other.

THE MODEL
We track the evolution of a pair of biallelic loci with arbitrary degree of physical linkage. Locus A is evolving under sexually antagonistic selection, with a male-beneficial allele (A m ) and female-beneficial allele (A f ). We assume that A f and A m are maintained by balancing selection, and ignore effects of mutation at the locus. Locus B segregates for two alleles: B 1 is a wild-type allele with relative fitness of one, and B 2 is a deleterious mutation that is maintained at a balance between recurrent mutation and purifying selection to remove it from the population. For simplicity, we assume that the per gamete mutation rate from B 1 to B 2 is u f in females and u m in males, with no back-mutation. We also assume that the pair of loci (A and B) affect fitness multiplicatively. The loci are separated by a recombinational distance of r (0 , r , 1/2), which reflects the probability of a recombination event between the loci, per meiosis. Generations are discrete, and follow the order of: (i) birth, (ii) viability selection, (iii) gametogenesis including mutation, and (iv) random mating and reproduction among the adults.
We consider two inheritance modelshaploidy and diploidythat arise often within dioecious animal and plant species. For the haploid modelapplicable, for example, in primitive plant species (see Jesson and Garnock-Jones 2012) and algaemeiosis occurs shortly after fertilization, and individuals spend most of their life cycle within the haploid stage. Consequently, selection occurs among haploid individuals, and genetic dominance does not factor into the evolutionary dynamics at either locus. In contrast, for the diploid model, which applies in many dioecious animals and higher plants, meiosis co-occurs with the production of gametes, selection acts among diploid individuals, and dominance becomes critically important in the dynamics of both sexually antagonistic and deleterious alleles (but see Immler et al. 2012;Otto et al. 2015).
For each model, we follow the evolution of four haplotypes that segregate (within zygotes) at frequencies: , and x 4 = [A f B 2 ]. Recursion equations follow from standard population genetics models of sex-specific selection (e.g., Kidwell et al. 1977;Connallon and Clark 2010), and are presented in full within the Appendix. The sex-specific fitness values for each of the resulting four haploid genotypes and nine diploid genotypes are presented in Table 1. Since we are primarily concerned with the overall frequency of deleterious mutations within the population, and the relative frequency of mutations on A f -bearing and A m -bearing haplotypes, our analysis focuses on each of the two mutant-bearing haplotypes: the A m B 2 haplotype at frequency x 2 , and the A f B 2 haplotype at frequency x 4 . The total frequency of the deleterious allele is simply their sum: x 2 + x 4 .
We adopt two approaches in analyzing the polymorphic equilibrium for each model of inheritance. We first develop a set of relatively simple analytical approximations for the case of tight linkage (r / 0; see the Appendix), and assuming weak mutation relative to the strength of purifying selection (i.e., u m + u f ,, t m + t f ). Second, we simulate exact equilibrium conditions under arbitrary mutation, selection, and recombination between loci. The case of tight linkage serves as a useful benchmark that provides an upper limit for the effect of sexually antagonistic polymorphism on the accumulation of deleterious mutations. As we later show, the effect of sexually antagonistic polymorphism on deleterious alleles predictably diminishes as the recombination rate increases. Ultimately, standard mutation-selection balance approximations (e.g., Haldane 1937;Crow and Kimura 1970) apply as linkage becomes loose (e.g., r / 1/2).

Data availability
The authors state that all data necessary for confirming the conclusions presented in the article are represented fully within the article.

RESULTS
We subdivide our results into three main sections. We first develop a general model that characterizes the qualitative effect of sexually antagonistic polymorphism on the evolution of a tightly linked locus evolving under purifying selection. Second, we develop haploid-and diploid-specific models of sexually antagonistic polymorphism to quantify deleterious mutation accumulation near sexually antagonistic genes; this once again assumes tight linkage between loci. Third, we present simulation results that evaluate the interaction between sexually antagonistic loci and partially linked deleterious alleles (e.g., under modest recombination rates between loci: 0 , r , 1/2).

Mutation accumulation under tight linkage to a sexually antagonistic locus
Under tight linkage between locus A and locus B, and weak mutation to the deleterious allele (r, u m , u f / 0), the equilibrium frequencies of B 2bearing haplotypes are: is the sex-averaged mutation rate to the B 2 allele,q = [A m ] is the equilibrium frequency of the male-benefit sexually antagonistic allele, d = [(A m | adult males) -(A m | adult females)] is the frequency difference of the A m allele between breeding males and females, and C is a constant that accounts for ploidy level (C = 1 for haploids; C = 1/h for diploids). Note that all five terms must take positive values (C, d, u, t, q . 0).
Sexually antagonistic polymorphism influences the cosegregation patterns of linked deleterious alleles. Equations (1a, b) can be used to quantify an individual's probability of inheriting a deleterious mutation, given their allelic status at a sexually antagonistic locus. The "relative risk" of inheriting a deleterious mutation, given inheritance of a malevs. a female-benefit allele at the sexually antagonistic locus, is: which follows from the similar concept of relative risk from epidemiology (see Jewell 2004, chapter 4). RR is scaled so that RR . 1 corresponds to an enrichment of the deleterious allele on haplotypes that carry the malebeneficial allele (A m ); RR , 1 corresponds to enrichment of deleterious alleles on the female-beneficial haplotype (A f ). The final approximation, to first order in d, shows that deleterious alleles tend to preferentially segregate on A m -bearing haplotypes when their fitness effects are more severe in females than males (t f . t m ), and when mutation rates are higher in males than females (u m . u f ). Taking both factors into account, deleterious alleles in which t f /t m . u f /u m will be enriched on A m -bearing haplotypes, and alleles in which t m /t f . u m /u f will be enriched on A fbearing haplotypes. When mutation rates are male-biased, the first of the two conditions is the most permissible, and under this scenario, deleterious mutations should preferentially accumulate on haplotypes that carry male-benefit sexually antagonistic alleles. Sexually antagonistic polymorphism also alters the total frequency of deleterious alleles in the population. In the absence of sexually antagonistic selection, the equilibrium frequency of the deleterious allele iŝ q B2 ¼x 2 þx 4 Cð u= tÞ, which corresponds to the standard result under mutation-selection balance. Linkage to a sexually antagonistic locus modifies the mutation-selection equilibrium to: Equation (3) provides insight into the effects of sexually antagonistic polymorphism on the accumulation of deleterious alleles. Linkage to a sexually antagonistic locus has no effect on the overall frequency of a deleterious allele unless the strength of purifying selection differs between the sexes (i.e., t m 6 ¼ t f ). When purifying selection is n Table 1 Sex-specific fitness for each genotype (w f and w m for females and males, respectively) For the sexually antagonistic locus (locus A), there are four parameters: sex-specific selection coefficients, s f and s m , against the deleterious allele for each sex, and dominance coefficients, k f and k m (0 , s m , s f , 1; 0 , k m , k f , 1). We use three parameters for the B locus: sex-specific selection coefficients, t m and t f , for the cost of the deleterious allele (B 2 ), and a single dominance coefficient (h) with respect to B 2 , which we assume does not differ between the sexes (0 , t m , t f , h , 1).
asymmetric between the sexes, the overall frequency of the deleterious allele can potentially increase or decrease relative to standard theory of mutation-selection balance. The deleterious allele decreases in frequency when (1t m /t f )(u m /u ft m /t f ) , 0, and it increases when Although sexual antagonism can, in principle, increase or decrease the total frequency of linked deleterious alleles, the former is a much more likely outcome than the latter. A decrease in frequency is only possible for mutations that carry a higher fitness cost in the sex with the higher mutation rate, but even in such cases, conditions remain restrictive for a net decrease of deleterious allele frequency (a net decrease requires that one of the following conditions is true: These conditions will be difficult to meet when mutation rates are similar between sexes, or when deleterious mutations have strongly asymmetric effects on male and female fitness. For example, mutations with sex-limited fitness effects always increase in frequency when linked to a sexually antagonistic locus. Sexual antagonism and mutation accumulation in haploid and diploid species Sexually antagonistic allele frequency differences between adult males and females favor the accumulation of deleterious alleles within a population, as demonstrated above [i.e., equation (3)]. Such allele frequency differences between sexes also give rise to biases in the mutation load linked to malevs. female-beneficial sexually antagonistic alleles [i.e., equation (2)]. The magnitude of these effects scales positively with the degree of the sexually antagonistic allele frequency differentiation between adult males and females, represented by d = [(A m | adult males) -(A m | adult females)]. To quantify effects of sexually antagonistic polymorphism on mutation accumulation, we develop expressions for d under haploid and diploid models of sexually antagonistic selection. We then use these expressions to identify scenarios of antagonistic polymorphism that are most favorable for deleterious mutation accumulation. We start with expressions for d for the haploid case, and then consider diploids. We then use each result to quantify effects of sexually antagonistic alleles on linked deleterious mutations.
Sexually antagonistic variation in haploids: Under a haploid model of balancing selection at a sexually antagonistic locus, the equilibrium frequency of the male-beneficial allele (in zygotes) is: The allele frequency difference between adults of each sex is: where w m ¼ 1 2 ð1 2q hap Þs m , w f ¼ 1 2q hap s f , andq j represents the frequency of A m within adults of sex j ( j = {m, f }). The maximum frequency difference between adults is d hap = 1, which occurs when s m = s f = 1. The allele frequency difference between adult females and males increases as the equilibrium for each allele approaches equality (q hap / 0.5, as s m / s f ), and as the strength of selection within each sex increases. Increasing d hap leads to increasingly pronounced differences in the mutation load between haplotypes that carry different sexually antagonistic alleles (Figure 1).

Sexually antagonistic variation in diploids:
In diploids, the pattern of genetic dominance at a sexually antagonistic locus greatly impacts the parameter space for balanced polymorphism (Kidwell et al. 1977;Fry 2010), along with the magnitude of allele frequency differences between adult females and males (as shown below). Following previous models (Fry 2010;Clark 2012, 2013), we focus on two idealized versions of dominance among sexually antagonistic alleles: (1) an additive effects model, where the fitness of heterozygotes of each sex is intermediate between homozygote fitnesses (k f = k m = 1/2 in Table 1); and (2) a strong dominance reversal model, in which the sexually antagonistic allele that is beneficial for each sex has a dominant fitness effect within that sex (k m = k f = 0 in Table 1; see Barson et al. 2015 for an empirical example).
For the additive model, the equilibrium for the male-beneficial allele under balancing selection is identical to that of the haploid model (q add ¼q hap ). The allele frequency difference between adult females and males is: which reduces to d add ¼ s fq hap 1 2q hap = w f ¼ s mq hap 1 2q hap = w m when sexually antagonistic selection is not too strong (e.g., when s f , s m , 1/2). The maximum allele frequency difference between diploid sexes under the additive model is d add = (O2 -1) 0.414.
The parameter space for balancing selection increases under dominance reversal conditions within diploids (e.g., k f = k m , 1/2; Kidwell et al. 1977;Fry 2010). For the complete dominance reversal scenario (k f = k m = 0), the equilibrium frequency of the male-beneficial allele is approximately:q DR s m s f þ s m ; (6a) and the allele frequency difference between adult males and females is: d DR 2s mq ð12qÞ 2 1 2 ð12qÞ 2 s m 2s fq 2 ð1 2qÞ 1 2q 2 s f : Equation (6b) is accurate unless sexually antagonistic selection is very strong (i.e., the approximation works well for selection coefficients up to roughly s m , s f 0.5), and otherwise tends to slightly overestimate allele frequency difference between adult sexes. The maximum for d DR is roughly one third, which falls below the maxima for both the haploid and additive diploid models of antagonistic selection.
Polymorphism in haploids vs. diploids: Allele frequency differences between adult sexes at the sexually antagonistic locus are often less pronounced in diploids than haploids. For the additive fitness model [contrasting equations (4b) and (5)], allele frequency differences between adult males and females will always be at least twofold higher in haploids than diploids for each parameter set that maintains a balanced polymorphism. Under weak selection, 2d add d hap , and the discrepancy between d add and d hap further increases with the strength of antagonistic selection. Thus, sex-biased mutation loads that accumulate near sexually antagonistic alleles will tend to be greater in haploids than diploids. Adult allele frequency differences at the sexually antagonistic locus are also greater in haploids than diploids with dominance reversal (k m = k f = 0) (Figure 2), unless the selection coefficients are asymmetric between the sexes. As the strength of sexually antagonistic selection becomes increasingly asymmetric, equilibrium heterozygosity at the sexually antagonistic locus declines rapidly in the haploid and additive diploid models (i.e., the equilibriaq hap andq dip approach zero or one as s ms f deviates from zero). In contrast, under asymmetric sexually antagonistic selection, substantial heterozygosity is still likely under dominance reversal conditions. When sexually antagonistic selection is weak (s m , s f ,, 1), d DR . d hap when min(q hap , 1 -q hap ) , 0.15, and stronger selection lowers this threshold (Figure 2). A qualitatively similar pattern is observed when contrasting allele frequency differences between the sexes under the additive and dominance reversal scenarios within diploids (Figure 2). In this case, however, d DR . d add is expected under most parameter conditions, with exceptions occurring when selection coefficients at the sexually antagonistic locus are nearly identical between the sexes.
Sexually antagonistic allele frequency differentiation between adult females and males predicts the degree of mutation accumulation at genes that are linked to sexually antagonistic loci (in terms of the overall frequency of B 2 , and the relative enrichment on haplotypes that carry male-or female-beneficial sexually antagonistic alleles). For example, when the strength of sexually antagonistic selection is similar between sexes (e.g., s f = s m ), sex-bias in genetic load (i.e., skew in the relative risk of carrying a deleterious allele between haplotypes) is more pronounced in haploids than diploids. Figure 3 illustrates this outcome both due to sex-limited deleterious effects and sex-biased mutation rates, and that these effects can increase strongly with the strength of sexually antagonistic selection under haploidy or diploidy. In contrast, as sexually antagonistic selection becomes increasingly asymmetric between the sexes, balancing selection at a sexually antagonistic locus, as well as mutation accumulation at nearby, linked loci, will be most pronounced under dominance reversal conditions, as apply within diploids. Figure 4 illustrates that the effects described above diminish rapidly with increasing recombination. Both strong sexually antagonistic selection (s m = s f = 0.5) and tight linkage to the deleterious locus are required to appreciably increase the frequency of the deleterious allele (the example in the figure exhibits a mild frequency increase of B 2 even under tight linkage; right panel). In contrast, when selection on sexually antagonistic alleles is sufficiently strong, the relative risk of inheriting deleterious alleles can substantially differ between haplotypes that carry malevs. female-benefit alleles, even under incomplete linkage between loci. Nevertheless, relative risk does decay rapidly away from the sexually antagonistic locus, with enrichment of deleterious alleles most likely at small recombinational distances (e.g., r = 0.05 or less; left panel).

DISCUSSION
Opportunities for maintaining sexually antagonistic polymorphisms are theoretically predicted to increase when sexually antagonistic selection is strong, when the strength of antagonistic selection is relatively symmetrical between the sexes, and when patterns of genetic dominance for fitness are sex-specific (as in the dominance reversal condition outlined above; see Kidwell et al. 1977;Patten and Haig 2009;Fry 2010;Jordan and Charlesworth 2012;Clark 2012, 2014;Arnqvist et al. 2014). Recent extensions of this single-locus theory suggest that genetic systems should also evolve to meet conditions that favor the maintenance of sexually antagonistic polymorphisms. For example, sexually antagonistic selection generates patterns of linkage disequilibrium that effectively increase the net strength of selection across tightly linked pairs of sexually antagonistic loci (e.g., Patten et al. 2010). Sexual antagonism also favors the evolution of genetic modifiers of sex-specific dominance, which stabilize balanced polymorphisms at sexually antagonistic loci (Spencer and Priest 2016).
Our results demonstrate that the evolutionary conditions that are most conducive for maintaining sexually antagonistic polymorphisms also promote the accumulation of deleterious mutations. Sexually antagonistic variation inflates the overall frequencies of deleterious mutations at genes that are tightly linked to sexually antagonistic loci, typically to greatest effect in haploid species. Moreover, and to a greater degree, haplotypes that carry male-benefit vs. female-benefit sexually antagonistic alleles will differ in their relative loads of sex-specific mutations. Male-beneficial haplotypes become enriched for mutations that preferentially harm female fitness, whereas female-benefit haplotypes harbor predominantly male-harming deleterious alleles. This effect amplifies the fitness costs of sexually antagonistic alleles when they are inherited by the "wrong" sex. Females (or males) that inherit malebenefit (female-benefit) alleles will simultaneously inherit an increased linked load of female-harming (male-harming) deleterious mutations. This symmetrical packaging of female-and male-harming alleles on distinct haplotypes may be somewhat offset by male-biased mutation rates, which unequivocally increase the load of mutations carried on haplotypes with male-benefit sexually antagonistic alleles.
Although we framed our models in terms of separate-sexed species, our results should equally apply to obligate outcrossing hermaphrodite species given assumptions regarding the fertility of matings between different genotypes (Bodmer 1965;Jordan and Connallon 2014;Tazzyman and Abbott 2015). In this case, both sexually antagonistic and deleterious alleles impact female and male sex functions within individuals (e.g., pollen or ovule production in plants). Self-fertilization in hermaphrodite populations can potentially impact our results in three ways: (1) by decreasing the potential for sexually antagonistic polymorphism (Jordan and Connallon 2014;Tazzyman and Abbott 2015); (2) by facilitating purging of partially or completely recessive deleterious alleles (Charlesworth and Charlesworth 1987); and (3) by decreasing the effective recombination rate between sexually antagonistic and deleterious loci. The net outcome of these contrasting forces for genetic load awaits investigation.
The vast majority of functional loci in the genome evolve under purifying selection (Eyre-Walker and Keightley 2007), and consequently, our predictions should apply to any genomic region that is linked to a polymorphic sexually antagonistic locus. Although strong sexually antagonistic selection is necessary to drastically shift the equilibrium frequency of a single deleterious mutation (i.e., relative to predictions from standard theory of mutation-selection balance; e.g., Haldane 1937;Crow and Kimura 1970), weak sexual antagonism may nevertheless generate modest shifts in the deleterious allele frequencies across many linked sites. Thus, the cumulative mutation burden on sexually antagonistic haplotypes may still be appreciable. Yet strong selection may often be plausible for individual cases of sexually antagonistic polymorphism. Instances of strong sexually antagonistic polymorphism have been implied in several study systems, including delicate skinks (Forsman and Shine 1995), collared flycatchers (Husby et al. 2015), Atlantic salmon (Barson et al. 2015), and most recently, an inversion polymorphism in the yellow monkey-flower (Lee et al. 2016). Such inversions may be particularly prone to the accumulation of deleterious alleles, because recombination can be suppressed between haplotypes comprising many genes (e.g., .1200 genes; Lee et al. 2016).
Our results demonstrate that deleterious mutations can, at least in some contexts, accumulate more readily in predominantly haploid than Figure 4 Mutation accumulation under diploid selection and arbitrary linkage between loci. The left panel shows the relative risk of a deleterious mutation on haplotypes that carry a male-beneficial vs. a female-beneficial sexually antagonistic allele (as in Figure 3). The right panel displays the equilibrium frequency of a deleterious allele relative to its expected frequency in the absence of sexually antagonistic selection at the A locus (as in Figure 1). Results are based on deterministic forward simulations, using recursions from the Appendix. Both panels show the accumulation of a partially recessive (h = 0.2) deleterious mutation with female-limited effect (t m = 0; t f = 0.1), and under sexually symmetric mutation (u f = u m = 10 25 ). The sexually antagonistic locus evolves under balancing selection with additive and symmetrical fitness effects between the sexes (s m = s f = 0.1, 0.25 or 0.5; k m = k f = 1/2 within Table 1).

Figure 3
The relative risk of a deleterious mutation on haplotypes that carry a male-beneficial vs. a female-beneficial sexually antagonistic allele. Results show the relative enrichment of a deleterious allele on haplotypes that carry the sexually antagonistic allele A m . The left panel shows the accumulation of a deleterious mutation with female-limited effect (t m = 0; t f = 0.1) under sexually symmetric mutation (u f = u m = 10 25 ). The right panel shows the accumulation of a symmetrically deleterious mutation (t m = t f = 0.1) under male-biased mutation (u f = 10 25 ; u m /u f = 8). Sexually antagonistic selection is symmetric between the sexes (s m = s f ; q hap = q add = q DR = 0.5). Each curve is based on the general equation (2)  in diploid species. This result runs counter to the traditional view in which purifying selection is enhanced in haploid relative to diploid species, with haploids carrying a reduced genetic load under most conditions of dominance (Crow and Kimura 1965;Otto and Goldstein 1992). It is currently difficult to evaluate the extent of sexually antagonistic genetic variation in haploid species, along with opportunities for the accumulation of deleterious alleles near sexually antagonistic genes. Most empirical work on sexual antagonism has, to date, focused on species with predominantly diploid life cycles, though many dioecious species do have a dominant haploid phase (e.g., Jesson and Garnock-Jones 2012). Such species represent intriguing systems for studying genetic conflicts and their evolutionary genetic consequences (e.g., Immler and Otto 2015). Predictions from the current study should hopefully provide additional motivation for empirical research on sexual antagonism in species with predominantly haploid life cycles.
The enrichment of deleterious alleles near sexually antagonistic loci represents one of several contexts of mutation accumulation that should differentially affect subregions of each sexual genome. Three forms of linked selection are particularly relevant to the evolutionary processes that we have emphasized. First, genes with sex-limited transmission patternsi.e., those on the W or Y chromosome, and genes encoded by the mitochondrial genomenaturally accumulate mutations with sexbiased fitness effects. Such loci primarily respond to selection in only one of the two sexes, and can fix mutations that harm the other (Frank and Hurst 1996;Charlesworth et al. 2014).
Second, obligate outcrossing plant species typically exhibit high levels of genetic diversity at self-incompatibility (SI) loci. These SI loci promote heterozygosity at sites to which they are physically linked, and thereby shelter recessive deleterious mutations from purifying selection (see Glémin et al. 2001;Stift et al. 2013; also see Strobeck 1972 for consideration in heterostylous plants, and van Oosterhout 2009, for comparable models at MHC loci). The ultimate consequence of balancing selection at SI loci is the accumulation of mutations at nearby genesa similar outcome to the one outlined here, though the underlying mechanisms of selection and genetic transmission differ.
Finally, tight physical linkage can itself promote the accumulation of deleterious alleles by reducing the overall efficacy of natural selection relative to genetic drift (e.g., due to Hill-Robertson interference between linked loci under selection; Hill and Robertson 1966;Charlesworth 2012). Genomic variation in recombination rates gives rise to heterogeneity across the genome for the efficacy of natural selection, and genomic regions with little or no recombination are particularly prone to functional decaythrough the accumulation of deleterious mutations and a failure to retain beneficial ones. Linkage disequilibrium, arising from selection among linked loci, can also promote the evolution of increased or decreased recombination (Otto and Lenormand 2002), including sexually dimorphic recombination rates (Lenormand 2003;Lenormand and Dutheil 2006). The form of linkage disequilibrium predicted hereof sexually antagonistic alleles cosegregating with mutations that carry sex-biased fitness costsprovides an intriguing context for the evolution of sexually dimorphic recombination rates (as observed in many animal and plant taxa; Lenormand and Dutheil 2006), for which sex differences in selection and haploidy have both been shown to play important roles (Lenormand 2003).

ACKNOWLEDGMENTS
We thank Deborah and Brian Charlesworth for providing key references to the related theory of mutation accumulation under linkage to a self-incompatibility locus. Stephen Wright and two anonymous reviewers made useful suggestions that improved the manuscript. This research was funded by the School of Biological Sciences at Monash University, and a grant from the Australian Research Council.