Small GTPases and Stress Responses of vvran1 in the Straw Mushroom Volvariella volvacea

Small GTPases play important roles in the growth, development and environmental responses of eukaryotes. Based on the genomic sequence of the straw mushroom Volvariella volvacea, 44 small GTPases were identified. A clustering analysis using human small GTPases as the references revealed that V. volvacea small GTPases can be grouped into five families: nine are in the Ras family, 10 are in the Rho family, 15 are in the Rab family, one is in the Ran family and nine are in the Arf family. The transcription of vvran1 was up-regulated upon hydrogen peroxide (H2O2) stress, and could be repressed by diphenyleneiodonium chloride (DPI), a NADPH oxidase-specific inhibitor. The number of vvran1 transcripts also increased upon cold stress. Diphenyleneiodonium chloride, but not the superoxide dismutase (SOD) inhibitor diethy dithiocarbamate (DDC), could suppress the up-regulation of vvran1 gene expression to cold stress. These results combined with the high correlations between gene expression and superoxide anion (O2−) generation indicated that vvran1 could be one of the candidate genes in the downstream of O2− mediated pathways that are generated by NADPH oxidase under low temperature and oxidative stresses.


Introduction
Small GTPases, a group of guanine nucleotide binding protein monomers with molecular weights ranging from 20 to 30 kDa, widely exist in eukaryotic cells [1]. Small GTPases (sometimes called Ras superfamily GTPases) can bind to GTP and GDP, switching between the activated (GTP binding) and non-activated (GDP binding) states depending on the binding of different substrates, thereby turning on or off relevant metabolic pathways [2]. The Ras superfamily can be further divided into five families (Ras, Rab, Rho, Ran and Arf) depending on the structures and functions of each constituent [3]. The Ras family is important for cell proliferation, differentiation, apoptosis and the regulation of gene expression [4]; the Rho family is involved in cytoskeleton formation, cell polarity, the cell cycle, the regulation of gene expression and the process of mating in yeast cells [5][6][7]; the Rab and Arf families play important roles in transport across cell membranes [8][9][10][11][12]; and the Ran family regulates the transportation of proteins and RNA molecules at the nuclear pore and is crucial for nuclear assembly, spindle formation and the regulation of mitosis [13][14][15].
To date, fungal small GTPases have been annotated and classified in only a few fungi, including Coprinopsis cinerea, Cryptococcus neoformans, Laccaria bicolor, Ganoderma lucidum, Phanerochaete chrysosporium, "_G" represent for GENSCAN Prediction; "_g" represent for GeneMark-ES Prediction; "_A" represent for Augustus Prediction; "_T" means the gene intron and exon regions were confirmed by transcriptom data, and the amino acid sequences were identified by ORF finder software.

Classification of V. volvacea Small GTPases
The amino acid sequences of 151 human small GTPases identified by Rojas et al. [3] were downloaded from the RCSB Protein Data Bank [26]. Among these sequences, there were 39 Ras proteins, 22 Rho proteins, 60 Rab proteins, one Ran protein and 29 Arf proteins. The phylogenic analysis of the 44 predicted V. volvacea small GTPases using the 151 human small GTPases as references showed that all human and V. volvacea small GTPases were clustered into four large clades ( Figure 1). The Ran and Rab family proteins were relatively close and were clustered together in the same clade. The 44 predicted V. volvacea small GTPases could be divided into the Ras, Rab, Rho, Ran and Arf families. There were one Ras proteins, 10 Rho proteins (of which one, VvMitRho, was a mitochondrial protein), 15 Rab proteins, one Ran protein and nine Arf proteins.
"_G" represent for GENSCAN Prediction; "_g" represent for GeneMark-ES Prediction; "_A" represent for Augustus Prediction; "_T" means the gene intron and exon regions were confirmed by transcriptom data, and the amino acid sequences were identified by ORF finder software.

Classification of V. volvacea Small GTPases
The amino acid sequences of 151 human small GTPases identified by Rojas et al. [3] were downloaded from the RCSB Protein Data Bank [26]. Among these sequences, there were 39 Ras proteins, 22 Rho proteins, 60 Rab proteins, one Ran protein and 29 Arf proteins. The phylogenic analysis of the 44 predicted V. volvacea small GTPases using the 151 human small GTPases as references showed that all human and V. volvacea small GTPases were clustered into four large clades ( Figure 1). The Ran and Rab family proteins were relatively close and were clustered together in the same clade. The 44 predicted V. volvacea small GTPases could be divided into the Ras, Rab, Rho, Ran and Arf families. There were one Ras proteins, 10 Rho proteins (of which one, VvMitRho, was a mitochondrial protein), 15 Rab proteins, one Ran protein and nine Arf proteins.  [3]. The branches were designated as protein name_uniprot ID. The clades of 44 V. volvacea small GTPases were labeled with the protein name_GenBank ID. The confidence levels of nodes were tested by bootstrapping 1000 times; scores ≥80% were denoted. Phanerochaete Chrysosporium RP78 (v2.0; February 2005; JGI), were also annotated by using the same method and chosen for the clustering analysis using the neighbor joining method. Results showed that each basidiomycete fungus contains all types of small GTPases ( Figure 2; Table 2). The numbers of GTPases in the Rab family are consistent between the different basidiomycetes. The Rab numbers are also relatively greater than those of the GTPases in other families of GTPases in the basidiomycetes except for Laccaria bicolor. Phanerochaete Chrysosporium RP78 (v2.0; February 2005; JGI), were also annotated by using the same method and chosen for the clustering analysis using the neighbor joining method. Results showed that each basidiomycete fungus contains all types of small GTPases ( Figure 2; Table 2). The numbers of GTPases in the Rab family are consistent between the different basidiomycetes. The Rab numbers are also relatively greater than those of the GTPases in other families of GTPases in the basidiomycetes except for Laccaria bicolor. Phylogenetic analysis of small GTPases of 5 basidiomycetes. The clades were labeled with the abbreviated Latin species name, followed by "_predicted gene ID" except for the proteins from Volvariella volvacea, which were labeled with "abbreviated Latin species name_protein name". The following abbreviations were used: Sccom for Schizophyllum commune; Labic for Laccaria bicolor; Cocin for Coprinopsis cinerea; Pcchr for Pchrysosporium chrysosporium; and Vovol for Volvariella volvacea. The confidence levels of the nodes were tested by bootstrapping 1000 times; scores ≥80% were denoted.

Figure 2.
Phylogenetic analysis of small GTPases of 5 basidiomycetes. The clades were labeled with the abbreviated Latin species name, followed by "_predicted gene ID" except for the proteins from Volvariella volvacea, which were labeled with "abbreviated Latin species name_protein name". The following abbreviations were used: Sccom for Schizophyllum commune; Labic for Laccaria bicolor; Cocin for Coprinopsis cinerea; Pcchr for Pchrysosporium chrysosporium; and Vovol for Volvariella volvacea. The confidence levels of the nodes were tested by bootstrapping 1000 times; scores ≥80% were denoted.

Phylogenic Analysis of VvRan1
The V. volvacea small GTPase in Ran family was termed VvRan1 in this study. To determine the relationship between VvRan1 and other Ran proteins, the amino acid sequences of VvRan1 (this study) and nine other Ran family proteins, which were downloaded from the RCSB Protein Data Bank were used for the sequence alignment. Figure 3 shows that the Ran family proteins were highly conserved in plants, animals and fungi. All of these sequences had five G box domains, two effecter regions [27] and an acidic C terminal sequence. The C-terminus of each plant Ran protein had two additional amino acid residues compared to those from animals and fungi, suggesting a close molecular evolutionary relationship between the Ran proteins of fungi and animals. To further determine the phylogenic relationship between VvRan1 and other fungal Ran GTPases, the protein sequences of VvRan1 and 22 Ran GTPases from 17 other fungal species were clustered using MEGA5.1 [28]. Figure 4 shows that the Ran GTPases clustered into two large clades. One clade contained only the Ran GTPases from Ascomycota, whereas the other clade contained only the Ran GTPases from Basidiomycota. In the Ascomycota clade, yeast Ran GTPases clustered together, while the Ran proteins of filamentous fungi clustered together. Compared with the Ran proteins of Rhodosporidium toruloides and Ustilago maydis, VvRan1 had a closer relationship with the Ran GTPases of Agaricales fungi. However, VvRan1 could not be clustered into the subclade of the Ran proteins of Agaricales fungi.

Phylogenic Analysis of VvRan1
The V. volvacea small GTPase in Ran family was termed VvRan1 in this study. To determine the relationship between VvRan1 and other Ran proteins, the amino acid sequences of VvRan1 (this study) and nine other Ran family proteins, which were downloaded from the RCSB Protein Data Bank were used for the sequence alignment. Figure 3 shows that the Ran family proteins were highly conserved in plants, animals and fungi. All of these sequences had five G box domains, two effecter regions [27] and an acidic C terminal sequence. The C-terminus of each plant Ran protein had two additional amino acid residues compared to those from animals and fungi, suggesting a close molecular evolutionary relationship between the Ran proteins of fungi and animals. To further determine the phylogenic relationship between VvRan1 and other fungal Ran GTPases, the protein sequences of VvRan1 and 22 Ran GTPases from 17 other fungal species were clustered using MEGA5.1 [28]. Figure 4 shows that the Ran GTPases clustered into two large clades. One clade contained only the Ran GTPases from Ascomycota, whereas the other clade contained only the Ran GTPases from Basidiomycota. In the Ascomycota clade, yeast Ran GTPases clustered together, while the Ran proteins of filamentous fungi clustered together. Compared with the Ran proteins of Rhodosporidium toruloides and Ustilago maydis, VvRan1 had a closer relationship with the Ran GTPases of Agaricales fungi. However, VvRan1 could not be clustered into the subclade of the Ran proteins of Agaricales fungi.

The Superoxide Anion (O2 − ) Signal Molecular Triggers vvran1 to Respone Stresses
Ran proteins are involved in the stress responses of plants [22], humans [23], and yeast [24]. To test the relationship between vvran1 gene and stresses in Volvariella volvacea, the mycelia were treated with 10 mmol/L hydrogen peroxide (H2O2, as oxidative stress) or 4 °C (as cold stress) for 1 h, respectively. RT-qPCR results showed that the vvran1 transcript levels were 4.7-and 6.6-fold increased, respectively. These results indicated that vvran1 gene could be regulated by both oxidative and cold stresses ( Figure 5A,B).
It is known that NADPH oxidase acts as the critical role in cellular stress responses [29]. It can produce the O2 − and then converted into H2O2 by superoxide dismutase (SOD), both of them can activate specific stresses signaling, also called "redox signaling" [30]. To further understand the vvran1 gene up-regulated mechanism, an NADPH oxidase-specific inhibitor, diphenyleneiodonium chloride (DPI) was used to inactivate the NADPH oxidases [31,32]. RT-qPCR results showed that DPI pre-incubating could completely repress the expression of vvran1 which was up-regulated by H2O2 and 4 °C stresses ( Figure 5A,B). These results could be taken as evidence that redox signaling via activation of NADPH oxidase is a reason for vvran1 up-regulated.
Moreover, there are two effects that exogenous H2O2 can cause. One is to permeate through the plasma membrane into cell and activate signaling stresses directly [33], the other is to act as oxidative stress and trigger the cell response by NADPH oxidase. Because the DPI can only block the NADPH oxidase path but cannot prevent H2O2 into cellular, the positive result of DPI treatment under H2O2 stress suggested that specific signaling stress activated by O2 − should be the reason for vvran1 up-regulated ( Figure 5A). To further confirm this mechanism, Diethy dithiocarbamate (DDC), a SOD specific inhibitor that can keep the high level of O2 − but the low level of H2O2 in the cell [34], was added to the incubation solution during cold stress. The results showed that both the intracellular O2 − concentration and the vvran1 expression were maintain at the high level ( Figure 5B,C). These results further confirmed that O2 − but not H2O2 signal mediates the vvran1 up-regulated expression under stress.  Ran proteins are involved in the stress responses of plants [22], humans [23], and yeast [24]. To test the relationship between vvran1 gene and stresses in Volvariella volvacea, the mycelia were treated with 10 mmol/L hydrogen peroxide (H 2 O 2 , as oxidative stress) or 4 • C (as cold stress) for 1 h, respectively. RT-qPCR results showed that the vvran1 transcript levels were 4.7-and 6.6-fold increased, respectively. These results indicated that vvran1 gene could be regulated by both oxidative and cold stresses ( Figure 5A,B).
It is known that NADPH oxidase acts as the critical role in cellular stress responses [29]. It can produce the O 2 − and then converted into H 2 O 2 by superoxide dismutase (SOD), both of them can activate specific stresses signaling, also called "redox signaling" [30]. To further understand the vvran1 gene up-regulated mechanism, an NADPH oxidase-specific inhibitor, diphenyleneiodonium chloride (DPI) was used to inactivate the NADPH oxidases [31,32]. RT-qPCR results showed that DPI pre-incubating could completely repress the expression of vvran1 which was up-regulated by H 2 O 2 and 4 • C stresses ( Figure 5A,B). These results could be taken as evidence that redox signaling via activation of NADPH oxidase is a reason for vvran1 up-regulated. Moreover, there are two effects that exogenous H 2 O 2 can cause. One is to permeate through the plasma membrane into cell and activate signaling stresses directly [33], the other is to act as oxidative stress and trigger the cell response by NADPH oxidase. Because the DPI can only block the NADPH oxidase path but cannot prevent H 2 O 2 into cellular, the positive result of DPI treatment under H 2 O 2 stress suggested that specific signaling stress activated by O 2 − should be the reason for Additionally, the Pearson correlation coefficient method suggested that the integrated optical density (IOD) mean value of nitroblue tetrazolium (NBT) straining showed in Figure 5C was correlated with vvran1 gene expression (r = 0.901, p < 0.05), which also indicated that the transcription of vvran1 regulated by oxidative and cold stresses may be mediated by O2 − signal molecules. relative expression levels of vvran1 to cold stress; and (C) the integrated optical density (IOD) mean of superoxide anion (O2 − ) detection by nitroblue tetrazolium (NBT) straining. The control, hyphae were incubated in phosphate buffered solution (PBS) buffer for 1.5 h; the H2O2 and H2O2 + DPI treatment, hyphae were incubated in PBS buffer or PBS buffer with 50 μmol/L DPI for 0.5 h, respectively, and then were switched to PBS buffer containing 10 mmol/L H2O2 and incubated for 1 h; the 4 °C and 4 °C + DPI treatment, hyphae were incubated in PBS buffer or PBS buffer with 50 μmol/L DPI for 0.5 h at 34 °C, respectively, followed by replacing the warm PBS buffer with cold PBS buffer (4 °C) and incubating the hyphae at 4 °C for 1 h; and the 4 °C + DDC treatment, hyphae were incubated for 0.5 h in PBS buffer at 34 °C, followed by exchanging the warm PBS buffer with cold PBS buffer (4 °C) containing DDC and incubating the hyphae at 4 °C for 1 h. Relative expression levels of vvran1 in (A,B) were calculated relative to the transcript level of vvran1 in the control and three independent experiments with nine independent replicates are shown by different color of shapes. The values of IOD mean in (C) are the means ± standard deviation (n = 3). Statistical testing of significance was performed using a one-way ANOVA and a Bonferroni's posttest, NS means no significance, and the different letters over the columns within a graph denote significant differences (p < 0.05).

Discussion
The cluster analysis using the neighbor joining (NJ) method grouped the human and V. volvacea small GTPases into five well-defined families, indicating that small GTPases within the same family are highly conserved between humans and microbes. However, there were also some highly variable amino acid residues despite the conservation of certain domains, such as the G box, Additionally, the Pearson correlation coefficient method suggested that the integrated optical density (IOD) mean value of nitroblue tetrazolium (NBT) straining showed in Figure 5C was correlated with vvran1 gene expression (r = 0.901, p < 0.05), which also indicated that the transcription of vvran1 regulated by oxidative and cold stresses may be mediated by O 2 − signal molecules.

Discussion
The cluster analysis using the neighbor joining (NJ) method grouped the human and V. volvacea small GTPases into five well-defined families, indicating that small GTPases within the same family are highly conserved between humans and microbes. However, there were also some highly variable amino acid residues despite the conservation of certain domains, such as the G box, between different families of small GTPases. This finding is consistent with previous results [3]. Because the yeast genome is small, it contains only 29 small GTPases [16]. A previous study showed that using only yeast as the reference for classification results in some of the small GTPases in other filamentous fungi being excluded from their appropriate families [17]. Because human small GTPases have been clearly and completely annotated and classified, we used human GTPases as the references and showed that the clustering analysis based on human small GTPases is suitable for the classification of small GTPases in filamentous fungi with large genomes.
Although Ran is the most abundant small GTPase in eukaryotes and is involved in many cellular metabolic activities, such as the assembly of nuclei, the formation of the spindle and the regulation of mitosis [4,[13][14][15], each basidiomycete contains only one gene encoding a Ran GTPase. In addition, each basidiomycete contains only one Cdc42, one Rac, and one mitochondrial Rho GTPase (Figure 2). The discovery of these genes will be useful for further studies. It is worth noting that Cdc42 and Rac are evolutionarily close, which might explain the similar cellular functions of these two genes [35][36][37].
Xu and Cai [38] found that in rice, low temperature stimulation could significantly up-regulate the expression of the Ran gene OsRAN1 and that the over-expression of OsRAN1 could effectively improve cold tolerance in rice. Our results showed that the vvran1 gene is also sensitive to the cold stress. Therefore, the over-expression of the vvran1 gene may be helpful in developing cold-tolerant mushrooms. Previous studies have shown that environmental stresses can activate membrane-bound NADPH oxidase to produce active ROS signaling molecules (such as O 2 − and H 2 O 2 ), thereby inducing relevant genes involved in protecting cells from stress [30,39]. Yan et al. [40] found that H 2 O 2 can up-regulate the expression of Ran/TC4 in benign mammary epithelial cells but not in malignant cells, suggesting that Ran/TC4 is involved in the antioxidant response of normal cells. Recent studies have suggested that oxidative stress can not only disrupt the distribution of Ran protein in the cell, but can also regulate the Ran-related cellular signal transduction pathways [41,42]. Based on the increased expression of vvran1 under low temperature and oxidative stresses and the high correlations between gene expression and O 2 − content, we propose that vvran1 could be one of candidate genes in the downstream of O 2 − mediated pathways which was produced by NADPH oxidase after stimulated by these abiotic stresses. According to many studies in the literature [39,[43][44][45][46], intracellular H 2 O 2 is an important signal molecular that regulates gene expression in the response to environmental stresses. However, the positive result of the experiment using the DPI under oxidative stress (H 2 O 2 stress) and the negative result of the DDC under cold stress suggested that vvran1 may not be regulated by intracellular H 2 O 2 . It has been reported that the expression of both the OsRAN1 and OsRAN2 genes in rice can be up-regulated by cold stresses [38,47]. These genes can maintain cell division and the progression of the cell cycle by promoting the formation of an intact nuclear envelope and promoting the export of intranuclear tubulin, thereby enhancing the cold tolerance of the cell [38,47]. Some other studies have suggested that the abiotic stresses such as free radical nitric oxide production and oxidative stress can mediate Ras guanine nucleotide dissociation; this decreases the levels of intracellular RanGTP and changes its cytoplasmic distribution, thereby leading to cell death [23,41,48]. Furthermore, classical nuclear protein import can be inhibited by oxidative and other forms of stress by reducing the GTP/GDP ratio in Saccharomyces cerevisiae [24]. Our results showed that the expression of vvran1 could be rapidly up-regulated by cold and oxidative stresses. This may promote nucleocytoplasmic transport, thereby enhancing the ability of the cells to tolerate stress. Needless to say, more studies are needed to reveal the detailed role of the vvran1 gene to abiotic stresses.

Strains
The dikaryotic strain H1521 (collection number: ACCC52633) was used in all experiments in this study. H1521 is a heterokaryon strain generated by crossing PYd15 (ACCC52631) with PYd21 (ACCC52632), two homokaryon strains with opposite mating types.

Genome Sequencing, Splicing and Prediction
De novo sequencing of the whole genome of PYd21 was performed on the Solexa/Illumina platform at the Shenzhen Huada Gene Research Institute (Shenzhen, China). The genome was assembled using a SOAPdenovo assembler [49]. The NCBI accession number PRJNA171553 was assigned to the genome. A total of 11,534 encoding genes and deduced amino acid sequences were obtained using GeneMark-ES (version 2.3, Atlanta, GA, USA) [50,51].

Annotation of Small GTPases
RAS-related protein sequences were downloaded from the RAS Oncogene Database [25] and locally compared to the 11,534 amino acid sequences using BLASTP after a standardized library was constructed. Amino acid sequences with identities ≥30% and e-values ≤1 × 10 −2 were extracted using Perl scripts and submitted to Pfam for domain prediction. Amino acid sequences containing PF00071 or PF00025 and e-values ≤1 × 10 −10 were defined as small GTPases.

Validation of Sequence Accuracy
To identify genes encoding small GTPases in V. volvacea, the DNA sequences of all V. volvacea small GTPase coding sequences along with the 1000 bp upstream and downstream sequences were used as references to map the reads in 500 bp read pools using the ZOOM software [52]. The reads were obtained from genome sequencing, and the paired end method was used to validate the accuracy of the sequences. The software parameter settings were as follows: the distance of adjacent paired reads was set at 1 to 2000 bp; the data were presented in the Illumina format; the number of allowed mismatch bases was set to 0; and other parameters were set at default. If all base pairs of a gene were covered by reads, the sequence was considered accurate; if some of the reference sequences did not have corresponding reads, the sequences were verified using Sanger sequencing [53]. All the corrected gene sequences were used as references to map against transcriptome sequence raw reads using ZOOM software, and the software parameter settings were the same with above but the number of allowed mismatch bases was set to 40 to identify the intron region, then, predicted the amino acid sequences by ORF finder online software. If the number of mapping reads was not enough for intron identify, we predicted the gene sequence again using GeneMark-ES [50].
Finally, the integrity and accuracy of the validated amino acid sequences were determined by submitting these sequences to NCBI for BLASTP analysis. The sequences that were obviously longer than the sequences in other species or the sequences that did not align to the sequences of other species were considered to be erroneously predicted, and their extra segments were removed. For sequences that lacked intact conserved domains, we used the GENSCAN (using Vertebrate as the reference species) and Augustus (with Laccaria bicolor as the reference species) websites to predict the alignment again [54,55].

Sequence Homology Comparison and Phylogenetic Tree Construction
After using the MUSCLE program to align the sequences, a neighbor-joining tree was constructed with MEGA5.1 [28]. The structures of the conserved amino acid sequences were colored using the GeneDoc software [56].

Preparation of Solutions and Stress Treatments
According to our previous research, 10 mmol/L H 2 O 2 or 4 • C treated for 1 h can significantly reduce but not completely inhibit the mycelium growth of Volvariella volvacea strain H1521, and hyphal growth was the fastest at pH 8 condition [57].
PBS buffer (0.02 mol/L, pH 8) containing 10 mmol/L of H 2 O 2 was used as an oxidative stress solution. DPI (Sigma, Saint Louis, MO, USA) and DDC (Sigma, Saint Louis, MO, USA) were dissolved in sterile water and diluted to 50 µmol/L and 1 mmol/L with PBS buffer (0.02 mol/L, pH 8), respectively. Solutions for cold treatments were kept at 4 • C while solutions for other treatments were kept at 34 • C before use. The temperature mentioned in this article were allowed within 0.5 • C fluctuation.
The V. volvacea strain H1521 was used to test the transcriptional models of vvran1 to different stresses. The mycelia were cultured in the solid potato dextrose agar (PDA) medium with glass papers on the surfaces of the PDA plates (Φ = 6 cm) for three days of incubation at 34 • C in the dark. To expose the mycelia to different stresses, about ten microliters of sterile PBS buffer (0.02 mol/L, pH 8) was first added to each plate to completely submerge the colonies in the buffer. The plates were incubated at 34 • C for an additional 0.5 h before stress treatment. For the DPI inhibition experiments, the colonies were pre-incubated with the DPI inhibitor solution that was added to the plates instead of PBS buffer. All stress treatments (except for cold stress treatments) were conducted at 34 • C in the dark for 1 h. For the cold stress treatments, the plates were incubated at 4 • C in the dark for 1 h.

vvran1 Transcript Analysis
After the stress treatment, the mycelia were quickly scraped, blot-dried and stored in a −80 • C freezer. The RNA was extracted using an E.Z.N.A.™ Plant RNA kit (Omega Bio-Tek, Norcross, GA, USA). The first strands of cDNA were synthesized using TransScript ® All-in-One First-Strand cDNA Synthesis SuperMix for qPCR (One-Step gDNA Removal) (TransGen Biotech, Beijing, China). Real-time fluorescent quantitative PCR (RT-qPCR) was carried out using TransStart Top Green qPCR SuperMix (TransGen Biotech, Beijing, China) on a CFX96 real-time fluorescence quantitative PCR machine (Bio-Rad, Hercules, CA, USA). The level of expression of the untreated control was used as the reference for calculating the relative expression levels using the 2 −∆∆Ct method [58]. Glyceraldehyde-3-dehydrogenase (GAPDH) and 18S ribosomal RNA gene (18S rRNA) were used as reference genes. The PCR primers are shown in Table 3. Table 3. Primers used in real-time quantitative PCR.

Primer
Sequence (5 -3 ) Solon, OH, USA) as substrate [59]. Briefly, the mycelia on the surfaces of the PDA plates after stresses treatment were killed quickly by liquid nitrogen, after ice melting, the mycelia were incubated with 0.05 mol/L PBS (pH 7.5) containing 0.05% NBT for 2 h at the ice-bath condition. The pictures were taken by Nikon P500 digital camera with the same exposure conditions. The computer-assisted genuine color image analysis system (imagepro-plus 6.0) was used to quantify the mean of integrated optical density.

Statistical Analysis
The significance of gene expression and superoxide anion content among different samples were analysed using the one-way ANOVA of variance with Bonferroni's multiple comparisons test, and the analysis was performed by GraphPad Prism version 5.01 (San Diego, CA, USA). The correlations of gene expression patterns and O 2 − contents were analyzed using the Pearson correlation coefficient method by SPSS Statistics v20 software with two-tailed test.