The tree versus the forest: the fungal tree of life and the topological diversity within the yeast phylome

Characterization of a selective, iron-chelating antifungal compound that disrupts fungal metabolism and synergizes with fluconazole

Fungal infections are a growing global health concern due to the limited number of available antifungal therapies as well as the emergence of fungi that are resistant to first-line antimicrobials, particularly azoles and echinocandins. Development of novel, selective antifungal therapies is challenging due to similarities between fungal and mammalian cells. An attractive source of potential antifungal treatments is provided by ecological niches co-inhabited by bacteria, fungi, and multicellular organisms, where complex relationships between multiple organisms have resulted in evolution of a wide variety of selective antimicrobials. Here, we characterized several analogs of one such natural compound, collismycin A. We show that NR-6226C has antifungal activity against several pathogenic Candida species, including C. albicans and C. glabrata, whereas it only has little toxicity against mammalian cells. Mechanistically, NR-6226C selectively chelates iron, which is a limiting factor for pathogenic fungi during infection. As a result, NR-6226C treatment causes severe mitochondrial dysfunction, leading to formation of reactive oxygen species, metabolic reprogramming, and a severe reduction in ATP levels. Using an in vivo model for fungal infections, we show that NR-6226C significantly increases survival of Candida-infected Galleria mellonella larvae. Finally, our data indicate that NR-6226C synergizes strongly with fluconazole in inhibition of C. albicans. Taken together, NR-6226C is a promising antifungal compound that acts by chelating iron and disrupting mitochondrial functions.IMPORTANCEDrug-resistant fungal infections are an emerging global threat, and pan-resistance to current antifungal therapies is an increasing problem. Clearly, there is a need for new antifungal drugs. In this study, we characterized a novel antifungal agent, the collismycin analog NR-6226C. NR-6226C has a favorable toxicity profile for human cells, which is essential for further clinical development. We unraveled the mechanism of action of NR-6226C and found that it disrupts iron homeostasis and thereby depletes fungal cells of energy. Importantly, NR-6226C strongly potentiates the antifungal activity of fluconazole, thereby providing inroads for combination therapy that may reduce or prevent azole resistance. Thus, NR-6226C is a promising compound for further development into antifungal treatment.

Copper Acts Synergistically With Fluconazole in Candida glabrata by Compromising Drug Efflux, Sterol Metabolism, and Zinc Homeostasis

The synergistic combinations of drugs are promising strategies to boost the effectiveness of current antifungals and thus prevent the emergence of resistance. In this work, we show that copper and the antifungal fluconazole act synergistically against Candida glabrata, an opportunistic pathogenic yeast intrinsically tolerant to fluconazole. Analyses of the transcriptomic profile of C. glabrata after the combination of copper and fluconazole showed that the expression of the multidrug transporter gene CDR1 was decreased, suggesting that fluconazole efflux could be affected. In agreement, we observed that copper inhibits the transactivation of Pdr1, the transcription regulator of multidrug transporters and leads to the intracellular accumulation of fluconazole. Copper also decreases the transcriptional induction of ergosterol biosynthesis (ERG) genes by fluconazole, which culminates in the accumulation of toxic sterols. Co-treatment of cells with copper and fluconazole should affect the function of proteins located in the plasma membrane, as several ultrastructural alterations, including irregular cell wall and plasma membrane and loss of cell wall integrity, were observed. Finally, we show that the combination of copper and fluconazole downregulates the expression of the gene encoding the zinc-responsive transcription regulator Zap1, which possibly, together with the membrane transporters malfunction, generates zinc depletion. Supplementation with zinc reverts the toxic effect of combining copper with fluconazole, underscoring the importance of this metal in the observed synergistic effect. Overall, this work, while unveiling the molecular basis that supports the use of copper to enhance the effectiveness of fluconazole, paves the way for the development of new metal-based antifungal strategies.

Sirtuins in Epigenetic Silencing and Control of Gene Expression in Model and Pathogenic Fungi

Fungi, including yeasts, molds, and mushrooms, proliferate on decaying matter and then adopt quiescent forms once nutrients are depleted. This review explores how fungi use sirtuin deacetylases to sense and respond appropriately to changing nutrients. Because sirtuins are NAD⁺-dependent deacetylases, their activity is sensitive to intracellular NAD⁺ availability. This allows them to transmit information about a cell's metabolic state on to the biological processes they influence. Fungal sirtuins are primarily known to deacetylate histones, repressing transcription and modulating genome stability. Their target genes include those involved in NAD⁺ homeostasis, metabolism, sporulation, secondary metabolite production, and virulence traits of pathogenic fungi. By targeting different genes over evolutionary time, sirtuins serve as rewiring points that allow organisms to evolve novel responses to low NAD⁺ stress by bringing relevant biological processes under the control of sirtuins. Expected final online publication date for the Annual Review of Microbiology, Volume 76 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

SpeciesRax: A Tool for Maximum Likelihood Species Tree Inference from Gene Family Trees under Duplication, Transfer, and Loss

Species tree inference from gene family trees is becoming increasingly popular because it can account for discordance between the species tree and the corresponding gene family trees. In particular, methods that can account for multiple-copy gene families exhibit potential to leverage paralogy as informative signal. At present, there does not exist any widely adopted inference method for this purpose. Here, we present SpeciesRax, the first maximum likelihood method that can infer a rooted species tree from a set of gene family trees and can account for gene duplication, loss, and transfer events. By explicitly modeling events by which gene trees can depart from the species tree, SpeciesRax leverages the phylogenetic rooting signal in gene trees. SpeciesRax infers species tree branch lengths in units of expected substitutions per site and branch support values via paralogy-aware quartets extracted from the gene family trees. Using both empirical and simulated data sets we show that SpeciesRax is at least as accurate as the best competing methods while being one order of magnitude faster on large data sets at the same time. We used SpeciesRax to infer a biologically plausible rooted phylogeny of the vertebrates comprising 188 species from 31,612 gene families in 1 h using 40 cores. SpeciesRax is available under GNU GPL at https://github.com/BenoitMorel/GeneRax and on BioConda.

Zap1 is required for Candida glabrata response to fluconazole

The increasing prevalence of fluconazole-resistant clinical isolates of Candida spp. strongly hinders the widespread use of the drug. To tackle this problem, great efforts have been made to fully understand the fungal response to fluconazole. In this work, we show that the role of Zap1 in Candida glabrata goes beyond regulating yeast adaptation to zinc deficiency. In line with our previous observation that deletion of ZAP1 makes yeast cells more sensitive to fluconazole, we found that the mutant CgΔzap1 accumulates higher levels of the drug, which correlates well with its lower levels of ergosterol. Surprisingly, Zap1 is a negative regulator of the drug efflux transporter gene CDR1 and of its regulator, PDR1. The apparent paradox of drug accumulation in cells where genes encoding transporters relevant for drug extrusion are being overexpressed led us to postulate that their activity could be impaired. In agreement, Zap1-depleted cells present, in addition to decreased ergosterol levels, an altered composition of membrane phospholipids, which together should impact membrane function and impair the detoxification of fluconazole. Overall, our study brings to light Zap1 as an important hub in Candida glabrata response to fluconazole.

Analysis of Volatile Molecules Present in the Secretome of the Fungal Pathogen Candida glabrata

Candida albicans, Candida glabrata, Candida parapsilosis and Candida tropicalis are the four most common human fungal pathogens isolated that can cause superficial and invasive infections. It has been shown that specific metabolites present in the secretomes of these fungal pathogens are important for their virulence. C. glabrata is the second most common isolate world-wide and has an innate resistance to azoles, xenobiotics and oxidative stress that allows this fungal pathogen to evade the immune response and persist within the host. Here, we analyzed and compared the C. glabrata secretome with those of C. albicans, C. parapsilosis, C. tropicalis and the non-pathogenic yeast Saccharomyces cerevisiae. In C. glabrata, we identified a different number of metabolites depending on the growth media: 12 in synthetic complete media (SC), 27 in SC-glutamic acid and 23 in rich media (YPD). C. glabrata specific metabolites are 1-dodecene (0.09 ± 0.11%), 2,5-dimethylundecane (1.01 ± 0.19%), 3,7-dimethyldecane (0.14 ± 0.15%), and octadecane (0.4 ± 0.53%). The metabolites that are shared with C. albicans, C. glabrata, C. parapsilosis, C. tropicalis and S. cerevisiae are phenylethanol, which is synthesized from phenylalanine, and eicosane and nonanoic acid (identified as trimethylsilyl ester), which are synthesized from fatty acid metabolism. Phenylethanol is the most abundant metabolite in all fungi tested: 26.36 ± 17.42% (C. glabrata), 46.77 ± 15.58% (C. albicans), 49.76 ± 18.43% (C. tropicalis), 5.72 ± 0.66% (C. parapsilosis.) and 44.58 ± 27.91% (S. cerevisiae). The analysis of C. glabrata's secretome will allow us to further our understanding of the possible role these metabolites could play in its virulence.

Broccoli: Combining Phylogenetic and Network Analyses for Orthology Assignment

Orthology assignment is a key step of comparative genomic studies, for which many bioinformatic tools have been developed. However, all gene clustering pipelines are based on the analysis of protein distances, which are subject to many artifacts. In this article, we introduce Broccoli, a user-friendly pipeline designed to infer, with high precision, orthologous groups, and pairs of proteins using a phylogeny-based approach. Briefly, Broccoli performs ultrafast phylogenetic analyses on most proteins and builds a network of orthologous relationships. Orthologous groups are then identified from the network using a parameter-free machine learning algorithm. Broccoli is also able to detect chimeric proteins resulting from gene-fusion events and to assign these proteins to the corresponding orthologous groups. Tested on two benchmark data sets, Broccoli outperforms current orthology pipelines. In addition, Broccoli is scalable, with runtimes similar to those of recent distance-based pipelines. Given its high level of performance and efficiency, this new pipeline represents a suitable choice for comparative genomic studies. Broccoli is freely available at https://github.com/rderelle/Broccoli.

Patterns and impacts of nonvertical evolution in eukaryotes: a paradigm shift

Evolution of eukaryotic species and their genomes has been traditionally understood as a vertical process in which genetic material is transmitted from parents to offspring along a lineage, and in which genetic exchange is restricted within species boundaries. However, mounting evidence from comparative genomics indicates that this paradigm is often violated. Horizontal gene transfer and mating between diverged lineages blur species boundaries and challenge the reconstruction of evolutionary histories of species and their genomes. Nonvertical evolution might be more restricted in eukaryotes than in prokaryotes, yet it is not negligible and can be common in certain groups. Recognition of such processes brings about the need to incorporate this complexity into our models, as well as to conceptually reframe eukaryotic diversity and evolution. Here, I review the recent work from genomics studies that supports the effects of nonvertical modes of evolution including introgression, hybridization, and horizontal gene transfer in different eukaryotic groups. I then discuss emerging patterns and effects, illustrated by specific examples, that support the conclusion that nonvertical processes are often at the root of important evolutionary transitions and adaptations. I will argue that a paradigm shift is needed to naturally accommodate nonvertical processes in eukaryotic evolution.

Biological and genomic analyses of a clinical isolate of Yarrowia galli from China

Infections caused by emerging fungal pathogens represent a new threat to human health. The yeast Yarrowia (Candida) galli was first described from chicken breast and liver in 2004 and has occasionally been isolated in clinical settings. In this study, we present the first report of a Y. galli isolate from a face granuloma of a woman. Y. galli is unable to grow at human physiological temperature (37 °C). Phenotypic analysis demonstrates that Y. galli can exist as several morphological types, namely fluffy, sticky, tight, and yeast forms, based on their cellular and colony appearances. Interestingly, Y. galli is able to undergo switching among different morphologies. These morphological changes are similar to the switching systems in pathogenic Candida species such as Candida albicans and Candida tropicalis. We further sequenced the genome of the Y. galli isolate. A comparative analysis with pathogenic yeast species indicated that a set of lipid metabolism genes were enriched in Y. galli. Domain enrichment analysis demonstrated that, similar to Candida clade species, the genome of Y. galli maintained several gene families required for virulence. Our biological and genomic analyses provide new insights into the understanding of the biology of Y. galli as either an environmental isolate or a potential human pathogen.

1,8-Substituted anthraquinones, anthrones and bianthrones as potential non-azole leads against fungal infections

The synthesis of a variety of 1,8-substituted anthraquinones, anthrones and bianthrones and their potential as antifungal agents is evaluated. Preliminary screening against Schizosaccharomyces pombe (S. pombe), a fission yeast, and Saccharomyces cerevisiae (S. cerevisiae), a budding yeast, is reported. Both these yeast species demonstrate close homologue to a number of pathogenic fungi.

The lichen symbiosis re-viewed through the genomes of Cladonia grayi and its algal partner Asterochloris glomerata

BACKGROUND

Lichens, encompassing 20,000 known species, are symbioses between specialized fungi (mycobionts), mostly ascomycetes, and unicellular green algae or cyanobacteria (photobionts). Here we describe the first parallel genomic analysis of the mycobiont Cladonia grayi and of its green algal photobiont Asterochloris glomerata. We focus on genes/predicted proteins of potential symbiotic significance, sought by surveying proteins differentially activated during early stages of mycobiont and photobiont interaction in coculture, expanded or contracted protein families, and proteins with differential rates of evolution.

RESULTS

A) In coculture, the fungus upregulated small secreted proteins, membrane transport proteins, signal transduction components, extracellular hydrolases and, notably, a ribitol transporter and an ammonium transporter, and the alga activated DNA metabolism, signal transduction, and expression of flagellar components. B) Expanded fungal protein families include heterokaryon incompatibility proteins, polyketide synthases, and a unique set of G-protein α subunit paralogs. Expanded algal protein families include carbohydrate active enzymes and a specific subclass of cytoplasmic carbonic anhydrases. The alga also appears to have acquired by horizontal gene transfer from prokaryotes novel archaeal ATPases and Desiccation-Related Proteins. Expanded in both symbionts are signal transduction components, ankyrin domain proteins and transcription factors involved in chromatin remodeling and stress responses. The fungal transportome is contracted, as are algal nitrate assimilation genes. C) In the mycobiont, slow-evolving proteins were enriched for components involved in protein translation, translocation and sorting.

CONCLUSIONS

The surveyed genes affect stress resistance, signaling, genome reprogramming, nutritional and structural interactions. The alga carries many genes likely transferred horizontally through viruses, yet we found no evidence of inter-symbiont gene transfer. The presence in the photobiont of meiosis-specific genes supports the notion that sexual reproduction occurs in Asterochloris while they are free-living, a phenomenon with implications for the adaptability of lichens and the persistent autonomy of the symbionts. The diversity of the genes affecting the symbiosis suggests that lichens evolved by accretion of many scattered regulatory and structural changes rather than through introduction of a few key innovations. This predicts that paths to lichenization were variable in different phyla, which is consistent with the emerging consensus that ascolichens could have had a few independent origins.

Genome-Wide Comparative Analysis of Phylogenetic Trees: The Prokaryotic Forest of Life

Genome-wide comparison of phylogenetic trees is becoming an increasingly common approach in evolutionary genomics, and a variety of approaches for such comparison have been developed. In this article we present several methods for comparative analysis of large numbers of phylogenetic trees. To compare phylogenetic trees taking into account the bootstrap support for each internal branch, the boot-split distance (BSD) method is introduced as an extension of the previously developed split distance (SD) method for tree comparison. The BSD method implements the straightforward idea that comparison of phylogenetic trees can be made more robust by treating tree splits differentially depending on the bootstrap support. Approaches are also introduced for detecting treelike and netlike evolutionary trends in the phylogenetic Forest of Life (FOL), i.e., the entirety of the phylogenetic trees for conserved genes of prokaryotes. The principal method employed for this purpose includes mapping quartets of species onto trees to calculate the support of each quartet topology and so to quantify the tree and net contributions to the distances between species. We describe the applications methods used to analyze the FOL and the results obtained with these methods. These results support the concept of the Tree of Life (TOL) as a central evolutionary trend in the FOL as opposed to the traditional view of the TOL as a "species tree."

The role of PCNA as a scaffold protein in cellular signaling is functionally conserved between yeast and humans

Proliferating cell nuclear antigen (PCNA), a member of the highly conserved DNA sliding clamp family, is an essential protein for cellular processes including DNA replication and repair. A large number of proteins from higher eukaryotes contain one of two PCNA-interacting motifs: PCNA-interacting protein box (PIP box) and AlkB homologue 2 PCNA-interacting motif (APIM). APIM has been shown to be especially important during cellular stress. PIP box is known to be functionally conserved in yeast, and here, we show that this is also the case for APIM. Several of the 84 APIM-containing yeast proteins are associated with cellular signaling as hub proteins, which are able to interact with a large number of other proteins. Cellular signaling is highly conserved throughout evolution, and we recently suggested a novel role for PCNA as a scaffold protein in cellular signaling in human cells. A cell-penetrating peptide containing the APIM sequence increases the sensitivity toward the chemotherapeutic agent cisplatin in both yeast and human cells, and both yeast and human cells become hypersensitive when the Hog1/p38 MAPK pathway is blocked. These results suggest that the interactions between APIM-containing signaling proteins and PCNA during the DNA damage response is evolutionary conserved between yeast and mammals and that PCNA has a role in cellular signaling also in yeast.

Adhesins in Candida glabrata

The human fungal pathogen Candida glabrata is causing more and more problems in hospitals, as this species shows an intrinsic antifungal drug resistance or rapidly becomes resistant when challenged with antifungals. C. glabrata only grows in the yeast form, so it is lacking a yeast-to-hyphae switch, which is one of the main virulence factors of C. albicans. An important virulence factor of C. glabrata is its capacity to strongly adhere to many different substrates. To achieve this, C. glabrata expresses a large number of adhesin-encoding genes and genome comparisons with closely related species, including the non-pathogenic S. cerevisiae, which revealed a correlation between the number of adhesin-encoding genes and pathogenicity. The adhesins are involved in the first steps during an infection; they are the first point of contact with the host. For several of these adhesins, their importance in adherence to different substrates and subsequent biofilm formation was demonstrated in vitro or in vivo. In this review, we provide an overview of the role of C. glabrata adhesins during adhesion and biofilm formation both, under in vitro and in vivo conditions.

Multiple Approaches to Phylogenomic Reconstruction of the Fungal Kingdom

Fungi are possibly the most diverse eukaryotic kingdom, with over a million member species and an evolutionary history dating back a billion years. Fungi have been at the forefront of eukaryotic genomics, and owing to initiatives like the 1000 Fungal Genomes Project the amount of fungal genomic data has increased considerably over the last 5 years, enabling large-scale comparative genomics of species across the kingdom. In this chapter, we first review fungal evolution and the history of fungal genomics. We then review in detail seven phylogenomic methods and reconstruct the phylogeny of 84 fungal species from 8 phyla using each method. Six methods have seen extensive use in previous fungal studies, while a Bayesian supertree method is novel to fungal phylogenomics. We find that both established and novel phylogenomic methods can accurately reconstruct the fungal kingdom. Finally, we discuss the accuracy and suitability of each phylogenomic method utilized.

Beyond the Whole-Genome Duplication: Phylogenetic Evidence for an Ancient Interspecies Hybridization in the Baker's Yeast Lineage

Whole-genome duplications have shaped the genomes of several vertebrate, plant, and fungal lineages. Earlier studies have focused on establishing when these events occurred and on elucidating their functional and evolutionary consequences, but we still lack sufficient understanding of how genome duplications first originated. We used phylogenomics to study the ancient genome duplication occurred in the yeast Saccharomyces cerevisiae lineage and found compelling evidence for the existence of a contemporaneous interspecies hybridization. We propose that the genome doubling was a direct consequence of this hybridization and that it served to provide stability to the recently formed allopolyploid. This scenario provides a mechanism for the origin of this ancient duplication and the lineage that originated from it and brings a new perspective to the interpretation of the origin and consequences of whole-genome duplications.

Evolutionary Advantage Conferred by an Eukaryote-to-Eukaryote Gene Transfer Event in Wine Yeasts

Although an increasing number of horizontal gene transfers have been reported in eukaryotes, experimental evidence for their adaptive value is lacking. Here, we report the recent transfer of a 158-kb genomic region between Torulaspora microellipsoides and Saccharomyces cerevisiae wine yeasts or closely related strains. This genomic region has undergone several rearrangements in S. cerevisiae strains, including gene loss and gene conversion between two tandemly duplicated FOT genes encoding oligopeptide transporters. We show that FOT genes confer a strong competitive advantage during grape must fermentation by increasing the number and diversity of oligopeptides that yeast can utilize as a source of nitrogen, thereby improving biomass formation, fermentation efficiency, and cell viability. Thus, the acquisition of FOT genes has favored yeast adaptation to the nitrogen-limited wine fermentation environment. This finding indicates that anthropic environments offer substantial ecological opportunity for evolutionary diversification through gene exchange between distant yeast species.

High variability of mitochondrial gene order among fungi

From their origin as an early alpha proteobacterial endosymbiont to their current state as cellular organelles, large-scale genomic reorganization has taken place in the mitochondria of all main eukaryotic lineages. So far, most studies have focused on plant and animal mitochondrial (mt) genomes (mtDNA), but fungi provide new opportunities to study highly differentiated mtDNAs. Here, we analyzed 38 complete fungal mt genomes to investigate the evolution of mtDNA gene order among fungi. In particular, we looked for evidence of nonhomologous intrachromosomal recombination and investigated the dynamics of gene rearrangements. We investigated the effect that introns, intronic open reading frames (ORFs), and repeats may have on gene order. Additionally, we asked whether the distribution of transfer RNAs (tRNAs) evolves independently to that of mt protein-coding genes. We found that fungal mt genomes display remarkable variation between and within the major fungal phyla in terms of gene order, genome size, composition of intergenic regions, and presence of repeats, introns, and associated ORFs. Our results support previous evidence for the presence of mt recombination in all fungal phyla, a process conspicuously lacking in most Metazoa. Overall, the patterns of rearrangements may be explained by the combined influences of recombination (i.e., most likely nonhomologous and intrachromosomal), accumulated repeats, especially at intergenic regions, and to a lesser extent, mobile element dynamics.

The genome of the generalist plant pathogen Fusarium avenaceum is enriched with genes involved in redox, signaling and secondary metabolism

Fusarium avenaceum is a fungus commonly isolated from soil and associated with a wide range of host plants. We present here three genome sequences of F. avenaceum, one isolated from barley in Finland and two from spring and winter wheat in Canada. The sizes of the three genomes range from 41.6–43.1 MB, with 13217–13445 predicted protein-coding genes. Whole-genome analysis showed that the three genomes are highly syntenic, and share>95% gene orthologs. Comparative analysis to other sequenced Fusaria shows that F. avenaceum has a very large potential for producing secondary metabolites, with between 75 and 80 key enzymes belonging to the polyketide, non-ribosomal peptide, terpene, alkaloid and indole-diterpene synthase classes. In addition to known metabolites from F. avenaceum, fuscofusarin and JM-47 were detected for the first time in this species. Many protein families are expanded in F. avenaceum, such as transcription factors, and proteins involved in redox reactions and signal transduction, suggesting evolutionary adaptation to a diverse and cosmopolitan ecology. We found that 20% of all predicted proteins were considered to be secreted, supporting a life in the extracellular space during interaction with plant hosts.

Complete DNA sequence of Kuraishia capsulata illustrates novel genomic features among budding yeasts (Saccharomycotina)

The numerous yeast genome sequences presently available provide a rich source of information for functional as well as evolutionary genomics but unequally cover the large phylogenetic diversity of extant yeasts. We present here the complete sequence of the nuclear genome of the haploid-type strain of Kuraishia capsulata (CBS1993^T), a nitrate-assimilating Saccharomycetales of uncertain taxonomy, isolated from tunnels of insect larvae underneath coniferous barks and characterized by its copious production of extracellular polysaccharides. The sequence is composed of seven scaffolds, one per chromosome, totaling 11.4 Mb and containing 6,029 protein-coding genes, ∼13.5% of which being interrupted by introns. This GC-rich yeast genome (45.7%) appears phylogenetically related with the few other nitrate-assimilating yeasts sequenced so far, Ogataea polymorpha, O. parapolymorpha, and Dekkera bruxellensis, with which it shares a very reduced number of tRNA genes, a novel tRNA sparing strategy, and a common nitrate assimilation cluster, three specific features to this group of yeasts. Centromeres were recognized in GC-poor troughs of each scaffold. The strain bears MAT alpha genes at a single MAT locus and presents a significant degree of conservation with Saccharomyces cerevisiae genes, suggesting that it can perform sexual cycles in nature, although genes involved in meiosis were not all recognized. The complete absence of conservation of synteny between K. capsulata and any other yeast genome described so far, including the three other nitrate-assimilating species, validates the interest of this species for long-range evolutionary genomic studies among Saccharomycotina yeasts.

Sporadic distribution of prion-forming ability of Sup35p from yeasts and fungi

Sup35p of Saccharomyces cerevisiae can form the [PSI+] prion, an infectious amyloid in which the protein is largely inactive. The part of Sup35p that forms the amyloid is the region normally involved in control of mRNA turnover. The formation of [PSI+] by Sup35p's from other yeasts has been interpreted to imply that the prion-forming ability of Sup35p is conserved in evolution, and thus of survival/fitness/evolutionary value to these organisms. We surveyed a larger number of yeast and fungal species by the same criteria as used previously and find that the Sup35p from many species cannot form prions. [PSI+] could be formed by the Sup35p from Candida albicans, Candida maltosa, Debaromyces hansenii, and Kluyveromyces lactis, but orders of magnitude less often than the S. cerevisiae Sup35p converts to the prion form. The Sup35s from Schizosaccharomyces pombe and Ashbya gossypii clearly do not form [PSI+]. We were also unable to detect [PSI+] formation by the Sup35ps from Aspergillus nidulans, Aspergillus fumigatus, Magnaporthe grisea, Ustilago maydis, or Cryptococcus neoformans. Each of two C. albicans SUP35 alleles can form [PSI+], but transmission from one to the other is partially blocked. These results suggest that the prion-forming ability of Sup35p is not a conserved trait, but is an occasional deleterious side effect of a protein domain conserved for another function.

Systematic phenotyping of a large-scale Candida glabrata deletion collection reveals novel antifungal tolerance genes

The opportunistic fungal pathogen Candida glabrata is a frequent cause of candidiasis, causing infections ranging from superficial to life-threatening disseminated disease. The inherent tolerance of C. glabrata to azole drugs makes this pathogen a serious clinical threat. To identify novel genes implicated in antifungal drug tolerance, we have constructed a large-scale C. glabrata deletion library consisting of 619 unique, individually bar-coded mutant strains, each lacking one specific gene, all together representing almost 12% of the genome. Functional analysis of this library in a series of phenotypic and fitness assays identified numerous genes required for growth of C. glabrata under normal or specific stress conditions, as well as a number of novel genes involved in tolerance to clinically important antifungal drugs such as azoles and echinocandins. We identified 38 deletion strains displaying strongly increased susceptibility to caspofungin, 28 of which encoding proteins that have not previously been linked to echinocandin tolerance. Our results demonstrate the potential of the C. glabrata mutant collection as a valuable resource in functional genomics studies of this important fungal pathogen of humans, and to facilitate the identification of putative novel antifungal drug target and virulence genes.

Clinical infections by the yeast-like pathogen Candida glabrata have been ever-increasing over the past years. Importantly, C. glabrata is one of the most prevalent causes of drug-refractory fungal infections in humans. We have generated a novel large-scale collection encompassing 619 bar-coded C. glabrata mutants, each lacking a single gene. Extensive profiling of phenotypes reveals a number of novel genes implicated in tolerance to antifungal drugs that interfere with proper cell wall function, as well as genes affecting fitness of C. glabrata both during normal growth and under environmental stress. This fungal deletion collection will be a valuable resource for the community to study mechanisms of virulence and antifungal drug tolerance in C. glabrata, which is particularly relevant in view of the increasing prevalence of infections caused by this important human fungal pathogen.

Independent expansion of zincin metalloproteinases in Onygenales fungi may be associated with their pathogenicity

To get a comprehensive view of fungal M35 family (deuterolysin) and M36 family (fungalysin) genes, we conducted genome-wide investigations and phylogenetic analyses of genes in these two families from 50 sequenced Ascomycota fungi with different life styles. Large variations in the number of M35 family and M36 family genes were found among different fungal genomes, indicating that these two gene families have been highly dynamic through fungal evolution. Moreover, we found obvious expansions of Meps in two families of Onygenales: Onygenaceae and Arthodermataceae, whereas species in family Ajellomycetace did not show expansion of these genes. The strikingly different gene duplication and loss patterns in Onygenales may be associated with the different pathogenicity of these species. Interestingly, likelihood ratio tests (LRT) of both M35 family and M36 family genes suggested that several branches leading to the duplicated genes in dermatophytic and Coccidioides fungi had signatures of positive selection, indicating that the duplicated Mep genes have likely diverged functionally to play important roles during the evolution of pathogenicity of dermatophytic and Coccidioides fungi. The potentially positively selected residues discovered by our analysis may have contributed to the development of new physiological functions of the duplicated Mep genes in dermatophytic fungi and Coccidioides species. Our study adds to the current knowledge of the evolution of Meps in fungi and also establishes a theoretical foundation for future experimental investigations.

A phylogenomics approach for selecting robust sets of phylogenetic markers

Reconstructing the evolutionary relationships of species is a major goal in biology. Despite the increasing number of completely sequenced genomes, a large number of phylogenetic projects rely on targeted sequencing and analysis of a relatively small sample of marker genes. The selection of these phylogenetic markers should ideally be based on accurate predictions of their combined, rather than individual, potential to accurately resolve the phylogeny of interest. Here we present and validate a new phylogenomics strategy to efficiently select a minimal set of stable markers able to reconstruct the underlying species phylogeny. In contrast to previous approaches, our methodology does not only rely on the ability of individual genes to reconstruct a known phylogeny, but it also explores the combined power of sets of concatenated genes to accurately infer phylogenetic relationships of species not previously analyzed. We applied our approach to two broad sets of cyanobacterial and ascomycetous fungal species, and provide two minimal sets of six and four genes, respectively, necessary to fully resolve the target phylogenies. This approach paves the way for the informed selection of phylogenetic markers in the effort of reconstructing the tree of life.

PhylomeDB v4: zooming into the plurality of evolutionary histories of a genome

Phylogenetic trees representing the evolutionary relationships of homologous genes are the entry point for many evolutionary analyses. For instance, the use of a phylogenetic tree can aid in the inference of orthology and paralogy relationships, and in the detection of relevant evolutionary events such as gene family expansions and contractions, horizontal gene transfer, recombination or incomplete lineage sorting. Similarly, given the plurality of evolutionary histories among genes encoded in a given genome, there is a need for the combined analysis of genome-wide collections of phylogenetic trees (phylomes). Here, we introduce a new release of PhylomeDB (http://phylomedb.org), a public repository of phylomes. Currently, PhylomeDB hosts 120 public phylomes, comprising >1.5 million maximum likelihood trees and multiple sequence alignments. In the current release, phylogenetic trees are annotated with taxonomic, protein-domain arrangement, functional and evolutionary information. PhylomeDB is also a major source for phylogeny-based predictions of orthology and paralogy, covering >10 million proteins across 1059 sequenced species. Here we describe newly implemented PhylomeDB features, and discuss a benchmark of the orthology predictions provided by the database, the impact of proteome updates and the use of the phylome approach in the analysis of newly sequenced genomes and transcriptomes.

Comparative genomics of emerging pathogens in the Candida glabrata clade

BACKGROUND

Candida glabrata follows C. albicans as the second or third most prevalent cause of candidemia worldwide. These two pathogenic yeasts are distantly related, C. glabrata being part of the Nakaseomyces, a group more closely related to Saccharomyces cerevisiae. Although C. glabrata was thought to be the only pathogenic Nakaseomyces, two new pathogens have recently been described within this group: C. nivariensis and C. bracarensis. To gain insight into the genomic changes underlying the emergence of virulence, we sequenced the genomes of these two, and three other non-pathogenic Nakaseomyces, and compared them to other sequenced yeasts.

RESULTS

Our results indicate that the two new pathogens are more closely related to the non-pathogenic N. delphensis than to C. glabrata. We uncover duplications and accelerated evolution that specifically affected genes in the lineage preceding the group containing N. delphensis and the three pathogens, which may provide clues to the higher propensity of this group to infect humans. Finally, the number of Epa-like adhesins is specifically enriched in the pathogens, particularly in C. glabrata.

CONCLUSIONS

Remarkably, some features thought to be the result of adaptation of C. glabrata to a pathogenic lifestyle, are present throughout the Nakaseomyces, indicating these are rather ancient adaptations to other environments. Phylogeny suggests that human pathogenesis evolved several times, independently within the clade. The expansion of the EPA gene family in pathogens establishes an evolutionary link between adhesion and virulence phenotypes. Our analyses thus shed light onto the relationships between virulence and the recent genomic changes that occurred within the Nakaseomyces.

SEQUENCE ACCESSION NUMBERS

Nakaseomyces delphensis: CAPT01000001 to CAPT01000179Candida bracarensis: CAPU01000001 to CAPU01000251Candida nivariensis: CAPV01000001 to CAPV01000123Candida castellii: CAPW01000001 to CAPW01000101Nakaseomyces bacillisporus: CAPX01000001 to CAPX01000186.

Measuring guide-tree dependency of inferred gaps in progressive aligners

MOTIVATION

Multiple sequence alignments are generally reconstructed using a progressive approach that follows a guide-tree. During this process, gaps are introduced at a cost to maximize residue pairing, but it is unclear whether inferred gaps reflect actual past events of sequence insertions or deletions. It has been found that patterns of inferred gaps in alignments contain information towards the true phylogeny, but it is as yet unknown whether gaps are simply reflecting information that was already present in the guide-tree.

RESULTS

We here develop a framework to disentangle the phylogenetic signal carried by gaps from that which is already present in the guide-tree. Our results indicate that most gaps are incorrectly inserted in patterns that, nevertheless, follow the guide-tree. Thus, most gap patterns in current alignments are not informative per se. This affects different programs to various degrees, PRANK being the most sensitive to the guide-tree.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Phylogenetic signal variation in the genomes of Medicago (Fabaceae)

Genome-scale data offer the opportunity to clarify phylogenetic relationships that are difficult to resolve with few loci, but they can also identify genomic regions with evolutionary history distinct from that of the species history. We collected whole-genome sequence data from 29 taxa in the legume genus Medicago, then aligned these sequences to the Medicago truncatula reference genome to confidently identify 87 596 variable homologous sites. We used this data set to estimate phylogenetic relationships among Medicago species, to investigate the number of sites needed to provide robust phylogenetic estimates and to identify specific genomic regions supporting topologies in conflict with the genome-wide phylogeny. Our full genomic data set resolves relationships within the genus that were previously intractable. Subsampling the data reveals considerable variation in phylogenetic signal and power in smaller subsets of the data. Even when sampling 5000 sites, no random sample of the data supports a topology identical to that of the genome-wide phylogeny. Phylogenetic relationships estimated from 500-site sliding windows revealed genome regions supporting several alternative species relationships among recently diverged taxa, consistent with the expected effects of deep coalescence or introgression in the recent history of Medicago.

Gene make-up: rapid and massive intron gains after horizontal transfer of a bacterial α-amylase gene to Basidiomycetes

BACKGROUND

Increasing genome data show that introns, a hallmark of eukaryotes, already existed at a high density in the last common ancestor of extant eukaryotes. However, intron content is highly variable among species. The tempo of intron gains and losses has been irregular and several factors may explain why some genomes are intron-poor whereas other are intron-rich.

RESULTS

We studied the dynamics of intron gains and losses in an α-amylase gene, whose product breaks down starch and other polysaccharides. It was transferred from an Actinobacterium to an ancestor of Agaricomycotina. This gene underwent further duplications in several species. The results indicate a high rate of intron insertions soon after the gene settled in the fungal genome. A number of these oldest introns, regularly scattered along the gene, remained conserved. Subsequent gains and losses were lineage dependent, with a majority of losses. Moreover, a few species exhibited a high number of both specific intron gains and losses in recent periods. There was little sequence conservation around insertion sites, then probably little information for splicing, whereas splicing sites, inside introns, showed typical and conserved patterns. There was little variation of intron size.

CONCLUSIONS

Since most Basidiomycetes have intron-rich genomes and this richness was ancestral in Fungi, long before the transfer event, we suggest that the new gene was shaped to comply with requirements of the splicing machinery, such as short exon and intron sizes, in order to be correctly processed.

Genome sequence of the necrotrophic fungus Penicillium digitatum, the main postharvest pathogen of citrus

BACKGROUND

Penicillium digitatum is a fungal necrotroph causing a common citrus postharvest disease known as green mold. In order to gain insight into the genetic bases of its virulence mechanisms and its high degree of host-specificity, the genomes of two P. digitatum strains that differ in their antifungal resistance traits have been sequenced and compared with those of 28 other Pezizomycotina.

RESULTS

The two sequenced genomes are highly similar, but important differences between them include the presence of a unique gene cluster in the resistant strain, and mutations previously shown to confer fungicide resistance. The two strains, which were isolated in Spain, and another isolated in China have identical mitochondrial genome sequences suggesting a recent worldwide expansion of the species. Comparison with the closely-related but non-phytopathogenic P. chrysogenum reveals a much smaller gene content in P. digitatum, consistent with a more specialized lifestyle. We show that large regions of the P. chrysogenum genome, including entire supercontigs, are absent from P. digitatum, and that this is the result of large gene family expansions rather than acquisition through horizontal gene transfer. Our analysis of the P. digitatum genome is indicative of heterothallic sexual reproduction and reveals the molecular basis for the inability of this species to assimilate nitrate or produce the metabolites patulin and penicillin. Finally, we identify the predicted secretome, which provides a first approximation to the protein repertoire used during invasive growth.

CONCLUSIONS

The complete genome of P. digitatum, the first of a phytopathogenic Penicillium species, is a valuable tool for understanding the virulence mechanisms and host-specificity of this economically important pest.

The Schistosoma mansoni phylome: using evolutionary genomics to gain insight into a parasite's biology

BACKGROUND

Schistosoma mansoni is one of the causative agents of schistosomiasis, a neglected tropical disease that affects about 237 million people worldwide. Despite recent efforts, we still lack a general understanding of the relevant host-parasite interactions, and the possible treatments are limited by the emergence of resistant strains and the absence of a vaccine. The S. mansoni genome was completely sequenced and still under continuous annotation. Nevertheless, more than 45% of the encoded proteins remain without experimental characterization or even functional prediction. To improve our knowledge regarding the biology of this parasite, we conducted a proteome-wide evolutionary analysis to provide a broad view of the S. mansoni's proteome evolution and to improve its functional annotation.

RESULTS

Using a phylogenomic approach, we reconstructed the S. mansoni phylome, which comprises the evolutionary histories of all parasite proteins and their homologs across 12 other organisms. The analysis of a total of 7,964 phylogenies allowed a deeper understanding of genomic complexity and evolutionary adaptations to a parasitic lifestyle. In particular, the identification of lineage-specific gene duplications pointed to the diversification of several protein families that are relevant for host-parasite interaction, including proteases, tetraspanins, fucosyltransferases, venom allergen-like proteins, and tegumental-allergen-like proteins. In addition to the evolutionary knowledge, the phylome data enabled us to automatically re-annotate 3,451 proteins through a phylogenetic-based approach rather than solely sequence similarity searches. To allow further exploitation of this valuable data, all information has been made available at PhylomeDB (http://www.phylomedb.org).

CONCLUSIONS

In this study, we used an evolutionary approach to assess S. mansoni parasite biology, improve genome/proteome functional annotation, and provide insights into host-parasite interactions. Taking advantage of a proteome-wide perspective rather than focusing on individual proteins, we identified that this parasite has experienced specific gene duplication events, particularly affecting genes that are potentially related to the parasitic lifestyle. These innovations may be related to the mechanisms that protect S. mansoni against host immune responses being important adaptations for the parasite survival in a potentially hostile environment. Continuing this work, a comparative analysis involving genomic, transcriptomic, and proteomic data from other helminth parasites, other parasites, and vectors will supply more information regarding parasite's biology as well as host-parasite interactions.

Candida species: new insights into biofilm formation

Biofilms of Candida albicans, Candida parapsilosis, Candida glabrata and Candida tropicalis are associated with high indices of hospital morbidity and mortality. Major factors involved in the formation and growth of Candida biofilms are the chemical composition of the medical implant and the cell wall adhesins responsible for mediating Candida-Candida, Candida-human host cell and Candida-medical device adhesion. Strategies for elucidating the mechanisms that regulate the formation of Candida biofilms combine tools from biology, chemistry, nanoscience, material science and physics. This review proposes the use of new technologies, such as synchrotron radiation, to study the mechanisms of biofilm formation. In the future, this information is expected to facilitate the design of new materials and antifungal compounds that can eradicate nosocomial Candida infections due to biofilm formation on medical implants. This will reduce dissemination of candidiasis and hopefully improve the quality of life of patients.

Systems biology of host-fungus interactions: turning complexity into simplicity

Modeling interactions between fungi and their hosts at the systems level requires a molecular understanding both of how the host orchestrates immune surveillance and tolerance, and how this activation, in turn, affects fungal adaptation and survival. The transition from the commensal to pathogenic state, and the co-evolution of fungal strains within their hosts, necessitates the molecular dissection of fungal traits responsible for these interactions. There has been a dramatic increase in publically available genome-wide resources addressing fungal pathophysiology and host–fungal immunology. The integration of these existing data and emerging large-scale technologies addressing host–pathogen interactions requires novel tools to connect genome-wide data sets and theoretical approaches with experimental validation so as to identify inherent and emerging properties of host–pathogen relationships and to obtain a holistic view of infectious processes. If successful, a better understanding of the immune response in health and microbial diseases will eventually emerge and pave the way for improved therapies.

Phylogenomics supports microsporidia as the earliest diverging clade of sequenced fungi

BACKGROUND

Microsporidia is one of the taxa that have experienced the most dramatic taxonomic reclassifications. Once thought to be among the earliest diverging eukaryotes, the fungal nature of this group of intracellular pathogens is now widely accepted. However, the specific position of microsporidia within the fungal tree of life is still debated. Due to the presence of accelerated evolutionary rates, phylogenetic analyses involving microsporidia are prone to methodological artifacts, such as long-branch attraction, especially when taxon sampling is limited.

RESULTS

Here we exploit the recent availability of six complete microsporidian genomes to re-assess the long-standing question of their phylogenetic position. We show that microsporidians have a similar low level of conservation of gene neighborhood with other groups of fungi when controlling for the confounding effects of recent segmental duplications. A combined analysis of thousands of gene trees supports a topology in which microsporidia is a sister group to all other sequenced fungi. Moreover, this topology received increased support when less informative trees were discarded. This position of microsporidia was also strongly supported based on the combined analysis of 53 concatenated genes, and was robust to filters controlling for rate heterogeneity, compositional bias, long branch attraction and heterotachy.

CONCLUSIONS

Altogether, our data strongly support a scenario in which microsporidia is the earliest-diverging clade of sequenced fungi.

Global analysis of the evolution and mechanism of echinocandin resistance in Candida glabrata

The evolution of drug resistance has a profound impact on human health. Candida glabrata is a leading human fungal pathogen that can rapidly evolve resistance to echinocandins, which target cell wall biosynthesis and are front-line therapeutics for Candida infections. Here, we provide the first global analysis of mutations accompanying the evolution of fungal drug resistance in a human host utilizing a series of C. glabrata isolates that evolved echinocandin resistance in a patient treated with the echinocandin caspofungin for recurring bloodstream candidemia. Whole genome sequencing identified a mutation in the drug target, FKS2, accompanying a major resistance increase, and 8 additional non-synonymous mutations. The FKS2-T1987C mutation was sufficient for echinocandin resistance, and associated with a fitness cost that was mitigated with further evolution, observed in vitro and in a murine model of systemic candidemia. A CDC6-A511G(K171E) mutation acquired before FKS2-T1987C(S663P), conferred a small resistance increase. Elevated dosage of CDC55, which acquired a C463T(P155S) mutation after FKS2-T1987C(S663P), ameliorated fitness. To discover strategies to abrogate echinocandin resistance, we focused on the molecular chaperone Hsp90 and downstream effector calcineurin. Genetic or pharmacological compromise of Hsp90 or calcineurin function reduced basal tolerance and resistance. Hsp90 and calcineurin were required for caspofungin-dependent FKS2 induction, providing a mechanism governing echinocandin resistance. A mitochondrial respiration-defective petite mutant in the series revealed that the petite phenotype does not confer echinocandin resistance, but renders strains refractory to synergy between echinocandins and Hsp90 or calcineurin inhibitors. The kidneys of mice infected with the petite mutant were sterile, while those infected with the HSP90-repressible strain had reduced fungal burden. We provide the first global view of mutations accompanying the evolution of fungal drug resistance in a human host, implicate the premier compensatory mutation mitigating the cost of echinocandin resistance, and suggest a new mechanism of echinocandin resistance with broad therapeutic potential.

The evolution of drug resistance poses a severe threat to human health. Candida glabrata is a leading cause of mortality due to fungal infections worldwide. It can rapidly evolve resistance to drugs such as echinocandins, which target the fungal cell wall and are front-line therapeutics for Candida infections. We harness whole genome sequencing to provide a global view of mutations that accumulate in C. glabrata during the evolution of echinocandin resistance in a human host. Nine non-synonymous mutations were identified, including one in the echinocandin target. A mutation in an additional gene conferred a small resistance increase and another was in a gene whose dosage mitigated the fitness cost of resistance. We further discovered that compromising function of the molecular chaperone Hsp90 abrogates drug resistance and reduces kidney fungal burden in a mouse model of infection. Hsp90 and its downstream effector calcineurin are required for induction of the drug target in response to drug. Thus, we reveal the first global portrait of antifungal resistance mutations that evolve in a human host, identify the first compensatory mutation that mitigates the cost of echinocandin resistance, and suggest a new mechanism of echinocandin resistance that can be exploited to treat life-threatening fungal infections.

Genome-wide comparative analysis of phylogenetic trees: the prokaryotic forest of life

Genome-wide comparison of phylogenetic trees is becoming an increasingly common approach in evolutionary genomics, and a variety of approaches for such comparison have been developed. In this article, we present several methods for comparative analysis of large numbers of phylogenetic trees. To compare phylogenetic trees taking into account the bootstrap support for each internal branch, the Boot-Split Distance (BSD) method is introduced as an extension of the previously developed Split Distance method for tree comparison. The BSD method implements the straightforward idea that comparison of phylogenetic trees can be made more robust by treating tree splits differentially depending on the bootstrap support. Approaches are also introduced for detecting tree-like and net-like evolutionary trends in the phylogenetic Forest of Life (FOL), i.e., the entirety of the phylogenetic trees for conserved genes of prokaryotes. The principal method employed for this purpose includes mapping quartets of species onto trees to calculate the support of each quartet topology and so to quantify the tree and net contributions to the distances between species. We describe the application of these methods to analyze the FOL and the results obtained with these methods. These results support the concept of the Tree of Life (TOL) as a central evolutionary trend in the FOL as opposed to the traditional view of the TOL as a "species tree."

A consistent phylogenetic backbone for the fungi

The kingdom of fungi provides model organisms for biotechnology, cell biology, genetics, and life sciences in general. Only when their phylogenetic relationships are stably resolved, can individual results from fungal research be integrated into a holistic picture of biology. However, and despite recent progress, many deep relationships within the fungi remain unclear. Here, we present the first phylogenomic study of an entire eukaryotic kingdom that uses a consistency criterion to strengthen phylogenetic conclusions. We reason that branches (splits) recovered with independent data and different tree reconstruction methods are likely to reflect true evolutionary relationships. Two complementary phylogenomic data sets based on 99 fungal genomes and 109 fungal expressed sequence tag (EST) sets analyzed with four different tree reconstruction methods shed light from different angles on the fungal tree of life. Eleven additional data sets address specifically the phylogenetic position of Blastocladiomycota, Ustilaginomycotina, and Dothideomycetes, respectively. The combined evidence from the resulting trees supports the deep-level stability of the fungal groups toward a comprehensive natural system of the fungi. In addition, our analysis reveals methodologically interesting aspects. Enrichment for EST encoded data—a common practice in phylogenomic analyses—introduces a strong bias toward slowly evolving and functionally correlated genes. Consequently, the generalization of phylogenomic data sets as collections of randomly selected genes cannot be taken for granted. A thorough characterization of the data to assess possible influences on the tree reconstruction should therefore become a standard in phylogenomic analyses.

Reconstructing the fungal tree of life using phylogenomics and a preliminary investigation of the distribution of yeast prion-like proteins in the fungal kingdom

We have used three independent phylogenomic approaches (concatenated alignments, single-, and multi-gene supertrees) to reconstruct the fungal tree of life (FTOL) using publicly available fungal genomes. This is the first time multi-gene families have been used in fungal supertree reconstruction and permits us to use up to 66% of the 1,001,217 genes in our fungal database. Our analyses show that different phylogenomic datasets derived from varying clustering criteria and alignment orientation do not have a major effect on phylogenomic supertree reconstruction. Overall the resultant phylogenomic trees are relatively congruent with one another and successfully recover the major fungal phyla, subphyla and classes. We find that where incongruences do occur, the inferences are usually poorly supported. Within the Ascomycota phylum, our phylogenies reconstruct monophyletic Saccharomycotina and Pezizomycotina subphyla clades and infer a sister group relationship between these to the exclusion of the Taphrinomycotina. Within the Pezizomycotina subphylum, all three phylogenies infer a sister group relationship between the Leotiomycetes and Sordariomycetes classes. However, there is conflict regarding the relationships with the Dothideomycetes and Eurotiomycetes classes. Within the Basidiomycota phylum, supertrees derived from single- and multi-gene families infer a sister group relationship between the Pucciniomycotina and Agaricomycotina subphyla while the concatenated phylogeny infers a poorly supported relationship between the Agaricomycotina and Ustilagomycotina. The reconstruction of a robust FTOL is important for future fungal comparative analyses. We illustrate this point by performing a preliminary investigation into the phyletic distribution of yeast prion-like proteins in the fungal kingdom.

Addressing inter-gene heterogeneity in maximum likelihood phylogenomic analysis: yeasts revisited

Phylogenomic approaches to the resolution of inter-species relationships have become well established in recent years. Often these involve concatenation of many orthologous genes found in the respective genomes followed by analysis using standard phylogenetic models. Genome-scale data promise increased resolution by minimising sampling error, yet are associated with well-known but often inappropriately addressed caveats arising through data heterogeneity and model violation. These can lead to the reconstruction of highly-supported but incorrect topologies. With the aim of obtaining a species tree for 18 species within the ascomycetous yeasts, we have investigated the use of appropriate evolutionary models to address inter-gene heterogeneities and the scalability and validity of supermatrix analysis as the phylogenetic problem becomes more difficult and the number of genes analysed approaches truly phylogenomic dimensions. We have extended a widely-known early phylogenomic study of yeasts by adding additional species to increase diversity and augmenting the number of genes under analysis. We have investigated sophisticated maximum likelihood analyses, considering not only a concatenated version of the data but also partitioned models where each gene constitutes a partition and parameters are free to vary between the different partitions (thereby accounting for variation in the evolutionary processes at different loci). We find considerable increases in likelihood using these complex models, arguing for the need for appropriate models when analyzing phylogenomic data. Using these methods, we were able to reconstruct a well-supported tree for 18 ascomycetous yeasts spanning about 250 million years of evolution.

SUS1 introns are required for efficient mRNA nuclear export in yeast

Efficient coupling between mRNA synthesis and export is essential for gene expression. Sus1/ENY2, a component of the SAGA and TREX-2 complexes, is involved in both transcription and mRNA export. While most yeast genes lack introns, we previously reported that yeast SUS1 bears two. Here we show that this feature is evolutionarily conserved and critical for Sus1 function. We determine that while SUS1 splicing is inefficient, it responds to cellular conditions, and intronic mutations either promoting or blocking splicing lead to defects in mRNA export and cell growth. Consistent with this, we find that an intron-less SUS1 only partially rescues sus1Δ phenotypes. Remarkably, splicing of each SUS1 intron is also affected by the presence of the other and by SUS1 exonic sequences. Moreover, by following SUS1 RNA and protein levels we establish that nonsense-mediated decay (NMD) pathway and the splicing factor Mud2 both play a role in SUS1 expression. Our data (and those of the accompanying work by Hossain et al.) provide evidence of the involvement of splicing, translation, and decay in the regulation of early events in mRNP biogenesis; and imply the additional requirement for a balance in splicing isoforms from a single gene.

Finding single copy genes out of sequenced genomes for multilocus phylogenetics in non-model fungi

Historically, fungal multigene phylogenies have been reconstructed based on a small number of commonly used genes. The availability of complete fungal genomes has given rise to a new wave of model organisms that provide large number of genes potentially useful for building robust gene genealogies. Unfortunately, cross-utilization of these resources to study phylogenetic relationships in the vast majority of non-model fungi (i.e. “orphan” species) remains an unexamined question. To address this problem, we developed a method coupled with a program named “PHYLORPH” (PHYLogenetic markers for ORPHans). The method screens fungal genomic databases (107 fungal genomes fully sequenced) for single copy genes that might be easily transferable and well suited for studies at low taxonomic levels (for example, in species complexes) in non-model fungal species. To maximize the chance to target genes with informative regions, PHYLORPH displays a graphical evaluation system based on the estimation of nucleotide divergence relative to substitution type. The usefulness of this approach was tested by developing markers in four non-model groups of fungal pathogens. For each pathogen considered, 7 to 40% of the 10–15 best candidate genes proposed by PHYLORPH yielded sequencing success. Levels of polymorphism of these genes were compared with those obtained for some genes traditionally used to build fungal phylogenies (e.g. nuclear rDNA, β-tubulin, γ-actin, Elongation factor EF-1α). These genes were ranked among the best-performing ones and resolved accurately taxa relationships in each of the four non-model groups of fungi considered. We envision that PHYLORPH will constitute a useful tool for obtaining new and accurate phylogenetic markers to resolve relationships between closely related non-model fungal species.

Molecular evolution of urea amidolyase and urea carboxylase in fungi

Background

Urea amidolyase breaks down urea into ammonia and carbon dioxide in a two-step process, while another enzyme, urease, does this in a one step-process. Urea amidolyase has been found only in some fungal species among eukaryotes. It contains two major domains: the amidase and urea carboxylase domains. A shorter form of urea amidolyase is known as urea carboxylase and has no amidase domain. Eukaryotic urea carboxylase has been found only in several fungal species and green algae. In order to elucidate the evolutionary origin of urea amidolyase and urea carboxylase, we studied the distribution of urea amidolyase, urea carboxylase, as well as other proteins including urease, across kingdoms.

Results

Among the 64 fungal species we examined, only those in two Ascomycota classes (Sordariomycetes and Saccharomycetes) had the urea amidolyase sequences. Urea carboxylase was found in many but not all of the species in the phylum Basidiomycota and in the subphylum Pezizomycotina (phylum Ascomycota). It was completely absent from the class Saccharomycetes (phylum Ascomycota; subphylum Saccharomycotina). Four Sordariomycetes species we examined had both the urea carboxylase and the urea amidolyase sequences. Phylogenetic analysis showed that these two enzymes appeared to have gone through independent evolution since their bacterial origin. The amidase domain and the urea carboxylase domain sequences from fungal urea amidolyases clustered strongly together with the amidase and urea carboxylase sequences, respectively, from a small number of beta- and gammaproteobacteria. On the other hand, fungal urea carboxylase proteins clustered together with another copy of urea carboxylases distributed broadly among bacteria. The urease proteins were found in all the fungal species examined except for those of the subphylum Saccharomycotina.

Conclusions

We conclude that the urea amidolyase genes currently found only in fungi are the results of a horizontal gene transfer event from beta-, gamma-, or related species of proteobacteria. The event took place before the divergence of the subphyla Pezizomycotina and Saccharomycotina but after the divergence of the subphylum Taphrinomycotina. Urea carboxylase genes currently found in fungi and other limited organisms were also likely derived from another ancestral gene in bacteria. Our study presented another important example showing plastic and opportunistic genome evolution in bacteria and fungi and their evolutionary interplay.

TreeKO: a duplication-aware algorithm for the comparison of phylogenetic trees

Comparisons of tree topologies provide relevant information in evolutionary studies. Most existing methods share the drawback of requiring a complete and exact mapping of terminal nodes between the compared trees. This severely limits the scope of genome-wide analyses, since trees containing duplications are pruned arbitrarily or discarded. To overcome this, we have developed treeKO, an algorithm that enables the comparison of tree topologies, even in the presence of duplication and loss events. To do so treeKO recursively splits gene trees into pruned trees containing only orthologs to subsequently compute a distance based on the combined analyses of all pruned tree comparisons. In addition treeKO, implements the possibility of computing phylome support values, and reconciliation-based measures such as the number of inferred duplication and loss events.

Hal: an automated pipeline for phylogenetic analyses of genomic data

The rapid increase in genomic and genome-scale data is resulting in unprecedented levels of discrete sequence data available for phylogenetic analyses. Major analytical impasses exist, however, prior to analyzing these data with existing phylogenetic software. Obstacles include the management of large data sets without standardized naming conventions, identification and filtering of orthologous clusters of proteins or genes, and the assembly of alignments of orthologous sequence data into individual and concatenated super alignments. Here we report the production of an automated pipeline, Hal that produces multiple alignments and trees from genomic data. These alignments can be produced by a choice of four alignment programs and analyzed by a variety of phylogenetic programs. In short, the Hal pipeline connects the programs BLASTP, MCL, user specified alignment programs, GBlocks, ProtTest and user specified phylogenetic programs to produce species trees. The script is available at sourceforge (http://sourceforge.net/projects/bio-hal/). The results from an example analysis of Kingdom Fungi are briefly discussed.

Regulation of Candida glabrata oxidative stress resistance is adapted to host environment

The human fungal pathogen Candida glabrata is related to Saccharomyces cerevisiae but has developed high resistance against reactive oxygen species. We find that induction of conserved genes encoding antioxidant functions is dependent on the transcription factors CgYap1 and CgSkn7 which cooperate for promoter recognition. Superoxide stress resistance of C. glabrata is provided by superoxide dismutase CgSod1, which is not dependent on CgYap1/Skn7. Only double mutants lacking both CgSod1 and CgYap1 were efficiently killed by primary mouse macrophages. Our results suggest that in C. glabrata the regulation of key genes providing stress protection is adopted to meet a host-pathogen situation.

MetaPhOrs: orthology and paralogy predictions from multiple phylogenetic evidence using a consistency-based confidence score

Reliable prediction of orthology is central to comparative genomics. Approaches based on phylogenetic analyses closely resemble the original definition of orthology and paralogy and are known to be highly accurate. However, the large computational cost associated to these analyses is a limiting factor that often prevents its use at genomic scales. Recently, several projects have addressed the reconstruction of large collections of high-quality phylogenetic trees from which orthology and paralogy relationships can be inferred. This provides us with the opportunity to infer the evolutionary relationships of genes from multiple, independent, phylogenetic trees. Using such strategy, we combine phylogenetic information derived from different databases, to predict orthology and paralogy relationships for 4.1 million proteins in 829 fully sequenced genomes. We show that the number of independent sources from which a prediction is made, as well as the level of consistency across predictions, can be used as reliable confidence scores. A webserver has been developed to easily access these data (http://orthology.phylomedb.org), which provides users with a global repository of phylogeny-based orthology and paralogy predictions.