Mean Time to Resolution of Gene Duplication

The mean time to resolution of gene duplication (T(r)) is studied in this paper under the double null recessive (DNR) and haplo-insufficient (HI) models within the same analytical and simulation framework. We show that when population size is not too small (more precisely Nmu > 0.1), T(r) for unlinked duplication is usually larger than that for linked and T(r) for unlinked duplication under the HI model might be greatly prolonged, which were consistent with previous observations. Furthermore, by analytical approach we here indicate the primary underlying mechanism is that the frequency of the original (or wild-type) chromosomal haplotype of the linked duplication decreases nearly exponential to zero with time while that of the unlinked decreases quickly to an quasi-equilibrium; and this phenomenon is particularly profound under the HI model, because the quasi-equilibrium frequency of the original chromosomal haplotype (x(0)) under the HI model is higher than that under the DNR model. These results suggest that recombination and HI model might jointly contribute to the marked prolongation of T(r) even in a modest population. The prolonged T(r) and higher quasi-equilibrium frequency of the original allele at both duplicated loci might have offered more opportunities for the emergence of novel genes.