Tempo and Mode of Genome Evolution in the Budding Yeast Subphylum.

Budding yeasts (subphylum Saccharomycotina) are found in every biome and are as genetically diverse as plants or animals. To understand budding yeast evolution, we analyzed the genomes of 332 yeast species, including 220 newly sequenced ones, which represent nearly one-third of all known budding yeast diversity. Here, we establish a robust genus-level phylogeny comprising 12 major clades, infer the timescale of diversification from the Devonian period to the present, quantify horizontal gene transfer (HGT), and reconstruct the evolution of 45 metabolic traits and the metabolic toolkit of the budding yeast common ancestor (BYCA). We infer that BYCA was metabolically complex and chronicle the tempo and mode of genomic and phenotypic evolution across the subphylum, which is characterized by very low HGT levels and widespread losses of traits and the genes that control them. More generally, our results argue that reductive evolution is a major mode of evolutionary diversification.

Pumping iron: A multi-omics analysis of two extremophilic algae reveals iron economy management.

Marine algae are responsible for half of the world's primary productivity, but this critical carbon sink is often constrained by insufficient iron. One species of marine algae, Dunaliella tertiolecta, is remarkable for its ability to maintain photosynthesis and thrive in low-iron environments. A related species, Dunaliella salina Bardawil, shares this attribute but is an extremophile found in hypersaline environments. To elucidate how algae manage their iron requirements, we produced high-quality genome assemblies and transcriptomes for both species to serve as a foundation for a comparative multiomics analysis. We identified a host of iron-uptake proteins in both species, including a massive expansion of transferrins and a unique family of siderophore-iron-uptake proteins. Complementing these multiple iron-uptake routes, ferredoxin functions as a large iron reservoir that can be released by induction of flavodoxin. Proteomic analysis revealed reduced investment in the photosynthetic apparatus coupled with remodeling of antenna proteins by dramatic iron-deficiency induction of TIDI1, which is closely related but identifiably distinct from the chlorophyll binding protein, LHCA3. These combinatorial iron scavenging and sparing strategies make Dunaliella unique among photosynthetic organisms.

Behavioural differences underlie toxicity and predation variation in blooms of Prymnesium parvum.

Much of the evolutionary ecology of toxic algal blooms (TABs) remains unclear, including the role of algal toxins in the adaptive 'strategies' of TAB-forming species. Most eukaryotic TABs are caused by mixotrophs that augment autotrophy with organic nutrient sources, including competing algae (intraguild predation). We leverage the standing diversity of TABs formed by the toxic, invasive mixotroph Prymnesium parvum to identify cell-level behaviours involved in toxin-assisted predation using direct observations as well as comparisons between genetically distinct low- and high-toxicity isolates. Our results suggest that P. parvum toxins are primarily delivered at close range and promote subsequent prey capture/consumption. Surprisingly, we find opposite chemotactic preferences for organic (prey-derived) and inorganic nutrients between differentially toxic isolates, respectively, suggesting behavioural integration of toxicity and phagotrophy. Variation in toxicity may, therefore, reflect broader phenotypic integration of key traits that ultimately contribute to the remarkable flexibility, diversity, and success of invasive populations.

Conservation and diversification of HAIRY MERISTEM gene family in land plants.

The shoot apical meristems (SAMs) of land plants are crucial for plant growth and organ formation. In several angiosperms, the HAIRY MERISTEM (HAM) genes function as key regulators that control meristem development and stem cell homeostasis. To date, the origin and evolutionary history of the HAM family in land plants remains unclear. Potentially shared and divergent functions of HAM family members from angiosperms and non-angiosperms are also not known. In constructing a comprehensive phylogeny of the HAM family, we show that HAM proteins are widely present in land plants and that HAM proteins originated prior to the divergence of bryophytes. The HAM family was duplicated in a common ancestor of angiosperms, leading to two distinct groups: type I and type II. Type-II HAM members are widely present in angiosperms, whereas type-I HAM members were independently lost in different orders of monocots. Furthermore, HAM members from angiosperms and non-angiosperms (including bryophytes, lycophytes, ferns and gymnosperms) are able to replace the role of the type-II HAM genes in Arabidopsis, maintaining established SAMs and promoting the initiation of new stem cell niches. Our results uncover the conserved functions of HAM family members and reveal the conserved regulatory mechanisms underlying HAM expression patterning in meristems, providing insight into the evolution of key stem cell regulators in land plants.

Ecological generalism drives hyperdiversity of secondary metabolite gene clusters in xylarialean endophytes.

Although secondary metabolites are typically associated with competitive or pathogenic interactions, the high bioactivity of endophytic fungi in the Xylariales, coupled with their abundance and broad host ranges spanning all lineages of land plants and lichens, suggests that enhanced secondary metabolism might facilitate symbioses with phylogenetically diverse hosts. Here, we examined secondary metabolite gene clusters (SMGCs) across 96 Xylariales genomes in two clades (Xylariaceae s.l. and Hypoxylaceae), including 88 newly sequenced genomes of endophytes and closely related saprotrophs and pathogens. We paired genomic data with extensive metadata on endophyte hosts and substrates, enabling us to examine genomic factors related to the breadth of symbiotic interactions and ecological roles. All genomes contain hyperabundant SMGCs; however, Xylariaceae have increased numbers of gene duplications, horizontal gene transfers (HGTs) and SMGCs. Enhanced metabolic diversity of endophytes is associated with a greater diversity of hosts and increased capacity for lignocellulose decomposition. Our results suggest that, as host and substrate generalists, Xylariaceae endophytes experience greater selection to diversify SMGCs compared with more ecologically specialised Hypoxylaceae species. Overall, our results provide new evidence that SMGCs may facilitate symbiosis with phylogenetically diverse hosts, highlighting the importance of microbial symbioses to drive fungal metabolic diversity.

Horizontal Transfer of Bacterial Cytolethal Distending Toxin B Genes to Insects.

Horizontal gene transfer events have played a major role in the evolution of microbial species, but their importance in animals is less clear. Here, we report horizontal gene transfer of cytolethal distending toxin B (cdtB), prokaryotic genes encoding eukaryote-targeting DNase I toxins, into the genomes of vinegar flies (Diptera: Drosophilidae) and aphids (Hemiptera: Aphididae). We found insect-encoded cdtB genes are most closely related to orthologs from bacteriophage that infect Candidatus Hamiltonella defensa, a bacterial mutualistic symbiont of aphids that confers resistance to parasitoid wasps. In drosophilids, cdtB orthologs are highly expressed during the parasitoid-prone larval stage and encode a protein with ancestral DNase activity. We show that cdtB has been domesticated by diverse insects and hypothesize that it functions in defense against their natural enemies.

A Global Coexpression Network Approach for Connecting Genes to Specialized Metabolic Pathways in Plants.

Plants produce diverse specialized metabolites (SMs), but the genes responsible for their production and regulation remain largely unknown, hindering efforts to tap plant pharmacopeia. Given that genes comprising SM pathways exhibit environmentally dependent coregulation, we hypothesized that genes within a SM pathway would form tight associations (modules) with each other in coexpression networks, facilitating their identification. To evaluate this hypothesis, we used 10 global coexpression data sets, each a meta-analysis of hundreds to thousands of experiments, across eight plant species to identify hundreds of coexpressed gene modules per data set. In support of our hypothesis, 15.3 to 52.6% of modules contained two or more known SM biosynthetic genes, and module genes were enriched in SM functions. Moreover, modules recovered many experimentally validated SM pathways, including all six known to form biosynthetic gene clusters (BGCs). In contrast, bioinformatically predicted BGCs (i.e., those lacking an associated metabolite) were no more coexpressed than the null distribution for neighboring genes. These results suggest that most predicted plant BGCs are not genuine SM pathways and argue that BGCs are not a hallmark of plant specialized metabolism. We submit that global gene coexpression is a rich, largely untapped resource for discovering the genetic basis and architecture of plant natural products.

Drivers of genetic diversity in secondary metabolic gene clusters within a fungal species.

Filamentous fungi produce a diverse array of secondary metabolites (SMs) critical for defense, virulence, and communication. The metabolic pathways that produce SMs are found in contiguous gene clusters in fungal genomes, an atypical arrangement for metabolic pathways in other eukaryotes. Comparative studies of filamentous fungal species have shown that SM gene clusters are often either highly divergent or uniquely present in one or a handful of species, hampering efforts to determine the genetic basis and evolutionary drivers of SM gene cluster divergence. Here, we examined SM variation in 66 cosmopolitan strains of a single species, the opportunistic human pathogen Aspergillus fumigatus. Investigation of genome-wide within-species variation revealed 5 general types of variation in SM gene clusters: nonfunctional gene polymorphisms; gene gain and loss polymorphisms; whole cluster gain and loss polymorphisms; allelic polymorphisms, in which different alleles corresponded to distinct, nonhomologous clusters; and location polymorphisms, in which a cluster was found to differ in its genomic location across strains. These polymorphisms affect the function of representative A. fumigatus SM gene clusters, such as those involved in the production of gliotoxin, fumigaclavine, and helvolic acid as well as the function of clusters with undefined products. In addition to enabling the identification of polymorphisms, the detection of which requires extensive genome-wide synteny conservation (e.g., mobile gene clusters and nonhomologous cluster alleles), our approach also implicated multiple underlying genetic drivers, including point mutations, recombination, and genomic deletion and insertion events as well as horizontal gene transfer from distant fungi. Finally, most of the variants that we uncover within A. fumigatus have been previously hypothesized to contribute to SM gene cluster diversity across entire fungal classes and phyla. We suggest that the drivers of genetic diversity operating within a fungal species shown here are sufficient to explain SM cluster macroevolutionary patterns.

Genome-wide analysis of Fusarium verticillioides reveals inter-kingdom contribution of horizontal gene transfer to the expansion of metabolism.

Horizontal gene transfer (HGT) is believed to shape genomes by facilitating the rapid acquisition of adaptive traits. We hypothesized that the economically important fungus Fusarium verticillioides is an excellent candidate for investigating the potential impact of HGT on the expansion of metabolic activities given its soilborne nature and versatile lifestyle as both a symptomless endophyte as well as a maize pathogen. To test this hypothesis, we used a phylogenomic pipeline followed by manual curation to perform a genome-wide identification of inter-kingdom derived HGT events. We found strong support for 36 genes in F. verticillioides putatively acquired from bacteria. Functional enrichment assessment of these 36 candidates suggested HGT potentially influenced several biochemical activities, including lysine, glycine and nitrogen metabolism. The expression of 25 candidate HGT genes was detected among RNA-Seq datasets from normal and various stress-related growth conditions, thus indicating potential functionality. FVEG_10494, one of the HGT candidates with homologs in only a few Fusarium species, was highly and specifically up-regulated under nitric oxide (NO) challenge. Functional analysis of FVEG_10494 suggests the gene moderately enhanced NO-triggered protective responses and suppressed expression of the F. verticillioides secondary metabolism gene cluster responsible for production of fusarin C. Overall, our global analysis of HGT events in F. verticillioides identified a well-supported set of transferred genes, providing further evidence that HGT offers a mechanism by which fungi can expand their metabolic capabilities, which in turn may enhance their adaptive strategies.

Horizontally acquired genes in early-diverging pathogenic fungi enable the use of host nucleosides and nucleotides.

Horizontal gene transfer (HGT) among bacteria, archaea, and viruses is widespread, but the extent of transfers from these lineages into eukaryotic organisms is contentious. Here we systematically identify hundreds of genes that were likely acquired horizontally from a variety of sources by the early-diverging fungal phyla Microsporidia and Cryptomycota. Interestingly, the Microsporidia have acquired via HGT several genes involved in nucleic acid synthesis and salvage, such as those encoding thymidine kinase (TK), cytidylate kinase, and purine nucleotide phosphorylase. We show that these HGT-derived nucleic acid synthesis genes tend to function at the interface between the metabolic networks of the host and pathogen. Thus, these genes likely play vital roles in diversifying the useable nucleic acid components available to the intracellular parasite, often through the direct capture of resources from the host. Using an in vivo viability assay, we also demonstrate that one of these genes, TK, encodes an enzyme that is capable of activating known prodrugs to their active form, which suggests a possible treatment route for microsporidiosis. We further argue that interfacial genes with well-understood activities, especially those horizontally transferred from bacteria or viruses, could provide medical treatments for microsporidian infections.

Comparative genomics of biotechnologically important yeasts.

Ascomycete yeasts are metabolically diverse, with great potential for biotechnology. Here, we report the comparative genome analysis of 29 taxonomically and biotechnologically important yeasts, including 16 newly sequenced. We identify a genetic code change, CUG-Ala, in Pachysolen tannophilus in the clade sister to the known CUG-Ser clade. Our well-resolved yeast phylogeny shows that some traits, such as methylotrophy, are restricted to single clades, whereas others, such as l-rhamnose utilization, have patchy phylogenetic distributions. Gene clusters, with variable organization and distribution, encode many pathways of interest. Genomics can predict some biochemical traits precisely, but the genomic basis of others, such as xylose utilization, remains unresolved. Our data also provide insight into early evolution of ascomycetes. We document the loss of H3K9me2/3 heterochromatin, the origin of ascomycete mating-type switching, and panascomycete synteny at the MAT locus. These data and analyses will facilitate the engineering of efficient biosynthetic and degradative pathways and gateways for genomic manipulation.

Examining the evolution of the regulatory circuit controlling secondary metabolism and development in the fungal genus Aspergillus.

Filamentous fungi produce diverse secondary metabolites (SMs) essential to their ecology and adaptation. Although each SM is typically produced by only a handful of species, global SM production is governed by widely conserved transcriptional regulators in conjunction with other cellular processes, such as development. We examined the interplay between the taxonomic narrowness of SM distribution and the broad conservation of global regulation of SM and development in Aspergillus, a diverse fungal genus whose members produce well-known SMs such as penicillin and gliotoxin. Evolutionary analysis of the 2,124 genes comprising the 262 SM pathways in four Aspergillus species showed that most SM pathways were species-specific, that the number of SM gene orthologs was significantly lower than that of orthologs in primary metabolism, and that the few conserved SM orthologs typically belonged to non-homologous SM pathways. RNA sequencing of two master transcriptional regulators of SM and development, veA and mtfA, showed that the effects of deletion of each gene, especially veA, on SM pathway regulation were similar in A. fumigatus and A. nidulans, even though the underlying genes and pathways regulated in each species differed. In contrast, examination of the role of these two regulators in development, where 94% of the underlying genes are conserved in both species showed that whereas the role of veA is conserved, mtfA regulates development in the homothallic A. nidulans but not in the heterothallic A. fumigatus. Thus, the regulation of these highly conserved developmental genes is divergent, whereas-despite minimal conservation of target genes and pathways-the global regulation of SM production is largely conserved. We suggest that the evolution of the transcriptional regulation of secondary metabolism in Aspergillus represents a novel type of regulatory circuit rewiring and hypothesize that it has been largely driven by the dramatic turnover of the target genes involved in the process.

Dynamic Evolution of Nitric Oxide Detoxifying Flavohemoglobins, a Family of Single-Protein Metabolic Modules in Bacteria and Eukaryotes.

Due to their functional independence, proteins that comprise standalone metabolic units, which we name single-protein metabolic modules, may be particularly prone to gene duplication (GD) and horizontal gene transfer (HGT). Flavohemoglobins (flavoHbs) are prime examples of single-protein metabolic modules, detoxifying nitric oxide (NO), a ubiquitous toxin whose antimicrobial properties many life forms exploit, to nitrate, a common source of nitrogen for organisms. FlavoHbs appear widespread in bacteria and have been identified in a handful of microbial eukaryotes, but how the distribution of this ecologically and biomedically important protein family evolved remains unknown. Reconstruction of the evolutionary history of 3,318 flavoHb protein sequences covering the family's known diversity showed evidence of recurrent HGT at multiple evolutionary scales including intrabacterial HGT, as well as HGT from bacteria to eukaryotes. One of the most striking examples of HGT is the acquisition of a flavoHb by the dandruff- and eczema-causing fungus Malassezia from Corynebacterium Actinobacteria, a transfer that growth experiments show is capable of mediating NO resistance in fungi. Other flavoHbs arose via GD; for example, many filamentous fungi possess two flavoHbs that are differentially targeted to the cytosol and mitochondria, likely conferring protection against external and internal sources of NO, respectively. Because single-protein metabolic modules such as flavoHb function independently, readily undergo GD and HGT, and are frequently involved in organismal defense and competition, we suggest that they represent "plug-and-play" proteins for ecological arms races.

The evolution of fungal metabolic pathways.

Fungi contain a remarkable range of metabolic pathways, sometimes encoded by gene clusters, enabling them to digest most organic matter and synthesize an array of potent small molecules. Although metabolism is fundamental to the fungal lifestyle, we still know little about how major evolutionary processes, such as gene duplication (GD) and horizontal gene transfer (HGT), have interacted with clustered and non-clustered fungal metabolic pathways to give rise to this metabolic versatility. We examined the synteny and evolutionary history of 247,202 fungal genes encoding enzymes that catalyze 875 distinct metabolic reactions from 130 pathways in 208 diverse genomes. We found that gene clustering varied greatly with respect to metabolic category and lineage; for example, clustered genes in Saccharomycotina yeasts were overrepresented in nucleotide metabolism, whereas clustered genes in Pezizomycotina were more common in lipid and amino acid metabolism. The effects of both GD and HGT were more pronounced in clustered genes than in their non-clustered counterparts and were differentially distributed across fungal lineages; specifically, GD, which was an order of magnitude more abundant than HGT, was most frequently observed in Agaricomycetes, whereas HGT was much more prevalent in Pezizomycotina. The effect of HGT in some Pezizomycotina was particularly strong; for example, we identified 111 HGT events associated with the 15 Aspergillus genomes, which sharply contrasts with the 60 HGT events detected for the 48 genomes from the entire Saccharomycotina subphylum. Finally, the impact of GD within a metabolic category was typically consistent across all fungal lineages, whereas the impact of HGT was variable. These results indicate that GD is the dominant process underlying fungal metabolic diversity, whereas HGT is episodic and acts in a category- or lineage-specific manner. Both processes have a greater impact on clustered genes, suggesting that metabolic gene clusters represent hotspots for the generation of fungal metabolic diversity.

Dinoflagellate genome evolution.

The dinoflagellates are an ecologically important group of microbial eukaryotes that have evolved many novel genomic characteristics. They possess some of the largest nuclear genomes among eukaryotes arranged on permanently condensed liquid-crystalline chromosomes. Recent advances have revealed the presence of genes arranged in tandem arrays, trans-splicing of messenger RNAs, and a reduced role for transcriptional regulation compared to other eukaryotes. In contrast, the mitochondrial and plastid genomes have the smallest gene content among functional eukaryotic organelles. Dinoflagellate biology and genome evolution have been dramatically influenced by lateral transfer of individual genes and large-scale transfer of genes through endosymbiosis. Next-generation sequencing technologies have only recently made genome-scale analyses of these organisms possible, and these new methods are helping researchers better understand the biology and evolution of this enigmatic group of eukaryotes.

Insights into transcriptional changes that accompany organelle sequestration from the stolen nucleus of Mesodinium rubrum.

BACKGROUND: Organelle retention is a form of mixotrophy that allows organisms to reap metabolic benefits similar to those of photoautotrophs through capture of algal prey and sequestration of their plastids. Mesodinium rubrum is an abundant and broadly distributed photosynthetic marine ciliate that steals organelles from cryptophyte algae, such as Geminigera cryophila. M. rubrum is unique from most other acquired phototrophs because it also steals a functional nucleus that facilitates genetic control of sequestered plastids and other organelles. We analyzed changes in G. cryophila nuclear gene expression and transcript abundance after its incorporation into the cellular architecture of M. rubrum as an initial step towards understanding this complex system. METHODS: We compared Illumina-generated transcriptomes of the cryptophyte Geminigera cryophila as a free-living cell and as a sequestered nucleus in M. rubrum to identify changes in protein abundance and gene expression. After KEGG annotation, proteins were clustered by functional categories, which were evaluated for over- or under-representation in the sequestered nucleus. Similarly, coding sequences were grouped by KEGG categories/pathways, which were then evaluated for over- or under-expression via read count strategies. RESULTS: At the time of sampling, the global transcriptome of M. rubrum was dominated (~58-62 %) by transcription from its stolen nucleus. A comparison of transcriptomes from free-living G. cryophila cells to those of the sequestered nucleus revealed a decrease in gene expression and transcript abundance for most functional protein categories within the ciliate. However, genes coding for proteins involved in photosynthesis, oxidative stress reduction, and several other metabolic pathways revealed striking exceptions to this general decline. CONCLUSIONS: Major changes in G. cryophila transcript expression after sequestration by M. rubrum and the ciliate's success as a photoautotroph imply some level of control or gene regulation by the ciliate and at the very least reflect a degree of coordination between host and foreign organelles. Intriguingly, cryptophyte genes involved in protein transport are significantly under-expressed in M. rubrum, implicating a role for the ciliate's endomembrane system in targeting cryptophyte proteins to plastid complexes. Collectively, this initial portrait of an acquired transcriptome within a dynamic and ecologically successful ciliate highlights the remarkable cellular and metabolic chimerism of this system.

Evolution of saxitoxin synthesis in cyanobacteria and dinoflagellates.

Dinoflagellates produce a variety of toxic secondary metabolites that have a significant impact on marine ecosystems and fisheries. Saxitoxin (STX), the cause of paralytic shellfish poisoning, is produced by three marine dinoflagellate genera and is also made by some freshwater cyanobacteria. Genes involved in STX synthesis have been identified in cyanobacteria but are yet to be reported in the massive genomes of dinoflagellates. We have assembled comprehensive transcriptome data sets for several STX-producing dinoflagellates and a related non-toxic species and have identified 265 putative homologs of 13 cyanobacterial STX synthesis genes, including all of the genes directly involved in toxin synthesis. Putative homologs of four proteins group closely in phylogenies with cyanobacteria and are likely the functional homologs of sxtA, sxtG, and sxtB in dinoflagellates. However, the phylogenies do not support the transfer of these genes directly between toxic cyanobacteria and dinoflagellates. SxtA is split into two proteins in the dinoflagellates corresponding to the N-terminal portion containing the methyltransferase and acyl carrier protein domains and a C-terminal portion with the aminotransferase domain. Homologs of sxtB and N-terminal sxtA are present in non-toxic strains, suggesting their functions may not be limited to saxitoxin production. Only homologs of the C-terminus of sxtA and sxtG were found exclusively in toxic strains. A more thorough survey of STX+ dinoflagellates will be needed to determine if these two genes may be specific to SXT production in dinoflagellates. The A. tamarense transcriptome does not contain homologs for the remaining STX genes. Nevertheless, we identified candidate genes with similar predicted biochemical activities that account for the missing functions. These results suggest that the STX synthesis pathway was likely assembled independently in the distantly related cyanobacteria and dinoflagellates, although using some evolutionarily related proteins. The biological role of STX is not well understood in either cyanobacteria or dinoflagellates. However, STX production in these two ecologically distinct groups of organisms suggests that this toxin confers a benefit to producers that we do not yet fully understand.

The impact of automated filtering of BLAST-determined homologs in the phylogenetic detection of horizontal gene transfer from a transcriptome assembly.

Phylomes (comprehensive sets of gene phylogenies for organisms) are built to investigate fundamental questions in genomics and evolutionary biology, such as those pertaining to the detection and characterization of horizontal gene transfer in microbes. To address these questions, phylome construction demands rigorous yet efficient phylogenetic methods. Currently, many sequence alignment and tree-building models can analyze several thousands of genes in a high-throughput manner. However, the phylogenetics is complicated by variability in sequence divergence and different taxon sampling among genes. In addition, homolog selection for automated approaches often relies on arbitrary sequence similarity thresholds that are likely inappropriate for all genes in a genome. To investigate the effects of automated homolog selection on the detection of horizontal gene transfer using phylogenomics, we constructed the phylome of a transcriptome assembly of Alexandrium tamarense, a microbial eukaryote with a history of horizontal and endosymbiotic gene transfer, using seven sequence similarity thresholds for selecting putative homologs to be included in phylogenetic analyses. We show that no single threshold recovered informative trees for the majority of A. tamarense unigenes compared to the pooled results from all pipeline iterations. As much as 29% of trees built could have misleading phylogenetic relationships that appear biased in favor of those otherwise indicative of horizontal gene transfer. Perhaps worse, nearly half of the unigenes were represented by a single tree built at just one threshold, making it difficult to assess the validity of phylogenetic relationships recovered in these cases. However, combining the results from several pipeline iterations maximizes the number of informative phylogenies. Moreover, when the same phylogenetic relationship for a given unigene is recovered in multiple pipeline iterations, conclusions regarding gene origin are more robust to methodological artifact. Using these methods, the majority of A. tamarense unigenes showed evolutionary relationships indicative of vertical inheritance. Nevertheless, many other unigenes revealed diverse phylogenetic associations, suggestive of possible gene transfer. This analysis suggests that caution should be used when interpreting the results from phylogenetic pipelines implementing a single similarity threshold. Our approach is a practical method to mitigate the problems associated with automated sequence selection in phylogenomics.

Candida glabrata maintains two Hap1 homologs, Zcf27 and Zcf4, for distinct roles in ergosterol gene regulation to mediate sterol homeostasis under azole and hypoxic conditions.

Candida glabrata exhibits innate resistance to azole antifungal drugs but also has the propensity to rapidly develop clinical drug resistance. Azole drugs, which target Erg11, is one of the three major classes of antifungals used to treat Candida infections. Despite their widespread use, the mechanism controlling azole-induced ERG gene expression and drug resistance in C. glabrata has primarily revolved around Upc2 and/or Pdr1. In this study, we determined the function of two zinc cluster transcription factors, Zcf27 and Zcf4, as direct but distinct regulators of ERG genes. Our phylogenetic analysis revealed C. glabrata Zcf27 and Zcf4 as the closest homologs to Saccharomyces cerevisiae Hap1. Hap1 is a known zinc cluster transcription factor in S. cerevisiae in controlling ERG gene expression under aerobic and hypoxic conditions. Interestingly, when we deleted HAP1 or ZCF27 in either S. cerevisiae or C. glabrata, respectively, both deletion strains showed altered susceptibility to azole drugs, whereas the strain deleted for ZCF4 did not exhibit azole susceptibility. We also determined that the increased azole susceptibility in a zcf27Δ strain is attributed to decreased azole-induced expression of ERG genes, resulting in decreased levels of total ergosterol. Surprisingly, Zcf4 protein expression is barely detected under aerobic conditions but is specifically induced under hypoxic conditions. However, under hypoxic conditions, Zcf4 but not Zcf27 was directly required for the repression of ERG genes. This study provides the first demonstration that Zcf27 and Zcf4 have evolved to serve distinct roles allowing C. glabrata to adapt to specific host and environmental conditions.