Coalescence and linkage disequilibrium in facultatively sexual diploids
AbstractUnder neutrality, linkage disequilibrium (LD) results from physically linked sites having non-independent coalescent histories. In obligately sexual organisms, meiotic recombination is the dominant force separating linked variants from one another, and thus in determining the decay of LD with physical distance. In facultatively sexual diploid organisms that principally reproduce clonally, mechanisms of mitotic exchange are expected to become relatively more important in shaping LD. Here we outline mathematical and computational models of a facultative-sex coalescent process that includes meiotic and mitotic recombination, via both crossovers and gene conversion, to determine how LD is affected with facultative sex. We demonstrate that the degree to which LD is broken down by meiotic recombination simply scales with the probability of sex if it is sufficiently high (much greater than 1/N for N the population size). However, with very rare sex (occurring with frequency on the order of 1/N), mitotic gene conversion plays a particularly important and complicated role because it both breaks down associations between sites and removes within-individual diversity. Strong population structure under rare sex leads to lower average LD values than in panmictic populations, due to the influence of low-frequency polymorphisms created by allelic sequence divergence acting in individual subpopulations. These analyses provide information on how to interpret observed LD patterns in facultative sexuals, and determine what genomic forces are likely to shape them.