Evolutionary and functional patterns of shared gene neighbourhood in fungi

The evolution and ecology of psilocybin in nature

Fungi produce diverse metabolites that can have antimicrobial, antifungal, antifeedant, or psychoactive properties. Among these metabolites are the tryptamine-derived compounds psilocybin, its precursors, and natural derivatives (collectively referred to as psiloids), which have played significant roles in human society and culture. The high allocation of nitrogen to psiloids in mushrooms, along with evidence of convergent evolution and horizontal transfer of psilocybin genes, suggest they provide a selective benefit to some fungi. However, no precise ecological roles of psilocybin have been experimentally determined. The structural and functional similarities of psiloids to serotonin, an essential neurotransmitter in animals, suggest that they may enhance the fitness of fungi through interference with serotonergic processes. However, other ecological mechanisms of psiloids have been proposed. Here, we review the literature pertinent to psilocybin ecology and propose potential adaptive advantages psiloids may confer to fungi.

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/.

Description and Genome Characterization of Three Novel Fungal Strains Isolated from Mars 2020 Mission-Associated Spacecraft Assembly Facility Surfaces-Recommendations for Two New Genera and One Species

National Aeronautics and Space Administration’s (NASA) spacecraft assembly facilities are monitored for the presence of any bacteria or fungi that might conceivably survive a transfer to an extraterrestrial environment. Fungi present a broad and diverse range of phenotypic and functional traits to adapt to extreme conditions, hence the detection of fungi and subsequent eradication of them are needed to prevent forward contamination for future NASA missions. During the construction and assembly for the Mars 2020 mission, three fungal strains with unique morphological and phylogenetic properties were isolated from spacecraft assembly facilities. The reconstruction of phylogenetic trees based on several gene loci (ITS, LSU, SSU, RPB, TUB, TEF1) using multi-locus sequence typing (MLST) and whole genome sequencing (WGS) analyses supported the hypothesis that these were novel species. Here we report the genus or species-level classification of these three novel strains via a polyphasic approach using phylogenetic analysis, colony and cell morphology, and comparative analysis of WGS. The strain FJI-L9-BK-P1 isolated from the Jet Propulsion Laboratory Spacecraft Assembly Facility (JPL-SAF) exhibited a putative phylogenetic relationship with the strain Aaosphaeria arxii CBS175.79 but showed distinct morphology and microscopic features. Another JPL-SAF strain, FJII-L3-CM-DR1, was phylogenetically distinct from members of the family Trichomeriaceae and exhibited morphologically different features from the genera Lithohypha and Strelitziana. The strain FKI-L1-BK-DR1 isolated from the Kennedy Space Center facility was identified as a member of Dothideomycetes incertae sedis and is closely related to the family Kirschsteiniotheliaceae according to a phylogenetic analysis. The polyphasic taxonomic approach supported the recommendation for establishing two novel genera and one novel species. The names Aaosphaeria pasadenensis (FJI-L9-BK-P1 = NRRL 64424 = DSM 114621), Pasadenomyces melaninifex (FJII-L3-CM-DR1 = NRRL 64433 = DSM 114623), and Floridaphiala radiotolerans (FKI-L1-BK-DR1 = NRRL 64434 = DSM 114624) are proposed as type species. Furthermore, resistance to ultraviolet-C and presence of specific biosynthetic gene cluster(s) coding for metabolically active compounds are unique to these strains.

Transcriptome and proteome profiling reveals complex adaptations of Candida parapsilosis cells assimilating hydroxyaromatic carbon sources

Many fungal species utilize hydroxyderivatives of benzene and benzoic acid as carbon sources. The yeast Candida parapsilosis metabolizes these compounds via the 3-oxoadipate and gentisate pathways, whose components are encoded by two metabolic gene clusters. In this study, we determine the chromosome level assembly of the C. parapsilosis strain CLIB214 and use it for transcriptomic and proteomic investigation of cells cultivated on hydroxyaromatic substrates. We demonstrate that the genes coding for enzymes and plasma membrane transporters involved in the 3-oxoadipate and gentisate pathways are highly upregulated and their expression is controlled in a substrate-specific manner. However, regulatory proteins involved in this process are not known. Using the knockout mutants, we show that putative transcriptional factors encoded by the genes OTF1 and GTF1 located within these gene clusters function as transcriptional activators of the 3-oxoadipate and gentisate pathway, respectively. We also show that the activation of both pathways is accompanied by upregulation of genes for the enzymes involved in β-oxidation of fatty acids, glyoxylate cycle, amino acid metabolism, and peroxisome biogenesis. Transcriptome and proteome profiles of the cells grown on 4-hydroxybenzoate and 3-hydroxybenzoate, which are metabolized via the 3-oxoadipate and gentisate pathway, respectively, reflect their different connection to central metabolism. Yet we find that the expression profiles differ also in the cells assimilating 4-hydroxybenzoate and hydroquinone, which are both metabolized in the same pathway. This finding is consistent with the phenotype of the Otf1p-lacking mutant, which exhibits impaired growth on hydroxybenzoates, but still utilizes hydroxybenzenes, thus indicating that additional, yet unidentified transcription factor could be involved in the 3-oxoadipate pathway regulation. Moreover, we propose that bicarbonate ions resulting from decarboxylation of hydroxybenzoates also contribute to differences in the cell responses to hydroxybenzoates and hydroxybenzenes. Finally, our phylogenetic analysis highlights evolutionary paths leading to metabolic adaptations of yeast cells assimilating hydroxyaromatic substrates.

Changes in the distribution of fitness effects and adaptive mutational spectra following a single first step towards adaptation

Historical contingency and diminishing returns epistasis have been typically studied for relatively divergent genotypes and/or over long evolutionary timescales. Here, we use Saccharomyces cerevisiae to study the extent of diminishing returns and the changes in the adaptive mutational spectra following a single first adaptive mutational step. We further evolve three clones that arose under identical conditions from a common ancestor. We follow their evolutionary dynamics by lineage tracking and determine adaptive outcomes using fitness assays and whole genome sequencing. We find that diminishing returns manifests as smaller fitness gains during the 2^nd step of adaptation compared to the 1^st step, mainly due to a compressed distribution of fitness effects. We also find that the beneficial mutational spectra for the 2^nd adaptive step are contingent on the 1^st step, as we see both shared and diverging adaptive strategies. Finally, we find that adaptive loss-of-function mutations, such as nonsense and frameshift mutations, are less common in the second step of adaptation than in the first step.

Analyses of both natural and experimental evolution suggest that adaptation depends on the evolutionary past and adaptive potential decreases over time. Here, by tracking yeast adaptation with DNA barcoding, the authors show that such evolutionary phenomena can be observed even after a single adaptive step.

Current Approaches for Advancement in Understanding the Molecular Mechanisms of Mycotoxin Biosynthesis

Filamentous fungi are able to synthesise a remarkable range of secondary metabolites, which play various key roles in the interaction between fungi and the rest of the biosphere, determining their ecological fitness. Many of them can have a beneficial activity to be exploited, as well as negative impact on human and animal health, as in the case of mycotoxins contaminating large quantities of food, feed, and agricultural products worldwide and posing serious health and economic risks. The elucidation of the molecular aspects of mycotoxin biosynthesis has been greatly sped up over the past decade due to the advent of next-generation sequencing technologies, which greatly reduced the cost of genome sequencing and related omic analyses. Here, we briefly highlight the recent progress in the use and integration of omic approaches for the study of mycotoxins biosynthesis. Particular attention has been paid to genomics and transcriptomic approaches for the identification and characterisation of biosynthetic gene clusters of mycotoxins and the understanding of the regulatory pathways activated in response to physiological and environmental factors leading to their production. The latest innovations in genome-editing technology have also provided a more powerful tool for the complete explanation of regulatory and biosynthesis pathways. Finally, we address the crucial issue of the interpretation of the combined omics data on the biology of the mycotoxigenic fungi. They are rapidly expanding and require the development of resources for more efficient integration, as well as the completeness and the availability of intertwined data for the research community.

Transcription Factor Repurposing Offers Insights into Evolution of Biosynthetic Gene Cluster Regulation

The fungal kingdom has provided advances in our ability to identify biosynthetic gene clusters (BGCs) and to examine how gene composition of BGCs evolves across species and genera. However, little is known about the evolution of specific BGC regulators that mediate how BGCs produce secondary metabolites (SMs). A bioinformatics search for conservation of the Aspergillus fumigatus xanthocillin BGC revealed an evolutionary trail of xan-like BGCs across Eurotiales species. Although the critical regulatory and enzymatic genes were conserved in Penicillium expansum, overexpression (OE) of the conserved xan BGC transcription factor (TF) gene, PexanC, failed to activate the putative xan BGC transcription or xanthocillin production in P. expansum, in contrast to the role of AfXanC in A. fumigatus. Surprisingly, OE::PexanC was instead found to promote citrinin synthesis in P. expansum via trans induction of the cit pathway-specific TF, ctnA, as determined by cit BGC expression and chemical profiling of ctnA deletion and OE::PexanC single and double mutants. OE::AfxanC results in significant increases of xan gene expression and metabolite synthesis in A. fumigatus but had no effect on either xanthocillin or citrinin production in P. expansum. Bioinformatics and promoter mutation analysis led to the identification of an AfXanC binding site, 5'-AGTCAGCA-3', in promoter regions of the A. fumigatus xan BGC genes. This motif was not in the ctnA promoter, suggesting a different binding site of PeXanC. A compilation of a bioinformatics examination of XanC orthologs and the presence/absence of the 5'-AGTCAGCA-3' binding motif in xan BGCs in multiple Aspergillus and Penicillium spp. supports an evolutionary divergence of XanC regulatory targets that we speculate reflects an exaptation event in the Eurotiales. IMPORTANCE Fungal secondary metabolites (SMs) are an important source of pharmaceuticals on one hand and toxins on the other. Efforts to identify the biosynthetic gene clusters (BGCs) that synthesize SMs have yielded significant insights into how variation in the genes that compose BGCs may impact subsequent metabolite production within and between species. However, the role of regulatory genes in BGC activation is less well understood. Our finding that the bZIP transcription factor XanC, located in the xanthocillin BGC of both Aspergillus fumigatus and Penicillium expansum, has functionally diverged to regulate different BGCs in these two species emphasizes that the diversification of BGC regulatory elements may sometimes occur through exaptation, which is the co-option of a gene that evolved for one function to a novel function. Furthermore, this work suggests that the loss/gain of transcription factor binding site targets may be an important mediator in the evolution of secondary-metabolism regulatory elements.

Colocality to Cofunctionality: Eukaryotic Gene Neighborhoods as a Resource for Function Discovery

Diverging from the classic paradigm of random gene order in eukaryotes, gene proximity can be leveraged to systematically identify functionally related gene neighborhoods in eukaryotes, utilizing techniques pioneered in bacteria. Current methods of identifying gene neighborhoods typically rely on sequence similarity to characterized gene products. However, this approach is not robust for nonmodel organisms like algae, which are evolutionarily distant from well-characterized model organisms. Here, we utilize a comparative genomic approach to identify evolutionarily conserved proximal orthologous gene pairs conserved across at least two taxonomic classes of green algae. A total of 317 gene neighborhoods were identified. In some cases, gene proximity appears to have been conserved since before the streptophyte–chlorophyte split, 1,000 Ma. Using functional inferences derived from reconstructed evolutionary relationships, we identified several novel functional clusters. A putative mycosporine-like amino acid, “sunscreen,” neighborhood contains genes similar to either vertebrate or cyanobacterial pathways, suggesting a novel mosaic biosynthetic pathway in green algae. One of two putative arsenic-detoxification neighborhoods includes an organoarsenical transporter (ArsJ), a glyceraldehyde 3-phosphate dehydrogenase-like gene, homologs of which are involved in arsenic detoxification in bacteria, and a novel algal-specific phosphoglycerate kinase-like gene. Mutants of the ArsJ-like transporter and phosphoglycerate kinase-like genes in Chlamydomonas reinhardtii were found to be sensitive to arsenate, providing experimental support for the role of these identified neighbors in resistance to arsenate. Potential evolutionary origins of neighborhoods are discussed, and updated annotations for formerly poorly annotated genes are presented, highlighting the potential of this strategy for functional annotation.

Comparative Genomics Used to Predict Virulence Factors and Metabolic Genes among Monilinia Species

Brown rot, caused by Monilinia spp., is among the most important diseases in stone fruits, and some pome fruits (mainly apples). This disease is responsible for significant yield losses, particularly in stone fruits, when weather conditions favorable for disease development appear. To achieve future sustainable strategies to control brown rot on fruit, one potential approach will be to characterize genomic variation among Monilinia spp. to define, among others, the capacity to infect fruit in this genus. In the present work, we performed genomic and phylogenomic comparisons of five Monilinia species and inferred differences in numbers of secreted proteins, including CAZy proteins and other proteins important for virulence. Duplications specific to Monilinia were sparse and, overall, more genes have been lost than gained. Among Monilinia spp., low variability in the CAZome was observed. Interestingly, we identified several secondary metabolism clusters based on similarity to known clusters, and among them was a cluster with homology to pyriculol that could be responsible for the synthesis of chloromonilicin. Furthermore, we compared sequences of all strains available from NCBI of these species to assess their MAT loci and heterokaryon compatibility systems. Our comparative analyses provide the basis for future studies into understanding how these genomic differences underlie common or differential abilities to interact with the host plant.

Ten Years of Collaborative Progress in the Quest for Orthologs

Accurate determination of the evolutionary relationships between genes is a foundational challenge in biology. Homology-evolutionary relatedness-is in many cases readily determined based on sequence similarity analysis. By contrast, whether or not two genes directly descended from a common ancestor by a speciation event (orthologs) or duplication event (paralogs) is more challenging, yet provides critical information on the history of a gene. Since 2009, this task has been the focus of the Quest for Orthologs (QFO) Consortium. The sixth QFO meeting took place in Okazaki, Japan in conjunction with the 67th National Institute for Basic Biology conference. Here, we report recent advances, applications, and oncoming challenges that were discussed during the conference. Steady progress has been made toward standardization and scalability of new and existing tools. A feature of the conference was the presentation of a panel of accessible tools for phylogenetic profiling and several developments to bring orthology beyond the gene unit-from domains to networks. This meeting brought into light several challenges to come: leveraging orthology computations to get the most of the incoming avalanche of genomic data, integrating orthology from domain to biological network levels, building better gene models, and adapting orthology approaches to the broad evolutionary and genomic diversity recognized in different forms of life and viruses.

PCR screening reveals abundance of bovicin-like bacteriocins among ruminal Streptococcus spp. isolated from beef and dairy cattle

AIMS

To investigate the inhibitory activity and the distribution of biosynthetic genes encoding bovicin-like bacteriocins among ruminal Streptococcus isolated from beef and dairy cattle.

METHODS AND RESULTS

Most isolates were classified as Streptococcus equinus and Streptococcus lutetiensis based on 16S rRNA sequencing. The antimicrobial activity of 150 ruminal streptococci isolated from beef and dairy cattle were tested by deferred inhibition assays and their genetic diversity was characterized by BOX-PCR. The frequency of biosynthetic genes associated with the biosynthesis of bovicin-like bacteriocins (bovicin HC5 and bovicin 255) was investigated by PCR screening. Approximately 33% of the ruminal streptococci isolated from Nellore heifers showed inhibitory activity in vitro with the majority harbouring genes for bacteriocin biosynthesis. In contrast, streptococci from Holstein cows showed limited inhibitory activity and a lower frequency of bacteriocin biosynthetic genes.

CONCLUSIONS

Streptococcus from the rumen of beef and dairy cattle exhibit remarkable differences in inhibitory activity and distribution of genes associated with the biosynthesis of prototypical bovicins (bovicin HC5 and bovicin 255).

SIGNIFICANCE AND IMPACT OF THE STUDY

Our findings demonstrate that bovicin HC5 is distributed among ruminal streptococci from different breeds of cattle. The high degree of conservation of the bovicin HC5 structural gene among strains of ruminal streptococci suggests that random genetic drift is not a dominant force in the evolution of this bacteriocin.

Expanding the Knowledge on the Skillful Yeast Cyberlindnera jadinii

Cyberlindnera jadinii is widely used as a source of single-cell protein and is known for its ability to synthesize a great variety of valuable compounds for the food and pharmaceutical industries. Its capacity to produce compounds such as food additives, supplements, and organic acids, among other fine chemicals, has turned it into an attractive microorganism in the biotechnology field. In this review, we performed a robust phylogenetic analysis using the core proteome of C. jadinii and other fungal species, from Asco- to Basidiomycota, to elucidate the evolutionary roots of this species. In addition, we report the evolution of this species nomenclature over-time and the existence of a teleomorph (C. jadinii) and anamorph state (Candida utilis) and summarize the current nomenclature of most common strains. Finally, we highlight relevant traits of its physiology, the solute membrane transporters so far characterized, as well as the molecular tools currently available for its genomic manipulation. The emerging applications of this yeast reinforce its potential in the white biotechnology sector. Nonetheless, it is necessary to expand the knowledge on its metabolism, regulatory networks, and transport mechanisms, as well as to develop more robust genetic manipulation systems and synthetic biology tools to promote the full exploitation of C. jadinii.

Mycena genomes resolve the evolution of fungal bioluminescence

Mushroom-forming fungi in the order Agaricales represent an independent origin of bioluminescence in the tree of life; yet the diversity, evolutionary history, and timing of the origin of fungal luciferases remain elusive. We sequenced the genomes and transcriptomes of five bonnet mushroom species (Mycena spp.), a diverse lineage comprising the majority of bioluminescent fungi. Two species with haploid genome assemblies ∼150 Mb are among the largest in Agaricales, and we found that a variety of repeats between Mycena species were differentially mediated by DNA methylation. We show that bioluminescence evolved in the last common ancestor of mycenoid and the marasmioid clade of Agaricales and was maintained through at least 160 million years of evolution. Analyses of synteny across genomes of bioluminescent species resolved how the luciferase cluster was derived by duplication and translocation, frequently rearranged and lost in most Mycena species, but conserved in the Armillaria lineage. Luciferase cluster members were coexpressed across developmental stages, with the highest expression in fruiting body caps and stipes, suggesting fruiting-related adaptive functions. Our results contribute to understanding a de novo origin of bioluminescence and the corresponding gene cluster in a diverse group of enigmatic fungal species.

Fungal evolution: cellular, genomic and metabolic complexity

The question of how phenotypic and genomic complexity are inter-related and how they are shaped through evolution is a central question in biology that historically has been approached from the perspective of animals and plants. In recent years, however, fungi have emerged as a promising alternative system to address such questions. Key to their ecological success, fungi present a broad and diverse range of phenotypic traits. Fungal cells can adopt many different shapes, often within a single species, providing them with great adaptive potential. Fungal cellular organizations span from unicellular forms to complex, macroscopic multicellularity, with multiple transitions to higher or lower levels of cellular complexity occurring throughout the evolutionary history of fungi. Similarly, fungal genomes are very diverse in their architecture. Deep changes in genome organization can occur very quickly, and these phenomena are known to mediate rapid adaptations to environmental changes. Finally, the biochemical complexity of fungi is huge, particularly with regard to their secondary metabolites, chemical products that mediate many aspects of fungal biology, including ecological interactions. Herein, we explore how the interplay of these cellular, genomic and metabolic traits mediates the emergence of complex phenotypes, and how this complexity is shaped throughout the evolutionary history of Fungi.

Genetic bases for variation in structure and biological activity of trichothecene toxins produced by diverse fungi

Trichothecenes are sesquiterpene toxins produced by diverse but relatively few fungal species in at least three classes of Ascomycetes: Dothideomycetes, Eurotiomycetes, and Sordariomycetes. Approximately 200 structurally distinct trichothecene analogs have been described, but a given fungal species typically produces only a small subset of analogs. All trichothecenes share a core structure consisting of a four-ring nucleus known as 12,13-epoxytrichothec-9-ene. This structure can be substituted at various positions with hydroxyl, acyl, or keto groups to give rise to the diversity of trichothecene structures that has been described. Over the last 30 years, the genetic and biochemical pathways required for trichothecene biosynthesis in several species of the fungi Fusarium and Trichoderma have been elucidated. In addition, phylogenetic and functional analyses of trichothecene biosynthetic (TRI) genes from fungi in multiple genera have provided insights into how acquisition, loss, and changes in functions of TRI genes have given rise to the diversity of trichothecene structures. These analyses also suggest both divergence and convergence of TRI gene function during the evolutionary history of trichothecene biosynthesis. What has driven trichothecene structural diversification remains an unanswered question. However, insight into the role of trichothecenes in plant pathogenesis of Fusarium species and into plant glucosyltransferases that detoxify the toxins by glycosylating them point to a possible driver. Because the glucosyltransferases can have substrate specificity, changes in trichothecene structures produced by a fungus could allow it to evade detoxification by the plant enzymes. Thus, it is possible that advantages conferred by evading detoxification have contributed to trichothecene structural diversification. KEY POINTS : • TRI genes have evolved by diverse processes: loss, acquisition and changes in function. • Some TRI genes have acquired the same function by convergent evolution. • Some other TRI genes have evolved divergently to have different functions. • Some TRI genes were acquired or resulted from diversification in function of other genes. • Substrate specificity of plant glucosyltransferases could drive trichothecene diversity.

Diversity of Secondary Metabolism in Aspergillus nidulans Clinical Isolates

The filamentous fungus Aspergillus nidulans has been a primary workhorse used to understand fungal genetics. Much of this work has focused on elucidating the genetics of biosynthetic gene clusters (BGCs) and the secondary metabolites (SMs) they produce. SMs are both niche defining in fungi and of great economic importance to humans. Despite the focus on A. nidulans, very little is known about the natural diversity in secondary metabolism within this species. We determined the BGC content and looked for evolutionary patterns in BGCs from whole-genome sequences of two clinical isolates and the A4 reference genome of A. nidulans Differences in BGC content were used to explain SM profiles determined using liquid chromatography-high-resolution mass spectrometry. We found that in addition to genetic variation of BGCs contained by all isolates, nine BGCs varied by presence/absence. We discovered the viridicatumtoxin BGC in A. nidulans and suggest that this BGC has undergone a horizontal gene transfer from the Aspergillus section Nigri lineage into Penicillium sometime after the sections Nigri and Nidulantes diverged. We identified the production of viridicatumtoxin and several other compounds previously not known to be produced by A. nidulans One isolate showed a lack of sterigmatocystin production even though it contained an apparently intact sterigmatocystin BGC, raising questions about other genes and processes known to regulate this BGC. Altogether, our work uncovers a large degree of intraspecies diversity in BGC and SM production in this genetic model species and offers new avenues to understand the evolution and regulation of secondary metabolism.IMPORTANCE Much of what we know about the genetics underlying secondary metabolite (SM) production and the function of SMs in the model fungus Aspergillus nidulans comes from a single reference genome. A growing body of research indicates the importance of biosynthetic gene cluster (BGC) and SM diversity within a species. However, there is no information about the natural diversity of secondary metabolism in A. nidulans We discovered six novel clusters that contribute to the considerable variation in both BGC content and SM production within A. nidulans We characterize a diverse set of mutations and emphasize how findings of single nucleotide polymorphisms (SNPs), deletions, and differences in evolutionary history encompass much of the variation observed in nonmodel systems. Our results emphasize that A. nidulans may also be a strong model to use within-species diversity to elucidate regulatory cross talk, fungal ecology, and drug discovery systems.

Evolutionary Histories of Type III Polyketide Synthases in Fungi

Type III polyketide synthases (PKSs) produce secondary metabolites with diverse biological activities, including antimicrobials. While they have been extensively studied in plants and bacteria, only a handful of type III PKSs from fungi has been characterized in the last 15 years. The exploitation of fungal type III PKSs to produce novel bioactive compounds requires understanding the diversity of these enzymes, as well as of their biosynthetic pathways. Here, phylogenetic and reconciliation analyses of 522 type III PKSs from 1,193 fungal genomes revealed complex evolutionary histories with massive gene duplications and losses, explaining their discontinuous distribution in the fungal tree of life. In addition, horizontal gene transfer events from bacteria to fungi and, to a lower extent, between fungi, could be inferred. Ancestral gene duplication events have resulted in the divergence of eight phylogenetic clades. Especially, two clades show ancestral linkage and functional co-evolution between a type III PKS and a reducing PKS genes. Investigation of the occurrence of protein domains in fungal type III PKS predicted gene clusters highlighted the diversity of biosynthetic pathways, likely reflecting a large chemical landscape. Type III PKS genes are most often located next to genes encoding cytochrome P450s, MFS transporters and transcription factors, defining ancestral core gene clusters. This analysis also allowed predicting gene clusters for the characterized fungal type III PKSs and provides working hypotheses for the elucidation of the full biosynthetic pathways. Altogether, our analyses provide the fundamental knowledge to motivate further characterization and exploitation of fungal type III PKS biosynthetic pathways.