Candida species are the most common cause of opportunistic fungal infection worldwide. We report the genome sequences of six Candida species and compare these and related pathogens and nonpathogens. There are significant expansions of cell wall, secreted, and transporter gene families in pathogenic species, suggesting adaptations associated with virulence. Large genomic tracts are homozygous in three diploid species, possibly resulting from recent recombination events. Surprisingly, key components of the mating and meiosis pathways are missing from several species. These include major differences at the Mating-type loci (MTL); Lodderomyces elongisporus lacks MTL, and components of the a1/alpha2 cell identity determinant were lost in other species, raising questions about how mating and cell types are controlled. Analysis of the CUG leucine to serine genetic code change reveals that 99% of ancestral CUG codons were erased and new ones arose elsewhere. Lastly, we revise the C. albicans gene catalog, identifying many new genes.

Background: To date, most fungal phylogenies have been derived from single gene comparisons, or from concatenated alignments of a small number of genes. The increase in fungal genome sequencing presents an opportunity to reconstruct evolutionary events using entire genomes. As a tool for future comparative, phylogenomic and phylogenetic studies, we used both supertrees and concatenated alignments to infer relationships between 42 species of fungi for which complete genome sequences are available.

Ace2 transcription factor family genes are found in many fungal genomes and are required for regulation of expression of genes involved in cell separation. We used transcriptional profiling to identify the targets of Ace2 in Candida albicans, and we show that these include several cell wall components, such as glucanases and glycosylphosphatidylinositol-anchored proteins. Expression is downregulated in ace2 deletion mutants in both yeast and hyphal cells. In addition, deleting ace2 results in dramatic changes in expression of metabolic pathways. Expression of glycolytic enzymes is reduced, while expression of respiratory genes (including those involved in the tricarboxylic acid cycle, oxidative phosphorylation, and ATP synthesis) is increased. Similar changes occur in both yeast and hyphal cells. In contrast, genes required for acetyl-coenzyme A and lipid metabolism are upregulated in an ace2 deletion mutant grown predominantly as yeast cells but are downregulated in hyphae. These results suggest that in wild-type strains, Ace2 acts to increase glycolysis and reduce respiration. This is supported by the observation that deleting ace2 results in increased resistance to antimycin A, a drug that inhibits respiration. We also show that Ace2 is required for filamentation in response to low oxygen concentrations (hypoxia). We suggest that filamentation is induced in wild-type cells by reducing respiration (using low oxygen or respiratory drugs) and that mutants with increased respiratory activity fail to undergo filamentation under these conditions.In Saccharomyces cerevisiae, ACE2 and SWI5 encode two of a set of nine transcription factors that control the mitotic cell cycle (40). The nine regulators are SBF (Swi4 and Swi6) and MBF (Swi6 and Mbp1), which control the expression of genes in G 1 and S phase (18); Ndd1, Fkh1, Fkh2, and Mcm1, which regulate the expression of genes at the G 2 /M border; and Swi5 and Ace2 (and Mcm1), which control the expression of genes in late M and early G 1 . The transcription factors act in a cascade: SBF and MBF regulate expression of NDD1, the G 2 activators control expression of ACE2 and SWI5, and these in turn are required for exit from mitosis and subsequent activation of SBF and MBF.Ace2 and Swi5 share many functional similarities in S. cerevisiae. The proteins are 37% identical and have the same DNA-binding sites in vitro (29). They also regulate the expression of many of the same genes (9). However, there are substantial differences. There is only ϳ22% overlap between the groups of genes regulated by the two factors (40). Swi5 remains cytoplasmic until the end of M phase, when it enters the nucleus (33). Ace2 also enters the nucleus at the end of mitosis (35) but is rapidly exported from (or degraded in) the nucleus of the mother cell and remains only in the nucleus of the daughter cell (7, 46). Localization is regulated by components of the RAM pathway (34, 38). Expression of many Ace2 targets is therefore restricted to the daughter cell, and one of their main functions is to e...

The G + C content at synonymous codon positions (GC3s) in genes varies along chromosomes in most eukaryotes. In Saccharomyces cerevisiae, regions of high GC3s are correlated with recombination hot spots, probably due to biased gene conversion. Here we examined how GC3s differs among groups of related yeast species in the Saccharomyces and Candida clades. The chromosomal locations of GC3s peaks and troughs are conserved among four Saccharomyces species, but we find that there have been highly consistent small shifts in their GC3s values. For instance, 84% of all S. cerevisiae genes have a lower GC3s value than their S. bayanus orthologs. There are extensive interspecies differences in the Candida clade both in the median value of GC3s (ranging from 22% to 49%) and in the variance of GC3s among genes. In three species—Candida lusitaniae, Pichia stipitis, and Yarrowia lipolytica—there is one region on each chromosome in which GC3s is markedly reduced. We propose that these GC-poor troughs indicate the positions of centromeres because in Y. lipolytica they coincide with the five experimentally identified centromeres. In P. stipitis, the troughs contain clusters of the retrotransposon Tps5. Likewise, in Debaryomyces hansenii, there is one cluster of the retrotransposon Tdh5 per chromosome, and all these clusters are located in GC-poor troughs. Locally reduced G + C content around centromeres is consistent with a model in which G + C content correlates with recombination rate, and recombination is suppressed around centromeres, although the troughs are unexpectedly wide (100–300 kb).

BackgroundTo date very few incidences of interdomain gene transfer into fungi have been identified. Here, we used the emerging genome sequences of Candida albicans WO-1, Candida tropicalis, Candida parapsilosis, Clavispora lusitaniae, Pichia guilliermondii, and Lodderomyces elongisporus to identify recent interdomain HGT events. We refer to these as CTG species because they translate the CTG codon as serine rather than leucine, and share a recent common ancestor.ResultsPhylogenetic and syntenic information infer that two C. parapsilosis genes originate from bacterial sources. One encodes a putative proline racemase (PR). Phylogenetic analysis also infers that there were independent transfers of bacterial PR enzymes into members of the Pezizomycotina, and protists. The second HGT gene in C. parapsilosis belongs to the phenazine F (PhzF) superfamily. Most CTG species also contain a fungal PhzF homolog. Our phylogeny suggests that the CTG homolog originated from an ancient HGT event, from a member of the proteobacteria. An analysis of synteny suggests that C. parapsilosis has lost the endogenous fungal form of PhzF, and subsequently reacquired it from a proteobacterial source. There is evidence that Schizosaccharomyces pombe and Basidiomycotina also obtained a PhzF homolog through HGT.ConclusionOur search revealed two instances of well-supported HGT from bacteria into the CTG clade, both specific to C. parapsilosis. Therefore, while recent interkingdom gene transfer has taken place in the CTG lineage, its occurrence is rare. However, our analysis will not detect ancient gene transfers, and we may have underestimated the global extent of HGT into CTG species.