Variation Among Biosynthetic Gene Clusters, Secondary Metabolite Profiles, and Cards of Virulence Across Aspergillus Species

Recent developments in Aspergillus fumigatus research: diversity, drugs, and disease

SUMMARYAdvances in modern medical therapies for many previously intractable human diseases have improved patient outcomes. However, successful disease treatment outcomes are often prevented due to invasive fungal infections caused by the environmental mold Aspergillus fumigatus. As contemporary antifungal therapies have not experienced the same robust advances as other medical therapies, defining mechanisms of A. fumigatus disease initiation and progression remains a critical research priority. To this end, the World Health Organization recently identified A. fumigatus as a research priority human fungal pathogen and the Centers for Disease Control has highlighted the emergence of triazole-resistant A. fumigatus isolates. The expansion in the diversity of host populations susceptible to aspergillosis and the complex and dynamic A. fumigatus genotypic and phenotypic diversity call for a reinvigorated assessment of aspergillosis pathobiological and drug-susceptibility mechanisms. Here, we summarize recent advancements in the field and discuss challenges in our understanding of A. fumigatus heterogeneity and its pathogenesis in diverse host populations.

Context matters: assessing the impacts of genomic background and ecology on microbial biosynthetic gene cluster evolution

Encoded within many microbial genomes, biosynthetic gene clusters (BGCs) underlie the synthesis of various secondary metabolites that often mediate ecologically important functions. Several studies and bioinformatics methods developed over the past decade have advanced our understanding of both microbial pangenomes and BGC evolution. In this minireview, we first highlight challenges in broad evolutionary analysis of BGCs, including delineation of BGC boundaries and clustering of BGCs across genomes. We further summarize key findings from microbial comparative genomics studies on BGC conservation across taxa and habitats and discuss the potential fitness effects of BGCs in different settings. Afterward, recent research showing the importance of genomic context on the production of secondary metabolites and the evolution of BGCs is highlighted. These studies draw parallels to recent, broader, investigations on gene-to-gene associations within microbial pangenomes. Finally, we describe mechanisms by which microbial pangenomes and BGCs evolve, ranging from the acquisition or origination of entire BGCs to micro-evolutionary trends of individual biosynthetic genes. An outlook on how expansions in the biosynthetic capabilities of some taxa might support theories that open pangenomes are the result of adaptive evolution is also discussed. We conclude with remarks about how future work leveraging longitudinal metagenomics across diverse ecosystems is likely to significantly improve our understanding on the evolution of microbial genomes and BGCs.

Gene regulatory network resource aids in predicting trans-acting regulators of biosynthetic gene clusters in Aspergillus fumigatus

The field of secondary metabolism has greatly benefitted from computational advances in recent years. This has been particularly true for fungal natural product studies. Strides in genome mining have led to the identification of an extraordinary number of secondary metabolite biosynthetic gene clusters (BGCs) across the fungal Kingdom and metabologenomic platforms can group BGCs into gene cluster families and link them to initial chemical structures. Missing are computational applications focused on identifying BGC regulatory networks. Here, we applied the new online gene regulatory network resource, GRAsp (Gene Regulation of Aspergillus fumigatus), to identify unknown and unpredictable BGC trans-acting transcriptional/metabolite production modules. GRAsp correctly predicted a two-component regulatory module composed of the transcription factors (TFs), RogA (regulation of gliotoxin) and HsfA, which negatively regulate the gliotoxin BGC and are also involved in gliotoxin self-protection. RogA functions through the repression of gliZ, the pathway-specific gliotoxin TF, and HsfA functions by activating rogA expression. Furthermore, GRAsp identified TFs that regulate the production of two BGCs lacking pathway-specific TFs, the helvolic acid and fumitremorgin BGCs, respectively. Finally, the known TF, NsdD, was predicted and found to regulate the hexadehydroastechrome BGC. These advances highlight the power of inference algorithms to uncover unpredictable networks in specialized metabolite synthesis.IMPORTANCEToxic secondary metabolites are virulence factors of the opportunistic fungal pathogen Aspergillus fumigatus, yet the transcriptional networks regulating secondary metabolite production remain elusive. Uncovering novel regulators without any prior information is challenging. Computational programs have gained prominence in the field of secondary metabolite research due to their accuracy and ability to handle vast amounts of data, including DNA, RNA, and protein data. In this study, a newly developed online computer platform, Gene Regulation of A. fumigatus, was used to identify five regulators involved in the production of several A. fumigatus toxins, including gliotoxin, helvolic acid, fumitremorgin, and hexadehydroastechrome. This work illustrates the potential for discovering new trans-acting regulators and mechanisms of secondary metabolite regulation through the examination of computational gene regulatory networks.

Fungal Metabolomics: A Comprehensive Approach to Understanding Pathogenesis in Humans and Identifying Potential Therapeutics

Metabolomics has emerged as a transformative tool in the study of microbes, including pathogenic fungi, facilitating the identification of unique metabolic profiles that elucidate their pathogenic mechanisms, host interactions, and treatment resistance. This review highlights key applications of metabolomics in understanding fungal metabolites essential for human virulence, such as mycotoxins produced by various fungal species, including Aspergillus fumigatus (gliotoxin, fumagillins) and Candida species (phenylethyl alcohol, TCA cycle metabolites), and secondary metabolites that contribute to pathogenicity. It also explores the metabolic adaptations of fungi in relation to drug resistance and biofilm formation, revealing alterations in key metabolic pathways during infection, as seen in C. albicans and C. auris. Furthermore, metabolomics aids in deciphering host-pathogen interactions, showcasing how fungi like Cryptococcus neoformans and Candida modify host metabolism to promote survival and evade immune responses. The study of antifungal resistance mechanisms has also benefited from metabolomic approaches, identifying specific metabolite patterns that signify resistance, such as in Candida albicans and Candidozyma (Candida) auris, and informing new therapeutic strategies. The integration of metabolomics with other omics technologies is paving the way for a comprehensive understanding of fungal biology and pathogenesis. Such multi-omics approaches are crucial for discovering new therapeutic targets and developing innovative antifungal treatments. Thus, the purpose of this review is to provide an overview of how metabolomics is revolutionizing our understanding of fungal pathogenesis, drug resistance, and host interactions, and to highlight its potential for identifying new therapeutic targets and improving antifungal strategies.

Aspergillus fumigatus conidial surface-associated proteome reveals factors for fungal evasion and host immunity modulation

Aspergillus fumigatus causes aspergillosis and relies on asexual spores (conidia) for initiating host infection. There is scarce information about A. fumigatus proteins involved in fungal evasion and host immunity modulation. Here we analysed the conidial surface proteome of A. fumigatus, two closely related non-pathogenic species, Aspergillus fischeri and Aspergillus oerlinghausenensis, as well as pathogenic Aspergillus lentulus, to identify such proteins. After identifying 62 proteins exclusively detected on the A. fumigatus conidial surface, we assessed null mutants for 42 genes encoding these proteins. Deletion of 33 of these genes altered susceptibility to macrophage, epithelial cells and cytokine production. Notably, a gene that encodes a putative glycosylasparaginase, modulating levels of the host proinflammatory cytokine IL-1β, is important for infection in an immunocompetent murine model of fungal disease. These results suggest that A. fumigatus conidial surface proteins are important for evasion and modulation of the immune response at the onset of fungal infection.

Evolutionary origin and population diversity of a cryptic hybrid pathogen

Cryptic fungal pathogens pose disease management challenges due to their morphological resemblance to known pathogens. Here, we investigated the genomes and phenotypes of 53 globally distributed isolates of Aspergillus section Nidulantes fungi and found 30 clinical isolates-including four isolated from COVID-19 patients-were A. latus, a cryptic pathogen that originated via allodiploid hybridization. Notably, all A. latus isolates were misidentified. A. latus hybrids likely originated via a single hybridization event during the Miocene and harbor substantial genetic diversity. Transcriptome profiling of a clinical isolate revealed that both parental subgenomes are actively expressed and respond to environmental stimuli. Characterizing infection-relevant traits-such as drug resistance and growth under oxidative stress-revealed distinct phenotypic profiles among A. latus hybrids compared to parental and closely related species. Moreover, we identified four features that could aid A. latus taxonomic identification. Together, these findings deepen our understanding of the origin of cryptic pathogens.

Strain heterogeneity in a non-pathogenic Aspergillus fungus highlights factors associated with virulence

Fungal pathogens exhibit extensive strain heterogeneity, including variation in virulence. Whether closely related non-pathogenic species also exhibit strain heterogeneity remains unknown. Here, we comprehensively characterized the pathogenic potentials (i.e., the ability to cause morbidity and mortality) of 16 diverse strains of Aspergillus fischeri, a non-pathogenic close relative of the major pathogen Aspergillus fumigatus. In vitro immune response assays and in vivo virulence assays using a mouse model of pulmonary aspergillosis showed that A. fischeri strains varied widely in their pathogenic potential. Furthermore, pangenome analyses suggest that A. fischeri genomic and phenotypic diversity is even greater. Genomic, transcriptomic, and metabolic profiling identified several pathways and secondary metabolites associated with variation in virulence. Notably, strain virulence was associated with the simultaneous presence of the secondary metabolites hexadehydroastechrome and gliotoxin. We submit that examining the pathogenic potentials of non-pathogenic close relatives is key for understanding the origins of fungal pathogenicity.

The airway mycobiome and interactions with immunity in health and chronic lung disease

The existence of commensal fungi that reside within the respiratory tract, termed the airway mycobiome, has only recently been discovered. Studies are beginning to characterize the spectrum of fungi that inhabit the human upper and lower respiratory tract but heterogeneous sampling and analysis techniques have limited the generalizability of findings to date. In this review, we discuss existing studies that have examined the respiratory mycobiota in healthy individuals and in those with inflammatory lung conditions such as asthma, chronic obstructive pulmonary disease and cystic fibrosis. Associations between specific fungi and features of disease pathogenesis are emerging but the precise functional consequences imparted by mycobiota upon the immune system remain poorly understood. It is imperative that further research is conducted in this important area as a more detailed understanding could facilitate the development of novel approaches to manipulating the mycobiome for therapeutic benefit.

Genome-wide patterns of noncoding and protein-coding sequence variation in the major fungal pathogen Aspergillus fumigatus

Aspergillus fumigatus is a deadly fungal pathogen, responsible for >400,000 infections/year and high mortality rates. A. fumigatus strains exhibit variation in infection-relevant traits, including in their virulence. However, most A. fumigatus protein-coding genes, including those that modulate its virulence, are shared between A. fumigatus strains and closely related nonpathogenic relatives. We hypothesized that A. fumigatus genes exhibit substantial genetic variation in the noncoding regions immediately upstream to the start codons of genes, which could reflect differences in gene regulation between strains. To begin testing this hypothesis, we identified 5,812 single-copy orthologs across the genomes of 263 A. fumigatus strains. In general, A. fumigatus noncoding regions showed higher levels of sequence variation compared with their corresponding protein-coding regions. Focusing on 2,482 genes whose protein-coding sequence identity scores ranged between 75 and 99%, we identified 478 total genes with signatures of positive selection only in their noncoding regions and 65 total genes with signatures only in their protein-coding regions. Twenty-eight of the 478 noncoding regions and 5 of the 65 protein-coding regions under selection are associated with genes known to modulate A. fumigatus virulence. Noncoding region variation between A. fumigatus strains included single-nucleotide polymorphisms and insertions or deletions of at least a few nucleotides. These results show that noncoding regions of A. fumigatus genes harbor greater sequence variation than protein-coding regions, raising the hypothesis that this variation may contribute to A. fumigatus phenotypic heterogeneity.

Screening and genetic engineering of marine-derived Aspergillus terreus for high-efficient production of lovastatin

BACKGROUND

Lovastatin has widespread applications thanks to its multiple pharmacological effects. Fermentation by filamentous fungi represents the major way of lovastatin production. However, the current lovastatin productivity by fungal fermentation is limited and needs to be improved.

RESULTS

In this study, the lovastatin-producing strains of Aspergillus terreus from marine environment were screened, and their lovastatin productions were further improved by genetic engineering. Five strains of A. terreus were isolated from various marine environments. Their secondary metabolites were profiled by metabolomics analysis using Ultra Performance Liquid Chromatography-Mass spectrometry (UPLC-MS) with Global Natural Products Social Molecular Networking (GNPS), revealing that the production of secondary metabolites was variable among different strains. Remarkably, the strain of A. terreus MJ106 could principally biosynthesize the target drug lovastatin, which was confirmed by High Performance Liquid Chromatography (HPLC) and gene expression analysis. By one-factor experiment, lactose was found to be the best carbon source for A. terreus MJ106 to produce lovastatin. To improve the lovastatin titer in A. terreus MJ106, genetic engineering was applied to this strain. Firstly, a series of strong promoters was identified by transcriptomic and green fluorescent protein reporter analysis. Then, three selected strong promoters were used to overexpress the transcription factor gene lovE encoding the major transactivator for lov gene cluster expression. The results revealed that compared to A. terreus MJ106, all lovE over-expression mutants exhibited significantly more production of lovastatin and higher gene expression. One of them, LovE-b19, showed the highest lovastatin productivity at a titer of 1512 mg/L, which represents the highest production level reported in A. terreus.

CONCLUSION

Our data suggested that combination of strain screen and genetic engineering represents a powerful tool for improving the productivity of fungal secondary metabolites, which could be adopted for large-scale production of lovastatin in marine-derived A. terreus.

The evolution of the gliotoxin biosynthetic gene cluster in Penicillium fungi

Fungi biosynthesize diverse secondary metabolites, small organic bioactive molecules with key roles in fungal ecology. Fungal secondary metabolites are often encoded by physically clustered genes known as biosynthetic gene clusters (BGCs). Fungi in the genus Penicillium produce a cadre of secondary metabolites, some of which are useful (e.g. the antibiotic penicillin and the cholesterol-lowering drug mevastatin) and others harmful (e.g. the mycotoxin patulin and the immunosuppressant gliotoxin) to human affairs. Fungal genomes often also encode resistance genes that confer protection against toxic secondary metabolites. Some Penicillium species, such as Penicillium decumbens, are known to produce gliotoxin, a secondary metabolite with known immunosuppressant activity. To investigate the evolutionary conservation of homologs of the gliotoxin BGC and of genes involved in gliotoxin resistance in Penicillium, we analyzed 35 Penicillium genomes from 23 species. Homologous, lesser fragmented gliotoxin BGCs were found in 12 genomes, mostly fragmented remnants of the gliotoxin BGC were found in 21 genomes, whereas the remaining 2 Penicillium genomes lacked the gliotoxin BGC altogether. In contrast, broad conservation of homologs of resistance genes that reside outside the BGC across Penicillium genomes was observed. Evolutionary rate analysis revealed that BGCs with higher numbers of genes evolve slower than BGCs with few genes, suggestive of constraint and potential functional significance or more recent decay. Gene tree-species tree reconciliation analyses suggested that the history of homologs in the gliotoxin BGC across the genus Penicillium likely involved multiple duplications, losses, and horizontal gene transfers. Our analyses suggest that genes encoded in BGCs can have complex evolutionary histories and be retained in genomes long after the loss of secondary metabolite biosynthesis.

NmrB (AN9181) expression is activated under oxidative stress conditions acting as a metabolic repressor of Aspergillus nidulans

Aspergilli comprise a diversity of species that have been extensively studied due to their catabolic diversity, biotechnological and ecological value, and pathogenicity. An impressive level of structural and functional conservation has been shown for aspergilli, regardless of many (yet) cryptic genomic elements. We have hypothesized the existence of conserved genes responsive to stress in aspergilli. To test the hypothesis of such conserved stress regulators in aspergilli, a straightforward computational strategy integrating well-established bioinformatic tools was used as the starting point. Specifically, five transcriptome-based datasets on exposure to organic compounds were used, covering three distinct Aspergillus species. Among the identified up-regulated genes, only one gene showed the same response in all conditions, AN9181. This gene encodes a protein containing a phenylcoumaran benzylic ether reductase-like domain and a Nitrogen metabolite repressor regulator domain (NmrA). Deletion of this gene caused significant phenotypic alterations compared to that of the parental strain across diverse conditions. Specifically, the deletion of AN9181 raised the mutant's metabolic activity in different nitrogen sources. The acquired data supports that AN9181 acts by repressing (slowing down) A. nidulans growth when exposed to aromatic compounds in a concentration dependent manner. The same phenotype was observed for amphotericin B. Finally, AN9181 underwent differential upregulation under oxidative stress conditions. Collectively, the data suggest that AN9181, herein assigned as NmrB (Nitrogen Metabolite Repression Regulator B), builds up the genetic machinery of perception of oxidative stress by negatively regulating growth under such conditions.

Phylogenomics reveals extensive misidentification of fungal strains from the genus Aspergillus

Modern taxonomic classification is often based on phylogenetic analyses of a few molecular markers, although single-gene studies are still common. Here, we leverage genome-scale molecular phylogenetics (phylogenomics) of species and populations to reconstruct evolutionary relationships in a dense data set of 710 fungal genomes from the biomedically and technologically important genus Aspergillus. To do so, we generated a novel set of 1,362 high-quality molecular markers specific for Aspergillus and provided profile Hidden Markov Models for each, facilitating their use by others. Examining the resulting phylogeny helped resolve ongoing taxonomic controversies, identified new ones, and revealed extensive strain misidentification (7.59% of strains were previously misidentified), underscoring the importance of population-level sampling in species classification. These findings were corroborated using the current standard, taxonomically informative loci. These findings suggest that phylogenomics of species and populations can facilitate accurate taxonomic classifications and reconstructions of the Tree of Life.IMPORTANCEIdentification of fungal species relies on the use of molecular markers. Advances in genomic technologies have made it possible to sequence the genome of any fungal strain, making it possible to use genomic data for the accurate assignment of strains to fungal species (and for the discovery of new ones). We examined the usefulness and current limitations of genomic data using a large data set of 710 publicly available genomes from multiple strains and species of the biomedically, agriculturally, and industrially important genus Aspergillus. Our evolutionary genomic analyses revealed that nearly 8% of publicly available Aspergillus genomes are misidentified. Our work highlights the usefulness of genomic data for fungal systematic biology and suggests that systematic genome sequencing of multiple strains, including reference strains (e.g., type strains), of fungal species will be required to reduce misidentification errors in public databases.

A Timeline of Biosynthetic Gene Cluster Discovery in Aspergillus fumigatus: From Characterization to Future Perspectives

In 1999, the first biosynthetic gene cluster (BGC), synthesizing the virulence factor DHN melanin, was characterized in Aspergillus fumigatus. Since then, 19 additional BGCs have been linked to specific secondary metabolites (SMs) in this species. Here, we provide a comprehensive timeline of A. fumigatus BGC discovery and find that initial advances centered around the commonly expressed SMs where chemical structure informed rationale identification of the producing BGC (e.g., gliotoxin, fumigaclavine, fumitremorgin, pseurotin A, helvolic acid, fumiquinazoline). Further advances followed the transcriptional profiling of a ΔlaeA mutant, which aided in the identification of endocrocin, fumagillin, hexadehydroastechrome, trypacidin, and fumisoquin BGCs. These SMs and their precursors are the commonly produced metabolites in most A. fumigatus studies. Characterization of other BGC/SM pairs required additional efforts, such as induction treatments, including co-culture with bacteria (fumicycline/neosartoricin, fumigermin) or growth under copper starvation (fumivaline, fumicicolin). Finally, four BGC/SM pairs were discovered via overexpression technologies, including the use of heterologous hosts (fumicycline/neosartoricin, fumihopaside, sphingofungin, and sartorypyrone). Initial analysis of the two most studied A. fumigatus isolates, Af293 and A1160, suggested that both harbored ca. 34-36 BGCs. However, an examination of 264 available genomes of A. fumigatus shows up to 20 additional BGCs, with some strains showing considerable variations in BGC number and composition. These new BGCs present a new frontier in the future of secondary metabolism characterization in this important species.

Strain heterogeneity in a non-pathogenic fungus highlights factors contributing to virulence

Fungal pathogens exhibit extensive strain heterogeneity, including variation in virulence. Whether closely related non-pathogenic species also exhibit strain heterogeneity remains unknown. Here, we comprehensively characterized the pathogenic potentials (i.e., the ability to cause morbidity and mortality) of 16 diverse strains of Aspergillus fischeri, a non-pathogenic close relative of the major pathogen Aspergillus fumigatus. In vitro immune response assays and in vivo virulence assays using a mouse model of pulmonary aspergillosis showed that A. fischeri strains varied widely in their pathogenic potential. Furthermore, pangenome analyses suggest that A. fischeri genomic and phenotypic diversity is even greater. Genomic, transcriptomic, and metabolomic profiling identified several pathways and secondary metabolites associated with variation in virulence. Notably, strain virulence was associated with the simultaneous presence of the secondary metabolites hexadehydroastechrome and gliotoxin. We submit that examining the pathogenic potentials of non-pathogenic close relatives is key for understanding the origins of fungal pathogenicity.

Predicting fungal secondary metabolite activity from biosynthetic gene cluster data using machine learning

Fungal secondary metabolites (SMs) contribute to the diversity of fungal ecological communities, niches, and lifestyles. Many fungal SMs have one or more medically and industrially important activities (e.g., antifungal, antibacterial, and antitumor). The genes necessary for fungal SM biosynthesis are typically located right next to each other in the genome and are known as biosynthetic gene clusters (BGCs). However, whether fungal SM bioactivity can be predicted from specific attributes of genes in BGCs remains an open question. We adapted machine learning models that predicted SM bioactivity from bacterial BGC data with accuracies as high as 80% to fungal BGC data. We trained our models to predict the antibacterial, antifungal, and cytotoxic/antitumor bioactivity of fungal SMs on two data sets: (i) fungal BGCs (data set comprised of 314 BGCs) and (ii) fungal (314 BGCs) and bacterial BGCs (1,003 BGCs). We found that models trained on fungal BGCs had balanced accuracies between 51% and 68%, whereas training on bacterial and fungal BGCs had balanced accuracies between 56% and 68%. The low prediction accuracy of fungal SM bioactivities likely stems from the small size of the data set; this lack of data, coupled with our finding that including bacterial BGC data in the training data did not substantially change accuracies currently limits the application of machine learning approaches to fungal SM studies. With >15,000 characterized fungal SMs, millions of putative BGCs in fungal genomes, and increased demand for novel drugs, efforts that systematically link fungal SM bioactivity to BGCs are urgently needed.IMPORTANCEFungi are key sources of natural products and iconic drugs, including penicillin and statins. DNA sequencing has revealed that there are likely millions of biosynthetic pathways in fungal genomes, but the chemical structures and bioactivities of >99% of natural products produced by these pathways remain unknown. We used artificial intelligence to predict the bioactivities of diverse fungal biosynthetic pathways. We found that the accuracies of our predictions were generally low, between 51% and 68%, likely because the natural products and bioactivities of only very few fungal pathways are known. With >15,000 characterized fungal natural products, millions of putative biosynthetic pathways present in fungal genomes, and increased demand for novel drugs, our study suggests that there is an urgent need for efforts that systematically identify fungal biosynthetic pathways, their natural products, and their bioactivities.

Genome-wide patterns of non-coding sequence variation in the major fungal pathogen Aspergillus fumigatus

A.fumigatus is a deadly fungal pathogen, responsible for >400,000 infections/year and high mortality rates. A. fumigatus strains exhibit variation in infection-relevant traits, including in their virulence. However, most A. fumigatus protein-coding genes, including those that modulate its virulence, are shared between A. fumigatus strains and closely related non-pathogenic relatives. We hypothesized that A. fumigatus genes exhibit substantial genetic variation in the non-coding regions immediately upstream to the start codons of genes, which could reflect differences in gene regulation between strains. To begin testing this hypothesis, we identified 5,812 single-copy orthologs across the genomes of 263 A. fumigatus strains. A. fumigatus non-coding regions showed higher levels of sequence variation compared to their corresponding protein-coding regions. Specifically, we found that 1,274 non-coding regions exhibited <75% nucleotide sequence similarity (compared to 928 protein-coding regions) and 3,721 non-coding regions exhibited between 75% and 99% similarity (compared to 2,482 protein-coding regions) across strains. Only 817 non-coding regions exhibited ≥99% sequence similarity compared to 2,402 protein-coding regions. By examining 2,482 genes whose protein-coding sequence identity scores ranged between 75% and 99%, we identified 478 total genes with signatures of positive selection only in their non-coding regions and 65 total genes with signatures only in their protein-coding regions. 28 of the 478 non-coding regions and 5 of the 65 protein-coding regions under selection are associated with genes known to modulate A. fumigatus virulence. Non-coding region variation between A. fumigatus strains included single nucleotide polymorphisms and insertions or deletions of at least a few nucleotides. These results show that non-coding regions of A. fumigatus genes harbor greater sequence variation than protein-coding regions, raising the hypothesis that this variation may contribute to A. fumigatus phenotypic heterogeneity.

Aspergillus fumigatus mitogen-activated protein kinase MpkA is involved in gliotoxin production and self-protection

Aspergillus fumigatus is a saprophytic fungus that can cause a variety of human diseases known as aspergillosis. Mycotoxin gliotoxin (GT) production is important for its virulence and must be tightly regulated to avoid excess production and toxicity to the fungus. GT self-protection by GliT oxidoreductase and GtmA methyltransferase activities is related to the subcellular localization of these enzymes and how GT can be sequestered from the cytoplasm to avoid increased cell damage. Here, we show that GliT:GFP and GtmA:GFP are localized in the cytoplasm and in vacuoles during GT production. The Mitogen-Activated Protein kinase MpkA is essential for GT production and self-protection, interacts physically with GliT and GtmA and it is necessary for their regulation and subsequent presence in the vacuoles. The sensor histidine kinase SlnASln1 is important for modulation of MpkA phosphorylation. Our work emphasizes the importance of MpkA and compartmentalization of cellular events for GT production and self-defense.

Use of a human small airway epithelial cell line to study the interactions of Aspergillus fumigatus with pulmonary epithelial cells

During the initiation of invasive aspergillosis, inhaled Aspergillus fumigatus conidia are deposited on the epithelial cells lining the bronchi, terminal bronchioles, and alveoli. While the interactions of A. fumigatus with bronchial and type II alveolar cell lines have been investigated in vitro, little is known about the interactions of this fungus with terminal bronchiolar epithelial cells. Using the HSAEC1-KT human small airway epithelial (HSAE) cell line, we developed an in vitro model to study the interaction of two strains of A. fumigatus with these cells. We then compared the interactions of A. fumigatus with the A549 type II alveolar epithelial cell line and the HSAE cell line. We found that A. fumigatus conidia were poorly endocytosed by A549 cells, but avidly endocytosed by HSAE cells. A. fumigatus germlings invaded both cell types by induced endocytosis, but not by active penetration. A549 cell endocytosis of A. fumigatus was independent of fungal viability, more dependent on host microfilaments than microtubules, and induced by A. fumigatus CalA interacting with host cell integrin α5β1. By contrast, HSAE cell endocytosis required fungal viability, was more dependent on microtubules than microfilaments, and did not require CalA or integrin α5β1. HSAE cells were more susceptible than A549 cells to damage caused by direct contact with killed A. fumigatus germlings and by secreted fungal products. In response to A. fumigatus infection, A549 cells secreted a broader profile of cytokines and chemokines than HSAE cells. Taken together, these results demonstrate that studies of HSAE cells provide complementary data to A549 cells and thus represent a useful model for probing the interactions of A. fumigatus with bronchiolar epithelial cells in vitro. Importance During the initiation of invasive aspergillosis, Aspergillus fumigatus interacts with the epithelial cells that line the airways and alveoli. Previous studies of A. fumigatus-epithelial cell interactions in vitro used either large airway epithelial cell lines or the A549 type II alveolar epithelial cell line; the interactions of fungi with terminal bronchiolar epithelial cells were not investigated. Using the TERT-immortalized human small airway epithelial HSAEC1-KT (HSAE) cell line, we developed an in vitro model of the interactions of A. fumigatus with bronchiolar epithelial cells. We discovered that A. fumigatus invades and damages A549 and HSAE cell lines by distinct mechanisms. Also, the proinflammatory responses of the cell lines to A. fumigatus are different. These results provide insight into how A. fumigatus interacts with different types of epithelial cells during invasive aspergillosis and demonstrate that HSAE cells are useful in vitro model for investigating the interactions of this fungus with bronchiolar epithelial cells.

Fungal and mycotoxin contaminants in cannabis and hemp flowers: implications for consumer health and directions for further research

Medicinal and recreational uses of Cannabis sativa, commonly known as cannabis or hemp, has increased following its legalization in certain regions of the world. Cannabis and hemp plants interact with a community of microbes (i.e., the phytobiome), which can influence various aspects of the host plant. The fungal composition of the C. sativa phytobiome (i.e., mycobiome) currently consists of over 100 species of fungi, which includes phytopathogens, epiphytes, and endophytes, This mycobiome has often been understudied in research aimed at evaluating the safety of cannabis products for humans. Medical research has historically focused instead on substance use and medicinal uses of the plant. Because several components of the mycobiome are reported to produce toxic secondary metabolites (i.e., mycotoxins) that can potentially affect the health of humans and animals and initiate opportunistic infections in immunocompromised patients, there is a need to determine the potential health risks that these contaminants could pose for consumers. This review discusses the mycobiome of cannabis and hemp flowers with a focus on plant-infecting and toxigenic fungi that are most commonly found and are of potential concern (e.g., Aspergillus, Penicillium, Fusarium, and Mucor spp.). We review current regulations for molds and mycotoxins worldwide and review assessment methods including culture-based assays, liquid chromatography, immuno-based technologies, and emerging technologies for these contaminants. We also discuss approaches to reduce fungal contaminants on cannabis and hemp and identify future research needs for contaminant detection, data dissemination, and management approaches. These approaches are designed to yield safer products for all consumers.

A heterologous expression platform in Aspergillus nidulans for the elucidation of cryptic secondary metabolism biosynthetic gene clusters: discovery of the Aspergillus fumigatus sartorypyrone biosynthetic pathway

Aspergillus fumigatus is a serious human pathogen causing life-threatening Aspergillosis in immunocompromised patients. Secondary metabolites (SMs) play an important role in pathogenesis, but the products of many SM biosynthetic gene clusters (BGCs) remain unknown. In this study, we have developed a heterologous expression platform in Aspergillus nidulans, using a newly created genetic dereplication strain, to express a previously unknown BGC from A. fumigatus and determine its products. The BGC produces sartorypyrones, and we have named it the spy BGC. Analysis of targeted gene deletions by HRESIMS, NMR, and microcrystal electron diffraction (MicroED) enabled us to identify 12 products from the spy BGC. Seven of the compounds have not been isolated previously. We also individually expressed the polyketide synthase (PKS) gene spyA and demonstrated that it produces the polyketide triacetic acid lactone (TAL), a potentially important biorenewable platform chemical. Our data have allowed us to propose a biosynthetic pathway for sartorypyrones and related natural products. This work highlights the potential of using the A. nidulans heterologous expression platform to uncover cryptic BGCs from A. fumigatus and other species, despite the complexity of their secondary metabolomes.

A phylogenetic approach to explore the Aspergillus fumigatus conidial surface-associated proteome and its role in pathogenesis

Aspergillus fumigatus, an important pulmonary fungal pathogen causing several diseases collectively called aspergillosis, relies on asexual spores (conidia) for initiating host infection. Here, we used a phylogenomic approach to compare proteins in the conidial surface of A. fumigatus, two closely related non-pathogenic species, Aspergillus fischeri and Aspergillus oerlinghausenensis, and the cryptic pathogen Aspergillus lentulus. After identifying 62 proteins uniquely expressed on the A. fumigatus conidial surface, we assessed null mutants for 42 genes encoding conidial proteins. Deletion of 33 of these genes altered susceptibility to macrophage killing, penetration and damage to epithelial cells, and cytokine production. Notably, a gene that encodes glycosylasparaginase, which modulates levels of the host pro-inflammatory cytokine IL-1β, is important for infection in an immunocompetent murine model of fungal disease. These results suggest that A. fumigatus conidial surface proteins and effectors are important for evasion and modulation of the immune response at the onset of fungal infection.

Predicting fungal secondary metabolite activity from biosynthetic gene cluster data using machine learning

Fungal secondary metabolites (SMs) play a significant role in the diversity of ecological communities, niches, and lifestyles in the fungal kingdom. Many fungal SMs have medically and industrially important properties including antifungal, antibacterial, and antitumor activity, and a single metabolite can display multiple types of bioactivities. The genes necessary for fungal SM biosynthesis are typically found in a single genomic region forming biosynthetic gene clusters (BGCs). However, whether fungal SM bioactivity can be predicted from specific attributes of genes in BGCs remains an open question. We adapted previously used machine learning models for predicting SM bioactivity from bacterial BGC data to fungal BGC data. We trained our models to predict antibacterial, antifungal, and cytotoxic/antitumor bioactivity on two datasets: 1) fungal BGCs (dataset comprised of 314 BGCs), and 2) fungal (314 BGCs) and bacterial BGCs (1,003 BGCs); the second dataset was our control since a previous study using just the bacterial BGC data yielded prediction accuracies as high as 80%. We found that the models trained only on fungal BGCs had balanced accuracies between 51-68%, whereas training on bacterial and fungal BGCs yielded balanced accuracies between 61-74%. The lower accuracy of the predictions from fungal data likely stems from the small number of BGCs and SMs with known bioactivity; this lack of data currently limits the application of machine learning approaches in studying fungal secondary metabolism. However, our data also suggest that machine learning approaches trained on bacterial and fungal data can predict SM bioactivity with good accuracy. With more than 15,000 characterized fungal SMs, millions of putative BGCs present in fungal genomes, and increased demand for novel drugs, efforts that systematically link fungal SM bioactivity to BGCs are urgently needed.

A phylogenetic approach to explore the Aspergillus fumigatus conidial surface-associated proteome and its role in pathogenesis

Aspergillus fumigatus , an important pulmonary fungal pathogen causing several diseases collectively called aspergillosis, relies on asexual spores or conidia for initiating host infection. Here, we used a phylogenomic approach to compare proteins in the conidial surface of A. fumigatus , two closely related non-pathogenic species, Aspergillus fischeri and Aspergillus oerlinghausenensis , and the cryptic pathogen Aspergillus lentulus . After identifying 62 proteins uniquely expressed on the A. fumigatus conidial surface, we deleted 42 genes encoding conidial proteins. We found deletion of 33 of these genes altered susceptibility to macrophage killing, penetration and damage to epithelial cells, and cytokine production. Notably, a gene that encodes glycosylasparaginase, which modulates levels of the host pro-inflammatory cytokine IL-1β, is important for infection in an immunocompetent murine model of fungal disease. These results suggest that A. fumigatus conidial surface proteins and effectors are important for evasion and modulation of the immune response at the onset of fungal infection.

Evolutionary origin, population diversity, and diagnostics for a cryptic hybrid pathogen

Cryptic fungal pathogens pose significant identification and disease management challenges due to their morphological resemblance to known pathogenic species while harboring genetic and (often) infectionrelevant trait differences. The cryptic fungal pathogen Aspergillus latus, an allodiploid hybrid originating from Aspergillus spinulosporus and an unknown close relative of Aspergillus quadrilineatus within section Nidulantes, remains poorly understood. The absence of accurate diagnostics for A. latus has led to misidentifications, hindering epidemiological studies and the design of effective treatment plans. We conducted an in-depth investigation of the genomes and phenotypes of 44 globally distributed isolates (41 clinical isolates and three type strains) from Aspergillus section Nidulantes. We found that 21 clinical isolates were A. latus; notably, standard methods of pathogen identification misidentified all A. latus isolates. The remaining isolates were identified as A. spinulosporus (8), A. quadrilineatus (1), or A. nidulans (11). Phylogenomic analyses shed light on the origin of A. latus, indicating one or two hybridization events gave rise to the species during the Miocene, approximately 15.4 to 8.8 million years ago. Characterizing the A. latus pangenome uncovered substantial genetic diversity within gene families and biosynthetic gene clusters. Transcriptomic analysis revealed that both parental genomes are actively expressed in nearly equal proportions and respond to environmental stimuli. Further investigation into infection-relevant chemical and physiological traits, including drug resistance profiles, growth under oxidative stress conditions, and secondary metabolite biosynthesis, highlight distinct phenotypic profiles of the hybrid A. latus compared to its parental and closely related species. Leveraging our comprehensive genomic and phenotypic analyses, we propose five genomic and phenotypic markers as diagnostics for A. latus species identification. These findings provide valuable insights into the evolutionary origin, genomic outcome, and phenotypic implications of hybridization in a cryptic fungal pathogen, thus enhancing our understanding of the underlying processes contributing to fungal pathogenesis. Furthermore, our study underscores the effectiveness of extensive genomic and phenotypic analyses as a promising approach for developing diagnostics applicable to future investigations of cryptic and emerging pathogens.

SWD1 epigenetically chords fungal morphogenesis, aflatoxin biosynthesis, metabolism, and virulence of Aspergillus flavus

As the main producer of aflatoxins, Aspergillus flavus is also one of the most important causes of invasive and non-invasive aspergillosis. Therefore, it is crucial to unravel the regulatory mechanisms of growth, metabolism, and pathogenicity of A. flavus. SWD1 is highly conserved across species for maintaining COMPASS methyltransferase activity, but the bio-function of SWD1 in A. flavus has not been explored. Through genetic analysis, this study revealed that SWD1 is involved in fungal morphogenesis and AFB1 biosynthesis by regulating the orthodox pathways through H3K4me1-3. Stresses sensitivity and crop models analysis revealed that SWD1 is a key regulator for the resistance of A. flavus to adapt to extreme adverse environments and to colonize crop kernels. It also revealed that the WD40 domain and 25 aa highly conserved sequence are indispensable for SWD1 in the regulation of mycotoxin bio-synthesis and fungal virulence. Metabolomic analysis inferred that SWD1 is crucial for the biosynthesis of numerous primary and secondary metabolites, regulates biological functions by reshaping the whole metabolic process, and may inhibit fungal virulence by inducing the apoptosis of mycelia through the inducer sphingosine. This study elucidates the epigenetic mechanism of SWD1 in regulating fungal pathogenicity and mycotoxin biosynthesis, and provides a potential novel target for controlling the virulence of A. flavus.

Use of a human small airway epithelial cell line to study the interactions of Aspergillus fumigatus with pulmonary epithelial cells

During the initiation of invasive aspergillosis, inhaled Aspergillus fumigatus conidia are deposited on the epithelial cells lining the bronchi, terminal bronchioles, and alveoli. While the interactions of A. fumigatus with bronchial and type II alveolar cell lines have been investigated in vitro , little is known about the interactions of this fungus with terminal bronchiolar epithelial cells. We compared the interactions of A. fumigatus with the A549 type II alveolar epithelial cell line and the HSAEC1-KT human small airway epithelial (HSAE) cell line. We found that A. fumigatus conidia were poorly endocytosed by A549 cells, but avidly endocytosed by HSAE cells. A. fumigatus germlings invaded both cell types by induced endocytosis, but not by active penetration. A549 cell endocytosis of A. fumigatus was independent of fungal viability, more dependent on host microfilaments than microtubules, and induced by A. fumigatus CalA interacting with host cell integrin α5β1. By contrast, HSAE cell endocytosis required fungal viability, was more dependent on microtubules than microfilaments, and did not require CalA or integrin α5β1. HSAE cells were more susceptible than A549 cells to damage caused by direct contact with killed A. fumigatus germlings and by secreted fungal products. In response to A. fumigatus infection, A549 cells secreted a broader profile of cytokines and chemokines than HSAE cells. Taken together, these results demonstrate that studies of HSAE cells provide complementary data to A549 cells and thus represent a useful model for probing the interactions of A. fumigatus with bronchiolar epithelial cells in vitro .

Importance

During the initiation of invasive aspergillosis, Aspergillus fumigatus invades, damages, and stimulates the epithelial cells that line the airways and alveoli. Previous studies of A. fumigatus - epithelial cell interactions in vitro have used either large airway epithelial cell lines or the A549 type II alveolar epithelial cell line. The interactions of fungi with terminal bronchiolar epithelial cells have not been investigated. Here, we compared the interactions of A. fumigatus with A549 cells and the Tert-immortalized human small airway epithelial HSAEC1-KT (HSAE) cell line. We discovered that A. fumigatus invades and damages these two cell lines by distinct mechanisms. Also, the proinflammatory responses of the cell lines to A. fumigatus are different. These results provide insight into how A. fumigatus interacts with different types of epithelial cells during invasive aspergillosis and demonstrate that HSAE cells are useful in vitro model for investigating the interactions of this fungus with bronchiolar epithelial cells.

The pan-genome of Aspergillus fumigatus provides a high-resolution view of its population structure revealing high levels of lineage-specific diversity driven by recombination

Aspergillus fumigatus is a deadly agent of human fungal disease where virulence heterogeneity is thought to be at least partially structured by genetic variation between strains. While population genomic analyses based on reference genome alignments offer valuable insights into how gene variants are distributed across populations, these approaches fail to capture intraspecific variation in genes absent from the reference genome. Pan-genomic analyses based on de novo assemblies offer a promising alternative to reference-based genomics with the potential to address the full genetic repertoire of a species. Here, we evaluate 260 genome sequences of A. fumigatus including 62 newly sequenced strains, using a combination of population genomics, phylogenomics, and pan-genomics. Our results offer a high-resolution assessment of population structure and recombination frequency, phylogenetically structured gene presence-absence variation, evidence for metabolic specificity, and the distribution of putative antifungal resistance genes. Although A. fumigatus disperses primarily via asexual conidia, we identified extraordinarily high levels of recombination with the lowest linkage disequilibrium decay value reported for any fungal species to date. We provide evidence for 3 primary populations of A. fumigatus, with recombination occurring only rarely between populations and often within them. These 3 populations are structured by both gene variation and distinct patterns of gene presence-absence with unique suites of accessory genes present exclusively in each clade. Accessory genes displayed functional enrichment for nitrogen and carbohydrate metabolism suggesting that populations may be stratified by environmental niche specialization. Similarly, the distribution of antifungal resistance genes and resistance alleles were often structured by phylogeny. Altogether, the pan-genome of A. fumigatus represents one of the largest fungal pan-genomes reported to date including many genes unrepresented in the Af293 reference genome. These results highlight the inadequacy of relying on a single-reference genome-based approach for evaluating intraspecific variation and the power of combined genomic approaches to elucidate population structure, genetic diversity, and putative ecological drivers of clinically relevant fungi.

Genome-wide association mapping reveals genes underlying population-level metabolome diversity in a fungal crop pathogen

Background

Fungi produce a wide range of specialized metabolites (SMs) involved in biotic interactions. Pathways for the production of SMs are often encoded in clusters of tightly arranged genes identified as biosynthetic gene clusters. Such gene clusters can undergo horizontal gene transfers between species and rapid evolutionary change within species. The acquisition, rearrangement, and deletion of gene clusters can generate significant metabolome diversity. However, the genetic basis underlying variation in SM production remains poorly understood.

Results

Here, we analyzed the metabolite production of a large population of the fungal pathogen of wheat, Zymoseptoria tritici. The pathogen causes major yield losses and shows variation in gene clusters. We performed untargeted ultra-high performance liquid chromatography-high resolution mass spectrometry to profile the metabolite diversity among 102 isolates of the same species. We found substantial variation in the abundance of the detected metabolites among isolates. Integrating whole-genome sequencing data, we performed metabolite genome-wide association mapping to identify loci underlying variation in metabolite production (i.e., metabolite-GWAS). We found that significantly associated SNPs reside mostly in coding and gene regulatory regions. Associated genes encode mainly transport and catalytic activities. The metabolite-GWAS identified also a polymorphism in the 3′UTR region of a virulence gene related to metabolite production and showing expression variation.

Conclusions

Taken together, our study provides a significant resource to unravel polymorphism underlying metabolome diversity within a species. Integrating metabolome screens should be feasible for a range of different plant pathogens and help prioritize molecular studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01422-z.

Aspergillus-bees: A dynamic symbiotic association

Besides representing one of the most relevant threats of fungal origin to human and animal health, the genus Aspergillus includes opportunistic pathogens which may infect bees (Hymenoptera, Apoidea) in all developmental stages. At least 30 different species of Aspergillus have been isolated from managed and wild bees. Some efficient behavioral responses (e.g., diseased brood removal) exerted by bees negatively affect the chance to diagnose the pathology, and may contribute to the underestimation of aspergillosis importance in beekeeping. On the other hand, bee immune responses may be affected by biotic and abiotic stresses and suffer from the loose co-evolutionary relationships with Aspergillus pathogenic strains. However, if not pathogenic, these hive mycobiota components can prove to be beneficial to bees, by affecting the interaction with other pathogens and parasites and by detoxifying xenobiotics. The pathogenic aptitude of Aspergillus spp. likely derives from the combined action of toxins and hydrolytic enzymes, whose effects on bees have been largely overlooked until recently. Variation in the production of these virulence factors has been observed among strains, even belonging to the same species. Toxigenic and non-toxigenic strains/species may co-exist in a homeostatic equilibrium which is susceptible to be perturbed by several external factors, leading to mutualistic/antagonistic switch in the relationships between Aspergillus and bees.

The histone code of the fungal genus Aspergillus uncovered by evolutionary and proteomic analyses

Chemical modifications of DNA and histone proteins impact the organization of chromatin within the nucleus. Changes in these modifications, catalysed by different chromatin-modifying enzymes, influence chromatin organization, which in turn is thought to impact the spatial and temporal regulation of gene expression. While combinations of different histone modifications, the histone code, have been studied in several model species, we know very little about histone modifications in the fungal genus Aspergillus, whose members are generally well studied due to their importance as models in cell and molecular biology as well as their medical and biotechnological relevance. Here, we used phylogenetic analyses in 94 Aspergilli as well as other fungi to uncover the occurrence and evolutionary trajectories of enzymes and protein complexes with roles in chromatin modifications or regulation. We found that these enzymes and complexes are highly conserved in Aspergilli, pointing towards a complex repertoire of chromatin modifications. Nevertheless, we also observed few recent gene duplications or losses, highlighting Aspergillus species to further study the roles of specific chromatin modifications. SET7 (KMT6) and other components of PRC2 (Polycomb Repressive Complex 2), which is responsible for methylation on histone H3 at lysine 27 in many eukaryotes including fungi, are absent in Aspergilli as well as in closely related Penicillium species, suggesting that these lost the capacity for this histone modification. We corroborated our computational predictions by performing untargeted MS analysis of histone post-translational modifications in Aspergillus nidulans. This systematic analysis will pave the way for future research into the complexity of the histone code and its functional implications on genome architecture and gene regulation in fungi.

Examination of Genome-Wide Ortholog Variation in Clinical and Environmental Isolates of the Fungal Pathogen Aspergillus fumigatus

Aspergillus fumigatus is both an environmental saprobe and an opportunistic human fungal pathogen. Knowledge of genomic variation across A. fumigatus isolates is essential for understanding the evolution of pathogenicity, virulence, and resistance to antifungal drugs. Here, we investigated 206 A. fumigatus isolates (133 clinical and 73 environmental isolates), aiming to identify genes with variable presence across isolates and test whether this variation was related to the clinical or environmental origin of isolates. The PanOrtho genome of A. fumigatus consists of 13,085 ortholog groups, of which 7,773 (59.4%) are shared by all isolates (core groups) and 5,312 (40.6%) vary in their gene presence across isolates (accessory groups plus singletons). Despite differences in the distribution of orthologs across all isolates, no significant differences were observed among clinical versus environmental isolates when phylogeny was accounted for. Orthologs that differ in their distribution across isolates tend to occur at low frequency and/or be restricted to specific isolates; thus, the degree of genomic conservation between orthologs of A. fumigatus is high. These results suggest that differences in the distribution of orthologs within A. fumigatus cannot be associated with the clinical or environmental origin of isolates. IMPORTANCE Aspergillus fumigatus is a cosmopolitan species of fungus responsible for thousands of cases of invasive disease annually. Clinical and environmental isolates of A. fumigatus exhibit extensive phenotypic differences, including differences related to virulence and antifungal drug resistance. A comprehensive survey of the genomic diversity present in A. fumigatus and its relationship to the clinical or environmental origin of isolates can contribute to the prediction of the mechanisms of evolution and infection of the species. Our results suggest that there is no significant variation in ortholog distribution between clinical and environmental isolates when accounting for evolutionary history. The work supports the hypothesis that environmental and clinical isolates of A. fumigatus do not differ in their gene contents.

Fungal-fungal co-culture: a primer for generating chemical diversity

Covering: 2002 to 2020In their natural environment, fungi must compete for resources. It has been hypothesized that this competition likely induces the biosynthesis of secondary metabolites for defence. In a quest to discover new chemical diversity from fungal cultures, a growing trend has been to recapitulate this competitive environment in the laboratory, essentially growing fungi in co-culture. This review covers fungal-fungal co-culture studies beginning with the first literature report in 2002. Since then, there has been a growing number of new secondary metabolites reported as a result of fungal co-culture studies. Specifically, this review discusses and provides insights into (1) rationale for pairing fungal strains, (2) ways to grow fungi for co-culture, (3) different approaches to screening fungal co-cultures for chemical diversity, (4) determining the secondary metabolite-producing strain, and (5) final thoughts regarding the fungal-fungal co-culture approach. Our goal is to provide a set of practical strategies for fungal co-culture studies to generate unique chemical diversity that the natural products research community can utilize.

A genomic journey in the secondary metabolite diversity of fungal plant and insect pathogens: from functional to population genomics

Fungal pathogens produce a broad array of secondary metabolites (SMs), which allow the fungus to thrive in its natural habitat and gain competitive advantage. Analysis of the genetically encoded blueprints for SM assembly highlighted that only a small portion of the SMs these fungi are capable of producing are known, and even fewer have been investigated for their natural function. Using molecular tools, a lot of progress has been made recently in identifying the blueprint products and linking them to their ecological purpose such as the peptide virulence factor fusaoctaxin A released by Fusarium graminearum during infection of wheat or the F. oxysporum polyketide bikaverin that provides competitive advantage against bacteria in tomato. In addition, population genomics have given particularly important insights into the species-specific plasticity of the SM blueprint arsenal, showcasing the ongoing evolution and adaptation of fungal pathogens. This approach holds promise in inferring roles in pathogenicity of many more fungal SMs.

Extensive non-coding sequence divergence between the major human pathogen Aspergillus fumigatus and its relatives

Invasive aspergillosis is a deadly fungal disease; more than 400,000 patients are infected worldwide each year and the mortality rate can be as high as 50-95%. Of the ~450 species in the genus Aspergillus only a few are known to be clinically relevant, with the major pathogen Aspergillus fumigatus being responsible for ~50% of all invasive mold infections. Genomic comparisons between A. fumigatus and other Aspergillus species have historically focused on protein-coding regions. However, most A. fumigatus genes, including those that modulate its virulence, are also present in other pathogenic and non-pathogenic closely related species. Our hypothesis is that differential gene regulation - mediated through the non-coding regions upstream of genes' first codon - contributes to A. fumigatus pathogenicity. To begin testing this, we compared non-coding regions upstream of the first codon of single-copy orthologous genes from the two A. fumigatus reference strains Af293 and A1163 and eight closely related Aspergillus section Fumigati species. We found that these non-coding regions showed extensive sequence variation and lack of homology across species. By examining the evolutionary rates of both protein-coding and non-coding regions in a subset of orthologous genes with highly conserved non-coding regions across the phylogeny, we identified 418 genes, including 25 genes known to modulate A. fumigatus virulence, whose non-coding regions exhibit a different rate of evolution in A. fumigatus. Examination of sequence alignments of these non-coding regions revealed numerous instances of insertions, deletions, and other types of mutations of at least a few nucleotides in A. fumigatus compared to its close relatives. These results show that closely related Aspergillus species that vary greatly in their pathogenicity exhibit extensive non-coding sequence variation and identify numerous changes in non-coding regions of A. fumigatus genes known to contribute to virulence.

Evolution of the human pathogenic lifestyle in fungi

Fungal pathogens cause more than a billion human infections every year, resulting in more than 1.6 million deaths annually. Understanding the natural history and evolutionary ecology of fungi is helping us understand how disease-relevant traits have repeatedly evolved. Different types and mechanisms of genetic variation have contributed to the evolution of fungal pathogenicity and specific genetic differences distinguish pathogens from non-pathogens. Insights into the traits, genetic elements, and genetic and ecological mechanisms that contribute to the evolution of fungal pathogenicity are crucial for developing strategies to both predict emergence of fungal pathogens and develop drugs to combat them.

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.

Heterogeneity in the transcriptional response of the human pathogen Aspergillus fumigatus to the antifungal agent caspofungin

Aspergillus fumigatus is the main causative agent of invasive pulmonary aspergillosis (IPA), a severe disease that affects immunosuppressed patients worldwide. The fungistatic drug caspofungin (CSP) is the second line of therapy against IPA but has increasingly been used against clinical strains that are resistant to azoles, the first line antifungal therapy. In high concentrations, CSP induces a tolerance phenotype with partial reestablishment of fungal growth called CSP paradoxical effect (CPE), resulting from a change in the composition of the cell wall. An increasing number of studies has shown that different isolates of A. fumigatus exhibit phenotypic heterogeneity, including heterogeneity in their CPE response. To gain insights into the underlying molecular mechanisms of CPE response heterogeneity, we analyzed the transcriptomes of two A. fumigatus reference strains, Af293 and CEA17, exposed to low and high CSP concentrations. We found that there is a core transcriptional response that involves genes related to cell wall remodeling processes, mitochondrial function, transmembrane transport, and amino acid and ergosterol metabolism, and a variable response related to secondary metabolite (SM) biosynthesis and iron homeostasis. Specifically, we show here that the overexpression of a SM pathway that works as an iron chelator extinguishes the CPE in both backgrounds, whereas iron depletion is detrimental for the CPE in Af293 but not in CEA17. We next investigated the function of the transcription factor CrzA, whose deletion was previously shown to result in heterogeneity in the CPE response of the Af293 and CEA17 strains. We found that CrzA constitutively binds to and modulates the expression of several genes related to processes involved in CSP tolerance and that crzA deletion differentially impacts the SM production and growth of Af293 and CEA17. As opposed to the ΔcrzACEA17 mutant, the ΔcrzAAf293 mutant fails to activate cell wall remodeling genes upon CSP exposure, which most likely severely affects its macrostructure and extinguishes its CPE. This study describes how heterogeneity in the response to an antifungal agent between A. fumigatus strains stems from heterogeneity in the function of a transcription factor and its downstream target genes.

Regulation of gliotoxin biosynthesis and protection in Aspergillus species

Aspergillus fumigatus causes a range of human and animal diseases collectively known as aspergillosis. A. fumigatus possesses and expresses a range of genetic determinants of virulence, which facilitate colonisation and disease progression, including the secretion of mycotoxins. Gliotoxin (GT) is the best studied A. fumigatus mycotoxin with a wide range of known toxic effects that impair human immune cell function. GT is also highly toxic to A. fumigatus and this fungus has evolved self-protection mechanisms that include (i) the GT efflux pump GliA, (ii) the GT neutralising enzyme GliT, and (iii) the negative regulation of GT biosynthesis by the bis-thiomethyltransferase GtmA. The transcription factor (TF) RglT is the main regulator of GliT and this GT protection mechanism also occurs in the non-GT producing fungus A. nidulans. However, the A. nidulans genome does not encode GtmA and GliA. This work aimed at analysing the transcriptional response to exogenous GT in A. fumigatus and A. nidulans, two distantly related Aspergillus species, and to identify additional components required for GT protection. RNA-sequencing shows a highly different transcriptional response to exogenous GT with the RglT-dependent regulon also significantly differing between A. fumigatus and A. nidulans. However, we were able to observe homologs whose expression pattern was similar in both species (43 RglT-independent and 11 RglT-dependent). Based on this approach, we identified a novel RglT-dependent methyltranferase, MtrA, involved in GT protection. Taking into consideration the occurrence of RglT-independent modulated genes, we screened an A. fumigatus deletion library of 484 transcription factors (TFs) for sensitivity to GT and identified 15 TFs important for GT self-protection. Of these, the TF KojR, which is essential for kojic acid biosynthesis in Aspergillus oryzae, was also essential for virulence and GT biosynthesis in A. fumigatus, and for GT protection in A. fumigatus, A. nidulans, and A. oryzae. KojR regulates rglT, gliT, gliJ expression and sulfur metabolism in Aspergillus species. Together, this study identified conserved components required for GT protection in Aspergillus species.

A. fumigatus secretes mycotoxins that are essential for its virulence and pathogenicity. Gliotoxin (GT) is a sulfur-containing mycotoxin, which is known to impair several aspects of the human immune response. GT is also toxic to different fungal species, which have evolved several GT protection strategies. To further decipher these responses, we used transcriptional profiling aiming to compare the response to GT in the GT producer A. fumigatus and the GT non-producer A. nidulans. This analysis allowed us to identify additional genes with a potential role in GT protection. We also identified 15 transcription factors (TFs) encoded in the A. fumigatus genome that are important for conferring resistance to exogenous gliotoxin. One of these TFs, KojR, which is essential for A. oryzae kojic acid production, is also important for virulence in A. fumigatus and GT protection in A. fumigatus, A. nidulans and A. oryzae. KojR regulates the expression of genes important for gliotoxin biosynthesis and protection, and sulfur metabolism. Together, this work identified conserved components required for gliotoxin protection in Aspergillus species.

Phylogenetic Distribution of csp1 Types in Aspergillus fumigatus and Their Correlates to Azole Antifungal Drug Resistance

In Aspergillus fumigatus, the repetitive region of the csp1 gene is one of the most frequently used loci for intraspecies typing of this human pathogenic mold. Using PCR amplification and Sanger sequencing of only a single marker, csp1 typing is readily available to most laboratories and highly reproducible. Here, I evaluate the usefulness of the csp1 marker for resistance detection and epidemiologic stratification among A. fumigatus isolates. After resolving nomenclature conflicts from published studies and adding novel csp1 types, the number of known types now adds up to 38. Their distribution mostly correlates with A. fumigatus population structure, and they are also meaningful for narrowly defined cases of azole resistance phenotypes. Isolates carrying the pandemic resistance allele TR34/L98H show signs of interclade crossing of strains with t02 or t04A, into the t11 clade. Furthermore, absolute differences in voriconazole MIC values between t02/t04B versus t11 TR34/L98H isolates indicate that the genetic background of resistance mutations may have a pivotal role in cross-resistance phenotypes and, thus, clinical outcome and environmental selection. Despite the general genetic similarity of isolates with identical csp1 types, outcrossing into other clades is also observed. The csp1 type alone, therefore, does not sufficiently discriminate genetic clades to be used as the sole marker in epidemiologic studies. IMPORTANCE Aspergillus fumigatus is a ubiquitously distributed saprophytic mold and a leading cause of invasive aspergillosis in human hosts. Pandemic azole-resistant strains have emerged on a global scale, which are thought to be propagated through use of azole-based fungicides in agriculture. To perform epidemiologic studies, genetic typing of large cohorts is key. Here, I evaluate the usefulness of the frequently used csp1 marker for resistance detection and epidemiologic stratification among A. fumigatus isolates. The phylogenetic distribution of csp1 types mostly correlates with A. fumigatus population structure and is also meaningful for narrowly defined cases of azole resistance phenotypes. Nevertheless, outcrossing of csp1 into other clades is also observed. The csp1 type alone, therefore, does not sufficiently discriminate genetic clades and should not be used as the sole marker in epidemiologic studies.

Aspergillus fumigatus pan-genome analysis identifies genetic variants associated with human infection

Aspergillus fumigatus is an environmental saprobe and opportunistic human fungal pathogen. Despite an estimated annual occurrence of more than 300,000 cases of invasive disease worldwide, a comprehensive survey of the genomic diversity present in A. fumigatus-including the relationship between clinical and environmental isolates and how this genetic diversity contributes to virulence and antifungal drug resistance-has been lacking. In this study we define the pan-genome of A. fumigatus using a collection of 300 globally sampled genomes (83 clinical and 217 environmental isolates). We found that 7,563 of the 10,907 unique orthogroups (69%) are core and present in all isolates and the remaining 3,344 show presence/absence of variation, representing 16-22% of the genome of each isolate. Using this large genomic dataset of environmental and clinical samples, we found an enrichment for clinical isolates in a genetic cluster whose genomes also contain more accessory genes, including genes coding for transmembrane transporters and proteins with iron-binding activity, and genes involved in both carbohydrate and amino-acid metabolism. Finally, we leverage the power of genome-wide association studies to identify genomic variation associated with clinical isolates and triazole resistance as well as characterize genetic variation in known virulence factors. This characterization of the genomic diversity of A. fumigatus allows us to move away from a single reference genome that does not necessarily represent the species as a whole and better understand its pathogenic versatility, ultimately leading to better management of these infections.

Metabolomics and genomics in natural products research: complementary tools for targeting new chemical entities

Covering: 2010 to 2021Organisms in nature have evolved into proficient synthetic chemists, utilizing specialized enzymatic machinery to biosynthesize an inspiring diversity of secondary metabolites. Often serving to boost competitive advantage for their producers, these secondary metabolites have widespread human impacts as antibiotics, anti-inflammatories, and antifungal drugs. The natural products discovery field has begun a shift away from traditional activity-guided approaches and is beginning to take advantage of increasingly available metabolomics and genomics datasets to explore undiscovered chemical space. Major strides have been made and now enable -omics-informed prioritization of chemical structures for discovery, including the prospect of confidently linking metabolites to their biosynthetic pathways. Over the last decade, more integrated strategies now provide researchers with pipelines for simultaneous identification of expressed secondary metabolites and their biosynthetic machinery. However, continuous collaboration by the natural products community will be required to optimize strategies for effective evaluation of natural product biosynthetic gene clusters to accelerate discovery efforts. Here, we provide an evaluative guide to scientific literature as it relates to studying natural product biosynthesis using genomics, metabolomics, and their integrated datasets. Particular emphasis is placed on the unique insights that can be gained from large-scale integrated strategies, and we provide source organism-specific considerations to evaluate the gaps in our current knowledge.

Fungal quinones: diversity, producers, and applications of quinones from Aspergillus, Penicillium, Talaromyces, Fusarium, and Arthrinium

Quinones represent an important group of highly structurally diverse, mainly polyketide-derived secondary metabolites widely distributed among filamentous fungi. Many quinones have been reported to have important biological functions such as inhibition of bacteria or repression of the immune response in insects. Other quinones, such as ubiquinones are known to be essential molecules in cellular respiration, and many quinones are known to protect their producing organisms from exposure to sunlight. Most recently, quinones have also attracted a lot of industrial interest since their electron-donating and -accepting properties make them good candidates as electrolytes in redox flow batteries, like their often highly conjugated double bond systems make them attractive as pigments. On an industrial level, quinones are mainly synthesized from raw components in coal tar. However, the possibility of producing quinones by fungal cultivation has great prospects since fungi can often be grown in industrially scaled bioreactors, producing valuable metabolites on cheap substrates. In order to give a better overview of the secondary metabolite quinones produced by and shared between various fungi, mainly belonging to the genera Aspergillus, Penicillium, Talaromyces, Fusarium, and Arthrinium, this review categorizes quinones into families such as emodins, fumigatins, sorbicillinoids, yanuthones, and xanthomegnins, depending on structural similarities and information about the biosynthetic pathway from which they are derived, whenever applicable. The production of these quinone families is compared between the different genera, based on recently revised taxonomy. KEY POINTS: • Quinones represent an important group of secondary metabolites widely distributed in important fungal genera such as Aspergillus, Penicillium, Talaromyces, Fusarium, and Arthrinium. • Quinones are of industrial interest and can be used in pharmacology, as colorants and pigments, and as electrolytes in redox flow batteries. • Quinones are grouped into families and compared between genera according to the revised taxonomy.

Examining Signatures of Natural Selection in Antifungal Resistance Genes Across Aspergillus Fungi

Certain Aspergillus fungi cause aspergillosis, a set of diseases that typically affect immunocompromised individuals. Most cases of aspergillosis are caused by Aspergillus fumigatus, which infects millions of people annually. Some closely related so-called cryptic species, such as Aspergillus lentulus, can also cause aspergillosis, albeit at lower frequencies, and they are also clinically relevant. Few antifungal drugs are currently available for treating aspergillosis and there is increasing worldwide concern about the presence of antifungal drug resistance in Aspergillus species. Furthermore, isolates from both A. fumigatus and other Aspergillus pathogens exhibit substantial heterogeneity in their antifungal drug resistance profiles. To gain insights into the evolution of antifungal drug resistance genes in Aspergillus, we investigated signatures of positive selection in 41 genes known to be involved in drug resistance across 42 susceptible and resistant isolates from 12 Aspergillus section Fumigati species. Using codon-based site models of sequence evolution, we identified ten genes that contain 43 sites with signatures of ancient positive selection across our set of species. None of the sites that have experienced positive selection overlap with sites previously reported to be involved in drug resistance. These results identify sites that likely experienced ancient positive selection in Aspergillus genes involved in resistance to antifungal drugs and suggest that historical selective pressures on these genes likely differ from any current selective pressures imposed by antifungal drugs.

Genomic and Phenotypic Analysis of COVID-19-Associated Pulmonary Aspergillosis Isolates of Aspergillus fumigatus

The ongoing global pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for coronavirus disease 2019 (COVID-19), first described in Wuhan, China. A subset of COVID-19 patients has been reported to have acquired secondary infections by microbial pathogens, such as opportunistic fungal pathogens from the genus Aspergillus. To gain insight into COVID-19-associated pulmonary aspergillosis (CAPA), we analyzed the genomes and characterized the phenotypic profiles of four CAPA isolates of Aspergillus fumigatus obtained from patients treated in the area of North Rhine-Westphalia, Germany. By examining the mutational spectrum of single nucleotide polymorphisms, insertion-deletion polymorphisms, and copy number variants among 206 genes known to modulate A. fumigatus virulence, we found that CAPA isolate genomes do not exhibit significant differences from the genome of the Af293 reference strain. By examining a number of factors, including virulence in an invertebrate moth model, growth in the presence of osmotic, cell wall, and oxidative stressors, secondary metabolite biosynthesis, and the MIC of antifungal drugs, we found that CAPA isolates were generally, but not always, similar to A. fumigatus reference strains Af293 and CEA17. Notably, CAPA isolate D had more putative loss-of-function mutations in genes known to increase virulence when deleted. Moreover, CAPA isolate D was significantly more virulent than the other three CAPA isolates and the A. fumigatus reference strains Af293 and CEA17, but similarly virulent to two other clinical strains of A. fumigatus. These findings expand our understanding of the genomic and phenotypic characteristics of isolates that cause CAPA. IMPORTANCE The global pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the etiological agent of coronavirus disease 2019 (COVID-19), has already killed millions of people. COVID-19 patient outcome can be further complicated by secondary infections, such as COVID-19-associated pulmonary aspergillosis (CAPA). CAPA is caused by Aspergillus fungal pathogens, but there is little information about the genomic and phenotypic characteristics of CAPA isolates. We conducted genome sequencing and extensive phenotyping of four CAPA isolates of Aspergillus fumigatus from Germany. We found that CAPA isolates were often, but not always, similar to other reference strains of A. fumigatus across 206 genetic determinants of infection-relevant phenotypes, including virulence. For example, CAPA isolate D was more virulent than other CAPA isolates and reference strains in an invertebrate model of fungal disease, but similarly virulent to two other clinical strains. These results expand our understanding of COVID-19-associated pulmonary aspergillosis.

Aspergillus fumigatus and aspergillosis: From basics to clinics

The airborne fungus Aspergillus fumigatus poses a serious health threat to humans by causing numerous invasive infections and a notable mortality in humans, especially in immunocompromised patients. Mould-active azoles are the frontline therapeutics employed to treat aspergillosis. The global emergence of azole-resistant A. fumigatus isolates in clinic and environment, however, notoriously limits the therapeutic options of mould-active antifungals and potentially can be attributed to a mortality rate reaching up to 100 %. Although specific mutations in CYP 51A are the main cause of azole resistance, there is a new wave of azole-resistant isolates with wild-type CYP 51A genotype challenging the efficacy of the current diagnostic tools. Therefore, applications of whole-genome sequencing are increasingly gaining popularity to overcome such challenges. Prominent echinocandin tolerance, as well as liver and kidney toxicity posed by amphotericin B, necessitate a continuous quest for novel antifungal drugs to combat emerging azole-resistant A. fumigatus isolates. Animal models and the tools used for genetic engineering require further refinement to facilitate a better understanding about the resistance mechanisms, virulence, and immune reactions orchestrated against A. fumigatus. This review paper comprehensively discusses the current clinical challenges caused by A. fumigatus and provides insights on how to address them.

Comparison of Two Aspergillus oryzae Genomes From Different Clades Reveals Independent Evolution of Alpha-Amylase Duplication, Variation in Secondary Metabolism Genes, and Differences in Primary Metabolism

Microbes (bacteria, yeasts, molds), in addition to plants and animals, were domesticated for their roles in food preservation, nutrition and flavor. Aspergillus oryzae is a domesticated filamentous fungal species traditionally used during fermentation of Asian foods and beverage, such as sake, soy sauce, and miso. To date, little is known about the extent of genome and phenotypic variation of A. oryzae isolates from different clades. Here, we used long-read Oxford Nanopore and short-read Illumina sequencing to produce a highly accurate and contiguous genome assemble of A. oryzae 14160, an industrial strain from China. To understand the relationship of this isolate, we performed phylogenetic analysis with 90 A. oryzae isolates and 1 isolate of the A. oryzae progenitor, Aspergillus flavus. This analysis showed that A. oryzae 14160 is a member of clade A, in comparison to the RIB 40 type strain, which is a member of clade F. To explore genome variation between isolates from distinct A. oryzae clades, we compared the A. oryzae 14160 genome with the complete RIB 40 genome. Our results provide evidence of independent evolution of the alpha-amylase gene duplication, which is one of the major adaptive mutations resulting from domestication. Synteny analysis revealed that both genomes have three copies of the alpha-amylase gene, but only one copy on chromosome 2 was conserved. While the RIB 40 genome had additional copies of the alpha-amylase gene on chromosomes III, and V, 14160 had a second copy on chromosome II and an third copy on chromosome VI. Additionally, we identified hundreds of lineage specific genes, and putative high impact mutations in genes involved in secondary metabolism, including several of the core biosynthetic genes. Finally, to examine the functional effects of genome variation between strains, we measured amylase activity, proteolytic activity, and growth rate on several different substrates. RIB 40 produced significantly higher levels of amylase compared to 14160 when grown on rice and starch. Accordingly, RIB 40 grew faster on rice, while 14160 grew faster on soy. Taken together, our analyses reveal substantial genome and phenotypic variation within A. oryzae.

An evolutionary genomic approach reveals both conserved and species-specific genetic elements related to human disease in closely related Aspergillus fungi

Aspergillosis is an important opportunistic human disease caused by filamentous fungi in the genus Aspergillus. Roughly 70% of infections are caused by Aspergillus fumigatus, with the rest stemming from approximately a dozen other Aspergillus species. Several of these pathogens are closely related to A. fumigatus and belong in the same taxonomic section, section Fumigati. Pathogenic species are frequently most closely related to nonpathogenic ones, suggesting Aspergillus pathogenicity evolved multiple times independently. To understand the repeated evolution of Aspergillus pathogenicity, we performed comparative genomic analyses on 18 strains from 13 species, including 8 species in section Fumigati, which aimed to identify genes, both ones previously connected to virulence as well as ones never before implicated, whose evolution differs between pathogens and nonpathogens. We found that most genes were present in all species, including approximately half of those previously connected to virulence, but a few genes were section- or species-specific. Evolutionary rate analyses identified over 1700 genes whose evolutionary rate differed between pathogens and nonpathogens and dozens of genes whose rates differed between specific pathogens and the rest of the taxa. Functional testing of deletion mutants of 17 transcription factor-encoding genes whose evolution differed between pathogens and nonpathogens identified eight genes that affect either fungal survival in a model of phagocytic killing, host survival in an animal model of fungal disease, or both. These results suggest that the evolution of pathogenicity in Aspergillus involved both conserved and species-specific genetic elements, illustrating how an evolutionary genomic approach informs the study of fungal disease.

Interspecies Genomic Variation and Transcriptional Activeness of Secondary Metabolism-Related Genes in Aspergillus Section Fumigati

Filamentous fungi produce various bioactive compounds that are biosynthesized by sets of proteins encoded in biosynthesis gene clusters (BGCs). For an unknown reason, many BGCs are transcriptionally silent in laboratory conditions, which has hampered the discovery of novel fungal compounds. The transcriptional reactiveness of fungal secondary metabolism is not fully understood. To gain the comprehensive view, we conducted comparative genomic and transcriptomic analyses of nine closely-related species of Aspergillus section Fumigati (A. fumigatus, A. fumigatiaffinis, A. novofumigatus, A. thermomutatus, A. viridinutans, A. pseudoviridinutans, A. lentulus, A. udagawae, and Neosartorya fischeri). For expanding our knowledge, we newly sequenced genomes of A. viridinutans and A. pseudoviridinutans, and reassembled and reannotated the previously released genomes of A. lentulus and A. udagawae. Between 34 and 84 secondary metabolite (SM) backbone genes were identified in the genomes of these nine respective species, with 8.7-51.2% being unique to the species. A total of 247 SM backbone gene types were identified in the nine fungi. Ten BGCs are shared by all nine species. Transcriptomic analysis using A. fumigatus, A. lentulus, A. udagawae, A. viridinutans, and N. fischeri was conducted to compare expression levels of all SM backbone genes in four different culture conditions; 32-83% of SM backbone genes in these species were not expressed in the tested conditions, which reconfirmed that large part of fungal SM genes are hard to be expressed. The species-unique SM genes of the five species were expressed with lower frequency (18.8% in total) than the SM genes that are conserved in all five species (56%). These results suggest that the expression tendency of BGCs is correlated with their interspecies distribution pattern. Our findings increase understanding of the evolutionary processes associated with the regulation of fungal secondary metabolism.