Models for the Evolution of GC Content in Asexual Fungi Candida albicans and C. dubliniensis

Although guanine–cytosine (GC)-biased gene conversion (gBGC) following meiotic recombination seems the most probable mechanism accounting for large-scale variations in GC content for many eukaryotes, it cannot explain such variations for organisms belonging to ancient asexual lineages, such as the pathogenic fungi Candida albicans and C. dubliniensis. Analysis of the substitution patterns for these two species reveals a strong anticorrelation between the synonymous transition rates at third codon positions. I propose two models that can account for this observation. According to the first model, the evolution of GC content is driven by gBGC linked to mitotic recombination, either associated with parasexuality or with damage repair. Variations in the GC content thus reflect variations in the strength of gBGC, presumably variations in the mitotic recombination rate. According to the second model, the evolution of GC content is driven by misincorporation errors during the process of DNA replication in S phase. This model proposes that variations in GC content are due to variations in the proportions of dCTPs and dGTPs at the time when sequences are replicated. Experimental data regarding mitotic recombination rates or the variations of dCTPs and dGTPs during S phase are required to validate definitively one of the two models, but in any case, the fit of the models to the data suggests that C. albicans and C. dubliniensis constitute so far unique examples of GC content evolution driven either by mitotic recombination or replicative errors.