Next Article in Journal
The Finnish New Variant of Chlamydia trachomatis with a Single Nucleotide Polymorphism in the 23S rRNA Target Escapes Detection by the Aptima Combo 2 Test
Next Article in Special Issue
Encapsulated Bdellovibrio Powder as a Potential Bio-Disinfectant against Whiteleg Shrimp-Pathogenic Vibrios
Previous Article in Journal
Wild Nutria (Myocastor coypus) Is a Potential Reservoir of Carbapenem-Resistant and Zoonotic Aeromonas spp. in Korea
Previous Article in Special Issue
Trends in Diatom Research Since 1991 Based on Topic Modeling
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 7
Issue 8
10.3390/microorganisms7080225
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Oxidative Stress Response of Aspergillus oryzae Induced by Hydrogen Peroxide and Menadione Sodium Bisulfite

by
Huanhuan Shao
1,2,
Yayi Tu
1,
Yijing Wang
1,
Chunmiao Jiang
1,
Long Ma
1,
Zhihong Hu
1,
Jiangfan Wang
1,
Bin Zeng
1,* and
Bin He
1,*
1
Jiangxi Key Laboratory of Bioprocess Engineering and Co-Innovation Center for In-vitro Diagnostic Reagents and Devices of Jiangxi Province, College of Life Sciences, Jiangxi Science & Technology Normal University, Nanchang 330013, China
2
College of Life Sciences, Sichuan Normal University, Chengdu 610101, China
*
Authors to whom correspondence should be addressed.
Microorganisms 2019, 7(8), 225; https://doi.org/10.3390/microorganisms7080225
Submission received: 17 June 2019 / Revised: 18 July 2019 / Accepted: 29 July 2019 / Published: 30 July 2019
(This article belongs to the Special Issue Recent Advances in Applied Microbiology)
Browse Figures
Versions Notes

Abstract

:
Oxidative stress response protects organisms from deleterious effects of reactive oxygen species (ROS), which can damage cellular components and cause disturbance of the cellular homeostasis. Although the defensive biochemical mechanisms have been extensively studied in yeast and other filamentous fungi, little information is available about Aspergillus oryzae. We investigated the effect of two oxidant agents (menadione sodium bisulfite, MSB, and hydrogen peroxide, H2O2) on cellular growth and antioxidant enzyme induction in A. oryzae. Results indicated severe inhibition of biomass and conidia production when high concentration of oxidants was used. Transcriptomic analysis showed an up-regulated expression of genes involved in oxidoreduction, such as catalase, glutathione peroxidase, and superoxide dismutase. In addition, it was observed that oxidative stress stimuli enhanced the expression of Yap1 and Skn7 transcription factors. Further, metabolomic analysis showed that glutathione content was increased in the oxidative treatments when compared with the control. Moreover, the content of unsaturated fatty acid decreased with oxidative treatment accompanying with the down-regulated expression of genes involved in linoleic acid biosynthesis. This study provided a global transcriptome characterization of oxidative stress response in A. oryzae, and can offer multiple target genes for oxidative tolerance improvement via genetic engineering.

1. Introduction

Aspergillus oryzae (A. oryzae) has been used for the production of traditional fermentation such as soy sauce, miso, and douchi, and demonstrated promising potential for primary and secondary metabolites [1,2]. In additional, A. oryzae has been listed as “Generally Recognized as Safe (GRAS)” by the Food and Drug Administration (FDA), and its safety has been recognized by the World Health Organization (WHO) [3]. Like other microorganisms, A. oryzae has been subjected to pressure from internal and external environment during the process of growth and fermentation. Under normal aerobic condition, the metabolic byproduct of A. oryzae may result in the formation of reactive oxygen species (ROS), mainly including superoxide anion, hydroxyl radical, and hydrogen peroxide. The accumulation of ROS can damage the cellular components due to their stronger reactivity than molecular oxygen [4]. As a consequence, ROS can cause disturbance of cellular homeostasis, which may lead to the inactivation of enzymes and disruption of unsaturated fatty acid double bond, even resulting in cell ageing or death [5]. Microorganisms have evolved several defense mechanisms against the deleterious effects of ROS. The most detailed studies on the response to oxidative stress have been reported in Saccharomyces cerevisiae [6]. However, the biochemical defense mechanisms of A. oryzae against oxidative stress have not been completely elucidated, which is essential for the survival and adaptation to an oxidative environment in fermentation and production process.
In S. cerevisiae, there are two intricate antioxidant defense systems: enzyme defense and non-enzyme defense system [6]. The antioxidant enzymes, mainly including catalase (CAT), glutathione peroxidase (GPX) and superoxide dismutase (SOD), played an important role in the defense against ROS by catalytically reduce it through electron transfer. There are two SOD isoenzymes in S. cerevisiae (SOD1, SOD2), which catalyze the two-step disproportionation of superoxide anion to hydrogen peroxide and oxygen molecule. Deletion of SOD genes caused an increased sensitivity to oxidative stress-inducing agents, or decreased the level of reduced glutathione (GSH) [7]. Moreover, the Δsod mutant exhibited oxidative stress-related variations in cells, low growth rate in air, and inability to grow in an atmosphere of 100% oxygen in a rich cultivation medium [8]. Furthermore, the loss of SOD dramatically reduced the chronological and replicative lifespans of the yeast [9]. CAT is an iron-containing enzyme, which catalyzes the disproportionation of H2O2 into O2 and H2O molecules. In yeast, there are two CAT isoenzymes, including CTT1 located in the cytosol and CTA1 located in peroxisomes and mitochondria [10]. However, the loss of either individual catalase or both CTT1 and CTA1 exhibited no effect on H2O2 sensitivity [11]. The results indicate that there are other peroxidases, such as GPX, that can effectively remove H2O2 in yeast cell. GPX can catalyze the reduction of H2O2 or phospholipid hydroperoxides to H2O or corresponding alcohol using reduced glutathione [12]. The biological purpose of oxidative stress response in the first place is to reinstate the cellular redox homeostasis. This requires the regulation of transcriptional changes via Yap1 and Skn7 transcription factors responsible for the expression of key genes coding enzymes that are required for ROS detoxification [13]. Yap1 is a major regulator of oxidative stress response in S. cerevisiae, which is distributed mainly in the cell cytosol [14]. Skn7 is a constitutive nuclear protein that usually functions in combination with Yap1 in the response to peroxides. Approximately half of the genes controlled by Yap1 are regulated by both Yap1 and Skn7 (including SOD1 and SOD2), while the other half is regulated in a Skn7-independent manner [15,16].
The non-enzyme defense system reinstates the cellular redox homeostasis through the modifications of cellular components, such as glutathione and lipid. Glutathione plays an important role in the response to oxidative stress that might be the target of oxidation against the exposure to an exogenous oxidant. Glutathione can protect the mitochondria from oxidant damage by reacting with oxidative agent [17,18]. When the level of ROS exceeds the antioxidant capacity of cells, lipids can react with ROS to reinstate the cellular redox homeostasis as well. Under oxidative stress, oxidants attack lipids containing C-C bond(s), particularly polyunsaturated fatty acids [19,20]. Since yeast does not contain polyunsaturated fatty acids, there is no lipid peroxidation process in its cells. Thus, it is imperative to study the biochemical defense mechanisms against oxidative stress in filamentous fungi, such as A. oryzae. Few filamentous fungi, such as Neurospora crassa, Aspergillus niger, Aspergillus fumigatus, and Aspergillus nidulans, have been investigated the defense mechanisms against oxidative stress [21,22,23,24]. Although it is an important industrial microorganism for the production of traditional fermentation and enzymes, there is a scarcity of information about the responses of A. oryzae against various oxidative stresses. The objective of this study was to understand the response of A. oryzae to oxidative stress. Thus, two types of ROS generating agents were examined for their effects on cellular growth and antioxidant enzyme induction in A. oryzae.

2. Materials and Methods

2.1. Effect of Oxidative Stress on Aspergillus oryzae Growth

Aspergillus oryzae (A. oryzae) 3.042 (China center of industrial culture collection (CICC) no. 40092), identified with morphological characteristics and ITS gene sequence, was selected to examine the response to oxidative stress. For inoculum preparation, cryopreserved conidia were activated in fresh potato dextrose agar (PDA) medium and incubated at 30 °C for 3 days. Conidia were then harvested and suspended in sterile water for quantitative determination of conidia concentration using a hemocytometer [25]. Oxidative stress treatments were performed on solid PDA medium supplemented with stress-initiating agents hydrogen peroxide (H2O2 at final concentration of 3 mM and 6 mM) and menadione sodium bisulfite (MSB at final concentration of 0.28 mM and 0.42 mM) [26]. PDA medium absent of oxidant was used as control. Freshly prepared suspension (1 × 107 conidia) was inoculated per plate covered by cellophane, and then all cultures were placed at 30 °C for 72 h in darkness. Then, fungal mycelia were dried overnight at 65 °C for the biomass determination and the density of spores was determined by hemocytometer. For the determination of reduced glutathione and fatty acids, samples were collected in the same way. Three replicates were used for each sample in above experiments and images of the most representative colonies were taken.

2.2. Determination of Reduced Glutathione Levels

For the determination of reduced glutathione (GSH), a previously described method (5,5’-dithiobis-(2-nitrobenzoic acid) (DTNB) method) was used [27]. In brief, cells were harvested and washed twice with 1 mM phosphate buffered saline (pH 7.4) to remove any traces of growth medium. Cells were dried overnight at 65 °C. Dry mycelia (100 mg) were mixed with GSH assay buffer and grind with liquid nitrogen. The grinded cells were centrifuged at 8000 g for 10 min at 4 °C and the suspension obtained was used to determine GSH levels. GSH quantification was performed by incubating the reaction mixture containing crude extract and 3 mM of DTNB at 30 °C for 5 min followed by NADPH addition (0.4 mM) and 2 μL of glutathione reductase (GR) enzyme. Samples were kept for another 2 min to allow the reaction completed and absorbance was measured at 412 nm. GSH levels in the samples were compared with the standard curve prepared by using various concentrations (50, 100, 150, 200 μg/mL) of GSH.

2.3. Measurement of Intracellular Fatty Acid Content

Total lipids were extracted and quantified following a previously described method [25]. Cells were harvested and incubated in chloroform with 2% H2SO4-MeOH solution at 70 °C for 2 h to obtain fatty acid methyl esters (FAMEs). FAMEs components were separated and analyzed using a coupled QP2010 gas chromatography–mass spectrometry (GC-MS) (Shimadzu, Kyoto, Japan). The system was equipped with Supelco SP-2340 fused silica capillary column (30 m × 0.25 mm i.d., film thickness of 0.2 μm; Bellefonte, PA, USA). FAMEs were identified by comparing their mass spectra with a spectrum database. Fatty acid peaks were identified by comparing retention times with external standards or similarity search. Fatty acid quantification was performed by using the peak area of the most intensive ion of each peak [25,28].

2.4. RNA Extraction and Illumina Sequencing

Total RNA was extracted from all samples using Trizol reagent (Roche, Basel, Switzerland) following the manufacturer’s instructions and genomic DNAs were digested by DNase I (Fermentas, Waltham, MA, USA). RNA concentration was analyzed by Qubit fluorometer (Invitrogen, Carlsbad, CA, USA) and its quality was determined using Agilent 2100 Bioanalyzer. Further, mRNA was enriched from the total RNA using Oligo (dT) magnetic beads and fragmented into about 200 nt using fragmentation buffer. These fragments were then reverse transcripted into cDNA with random primers and the second-strand cDNA was synthesized using DNA polymerase I and RNase H. The resulting double-stranded cDNA was purified, end repaired and ligated to Illumina sequencing adapters. The cDNA sequencing was performed using Illumina HiSeq 2500 (Biomarker Biotechnology Co., Beijing, China).

2.5. Differentially Expressed Genes Analysis

After the removal of low-quality reads, adaptor sequences and reads with N (unknown nucleotides) exceeding 10%, clean-read datasets were obtained and aligned to A. oryzae 3.042 reference genome using Tophats2 (v2.0.3.12) [29]. Gene abundances were quantified by RSEM software and the quantification of gene expression level was normalized using FPKM (Fragments Per Kilobase of transcript per Million mapped reads) method [30]. Differentially expressed genes (DEGs) in all samples were identified using DESeq R package [31]. The log2 (fold change) over 1 and false discovery rate (FDR) within 0.05 were set as the threshold for significant DEGs [25]. DEGs were used for the enrichment analysis of GO functions and KEGG pathways.

2.6. Calculations and Statistical Analyses

Three independent experiments were performed and the values shown in this study were presented as mean ± standard error (SE). The effects of treatments were compared by one-way nested analysis of variance (ANOVA), followed by least significant difference test (LSD) for mean comparison. All statistical analysis was performed using SAS 9.20 software (SAS Institute Inc., Cary, NC, USA) at p < 0.05.

3. Results

3.1. Effect of Oxidant Species and Concentration on A. oryzae Cell Growth

To study the effect of oxidant species and its concentration on A. oryzae growth and development, the density of spores, dry biomass, and phenotype analysis were determined and performed. When A. oryzae cells were supplemented with 3 mM H2O2 and 0.28 mM MSB, no significant effects were observed on the spore density between control and oxidative treatments, while the dry biomass was decreased when compared to control (Figure 1A,B). Under high H2O2 (6 mM) and MSB (0.42 mM) concentration, the density of spores and dry biomass of A. oryzae showed a significant difference between the oxidative treatments and control. During the cultivation period, the cell growth and conidia formation was notably inhibited with high concentration of oxidant (6 mM H2O2 and 0.42 mM MSB) compared to others, indicating severe inhibition of biomass and conidia production (Figure 1C).

3.2. Global Analysis of A. oryzae Transcriptome under Oxidative Stress

To investigate the oxidative stress response of A. oryzae induced by H2O2 and MSB, a transcriptome analysis based on RNA sequencing of A. oryzae that was subjected to various oxidative stress treatments (no oxidant, 3 mM and 6 mM H2O2, 0.28 mM and 0.42 mM MSB) was performed. This resulted into the generation of 40.79, 41.41, 41.10, 43.61 and 39.90 million clean reads per library, respectively. The Q30 values (99.9% accuracy of bases) were higher than 94% for all samples, indicating the good quality of sequencing data. The clean reads were aligned to A. oryzae 3042 genome sequence, and more than 84% of the clean reads for each sample were uniquely mapped to the genome (Table 1).

3.3. Differentially Expressed genes Analysis of A. oryzae Transcriptome under Oxidative Stress

The relative expression level of each gene in terms of FPKM was estimated by mapping all clean reads from each library back to the reference genome. Compared to the control, a total of 880 genes were identified as DEGs at 3 mM H2O2, which comprised 410 up-regulated genes (accounting for 47% of all significant DEGs) and 470 down-regulated genes (accounting for 53%). Between control and 6 mM H2O2 treatment, 1030 DEGs were identified, including 456 up-regulated genes (44%) and 574 down-regulated genes (56%). Under 0.28 and 0.42 mM MSB treatments, 882 and 649 DEGs were identified (Figure 2A). A Venn diagram of DEGs distribution was constructed (Figure 2B), and it was observed that only 56 DEGs were commonly shared among the four groups. The numbers of specific DEGs between 6 mM H2O2 treatments and control (408) as well as 0.42 mM MSB treatments and control (355) were remarkably greater than that of the other two groups (Figure 2B). This suggested the involvement of complex developmental events at high concentration of oxidant treatments (Figure 2B). In addition, we found that more common DEGs were found between WT-vs-H2O2_3 and WT-vs-H2O2_6 or WT-vs-MSB_28 and WT-vs-MSB_42 when compared to that between WT-vs-H2O2 and WT-vs-MSB, suggesting that the response of A. oryzae to H2O2 and MSB might be different.

3.4. Expression Analysis of Genes Involved in Oxidoreduction under Oxidative Stress

Under oxidative stress, fungus could eliminate the effects of reactive oxygen species (ROS) reactions through the expression genes involved in oxidoreduction, including CAT, GPX, and SOD [32]. Therefore, the genes encoding CAT, GPX and SOD were identified in A. oryzae transcriptome and their expression profile was analyzed. As shown in Figure 3, six CAT genes, including gene8009, gene3678, gene2938, gene9633, gene7236, and gene1032 were identified in A. oryzae. The results of expression analysis showed that only two CAT genes (gene8009 and gene3678) were up-regulated in MSB treatments, while other four CAT genes (gene2938, gene9633, gene7236, and gene1032) were significantly up-regulated in H2O2 treatments. Furthermore, only one GPX gene was identified in A. oryzae and notably up-regulated at MSB treatments. The four SOD genes showed a different expression profile under H2O2 and MSB treatments. Two SOD genes (gene2883 and gene4586) exhibited a higher expression level in samples cultivated in solid PDA medium supplemented with H2O2, while the remaining two SOD genes (gene9213 and gene5567) exhibited highest expression level under MSB treatments. The different expression patterns of 11 genes involved in oxidoreduction under H2O2 and MSB treatments suggested that the mechanism of A. oryzae in response to various oxidants is complex and diverse.

3.5. Expression Analysis of Transcription Factors Possibly Involved in Oxidative Stress Response

Previous studies have reported that several transcription factors, such as Yap1 and Skn7, responded to oxidative stress through maintaining thiol redox homeostasis [15]. Thus, we identified and analyzed the expression profile of Yap1 and Skn7 transcription factors. A total of three Yap1 transcription factors (gene278, gene10610, and gene6095) were identified in A. oryzae, and they were highly expressed in the samples cultivated in solid PDA medium supplemented with H2O2 (Figure 4). Similarly, 12 Skn7 transcription factors were identified in A. oryzae transcriptome. Seven Skn7 transcription factors exhibited significant increase in the expression levels when treated with H2O2, while the remaining five Skn7 transcription factors showed up-regulated expressions in MSB treatments. These results suggested that the response of A. oryzae to H2O2 and MSB treatment might be different.

3.6. Effects of Oxidant on A. oryzae GSH

GSH plays an important role in the response to oxidative stress as the most important antioxidant through the enzymes involved in oxidoreduction and Yap1 and Skn7 transcription factors. We determined the content of A. oryzae GSH in the oxidative treatments through DTNB method. As we expected, the content of GSH was found to be increased in oxidative treatments compared to control (Figure 5). In addition, GSH content under MSB stress was higher than those under H2O2 stress. To further confirm the increase of glutathione level, the expression profile of genes involved in glutathione biosynthesis, including glutamate-cysteine ligase (Gcl) and glutathione synthase (Gss), was investigated. The results revealed that two genes encoding Gss and the gene encoding Gcl both exhibited highest expression levels in 0.42 mM MSB treatments, which were supported the increase of GSH content.

3.7. Change in Linoleic acid Biosynthesis under Oxidative Stress

According to KEGG pathway analysis, four genes involved in A. oryzae linoleic acid biosynthesis were identified. They encoded delta 9 fatty acid desaturases (D9D) and delta 12 fatty acid desaturases (D12D). D9D and D12D were rate limiting enzymes in the biosynthesis of unsaturated fatty acid and polyunsaturated fatty acids, respectively. Analysis of gene expression profiling for D9D and D12D revealed that the four genes involved in linoleic acid biosynthesis were down-regulated in the response to oxidative stress, suggesting that unsaturated fatty acid synthesis was inhibited in A. oryzae cells in response to oxidative stress (Figure 6A). This was further supported by the determination of fatty acid content, showing a decrease in unsaturated fatty acid (Figure 6B).

4. Discussion

Environmental stress responses of A. oryzae cells are not isolated in linear metabolic or signaling pathways. Global responses have been associated with a much larger stress signaling network that integrates the information from various pathways and enzymes [6]. Therefore, a systematic analysis method is required to observe the intracellular global regulation patterns. Genome-wide approaches are considered as a convenient and efficient method to characterize gene expression profiles and elucidate the whole regulation network in plants, animals and microorganisms. The genome-wide projects of A. oryzae involved a lot of effort, and the complete genome sequence of A. oryzae was obtained through whole-genome shotgun approach in 2005 [2,33]. Wang et al. analyzed the transcriptomes of A. oryzae grown on solid-state culture, and in liquid culture with and without DTT treatment to understand the complex mechanism in A. oryzae at whole genome level [34]. A comprehensive transcriptome analysis of A. oryzae was performed to understand the possible biological processes across various growth stages [35]. In this study, genome-wide response of gene expressions in A. oryzae under oxidative conditions was analyzed to better understand the oxidative stress response of A. oryzae induced by H2O2 and MSB.
The results of this study showed that exposure of A. oryzae cells to H2O2 and MSB led to a substantial transcriptional regulation. However, the response of A. oryzae cells to H2O2 and MSB treatment was different. Only 56 DEGs were found to be commonly shared among the different concentration of H2O2 and MSB treatments (Figure 2B). In addition, there were less than 200 common DEGs shared between H2O2 and MSB treatments, whereas 525 DEGs were commonly shared between 3 mM and 6 mM H2O2 and 378 DEGs were commonly shared between 0.28 mM and 0.42 mM MSB treatment. Moreover, oxidative stress enhanced the expression of genes involved in oxidoreduction, including CAT, GPX, and SOD, and transcription factors (Yap1 and Skn7). The critical roles of these genes involved in oxidoreduction and transcription factors in oxidative stress have been reported in many microorganisms, such as yeast, Aspergillus fumigatus and Alternaria alternata [36,37]. In this study, the up-regulated expression of genes involved in oxidoreduction and two transcription factors were observed in A. oryzae cells in response to H2O2 and MSB treatment. However, the expression levels were different. MSB treatment increased GPX, whereas exogenous H2O2 enhanced Yap1 transcription factor. Further, seven Skn7 transcription factors exhibited significant increase in the gene expressions when treated with H2O2 and the remaining five Skn7 transcription factors showed up-regulated expressions at MSB treatments. These results demonstrated the different response of A. oryzae cells to H2O2 and MSB treatment. Previous studies on oxidative stress response of filamentous fungi induced by H2O2 and paraquat also showed that cell responses were different with different treatment, in which paraquat increased mainly SOD, whereas exogenous H2O2 increased CAT [38]. Although the expression patterns of these genes involved in oxidoreduction and transcription factors were different, GSH content was found to be high in both treatments. GSH is one of the major antioxidant molecules in yeast cells, and is considered to play a vital role in cell defense mechanism against ROS [39]. Glutathione can react with oxidants non-enzymically, or the reaction can be catalyzed by glutathione peroxidase, and oxidized glutathiol (GSSG) could be recycled by glutathione reductase (GR) at the expense of NADPH [22]. Mutants of S. cerevisiae deficient in glutathione synthesis or glutathione recycling showed increased sensitivity to H2O2 [40]. In A. niger, GSH concentration increased significantly following the addition of H2O2 when compared to control [24]. These reports supported the results in this study that GSH has essential role in oxidative defense. Besides GSH, the composition of lipids in cells has close relationship with the response to oxidative stress. Fatty acids of yeast cells are more resistant to oxidative attack than those of mammalian cells due to the lack of enzymes capable of synthesizing polyunsaturated fatty acids (PUFAs) [41]. Cipak et al. observed that the PUFA-producing yeast was initially more sensitive to oxidative stress than the wild-type strain because induced PUFA production increased the levels of endogenous ROS [42]. However, the changes of lipids in A. niger were different from that observed in yeasts, indicating that H2O2 treatment significantly increased the unsaturation degree of fatty acids [43]. In this study, it was observed that the changes of fatty acid compositions in A. oryzae under oxidative stress were highly similar with yeast, which exhibited an increased content of saturated fatty acids.

Author Contributions

Conceptualization, B.H. and Y.T.; methodology, C.J., H.S. and Y.W.; validation, Y.W., J.W. and L.M.; writing—original draft preparation, Z.H. and B.H.; supervision, B.Z.; project administration, H.S.; funding acquisition, B.Z. and B.H.

Funding

This research was funded by National Natural Science Foundation of China (NSFC) (Grant no. 31460447 and 31700068), International S&T Cooperation Project of Jiangxi Provincial (Grant no. 20142BDH80003), Doctor and master specific projects of Honghe University (XJ17B09), the Science Funds of Natural Science Foundation of Jiangxi Province (20114BAB205039).

Acknowledgments

The authors thank Wei Zhang and Jia Qin for critical reading of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. He, B.; Li, H.; Hu, Z.; Zhang, Y.; Sun, M.; Qiu, S.; Zeng, B. Difference in microbial community and taste compounds between Mucor-type and Aspergillus-type Douchi during koji-making. Food Res. Int. (Ott. Ont.) 2019, 121, 136–143. [Google Scholar] [CrossRef] [PubMed]
  2. Machida, M.; Asai, K.; Sano, M.; Tanaka, T.; Kumagai, T.; Terai, G.; Kusumoto, K.-I.; Arima, T.; Akita, O.; Kashiwagi, Y.; et al. Genome sequencing and analysis of Aspergillus oryzae. Nature 2005, 438, 1157. [Google Scholar] [CrossRef] [PubMed]
  3. He, B.; Tu, Y.Y.; Jiang, C.M.; Zhang, Z.; Li, Y.K.; Zeng, B. Functional Genomics of Aspergillus oryzae: Strategies and Progress. Microorganisms 2019, 7, 103. [Google Scholar] [CrossRef] [PubMed]
  4. Zadrąg-Tęcza, R.; Maślanka, R.; Bednarska, S.; Kwolek-Mirek, M. Response Mechanisms to Oxidative Stress in Yeast and Filamentous Fungi. In Stress Response Mechanisms in Fungi: Theoretical and Practical Aspects; Skoneczny, M., Ed.; Springer International Publishing: Cham, Switzerland, 2018. [Google Scholar]
  5. Kreiner, M.; Harveya, L.M.; Mcneila, B. Oxidative stress response of a recombinant Aspergillus niger to exogenous menadione and H2O2 addition. Enzym. Microb. Technol. 2002, 30, 346–353. [Google Scholar] [CrossRef]
  6. Zhao, H.; Chen, J.; Liu, J.; Han, B. Transcriptome analysis reveals the oxidative stress response in Saccharomyces cerevisiae. Rsc Adv. 2015, 5, 22923–22934. [Google Scholar] [CrossRef]
  7. Lópezmirabal, H.R.; Thorsen, M.; Kiellandbrandt, M.C.; Toledano, M.B.; Winther, J.R. Cytoplasmic glutathione redox status determines survival upon exposure to the thiol-oxidant 4,4′-dipyridyl disulfide. Fems Yeast Res. 2010, 7, 391–403. [Google Scholar] [CrossRef] [PubMed]
  8. Suo, Y.; Liu, Y.; Zhou, X.; Huang, Y.; Shi, C.; Matthews, K.; Shi, X. Impact of sod on the expression of stress-related genes in Listeria monocytogenes 4b G with/without paraquat treatment. J. Food Sci. 2015, 79, M1745–M1749. [Google Scholar] [CrossRef] [PubMed]
  9. Barker, M.G.; Brimage, L.J.E.; Smart, K.A. Effect of Cu,Zn superoxide dismutase disruption mutation on replicative senescence in Saccharomyces cerevisiae. Fems Microbiol. Lett. 2010, 177, 199–204. [Google Scholar] [CrossRef] [PubMed]
  10. Nishimoto, T.; Watanabe, T.; Furuta, M.; Kataoka, M.; Kishida, M. Roles of Catalase and Trehalose in the Protection from Hydrogen Peroxide Toxicity in Saccharomyces cerevisiae. Biocontrol Sci. 2016, 21, 179–182. [Google Scholar] [CrossRef] [PubMed]
  11. Grant, C.M.; Perrone, G.; Dawes, I.W. Glutathione and Catalase Provide Overlapping Defenses for Protection against Hydrogen Peroxide in the Yeast Saccharomyces cerevisiae. Biochem. Biophys. Res. Commun. 1998, 253, 893–898. [Google Scholar] [CrossRef] [PubMed]
  12. Dalto, D.B.; Matte, J.J. Pyridoxine (Vitamin B6) and the Glutathione Peroxidase System; a Link between One-Carbon Metabolism and Antioxidation. Nutrients 2017, 9, 189. [Google Scholar] [CrossRef] [PubMed]
  13. Mulford, K.E.; Fassler, J.S. Association of the Skn7 and Yap1 Transcription Factors in the Saccharomyces cerevisiae Oxidative Stress Response. Eukaryot. Cell 2011, 10, 761–769. [Google Scholar] [CrossRef] [PubMed]
  14. Carmel-Harel, O.; Stearman, R.; Gasch, A.P.; Botstein, D.; Brown, P.O.; Storz, G. Role of thioredoxin reductase in the Yap1p-dependent response to oxidative stress in Saccharomyces cerevisiae. Mol. Microbiol. 2010, 39, 595–605. [Google Scholar] [CrossRef] [PubMed]
  15. Chen, L.H.; Tsai, H.C.; Yu, P.L.; Chung, K.R. A Major Facilitator Superfamily Transporter-Mediated Resistance to Oxidative Stress and Fungicides Requires Yap1, Skn7, and MAP Kinases in the Citrus Fungal Pathogen Alternaria alternata. PLoS ONE 2017, 12, e0169103. [Google Scholar] [CrossRef] [PubMed]
  16. Yi, D.G.; Kim, M.J.; Ji, E.C.; Lee, J.; Jung, J.; Huh, W.K.; Chung, W.H. Yap1 and Skn7 genetically interact with Rad51 in response to oxidative stress and DNA double-strand break in Saccharomyces cerevisiae. Free Radic. Biol. Med. 2016, 101, 424–433. [Google Scholar] [CrossRef]
  17. Ilyas, S.; Rehman, A. Oxidative stress, glutathione level and antioxidant response to heavy metals in multi-resistant pathogen, Candida tropicalis. Environ. Monit. Assess. 2015, 187, 1–7. [Google Scholar] [CrossRef] [PubMed]
  18. Grant, C.M.; Maciver, F.H.; Dawes, I.W. Glutathione is an essential metabolite required for resistance to oxidative stress in the yeast Saccharomyces cerevisiae. Curr. Genet. 1996, 29, 511–515. [Google Scholar] [CrossRef]
  19. Kikugawa, K.; Yasuhara, Y.; Ando, K.; Koyama, K.; Hiramoto, K.; Suzuki, M. Protective Effect of Supplementation of Fish Oil with High n−3 Polyunsaturated Fatty Acids against Oxidative Stress-Induced DNA Damage of Rat Liver in Vivo. J. Agric. Food Chem. 2003, 51, 6073–6079. [Google Scholar] [CrossRef]
  20. Lagarde, M.; Guichardant, M.; Bernoudhubac, N.; Calzada, C.; Véricel, E. Oxygenation of polyunsaturated fatty acids and oxidative stress within blood platelets. Biochim. Et Biophys. Acta (Bba) - Mol. Cell Biol. Lipids 2018, 1863, 651–656. [Google Scholar] [CrossRef]
  21. Belozerskaya, T.A.; Gessler, N.N. Oxidative stress and differentiation in Neurospora crassa. Microbiology 2006, 75, 427–431. [Google Scholar] [CrossRef]
  22. Li, Q.; Mcneil, B.; Harvey, L.M. Adaptive response to oxidative stress in the filamentous fungus Aspergillus niger B1-D. Free Radic. Biol. Med. 2008, 44, 394–402. [Google Scholar] [CrossRef]
  23. Emri, T.; Szarvas, V.; Orosz, E.; Antal, K.; Park, H.S.; Han, K.H.; Yu, J.H.; Pócsi, I. Core oxidative stress response in Aspergillus nidulans. BMC Genom. 2015, 16, 478. [Google Scholar] [CrossRef]
  24. Emilie, F.; Cécile-Marie, A.D.; Laurence, D.; David, H.; Eduardo, D.C. Transcriptomic insights into the oxidative response of stress-exposed Aspergillus fumigatus. Curr. Pharm. Des. 2013, 19, 3713–3737. [Google Scholar]
  25. He, B.; Long, M.; Hu, Z.; Li, H.; Ai, M.; Long, C.; Zeng, B. Deep sequencing analysis of transcriptomes in Aspergillus oryzae in response to salinity stress. Appl. Microbiol. Biotechnol. 2018, 102, 897–906. [Google Scholar] [CrossRef]
  26. Vries, R.P.D.; Riley, R.; Wiebenga, A.; Aguilar-Osorio, G.; Amillis, S.; Uchima, C.A.; Anderluh, G.; Asadollahi, M.; Askin, M.; Barry, K. 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]
  27. Rahman, I.; Kode, A.; Biswas, S. Assay for quantitative determination of glutathione and glutathione disulfide levels using enzymatic recycling method. Nat. Protoc. 2007, 1, 3159–3165. [Google Scholar] [CrossRef]
  28. Wu, X.; Tong, Y.; Shankar, K.; Baumgardner, J.N.; Kang, J.; Badeaux, J.; Badger, T.M.; Ronis, M.J. Lipid fatty acid profile analyses in liver and serum in rats with nonalcoholic steatohepatitis using improved gas chromatography-mass spectrometry methodology. J. Agric. Food Chem. 2011, 59, 747–754. [Google Scholar] [CrossRef]
  29. Kim, D.; Salzberg, S.L. TopHat-Fusion: An algorithm for discovery of novel fusion transcripts. Genome Biol. 2011, 12, R72. [Google Scholar] [CrossRef]
  30. Li, B.; Dewey, C.N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011, 12, 323. [Google Scholar] [CrossRef]
  31. Zuo, Z.; Yang, L.; Chen, S.; Ye, C.; Han, Y.; Wang, S.; Ma, Y. Effects of nitrogen nutrients on the volatile organic compound emissions from Microcystis aeruginosa. Ecotoxicol. Environ. Saf. 2018, 161, 214–220. [Google Scholar] [CrossRef]
  32. Bai, Z.; Harvey, L.; Mcneil, B. Physiological responses of chemostat cultures of Aspergillus niger (B1-D) to simulated and actual oxidative stress. Biotechnol. Bioeng. 2003, 82, 691–701. [Google Scholar] [CrossRef]
  33. He, B.; Tu, Y.; Hu, Z.; Ma, L.; Dai, J.; Cheng, X.; Li, H.; Liu, L.; Zeng, B. Genome-wide identification and expression profile analysis of the HOG gene family in Aspergillus oryzae. World J. Microb. Biot. 2018, 34, 35. [Google Scholar] [CrossRef]
  34. Bin, W.; Guangwu, G.; Chao, W.; Ying, L.; Xiaoning, W.; Mouming, Z.; Yong, G.; Minghui, H.; Yong, Z.; Li, P. Survey of the transcriptome of Aspergillus oryzae via massively parallel mRNA sequencing. Nucleic Acids Res. 2010, 38, 5075–5087. [Google Scholar]
  35. He, B.; Hu, Z.; Long, M.; Li, H.; Ai, M.; Han, J.; Zeng, B. Transcriptome analysis of different growth stages of Aspergillus oryzae reveals dynamic changes of distinct classes of genes during growth. BMC Microbiol. 2018, 18, 12. [Google Scholar] [CrossRef]
  36. Rocha, M.C.; Godoy, K.F.D.; Bannitzfernandes, R.; Fabri, J.H.T.M.; Barbosa, M.M.F.; Castro, P.A.D.; Almeida, F.; Goldman, G.H.; Cunha, A.F.D.; Netto, L.E.S. Analyses of the three 1-Cys Peroxiredoxins from Aspergillus fumigatus reveal that cytosolic Prx1 is central to H2O2 metabolism and virulence. Sci. Rep. 2018, 8, 12314. [Google Scholar] [CrossRef]
  37. Wang, M.; Sun, X.; Yu, D.; Xu, J.; Chung, K.; Li, H. Genomic and transcriptomic analyses of the tangerine pathotype of Alternaria alternata in response to oxidative stress. Sci. Rep. 2016, 6, 32437. [Google Scholar] [CrossRef]
  38. Angelova, M.B.; Pashova, S.B.; Spasova, B.K.; Vassilev, S.V.; Slokoska, L.S. Oxidative stress response of filamentous fungi induced by hydrogen peroxide and paraquat. Mycol. Res. 2005, 109, 150–158. [Google Scholar] [CrossRef]
  39. Mitchell, D.A.; Marshall, T.K.; Deschenes, R.J. Vectors for the inducible overexpression of glutathione S-transferase fusion proteins in yeast. Yeast 2010, 9, 715–722. [Google Scholar] [CrossRef]
  40. Lee, J.C.; Straffon, M.J.; Jang, T.Y.; Higgins, V.J.; Grant, C.M.; Dawes, I.W. The essential and ancillary role of glutathione in Saccharomyces cerevisiae analysed using a grande gsh1 disruptant strain. Fems Yeast Res. 2001, 1, 57–65. [Google Scholar]
  41. Andrišić, L.; Collinson, E.J.; Tehlivets, O.; Perak, E.; Žarković, T.; Dawes, I.W.; Žarković, N.; Gasparović, A.Č. Transcriptional and antioxidative responses to endogenous polyunsaturated fatty acid accumulation in yeast. Mol. Cell. Biochem. 2015, 399, 27–37. [Google Scholar] [CrossRef]
  42. Cipak, A.; Jaganjac, M.; Tehlivets, O.; Kohlwein, S.D.; Zarkovic, N. Adaptation to oxidative stress induced by polyunsaturated fatty acids in yeast. Bba - Mol. Cell Biol. Lipids 2008, 1781, 283–287. [Google Scholar] [CrossRef] [PubMed]
  43. Ianutsevich, E.A.; Danilova, O.A.; Groza, N.V.; Tereshina, V.M. Membrane lipids and cytosol carbohydrates in Aspergillus niger under osmotic, oxidative, and cold impact. Microbiology 2016, 85, 302–310. [Google Scholar] [CrossRef]
Microorganisms 07 00225 g001 550
Figure 1. The density of spores, dry biomass and phenotype as affected by oxidative treatment. (A) The density of spores in the oxidative treatments (H2O2 and MSB) and the control; (B) The dry biomass in the oxidative treatments and the control. The mycelia were collected by peeling them off from the plates and dried overnight for the determination of biomass; The bars represent the average (±SE) of biological repeats. (C) The phenotype of A. oryzae in the control and the oxidative treatment at 72 h. The colony growth was retarded and the conidia production was inhibited. WT: the control; H2O2: the H2O2 treatments (3 and 6 mM); MSB: the MSB treatments (0.28 and 0.42 mM).
Figure 1. The density of spores, dry biomass and phenotype as affected by oxidative treatment. (A) The density of spores in the oxidative treatments (H2O2 and MSB) and the control; (B) The dry biomass in the oxidative treatments and the control. The mycelia were collected by peeling them off from the plates and dried overnight for the determination of biomass; The bars represent the average (±SE) of biological repeats. (C) The phenotype of A. oryzae in the control and the oxidative treatment at 72 h. The colony growth was retarded and the conidia production was inhibited. WT: the control; H2O2: the H2O2 treatments (3 and 6 mM); MSB: the MSB treatments (0.28 and 0.42 mM).
Microorganisms 07 00225 g001
Microorganisms 07 00225 g002 550
Figure 2. Distribution of differentially expressed genes (DEGs) in the five samples. (A) DEGs’ distribution between control and oxidative treatment; (B) Venn diagram exhibiting the DEGs’ distribution in five samples. WT: the control; H2O2: the H2O2 treatments (3 and 6 mM); MSB: the MSB treatments (0.28 and 0.42 mM).
Figure 2. Distribution of differentially expressed genes (DEGs) in the five samples. (A) DEGs’ distribution between control and oxidative treatment; (B) Venn diagram exhibiting the DEGs’ distribution in five samples. WT: the control; H2O2: the H2O2 treatments (3 and 6 mM); MSB: the MSB treatments (0.28 and 0.42 mM).
Microorganisms 07 00225 g002
Microorganisms 07 00225 g003 550
Figure 3. The expression profile of catalase (CAT), GPX1 (glutathione peroxidase) and superoxide dismutase (SOD) in the five samples. WT: the control; H2O2: the H2O2 treatments (3 and 6 mM); MSB: the MSB treatments (0.28 and 0.42 mM).
Figure 3. The expression profile of catalase (CAT), GPX1 (glutathione peroxidase) and superoxide dismutase (SOD) in the five samples. WT: the control; H2O2: the H2O2 treatments (3 and 6 mM); MSB: the MSB treatments (0.28 and 0.42 mM).
Microorganisms 07 00225 g003
Microorganisms 07 00225 g004 550
Figure 4. The expression profile of Yap1 and Skn7 transcription factors in the five samples. WT: the control; H2O2: the H2O2 treatments (3 and 6 mM); MSB: the MSB treatments (0.28 and 0.42 mM).
Figure 4. The expression profile of Yap1 and Skn7 transcription factors in the five samples. WT: the control; H2O2: the H2O2 treatments (3 and 6 mM); MSB: the MSB treatments (0.28 and 0.42 mM).
Microorganisms 07 00225 g004
Microorganisms 07 00225 g005 550
Figure 5. The content of glutathione (GSH) and expression profile of genes involved in glutathione biosynthesis in the oxidative treatments (H2O2 and MSB) and the control. (A) The standard curve of GSH concentration. (B) The changes of GSH level in response to oxidative stress. The bars represent the average (±SE) of biological repeats. The content of GSH increased with the increasing concentration of H2O2 and MSB. (C) The expression pattern of DEGs involved in the linoleic acid biosynthesis pathways (from left to right: the control, 3 mM and 6 mM H2O2, 0.28 mM and 0.42 mM MSB).
Figure 5. The content of glutathione (GSH) and expression profile of genes involved in glutathione biosynthesis in the oxidative treatments (H2O2 and MSB) and the control. (A) The standard curve of GSH concentration. (B) The changes of GSH level in response to oxidative stress. The bars represent the average (±SE) of biological repeats. The content of GSH increased with the increasing concentration of H2O2 and MSB. (C) The expression pattern of DEGs involved in the linoleic acid biosynthesis pathways (from left to right: the control, 3 mM and 6 mM H2O2, 0.28 mM and 0.42 mM MSB).
Microorganisms 07 00225 g005
Microorganisms 07 00225 g006 550
Figure 6. DEGs involved in the linoleic acid biosynthesis pathways and the changes of C18 fatty acid in response to oxidative stress. (A) The expression pattern of DEGs involved in the linoleic acid biosynthesis pathways (from left to right: the control, 3 mM and 6 mM H2O2, 0.28 mM and 0.42 mM MSB); (B) The changes of C18 fatty acid in response to oxidative stress. The bars represent the average (±SE) of biological repeats. WT: the control; H2O2: the H2O2 treatments (3 and 6 mM); MSB: the MSB treatments (0.28 and 0.42 mM).
Figure 6. DEGs involved in the linoleic acid biosynthesis pathways and the changes of C18 fatty acid in response to oxidative stress. (A) The expression pattern of DEGs involved in the linoleic acid biosynthesis pathways (from left to right: the control, 3 mM and 6 mM H2O2, 0.28 mM and 0.42 mM MSB); (B) The changes of C18 fatty acid in response to oxidative stress. The bars represent the average (±SE) of biological repeats. WT: the control; H2O2: the H2O2 treatments (3 and 6 mM); MSB: the MSB treatments (0.28 and 0.42 mM).
Microorganisms 07 00225 g006
Table
Table 1. Summary of the sequencing data and gene numbers of A. oryzae transcriptome under oxidative stress.
Table 1. Summary of the sequencing data and gene numbers of A. oryzae transcriptome under oxidative stress.
SamplesWT3 mM H2O26 mM H2O20.28 mM MSB0.42 mM MSB
Clean Reads40,788,83241,412,22841,102,85043,613,78239,908,124
GC Content52.28%52.47%52.26%51.64%51.88%
% ≥ Q3094.27%94.33%94.40%94.36%94.35%
Mapped reads36,552,991 (89.62%)34,900,808 (84.28%)34,852,216 (84.79%)38,343,753 (87.92%)35,212,300 (88.23%)
Unique mapped36,452,891 (89.37%)34,787,458 (84.00%)34,722,583 (84.48%)38,241,527 (87.68%)35,122,814 (88.01%)
Multiple mapped100,100 (0.25%)113,350 (0.27%)129,633 (0.32%)102,226 (0.23%)89,486 (0.22%)

Share and Cite

MDPI and ACS Style

Shao, H.; Tu, Y.; Wang, Y.; Jiang, C.; Ma, L.; Hu, Z.; Wang, J.; Zeng, B.; He, B. Oxidative Stress Response of Aspergillus oryzae Induced by Hydrogen Peroxide and Menadione Sodium Bisulfite. Microorganisms 2019, 7, 225. https://doi.org/10.3390/microorganisms7080225

AMA Style

Shao H, Tu Y, Wang Y, Jiang C, Ma L, Hu Z, Wang J, Zeng B, He B. Oxidative Stress Response of Aspergillus oryzae Induced by Hydrogen Peroxide and Menadione Sodium Bisulfite. Microorganisms. 2019; 7(8):225. https://doi.org/10.3390/microorganisms7080225

Chicago/Turabian Style

Shao, Huanhuan, Yayi Tu, Yijing Wang, Chunmiao Jiang, Long Ma, Zhihong Hu, Jiangfan Wang, Bin Zeng, and Bin He. 2019. "Oxidative Stress Response of Aspergillus oryzae Induced by Hydrogen Peroxide and Menadione Sodium Bisulfite" Microorganisms 7, no. 8: 225. https://doi.org/10.3390/microorganisms7080225

APA Style

Shao, H., Tu, Y., Wang, Y., Jiang, C., Ma, L., Hu, Z., Wang, J., Zeng, B., & He, B. (2019). Oxidative Stress Response of Aspergillus oryzae Induced by Hydrogen Peroxide and Menadione Sodium Bisulfite. Microorganisms, 7(8), 225. https://doi.org/10.3390/microorganisms7080225

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