Discovery of CO2 tolerance genes associated with virulence in the fungal pathogen Cryptococcus neoformans.

Cryptococcus neoformans is a ubiquitous soil fungus and airborne pathogen that causes over 180,000 deaths each year. Cryptococcus must adapt to host CO2 levels to cause disease, but the genetic basis for this adaptation is unknown. We utilized quantitative trait loci mapping with 374 progeny from a cross between a CO2-tolerant clinical isolate and a CO2-sensitive environmental isolate to identify genetic regions regulating CO2 tolerance. To identify specific quantitative trait genes, we applied fine mapping through bulk segregant analysis of near-isogenic progeny with distinct tolerance levels to CO2. We found that virulence among near-isogenic strains in a murine model of cryptococcosis correlated with CO2 tolerance. Moreover, we discovered that sensitive strains may adapt in vivo to become more CO2 tolerant and more virulent. These findings highlight the underappreciated role of CO2 tolerance and its importance in the ability of an opportunistic environmental pathogen to cause disease.

The identification of a key gene highlights macrocyclic ring's role in trichothecene toxicity.

Trichothecenes are toxins produced by certain species from several fungal genera, including Aspergillus, Fusarium, Isaria, Paramyrothecium, Stachybotrys, Trichoderma, and Trichothecium. These toxins are of interest because they contribute to the toxigenicity, plant pathogenicity, and/or biological control activities of some fungi. All trichothecenes have the same core (12,13-epoxytrichothec-9-ene or EPT) structure but can differ from one another by the presence or absence of a macrocyclic ring formed from polyketide and isoprenoid substituents esterified to carbon atoms 4 and 15 of EPT, respectively. Genes required for formation and some modifications of EPT have been elucidated, but almost nothing is known about genes specific to the formation of the macrocyclic ring. Therefore, we used genomic, transcriptomic, metabolomic, and gene deletion analyses to identify genes that are required specifically for the formation of the macrocyclic ring. These analyses identified one gene, TRI24, that is predicted to encode an acyltransferase and that is required for macrocyclic ring formation during biosynthesis of macrocyclic trichothecenes by the fungus Paramyrothecium roridum. In addition, a TRI24 deletion mutant of P. roridum caused less severe disease symptoms on common bean and had less antifungal activity than its wild-type progenitor strain. We propose that the reduced aggressiveness and antifungal activity of the mutant resulted from its inability to produce trichothecenes with a macrocyclic ring. To our knowledge, this is the first report of a gene required specifically for the formation of the macrocyclic ring of trichothecenes and that loss of the macrocyclic ring of trichothecenes can alter the biological activities of a fungus. KEY POINTS: • TRI24 gene is found in all known macrocyclic trichothecene-producing fungi. • A tri24-deletion mutant exhibits a reduction in antifungal and plant disease activities. • TRI24 is the first described gene specific to macrocyclic trichothecene biosynthesis.

Virulence of banana wilt-causing fungal pathogen Fusarium oxysporum tropical race 4 is mediated by nitric oxide biosynthesis and accessory genes.

Fusarium wilt of banana, caused by Fusarium oxysporum f. sp. cubense (Foc), is one of the most damaging plant diseases known. Foc race 1 (R1) decimated the Gros Michel-based banana (Musa acuminata) trade, and now Foc tropical race 4 (TR4) threatens global production of its replacement, the Cavendish banana. Here population genomics revealed that all Cavendish banana-infecting Foc race 4 strains share an evolutionary origin distinct from that of R1 strains. Although TR4 lacks accessory chromosomes, it contains accessory genes at the ends of some core chromosomes that are enriched for virulence and mitochondria-related functions. Meta-transcriptomics revealed the unique induction of the entire mitochondrion-localized nitric oxide (NO) biosynthesis pathway upon TR4 infection. Empirically, we confirmed the unique induction of a NO burst in TR4, suggesting that nitrosative pressure may contribute to virulence. Targeted mutagenesis demonstrated the functional importance of fungal NO production and the accessory gene SIX4 as virulence factors.

Screening of pathogenicity-deficient Penicillium italicum mutants established by Agrobacterium tumefaciens-mediated transformation.

Blue mold, caused by Penicillium italicum, is one of the main postharvest diseases of citrus fruits during storage and marketing. The pathogenic mechanism remains largely unclear. To explore the potential pathogenesis-related genes of this pathogen, a T-DNA insertion library of P. italicum PI5 was established via Agrobacterium tumefaciens-mediated transformation (ATMT). The system yielded 200-250 transformants per million conidia, and the transformants were genetically stable after five generations of successive subcultures on hygromycin-free media. 2700 transformants were obtained to generate a T-DNA insertion library of P. italicum. Only a few of the 200 randomly selected mutants exhibited significantly weakened virulence on citrus fruits, with two mutants displaying attenuated sporulation. The T-DNA in the two mutants existed as a single copy. Moreover, the mutant genes PiBla (PITC_048370) and PiFTF1 (PITC_077280) identified may be involved in conidia production by regulating expressions of the key regulatory components for conidiogenesis. These results demonstrated that the ATMT system is useful to obtain mutants of P. italicum for further investigation of the molecular mechanisms of pathogenicity and the obtained two pathogenesis-related genes might be novel loci associated with pathogenesis and conidia production.

High rate of gene family evolution in proximity to the origin of ectomycorrhizal symbiosis in Inocybaceae.

The genomes of ectomycorrhizal (ECM) fungi have a reduced number of genes encoding Carbohydrate-Active EnZymes (CAZymes), expansions in transposable elements (TEs) and small secreted proteins (SSPs) compared with saprotrophs. Fewer genes for specific peptidases and lipases in ECM fungi are also reported. It is unclear whether these changes occur at the shift to the ECM habit or are more gradual throughout the evolution of ECM lineages. We generated a genomic dataset of 20 species in the ECM lineage Inocybaceae and compared them with six saprotrophic species. Inocybaceae genomes have fewer CAZymes, peptidases, lipases, secondary metabolite clusters and SSPs and higher TE content than their saprotrophic relatives. There was an increase in the rate of gene family evolution along the branch with the transition to the ECM lifestyle. This branch had very high rate of evolution in CAZymes and had the largest number of contractions. Other significant changes along this branch included expansions in transporters, transposons-related genes and communication genes such as fungal kinases. There is a high concentration of changes in proximity to the transition to the ECM lifestyle, which correspond to the identified key changes for the gain of this lifestyle.

CLOCI: unveiling cryptic fungal gene clusters with generalized detection.

Gene clusters are genomic loci that contain multiple genes that are functionally and genetically linked. Gene clusters collectively encode diverse functions, including small molecule biosynthesis, nutrient assimilation, metabolite degradation, and production of proteins essential for growth and development. Identifying gene clusters is a powerful tool for small molecule discovery and provides insight into the ecology and evolution of organisms. Current detection algorithms focus on canonical 'core' biosynthetic functions many gene clusters encode, while overlooking uncommon or unknown cluster classes. These overlooked clusters are a potential source of novel natural products and comprise an untold portion of overall gene cluster repertoires. Unbiased, function-agnostic detection algorithms therefore provide an opportunity to reveal novel classes of gene clusters and more precisely define genome organization. We present CLOCI (Co-occurrence Locus and Orthologous Cluster Identifier), an algorithm that identifies gene clusters using multiple proxies of selection for coordinated gene evolution. Our approach generalizes gene cluster detection and gene cluster family circumscription, improves detection of multiple known functional classes, and unveils non-canonical gene clusters. CLOCI is suitable for genome-enabled small molecule mining, and presents an easily tunable approach for delineating gene cluster families and homologous loci.

Gene function characterization in Aspergillus niger using a dual resistance marker transformation system mediated by Agrobacterium tumefaciens.

Aspergillus niger is a well-known workhorse for the industrial production of enzymes and organic acids. This fungus can also cause postharvest diseases in fruits. Although Agrobacterium tumefaciens-mediated transformation (ATMT) based on antibiotic resistance markers has been effectively exploited for inspecting functions of target genes in wild-type fungi, it still needs to be further improved in A. niger. In the present study, we re-examined the ATMT in the wild-type A. niger strains using the hygromycin resistance marker and introduced the nourseothricin resistance gene as a new selection marker for this fungus. Unexpectedly, our results revealed that the ATMT method using the resistance markers in A. niger led to numerous small colonies as false-positive transformants on transformation plates. Using the top agar overlay technique to restrict false positive colonies, a transformation efficiency of 87 ± 18 true transformants could be achieved for 106 conidia. With two different selection markers, we could perform both the deletion and complementation of a target gene in a single wild-type A. niger strain. Our results also indicated that two key regulatory genes (laeA and veA) of the velvet complex are required for A. niger to infect apple fruits. Notably, we demonstrated for the first time that a laeA homologous gene from the citrus postharvest pathogen Penicillium digitatum was able to restore the acidification ability and pathogenicity of the A. niger ΔlaeA mutant. The dual resistance marker ATMT system from our work represents an improved genetic tool for gene function characterization in A. niger.

Ancestral Sequence Reconstruction as a Tool to Detect and Study De Novo Gene Emergence.

New protein-coding genes can evolve from previously noncoding genomic regions through a process known as de novo gene emergence. Evidence suggests that this process has likely occurred throughout evolution and across the tree of life. Yet, confidently identifying de novo emerged genes remains challenging. Ancestral sequence reconstruction is a promising approach for inferring whether a gene has emerged de novo or not, as it allows us to inspect whether a given genomic locus ancestrally harbored protein-coding capacity. However, the use of ancestral sequence reconstruction in the context of de novo emergence is still in its infancy and its capabilities, limitations, and overall potential are largely unknown. Notably, it is difficult to formally evaluate the protein-coding capacity of ancestral sequences, particularly when new gene candidates are short. How well-suited is ancestral sequence reconstruction as a tool for the detection and study of de novo genes? Here, we address this question by designing an ancestral sequence reconstruction workflow incorporating different tools and sets of parameters and by introducing a formal criterion that allows to estimate, within a desired level of confidence, when protein-coding capacity originated at a particular locus. Applying this workflow on ∼2,600 short, annotated budding yeast genes (<1,000 nucleotides), we found that ancestral sequence reconstruction robustly predicts an ancient origin for the most widely conserved genes, which constitute "easy" cases. For less robust cases, we calculated a randomization-based empirical P-value estimating whether the observed conservation between the extant and ancestral reading frame could be attributed to chance. This formal criterion allowed us to pinpoint a branch of origin for most of the less robust cases, identifying 49 genes that can unequivocally be considered de novo originated since the split of the Saccharomyces genus, including 37 Saccharomyces cerevisiae-specific genes. We find that for the remaining equivocal cases we cannot rule out different evolutionary scenarios including rapid evolution, multiple gene losses, or a recent de novo origin. Overall, our findings suggest that ancestral sequence reconstruction is a valuable tool to study de novo gene emergence but should be applied with caution and awareness of its limitations.

Revealing Hidden Genes in Botrytis cinerea: New Insights into Genes Involved in the Biosynthesis of Secondary Metabolites.

Utilizing bioinformatics tools, this study expands our understanding of secondary metabolism in Botrytis cinerea, identifying novel genes within polyketide synthase (PKS), non-ribosomal peptide synthetase (NRPS), sesquiterpene cyclase (STC), diterpene cyclase (DTC), and dimethylallyltryptophan synthase (DMATS) families. These findings enrich the genetic framework associated with B. cinerea's pathogenicity and ecological adaptation, offering insights into uncharted metabolic pathways. Significantly, the discovery of previously unannotated genes provides new molecular targets for developing targeted antifungal strategies, promising to enhance crop protection and advance our understanding of fungal biochemistry. This research not only broadens the scope of known secondary metabolites but also opens avenues for future exploration into B. cinerea's biosynthetic capabilities, potentially leading to novel antifungal compounds. Our work underscores the importance of integrating bioinformatics and genomics for fungal research, paving the way for sustainable agricultural practices by pinpointing precise molecular interventions against B. cinerea. This study sets a foundation for further investigations into the fungus's secondary metabolism, with implications for biotechnology and crop disease management.

Holistic transcriptional responses of Cordyceps militaris to different culture temperatures.

Cordyceps militaris is a medicinal entomopathogenic fungus containing valuable biometabolites for pharmaceutical applications. Its genetic inheritance and environmental factors play a crucial role in the production of biomass enriched with cordycepin. While temperature is a crucial controlled parameter for fungal cultivation, its impacts on growth and metabolite biosynthesis remains poorly characterized. This study aimed to investigate the metabolic responses and cordycepin production of C. militaris strain TBRC6039 under various temperature conditions through transcriptome analysis. Among 9599 expressed genes, 576 genes were significantly differentially expressed at culture temperatures of 15 and 25 °C. The changes in the transcriptional responses induced by these temperatures were found in several metabolisms involved in nutrient assimilation and energy source, including amino acids metabolism (e.g., glycine, serine and threonine metabolism) and lipid metabolism (e.g., biosynthesis of unsaturated fatty acids and steroid biosynthesis). At the lower temperature (15 °C), the biosynthetic pathways of lipids, specifically ergosterol and squalene, were the target for maintaining membrane function by transcriptional upregulation. Our study revealed the responsive mechanisms of C. militaris in acclimatization to temperature conditions that provide an insight on physiological manipulation for the production of metabolites by C. militaris.

Aspergillus flavus pangenome (AflaPan) uncovers novel aflatoxin and secondary metabolite associated gene clusters.

BACKGROUND: Aspergillus flavus is an important agricultural and food safety threat due to its production of carcinogenic aflatoxins. It has high level of genetic diversity that is adapted to various environments. Recently, we reported two reference genomes of A. flavus isolates, AF13 (MAT1-2 and highly aflatoxigenic isolate) and NRRL3357 (MAT1-1 and moderate aflatoxin producer). Where, an insertion of 310 kb in AF13 included an aflatoxin producing gene bZIP transcription factor, named atfC. Observations of significant genomic variants between these isolates of contrasting phenotypes prompted an investigation into variation among other agricultural isolates of A. flavus with the goal of discovering novel genes potentially associated with aflatoxin production regulation. Present study was designed with three main objectives: (1) collection of large number of A. flavus isolates from diverse sources including maize plants and field soils; (2) whole genome sequencing of collected isolates and development of a pangenome; and (3) pangenome-wide association study (Pan-GWAS) to identify novel secondary metabolite cluster genes. RESULTS: Pangenome analysis of 346 A. flavus isolates identified a total of 17,855 unique orthologous gene clusters, with mere 41% (7,315) core genes and 59% (10,540) accessory genes indicating accumulation of high genomic diversity during domestication. 5,994 orthologous gene clusters in accessory genome not annotated in either the A. flavus AF13 or NRRL3357 reference genomes. Pan-genome wide association analysis of the genomic variations identified 391 significant associated pan-genes associated with aflatoxin production. Interestingly, most of the significantly associated pan-genes (94%; 369 associations) belonged to accessory genome indicating that genome expansion has resulted in the incorporation of new genes associated with aflatoxin and other secondary metabolites. CONCLUSION: In summary, this study provides complete pangenome framework for the species of Aspergillus flavus along with associated genes for pathogen survival and aflatoxin production. The large accessory genome indicated large genome diversity in the species A. flavus, however AflaPan is a closed pangenome represents optimum diversity of species A. flavus. Most importantly, the newly identified aflatoxin producing gene clusters will be a new source for seeking aflatoxin mitigation strategies and needs new attention in research.

Development of a telomere vector-based approach to overcome limitations caused by lethal phenotypes in the study of essential genes in Magnaporthe oryzae.

Reverse genetic approaches are common tools in genomics for elucidating gene functions, involving techniques such as gene deletion followed by screening for aberrant phenotypes. If the generation of gene deletion mutants fails, the question arises whether the failure stems from technical issues or because the gene of interest (GOI) is essential, meaning that the deletion causes lethality. In this report, we introduce a novel method for assessing gene essentiality using the phytopathogenic ascomycete Magnaporthe oryzae. The method is based on the observation that telomere vectors are lost in transformants during cultivation without selection pressure. We tested the hypothesis that essential genes can be identified in deletion mutants co-transformed with a telomere vector. The M. oryzae gene MoPKC, described in literature as essential, was chosen as GOI. Using CRISPR/Cas9 technology transformants with deleted GOI were generated and backed up by a telomere vector carrying a copy of the GOI and conferring fenhexamid resistance. Transformants in which the GOI deletion in the genome was not successful lost the telomere vector on media without fenhexamid. In contrast, transformants with confirmed GOI deletion retained the telomere vector even in absence of fenhexamid selection. In the latter case, the maintenance of the telomere indicates that the GOI is essential for the surveillance of the fungi, as it would have been lost otherwise. The method presented here allows to test for essentiality of genes when no mutants can be obtained from gene deletion approaches, thereby expanding the toolbox for studying gene function in ascomycetes.

Selection of stable reference genes for qPCR expression of Colletotrichum lindemuthianum, the bean anthracnose pathogen.

Phaseolus vulgaris L., commonly known as the common bean, is a highly nutritious crop often called the "poor man's meat". However, it is susceptible to various diseases throughout the cropping season, with anthracnose caused by Colletotrichum lindemuthianum being a significant threat that leads to substantial losses. There is still a lack of understanding about the molecular basis of C. lindemuthianum pathogenicity. The first step in understanding this is to identify pathogenicity genes that express more during infection of common beans. A reverse transcription quantitative real-time PCR (qPCR) method can be used for virulence gene expression. However, this approach requires selecting appropriate reference genes to normalize relative gene expression data. Currently, there is no reference gene available for C. lindemuthianum. In this study, we selected eight candidate reference genes from the available genome of C. lindemuthianum to bridge the gap. These genes were ACT (Actin), ß-tub (ß-tubulin), EF (Elongation Factor), Cyt C (Cytochrome C), His H3 (Histone H3), CHS1 (Chitin synthetase), GAPDH (Glyceraldehyde-3-phosphate dehydrogenase) and abfA (Alpha-l-Arabinofuranosidase A). The primers for these candidate reference genes were able to amplify cDNA only from the pathogen, demonstrating their specificity. The qPCR efficiency of the primers ranged from 80% to 103%. We analyzed the stability of gene expression in C. lindemuthianum by exposing the mycelium to nine different stress conditions. We employed algorithms, such as GeNorm, NormFinder, BestKeeper, and RefFinder tools, to identify the most stable gene. The analysis using these tools revealed that EF, GAPDH, and ß-tub most stable genes, while ACT and CHS1 showed relatively low expression stability. A large number of potential effector genes have been identified through bioinformatics analysis in C. lindemuthianum. The stable genes for qPCR (EF and GAPDH) discovered in this study will aid the scientific community in determining the relative expression of C. lindemuthianum effector genes.

The Narrow Footprint of Ancient Balancing Selection Revealed by Heterokaryon Incompatibility Genes in Aspergillus fumigatus.

In fungi, fusion between individuals leads to localized cell death, a phenomenon termed heterokaryon incompatibility. Generally, the genes responsible for this incompatibility are observed to be under balancing selection resulting from negative frequency-dependent selection. Here, we assess this phenomenon in Aspergillus fumigatus, a human pathogenic fungus with a very low level of linkage disequilibrium as well as an extremely high crossover rate. Using complementation of auxotrophic mutations as an assay for hyphal compatibility, we screened sexual progeny for compatibility to identify genes involved in this process, called het genes. In total, 5/148 (3.4%) offspring were compatible with a parent and 166/2,142 (7.7%) sibling pairs were compatible, consistent with several segregating incompatibility loci. Genetic mapping identified five loci, four of which could be fine mapped to individual genes, of which we tested three through heterologous expression, confirming their causal relationship. Consistent with long-term balancing selection, trans-species polymorphisms were apparent across several sister species, as well as equal allele frequencies within A. fumigatus. Surprisingly, a sliding window genome-wide population-level analysis of an independent dataset did not show increased Tajima's D near these loci, in contrast to what is often found surrounding loci under balancing selection. Using available de novo assemblies, we show that these balanced polymorphisms are restricted to several hundred base pairs flanking the coding sequence. In addition to identifying the first het genes in an Aspergillus species, this work highlights the interaction of long-term balancing selection with rapid linkage disequilibrium decay.

A Rapid, Equipment-Free Method for Detecting Avirulence Genes of Pyricularia oryzae Using a Lateral Flow Strip-Based RPA Assay.

Rice blast, caused by Pyricularia oryzae, is one of the most destructive rice diseases worldwide. Using resistant rice varieties is the most cost-effective way to control rice blast. Consequently, it is critical to monitor the distribution frequency of avirulence (Avr) genes in rice planting fields to facilitate the breeding of resistant rice varieties. In this study, we established a rapid recombinase polymerase amplification-lateral flow dipstick (RPA-LFD) detection system for the identification of AvrPik, Avr-Piz-t, and Avr-Pi9. The optimized reaction temperature and duration were 37°C and 20 min, indicating that the reaction system could be initiated by body temperature without relying on any precision instruments. Specificity analysis showed that the primer and probe combinations targeting the three Avr genes exhibited a remarkable specificity at genus-level detection. Under the optimized condition, the lower detected thresholds of AvrPik, Avr-Piz-t, and Avr-Pi9 were 10 fg/µl, 100 fg/µl, and 10 pg/µl, respectively. Notably, the detection sensitivity of the three Avr genes was much higher than that of PCR. In addition, we also successfully detected the presence of AvrPik, Avr-Piz-t, and Avr-Pi9 in the leaf and panicle blast lesions with the RPA-LFD detection system. In particular, the genomic DNA was extracted using the simpler PEG-NaOH rapid extraction method. In summary, we developed an RPA detection system for AvrPik, Avr-Pi9, and Avr-Piz-t, combined with the PEG-NaOH rapid DNA extraction method. The innovative approach achieved rapid, real-time, and accurate detection of the three Avr genes in the field, which is helpful to understand the distribution frequency of the three Avr genes in the field and provide theoretical reference for the scientific layout of resistant rice varieties.

Identification of two genes essential for basidiospore formation during the postmeiotic stages in Pleurotus ostreatus.

A sporeless strain is an important breeding target in the mushroom industry. However, basidiospore production in the oyster mushroom Pleurotus ostreatus has been shown to be impaired by single-gene mutations in only two meiosis-related genes, mer3 and msh4. This study proposed a strategy for identifying the genes essential for basidiospore formation after meiotic division to determine new targets for molecular breeding. RNA-seq analysis was performed to identify P. ostreatus genes that are specifically expressed in the gill tissue of fruiting bodies, where basidiospore formation occurs. Transcriptome data during fruiting development of Coprinopsis cinerea, in which the meiotic steps progress synchronously, were then used to identify genes that are active in the postmeiotic stages. Based on these comparative analyses, five P. ostreatus genes were identified. Plasmids containing expression cassettes for hygromycin B-resistance screening, Cas9, and single-guide RNA targeting each gene were introduced into the protoplasts of dikaryotic strain, PC9×#64, to generate dikaryotic gene disruptants. Among the obtained transformants, three dikaryotic pcl1 disruptants and two cro6c disruptants did not produce basidiospores. Microscopic analyses indicated that spore formation was arrested at particular stages in these gene disruptants. These results indicate that these two genes are essential for mature spore formation in this fungus.

FGCD: a database of fungal gene clusters related to secondary metabolism.

Fungal secondary metabolites are not necessary for growth, but they are important for fungal metabolism and ecology because they provide selective advantages for competition, survival and interactions with the environment. These various metabolites are widely used as medicinal precursors and insecticides. Secondary metabolism genes are commonly arranged in clusters along chromosomes, which allow for the coordinate control of complete pathways. In this study, we created the Fungal Gene Cluster Database to store, retrieve, and visualize secondary metabolite gene cluster information across fungal species. The database was created by merging data from RNA sequencing, Basic Local Alignment Search Tool, genome browser, enrichment analysis and the R Shiny web framework to visualize and query putative gene clusters. This database facilitated the rapid and thorough examination of significant gene clusters across fungal species by detecting, defining and graphically displaying the architecture, organization and expression patterns of secondary metabolite gene clusters. In general, this genomic resource makes use of the tremendous chemical variety of the products of these ecologically and biotechnologically significant gene clusters to our further understanding of fungal secondary metabolism. Database URL: https://www.hebaubioinformatics.cn/FungalGeneCluster/.

Transcriptional study of genes involved in the passage from teliospore to hyphae stage in the fungus Thecaphora frezii, the causal agent of peanut smut.

Peanuts (Arachis hypogaea L.) are among the most important leguminous crops in Argentina. During the growing season, they are frequently attacked by fungal diseases, including Thecaphora frezii. The spores of T. frezii are structures that confer resistance to this phytopathogen. The transition from teliospore to hypha is a characteristic process of some fungi, which is essential for completing their life cycle. Using the transcriptomes of teliospores and hyphae of T. frezii, we aimed to identify genes that were differentially expressed during this transition, and we found 134 up-regulated and 66 down-regulated genes, which would participate in different cellular processes such as: (a) cell cycle and DNA processing; (b) cell fate; (c) rescue, defense and cellular virulence; (d) detoxification by CYP450; (e) energy; (f) nutrient interaction and nutritional adaptation; (g) metabolism; (g) proteins with binding functions or cofactor requirements; (h) stress, cell differentiation and biogenesis of cell components; and (i) transport, cell communication and transcription. The identification of genes in T. frezii and their expression levels during different stages of differentiation could contribute to our understanding of the biological mechanisms in this fungus.

Functional and Evolutionary Integration of a Fungal Gene With a Bacterial Operon.

Siderophores are crucial for iron-scavenging in microorganisms. While many yeasts can uptake siderophores produced by other organisms, they are typically unable to synthesize siderophores themselves. In contrast, Wickerhamiella/Starmerella (W/S) clade yeasts gained the capacity to make the siderophore enterobactin following the remarkable horizontal acquisition of a bacterial operon enabling enterobactin synthesis. Yet, how these yeasts absorb the iron bound by enterobactin remains unresolved. Here, we demonstrate that Enb1 is the key enterobactin importer in the W/S-clade species Starmerella bombicola. Through phylogenomic analyses, we show that ENB1 is present in all W/S clade yeast species that retained the enterobactin biosynthetic genes. Conversely, it is absent in species that lost the ent genes, except for Starmerella stellata, making this species the only cheater in the W/S clade that can utilize enterobactin without producing it. Through phylogenetic analyses, we infer that ENB1 is a fungal gene that likely existed in the W/S clade prior to the acquisition of the ent genes and subsequently experienced multiple gene losses and duplications. Through phylogenetic topology tests, we show that ENB1 likely underwent horizontal gene transfer from an ancient W/S clade yeast to the order Saccharomycetales, which includes the model yeast Saccharomyces cerevisiae, followed by extensive secondary losses. Taken together, these results suggest that the fungal ENB1 and bacterial ent genes were cooperatively integrated into a functional unit within the W/S clade that enabled adaptation to iron-limited environments. This integrated fungal-bacterial circuit and its dynamic evolution determine the extant distribution of yeast enterobactin producers and cheaters.