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

Machine learning enables identification of an alternative yeast galactose utilization pathway

How genomic differences contribute to phenotypic differences is a major question in biology. The recently characterized genomes, isolation environments, and qualitative patterns of growth on 122 sources and conditions of 1,154 strains from 1,049 fungal species (nearly all known) in the yeast subphylum Saccharomycotina provide a powerful, yet complex, dataset for addressing this question. We used a random forest algorithm trained on these genomic, metabolic, and environmental data to predict growth on several carbon sources with high accuracy. Known structural genes involved in assimilation of these sources and presence/absence patterns of growth in other sources were important features contributing to prediction accuracy. By further examining growth on galactose, we found that it can be predicted with high accuracy from either genomic (92.2%) or growth data (82.6%) but not from isolation environment data (65.6%). Prediction accuracy was even higher (93.3%) when we combined genomic and growth data. After the GALactose utilization genes, the most important feature for predicting growth on galactose was growth on galactitol, raising the hypothesis that several species in two orders, Serinales and Pichiales (containing the emerging pathogen Candida auris and the genus Ogataea, respectively), have an alternative galactose utilization pathway because they lack the GAL genes. Growth and biochemical assays confirmed that several of these species utilize galactose through an alternative oxidoreductive D-galactose pathway, rather than the canonical GAL pathway. Machine learning approaches are powerful for investigating the evolution of the yeast genotype-phenotype map, and their application will uncover novel biology, even in well-studied traits.

Genomic insights into the evolution and adaptation of secondary metabolite gene clusters in fungicolous species Cladobotryum mycophilum ATHUM6906

Mycophilic or fungicolous fungi can be found wherever fungi exist since they are able to colonize other fungi, which occupy a diverse range of habitats. Some fungicolous species cause important diseases on Basidiomycetes, and thus, they are the main reason for the destruction of mushroom cultivations. Nonetheless, despite their ecological significance, their genomic data remain limited. Cladobotryum mycophilum is one of the most aggressive species of the genus, destroying the economically important Agaricus bisporus cultivations. The 40.7 Mb whole genome of the Greek isolate ATHUM6906 is assembled in 16 fragments, including the mitochondrial genome and 2 small circular mitochondrial plasmids, in this study. This genome includes a comprehensive set of 12,282 protein coding, 56 rRNA, and 273 tRNA genes. Transposable elements, CAZymes, and pathogenicity related genes were also examined. The genome of C. mycophilum contained a diverse arsenal of genes involved in secondary metabolism, forming 106 biosynthetic gene clusters, which renders this genome as one of the most BGC abundant among fungicolous species. Comparative analyses were performed for genomes of species of the family Hypocreaceae. Some BGCs identified in C. mycophilum genome exhibited similarities to clusters found in the family Hypocreaceae, suggesting vertical heritage. In contrast, certain BGCs showed a scattered distribution among Hypocreaceae species or were solely found in Cladobotryum genomes. This work provides evidence of extensive BGC losses, horizontal gene transfer events, and formation of novel BGCs during evolution, potentially driven by neutral or even positive selection pressures. These events may increase Cladobotryum fitness under various environmental conditions and potentially during host-fungus interaction.

Transcription Factor Condensates Mediate Clustering of MET Regulon and Enhancement in Gene Expression

Some transcription factors (TFs) can form liquid-liquid phase separated (LLPS) condensates. However, the functions of these TF condensates in 3D genome organization and gene regulation remain elusive. In response to methionine (met) starvation, budding yeast TF Met4 and a few co-activators, including Met32, induce a set of genes involved in met biosynthesis. Here, we show that the endogenous Met4 and Met32 form co-localized puncta-like structures in yeast nuclei upon met depletion. Recombinant Met4 and Met32 form mixed droplets with LLPS properties in vitro. In relation to chromatin, Met4 puncta co-localize with target genes, and at least a subset of these target genes are clustered in 3D in a Met4-dependent manner. A MET3pr-GFP reporter inserted near several native Met4 binding sites becomes co-localized with Met4 puncta and displays enhanced transcriptional activity. A Met4 variant with a partial truncation of an intrinsically disordered region (IDR) shows less puncta formation, and this mutant selectively reduces the reporter activity near Met4 binding sites to the basal level. Overall, these results support a model where Met4 and co-activators form condensates to bring multiple target genes into a vicinity with higher local TF concentrations, which facilitates a strong response to methionine depletion.

Functional Clustering of Metabolically Related Genes Is Conserved across Dikarya

Transcriptional regulation is vital for organismal survival, with many layers and mechanisms collaborating to balance gene expression. One layer of this regulation is genome organization, specifically the clustering of functionally related, co-expressed genes along the chromosomes. Spatial organization allows for position effects to stabilize RNA expression and balance transcription, which can be advantageous for a number of reasons, including reductions in stochastic influences between the gene products. The organization of co-regulated gene families into functional clusters occurs extensively in Ascomycota fungi. However, this is less characterized within the related Basidiomycota fungi despite the many uses and applications for the species within this clade. This review will provide insight into the prevalence, purpose, and significance of the clustering of functionally related genes across Dikarya, including foundational studies from Ascomycetes and the current state of our understanding throughout representative Basidiomycete species.

EvolClustDB: Exploring Eukaryotic Gene Clusters with Evolutionarily Conserved Genomic Neighbourhoods

Conservation of gene neighbourhood over evolutionary distances is generally indicative of shared regulation or functional association among genes. This concept has been broadly exploited in prokaryotes but its use on eukaryotic genomes has been limited to specific functional classes, such as biosynthetic gene clusters. We here used an evolutionary-based gene cluster discovery algorithm (EvolClust) to pre-compute evolutionarily conserved gene neighbourhoods, which can be searched, browsed and downloaded in EvolClustDB. We inferred ∼35,000 cluster families in 882 different species in genome comparisons of five taxonomically broad clades: Fungi, Plants, Metazoans, Insects and Protists. EvolClustDB allows browsing through the cluster families, as well as searching by protein, species, identifier or sequence. Visualization allows inspecting gene order per species in a phylogenetic context, so that relevant evolutionary events such as gain, loss or transfer, can be inferred. EvolClustDB is freely available, without registration, at http://evolclustdb.org/.

Recognition of galactose by a scaffold protein recruits a transcriptional activator for the GAL regulon induction in Candida albicans

The GAL pathway of yeasts has long served as a model system for understanding of how regulatory mode of eukaryotic metabolic pathways evolves. While Gal4 mode has been well-characterized in Saccharomycetaceae clade, little is known about the regulation of the GAL pathway in other yeasts. Here, we find that Rep1, a Ndt80-like family transcription factor, serves as a galactose sensor in the commensal-pathogenic fungus Candida albicans. It is presented at the GAL gene promoters independent of the presence of galactose. Rep1 recognizes galactose via a direct physical interaction. The net result of this interaction is the recruitment of a transcriptional activator Cga1 (Candida galactose gene activator, orf19.4959) and transcription of the GAL genes proceeds. Rep1 and Cga1 are conserved across the CTG species. Rep1 itself does not possess transcriptional activity. Instead, it provides a scaffold to recruit different factors for transcriptional regulation. Rep1-Cga1 mode of regulation represents a new example of network rewiring in fungi, which provides insight into how C. albicans evolves transcriptional programs to colonize diverse host niches.

Uncovering a miltiradiene biosynthetic gene cluster in the Lamiaceae reveals a dynamic evolutionary trajectory

The spatial organization of genes within plant genomes can drive evolution of specialized metabolic pathways. Terpenoids are important specialized metabolites in plants with diverse adaptive functions that enable environmental interactions. Here, we report the genome assemblies of Prunella vulgaris, Plectranthus barbatus, and Leonotis leonurus. We investigate the origin and subsequent evolution of a diterpenoid biosynthetic gene cluster (BGC) together with other seven species within the Lamiaceae (mint) family. Based on core genes found in the BGCs of all species examined across the Lamiaceae, we predict a simplified version of this cluster evolved in an early Lamiaceae ancestor. The current composition of the extant BGCs highlights the dynamic nature of its evolution. We elucidate the terpene backbones generated by the Callicarpa americana BGC enzymes, including miltiradiene and the terpene (+)-kaurene, and show oxidization activities of BGC cytochrome P450s. Our work reveals the fluid nature of BGC assembly and the importance of genome structure in contributing to the origin of metabolites.

Contrasting modes of macro and microsynteny evolution in a eukaryotic subphylum

Examination of the changes in order and arrangement of homologous genes is key for understanding the mechanisms of genome evolution in eukaryotes. Previous comparisons between eukaryotic genomes have revealed considerable conservation across species that diverged hundreds of millions of years ago (e.g., vertebrates,1,2,3 bilaterian animals,4,5 and filamentous fungi6). However, understanding how genome organization evolves within and between eukaryotic major lineages remains underexplored. We analyzed high-quality genomes of 120 representative budding yeast species (subphylum Saccharomycotina) spanning ∼400 million years of eukaryotic evolution to examine how their genome organization evolved and to compare it with the evolution of animal and plant genome organization.7 We found that the decay of both macrosynteny (the conservation of homologous chromosomes) and microsynteny (the conservation of local gene content and order) was strongly associated with evolutionary divergence across budding yeast major clades. However, although macrosynteny decayed very fast, within ∼100 million years, the microsynteny of many genes-especially genes in metabolic clusters (e.g., in the GAL gene cluster8)-was much more deeply conserved both within major clades and across the subphylum. We further found that when genomes with similar evolutionary divergence times were compared, budding yeasts had lower macrosynteny conservation than animals and filamentous fungi but higher conservation than angiosperms. In contrast, budding yeasts had levels of microsynteny conservation on par with mammals, whereas angiosperms exhibited very low conservation. Our results provide new insight into the tempo and mode of the evolution of gene and genome organization across an entire eukaryotic subphylum.

Integrative transcriptome analysis revealed the pathogenic molecular basis of Rhizoctonia solani AG-3 TB at three progressive stages of infection

Rhizoctonia solani has a broad host range and results in significant losses in agricultural production. Here, an integrated transcriptomic analysis was performed to reveal the critical genes responsible for the pathogenesis of R. solani AG-3 TB on Nicotiana tabacum at different infection stages. The results showed that various differential expressed genes (DEGs) were enriched in fatty acid metabolism, amino sugar, carbon metabolism, and cellular carbohydrate biosynthetic process at the early (6–12 hpi), middle (24–36 hpi), and late stage (48–72 hpi) of infection. Specifically, several critical genes such as shikimate kinase that were involved in the biosynthesis of an important fungal toxin, phenylacetic acid (PAA) showed markedly increase at 24 hpi. Additionally, the genes expression levels of carbohydrate-active enzymes (CAZymes) and cell wall degrading enzymes (CWDEs) were significantly increased at the late infection stage. Furthermore, we identified 807 potential secreted proteins and 78 small cysteine-rich proteins, which may function as fungal effectors and involved in the pathogenicity. These results provide valuable insights into critical and potential genes as well as the pathways involved in the pathogenesis of R. solani AG-3 TB.

Horizontal gene transfer in yeasts

Horizontal gene transfer (HGT), defined as the exchange of genetic material other than from parent to progeny, is very common in bacteria and appears to constitute the most important mechanism contributing to enlarge a species gene pool. However, in eukaryotes, HGT is certainly much less common and some early insufficiently consubstantiated cases involving bacterial donors led some to consider that it was unlikely to occur in eukaryotes outside the host/endosymbiont relationship. More recently, plenty of reports of interdomain HGT have seen the light based on the strictest criteria, many concerning filamentous fungi and yeasts. Here, we attempt to summarize the most prominent instances of HGT reported in yeasts as well as what we have been able to learn so far concerning frequency and distribution, mechanisms, barriers, function of horizontally acquired genes, and the role of HGT in domestication.

Neighbours matter: Effects of genomic organization on gene expression plasticity in response to environmental stresses during biological invasions

Gene expression regulation has been widely recognized as an important molecular mechanism underlying phenotypic plasticity in environmental adaptation. However, it remains largely unexplored on the effects of genomic organization on gene expression plasticity under environmental stresses during biological invasions. Here, we use an invasive model ascidian, Ciona robusta, to investigate how genomic organization affects gene expression in response to salinity stresses during range expansions. Our study showed that neighboring genes were co-expressed and approximately 30% of stress responsive genes were physically clustered on chromosomes. Such coordinated expression was substantially affected by the physical distance and orientation of genes. Interestingly, the overall expression correlation of neighboring genes was significantly decreased under high salinity stresses, illustrating that the co-expression regulation could be disrupted by salinity challenges. Furthermore, the clustering of genes was associated with their function constraints and expression patterns - operon genes enriched in gene expression machinery had the highest transcriptional activity and expression stability. Notably, our analyses showed that the tail-to-tail genes, mainly involved in biological functions related to phosphorylation, homeostatic process, and ion transport, exhibited higher intrinsic expression variability and greater response to salinity challenges. Altogether, the results obtained here provide new insights into the effects of gene organization on gene expression plasticity under environmental challenges, hence improving our knowledge on mechanisms of rapid environmental adaptation during biological invasions.

Horizontal transfer and evolution of the biosynthetic gene cluster for benzoxazinoids in plants

Benzoxazinoids are a class of protective and allelopathic plant secondary metabolites that have been identified in multiple grass species and are encoded by the Bx biosynthetic gene cluster (BGC) in maize. Data mining of 41 high-quality grass genomes identified complete Bx clusters (containing genes Bx1-Bx5 and Bx8) in three genera (Zea, Echinochloa, and Dichanthelium) of Panicoideae and partial clusters in Triticeae. The Bx cluster probably originated from gene duplication and chromosomal translocation of native homologs of Bx genes. An ancient Bx cluster that included additional Bx genes (e.g., Bx6) is presumed to have been present in ancestral Panicoideae. The ancient Bx cluster was putatively gained by the Triticeae ancestor via horizontal transfer (HT) from the ancestral Panicoideae and later separated into multiple segments on different chromosomes. Bx6 appears to have been under less constrained selection compared with the Bx cluster during the evolution of Panicoideae, as evidenced by the fact that it was translocated away from the Bx cluster in Zea mays, moved to other chromosomes in Echinochloa, and even lost in Dichanthelium. Further investigations indicate that purifying selection and polyploidization have shaped the evolutionary trajectory of Bx clusters in the grass family. This study provides the first candidate case of HT of a BGC between plants and sheds new light on the evolution of BGCs.

Tracking of Diversity and Evolution in the Brown Rot Fungi Monilinia fructicola, Monilinia fructigena, and Monilinia laxa

Monilinia species are among the most devastating fungi worldwide as they cause brown rot and blossom blight on fruit trees. To understand the molecular bases of their pathogenic lifestyles, we compared the newly assembled genomes of single strains of Monilinia fructicola, M. fructigena and M. laxa, with those of Botrytis cinerea and Sclerotinia sclerotiorum, as the closest species within Sclerotiniaceae. Phylogenomic analysis of orthologous proteins and syntenic investigation suggest that M. laxa is closer to M. fructigena than M. fructicola, and is closest to the other investigated Sclerotiniaceae species. This indicates that M. laxa was the earliest result of the speciation process. Distinct evolutionary profiles were observed for transposable elements (TEs). M. fructicola and M. laxa showed older bursts of TE insertions, which were affected (mainly in M. fructicola) by repeat-induced point (RIP) mutation gene silencing mechanisms. These suggested frequent occurrence of the sexual process in M. fructicola. More recent TE expansion linked with low RIP action was observed in M. fructigena, with very little in S. sclerotiorum and B. cinerea. The detection of active non-syntenic TEs is indicative of horizontal gene transfer and has resulted in alterations in specific gene functions. Analysis of candidate effectors, biosynthetic gene clusters for secondary metabolites and carbohydrate-active enzymes, indicated that Monilinia genus has multiple virulence mechanisms to infect host plants, including toxins, cell-death elicitor, putative virulence factors and cell-wall-degrading enzymes. Some species-specific pathogenic factors might explain differences in terms of host plant and organ preferences between M. fructigena and the other two Monilinia species.

The evolution of the GALactose utilization pathway in budding yeasts

The Leloir galactose utilization or GAL pathway of budding yeasts, including that of the baker's yeast Saccharomyces cerevisiae and the opportunistic human pathogen Candida albicans, breaks down the sugar galactose for energy and biomass production. The GAL pathway has long served as a model system for understanding how eukaryotic metabolic pathways, including their modes of regulation, evolve. More recently, the physical linkage of the structural genes GAL1, GAL7, and GAL10 in diverse budding yeast genomes has been used as a model for understanding the evolution of gene clustering. In this review, we summarize exciting recent work on three different aspects of this iconic pathway's evolution: gene cluster organization, GAL gene regulation, and the population genetics of the GAL pathway.

Three chromosome-scale Papaver genomes reveal punctuated patchwork evolution of the morphinan and noscapine biosynthesis pathway

For millions of years, plants evolve plenty of structurally diverse secondary metabolites (SM) to support their sessile lifestyles through continuous biochemical pathway innovation. While new genes commonly drive the evolution of plant SM pathway, how a full biosynthetic pathway evolves remains poorly understood. The evolution of pathway involves recruiting new genes along the reaction cascade forwardly, backwardly, or in a patchwork manner. With three chromosome-scale Papaver genome assemblies, we here reveal whole-genome duplications (WGDs) apparently accelerate chromosomal rearrangements with a nonrandom distribution towards SM optimization. A burst of structural variants involving fusions, translocations and duplications within 7.7 million years have assembled nine genes into the benzylisoquinoline alkaloids gene cluster, following a punctuated patchwork model. Biosynthetic gene copies and their total expression matter to morphinan production. Our results demonstrate how new genes have been recruited from a WGD-induced repertoire of unregulated enzymes with promiscuous reactivities to innovate efficient metabolic pathways with spatiotemporal constraint.

Genome organization and evolution of a eukaryotic nicotinate co-inducible pathway

In Aspergillus nidulans a regulon including 11 hxn genes (hxnS, T, R, P, Y, Z, X, W, V, M and N) is inducible by a nicotinate metabolic derivative, repressible by ammonium and under stringent control of the nitrogen-state-sensitive GATA factor AreA and the specific transcription factor HxnR. This is the first report in a eukaryote of the genomic organization of a possibly complete pathway of nicotinate utilization. In A. nidulans the regulon is organized in three distinct clusters, this organization is variable in the Ascomycota. In some Pezizomycotina species all 11 genes map in a single cluster; in others they map in two clusters. This variable organization sheds light on cluster evolution. Instances of gene duplication followed by or simultaneous with integration in the cluster, partial or total cluster loss, and horizontal gene transfer of several genes (including an example of whole cluster re-acquisition in Aspergillus of section Flavi) were detected, together with the incorporation in some clusters of genes not found in the A. nidulans co-regulated regulon, which underlie both the plasticity and the reticulate character of metabolic cluster evolution. This study provides a comprehensive phylogeny of six members of the cluster across representatives of all Ascomycota classes.

Repeated horizontal gene transfer of GALactose metabolism genes violates Dollo's law of irreversible loss

Dollo's law posits that evolutionary losses are irreversible, thereby narrowing the potential paths of evolutionary change. While phenotypic reversals to ancestral states have been observed, little is known about their underlying genetic causes. The genomes of budding yeasts have been shaped by extensive reductive evolution, such as reduced genome sizes and the losses of metabolic capabilities. However, the extent and mechanisms of trait reacquisition after gene loss in yeasts have not been thoroughly studied. Here, through phylogenomic analyses, we reconstructed the evolutionary history of the yeast galactose utilization pathway and observed widespread and repeated losses of the ability to utilize galactose, which occurred concurrently with the losses of GALactose (GAL) utilization genes. Unexpectedly, we detected multiple galactose-utilizing lineages that were deeply embedded within clades that underwent ancient losses of galactose utilization. We show that at least two, and possibly three, lineages reacquired the GAL pathway via yeast-to-yeast horizontal gene transfer. Our results show how trait reacquisition can occur tens of millions of years after an initial loss via horizontal gene transfer from distant relatives. These findings demonstrate that the losses of complex traits and even whole pathways are not always evolutionary dead-ends, highlighting how reversals to ancestral states can occur.

Plant Metabolic Gene Clusters: Evolution, Organization, and Their Applications in Synthetic Biology

Plants are a remarkable source of high-value specialized metabolites having significant physiological and ecological functions. Genes responsible for synthesizing specialized metabolites are often clustered together for a coordinated expression, which is commonly observed in bacteria and filamentous fungi. Similar to prokaryotic gene clustering, plants do have gene clusters encoding enzymes involved in the biosynthesis of specialized metabolites. More than 20 gene clusters involved in the biosynthesis of diverse metabolites have been identified across the plant kingdom. Recent studies demonstrate that gene clusters are evolved through gene duplications and neofunctionalization of primary metabolic pathway genes. Often, these clusters are tightly regulated at nucleosome level. The prevalence of gene clusters related to specialized metabolites offers an attractive possibility of an untapped source of highly useful biomolecules. Accordingly, the identification and functional characterization of novel biosynthetic pathways in plants need to be worked out. In this review, we summarize insights into the evolution of gene clusters and discuss the organization and importance of specific gene clusters in the biosynthesis of specialized metabolites. Regulatory mechanisms which operate in some of the important gene clusters have also been briefly described. Finally, we highlight the importance of gene clusters to develop future metabolic engineering or synthetic biology strategies for the heterologous production of novel metabolites.

Application of an optimized annotation pipeline to the Cryptococcus deuterogattii genome reveals dynamic primary metabolic gene clusters and genomic impact of RNAi loss

Evaluating the quality of a de novo annotation of a complex fungal genome based on RNA-seq data remains a challenge. In this study, we sequentially optimized a Cufflinks-CodingQuary-based bioinformatics pipeline fed with RNA-seq data using the manually annotated model pathogenic yeasts Cryptococcus neoformans and Cryptococcus deneoformans as test cases. Our results show that the quality of the annotation is sensitive to the quantity of RNA-seq data used and that the best quality is obtained with 5–10 million reads per RNA-seq replicate. We also showed that the number of introns predicted is an excellent a priori indicator of the quality of the final de novo annotation. We then used this pipeline to annotate the genome of the RNAi-deficient species Cryptococcus deuterogattii strain R265 using RNA-seq data. Dynamic transcriptome analysis revealed that intron retention is more prominent in C. deuterogattii than in the other RNAi-proficient species C. neoformans and C. deneoformans. In contrast, we observed that antisense transcription was not higher in C. deuterogattii than in the two other Cryptococcus species. Comparative gene content analysis identified 21 clusters enriched in transcription factors and transporters that have been lost. Interestingly, analysis of the subtelomeric regions in these three annotated species identified a similar gene enrichment, reminiscent of the structure of primary metabolic clusters. Our data suggest that there is active exchange between subtelomeric regions, and that other chromosomal regions might participate in adaptive diversification of Cryptococcus metabolite assimilation potential.

Signatures of optimal codon usage in metabolic genes inform budding yeast ecology

Reverse ecology is the inference of ecological information from patterns of genomic variation. One rich, heretofore underutilized, source of ecologically relevant genomic information is codon optimality or adaptation. Bias toward codons that match the tRNA pool is robustly associated with high gene expression in diverse organisms, suggesting that codon optimization could be used in a reverse ecology framework to identify highly expressed, ecologically relevant genes. To test this hypothesis, we examined the relationship between optimal codon usage in the classic galactose metabolism (GAL) pathway and known ecological niches for 329 species of budding yeasts, a diverse subphylum of fungi. We find that optimal codon usage in the GAL pathway is positively correlated with quantitative growth on galactose, suggesting that GAL codon optimization reflects increased capacity to grow on galactose. Optimal codon usage in the GAL pathway is also positively correlated with human-associated ecological niches in yeasts of the CUG-Ser1 clade and with dairy-associated ecological niches in the family Saccharomycetaceae. For example, optimal codon usage of GAL genes is greater than 85% of all genes in the genome of the major human pathogen Candida albicans (CUG-Ser1 clade) and greater than 75% of genes in the genome of the dairy yeast Kluyveromyces lactis (family Saccharomycetaceae). We further find a correlation between optimization in the GALactose pathway genes and several genes associated with nutrient sensing and metabolism. This work suggests that codon optimization harbors information about the metabolic ecology of microbial eukaryotes. This information may be particularly useful for studying fungal dark matter-species that have yet to be cultured in the lab or have only been identified by genomic material.

Fungal natural products galaxy: Biochemistry and molecular genetics toward blockbuster drugs discovery

Secondary metabolites synthesized by fungi have become a precious source of inspiration for the design of novel drugs. Indeed, fungi are prolific producers of fascinating, diverse, structurally complex, and low-molecular-mass natural products with high therapeutic leads, such as novel antimicrobial compounds, anticancer compounds, immunosuppressive agents, among others. Given that these microorganisms possess the extraordinary capacity to secrete diverse chemical scaffolds, they have been highly exploited by the giant pharma companies to generate small molecules. This has been made possible because the isolation of metabolites from fungal natural sources is feasible and surpasses the organic synthesis of compounds, which otherwise remains a significant bottleneck in the drug discovery process. Here in this comprehensive review, we have discussed recent studies on different fungi (pathogenic, non-pathogenic, commensal, and endophytic/symbiotic) from different habitats (terrestrial and marines), the specialized metabolites they biosynthesize, and the drugs derived from these specialized metabolites. Moreover, we have unveiled the logic behind the biosynthesis of vital chemical scaffolds, such as NRPS, PKS, PKS-NRPS hybrid, RiPPS, terpenoids, indole alkaloids, and their genetic mechanisms. Besides, we have provided a glimpse of the concept behind mycotoxins, virulence factor, and host immune response based on fungal infections.

Giant GAL gene clusters for the melibiose-galactose pathway in Torulaspora

In many yeast species, the three genes at the centre of the galactose catabolism pathway, GAL1, GAL10 and GAL7, are neighbours in the genome and form a metabolic gene cluster. We report here that some yeast strains in the genus Torulaspora have much larger GAL clusters that include genes for melibiase (MEL1), galactose permease (GAL2), glucose transporter (HGT1), phosphoglucomutase (PGM1) and the transcription factor GAL4, in addition to GAL1, GAL10, and GAL7. Together, these eight genes encode almost all the steps in the pathway for catabolism of extracellular melibiose (a disaccharide of galactose and glucose). We show that a progenitor 5-gene cluster containing GAL 7-1-10-4-2 was likely present in the common ancestor of Torulaspora and Zygotorulaspora. It added PGM1 and MEL1 in the ancestor of most Torulaspora species. It underwent further expansion in the T. pretoriensis clade, involving the fusion of three progenitor clusters in tandem and the gain of HGT1. These giant GAL clusters are highly polymorphic in structure, and subject to horizontal transfers, pseudogenization and gene losses. We identify recent horizontal transfers of complete GAL clusters from T. franciscae into one strain of T. delbrueckii, and from a relative of T. maleeae into one strain of T. globosa. The variability and dynamic evolution of GAL clusters in Torulaspora indicates that there is strong natural selection on the GAL pathway in this genus.

The Path towards Predicting Evolution as Illustrated in Yeast Cell Polarity

A bottom-up route towards predicting evolution relies on a deep understanding of the complex network that proteins form inside cells. In a rapidly expanding panorama of experimental possibilities, the most difficult question is how to conceptually approach the disentangling of such complex networks. These can exhibit varying degrees of hierarchy and modularity, which obfuscate certain protein functions that may prove pivotal for adaptation. Using the well-established polarity network in budding yeast as a case study, we first organize current literature to highlight protein entrenchments inside polarity. Following three examples, we see how alternating between experimental novelties and subsequent emerging design strategies can construct a layered understanding, potent enough to reveal evolutionary targets. We show that if you want to understand a cell’s evolutionary capacity, such as possible future evolutionary paths, seemingly unimportant proteins need to be mapped and studied. Finally, we generalize this research structure to be applicable to other systems of interest.

Genomic clustering within functionally related gene families in Ascomycota fungi

Multiple mechanisms collaborate for proper regulation of gene expression. One layer of this regulation is through the clustering of functionally related genes at discrete loci throughout the genome. This phenomenon occurs extensively throughout Ascomycota fungi and is an organizing principle for many gene families whose proteins participate in diverse molecular functions throughout the cell. Members of this phylum include organisms that serve as model systems and those of interest medically, pharmaceutically, and for industrial and biotechnological applications. In this review, we discuss the prevalence of functional clustering through a broad range of organisms within the phylum. Position effects on transcription, genomic locations of clusters, transcriptional regulation of clusters, and selective pressures contributing to the formation and maintenance of clusters are addressed, as are common methods to identify and characterize clusters.

Horizontal Gene Transfer and Endophytes: An Implication for the Acquisition of Novel Traits

Horizontal gene transfer (HGT), an important evolutionary mechanism observed in prokaryotes, is the transmission of genetic material across phylogenetically distant species. In recent years, the availability of complete genomes has facilitated the comprehensive analysis of HGT and highlighted its emerging role in the adaptation and evolution of eukaryotes. Endophytes represent an ecologically favored association, which highlights its beneficial attributes to the environment, in agriculture and in healthcare. The HGT phenomenon in endophytes, which features an important biological mechanism for their evolutionary adaptation within the host plant and simultaneously confers "novel traits" to the associated microbes, is not yet completely understood. With a focus on the emerging implications of HGT events in the evolution of biological species, the present review discusses the occurrence of HGT in endophytes and its socio-economic importance in the current perspective. To our knowledge, this review is the first report that provides a comprehensive insight into the impact of HGT in the adaptation and evolution of endophytes.

Complex Evolutionary Origins of Specialized Metabolite Gene Cluster Diversity among the Plant Pathogenic Fungi of the Fusarium graminearum Species Complex

Fungal genomes encode highly organized gene clusters that underlie the production of specialized (or secondary) metabolites. Gene clusters encode key functions to exploit plant hosts or environmental niches. Promiscuous exchange among species and frequent reconfigurations make gene clusters some of the most dynamic elements of fungal genomes. Despite evidence for high diversity in gene cluster content among closely related strains, the microevolutionary processes driving gene cluster gain, loss, and neofunctionalization are largely unknown. We analyzed the Fusarium graminearum species complex (FGSC) composed of plant pathogens producing potent mycotoxins and causing Fusarium head blight on cereals. We de novo assembled genomes of previously uncharacterized FGSC members (two strains of F. austroamericanum, F. cortaderiae, and F. meridionale). Our analyses of 8 species of the FGSC in addition to 15 other Fusarium species identified a pangenome of 54 gene clusters within FGSC. We found that multiple independent losses were a key factor generating extant cluster diversity within the FGSC and the Fusarium genus. We identified a modular gene cluster conserved among distantly related fungi, which was likely reconfigured to encode different functions. We also found strong evidence that a rare cluster in FGSC was gained through an ancient horizontal transfer between bacteria and fungi. Chromosomal rearrangements underlying cluster loss were often complex and were likely facilitated by an enrichment in specific transposable elements. Our findings identify important transitory stages in the birth and death process of specialized metabolism gene clusters among very closely related species.

Multilayered horizontal operon transfers from bacteria reconstruct a thiamine salvage pathway in yeasts

Horizontal acquisition of bacterial genes is presently recognized as an important contribution to the adaptation and evolution of eukaryotic genomes. However, the mechanisms underlying expression and consequent selection and fixation of the prokaryotic genes in the new eukaryotic setting are largely unknown. Here we show that genes composing the pathway for the synthesis of the essential vitamin B1 (thiamine) were lost in an ancestor of a yeast lineage, the Wickerhamiella/Starmerella (W/S) clade, known to harbor an unusually large number of genes of alien origin. The thiamine pathway was subsequently reassembled, at least twice, by multiple HGT events from different bacterial donors involving both single genes and entire operons. In the W/S-clade species Starmerella bombicola we obtained direct genetic evidence that all bacterial genes of the thiamine pathway are functional. The reconstructed pathway is composed by yeast and bacterial genes operating coordinately to scavenge thiamine derivatives from the environment. The adaptation of the newly acquired operons to the eukaryotic setting involved a repertoire of mechanisms until now only sparsely documented, namely longer intergenic regions, post-horizontal gene transfer (HGT) gene fusions fostering coordinated expression, gene relocation, and possibly recombination generating mosaic genes. The results provide additional evidence that HGT occurred recurrently in this yeast lineage and was crucial for the reestablishment of lost functions and that similar mechanisms are used across a broad range of eukaryotic microbes to promote adaptation of prokaryotic genes to their new environment.

Synchronization of stochastic expressions drives the clustering of functionally related genes

Functionally related genes tend to be chromosomally clustered in eukaryotic genomes even after the exclusion of tandem duplicates, but the biological significance of this widespread phenomenon is unclear. We propose that stochastic expression fluctuations of neighboring genes resulting from chromatin dynamics are more or less synchronized such that their expression ratio is more stable than that for unlinked genes. Consequently, chromosomal clustering could be advantageous when the expression ratio of the clustered genes needs to stay constant, for example, because of the accumulation of toxic compounds when this ratio is altered. Evidence from manipulative experiments on the yeast GAL cluster, comprising three chromosomally adjacent genes encoding enzymes catalyzing consecutive reactions in galactose catabolism, unequivocally supports this hypothesis and elucidates how disorder in one biological phenomenon-gene expression noise-could prompt the emergence of order in another-genome organization.

Evolutionary and functional patterns of shared gene neighbourhood in fungi

Gene clusters comprise genomically co-localized and potentially co-regulated genes that tend to be conserved across species. In eukaryotes, multiple examples of metabolic gene clusters are known, particularly among fungi and plants. However, little is known about how gene clustering patterns vary among taxa or with respect to functional roles. Furthermore, mechanisms of the formation, maintenance and evolution of gene clusters remain unknown. We surveyed 341 fungal genomes to discover gene clusters shared by different species, independently of their functions. We inferred 12,120 cluster families, which comprised roughly one third of the gene space and were enriched in genes associated with diverse cellular functions. Additionally, most clusters did not encode transcription factors, suggesting that they are regulated distally. We used phylogenomics to characterize the evolutionary history of these clusters. We found that most clusters originated once and were transmitted vertically, coupled to differential loss. However, convergent evolution-that is, independent appearance of the same cluster-was more prevalent than anticipated. Finally, horizontal gene transfer of entire clusters was somewhat restricted, with the exception of those associated with secondary metabolism. Altogether, our results provide insights on the evolution of gene clustering as well as a broad catalogue of evolutionarily conserved gene clusters whose function remains to be elucidated.

The birth, evolution and death of metabolic gene clusters in fungi

Fungi contain a remarkable diversity of both primary and secondary metabolic pathways involved in ecologically specialized or accessory functions. Genes in these pathways are frequently physically linked on fungal chromosomes, forming metabolic gene clusters (MGCs). In this Review, we describe the diversity in the structure and content of fungal MGCs, their population-level and species-level variation, the evolutionary mechanisms that underlie their formation, maintenance and decay, and their ecological and evolutionary impact on fungal populations. We also discuss MGCs from other eukaryotes and the reasons for their preponderance in fungi. Improved knowledge of the evolutionary life cycle of MGCs will advance our understanding of the ecology of specialized metabolism and of the interplay between the lifestyle of an organism and genome architecture.

Extensive loss of cell-cycle and DNA repair genes in an ancient lineage of bipolar budding yeasts

Cell-cycle checkpoints and DNA repair processes protect organisms from potentially lethal mutational damage. Compared to other budding yeasts in the subphylum Saccharomycotina, we noticed that a lineage in the genus Hanseniaspora exhibited very high evolutionary rates, low Guanine-Cytosine (GC) content, small genome sizes, and lower gene numbers. To better understand Hanseniaspora evolution, we analyzed 25 genomes, including 11 newly sequenced, representing 18/21 known species in the genus. Our phylogenomic analyses identify two Hanseniaspora lineages, a faster-evolving lineage (FEL), which began diversifying approximately 87 million years ago (mya), and a slower-evolving lineage (SEL), which began diversifying approximately 54 mya. Remarkably, both lineages lost genes associated with the cell cycle and genome integrity, but these losses were greater in the FEL. E.g., all species lost the cell-cycle regulator WHIskey 5 (WHI5), and the FEL lost components of the spindle checkpoint pathway (e.g., Mitotic Arrest-Deficient 1 [MAD1], Mitotic Arrest-Deficient 2 [MAD2]) and DNA-damage-checkpoint pathway (e.g., Mitosis Entry Checkpoint 3 [MEC3], RADiation sensitive 9 [RAD9]). Similarly, both lineages lost genes involved in DNA repair pathways, including the DNA glycosylase gene 3-MethylAdenine DNA Glycosylase 1 (MAG1), which is part of the base-excision repair pathway, and the DNA photolyase gene PHotoreactivation Repair deficient 1 (PHR1), which is involved in pyrimidine dimer repair. Strikingly, the FEL lost 33 additional genes, including polymerases (i.e., POLymerase 4 [POL4] and POL32) and telomere-associated genes (e.g., Repressor/activator site binding protein-Interacting Factor 1 [RIF1], Replication Factor A 3 [RFA3], Cell Division Cycle 13 [CDC13], Pbp1p Binding Protein [PBP2]). Echoing these losses, molecular evolutionary analyses reveal that, compared to the SEL, the FEL stem lineage underwent a burst of accelerated evolution, which resulted in greater mutational loads, homopolymer instabilities, and higher fractions of mutations associated with the common endogenously damaged base, 8-oxoguanine. We conclude that Hanseniaspora is an ancient lineage that has diversified and thrived, despite lacking many otherwise highly conserved cell-cycle and genome integrity genes and pathways, and may represent a novel, to our knowledge, system for studying cellular life without them.

A putative antimicrobial peptide from Hymenoptera in the megaplasmid pSCL4 of Streptomyces clavuligerus ATCC 27064 reveals a singular case of horizontal gene transfer with potential applications

Streptomyces clavuligerus is a Gram-positive bacterium that is a high producer of secondary metabolites with industrial applications. The production of antibiotics such as clavulanic acid or cephamycin has been extensively studied in this species; nevertheless, other aspects, such as evolution or ecology, have received less attention. Furthermore, genes that arise from ancient events of lateral transfer have been demonstrated to be implicated in important functions of host species. This approximation discovered relevant genes that genomic analyses overlooked. Thus, we studied the impact of horizontal gene transfer in the S. clavuligerus genome. To perform this task, we applied whole-genome analysis to identify a laterally transferred sequence from different domains. The most relevant result was a putative antimicrobial peptide (AMP) with a clear origin in the Hymenoptera order of insects. Next, we determined that two copies of these genes were present in the megaplasmid pSCL4 but absent in the S. clavuligerus ATCC 27064 chromosome. Additionally, we found that these sequences were exclusive to the ATCC 27064 strain (and so were not present in any other bacteria) and we also verified the expression of the genes using RNAseq data. Next, we used several AMP predictors to validate the original annotation extracted from Hymenoptera sequences and explored the possibility that these proteins had post-translational modifications using peptidase cleavage prediction. We suggest that Hymenoptera AMP-like proteins of S. clavuligerus ATCC 27064 may be useful for both species adaptation and as an antimicrobial molecule with industrial applications.

Eukaryotic Acquisition of a Bacterial Operon

Operons are a hallmark of bacterial genomes, where they allow concerted expression of functionally related genes as single polycistronic transcripts. They are rare in eukaryotes, where each gene usually drives expression of its own independent messenger RNAs. Here, we report the horizontal operon transfer of a siderophore biosynthesis pathway from relatives of Escherichia coli into a group of budding yeast taxa. We further show that the co-linearly arranged secondary metabolism genes are expressed, exhibit eukaryotic transcriptional features, and enable the sequestration and uptake of iron. After transfer, several genetic changes occurred during subsequent evolution, including the gain of new transcription start sites that were sometimes within protein-coding sequences, acquisition of polyadenylation sites, structural rearrangements, and integration of eukaryotic genes into the cluster. We conclude that the genes were likely acquired as a unit, modified for eukaryotic gene expression, and maintained by selection to adapt to the highly competitive, iron-limited environment.

Pan-genome analyses of model fungal species

The concept of the species ‘pan-genome’, the union of ‘core’ conserved genes and all ‘accessory’ non-conserved genes across all strains of a species, was first proposed in prokaryotes to account for intraspecific variability. Species pan-genomes have been extensively studied in prokaryotes, but evidence of species pan-genomes has also been demonstrated in eukaryotes such as plants and fungi. Using a previously published methodology based on sequence homology and conserved microsynteny, in addition to bespoke pipelines, we have investigated the pan-genomes of four model fungal species: Saccharomyces cerevisiae, Candida albicans, Cryptococcus neoformans var. grubii and Aspergillus fumigatus. Between 80 and 90 % of gene models per strain in each of these species are core genes that are highly conserved across all strains of that species, many of which are involved in housekeeping and conserved survival processes. In many of these species, the remaining ‘accessory’ gene models are clustered within subterminal regions and may be involved in pathogenesis and antimicrobial resistance. Analysis of the ancestry of species core and accessory genomes suggests that fungal pan-genomes evolve by strain-level innovations such as gene duplication as opposed to wide-scale horizontal gene transfer. Our findings lend further supporting evidence to the existence of species pan-genomes in eukaryote taxa.

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.

Metabolic Gene Clusters in Eukaryotes

In bacteria, more than half of the genes in the genome are organized in operons. In contrast, in eukaryotes, functionally related genes are usually dispersed across the genome. There are, however, numerous examples of functional clusters of nonhomologous genes for metabolic pathways in fungi and plants. Despite superficial similarities with operons (physical clustering, coordinate regulation), these clusters have not usually originated by horizontal gene transfer from bacteria, and (unlike operons) the genes are typically transcribed separately rather than as a single polycistronic message. This clustering phenomenon raises intriguing questions about the origins of clustered metabolic pathways in eukaryotes and the significance of clustering for pathway function. Here we review metabolic gene clusters from fungi and plants, highlight commonalities and differences, and consider how these clusters form and are regulated. We also identify opportunities for future research in the areas of large-scale genomics, synthetic biology, and experimental evolution.

Modeling three-dimensional genomic organization in evolution and pathogenesis

The regulation of gene expression is mediated via the complex three-dimensional (3D) conformation of the genetic material and its interactions with various intracellular factors. Various experimental and computational approaches have been developed in recent years for understating the relation between the 3D conformation of the genome and the phenotypes of cells in normal condition and diseases. In this review, we will discuss novel approaches for analyzing and modeling the 3D genomic conformation, focusing on deciphering disease-causing mutations that affect gene expression. We conclude that as this is a very challenging mission, an important direction should involve the comparative analysis of various 3D models from various organisms or cells.

Genetic Interaction Network as an Important Determinant of Gene Order in Genome Evolution

Although it is generally accepted that eukaryotic gene order is not random, the basic principles of gene arrangement on a chromosome remain poorly understood. Here, we extended existing population genetics theories that were based on two-locus models and proposed a hypothesis that genetic interaction networks drive the evolution of eukaryotic gene order. We predicted that genes with positive epistasis would move toward each other in evolution, during which a negative correlation between epistasis and gene distance formed. We tested and confirmed our prediction with computational simulations and empirical data analyses. Importantly, we demonstrated that gene order in the budding yeast could be successfully predicted from the genetic interaction network. Taken together, our study reveals the role of the genetic interaction network in the evolution of gene order, extends our understanding of the encoding principles in genomes, and potentially offers new strategies to improve synthetic biology.

Genome Mining of Non-Conventional Yeasts: Search and Analysis of MAL Clusters and Proteins

Genomic clustering of functionally related genes is rare in yeasts and other eukaryotes with only few examples available. Here, we summarize our data on a nontelomeric MAL cluster of a non-conventional methylotrophic yeast Ogataea (Hansenula) polymorpha containing genes for α-glucosidase MAL1, α-glucoside permease MAL2 and two hypothetical transcriptional activators. Using genome mining, we detected MAL clusters of varied number, position and composition in many other maltose-assimilating non-conventional yeasts from different phylogenetic groups. The highest number of MAL clusters was detected in Lipomyces starkeyi while no MAL clusters were found in Schizosaccharomyces pombe and Blastobotrys adeninivorans. Phylograms of α-glucosidases and α-glucoside transporters of yeasts agreed with phylogenesis of the respective yeast species. Substrate specificity of unstudied α-glucosidases was predicted from protein sequence analysis. Specific activities of Scheffersomycesstipitis α-glucosidases MAL7, MAL8, and MAL9 heterologously expressed in Escherichia coli confirmed the correctness of the prediction-these proteins were verified promiscuous maltase-isomaltases. α-Glucosidases of earlier diverged yeasts L. starkeyi, B. adeninivorans and S. pombe showed sequence relatedness with α-glucosidases of filamentous fungi and bacilli.

Genome sequence and physiological analysis of Yamadazyma laniorum f.a. sp. nov. and a reevaluation of the apocryphal xylose fermentation of its sister species, Candida tenuis

Xylose fermentation is a rare trait that is immensely important to the cellulosic biofuel industry, and Candida tenuis is one of the few yeasts that has been reported with this trait. Here we report the isolation of two strains representing a candidate sister species to C. tenuis. Integrated analysis of genome sequence and physiology suggested the genetic basis of a number of traits, including variation between the novel species and C. tenuis in lactose metabolism due to the loss of genes encoding lactose permease and β-galactosidase in the former. Surprisingly, physiological characterization revealed that neither the type strain of C. tenuis nor this novel species fermented xylose in traditional assays. We reexamined three xylose-fermenting strains previously identified as C. tenuis and found that these strains belong to the genus Scheffersomyces and are not C. tenuis. We propose Yamadazyma laniorum f.a. sp. nov. to accommodate our new strains and designate its type strain as yHMH7 (=CBS 14780 = NRRL Y-63967T). Furthermore, we propose the transfer of Candida tenuis to the genus Yamadazyma as Yamadazyma tenuis comb. nov. This approach provides a roadmap for how integrated genome sequence and physiological analysis can yield insight into the mechanisms that generate yeast biodiversity.

Evidence for loss and reacquisition of alcoholic fermentation in a fructophilic yeast lineage

Fructophily is a rare trait that consists of the preference for fructose over other carbon sources. Here, we show that in a yeast lineage (the Wickerhamiella/Starmerella, W/S clade) comprised of fructophilic species thriving in the high-sugar floral niche, the acquisition of fructophily is concurrent with a wider remodeling of central carbon metabolism. Coupling comparative genomics with biochemical and genetic approaches, we gathered ample evidence for the loss of alcoholic fermentation in an ancestor of the W/S clade and subsequent reinstatement through either horizontal acquisition of homologous bacterial genes or modification of a pre-existing yeast gene. An enzyme required for sucrose assimilation was also acquired from bacteria, suggesting that the genetic novelties identified in the W/S clade may be related to adaptation to the high-sugar environment. This work shows how even central carbon metabolism can be remodeled by a surge of HGT events.

Cells build their components, such as the molecular machinery that helps them obtain energy from their environment, by following the instructions contained in genes. This genetic information is usually transferred from parents to offspring. Over the course of several generations, genes can accumulate small changes and the molecules they code for can acquire new roles: yet, this process is normally slow. However, certain organisms can also obtain completely new genes by ‘stealing’ them from other species. For example, yeasts, such as the ones used to make bread and beer, can take genes from nearby bacteria. This ‘horizontal gene transfer’ helps organisms to rapidly gain new characteristics, which is particularly useful if the environment changes quickly.

One way that yeasts get the energy they need is by breaking down sugars through a process called alcoholic fermentation. To do this, most yeast species prefer to use a sugar called glucose, but a small group of ‘fructophilic’ species instead favors a type of sugar known as fructose. Scientists do not know exactly how fructophilic yeasts came to be, but there is some evidence horizontal gene transfers may have been involved in the process.

Now, Gonçalves et al. have compared the genetic material of fructophilic yeasts with that of other groups of yeasts . Comparing genetic material helps scientists identify similarities and differences between species, and gives clues about why specific genetic features first evolved.

The experiments show that, early in their history, fructophilic yeasts lost the genes that allowed them to do alcoholic fermentation, probably since they could obtain energy in a different way. However, at a later point in time, these yeasts had to adapt to survive in flower nectar, an environment rich in sugar. They then favored fructose as their source of energy, possibly because this sugar can compensate more effectively for the absence of alcoholic fermentation. Later, the yeasts acquired a gene from nearby bacteria, which allowed them to do alcoholic fermentation again: this improved their ability to use the other sugars present in flower nectars.

When obtaining energy, yeasts and other organisms produce substances that are relevant to industry. Studying natural processes of evolution can help scientists understand how organisms can change the way they get their energy and adapt to new challenges. In turn, this helps to engineer yeasts into ‘cell factories’ that produce valuable chemicals in environmentally friendly and cost-effective ways.

Galactose metabolism and toxicity in Ustilago maydis

In most organisms, galactose is metabolized via the Leloir pathway, which is conserved from bacteria to mammals. Utilization of galactose requires a close interplay of the metabolic enzymes, as misregulation or malfunction of individual components can lead to the accumulation of toxic intermediate compounds. For the phytopathogenic basidiomycete Ustilago maydis, galactose is toxic for wildtype strains, i.e. leads to growth repression despite the presence of favorable carbon sources as sucrose. The galactose sensitivity can be relieved by two independent modifications: (1) by disruption of Hxt1, which we identify as the major transporter for galactose, and (2) by a point mutation in the gene encoding the galactokinase Gal1, the first enzyme of the Leloir pathway. The mutation in gal1(Y67F) leads to reduced enzymatic activity of Gal1 and thus may limit the formation of putatively toxic galactose-1-phosphate. However, systematic deletions and double deletions of different genes involved in galactose metabolism point to a minor role of galactose-1-phosphate in galactose toxicity. Our results show that molecular triggers for galactose toxicity in U. maydis differ from yeast and mammals.

Mapping Causal Variants with Single-Nucleotide Resolution Reveals Biochemical Drivers of Phenotypic Change

Understanding the sequence determinants that give rise to diversity among individuals and species is the central challenge of genetics. However, despite ever greater numbers of sequenced genomes, most genome-wide association studies cannot distinguish causal variants from linked passenger mutations spanning many genes. We report that this inherent challenge can be overcome in model organisms. By pushing the advantages of inbred crossing to its practical limit in Saccharomyces cerevisiae, we improved the statistical resolution of linkage analysis to single nucleotides. This "super-resolution" approach allowed us to map 370 causal variants across 26 quantitative traits. Missense, synonymous, and cis-regulatory mutations collectively gave rise to phenotypic diversity, providing mechanistic insight into the basis of evolutionary divergence. Our data also systematically unmasked complex genetic architectures, revealing that multiple closely linked driver mutations frequently act on the same quantitative trait. Single-nucleotide mapping thus complements traditional deletion and overexpression screening paradigms and opens new frontiers in quantitative genetics.

Identification of a mutarotase gene involved in D-galactose utilization in Aspergillus nidulans

Aldose 1-epimerases or mutarotases (EC 5.1.3.3) are catalyzing the interconversion of α- and β-anomers of hemiacetals of aldose sugars such as D-glucose and D-galactose, and are presumed to play an auxiliary role in carbohydrate metabolism as mutarotation occurs spontaneously in watery solutions. The first step in the Leloir pathway of D-galactose breakdown is preceded by accelerated conversion of β-D-galactopyranose into the α-anomer, the substrate of the anomer-specific D-galactose 1-kinase. Here, we identified two putative aldose-1-epimerase genes (galmA and galmB) in the model organism Aspergillus nidulans, and characterized them upon generation of single- and double deletion mutant strains, as well as overexpressing mutants carrying multiple copies of either. Assaying cell-free extracts from the galmB single- and galm double mutants, we observed that the mutarotation hardly exceeded spontaneous anomer conversion, while galmB multicopy strains displayed higher activities than the wild type, increasing with the copy number. When grown on D-galactose in submerged cultures, biomass formation and D-galactose uptake rates in mutants lacking galmB were considerably reduced. None such effects were observed studying galmA deletion mutants, which consistently behave like the wild type. We conclude that GalmB is the physiologically relevant mutarotase for the utilization of D-galactose in A. nidulans.

Fungal Gene Cluster Diversity and Evolution

Metabolic gene clusters (MGCs) have provided some of the earliest glimpses at the biochemical machinery of yeast and filamentous fungi. MGCs encode diverse genetic mechanisms for nutrient acquisition and the synthesis/degradation of essential and adaptive metabolites. Beyond encoding the enzymes performing these discrete anabolic or catabolic processes, MGCs may encode a range of mechanisms that enable their persistence as genetic consortia; these include enzymatic mechanisms to protect their host fungi from their inherent toxicities, and integrated regulatory machinery. This modular, self-contained nature of MGCs contributes to the metabolic and ecological adaptability of fungi. The phylogenetic and ecological patterns of MGC distribution reflect the broad diversity of fungal life cycles and nutritional modes. While the origins of most gene clusters are enigmatic, MGCs are thought to be born into a genome through gene duplication, relocation, or horizontal transfer, and analyzing the death and decay of gene clusters provides clues about the mechanisms selecting for their assembly. Gene clustering may provide inherent fitness advantages through metabolic efficiency and specialization, but experimental evidence for this is currently limited. The identification and characterization of gene clusters will continue to be powerful tools for elucidating fungal metabolism as well as understanding the physiology and ecology of fungi.

Transcriptome analysis reveals the host selection fitness mechanisms of the Rhizoctonia solani AG1IA pathogen

Rhizoctonia solani AG1IA is a major generalist pathogen that causes sheath blight. Its genome, which was the first to be sequenced from the Rhizoctonia genus, may serve as a model for studying pathogenic mechanisms. To explore the pathogen-host fitness mechanism of sheath-blight fungus, a comprehensive comparative transcriptome ecotype analysis of R. solani AG1IA isolated from rice, soybean and corn during infection was performed. Special characteristics in gene expression, gene ontology terms and expression of pathogenesis-associated genes, including genes encoding secreted proteins, candidate effectors, hydrolases, and proteins involved in secondary metabolite production and the MAPK pathway, were revealed. Furthermore, as an important means of pathogenic modulation, diverse alternative splicing of key pathogenic genes in Rhizoctonia solani AG1IA during infections of the abovementioned hosts was uncovered for the first time. These important findings of key factors in the pathogenicity of R. solani AG1IA ecotypes during infection of various hosts explain host preference and provide novel insights into the pathogenic mechanisms and host-pathogen selection. Furthermore, they provide information on the fitness of Rhizoctonia, a severe pathogen with a wide host range.

Reconstructing the Backbone of the Saccharomycotina Yeast Phylogeny Using Genome-Scale Data

Understanding the phylogenetic relationships among the yeasts of the subphylum Saccharomycotina is a prerequisite for understanding the evolution of their metabolisms and ecological lifestyles. In the last two decades, the use of rDNA and multilocus data sets has greatly advanced our understanding of the yeast phylogeny, but many deep relationships remain unsupported. In contrast, phylogenomic analyses have involved relatively few taxa and lineages that were often selected with limited considerations for covering the breadth of yeast biodiversity. Here we used genome sequence data from 86 publicly available yeast genomes representing nine of the 11 known major lineages and 10 nonyeast fungal outgroups to generate a 1233-gene, 96-taxon data matrix. Species phylogenies reconstructed using two different methods (concatenation and coalescence) and two data matrices (amino acids or the first two codon positions) yielded identical and highly supported relationships between the nine major lineages. Aside from the lineage comprised by the family Pichiaceae, all other lineages were monophyletic. Most interrelationships among yeast species were robust across the two methods and data matrices. However, eight of the 93 internodes conflicted between analyses or data sets, including the placements of: the clade defined by species that have reassigned the CUG codon to encode serine, instead of leucine; the clade defined by a whole genome duplication; and the species Ascoidea rubescens These phylogenomic analyses provide a robust roadmap for future comparative work across the yeast subphylum in the disciplines of taxonomy, molecular genetics, evolutionary biology, ecology, and biotechnology. To further this end, we have also provided a BLAST server to query the 86 Saccharomycotina genomes, which can be found at http://y1000plus.org/blast.

Several steps of lateral gene transfer followed by events of 'birth-and-death' evolution shaped a fungal sorbicillinoid biosynthetic gene cluster

Background

Sorbicillinoids are a family of complex cyclic polyketides produced by only a small number of distantly related ascomycete fungi such as Trichoderma (Sordariomycetes) and Penicillium (Eurotiomycetes). In T. reesei, they are synthesized by a gene cluster consisting of eight genes including two polyketide synthases (PKS). To reconstruct the evolutionary origin of this gene cluster, we examined the occurrence of these eight genes in ascomycetes.

Results

A cluster comprising at least six of them was only found in Hypocreales (Acremonium chrysogenum, Ustilaginoidea virens, Trichoderma species from section Longibrachiatum) and in Penicillium rubens (Eurotiales). In addition, Colletotrichum graminicola contained the two pks (sor1 and sor2), but not the other sor genes. A. chrysogenum was the evolutionary eldest species in which sor1, sor2, sor3, sor4 and sor6 were present. Sor5 was gained by lateral gene transfer (LGT) from P. rubens. In the younger Hypocreales (U. virens, Trichoderma spp.), the cluster evolved by vertical transfer, but sor2 was lost and regained by LGT from C. graminicola. SorB (=sor2) and sorD (=sor4) were symplesiomorphic in P. rubens, whereas sorA, sorC and sorF were obtained by LGT from A. chrysogenum, and sorE by LGT from Pestalotiopsis fici (Xylariales). The sorbicillinoid gene cluster in Trichoderma section Longibrachiatum is under strong purifying selection. The T. reesei sor genes are expressed during fast vegetative growth, during antagonism of other fungi and regulated by the secondary metabolism regulator LAE1.

Conclusions

Our findings pinpoint the evolution of the fungal sorbicillinoid biosynthesis gene cluster. The core cluster arose in early Hypocreales, and was complemented by LGT. During further speciation in the Hypocreales, it became subject to birth and death evolution in selected lineages. In P. rubrens (Eurotiales), two cluster genes were symplesiomorphic, and the whole cluster formed by LGT from at least two different fungal donors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0834-6) contains supplementary material, which is available to authorized users.