The exact distributions of F(IS) under partial asexuality in small finite populations with mutation

Reproductive systems like partial asexuality participate to shape the evolution of genetic diversity within populations, which is often quantified by the inbreeding coefficient F IS. Understanding how those mating systems impact the possible distributions of F IS values in theoretical populations helps to unravel forces shaping the evolution of real populations. We proposed a population genetics model based on genotypic states in a finite population with mutation. For populations with less than 400 individuals, we assessed the impact of the rates of asexuality on the full exact distributions of F IS, the probabilities of positive and negative F IS, the probabilities of fixation and the probabilities to observe changes in the sign of F IS over one generation. After an infinite number of generations, we distinguished three main patterns of effects of the rates of asexuality on genetic diversity that also varied according to the interactions of mutation and genetic drift. Even rare asexual events in mainly sexual populations impacted the balance between negative and positive F IS and the occurrence of extreme values. It also drastically modified the probability to change the sign of F IS value at one locus over one generation. When mutation prevailed over genetic drift, increasing rates of asexuality continuously increased the variance of F IS that reached its highest value in fully asexual populations. In consequence, even ancient asexual populations showed the entire F IS spectrum, including strong positive F IS. The prevalence of heterozygous loci only occurred in full asexual populations when genetic drift dominated.