Jason C. Slot

Reconstructing the early evolution of Fungi using a six-gene phylogeny

The ancestors of fungi are believed to be simple aquatic forms with flagellated spores, similar to members of the extant phylum Chytridiomycota (chytrids). Current classifications assume that chytrids form an early-diverging clade within the kingdom Fungi and imply a single loss of the spore flagellum, leading to the diversification of terrestrial fungi. Here we develop phylogenetic hypotheses for Fungi using data from six gene regions and nearly 200 species. Our results indicate that there may have been at least four independent losses of the flagellum in the kingdom Fungi. These losses of swimming spores coincided with the evolution of new mechanisms of spore dispersal, such as aerial dispersal in mycelial groups and polar tube eversion in the microsporidia (unicellular forms that lack mitochondria). The enigmatic microsporidia seem to be derived from an endoparasitic chytrid ancestor similar to Rozella allomycis, on the earliest diverging branch of the fungal phylogenetic tree.

The Paleozoic Origin of Enzymatic Lignin Decomposition Reconstructed from 31 Fungal Genomes

Dating Wood Rot Specific lineages within the basidiomycete fungi, white rot species, have evolved the ability to break up a major structural component of woody plants, lignin, relative to their non–lignin-decaying brown rot relatives. Through the deep phylogenetic sampling of fungal genomes, Floudas et al. (p. 1715; see the Perspective by Hittinger) mapped the detailed evolution of wood-degrading enzymes. A key peroxidase and other enzymes involved in lignin decay were present in the common ancestor of the Agaricomycetes. These genes then expanded through gene duplications in parallel, giving rise to white rot lineages. The enzyme family that enables fungi to digest lignin expanded around the end of the coal-forming Carboniferous period. Wood is a major pool of organic carbon that is highly resistant to decay, owing largely to the presence of lignin. The only organisms capable of substantial lignin decay are white rot fungi in the Agaricomycetes, which also contains non–lignin-degrading brown rot and ectomycorrhizal species. Comparative analyses of 31 fungal genomes (12 generated for this study) suggest that lignin-degrading peroxidases expanded in the lineage leading to the ancestor of the Agaricomycetes, which is reconstructed as a white rot species, and then contracted in parallel lineages leading to brown rot and mycorrhizal species. Molecular clock analyses suggest that the origin of lignin degradation might have coincided with the sharp decrease in the rate of organic carbon burial around the end of the Carboniferous period.

Multiple GAL pathway gene clusters evolved independently and by different mechanisms in fungi

A notable characteristic of fungal genomes is that genes involved in successive steps of a metabolic pathway are often physically linked or clustered. To investigate how such clusters of functionally related genes are assembled and maintained, we examined the evolution of gene sequences and order in the galactose utilization (GAL) pathway in whole-genome data from 80 diverse fungi. We found that GAL gene clusters originated independently and by different mechanisms in three unrelated yeast lineages. Specifically, the GAL cluster found in Saccharomyces and Candida yeasts originated through the relocation of native unclustered genes, whereas the GAL cluster of Schizosaccharomyces yeasts was acquired through horizontal gene transfer from a Candida yeast. In contrast, the GAL cluster of Cryptococcus yeasts was assembled independently from the Saccharomyces/Candida and Schizosaccharomyces GAL clusters and coexists in the Cryptococcus genome with unclustered GAL paralogs. These independently evolved GAL clusters represent a striking example of analogy at the genomic level. We also found that species with GAL clusters exhibited significantly higher rates of GAL pathway loss than species with unclustered GAL genes. These results suggest that clustering of metabolic genes might facilitate fungal adaptation to changing environments both through the acquisition and loss of metabolic capacities.

Comparative genomics of biotechnologically important yeasts

Significance The highly diverse Ascomycete yeasts have enormous biotechnological potential. Collectively, these yeasts convert a broad range of substrates into useful compounds, such as ethanol, lipids, and vitamins, and can grow in extremes of temperature, salinity, and pH. We compared 29 yeast genomes with the goal of correlating genetics to useful traits. In one rare species, we discovered a genetic code that translates CUG codons to alanine rather than canonical leucine. Genome comparison enabled correlation of genes to useful metabolic properties and showed the synteny of the mating-type locus to be conserved over a billion years of evolution. Our study provides a roadmap for future biotechnological exploitations. 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.

Horizontal Transfer and Death of a Fungal Secondary Metabolic Gene Cluster

A cluster composed of four structural and two regulatory genes found in several species of the fungal genus Fusarium (class Sordariomycetes) is responsible for the production of the red pigment bikaverin. We discovered that the unrelated fungus Botrytis cinerea (class Leotiomycetes) contains a cluster of five genes that is highly similar in sequence and gene order to the Fusarium bikaverin cluster. Synteny conservation, nucleotide composition, and phylogenetic analyses of the cluster genes indicate that the B. cinerea cluster was acquired via horizontal transfer from a Fusarium donor. Upon or subsequent to the transfer, the B. cinerea gene cluster became inactivated; one of the four structural genes is missing, two others are pseudogenes, and the fourth structural gene shows an accelerated rate of nonsynonymous substitutions along the B. cinerea lineage, consistent with relaxation of selective constraints. Interestingly, the bik4 regulatory gene is still intact and presumably functional, whereas bik5, which is a pathway-specific regulator, also shows a mild but significant acceleration of evolutionary rate along the B. cinerea lineage. This selective preservation of the bik4 regulator suggests that its conservation is due to its likely involvement in other non–bikaverin-related biological processes in B. cinerea. Thus, in addition to novel metabolism, horizontal transfer of wholesale metabolic gene clusters might also be contributing novel regulation.

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.

Physical linkage of metabolic genes in fungi is an adaptation against the accumulation of toxic intermediate compounds

Genomic analyses have proliferated without being tied to tangible phenotypes. For example, although coordination of both gene expression and genetic linkage have been offered as genetic mechanisms for the frequently observed clustering of genes participating in fungal metabolic pathways, elucidation of the phenotype(s) favored by selection, resulting in cluster formation and maintenance, has not been forthcoming. We noted that the cause of certain well-studied human metabolic disorders is the accumulation of toxic intermediate compounds (ICs), which occurs when the product of an enzyme is not used as a substrate by a downstream neighbor in the metabolic network. This raises the hypothesis that the phenotype favored by selection to drive gene clustering is the mitigation of IC toxicity. To test this, we examined 100 diverse fungal genomes for the simplest type of cluster, gene pairs that are both metabolic neighbors and chromosomal neighbors immediately adjacent to each other, which we refer to as “double neighbor gene pairs” (DNGPs). Examination of the toxicity of their corresponding ICs shows that, compared with chromosomally nonadjacent metabolic neighbors, DNGPs are enriched for ICs that have acutely toxic LD50 doses or reactive functional groups. Furthermore, DNGPs are significantly more likely to be divergently oriented on the chromosome; remarkably, ∼40% of these DNGPs have ICs known to be toxic. We submit that the structure of synteny in metabolic pathways of fungi is a signature of selection for protection against the accumulation of toxic metabolic intermediates.

In Silico Characterization and Molecular Evolutionary Analysis of a Novel Superfamily of Fungal Effector Proteins

Most fungal plant pathogens secrete effector proteins during pathogenesis to manipulate their hosts defense and promote disease. These are so highly diverse in sequence and distribution, they are essentially considered as species-specific. However, we have recently shown the presence of homologous effectors in fungal species of the Dothideomycetes class. One such example is Ecp2, an effector originally described in the tomato pathogen Cladosporium fulvum but later detected in the plant pathogenic fungi Mycosphaerella fijiensis and Mycosphaerella graminicola as well. Here, using in silico sequence-similarity searches against a database of 135 fungal genomes and GenBank, we extend our queries for homologs of Ecp2 to the fungal kingdom and beyond, and further study their history of diversification. Our analyses show that Ecp2 homologs are members of an ancient and widely distributed superfamily of putative fungal effectors, which we term Hce2 for Homologs of C. fulvum Ecp2. Molecular evolutionary analyses show that the superfamily originated and diversified within the fungal kingdom, experiencing multiple lineage-specific expansions and losses that are consistent with the birth-and-death model of gene family evolution. Newly formed paralogs appear to be subject to diversification early after gene duplication events, whereas at later stages purifying selection acts to preserve diversity and the newly evolved putative functions. Some members of the Hce2 superfamily are fused to fungal Glycoside Hydrolase family 18 chitinases that show high similarity to the Zymocin killer toxin from the dairy yeast Kluyveromyces lactis, suggesting an analogous role in antagonistic interactions. The observed high rates of gene duplication and loss in the Hce2 superfamily, combined with diversification in both sequence and possibly functions within and between species, suggest that Hce2s are involved in adaptation to stresses and new ecological niches. Such findings address the need to rationalize effector biology and evolution beyond the perspective of solely host-microbe interactions.

Dimensions of Horizontal Gene Transfer in Eukaryotic Microbial Pathogens

Comparative genomic studies of microorganisms have disrupted the paradigm of vertical inheritance with modification. First in bacteria, and more recently in microscopic and even multicellular eukaryotes, horizontal gene transfer (HGT) has been implicated in genomic and ecological evolution. HGT is the exchange of genetic material between organisms that occurs independently of meiotic and mitotic recombination between mating or hybridizing individuals. HGT occurs as viral and plasmid-mediated transfer, and transformation by environmental DNA via known or yet-unknown mechanisms [1]. The existence of environmental gene pools and pan-genomes is supported by decades of functional and phylogenetic studies in bacteria that have highlighted the exchange and proliferation of virulence factors and antibiotic resistance mechanisms [2–4]. Presently, accumulating reports of HGT in eukaryotes raise similar questions of how exposure to such gene pools has impacted the evolution of eukaryotic microbes, and whether or not human activities influence HGT dynamics. Here, we describe the evidence supporting HGT in eukaryotic microbial pathogens from divergent lineages that impact human, animal, and plant health (S1 Table). We consider three interacting dimensions affecting the prevalence of HGT (genetic network structure, selectable functions, and opportunity for contact) in order to better understand how HGT manifests in this important group of organisms.

Repeated loss of an anciently horizontally transferred gene cluster in Botrytis

At least five of the six genes of the bikaverin secondary metabolic gene cluster were shown to have undergone horizontal transfer (HGT) from a Fusarium donor to the Botrytis lineage. Of these five, two enzyme-encoding genes are found as pseudogenes in B. cinerea whereas two regulatory genes and the transporter remain intact. To reconstruct the evolutionary events leading to decay of this gene cluster and infer a more precise timing of its transfer, we examined the genomes of nine additional broadly sampled Botrytis species. We found evidence that a Botrytis ancestor acquired the entire gene cluster through an ancient HGT that occurred before the diversification of the genus. During the subsequent evolution and diversification of the genus, four of the 10 genomes appear to have lost the gene cluster, while in the other six the cluster is in various stages of degeneration. Across the Botrytis genomes, the modes of gene decay in the cluster differed between enzyme-encoding genes, which had higher rates of transition to or retention of pseudogenes and were universally inactivated, and regulatory genes (particularly the non-pathway-specific regulator bik4), which more frequently appeared intact. Consistent with these results, the regulatory genes bik4 and bik5 showed stronger evidence of transcriptional expression than other bikaverin genes under multiple conditions in B. cinerea. These results could be explained by pleiotropy in the bikaverin regulatory genes either through rewiring or their interaction with more central pathways or by constraints on the order of gene loss driven by the intrinsic toxicity of the pathway. Our finding that most of the bikaverin pathway genes have been lost or pseudogenized in these Botrytis genomes suggests that the incidence of HGT of gene cluster-encoded metabolic pathways might be higher than what is possible to be inferred from isolated genome analyses.