Fixation of a deleterious allele at one of two "duplicate" loci by mutation pressure and random drift

We consider a diploid population and assume two gene loci with two alleles each, A and a at one locus and B and b at the second locus. Mutation from wild-type alleles A and B to deleterious alleles a and b occurs with mutation rates va and vb, respectively. We assume that alleles are completely recessive and that only the double recessive genotype aabb shows a deleterious effect with relative fitness 1-epsilon. Then, it can be shown that if va greater than vb mutant a becomes fixed in the population by mutation pressure and a mutation-selection balance is ultimately attained with respect to the B/b locus alone. The main aim of this paper is to investigate the situation in which va = vb exactly. In this case a neutral equilibrium is attained and either locus can drift to fixation for the mutant allele. Diffusion models are developed to treat the stochastic process involved whereby the deleterious mutant eventually becomes fixed in one of the two duplicated loci by random sampling drift in finite populations. In particular, the equation for the average time until fixation of mutant a or b is derived, and this is solved numerically for some combinations of parameters 4Nev and 4Ne epsilon, where v is the mutation rate (va = vb = v) and Ne is the effective size of the population. Monte Carlo experiments have been performed (using a device termed "pseudo sampling variable") to supplement the numerical analysis.