Next Article in Journal
Microbiota Modulates the Immunomodulatory Effects of Filifolinone on Atlantic Salmon
Next Article in Special Issue
Colistin and Isavuconazole Interact Synergistically In Vitro against Aspergillus nidulans and Aspergillus niger
Previous Article in Journal
Bacterial Populations in International Artisanal Kefirs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome Wide Analysis Reveals the Role of VadA in Stress Response, Germination, and Sterigmatocystin Production in Aspergillus nidulans Conidia

1
School of Food Science and Biotechnology, Kyungpook National University, Daegu 41566, Korea
2
Department of Integrative Biology, Kyungpook National University, Daegu 41566, Korea
*
Author to whom correspondence should be addressed.
Microorganisms 2020, 8(9), 1319; https://doi.org/10.3390/microorganisms8091319
Submission received: 22 July 2020 / Revised: 27 August 2020 / Accepted: 28 August 2020 / Published: 30 August 2020
(This article belongs to the Special Issue Aspergillus and Health 1.0)

Abstract

:
In the Aspergillus species, conidia are asexual spores that are infectious particles responsible for propagation. Conidia contain various mycotoxins that can have detrimental effects in humans. Previous study demonstrated that VadA is required for fungal development and spore viability in the model fungus Aspergillus nidulans. In the present study, vadA transcriptomic analysis revealed that VadA affects the mRNA expression of a variety of genes in A. nidulans conidia. The genes that were primarily affected in conidia were associated with trehalose biosynthesis, cell-wall integrity, stress response, and secondary metabolism. Genetic changes caused by deletion of vadA were related to phenotypes of the vadA deletion mutant conidia. The deletion of vadA resulted in increased conidial sensitivity against ultraviolet stress and induced germ tube formation in the presence and absence of glucose. In addition, most genes in the secondary metabolism gene clusters of sterigmatocystin, asperfuranone, monodictyphenone, and asperthecin were upregulated in the mutant conidia with vadA deletion. The deletion of vadA led to an increase in the amount of sterigmatocystin in the conidia, suggesting that VadA is essential for the repression of sterigmatocystin production in conidia. These results suggest that VadA coordinates conidia maturation, stress response, and secondary metabolism in A. nidulans conidia.

1. Introduction

The Aspergillus species produce many spores for propagation under a variety of environmental conditions [1,2]. Most species reproduce asexual spores, known as conidia, which are light and small and can be easily dispersed in the environment [3,4]. These conidia can be the main infectious agent in humans, and can lead to aspergillosis in immunocompromised patients [5], or to allergic bronchopulmonary aspergillosis [6]. Additionally, fungal conidia contain many mycotoxins, including aflatoxins and sterigmatocystin, which can also lead to carcinogenesis [7]. Conidia have a hard cell wall with a protective layer in order to protect against a variety of aggressive environments, such as osmotic, oxidative, temperature, and ultraviolet (UV) stress [8]. During spore wall formation, a myriad of spore-specific regulators alter the composition of polysaccharides in the cell wall, including α-1,3 glucan, β-1,3 glucan, and chitin [9,10,11,12]. Research assessing cell wall composition has primarily focused on Aspergillus species, such as Aspergillus nidulans and A. fumigatus [13,14,15].
A. nidulans is the model fungus for asexual development, as it produces asexual specialized structures, called conidiophores, which bear conidia [16,17]. A conidiophore is a multicellular structure consisting of a vesicle, metulae, phialides, and conidia, and the formation of conidiophores is regulated by three key transcription factors, including BrlA, AbaA, and WetA [2,18,19]. Among these transcription factors, WetA is a key regulator that plays a central role in conidia formation during the late phase of asexual development [20,21]. Two velvet regulators, VosA and VelB, are involved in conidial formation, maturation, and dormancy in Aspergillus species [22,23]. These three transcription factors regulate the mRNA expression a thousand conidia-specific genes [21,24,25]. A variety of studies have focused on these three transcription factors; however, other spore-specific genes have not been investigated in detail.
vadA is one of the VosA/VelB-activated developmental genes that is specifically expressed during formation of conidia in A. nidulans [26]. Previous studies have shown that vadA is required for the balance between asexual and sexual development in A. nidulans. VadA has also been shown to be involved in conidial viability, conidial stress response, β-glucan biosynthesis, and sterigmatocystin production, suggesting that VadA is a multi-functional regulator in A. nidulans. Although the roles of VadA have been investigated, detailed regulatory mechanisms of VadA have not been elucidated. In the present study, a transcriptome analysis identified that deletion of vadA leads to significant alterations in the mRNA levels of more than 3000 genes in the conidia of A. nidulans. Several secondary metabolite gene clusters and genes associated with polysaccharide metabolic processes are affected by the loss of vadA, suggesting that VadA is a key regulator of conidia formation in A. nidulans.

2. Materials and Methods

2.1. Strains and Media

Control (THS30, pyrG89; AfupyrG+; veA+) [24], vadA deletion mutant (ΔvadA, THS33, pyrG89; pyroA4; ΔvadA::AfupyrG+; veA+), and vadA complemented (C’ vadA, THS34, pyrG89; pyroA::vadA(p)::vadA::FLAG3x::pyroA; ΔvadA::AfupyrG+; veA+) strains [26] were grown on solid minimal medium (MM) with 1% glucose (MMG) [27].

2.2. RNA-Sequencing (RNA-Seq) Analysis

Fungal conidia were inoculated onto solid MMG and incubated at 37 °C for 2 days. Conidia were collected from the plates after the two days of incubation, washed using ddH2O, and filtrated through Miracloth (Calbiochem, San Diego, CA, USA). Conidia were then collected for total RNA extraction. Total RNA was extracted from fungal conidia using Trizol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions, with some modifications. RNA samples were treated with DNase I (Promega, Madison, WI, USA) to remove DNA contamination, and then purified using the RNeasy Mini Kit (Qiagen, Germantown, MD, USA). Library preparation and sequencing was performed by Theragen Bio (Suwon, South Korea).
RNA quality was assessed via analysis of the rRNA band integrity using the Agilent RNA 6000 Nano Kit (Agilent Technologies, Santa Clara, CA, USA). Before cDNA library construction, 2 µg of total RNA and magnetic beads with oligos (dT) were used for poly (A) enrichment. The purified mRNA was disrupted into short fragments and double-stranded cDNAs were immediately synthesized. The cDNAs were subjected to end-repair, poly (A) addition, and connected with sequencing adapters using the TruSeq Stranded mRNA Sample Prep Kit (Illumina, San Diego, CA, USA). Suitable fragments were automatically purified using the BluePippin 2% agarose gel cassette (Sage Science, Beverly, MA, USA), and were selected as templates for PCR amplification. The final library sizes and qualities were evaluated electrophoretically using the Agilent High Sensitivity DNA Kit (Agilent Technologies, Santa Clara, CA, USA). Fragments were between 350–450 bp. The library was sequenced using an Illumina HiSeq2500 sequencer (Illumina, San Diego, CA, USA).2.3. RNA-seq data analysis
RNA-seq data analysis and the filtering process were performed as described previously [28]. Filtered reads were mapped to the A. nidulans A4 transcriptome [29] using the aligner STAR v.2.3.0e software [30]. Gene expression levels were measured using Cufflinks v2.1.1 [31] with the gene annotation database from the Aspergillus Genome Database (AspGD) [29]. For differential expression analysis, gene level count data were generated using the HTSeq-count v0.5.4p3 [32] tool with the option “-m intersection-nonempty” and “-r option considering paired-end sequence.” Based on the read count data that was calculated, differentially expressed genes (DEG) were identified using the R package TCC [33]. Normalization factors were calculated using the iterative DEGES/edgeR method. DEGs were identified based on a q-value threshold less than 0.05.

2.3. Gene Ontology Enrichment Analysis

Gene ontology enrichment analyses were carried out using at AspGD [32] and FungiFun [34] database. To characterize the genes identified from the DEG analysis, gene ontology (GO)-based trend tests were performed using the Fisher’s exact test. p-values < 0.001 were considered statistically significant.

2.4. Data Availability

All RNA-seq data files are available from the NCBI BioProject database (PRJNA646843).

2.5. UV Stress Tolerance Test

The UV stress tolerance test was carried out as described previously [35]. Two-day old conidia (about 100 conidia per plate) were spread on MMG plates. The plates were UV -irradiated using a UV crosslinker (Spectrolinke XL-1000 UV crosslinker, Thomas Scientific, Swedesboro, NJ, USA), and the irradiated plates were further incubated at 37 °C for 48 h. The colony numbers were counted and calculated as the survival ratio to that of the untreated control, in triplicates.

2.6. Conidial Germination Test

The germination of conidia was measured as described previously [36]. Briefly, conidia (107 per plate) were spread onto solid MM (containing 0.6% (w/v) sodium nitrate as nitrogen source) with or without glucose, and incubated at 37 °C. The germination rates of the conidia were assessed every hour after inoculation using a Zeiss Lab.A1 microscope (Jena, Germany) equipped with AxioCam 105 and AxioVision (Rel. 4.9) digital imaging software.

2.7. Sterigmatocystin Extraction from Conidia and Thin-Layer Chromatography (TLC) Analysis

Sterigmatocystin was extracted from conidia as described previously, with some modifications [28,37]. Conidia were cultured for 5 days and 109 conidia were mixed with 2 mL of CHCl3 and 0.5 mm Zirconia/Silica beads (RPI, Mt. Prospect, IL, USA). Conidia were ruptured using a Mini-Beadbeater (BioSpec Products Inc, Bartlesville, OK, USA) and the organic phase was separated by centrifugation and transferred to new vials. After transfer, the organic phase was evaporated and resuspended in 0.1 mL CHCl3. Samples were spotted onto a thin-layer chromatography (TLC) silica plate (Kiesel gel 60, 0.25 mm; Merck) and resolved in toluene:ethyl acetate:acetic acid (8:1:1 v/v). TLC plates were treated with 1% aluminum hydroxide hydrate (Sigma, St. Louis, MO, USA). Images of TLC plates were captured following UV exposure (366 nm). Quantification of sterigmatocystin and aflatoxin spot intensities was calculated using Image J software. Experiments were performed in triplicate.

2.8. Real-Time PCR Analysis

Real-time PCR analysis was conducted as described previously [28]. The RNA was extracted as described above. For cDNA synthesis, the GoScript Reverse Transcription System (Promega, Madison, WI, USA) was used. Gene-specific primers (Table S1) and iTaq Universal SYBR Green Supermix (Bio-Rad, Hercules, CA, USA) were used for quantitative real-time PCR. To calculate the expression levels of the target genes, the 2−∆∆CT method was used. β-actin was used as an endogenous control. All experiments were carried out in triplicate.

2.9. Statistical Analysis

Statistical differences between the control and mutant strains were evaluated using unpaired Student’s t-tests. Data shown represent the mean ± standard deviation.

3. Results

3.1. VadA Controls Transcript Expression of a Variety of Genes in Conidia

A previous study showed that loss of vadA results in reduced spore viability, reduced trehalose content, reduced tolerance to oxidative stress, and increased β (1,3)-glucan levels and sterigmatocystin production in conidia [26]. To further assess the regulatory role of VadA in conidia, transcriptomic analysis of control conidia and conidia with deletion of the vadA gene (ΔvadA) was performed. As shown Figure S1, deletion of vadA led to differential expression of 3678 genes in A. nidulans (fold change > 2.0; q-value < 0.05). 1489 genes were downregulated, and 2189 genes were upregulated in the conidia of ΔvadA mutants.
To further identify the role of VadA in conidia, gene ontology (GO) analyses were performed [38] using the ASPGD and FungiFun platforms. A large number of genes associated with responses to stimuli, including heat, oxidative, and UV, and small molecule metabolic processes were downregulated in ΔvadA conidia (Figure 1). Whereas the genes involved in the transmembrane transport, secondary metabolism, anatomical structure morphogenesis, and polysaccharide metabolic process were upregulated in the ΔvadA conidia (Figure 1 and Figure S2). These findings suggest that VadA is involved in polysaccharide metabolic processes, secondary metabolic processes, trehalose formation, and stress response.

3.2. VadA Is Required for Expression of Genes Involved in Conidial Wall Integrity

Previous findings have suggested that ΔvadA conidia have reduced levels of trehalose, but increased levels of β-glucan, as compared to wild type conidia. VadA was shown to affect the expression of key genes involved in trehalose (tpsA and orlA) and β-glucan (fksA) biosynthesis [26]. In the present study, RNA-seq data revealed that the mRNA expression of a variety of genes associated with trehalose, β-glucan, and chitin biogenesis were affected by deletion of vadA (Figure 2). The expression of genes related to trehalose synthesis, such as tpsA, orlA, ccg9, and tpsC, were decreased, and the mRNA expression of treA (encoding alpha, alpha-trehalase) was increased in the ΔvadA conidia (Figure 2A). Interestingly, the mRNA levels of many genes associated with chitin and β-glucan biosynthesis or degradation were increased in the mutant conidia (Figure 2B–F). These data suggest that VadA is required for the gene expression of genes involved in cell wall integrity, as the deletion of vadA results in alterations in the content of these carbohydrates.

3.3. VadA Is Required for Conidial Germination

During spore germination and hyphal growth, the dynamic gene expression for the genes involved in the biosynthesis of chitin when β-glucan biosynthesis is occurred [39,40]. As we observed alterations in the gene expression of genes involved in these pathways in the ΔvadA conidia, we hypothesized that vadA deletion may affect conidial germination. To test this hypothesis, the conidial germination rate was assessed. As shown Figure 3, the germination rate of the ΔvadA conidia was faster than that of the control conidia and complemented strains. The ΔvadA conidia, but not the control or complemented conidia, formed germ tubes, suggesting that VadA is required for proper conidial germination.

3.4. Deletion of vadA Affects for UV Stress Response

Previous studies have also shown that ΔvadA conidia are highly susceptible to oxidative stress [26]. In the present study, functional analysis showed that the expression of stress response-related genes was reduced in the ΔvadA conidia (Figure 1). Importantly, the mRNA levels of genes associated with response to radiation were decreased (Figure 4A). Therefore, we assessed the tolerance of the ΔvadA conidia against UV stress. As shown Figure 4B, the ΔvadA conidia were more sensitive to UV stress than the control conidia. Taken together, these results, which are consistent with previous findings, indicate that VadA is important for UV and oxidative stress response in A. nidulans conidia.

3.5. Role of VadA in Secondary Metabolite Gene Clusters in A. nidulans Conidia

Our previous study identified that VadA is required for proper sterigmatocystin production during development [26]. In the present study, we used RNA-seq data to obtain further insight into vadA regulation of secondary metabolism in A. nidulans conidia. Among the 25 secondary metabolite gene clusters, the mRNA expression of genes related to six secondary metabolites gene clusters was upregulated in the ΔvadA conidia (Table 1). In the sterigmatocystin gene cluster, the mRNA levels of a majority of the genes were increased in the mutant conidia (Figure 5A,B). In addition, the amount of sterigmatocystin in the ΔvadA conidia was increased as compared to that of control conidia, suggesting that VadA is negatively associated with sterigmatocystin production in conidia (Figure 5C,D). In addition to sterigmatocystin biosynthesis, other secondary metabolites gene clusters involved in the production of asperfuranone, monodictyphenone, emericellamide, and asperthecin were also upregulated with the deletion of vadA (Figure 6). Overall, these results indicate that VadA is important for proper secondary metabolism in A. nidulans conidia.

4. Discussion

Conidia production is tightly regulated by spore-specific transcription factors. In particular, WetA, VosA, and VelB have been identified as key transcription factors in Aspergillus species [21,24]. Previous studies identified a variety of target genes of the VosA-VelB complex in conidia in A. nidulans [24,25]. VadA is one VosA/VelB-activated developmental gene that is a conidial specific gene in the Aspergillus species [26]. Deletion of vadA affects the balance between asexual and sexual development, conidial viability, conidial trehalose and β-glucan content, and sterigmatocystin production, suggesting that VadA is a multifunctional protein involved in fungal development and conidial maturation in A. nidulans.
To further assess the mode of action of VadA in A. nidulans, RNA-sequencing was performed in asexual spores from wildtype and vadA deletion mutants. Our analyses revealed that the absence of vadA affects the expression of more than 30% of the genes in the A. nidulans genome. The genetic changes caused by the deletion of vadA are consistent with previous phenotypic findings [26]. The major phenotypes observed in the ΔvadA conidia were loss of conidial viability, reduced trehalose content in conidia, and increased conidial sensitivity against UV and oxidative stresses. Our transcriptome analysis indicated that loss of vadA led to decreased mRNA expression of genes associated with spore viability (vosA, velB) [36,41], oxidative stress response (atfA, catA, trxA, msrA, napA), and UV stress response (uvsD, uvsF) [42] indicating that VadA is a key regulator of stress response in tolerance in conidia.
Another feature of ΔvadA conidia is that they have higher sterigmatocystin content than control conidia (Figure 5). Among the 23 genes in the sterigmatocystin gene cluster, the mRNA expression of 21 genes (nine genes were DEG) were increased in the ΔvadA strains. This suggests that upregulation of these genes may affect sterigmatocystin production. In addition, transcripts of genes in gene clusters involved in the production of asperfuranone, asperthecin, emericellamide, and monodictyphenone were also increased in the ΔvadA conidia. Several studies have aimed to identify the gene clusters of these secondary metabolites; however, the function of these metabolites in fungal growth and development has not yet been studied. Therefore, more research is needed to elucidate the functional association of vadA with these pathways.
Chitin and β-glucan are the primary polysaccharides involved in cell wall integrity in most fungi [9,43]. There are dynamic levels of chitin and β-glucan biosynthesis during hyphal tip extension [44,45]. In conidia, the VosA-VelB complex negatively regulates the expression of various genes, including cell wall-related genes, resulting in the repression of chitin and β-glucan biosynthesis [24]. In the ΔvadA conidia, the β-glucan content and the expression of genes involved in β-glucan biosynthesis were increased, suggesting that VadA is also necessary for proper repression of β-glucan biosynthesis in conidia. One possible explanation for the role of VadA in β-glucan biosynthesis in conidia is that VadA indirectly regulates β-glucan biosynthesis through the VosA-VelB complex. Deletion of vadA results in decreased mRNA expression of vosA and velB, two key repressors of β-glucan biosynthesis. However, as the molecular function of VadA has not been studied in detail, further work is necessary to elucidate the molecular roles of VadA.
The new finding in this study is that deletion of vadA led to increased germination rates (Figure 3). Several factors are involved in the gemmation, one of the important components is trehalose [40]. This is because trehalose acts as not only a protectant against oxidative stress in spores, but also an energy source during spore germination [46]. A previous study found that the ΔvadA conidia contain less amount of trehalose than WT conidia [26] The reduced amount of trehalose in the ΔvadA conidia may induce spore germination. Another speculation about the acceleration of the conidial germination rate of the ΔvadA conidia is trehalose hydrolysis. A recent metabolomic study proposes that trehalose hydrolysis plays an important role in conidia dormancy breakage and germination [47]. RNA-seq results showed that mRNA expression of treA [48], was increased in the ΔvadA conidia, speculating that trehalose hydrolysis can be induced, due to the deletion of vadA. Further study is needed to clarify the role of VadA in conidial germination.
Based on the previous phenotypic findings and the genome-wide data obtained in the present study, we propose a model of VadA-related regulation in A. nidulans conidia (Figure 7). In conidia, VadA primarily affects a variety of biological processes. First, VadA is required for the regulation of mRNA expression of secondary metabolites gene clusters, such as sterigmatocystin, asperfuranone, asperthecin, emericellamide, and monodictyphenone. In particular, VadA may play a role in repressing the expression of sterigmatocystin related genes and the production of sterigmatocystin. Second, VadA represses the mRNA expression of genes associated with cell-wall integrity, thereby affecting the production of polysaccharides, such as β-glucan and chitin. Third, VadA involves in the regulation of mRNA expression of genes associated with trehalose biosynthesis, which then affects oxidative stress response and conidial germination. Forth, VadA is involved in the oxidative and UV stress response in conidia. Overall, these findings suggest that VadA is a conidia-specific regulator that governs metabolic processes and conidial maturation in A. nidulans. This study provides insights about how to regulate conidia viability and mycotoxin metabolism in conidia in other pathogenic fungi, such as A. fumigatus and A. flavus.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2607/8/9/1319/s1. Figure S1: Heat map showing the clustering of the differentially expressed genes (DEGs) between control and ΔvadA conidia (fold change > 2.0 and q-value < 0.05) is presented. Figure S2: Functional enrichment analysis in ΔvadA conidia. The charts show categories that were significantly enriched for down-regulated (A) or up-regulated (B) genes. The analysis was performed using the FungiFun database. Table S1: Oligonucleotides used in this study.

Author Contributions

Conceptualization, H.-S.P.; methodology, validation, formal analysis, and data curation, Y.-E.S.; writing—original draft preparation, view and editing, Y.-E.S. and H.-S.P.; supervision, project administration, and funding acquisition, H.-S.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Research Foundation of Korea (NRF) grant to HSP funded by the Korean government (NRF-2016R1C1B2010945 and NRF-2020R1C1C1004473).

Acknowledgments

We thank lab members for critically reading this manuscript.

Conflicts of Interest

The authors have no conflict of interest to declare.

References

  1. Jung, B.; Kim, S.; Lee, J. Microcyle conidiation in filamentous fungi. Mycobiology 2014, 42, 1–5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Park, H.-S.; Yu, J.-H. Genetic control of asexual sporulation in filamentous fungi. Curr. Opin. Microbiol. 2012, 15, 669–677. [Google Scholar] [CrossRef] [PubMed]
  3. Wyatt, T.T.; Wosten, H.A.; Dijksterhuis, J. Fungal spores for dispersion in space and time. Adv. Appl. Microbiol. 2013, 85, 43–91. [Google Scholar] [CrossRef] [PubMed]
  4. Ebbole, D.J. The conidium. In Cellular and Molecular Biology of Filamentous Fungi; American Society for Microbiology Press: Washington, DC, USA, 2010; pp. 577–590. [Google Scholar]
  5. Latge, J.P. Aspergillus fumigatus and aspergillosis. Clin. Microbiol. Rev. 1999, 12, 310–350. [Google Scholar] [CrossRef] [Green Version]
  6. Tracy, M.C.; Okorie, C.U.A.; Foley, E.A.; Moss, R.B. Allergic Bronchopulmonary Aspergillosis. J. Fungi 2016, 2, 17. [Google Scholar] [CrossRef]
  7. Huang, X.H.; Zhang, X.H.; Li, Y.H.; Wang, J.L.; Yan, X.; Xing, L.X.; Wang, F.R. Carcinogenic effects of sterigmatocystin and deoxynivalenol in NIH mice. Zhonghua Zhong Liu Za Zhi 2004, 26, 705–708. [Google Scholar]
  8. Dijksterhuis, J. Fungal spores: Highly variable and stress-resistant vehicles for distribution and spoilage. Food Microbiol. 2019, 81, 2–11. [Google Scholar] [CrossRef]
  9. Beauvais, A.; Latge, J.P. Special Issue: Fungal Cell Wall. J. Fungi 2018, 4, 91. [Google Scholar] [CrossRef] [Green Version]
  10. Bowman, S.M.; Free, S.J. The structure and synthesis of the fungal cell wall. Bioessays 2006, 28, 799–808. [Google Scholar] [CrossRef]
  11. Aimanianda, V.; Simenel, C.; Garnaud, C.; Clavaud, C.; Tada, R.; Barbin, L.; Mouyna, I.; Heddergott, C.; Popolo, L.; Ohya, Y.; et al. The dual activity responsible for the elongation and branching of beta-(1,3)-glucan in the fungal cell wall. mBio 2017, 8. [Google Scholar] [CrossRef] [Green Version]
  12. Gow, N.A.R.; Latge, J.P.; Munro, C.A. The fungal cell wall: Structure, biosynthesis, and function. Microbiol. Spectr. 2017, 5. [Google Scholar] [CrossRef] [Green Version]
  13. Latge, J.P.; Beauvais, A.; Chamilos, G. The cell wall of the human fungal pathogen aspergillus fumigatus: Biosynthesis, organization, immune response, and virulence. Annu. Rev. Microbiol. 2017, 71, 99–116. [Google Scholar] [CrossRef] [PubMed]
  14. Yoshimi, A.; Miyazawa, K.; Abe, K. Cell wall structure and biogenesis in Aspergillus species. Biosci. Biotechnol. Biochem. 2016, 80, 1700–1711. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Garcia-Rubio, R.; de Oliveira, H.C.; Rivera, J.; Trevijano-Contador, N. The Fungal Cell Wall: Candida, Cryptococcus, and Aspergillus Species. Front. Microbiol. 2019, 10, 2993. [Google Scholar] [CrossRef]
  16. De Vries, R.P.; Riley, R.; Wiebenga, A.; Aguilar-Osorio, G.; Amillis, S.; Uchima, C.A.; Anderluh, G.; Asadollahi, M.; Askin, M.; Barry, K.; et al. Comparative genomics reveals high biological diversity and specific adaptations in the industrially and medically important fungal genus Aspergillus. Genome Biol. 2017, 18, 28. [Google Scholar] [CrossRef] [Green Version]
  17. Etxebeste, O.; Espeso, E.A. Aspergillus nidulans in the post-genomic era: A top-model filamentous fungus for the study of signaling and homeostasis mechanisms. Int. Microbiol. 2019. [Google Scholar] [CrossRef]
  18. Yu, J.H. Regulation of development in Aspergillus nidulans and Aspergillus fumigatus. Mycobiology 2010, 38, 229–237. [Google Scholar] [CrossRef] [Green Version]
  19. Ojeda-Lopez, M.; Chen, W.; Eagle, C.E.; Gutierrez, G.; Jia, W.L.; Swilaiman, S.S.; Huang, Z.; Park, H.S.; Yu, J.H.; Canovas, D.; et al. Evolution of asexual and sexual reproduction in the Aspergilli. Stud. Mycol. 2018, 91, 37–59. [Google Scholar] [CrossRef]
  20. Wu, M.Y.; Mead, M.E.; Kim, S.C.; Rokas, A.; Yu, J.H. WetA bridges cellular and chemical development in Aspergillus flavus. PLoS ONE 2017, 12, e0179571. [Google Scholar] [CrossRef] [Green Version]
  21. Wu, M.Y.; Mead, M.E.; Lee, M.K.; Ostrem Loss, E.M.; Kim, S.C.; Rokas, A.; Yu, J.H. Systematic dissection of the evolutionarily conserved weta developmental regulator across a genus of filamentous fungi. mBio 2018, 9. [Google Scholar] [CrossRef] [Green Version]
  22. Park, H.S.; Yu, J.H. Velvet regulators in Aspergillus spp. Microbiol. Biotechnol. Lett. 2017, 44, 409–419. [Google Scholar] [CrossRef]
  23. Bayram, O.; Braus, G.H. Coordination of secondary metabolism and development in fungi: The velvet family of regulatory proteins. FEMS Microbiol. Rev. 2012, 36, 1–24. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Park, H.-S.; Yu, Y.M.; Lee, M.-K.; Maeng, P.J.; Kim, S.C.; Yu, J.-H. Velvet-mediated repression of beta-glucan synthesis in Aspergillus nidulans spores. Sci. Rep. 2015, 5, 10199. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Ahmed, Y.L.; Gerke, J.; Park, H.S.; Bayram, O.; Neumann, P.; Ni, M.; Dickmanns, A.; Kim, S.C.; Yu, J.H.; Braus, G.H.; et al. The velvet family of fungal regulators contains a DNA-binding domain structurally similar to NF-kappaB. PLoS Biol. 2013, 11, e1001750. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Park, H.S.; Lee, M.K.; Kim, S.C.; Yu, J.H. The role of VosA/VelB-activated developmental gene vadA in Aspergillus nidulans. PLoS ONE 2017, 12, e0177099. [Google Scholar] [CrossRef] [Green Version]
  27. Park, H.-S.; Yu, J.-H. Multi-copy genetic screen in Aspergillus nidulans. Methods Mol. Biol. 2012, 944, 183–190. [Google Scholar] [CrossRef]
  28. Kim, M.J.; Lee, M.K.; Pham, H.Q.; Gu, M.J.; Zhu, B.; Son, S.H.; Hahn, D.; Shin, J.H.; Yu, J.H.; Park, H.S.; et al. The velvet regulator VosA governs survival and secondary metabolism of sexual spores in Aspergillus nidulans. Genes 2020, 11, 103. [Google Scholar] [CrossRef] [Green Version]
  29. Arnaud, M.B.; Chibucos, M.C.; Costanzo, M.C.; Crabtree, J.; Inglis, D.O.; Lotia, A.; Orvis, J.; Shah, P.; Skrzypek, M.S.; Binkley, G.; et al. The Aspergillus genome database. A curated comparative genomics resource for gene, protein and sequence information for the Aspergillus research community. Nucl. Acids Res. 2010, 38, D420–D427. [Google Scholar] [CrossRef] [Green Version]
  30. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  31. Trapnell, C.; Williams, B.A.; Pertea, G.; Mortazavi, A.; Kwan, G.; van Baren, M.J.; Salzberg, S.L.; Wold, B.J.; Pachter, L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 2010, 28, 511–515. [Google Scholar] [CrossRef] [Green Version]
  32. Anders, S.; Pyl, P.T.; Huber, W. HTSeq—A Python framework to work with high-throughput sequencing data. Bioinformatics 2015, 31, 166–169. [Google Scholar] [CrossRef] [PubMed]
  33. Sun, J.; Nishiyama, T.; Shimizu, K.; Kadota, K. TCC: An R package for comparing tag count data with robust normalization strategies. BMC Bioinf. 2013, 14, 219. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Priebe, S.; Kreisel, C.; Horn, F.; Guthke, R.; Linde, J. FungiFun2: A comprehensive online resource for systematic analysis of gene lists from fungal species. Bioinformatics 2015, 31, 445–446. [Google Scholar] [CrossRef] [Green Version]
  35. Sarikaya Bayram, O.; Bayram, O.; Valerius, O.; Park, H.-S.; Irniger, S.; Gerke, J.; Ni, M.; Han, K.H.; Yu, J.-H.; Braus, G.H. LaeA control of velvet family regulatory proteins for light-dependent development and fungal cell-type specificity. PLoS Genet. 2010, 6, e1001226. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Park, H.S.; Ni, M.; Jeong, K.C.; Kim, Y.H.; Yu, J.H. The role, interaction and regulation of the velvet regulator VelB in Aspergillus nidulans. PLoS ONE 2012, 7, e45935. [Google Scholar] [CrossRef] [Green Version]
  37. Son, Y.E.; Cho, H.J.; Chen, W.; Son, S.H.; Lee, M.K.; Yu, J.H.; Park, H.S. The role of the VosA-repressed dnjA gene in development and metabolism in Aspergillus species. Curr. Genet. 2020, 66, 621–633. [Google Scholar] [CrossRef] [PubMed]
  38. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T.; et al. Gene ontology: Tool for the unification of biology. The gene ontology consortium. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef] [Green Version]
  39. Van Leeuwen, M.R.; Krijgsheld, P.; Bleichrodt, R.; Menke, H.; Stam, H.; Stark, J.; Wosten, H.A.; Dijksterhuis, J. Germination of conidia of Aspergillus niger is accompanied by major changes in RNA profiles. Stud. Mycol. 2013, 74, 59–70. [Google Scholar] [CrossRef]
  40. Baltussen, T.J.H.; Zoll, J.; Verweij, P.E.; Melchers, W.J.G. Molecular mechanisms of conidial germination in Aspergillus spp. Microbiol. Mol. Biol. Rev. 2020, 84. [Google Scholar] [CrossRef]
  41. Ni, M.; Yu, J.H. A novel regulator couples sporogenesis and trehalose biogenesis in Aspergillus nidulans. PLoS ONE 2007, 2, e970. [Google Scholar] [CrossRef] [Green Version]
  42. Kafer, E.; Mayor, O. Genetic analysis of DNA repair in Aspergillus: Evidence for different types of MMS-sensitive hyperrec mutants. Mutat. Res. 1986, 161, 119–134. [Google Scholar] [CrossRef]
  43. Free, S.J. Fungal cell wall organization and biosynthesis. Adv. Genet. 2013, 81, 33–82. [Google Scholar] [CrossRef] [PubMed]
  44. Lenardon, M.D.; Munro, C.A.; Gow, N.A. Chitin synthesis and fungal pathogenesis. Curr. Opin. Microbiol. 2010, 13, 416–423. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Patel, P.K.; Free, S.J. The genetics and biochemistry of cell wall structure and synthesis in Neurospora crassa, a model filamentous fungus. Front. Microbiol. 2019, 10, 2294. [Google Scholar] [CrossRef] [PubMed]
  46. Fillinger, S.; Chaveroche, M.K.; van Dijck, P.; de Vries, R.; Ruijter, G.; Thevelein, J.; d’Enfert, C. Trehalose is required for the acquisition of tolerance to a variety of stresses in the filamentous fungus Aspergillus nidulans. Microbiology 2001, 147, 1851–1862. [Google Scholar] [CrossRef] [Green Version]
  47. Kalampokis, I.F.; Erban, A.; Amillis, S.; Diallinas, G.; Kopka, J.; Aliferis, K.A. Untargeted metabolomics as a hypothesis-generation tool in plant protection product discovery: Highlighting the potential of trehalose and glycerol metabolism of fungal conidiospores as novel targets. Metabolomics 2020, 16, 79. [Google Scholar] [CrossRef]
  48. D’Enfert, C.; Fontaine, T. Molecular characterization of the Aspergillus nidulans treA gene encoding an acid trehalase required for growth on trehalose. Mol. Microbiol. 1997, 24, 203–216. [Google Scholar] [CrossRef]
Figure 1. Transcriptomic analyses of ΔvadA conidia. GO analysis of DEGs in ΔvadA conidia.
Figure 1. Transcriptomic analyses of ΔvadA conidia. GO analysis of DEGs in ΔvadA conidia.
Microorganisms 08 01319 g001
Figure 2. mRNA expression of cell-wall integrity-related genes in ΔvadA conidia. (AE) Heat map showing the differentially expressed genes (DEGs) of trehalose biogenesis (A), chitin degradation (B), chitin biosynthesis (C), β-glucan degradation and β-glucan biosynthesis (C) in control and ΔvadA conidia. (F) mRNA expression of gelA and gelB in control (THS30), ΔvadA (THS33), and C’vadA (THS34) conidia were assessed by qRT-PCR. Error bars indicates the standard errors from the mean of triplicates.
Figure 2. mRNA expression of cell-wall integrity-related genes in ΔvadA conidia. (AE) Heat map showing the differentially expressed genes (DEGs) of trehalose biogenesis (A), chitin degradation (B), chitin biosynthesis (C), β-glucan degradation and β-glucan biosynthesis (C) in control and ΔvadA conidia. (F) mRNA expression of gelA and gelB in control (THS30), ΔvadA (THS33), and C’vadA (THS34) conidia were assessed by qRT-PCR. Error bars indicates the standard errors from the mean of triplicates.
Microorganisms 08 01319 g002
Figure 3. Role of vadA in germination. (A) Conidial germination rate in control (THS30), ΔvadA (THS33), and C’ vadA (THS34) conidia inoculated on solid minimal medium with glucose. (B) Control (THS30), ΔvadA (THS33), and C’ vadA (THS34) conidia were inoculated on solid MM (1% glucose), with or without a glucose source and incubated for 4 or 5 h, respectively (bar = 25 µm). Arrows indicate germ tubes.
Figure 3. Role of vadA in germination. (A) Conidial germination rate in control (THS30), ΔvadA (THS33), and C’ vadA (THS34) conidia inoculated on solid minimal medium with glucose. (B) Control (THS30), ΔvadA (THS33), and C’ vadA (THS34) conidia were inoculated on solid MM (1% glucose), with or without a glucose source and incubated for 4 or 5 h, respectively (bar = 25 µm). Arrows indicate germ tubes.
Microorganisms 08 01319 g003
Figure 4. Role of vadA in the UV-mediated stress response. (A) Heat map showing the DEGs with UV stress response. (B) The tolerance of control (THS30), ΔvadA (THS33), and C’ vadA (THS34) conidia against UV irradiation. About 102 conidia were irradiated with 50 or 75 J/m2 of UV at 260 nm (triplicate measurements). Error bars represent the standard deviation (differences between the vadA mutants and control strains). ** p < 0.01; * p < 0.05.
Figure 4. Role of vadA in the UV-mediated stress response. (A) Heat map showing the DEGs with UV stress response. (B) The tolerance of control (THS30), ΔvadA (THS33), and C’ vadA (THS34) conidia against UV irradiation. About 102 conidia were irradiated with 50 or 75 J/m2 of UV at 260 nm (triplicate measurements). Error bars represent the standard deviation (differences between the vadA mutants and control strains). ** p < 0.01; * p < 0.05.
Microorganisms 08 01319 g004
Figure 5. The role of VadA in sterigmatocystin production in conidia. (A) Heat map diagram of genes in sterigmatocystin gene clusters. (B) mRNA expression of stcE and stcU in control (THS30), ΔvadA (THS33), and C’vadA (THS34) conidia were assessed by qRT-PCR. (C) TLC image of sterigmatocystin from control (THS30), ΔvadA (THS33), and C’ vadA (THS34) conidia. (D) Relative production of sterigmatocystin shown in (C). Error bars represent standard deviation (differences between the vadA mutants and control strains. *** p < 0.001).
Figure 5. The role of VadA in sterigmatocystin production in conidia. (A) Heat map diagram of genes in sterigmatocystin gene clusters. (B) mRNA expression of stcE and stcU in control (THS30), ΔvadA (THS33), and C’vadA (THS34) conidia were assessed by qRT-PCR. (C) TLC image of sterigmatocystin from control (THS30), ΔvadA (THS33), and C’ vadA (THS34) conidia. (D) Relative production of sterigmatocystin shown in (C). Error bars represent standard deviation (differences between the vadA mutants and control strains. *** p < 0.001).
Microorganisms 08 01319 g005
Figure 6. Expression of secondary metabolite gene clusters in ΔvadA conidia. (AD) Heat map diagram of genes related to emericellamide (A), monodictyphenone (B), asperthecin (C), and asperfuranone (D) gene clusters. (E) mRNA expression of afoA and afoB in control (THS30), ΔvadA (THS33), and C’vadA (THS34) conidia were assessed by qRT-PCR.
Figure 6. Expression of secondary metabolite gene clusters in ΔvadA conidia. (AD) Heat map diagram of genes related to emericellamide (A), monodictyphenone (B), asperthecin (C), and asperfuranone (D) gene clusters. (E) mRNA expression of afoA and afoB in control (THS30), ΔvadA (THS33), and C’vadA (THS34) conidia were assessed by qRT-PCR.
Microorganisms 08 01319 g006
Figure 7. Genetic model for VadA-mediated regulation in A. nidulans conidia. A proposed model for the VadA-mediated conidiogenesis is presented.
Figure 7. Genetic model for VadA-mediated regulation in A. nidulans conidia. A proposed model for the VadA-mediated conidiogenesis is presented.
Microorganisms 08 01319 g007
Table 1. Upregulated genes in secondary metabolite genes clusters in the ΔvadA conidia.
Table 1. Upregulated genes in secondary metabolite genes clusters in the ΔvadA conidia.
Secondary Metabolite Gene Cluster
(Number of Genes in the Cluster)
Upregulated Genes in the ΔvadA Conidia
Asperfuranone (7)afoA, afoB, afoC, afoD, afoE, afoF, afoG
Asperthecin (3)aptA, aptB, aptC
Aspyridone (8)apdB, apdA, apdF, apdG
Austinol cluster 1(10)ausH, ausI, ausJ, ausK
Austinol cluster 2 (4)ausC, ausA, ausD
Derivative of Benzaldehyde1 and F9775 hybrid cluster 1 (9)dbaA, dbaB, dbaC, dbaG, dbaH, dbaI
Derivative of Benzaldehyde1 and F9775 hybrid cluster 2 (3)orsA, orsB, orsC
Emericellamide cluster (4)easA, easB, easC, easD
inp cluster (6)inpC, inpD, scpR, inpE
ivo cluster (2)_
Microperfuranone cluster (3)AN3394, AN3395, micA
Monodictyphenone cluster (12)mdpA, mdpB, mdpC, mdpD, mdpE, mdpF, mdpG, mdpH, mdpI, mdpJ, mdpK, mdpL
Penicillin cluster (3)acvA, ipnA, aatA
Nidulanin A cluster (2)AN11080
Sterigmatocystin cluster (23)stcW, stcV, stcU, stcS, AN7809, stcQ, stcO, stcN, stcL, AN7815, stcI, AN7817, stcF, AN7819, aflR, stcE, AN7822, stcC, stcA, stcK, AN11021
Terriquinone cluster (4)tdiA, tdiC
pkb cluster (8)cicC, cicE, pkbA, AN6449, AN6450, AN6451
pkdA cluster (11)pkdA, AN0524, AN0525, AN0526, AN0527, AN0528, AN0529, AN0530, AN0531, AN0533
pkf cluster (6)AN3225, pkfA, pkfC, pkfD, pkfE, pkfF
pkg cluster (8)AN7072, AN7074, AN7075, AN10884
pkh cluster (8)AN2031, AN2033, AN2034, pkhB, AN2036, AN2037, AN2038
pki cluster (8)AN3379, pkiB, pkiC, salA, AN3383, AN3385, AN3386
xptA-containing cluster (10)xptA, AN6785, AN6786, AN6789, AN6791
xptB-containing cluster (4)AN7999, xptB, AN9467
sidC cluster (4)AN0606, sidC, AN0608

Share and Cite

MDPI and ACS Style

Son, Y.-E.; Park, H.-S. Genome Wide Analysis Reveals the Role of VadA in Stress Response, Germination, and Sterigmatocystin Production in Aspergillus nidulans Conidia. Microorganisms 2020, 8, 1319. https://doi.org/10.3390/microorganisms8091319

AMA Style

Son Y-E, Park H-S. Genome Wide Analysis Reveals the Role of VadA in Stress Response, Germination, and Sterigmatocystin Production in Aspergillus nidulans Conidia. Microorganisms. 2020; 8(9):1319. https://doi.org/10.3390/microorganisms8091319

Chicago/Turabian Style

Son, Ye-Eun, and Hee-Soo Park. 2020. "Genome Wide Analysis Reveals the Role of VadA in Stress Response, Germination, and Sterigmatocystin Production in Aspergillus nidulans Conidia" Microorganisms 8, no. 9: 1319. https://doi.org/10.3390/microorganisms8091319

APA Style

Son, Y. -E., & Park, H. -S. (2020). Genome Wide Analysis Reveals the Role of VadA in Stress Response, Germination, and Sterigmatocystin Production in Aspergillus nidulans Conidia. Microorganisms, 8(9), 1319. https://doi.org/10.3390/microorganisms8091319

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop