AtfA bZIP-type transcription factor regulates oxidative and osmotic stress responses in Aspergillus nidulans

HacA, a key transcription factor for the unfolded protein response, is required for fungal development, aflatoxin biosynthesis and pathogenicity of Aspergillus flavus

Aspergillus flavus is a fungus notorious for contaminating food and feed with aflatoxins. As a saprophytic fungus, it secretes large amounts of enzymes to access nutrients, making endoplasmic reticulum (ER) homeostasis important for protein folding and secretion. The role of HacA, a key transcription factor in the unfolded protein response pathway, remains poorly understood in A. flavus. In this study, the hacA gene in A. flavus was knockout. Results showed that the absence of hacA led to a decreased pathogenicity of the strain, as it failed to colonize intact maize kernels. This may be due to retarded vegetable growth, especially the abnormal development of swollen tips and shorter hyphal septa. Deletion of hacA also hindered conidiogenesis and sclerotial development. Notably, the mutant strain failed to produce aflatoxin B1. Moreover, compared to the wild type, the mutant strain showed increased sensitivity to ER stress inducer such as Dithiothreitol (DTT), and heat stress. It also displayed heightened sensitivity to other environmental stresses, including cell wall, osmotic, and pH stresses. Further transcriptomic analysis revealed the involvement of the hacA in numerous biological processes, including filamentous growth, asexual reproduction, mycotoxin biosynthetic process, signal transduction, budding cell apical bud growth, invasive filamentous growth, response to stimulus, and so on. Taken together, HacA plays a vital role in fungal development, pathogenicity and aflatoxins biosynthesis. This highlights the potential of targeting hacA as a novel approach for early prevention of A. flavus contamination.

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

BACKGROUND

Aspergillus flavus is an important agricultural and food safety threat due to its production of carcinogenic aflatoxins. It has high level of genetic diversity that is adapted to various environments. Recently, we reported two reference genomes of A. flavus isolates, AF13 (MAT1-2 and highly aflatoxigenic isolate) and NRRL3357 (MAT1-1 and moderate aflatoxin producer). Where, an insertion of 310 kb in AF13 included an aflatoxin producing gene bZIP transcription factor, named atfC. Observations of significant genomic variants between these isolates of contrasting phenotypes prompted an investigation into variation among other agricultural isolates of A. flavus with the goal of discovering novel genes potentially associated with aflatoxin production regulation. Present study was designed with three main objectives: (1) collection of large number of A. flavus isolates from diverse sources including maize plants and field soils; (2) whole genome sequencing of collected isolates and development of a pangenome; and (3) pangenome-wide association study (Pan-GWAS) to identify novel secondary metabolite cluster genes.

RESULTS

Pangenome analysis of 346 A. flavus isolates identified a total of 17,855 unique orthologous gene clusters, with mere 41% (7,315) core genes and 59% (10,540) accessory genes indicating accumulation of high genomic diversity during domestication. 5,994 orthologous gene clusters in accessory genome not annotated in either the A. flavus AF13 or NRRL3357 reference genomes. Pan-genome wide association analysis of the genomic variations identified 391 significant associated pan-genes associated with aflatoxin production. Interestingly, most of the significantly associated pan-genes (94%; 369 associations) belonged to accessory genome indicating that genome expansion has resulted in the incorporation of new genes associated with aflatoxin and other secondary metabolites.

CONCLUSION

In summary, this study provides complete pangenome framework for the species of Aspergillus flavus along with associated genes for pathogen survival and aflatoxin production. The large accessory genome indicated large genome diversity in the species A. flavus, however AflaPan is a closed pangenome represents optimum diversity of species A. flavus. Most importantly, the newly identified aflatoxin producing gene clusters will be a new source for seeking aflatoxin mitigation strategies and needs new attention in research.

The bZIP transcription factor ATF1 regulates blue light and oxidative stress responses in Trichoderma guizhouense

In several filamentous fungi, incident light and environmental stress signaling share the mitogen-activated protein kinase (MAPK) HOG (SAK) pathway. It has been revealed that short-term illumination with blue light triggers the activation of the HOG pathway in Trichoderma spp. In this study, we demonstrate the crucial role of the basic leucine zipper transcription factor ATF1 in blue light responses and signaling downstream of the MAPK HOG1 in Trichoderma guizhouense. The lack of ATF1 severely impaired photoconidiation and delayed vegetative growth and conidial germination. Upon blue light or H2O2 stimuli, HOG1 interacted with ATF1 in the nucleus. Genome-wide transcriptome analyses revealed that 61.8% (509 out of 824) and 85.2% (702 out of 824) of blue light-regulated genes depended on ATF1 and HOG1, respectively, of which 58.4% (481 out of 824) were regulated by both of them. Our results also show that blue light promoted conidial germination and HOG1 and ATF1 played opposite roles in controlling conidial germination in the dark. Additionally, the lack of ATF1 led to reduced oxidative stress resistance, probably because of the downregulation of catalase-encoding genes. Overall, our results demonstrate that ATF1 is the downstream component of HOG1 and is responsible for blue light responses, conidial germination, vegetative growth, and oxidative stress resistance in T. guizhouense.

Global Transcriptomic Changes Elicited by sodB Deletion and Menadione Exposure in Aspergillus nidulans

Manganese superoxide dismutases (MnSODs) play a pivotal role in the preservation of mitochondrial integrity and function in fungi under various endogenous and exogenous stresses. Deletion of Aspergillus nidulans mnSOD/SodB increased oxidative stress sensitivity and apoptotic cell death rates as well as affected antioxidant enzyme and sterigmatocystin productions, respiration, conidiation and the stress tolerance of conidiospores. The physiological consequences of the lack of sodB were more pronounced during carbon starvation than in the presence of glucose. Lack of SodB also affected the changes in the transcriptome, recorded by high-throughput RNA sequencing, in menadione sodium bisulfite (MSB)-exposed, submerged cultures supplemented with glucose. Surprisingly, the difference between the global transcriptional changes of the ΔsodB mutant and the control strain were relatively small, indicating that the SodB-dependent maintenance of mitochondrial integrity was not essential under these experimental conditions. Owing to the outstanding physiological flexibility of the Aspergilli, certain antioxidant enzymes and endogenous antioxidants together with the reduction in mitochondrial functions compensated well for the lack of SodB. The lack of sodB reduced the growth of surface cultures more than of the submerged culture, which should be considered in future development of fungal disinfection methods.

Host-derived reactive oxygen species trigger activation of the Candida albicans transcription regulator Rtg1/3

The signals that denote mammalian host environments and dictate the activation of signaling pathways in human-associated microorganisms are often unknown. The transcription regulator Rtg1/3 in the human fungal pathogen Candida albicans is a crucial determinant of host colonization and pathogenicity. Rtg1/3's activity is controlled, in part, by shuttling the regulator between the cytoplasm and nucleus of the fungus. The host signal(s) that Rtg1/3 respond(s) to, however, have remained unclear. Here we report that neutrophil-derived reactive oxygen species (ROS) direct the subcellular localization of this C. albicans transcription regulator. Upon engulfment of Candida cells by human or mouse neutrophils, the regulator shuttles to the fungal nucleus. Using genetic and chemical approaches to disrupt the neutrophils' oxidative burst, we establish that the oxidants produced by the NOX2 complex-but not the oxidants generated by myeloperoxidase-trigger Rtg1/3's migration to the nucleus. Furthermore, screening a collection of C. albicans kinase deletion mutants, we implicate the MKC1 signaling pathway in the ROS-dependent regulation of Rtg1/3 in this fungus. Finally, we show that Rtg1/3 contributes to C. albicans virulence in the nematode Caenorhabditis elegans in an ROS-dependent manner as the rtg1 and rtg3 mutants display virulence defects in wild-type but not in ROS deficient worms. Our findings establish NOX2-derived ROS as a key signal that directs the activity of the pleiotropic fungal regulator Rtg1/3.

Species-specific effects of the introduction of Aspergillus nidulans gfdB in osmophilic aspergilli

Industrial fungi need a strong environmental stress tolerance to ensure acceptable efficiency and yields. Previous studies shed light on the important role that Aspergillus nidulans gfdB, putatively encoding a NAD⁺-dependent glycerol-3-phosphate dehydrogenase, plays in the oxidative and cell wall integrity stress tolerance of this filamentous fungus model organism. The insertion of A. nidulans gfdB into the genome of Aspergillus glaucus strengthened the environmental stress tolerance of this xerophilic/osmophilic fungus, which may facilitate the involvement of this fungus in various industrial and environmental biotechnological processes. On the other hand, the transfer of A. nidulans gfdB to Aspergillus wentii, another promising industrial xerophilic/osmophilic fungus, resulted only in minor and sporadic improvement in environmental stress tolerance and meanwhile partially reversed osmophily. Because A. glaucus and A. wentii are phylogenetically closely related species and both fungi lack a gfdB ortholog, these results warn us that any disturbance of the stress response system of the aspergilli may elicit rather complex and even unforeseeable, species-specific physiological changes. This should be taken into consideration in any future targeted industrial strain development projects aiming at the fortification of the general stress tolerance of these fungi. KEY POINTS: • A. wentii c' gfdB strains showed minor and sporadic stress tolerance phenotypes. • The osmophily of A. wentii significantly decreased in the c' gfdB strains. • Insertion of gfdB caused species-specific phenotypes in A. wentii and A. glaucus.

Genome-Wide Gene Expression Analyses of the AtfA/AtfB-Mediated Menadione Stress Response in Aspergillus nidulans

The bZIP transcription factors (TFs) govern regulation of development, secondary metabolism, and various stress responses in filamentous fungi. In this work, we carried out genome-wide expression studies employing Illumina RNAseq to understand the roles of the two bZIP transcription factors AtfA and AtfB in Aspergillus nidulans. Comparative analyses of transcriptomes of control, ΔatfA, ΔatfB, and ΔatfAΔatfB mutant strains were performed. Dependence of a gene on AtfA (AtfB) was decided by its differential downregulation both between the reference and ΔatfA (ΔatfB) strains and between the ΔatfB (ΔatfA) and the ΔatfAΔatfB strains in vegetatively grown cells (mycelia) and asexual spores (conidia) of menadione sodium bisulfite (MSB)-treated or untreated cultures. As AtfA is the primary bZIP TF governing stress-response in A. nidulans, the number of differentially expressed genes for ΔatfA was significantly higher than for ΔatfB in both mycelial and conidial samples, and most of the AtfB-dependent genes showed AtfA dependence, too. Moreover, the low number of genes depending on AtfB but not on AtfA can be a consequence of ΔatfA leading to downregulation of atfB expression. Conidial samples showed much higher abundance of atfA and atfB mRNAs and more AtfA- and AtfB-affected genes than mycelial samples. In the presence of MSB, the number of AtfB- (but not of AtfA-) affected genes decreased markedly, which was accompanied with decreased mRNA levels of atfB in MSB-treated mycelial (reference strain) and conidial (ΔatfA mutant) samples. In mycelia, the overlap between the AtfA-dependent genes in MSB-treated and in untreated samples was low, demonstrating that distinct genes can be under AtfA control under different conditions. Carbohydrate metabolism genes were enriched in the set of AtfA-dependent genes. Among them, AtfA-dependence of glycolytic genes in conidial samples was the most notable. Levels of transcripts of certain secondary metabolitic gene clusters, such as the Emericellamide cluster, also showed AtfA-dependent regulation. Genes encoding catalase and histidine-containing phosphotransfer proteins showed AtfA-dependence under all experimental conditions. There were 23 AtfB-dependent genes that did not depend on AtfA under any of our experimental conditions. These included a putative α-glucosidase (agdB), a putative α-amylase, calA, which is involved in early conidial germination, and an alternative oxidase. In summary, in A. nidulans there is a complex interaction between the two bZIP transcription factors, where AtfA plays the primary regulatory role.

Genetic Interaction of Global Regulators AflatfA and AflatfB Mediating Development, Stress Response and Aflatoxins B1 Production in Aspergillus flavus

Aspergillus flavus produces carcinogenic and mutagenic aflatoxins, which cause economic losses and risk of food safety by contaminating grains, food and feed. In this study, we characterized two bZIP transcription factors, AflatfA and AflatfB, and their genetic interaction. Compared to the wild type (WT), AflatfA deletion and AflatfA and AflatfB double deletion both caused retarded vegetative growth of mycelia. Relative to WT, the AflatfA deletion strain (ΔAflatfA) and AflatfA and AflatfB double deletion strain (ΔAflatfAΔAflatfB) produced more sclerotia, whereas the AflatfB deletion strain (ΔAflatfB) produced less sclerotia. After 4 °C preservation and incubation at 50 °C, conidia viability dramatically decreased in the ΔAflatfA and ΔAflatfAΔAflatfB but ΔAflatfB mutants, whereas conidia viability of the ΔAflatfAΔAflatfB strain was higher after storage at 4 °C than in AflatfA mutant. Conidia of ΔAflatfA, ΔAflatfB and ΔAflatfAΔAflatfB strains significantly increased in sensitivity to H2O2 in comparison with WT. Compared to WT, the mycelium of ΔAflatfA and ΔAflatfB strains were more sensitive to H2O2; conversely, the ΔAflatfAΔAflatfB strain showed less sensitivity to H2O2. ΔAflatfA and ΔAflatfAΔAflatfB strains displayed less sensitivity to the osmotic reagents NaCl, KCl and Sorbitol, in comparison with WT and ΔAflatfB strains. When on YES medium and hosts corn and peanut, ΔAflatfA and ΔAflatfAΔAflatfB strains produced less aflatoxin B1 (AFB1) than ΔAflatfB, and the AFB1 yield of ΔAflatfB was higher than that of WT. When WT and mutants were inoculated on corn and peanut, the ΔAflatfA and ΔAflatfAΔAflatfB but not ΔAflatfB mutants produced less conidia than did WT. Taken together, this study reveals that AflatfA controls more cellular processes, and the function of AflatfA is stronger than that of AflatfB when of the same process is regulated, except the response to H2O2, which might result from the effect of AflatfA on the transcriptional level of AflatfB.

Transcriptome profiling of transcription factors in Ganoderma lucidum in response to methyl jasmonate

Ganoderma lucidum is a traditional Chinese medicine and its major active ingredients are ganoderma triterpenoids (GTs). To screen for transcription factors (TFs) that involved in the biosynthetic pathway of GTs in G. lucidum, the chemical composition in mycelia, primordium and fruiting body were analyzed, and the transcriptomes of mycelia induced by methyl jasmonate (MeJA) were analyzed. In addition, the expression level data of MeJA-responsive TFs in mycelia, primordia and fruiting body were downloaded from the database, and the correlation analysis was carried out between their expression profiles and the content of total triterpenoids. The results showed that a total of 89 components were identified, and the content of total triterpenoids was the highest in primordium, followed by fruiting body and mycelia. There were 103 differentially expressed TFs that response to MeJA-induction including 95 upregulated and 8 downregulated genes. These TFs were classified into 22 families including C2H2 (15), TFII-related (12), HTH (9), fungal (8), bZIP (6), HMG (5), DADS (2), etc. Correlation analysis showed that the expression level of GL23559 (MADS), GL26472 (HTH), and GL31187 (HMG) showed a positive correlation with the GTs content, respectively. While the expression level of GL25628 (fungal) and GL26980 (PHD) showed a negative correlation with the GTs content, respectively. Furthermore, the over expression of the Glmhr1 gene (GL25628) in Pichia pastoris GS115 indicated that it might be a negative regulator of GT biosynthesis through decreasing the production of lanosterol. This study provided useful information for a better understanding of the regulation of TFs involved in GT biosynthesis and fungal growth in G. lucidum.

GCN4 Enhances the Transcriptional Regulation of AreA by Interacting with SKO1 To Mediate Nitrogen Utilization in Ganoderma lucidum

Nitrogen is an essential nutrient for cell growth and proliferation. Limitations of nitrogen availability in organisms elicit a series of rapid transcriptional reprogramming mechanisms, which involve the participation of many transcription factors.

Potential Original Drug for Aspergillosis: In Vitro and In Vivo Effects of 1-N,N-Dimethylamino-5-Isocyanonaphthalene (DIMICAN) on Aspergillus fumigatus

As the recent outbreak of coronavirus disease 2019 (COVID-19) has shown, viral infections are prone to secondary complications like invasive aspergillosis with a high mortality rate, and therefore the development of novel, effective antifungals is of paramount importance. We have previously demonstrated that 1-amino-5-isocyanonaphthalene (ICAN) derivatives are promising original drug candidates against Candida strains (Patent pending), even against fluconazole resistant C. albicans. Consequently, in this study ICANs were tested on Aspergillus fumigatus, an opportunistic pathogen, which is the leading cause of invasive and systematic pulmonary aspergillosis in immunosuppressed, transplanted and cancer- or COVID-19 treated patients. We have tested several N-alkylated ICANs, a well as 1,5-naphthalene-diisocyanide (DIN) with the microdilution method against Aspergillus fumigatus strains. The results revealed that the diisocyanide (DIN) was the most effective with a minimum inhibitory concentration (MIC) value as low as 0.6 µg mL⁻¹ (3.4 µM); however, its practical applicability is limited by its poor water solubility, which needs to be overcome by proper formulation. The other alkylated derivatives also have in vitro and in vivo anti-Aspergillus fumigatus effects. For animal experiments the second most effective derivative 1-N, N-dimethylamino-5-isocyanonaphthalene (DIMICAN, MIC: 7–8 µg mL⁻¹, 36–41 µM) was selected, toxicity tests were made with mice, and then the antifungal effect of DIMICAN was tested in a neutropenic aspergillosis murine model. Compared to amphotericin B (AMB), a well-known antifungal, the antifungal effect of DIMICAN in vivo turned out to be much better (40% vs. 90% survival after eight days), indicating its potential as a clinical drug candidate.

Functional analysis of the bZIP-type transcription factors AtfA and AtfB in Aspergillus nidulans

Transcription factors (TFs) with the basic leucin zipper domain are key elements of the stress response pathways in filamentous fungi. In this study, we functionally characterized the two bZIP type TFs AtfA and AtfB by deletion (Δ) and overexpression (OE) of their encoding genes in all combination: ΔatfA, ΔatfB, ΔatfAΔatfB, ΔatfAatfBOE, ΔatfBatfAOE, atfAOE, atfBOE and atfAOEatfBOE in Aspergillus nidulans. Based on our previous studies, ΔatfA increased the sensitivity of the fungus to oxidative stress mediated by menadione sodium bisulfite (MSB) and tert-butylhydroperoxide (tBOOH), while ΔatfB was not sensitive to any oxidative stress generating agents, namely MSB, tBOOH and diamide at all. Contrarily, the ΔatfB mutant was sensitive to NaCl, but tolerant to sorbitol. Overexpression of atfB was able to compensate the MSB sensitivity of the ΔatfA mutant. Heavy metal stress elicited by CdCl2 reduced diameter of the atfBOE and atfAOEatfBOE mutant colonies to about 50% of control colony, while the cell wall stress generating agent CongoRed increased the tolerance of the ΔatfA mutant. When we tested the heat stress sensitivity of the asexual spores (conidiospores) of the mutants, we found that conidiospores of ΔatfAatfBOE and ΔatfBatfAOE showed nearly 100% tolerance to heat stress. Asexual development was negatively affected by ΔatfA, while atfAOE and atfAOE coupled with ΔatfB increased the number of conidiospores of the fungus approximately 150% compared to the control. Overexpression of atfB led to a 25% reduction in the number of conidiospores, but increased levels of abaA mRNA and size of conidiospores. Sexual fruiting body (cleistothecium) formation was diminished in the ΔatfA and the ΔatfAΔatfB mutants, while relatively elevated in the ΔatfB and the ΔatfBatfAOE mutants. Production of the mycotoxin sterigmatocystin (ST) was decreased to undetectable levels in the ΔatfA mutant, yet ST production was restored in the ΔatfAΔatfB mutant, suggesting that ΔatfB can suppress ST production defect caused by ΔatfA. Levels of ST were also significantly decreased in the ΔatfAatfBOE, ΔatfBatfAOE and atfAOEatfBOE mutants.

Genome-Wide Identification and Functional Analysis of the bZIP Transcription Factor Family in Rice Bakanae Disease Pathogen, Fusarium fujikuroi

Fungal basic leucine zipper (bZIP) proteins play a vital role in biological processes such as growth, biotic/abiotic stress responses, nutrient utilization, and invasion. In this study, genome-wide identification of bZIP genes in the fungus Fusarium fujikuroi, the pathogen of bakanae disease, was carried out. Forty-four genes encoding bZIP transcription factors (TFs) from the genome of F. fujikuroi (FfbZIP) were identified and functionally characterized. Structures, domains, and phylogenetic relationships of the sequences were analyzed by bioinformatic approaches. Based on the phylogenetic relationships with the FfbZIP proteins of eight other fungi, the bZIP genes can be divided into six groups (A–F). The additional conserved motifs have been identified and their possible functions were predicted. To analyze functions of the bZIP genes, 11 FfbZIPs were selected according to different motifs they contained and were knocked out by genetic recombination. Results of the characteristic studies revealed that these FfbZIPs were involved in oxygen stress, osmotic stress, cell wall selection pressure, cellulose utilization, cell wall penetration, and pathogenicity. In conclusion, this study enhanced understandings of the evolution and regulatory mechanism of the FfbZIPs in fungal growth, abiotic/biotic stress resistance, and pathogenicity, which could be the reference for other fungal bZIP studies.

Drosophila attack inhibits hyphal regeneration and defense mechanisms activation for the fungus Trichoderma atroviride

The capability to respond to wounding is a process shared by organisms of different kingdoms that can result in the regeneration of whole-body parts or lost structures or organs. Filamentous fungi constitute a rich food source that ensures survival and reproduction of their predators and are therefore continuously exposed to mechanical damage. Nevertheless, our understanding of how fungi respond to wounding and predators is scarce. Fungi like plants and animals respond to injury recognizing Damage- and Microbe-Associated Molecular Patterns (DAMPs/MAMPs) that activate Ca2+ and Mitogen-Activated Protein Kinase dependent signaling for the activation of defense mechanisms. During herbivory, plants, in addition to activating pathways related to injury, activate specific responses to combat their predators. Using a transcriptional approach, we studied the capacity of the filamentous fungus Trichoderma atroviride to activate specific responses to injury and attack by different arthropods. Attack by Drosophila melanogaster inhibited the transcriptional activation of genes required for hyphal regeneration, and the fungal innate immune and chemical defense responses. We also provide mechanistic insight of this inhibition involving components of the D. melanogaster salivary glands that repress the expression of a set of genes and block hyphal regeneration.

Distinct and essential roles of bZIP transcription factors in the stress response and pathogenesis in Alternaria alternata

The ability to cope with environmental abiotic stress and biotic stress is crucial for the survival of plants and microorganisms, which enable them to occupy multiple niches in the environment. Previous studies have shown that transcription factors play crucial roles in regulating various biological processes including multiple stress tolerance and response in eukaryotes. This work identified multiple critical transcription factor genes, metabolic pathways and gene ontology (GO) terms related to abiotic stress response were broadly activated by analyzing the transcriptome of phytopathogenic fungus Alternaria alternata under metal ions stresses, oxidative stress, salt stresses, and host-pathogen interaction. We investigated the biological functions and regulatory roles of the bZIP transcriptional factor (TF) genes in the phytopathogenic fungus A.alternata by analyzing targeted gene disrupted mutants. Morphological analysis provides evidence that the bZIP transcription factors (Gcn4, MeaB, Atf1, the ER stress regulator Hac1, and the all development altered-1 gene Ada1) are required for morphogenesis as the colony morphology of these gene deletion mutants was significantly different from that of the wild-type. In addition, bZIPs are involved in the resistance to multiple stresses such as oxidative stress (Ada1, Yap1, MetR) and virulence (Hac1, MetR, Yap1, Ada1) at varying degrees. Transcriptome data demonstrated that the inactivation of bZIPs (Hac1, Atf1, Ada1 and Yap1) significantly affected many genes in multiple critical metabolism pathways and gene ontology (GO) terms. Moreover,the ΔHac1 mutants displayed reduced aerial hypha and are hypersensitivity to endoplasmic reticulum disruptors such as tunicamycin and dithiothreitol. Transcriptome analysis showed that inactivation of Hac1 significantly affected the proteasome process and its downstream unfolded protein binding, indicating that Hac1 participates in the endoplasmic reticulum stress response through the conserved unfolded protein response. Taken together, our findings reveal that bZIP transcription factors function as key regulators of fungal morphogenesis, abiotic stress response and pathogenesis, and expand our understanding of how microbial pathogens utilize these genes to deal with environmental stresses and achieve successful infection in the host plant.

Study on the bZIP-Type Transcription Factors NapA and RsmA in the Regulation of Intracellular Reactive Species Levels and Sterigmatocystin Production of Aspergillus nidulans

Basic leucine zipper (bZIP) transcription factors play a crucial role in the environmental stress response of eukaryotes. In this work, we studied the effect of gene manipulations, including both deletions and overexpressions, of two selected bZIP transcription factors, NapA and RsmA, in the oxidative stress response and sterigmatocystin production of Aspergillus nidulans. We found that NapA was important in the oxidative stress response by negatively regulating intracellular reactive species production and positively regulating catalase activities, whereas RsmA slightly negatively regulated catalase activities. Concerning sterigmatocystin production, the highest concentration was measured in the ΔrsmAΔnapA double deletion mutant, but elevated sterigmatocystin production was also found in the OErsmA OEnapA strain. Our results indicate that NapA influences sterigmatocystin production via regulating reactive species level whereas RsmA modulates toxin production independently of the redox regulation of the cells.

Profiling of chromatin accessibility identifies transcription factor binding sites across the genome of Aspergillus species

Background

The identification of open chromatin regions and transcription factor binding sites (TFBs) is an important step in understanding the regulation of gene expression in diverse species. ATAC-seq is a technique used for such purpose by providing high-resolution measurements of chromatin accessibility revealed through integration of Tn5 transposase. However, the existence of cell walls in filamentous fungi and associated difficulty in purifying nuclei have precluded the routine application of this technique, leading to a lack of experimentally determined and computationally inferred data on the identity of genome-wide cis-regulatory elements (CREs) and TFBs. In this study, we constructed an ATAC-seq platform suitable for filamentous fungi and generated ATAC-seq libraries of Aspergillus niger and Aspergillus oryzae grown under a variety of conditions.

Results

We applied the ATAC-seq assay for filamentous fungi to delineate the syntenic orthologue and differentially changed chromatin accessibility regions among different Aspergillus species, during different culture conditions, and among specific TF-deleted strains. The syntenic orthologues of accessible regions were responsible for the conservative functions across Aspergillus species, while regions differentially changed between culture conditions and TFs mutants drove differential gene expression programs. Importantly, we suggest criteria to determine TFBs through the analysis of unbalanced cleavage of distinct TF-bound DNA strands by Tn5 transposase. Based on this criterion, we constructed data libraries of the in vivo genomic footprint of A. niger under distinct conditions, and generated a database of novel transcription factor binding motifs through comparison of footprints in TF-deleted strains. Furthermore, we validated the novel TFBs in vivo through an artificial synthetic minimal promoter system.

Conclusions

We characterized the chromatin accessibility regions of filamentous fungi species, and identified a complete TFBs map by ATAC-seq, which provides valuable data for future analyses of transcriptional regulation in filamentous fungi.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01114-0.

The impact of bZIP Atf1ortholog global regulators in fungi

Regulation of signal transduction pathways is crucial for the maintenance of cellular homeostasis and organismal development in fungi. Transcription factors are key elements of this regulatory network. The basic-region leucine zipper (bZIP) domain of the bZIP-type transcription factors is responsible for DNA binding while their leucine zipper structural motifs are suitable for dimerization with each other facilitiating the formation of homodimeric or heterodimeric bZIP proteins. This review highlights recent knowledge on the function of fungal orthologs of the Schizosaccharomyces pombe Atf1, Aspergillus nidulans AtfA, and Fusarium verticillioides FvAtfA, bZIP-type transcription factors with a special focus on pathogenic species. We demonstrate that fungal Atf1-AtfA-FvAtfA orthologs play an important role in vegetative growth, sexual and asexual development, stress response, secondary metabolite production, and virulence both in human pathogens, including Aspergillus fumigatus, Mucor circinelloides, Penicillium marneffei, and Cryptococcus neoformans and plant pathogens, like Fusarium ssp., Magnaporthe oryzae, Claviceps purpurea, Botrytis cinerea, and Verticillium dahliae. KEY POINTS: • Atf1 orthologs play crucial role in the growth and development of fungi. • Atf1 orthologs orchestrate environmental stress response of fungi. • Secondary metabolite production and virulence are coordinated by Atf1 orthologs.

AtfA-Independent Adaptation to the Toxic Heavy Metal Cadmium in Aspergillus nidulans

Cadmium is an exceptionally toxic industrial and environmental pollutant classified as a human carcinogen. In order to provide insight into how we can keep our environment safe from cadmium contamination and prevent the accumulation of it in the food chain, we aim to elucidate how Aspergillus nidulans, one of the most abundant fungi in soil, survives and handles cadmium stress. As AtfA is the main transcription factor governing stress responses in A. nidulans, we examined genome-wide expression responses of wild-type and the atfA null mutant exposed to CdCl₂. Both strains showed up-regulation of the crpA Cu²⁺/Cd²⁺ pump gene and AN7729 predicted to encode a putative bis(glutathionato)-cadmium transporter, and transcriptional changes associated with elevated intracellular Cys availability leading to the efficient adaptation to Cd²⁺. Although the deletion of atfA did not alter the cadmium tolerance of the fungus, the cadmium stress response of the mutant differed from that of a reference strain. Promoter and transcriptional analyses of the “Phospho-relay response regulator” genes suggest that the AtfA-dependent regulation of these genes can be relevant in this phenomenon. We concluded that the regulatory network of A. nidulans has a high flexibility allowing the fungus to adapt efficiently to stress both in the presence and absence of this important transcription factor.

Mimosa tenuiflora Aqueous Extract: Role of Condensed Tannins in Anti-Aflatoxin B1 Activity in Aspergillus flavus

Aflatoxin B1 (AFB1) is a potent carcinogenic mycotoxin that contaminates numerous crops pre- and post-harvest. To protect foods and feeds from such toxins without resorting to pesticides, the use of plant extracts has been increasingly studied. The most interesting candidate plants are those with strong antioxidative activity because oxidation reactions may interfere with AFB1 production. The present study investigates how an aqueous extract of Mimosa tenuiflora bark affects both the growth of Aspergillus flavus and AFB1 production. The results reveal a dose-dependent inhibition of toxin synthesis with no impact on fungal growth. AFB1 inhibition is related to a down-modulation of the cluster genes of the biosynthetic pathway and especially to the two internal regulators aflR and aflS. Its strong anti-oxidative activity also allows the aqueous extract to modulate the expression of genes involved in fungal oxidative-stress response, such as msnA, mtfA, atfA, or sod1. Finally, a bio-guided fractionation of the aqueous extract demonstrates that condensed tannins play a major role in the anti-aflatoxin activity of Mimosa tenuiflora bark.

Basic leucine zipper (bZIP) transcription factor genes and their responses to drought stress in ginseng, Panax ginseng C.A. Meyer

BACKGROUND

Ginseng is an important medicinal herb in Asia and Northern America. The basic leucine zipper (bZIP) transcription factor genes play important roles in many biological processes and plant responses to abiotic and biotic stresses, such as drought stress. Nevertheless, the genes remain unknown in ginseng.

RESULTS

Here, we report 91 bZIP genes identified from ginseng, designated PgbZIP genes. These PgbZIP genes were alternatively spliced into 273 transcripts. Phylogenetic analysis grouped the PgbZIP genes into ten groups, including A, B, C, D, E, F, G, H, I and S. Gene Ontology (GO) categorized the PgbZIP genes into five functional subcategories, suggesting that they have diversified in functionality, even though their putative proteins share a number of conserved motifs. These 273 PgbZIP transcripts expressed differentially across 14 tissues, the roots of different ages and the roots of different genotypes. However, the transcripts of the genes expressed coordinately and were more likely to form a co-expression network. Furthermore, we studied the responses of the PgbZIP genes to drought stress in ginseng using a random selection of five PgbZIP genes, including PgbZIP25, PgbZIP38, PgbZIP39, PgbZIP53 and PgbZIP54. The results showed that all five PgbZIP genes responded to drought stress in ginseng, indicating that the PgbZIP genes play important roles in ginseng responses to drought stress.

CONCLUSIONS

These results provide knowledge and gene resources for deeper functional analysis of the PgbZIP genes and molecular tools for enhanced drought tolerance breeding in ginseng.

FvmnSOD is involved in oxidative stress defence, mitochondrial stability and apoptosis prevention in Fusarium verticillioides

Superoxide dismutases are key enzymes in elimination of the superoxide anion radical (O₂ ^•- ) generated intracellularly or by exogenous oxidative stress eliciting agents, like menadione. In this study, we investigated the physiological role of the manganese superoxide dismutase-encoding gene in Fusarium verticillioides via the construction of a gene deletion mutant, ΔFvmnSOD and comparing its phenotype with that of the wild-type parental strain and a ΔFvmnSOD' C strain, complemented with the functional manganese superoxide dismutase gene. Deletion of FvmnSOD had no effect on the relative intracellular superoxide ratio but increased the sensitivity of the fungus to menadione sodium bisulphite on Czapek-Dox stress agar plates. The lack of FvmnSOD caused changes in mitochondrial morphology and physiology: The volumetric ratio of these cell organelles in the second hyphal segment, as well as the total, the KCN-sensitive cytochrome c-dependent and the KCN+SHAM (salicylhidroxamic acid)-resistant residual respiration rates, were higher in the mutant as compared to the wild-type and the complemented strains. Nevertheless, changes in the respiration rates were attributable to the higher volumetric ratio of mitochondria found in the gene deletion mutant. Changes in the mitochondrial functions also brought about higher sensitivity to apoptotic cell death elicited by the Penicillium chrysogenum antifungal protein. The gene deletion mutant developed significantly thinner hyphae in comparison to the wild-type strain. Deletion of FvmnSOD had no effect on fumonisin B₁ and B₂ production of the fungus grown in Myro medium as a static culture.

Two New Aspergillus flavus Reference Genomes Reveal a Large Insertion Potentially Contributing to Isolate Stress Tolerance and Aflatoxin Production

Efforts in genome sequencing in the Aspergillus genus have led to the development of quality reference genomes for several important species including A. nidulans, A. fumigatus, and A. oryzae However, less progress has been made for A. flavus As part of the effort of the USDA-ARS Annual Aflatoxin Workshop Fungal Genome Project, the isolate NRRL3357 was sequenced and resulted in a scaffold-level genome released in 2005. Our goal has been biologically driven, focusing on two areas: isolate variation in aflatoxin production and drought stress exacerbating aflatoxin production by A. flavus Therefore, we developed two reference pseudomolecule genome assemblies derived from chromosome arms for two isolates: AF13, a MAT1-2, highly stress tolerant, and highly aflatoxigenic isolate; and NRRL3357, a MAT1-1, less stress tolerant, and moderate aflatoxin producer in comparison to AF13. Here, we report these two reference-grade assemblies for these isolates through a combination of PacBio long-read sequencing and optical mapping, and coupled them with comparative, functional, and phylogenetic analyses. This analysis resulted in the identification of 153 and 45 unique genes in AF13 and NRRL3357, respectively. We also confirmed the presence of a unique 310 Kb insertion in AF13 containing 60 genes. Analysis of this insertion revealed the presence of a bZIP transcription factor, named atfC, which may contribute to isolate pathogenicity and stress tolerance. Phylogenomic analyses comparing these and other available assemblies also suggest that the species complex of A. flavus is polyphyletic.

Tolerance to alkaline ambient pH in Aspergillus nidulans depends on the activity of ENA proteins

Tolerance of microorganisms to abiotic stress is enabled by regulatory mechanisms that coordinate the expression and activity of resistance genes. Alkalinity and high salt concentrations are major environmental physicochemical stresses. Here, we analyzed the roles of sodium-extrusion family (ENA) transporters EnaA, EnaB and EnaC in the response to these stress conditions in the filamentous fungus Aspergillus nidulans. While EnaC has a minor role, EnaB is a key element for tolerance to Na⁺ and Li⁺ toxicity. Adaptation to alkaline pH requires the concerted action of EnaB with EnaA. Accordingly, expression of enaA and enaB was induced by Na⁺, Li⁺ and pH 8. These expression patterns are altered in a sltAΔ background and completely inhibited in a mutant expressing non-functional PacC protein (palH72). However, a constitutively active PacC form was not sufficient to restore maximum enaA expression. In agreement with their predicted role as membrane ATPases, EnaA localized to the plasma membrane while EnaB accumulated at structures resembling the endoplasmic reticulum. Overall, results suggest different PacC- and SltA-dependent roles for EnaB in pH and salt homeostasis, acting in coordination with EnaA at pH 8 but independently under salt stress.

FvatfA regulates growth, stress tolerance as well as mycotoxin and pigment productions in Fusarium verticillioides

FvatfA from the maize pathogen Fusarium verticillioides putatively encodes the Aspergillus nidulans AtfA and Schizasaccharomyces pombe Atf1 orthologous bZIP-type transcription factor, FvAtfA. In this study, a ΔFvatfA deletion mutant was constructed and then genetically complemented with the fully functional FvatfA gene. Comparing phenotypic features of the wild-type parental, the deletion mutant and the restored strains shed light on the versatile regulatory functions played by FvAtfA in (i) the maintenance of vegetative growth on Czapek-Dox and Potato Dextrose agars and invasive growth on unwounded tomato fruits, (ii) the preservation of conidiospore yield and size, (iii) the orchestration of oxidative (H₂O₂, menadione sodium bisulphite) and cell wall integrity (Congo Red) stress defences and (iv) the regulation of mycotoxin (fumonisins) and pigment (bikaverin, carotenoid) productions. Expression of selected biosynthetic genes both in the fumonisin (fum1, fum8) and the carotenoid (carRA, carB) pathways were down-regulated in the ΔFvatfA strain resulting in defected fumonisin production and considerably decreased carotenoid yields. The expression of bik1, encoding the polyketide synthase needed in bikaverin biosynthesis, was not up-regulated by the deletion of FvatfA meanwhile the ΔFvatfA strain produced approximately ten times more bikaverin than the wild-type or the genetically complemented strains. The abolishment of fumonisin production of the ΔFvatfA strain may lead to the development of new-type, biology-based mycotoxin control strategies. The novel information gained on the regulation of pigment production by this fungus can be interesting for experts working on new, Fusarium-based biomass and pigment production technologies.

Key points • FvatfA regulates vegetative and invasive growths of F. verticillioides. • FvatfA also orchestrates oxidative and cell wall integrity stress defenses. • The ΔFvatfA mutant was deficient in fumonisin production. • FvatfA deletion resulted in decreased carotenoid and increased bikaverin yields.

Genome-Wide Identification and Expression Analysis of the bZIP Transcription Factors in the Mycoparasite Coniothyrium minitans

The basic leucine zipper (bZIP) proteins family is one of the largest and most diverse transcription factors, widely distributed in eukaryotes. However, no information is available regarding the bZIP gene family in Coniothyrium minitans, an important biocontrol agent of the plant pathogen Sclerotinia sclerotiorum. In this study, we identified 34 bZIP genes from the C. minitans genome, which were classified into 8 groups based on their phylogenetic relationships. Intron analysis showed that 28 CmbZIP genes harbored a variable number of introns, and 15 of them shared a feature that intron inserted into the bZIP domain. The intron position in bZIP domain was highly conserved, which was related to recognize the arginine (R) and could be treated as a genomic imprinting. Expression analysis of the CmbZIP genes in response to abiotic stresses indicated that they might play distinct roles in abiotic stress responses. Results showed that 22 CmbZIP genes were upregulated during the later stage of conidial development. Furthermore, transcriptome analysis indicated that CmbZIP genes are involved in different stages of mycoparasitism. Among deletion mutants of four CmbZIPs (CmbZIP07, -09, -13, and -16), only ΔCmbZIP16 mutants significantly reduced its tolerance to the oxidative stress. The other mutants exhibited no significant effects on colony morphology, mycelial growth, conidiation, and mycoparasitism. Taken together, our results suggested that CmbZIP genes play important roles in the abiotic stress responses, conidial development, and mycoparasitism. These results provide comprehensive information of the CmbZIP gene family and lay the foundation for further research on the bZIP gene family regarding their biological functions and evolutionary history.

Supplementation of Aspergillus glaucus with gfdB gene encoding a glycerol 3-phosphate dehydrogenase in Aspergillus nidulans

In Aspergillus nidulans, there are two putative glycerol 3-phosphate dehydrogenases encoded by the genes gfdA and gfdB, while the genome of the osmophilic Aspergillus glaucus harbors only the ortholog of the A. nidulans gfdA gene. Our aim was to insert the gfdB gene into the genome of A. glaucus, and we reached this goal with the adaptation of the Agrobacterium tumefaciens-mediated transformation method. We tested the growth of the gfdB-complemented A. glaucus strains on a medium containing 2 mol l^-1 sorbitol in the presence of oxidative stress generating agents such as tert-butyl hydroperoxide, H₂ O₂ , menadione sodium bisulfite, as well as the cell wall integrity stress-inducing agent Congo Red and the heavy metal stress eliciting CdCl₂ . The growth of the complemented strains was significantly higher than that of the wild-type strain on media supplemented with these stress generating agents. The A. nidulans ΔgfdB mutant was also examined under the same conditions and resulted in a considerably lower growth than that of the control strain in all stress exposure experiments. Our results shed light on the fact that the gfdB gene from A. nidulans was also involved in the stress responses of the complemented A. glaucus strains supporting our hypothesis on the antioxidant function of GfdB in the Aspergilli. Nevertheless, the osmotolerant nature of A. glaucus could not be explained by the lack of the gfdB gene in A. glaucus, as we hypothesized earlier.

Increased Cd2+ biosorption capability of Aspergillus nidulans elicited by crpA deletion

The P-type ATPase CrpA is an important Cu²⁺ /Cd²⁺ pump in the Aspergilli, significantly contributing to the heavy metal stress tolerance of these ascomycetous fungi. As expected, the deletion of crpA resulted in Cu²⁺ /Cd²⁺ -sensitive phenotypes in Aspergillus nidulans on stress agar plates inoculated with conidia. Nevertheless, paradoxical growth stimulations were observed with the ΔcrpA strain in both standard Cu²⁺ stress agar plate experiments and cellophane colony harvest (CCH) cultures, when exposed to Cd²⁺ . These observations reflect efficient compensatory mechanisms for the loss of CrpA operating under these experimental conditions. It is remarkable that the ΔcrpA strain showed a 2.7 times higher Cd biosorption capacity in CCH cultures, which may facilitate the development of new, fungal biomass-based bioremediation technologies to extract harmful Cd²⁺ ions from the environment. The nullification of crpA also significantly changed the spatial distribution of Cu and Cd in CCH cultures, as demonstrated by the combined particle-induced X-ray emission and scanning transmission ion microscopy technique. Most important, the centers of gravity for Cu and Cd accumulations of the ΔcrpA colonies shifted toward the older regions as compared with wild-type surface cultures.

The bZIP Transcription Factor AflRsmA Regulates Aflatoxin B1 Biosynthesis, Oxidative Stress Response and Sclerotium Formation in Aspergillus flavus

Fungal secondary metabolites play important roles not only in fungal ecology but also in humans living as beneficial medicine or harmful toxins. In filamentous fungi, bZIP-type transcription factors (TFs) are associated with the proteins involved in oxidative stress response and secondary metabolism. In this study, a connection between a bZIP TF and oxidative stress induction of secondary metabolism is uncovered in an opportunistic pathogen Aspergillus flavus, which produces carcinogenic and mutagenic aflatoxins. The bZIP transcription factor AflRsmA was identified by a homology research of A. flavus genome with the bZIP protein RsmA, involved in secondary metabolites production in Aspergillusnidulans. The AflrsmA deletion strain (ΔAflrsmA) displayed less sensitivity to the oxidative reagents tert-Butyl hydroperoxide (tBOOH) in comparison with wild type (WT) and AflrsmA overexpression strain (AflrsmA^OE), while AflrsmA^OE strain increased sensitivity to the oxidative reagents menadione sodium bisulfite (MSB) compared to WT and ΔAflrsmA strains. Without oxidative treatment, aflatoxin B₁ (AFB₁) production of ΔAflrsmA strains was consistent with that of WT, but AflrsmA^OE strain produced more AFB₁ than WT; tBOOH and MSB treatment decreased AFB₁ production of ΔAflrsmA compared to WT. Besides, relative to WT, ΔAflrsmA strain decreased sclerotia, while AflrsmA^OE strain increased sclerotia. The decrease of AFB₁ by ΔAflrsmA but increase of AFB₁ by AflrsmA^OE was on corn. Our results suggest that AFB₁ biosynthesis is regulated by AflRsmA by oxidative stress pathways and provide insights into a possible function of AflRsmA in mediating AFB₁ biosynthesis response host defense in pathogen A. flavus.

Diversity of Secondary Metabolism in Aspergillus nidulans Clinical Isolates

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

The Molecular Mechanism of Perillaldehyde Inducing Cell Death in Aspergillus flavus by Inhibiting Energy Metabolism Revealed by Transcriptome Sequencing

Perillaldehyde (PAE), an essential oil in Perilla plants, serves as a safe flavor ingredient in foods, and shows an effectively antifungal activity. Reactive oxygen species (ROS) accumulation in Aspergillus flavus plays a critical role in initiating a metacaspase-dependent apoptosis. However, the reason for ROS accumulation in A. flavus is not yet clear. Using transcriptome sequencing of A. flavus treated with different concentrations of PAE, our data showed that the ROS accumulation might have been as a result of an inhibition of energy metabolism with less production of reducing power. By means of GO and KEGG enrichment analysis, we screened four key pathways, which were divided into two distinct groups: a downregulated group that was made up of the glycolysis and pentose phosphate pathway, and an upregulated group that consisted of MAPK signaling pathway and GSH metabolism pathway. The inhibition of dehydrogenase gene expression in two glycometabolism pathways might play a crucial role in antifungal mechanism of PAE. Also, in our present study, we systematically showed a gene interaction network of how genes of four subsets are effected by PAE stress on glycometabolism, oxidant damage repair, and cell cycle control. This research may contribute to explaining an intrinsic antifungal mechanism of PAE against A. flavus.

rmtA-Dependent Transcriptome and Its Role in Secondary Metabolism, Environmental Stress, and Virulence in Aspergillus flavus

Aspergillus flavus colonizes numerous oil seed crops such as maize, peanuts, treenuts and cottonseed worldwide, contaminating them with aflatoxins and other harmful toxins. Previously our lab characterized the gene rmtA, which encodes an arginine methyltransferase in A. flavus, and demonstrated its role governing the expression of regulators in the aflatoxin gene cluster and subsequent synthesis of toxin. Furthermore, our studies revealed that rmtA also controls conidial and sclerotial development implicating it as an epigenetic regulator in A. flavus. To confirm this, we performed a RNA sequencing analysis to ascertain the extent of rmtA’s influence on the transcriptome of A. flavus. In this analysis we identified over 2000 genes that were rmtA-dependent, including over 200 transcription factor genes, as well as an uncharacterized secondary metabolite gene cluster possibly responsible for the synthesis of an epidithiodiketopiperazine-like compound. Our results also revealed rmtA-dependent genes involved in multiple types of abiotic stress response in A. flavus. Importantly, hundreds of genes active during maize infection were also regulated by rmtA. In addition, in the animal infection model, rmtA was dispensable for virulence, however forced overexpression of rmtA increased mortality with respect to the wild type.

Transcription Factor Atf1 Regulates Expression of Cellulase and Xylanase Genes during Solid-State Fermentation of Ascomycetes

Transcriptional regulation of cellulolytic and xylolytic genes in ascomycete fungi is controlled by specific carbon sources in different external environments. Here, comparative transcriptomic analyses of Penicillium oxalicum grown on wheat bran (WB), WB plus rice straw (WR), or WB plus Avicel (WA) as the sole carbon source under solid-state fermentation (SSF) revealed that most of the differentially expressed genes (DEGs) were involved in metabolism, specifically, carbohydrate metabolism. Of the DEGs, the basic core carbohydrate-active enzyme-encoding genes which responded to the plant biomass resources were identified in P. oxalicum, and their transcriptional levels changed to various extents depending on the different carbon sources. Moreover, this study found that three deletion mutants of genes encoding putative transcription factors showed significant alterations in filter paper cellulase production compared with that of a parental P. oxalicum strain with a deletion of Ku70 (ΔPoxKu70 strain) when grown on WR under SSF. Importantly, the ΔPoxAtf1 mutant (with a deletion of P. oxalicum Atf1, also called POX03016) displayed 46.1 to 183.2% more cellulase and xylanase production than a ΔPoxKu70 mutant after 2 days of growth on WR. RNA sequencing and quantitative reverse transcription-PCR revealed that PoxAtf1 dynamically regulated the expression of major cellulase and xylanase genes under SSF. PoxAtf1 bound to the promoter regions of the key cellulase and xylanase genes in vitro This study provides novel insights into the regulatory mechanism of fungal cellulase and xylanase gene expression under SSF.IMPORTANCE The transition to a more environmentally friendly economy encourages studies involving the high-value-added utilization of lignocellulosic biomass. Solid-state fermentation (SSF), that simulates the natural habitat of soil microorganisms, is used for a variety of applications such as biomass biorefinery. Prior to the current study, our understanding of genome-wide gene expression and of the regulation of gene expression of lignocellulose-degrading enzymes in ascomycete fungi during SSF was limited. Here, we employed RNA sequencing and genetic analyses to investigate transcriptomes of Penicillium oxalicum strain EU2101 cultured on medium containing different carbon sources and to identify and characterize transcription factors for regulating the expression of cellulase and xylanase genes during SSF. The results generated will provide novel insights into genetic engineering of filamentous fungi to further increase enzyme production.

General stress response or adaptation to rapid growth in Aspergillus nidulans?

Genome-wide transcriptional changes in Aspergillus nidulans induced by nine different stress conditions were evaluated to reveal the general environmental stress response gene set showing unidirectional expressional changes under various types of stress. Clustering the genes by their transcriptional changes was a useful technique for identifying large groups of co-regulated genes. Altogether, 1642 co-upregulated and 3916 co-downregulated genes were identified. Nevertheless, the co-regulated genes describe the difference between the transcriptomes recorded under the stress conditions tested and one chosen reference culture condition which is designated as the "unstressed" condition. Obviously, the corresponding transcriptional differences may be attributed to either the general stress response or the reference condition. Accordingly, reduced growth and increased transcription of certain antioxidative enzymes observed under stress may be interpreted as elements of the general stress response or as a feature of the "optimal growth" reference condition and decreased antioxidative protection due to "rapid growth" stress. Reversing the many to one comparison underlying the identification of co-regulated gene sets allows the same procedure to highlight changes under a single condition with respect to a set of other "background" conditions. As an example, we compared menadione treatment to our other conditions and identified downregulation of endoplasmic reticulum dependent processes and upregulation of iron-sulfur cluster assembly as well as glutathione-S-transferase genes as changes characteristic of MSB-treated cultures. Deletion of the atfA gene markedly altered the co-regulated gene sets primarily by changing the reference transcriptome; not by changing the stress responsiveness of genes. The functional characterization of AtfA-dependent co-regulated genes demonstrated the involvement of AtfA in the regulation of both vegetative growth and conidiogenesis in untreated cultures. Our data also suggested that the diverse effects of atfA gene deletion on the transcriptome under different stress conditions were the consequence of the altered transcription of several phosphorelay signal transduction system genes, including fphA, nikA, phkA, srrB, srrC, sskA and tcsB. Hopefully, this study will draw further attention to the importance of the proper selection of reference cultures in fungal transcriptomics studies especially when elements of specific stress responses are mapped.

Characterization of gfdB, putatively encoding a glycerol 3-phosphate dehydrogenase in Aspergillus nidulans

The genome of Aspergillus nidulans accommodates two glycerol 3-phosphate dehydrogenase genes, gfdA and gfdB. Previous studies confirmed that GfdA is involved in the osmotic stress defence of the fungus. In this work, the physiological role of GfdB was characterized via the construction and functional characterization of the gene deletion mutant ΔgfdB. Unexpectedly, ΔgfdB strains showed oxidative stress sensitivity in the presence of a series of well-known oxidants including tert-butyl-hydroperoxide (tBOOH), diamide as well as hydrogen peroxide. Moderate sensitivity of the mutant towards the cell wall stress inducing agent CongoRed was also observed. Hence, both Gfd isoenzymes contributed to the environmental stress defence of the fungus but their functions were stress-type-specific. Furthermore, the specific activities of certain antioxidant enzymes, like catalase and glutathione peroxidase, were lower in ΔgfdB hyphae than those recorded in the control strain. As a consequence, mycelia from ΔgfdB cultures accumulated reactive species at higher levels than the control. On the other hand, the specific glutathione reductase activity was higher in the mutant, most likely to compensate for the elevated intracellular oxidative species concentrations. Nevertheless, the efficient control of reactive species failed in ΔgfdB cultures, which resulted in reduced viability and, concomitantly, early onset of programmed cell death in mutant hyphae. Inactivation of gfdB brought about higher mannitol accumulation in mycelia meanwhile the erythritol production was not disturbed in unstressed cultures. After oxidative stress treatment with tBOOH, only mannitol was detected in both mutant and control mycelia and the accumulation of mannitol even intensified in the ΔgfdB strain.

Effect of amphotericin B and voriconazole on the outgrowth of conidia of Aspergillus fumigatus followed by time-lapse microscopy

Studies of morphological measurements from the outgrowth of cells to a network of hyphae have been extended from Candida albicans (Nagy et al. in Appl Microbiol Biotechnol 98(11):5185-5194. https://doi.org/10.1007/s00253-014-5696-5 , 2014) to invasive conidiospores of Aspergillus fumigatus upon treatment with antifungal agents. The understanding of mycelial processes is important to optimize industrial processes such as fermentation and contributes to the fight against pathogenic fungi. This brief study combines TLS with digital image analysis. The TLS system was adapted to get information related to the adherence and growth dynamics of filamentous fungi. This approach was used earlier to distinguish among subphases of bacterial and fungal infections of mammal cells by detecting Mycoplasma infection in cell cultures causing serious damages in cell cultures. We describe changes in adherence, germination of spores, and hyphal growth of A. fumigatus, taking place in the absence and presence of amphotericin B (AMB) and voriconazole (VRC). These growth parameters were measured by TLS in CO₂ incubator under physiological Photomicrography by TLS and extended for a longer period of time up to several weeks combined with image analysis represents a comfortable and reliable means to characterize the growth dynamism of A. fumigatus. The most important observation of medical importance related to the pathomechanism of VRC was that it did not adhere to conidiospores, i.e. that it did not contribute to the attachment of spores to the growth surface, and did not prevent germination but delayed hypha protrusion and elongation. In contrast AMB adhered to conidia, inhibited germination, hypha elongation and branching. It was concluded that AMB was efficient against the therapy of growth but not against the prevention of fungal infection.

Oxidative stress and the presence of bacteria increase gene expression of the antimicrobial peptide aclasin, a fungal CSαβ defensin in Aspergillus clavatus

Background

Antimicrobial peptides (AMPs) represent a broad class of naturally occurring antimicrobial compounds. Plants, invertebrates and fungi produce various AMPs as, for example, defensins. Most of these defensins are characterised by the presence of a cysteine-stabilised α-helical and β-sheet (CSαβ) motif. The changes in gene expression of a fungal CSαβ defensin by stress conditions were investigated in Aspergillus clavatus. A. clavatus produces the CSαβ defensin Aclasin, which is encoded by the aclasin gene.

Methods

Aclasin expression was evaluated in submerged mycelium cultures under heat shock, osmotic stress, oxidative stress and the presence of bacteria by quantitative real-time PCR.

Results

Aclasin expression increased two fold under oxidative stress conditions and in the presence of viable and heat-killed Bacillus megaterium. Under heat shock and osmotic stress, aclasin expression decreased.

Discussion

The results suggest that oxidative stress and the presence of bacteria might regulate fungal defensin expression. Moreover, fungi might recognise microorganisms as plants and animals do.

Fungal Stress Database (FSD)--a repository of fungal stress physiological data

The construction of the Fungal Stress Database (FSD) was initiated and fueled by two major goals. At first, some outstandingly important groups of filamentous fungi including the aspergilli possess remarkable capabilities to adapt to a wide spectrum of environmental stress conditions but the underlying mechanisms of this stress tolerance have remained yet to be elucidated. Furthermore, the lack of any satisfactory interlaboratory standardization of stress assays, e.g. the widely used stress agar plate experiments, often hinders the direct comparison and discussion of stress physiological data gained for various fungal species by different research groups. In order to overcome these difficulties and to promote multilevel, e.g. combined comparative physiology-based and comparative genomics-based, stress research in filamentous fungi, we constructed FSD, which currently stores 1412 photos taken on Aspergillus colonies grown under precisely defined stress conditions. This study involved altogether 18 Aspergillus strains representing 17 species with two different strains for Aspergillus niger and covered six different stress conditions. Stress treatments were selected considering the frequency of various stress tolerance studies published in the last decade in the aspergilli and included oxidative (H₂O₂, menadione sodium bisulphite), high-osmolarity (NaCl, sorbitol), cell wall integrity (Congo Red) and heavy metal (CdCl₂) stress exposures. In the future, we would like to expand this database to accommodate further fungal species and stress treatments.

URL: http://www.fung-stress.org/

Duplications and losses of genes encoding known elements of the stress defence system of the Aspergilli contribute to the evolution of these filamentous fungi but do not directly influence their environmental stress tolerance

The contribution of stress protein duplication and deletion events to the evolution of the Aspergilli was studied. We performed a large-scale homology analysis of stress proteins and generated and analysed three stress defence system models based on Saccharomyces cerevisiae, Schizosaccharomyces pombe and Aspergillus nidulans. Although both yeast-based and A. nidulans-based models were suitable to trace evolutionary changes, the A. nidulans-based model performed better in mapping stress protein radiations. The strong Mantel correlation found between the positions of species in the phylogenetic tree on the one hand and either in the A. nidulans-based or S. cerevisiae-based models on the other hand demonstrated that stress protein expansions and reductions contributed significantly to the evolution of the Aspergilli. Interestingly, stress tolerance attributes correlated well with the number of orthologs only for a few stress proteins. Notable examples are Ftr1 iron permease and Fet3 ferro-O2-oxidoreductase, elements of the reductive iron assimilation pathway, in the S. cerevisiae-based model, as well as MpkC, a HogA-like mitogen activated protein kinase in the A. nidulans-based model. In the case of the iron assimilation proteins, the number of orthologs showed a positive correlation with H2O2-induced stress tolerance while the number of MpkC orthologs correlated positively with Congo Red induced cell wall stress, sorbitol induced osmotic stress and H2O2 induced oxidative stress tolerances. For most stress proteins, changes in the number of orthologs did not correlate well with any stress tolerance attributes. As a consequence, stress tolerance patterns of the studied Aspergilli did not correlate with either the sets of stress response proteins in general or with the phylogeny of the species studied. These observations suggest that stress protein duplication and deletion events significantly contributed to the evolution of stress tolerance attributes of Aspergilli. In contrast, there are other processes, which may counterbalance the effects of stress gene duplications or deletions including (i) alterations in the structures of stress proteins leading to changes in their biological activities, (ii) varying biosynthesis of stress proteins, (iii) rewiring stress response regulatory networks or even (iv) acquiring new stress response genes by horizontal gene transfer. All these multilevel changes are indispensable for the successful adaptation of filamentous fungi to altering environmental conditions, especially when these organisms are entering new ecological niches.

The Antioxidant Gallic Acid Inhibits Aflatoxin Formation in Aspergillus flavus by Modulating Transcription Factors FarB and CreA

Aflatoxin biosynthesis is correlated with oxidative stress and is proposed to function as a secondary defense mechanism to redundant intracellular reactive oxygen species (ROS). We find that the antioxidant gallic acid inhibits aflatoxin formation and growth in Aspergillus flavus in a dose-dependent manner. Global expression analysis (RNA-Seq) of gallic acid-treated A. flavus showed that 0.8% (w/v) gallic acid revealed two possible routes of aflatoxin inhibition. Gallic acid significantly inhibited the expression of farB, encoding a transcription factor that participates in peroxisomal fatty acid β-oxidation, a fundamental contributor to aflatoxin production. Secondly, the carbon repression regulator encoding gene, creA, was significantly down regulated by gallic acid treatment. CreA is necessary for aflatoxin synthesis, and aflatoxin biosynthesis genes were significantly downregulated in ∆creA mutants. In addition, the results of antioxidant enzyme activities and the lipid oxidation levels coupled with RNA-Seq data of antioxidant genes indicated that gallic acid may reduce oxidative stress through the glutathione- and thioredoxin-dependent systems in A. flavus.

Murine model to follow hyphal development in invasive pulmonary aspergillosis

Aspergillus fumigatus is an opportunistic pathogen, the leading cause of invasive and disseminated aspergillosis in systemic immunocompromised patients, and an important cause of mortality. The aim of the present study was to adapt a pulmonary aspergillosis murine model, to determine pathodynamical parameters quantitatively, and to follow the progression of fungal infection in vivo. The nasal inoculation of Aspergillus conidia in mice previously subjected to immunosuppression with cyclophosphamide (CP) turned out to be a more suitable model than that of immunosuppressed with hydrocortisone (HC). The following parameters were found to correlate quantitatively with the progress of the infection: (i) survival rate, (ii) weight loss of mice, (iii) infected focal plaque size, (iv) hyphal density, (v) hyphal length distribution of A. fumigatus, and the (vi) the histopathological status and scores. These parameters will be essential elements for the development of antifungal drugs and therapies, and important for the investigation of the pathogenicity in different strains of A. fumigatus.

Genome-wide identification and analysis of the basic leucine zipper (bZIP) transcription factor gene family in Ustilaginoidea virens

The basic leucine zipper (bZIP) transcription factor (TF) family is one of the largest and most diverse TF families widely distributed across the eukaryotes. The bZIP TF family plays an important role in growth, development, and response to abiotic or biotic stresses, which have been well characterized in plants, but not in plant pathogenic fungi. In this study, we performed genome-wide and systematic bioinformatics analysis of bZIP genes in the fungus Ustilaginoidea virens, the causal agent of rice false smut disease. We identified 28 bZIP family members in the U. virens genome by searching for the bZIP domain in predicted genes. The gene structures, motifs, and phylogenetic relationships were analyzed for bZIP genes in U. virens (UvbZIP). Together with bZIP proteins from two other fungi, the bZIP genes can be divided into eight groups according to their phylogenetic relationships. Based on RNA-Seq data, the expression profiles of UvbZIP genes at different infection stages were evaluated. Results showed that 17 UvbZIP genes were up-regulated during the infection period. Furthermore, 11 infection-related UvbZIP genes were investigated under H2O2 stress and the expression level of eight genes were changed, which confirmed their role in stress tolerance and pathogenicity. In summary, our genome-wide systematic characterization and expression analysis of UvbZIP genes provided insight into the molecular function of these genes in U. virens and provides a reference for other pathogens.

Transcriptome-Based Modeling Reveals that Oxidative Stress Induces Modulation of the AtfA-Dependent Signaling Networks in Aspergillus nidulans

To better understand the molecular functions of the master stress-response regulator AtfA in Aspergillus nidulans, transcriptomic analyses of the atfA null mutant and the appropriate control strains exposed to menadione sodium bisulfite- (MSB-), t-butylhydroperoxide- and diamide-induced oxidative stresses were performed. Several elements of oxidative stress response were differentially expressed. Many of them, including the downregulation of the mitotic cell cycle, as the MSB stress-specific upregulation of FeS cluster assembly and the MSB stress-specific downregulation of nitrate reduction, tricarboxylic acid cycle, and ER to Golgi vesicle-mediated transport, showed AtfA dependence. To elucidate the potential global regulatory role of AtfA governing expression of a high number of genes with very versatile biological functions, we devised a model based on the comprehensive transcriptomic data. Our model suggests that an important function of AtfA is to modulate the transduction of stress signals. Although it may regulate directly only a limited number of genes, these include elements of the signaling network, for example, members of the two-component signal transduction systems. AtfA acts in a stress-specific manner, which may increase further the number and diversity of AtfA-dependent genes. Our model sheds light on the versatility of the physiological functions of AtfA and its orthologs in fungi.

Functional characterization of two bZIP transcription factors in Verticillium dahliae

bZIP transcription factors play various biological roles in stress responses, conidiation, and pathogenicity in pathogenic fungi. Here, we report two bZIP transcription factors (VDAG_08640 and VDAG_08676) of Verticillium dahliae, which were differentially expressed during microsclerotia development and induced by hydrogen peroxide as well. We find that deletion of either gene does not affect microsclerotia formation and the sensitivity to hydrogen peroxide; however, the mutants manifest decreased activity of extracellular peroxidase and laccase. Other phenotypic characterization reveals that VDAG_08676 disruption results in significant reduction of conidial production and virulence, while VDAG_08640 disruption does not lead to observable phenotypic variances compared with the wild-type strain. To elucidate whether they exhibit functional redundancy, double deletion mutants were generated. The double deletion mutants show remarkably increased sensitivity to hydrogen peroxide stress, whereas the two genes are not involved in microsclerotia formation. Taken together, our data demonstrate that a bZIP transcription factor gene VDAG_08676 is involved in the conidial production, oxidative stress response and virulence which may lay a foundation for further analysis of other bZIP transcription factors in V. dahliae.

Genome-Wide Identification of bZIP Family Genes Involved in Drought and Heat Stresses in Strawberry (Fragaria vesca)

Basic leucine zipper (bZIP) genes are known to play a crucial role in response to various processes in plant as well as abiotic or biotic stress challenges. We have performed an identification and characterization of 50 bZIP genes across the woodland strawberry (Fragaria vesca) genome, which were divided into 10 clades according to the phylogenetic relationship of the strawberry bZIP proteins with those in Arabidopsis and rice. Five categories of intron patterns were observed within basic and hinge regions of the bZIP domains. Some additional conserved motifs have been found with the group specificity. Further, we predicted DNA-binding specificity of the basic and hinge regions as well as dimerization properties of leucine zipper regions, which was consistent with our phylogenetic clade and classified into 20 subfamilies. Across the different developmental stages of 15 organs and two types of fruits, the clade A bZIP members showed different tissue-specific expression patterns and the duplicated genes were differentially regulated, indicating a functional diversification coupled with the expansion of this gene family in strawberry. Under normal growth conditions, mrna11837 and mrna30280 of clade A showed very weak expression levels in organs and fruits, respectively; but higher expression was observed with different set of genes following drought and heat treatment, which may be caused by the separate response pathway between drought and heat treatments.

Genomics of Compensatory Adaptation in Experimental Populations of Aspergillus nidulans

Knowledge of the number and nature of genetic changes responsible for adaptation is essential for understanding and predicting evolutionary trajectories. Here, we study the genomic basis of compensatory adaptation to the fitness cost of fungicide resistance in experimentally evolved strains of the filamentous fungus Aspergillus nidulans The original selection experiment tracked the fitness recovery of lines founded by an ancestral strain that was resistant to fludioxonil, but paid a fitness cost in the absence of the fungicide. We obtained whole-genome sequence data for the ancestral A. nidulans strain and eight experimentally evolved strains. We find that fludioxonil resistance in the ancestor was likely conferred by a mutation in histidine kinase nikA, part of the two-component signal transduction system of the high-osmolarity glycerol (HOG) stress response pathway. To compensate for the pleiotropic negative effects of the resistance mutation, the subsequent fitness gains observed in the evolved lines were likely caused by secondary modification of HOG pathway activity. Candidate genes for the compensatory fitness increases were significantly overrepresented by stress response functions, and some were specifically associated with the HOG pathway itself. Parallel evolution at the gene level was rare among evolved lines. There was a positive relationship between the predicted number of adaptive steps, estimated from fitness data, and the number of genomic mutations, determined by whole-genome sequencing. However, the number of genomic mutations was, on average, 8.45 times greater than the number of adaptive steps inferred from fitness data. This research expands our understanding of the genetic basis of adaptation in multicellular eukaryotes and lays out a framework for future work on the genomics of compensatory adaptation in A. nidulans.

Genome-wide transcriptome analysis of Aspergillus fumigatus exposed to osmotic stress reveals regulators of osmotic and cell wall stresses that are SakAHOG1 and MpkC dependent

Invasive aspergillosis is predominantly caused by Aspergillus fumigatus, and adaptations to stresses experienced within the human host are a prerequisite for the survival and virulence strategies of the pathogen. The central signal transduction pathway operating during hyperosmotic stress is the high osmolarity glycerol mitogen-activated protein kinase cascade. A. fumigatus MpkC and SakA, orthologues of the Saccharomyces cerevisiae Hog1p, constitute the primary regulator of the hyperosmotic stress response. We compared A. fumigatus wild-type transcriptional response to osmotic stress with the ΔmpkC, ΔsakA, and ΔmpkC ΔsakA strains. Our results strongly indicate that MpkC and SakA have independent and collaborative functions during the transcriptional response to transient osmotic stress. We have identified and characterized null mutants for four A. fumigatus basic leucine zipper proteins transcription factors. The atfA and atfB have comparable expression levels with the wild-type in ΔmpkC but are repressed in ΔsakA and ΔmpkC ΔsakA post-osmotic stress. The atfC and atfD have reduced expression levels in all mutants post-osmotic stress. The atfA-D null mutants displayed several phenotypes related to osmotic, oxidative, and cell wall stresses. The ΔatfA and ΔatfB were shown to be avirulent and to have attenuated virulence, respectively, in both Galleria mellonella and a neutropenic murine model of invasive pulmonary aspergillosis.

Genome-Wide Chromatin Immunoprecipitation Sequencing Analysis of the Penicillium chrysogenum Velvet Protein PcVelA Identifies Methyltransferase PcLlmA as a Novel Downstream Regulator of Fungal Development

Penicillium chrysogenum is the sole industrial producer of the β-lactam antibiotic penicillin, which is the most commonly used drug for treating bacterial infections. In P. chrysogenum and other filamentous fungi, secondary metabolism and morphogenesis are controlled by the highly conserved multisubunit velvet complex. Here we present the first chromatin immunoprecipitation next-generation sequencing (ChIP-seq) analysis of a fungal velvet protein, providing experimental evidence that a velvet homologue in P. chrysogenum (PcVelA) acts as a direct transcriptional regulator at the DNA level in addition to functioning as a regulator at the protein level in P. chrysogenum, which was previously described. We identified many target genes that are related to processes known to be dependent on PcVelA, e.g., secondary metabolism as well as asexual and sexual development. We also identified seven PcVelA target genes that encode putative methyltransferases. Yeast two-hybrid and bimolecular fluorescence complementation analyses showed that one of the putative methyltransferases, PcLlmA, directly interacts with PcVelA. Furthermore, functional characterization of PcLlmA demonstrated that this protein is involved in the regulation of conidiosporogenesis, pellet formation, and hyphal morphology, all traits with major biotechnological relevance. IMPORTANCE Filamentous fungi are of major interest for biotechnological and pharmaceutical applications. This is due mainly to their ability to produce a wide variety of secondary metabolites, many of which are relevant as antibiotics. One of the most prominent examples is penicillin, a β-lactam antibiotic that is produced on the industrial scale by fermentation of P. chrysogenum. In recent years, the multisubunit protein complex velvet has been identified as one of the key regulators of fungal secondary metabolism and development. However, until recently, only a little has been known about how velvet mediates regulation at the molecular level. To address this issue, we performed ChIP-seq (chromatin immunoprecipitation in combination with next-generation sequencing) on and follow-up analysis of PcVelA, the core component of the velvet complex in P. chrysogenum. We demonstrate direct involvement of velvet in transcriptional control and present the putative methyltransferase PcLlmA as a new downstream factor and interaction partner of PcVelA.