Deep sequencing analysis of transcriptomes in Aspergillus flavus in response to resveratrol

Background Resveratrol has been reported as a natural phytoalexin that inhibits infection or the growth of certain fungi including Aspergillus flavus. Our previous research revealed that aflatoxin production in A. flavus was reduced in medium with resveratrol. To understand the molecular mechanism of the A. flavus response to resveratrol treatment, the high-throughput paired-end RNA-Seq was applied to analyze the transcriptomic profiles of A. flavus. Results In total, 366 and 87 genes of A. flavus were significantly up- and down- regulated, respectively, when the fungus was treated with resveratrol. Gene Ontology (GO) functional enrichment analysis revealed that 48 significantly differentially expressed genes were involved in 6 different terms. Most genes in the aflatoxin biosynthetic pathway genes cluster (#54) did not show a significant change when A. flavus was treated with resveratrol, but 23 of the 30 genes in the #54 cluster were down-regulated. The transcription of aflA and aflB was significantly suppressed under resveratrol treatment, resulting in an insufficient amount of the starter unit hexanoate for aflatoxin biosynthesis. In addition, resveratrol significantly increased the activity of antioxidative enzymes that destroy radicals, leading to decreased aflatoxin production. Moreover, stuA, fluG, flbC, and others genes involved in mycelial and conidial development were down-regulated, which disrupted the cell’s orderly differentiation and blocked conidia formation and mycelia development. The transcripts of laeA and veA were slightly inhibited by resveratrol, which may partly decrease aflatoxin production and depress conidia formation. Conclusions Resveratrol can affect the expression of A. flavus genes that are related to developmental and secondary metabolic processes, resulting in decreased aflatoxin production and conidia formation and could also cause abnormal mycelia development. These results provide insight into the transcriptome of A. flavus in response to resveratrol and a new clew for further study in regulation of aflatoxin biosynthesis in A. flavus. Electronic supplementary material The online version of this article (doi:10.1186/s12866-015-0513-6) contains supplementary material, which is available to authorized users.


Background
The widespread fungus Aspergillus flavus is a saprophytic and opportunistic pathogen that is capable of causing diseases in agricultural crops [1], producing toxic secondary metabolites, and bringing about health problems in animals and humans [2]. Under favorable conditions, A. flavus produces polyketide-derived carcinogenic and mutagenic secondary metabolites, aflatoxins [3], which are extremely hepatotoxic, immunosuppressive, and antinutritional mycotoxins associated with both acute and chronic toxicity in humans and animals [4]. Besides the health implications in humans and animals, A. flavus colonization in crops causes considerable economic losses because of reduced utilization and lower price of aflatoxincontaminated grains.
To develop effective means of combating aflatoxin contamination, it is of vital importance to investigate the molecular mechanisms of development and secondary metabolism in A. flavus. Previous studies identified numerous diverse biochemical substances and extracts with adverse effects on the development and/or secondary metabolism of A. flavus [5][6][7][8][9][10]. Resveratrol (3,5,4-trihydroxystilbene) is a phytoalexin with antifungal activity and antioxidant capacity, which can enhance host plant resistance against biotic and abiotic stresses [11]. Our previous study found that resveratrol is highly related to resistance to aflatoxin production in peanut seeds and could inhibit aflatoxin production in medium [12]. In addition, resveratrol affected mycelial morphologic characteristics and inhibited conidia formation of A. flavus in medium. However, how resveratrol affects the development and secondary metabolism of A. flavus remains unknown.
RNA-sequencing (RNA-seq) is a powerful and costefficient high-throughput technology for transcriptomic profiling analysis that has been successfully used to interrogate transcriptomes of numerous fungi [10,[13][14][15][16][17]. The application of this technology has greatly accelerated our understanding of the complexity of gene expression, regulation, and networks [18]. Gene expression can be more accurately quantified using RNA-seq approaches than conventional transcriptomic analysis [18]. With higher sensitivity, RNA-seq can also be more efficient in detecting a larger range of dynamically expressed genes than microarrays. Furthermore, RNAseq has been used to survey sequence variations and complex transcriptomes with low false-positive rates and high sensitivity and reproducibility [14,19]. For an organism with a well-annotated genome, mapping read sequences to the corresponding reference genome/transcriptome is usually the first and essential step for RNA-seq data analysis [10]. Moreover, the "referencebased" RNA-seq approach was proved to be more sensitive and have higher computational efficiency than a de novo assembly approach [20]. The construction of expressed sequences tags in A. flavus was reported by several researchers [21], and the whole-genome sequencing of the fungus was completed [22]. Annotation of the A. flavus genome revealed numerous genes and gene clusters that are potentially involved in conidia formation and secondary metabolic processes including aflatoxin biosynthesis [23,24]. The transcriptomic analysis of aflatoxin biosynthesis and mycelia development in A. flavus in response to 5-azacytidine [10,25], decanal [17], water activity [26], and temperature [27] also used RNA-seq technology.
To comprehensively understand the effect of resveratrol on the development and secondary metabolism of A. flavus, an RNA-seq approach was applied in this study to obtain transcriptomic profiles of A. flavus at the whole-genome level. Differentially expressed genes (DEGs) associated with mycelial growth, conidia formation, and aflatoxin biosynthesis were revealed by comparing resveratrol-treated and untreated samples of A. flavus. This study would be meaningful for further annotating the genome of A. flavus, an understanding of aflatoxin biosynthesis and improving aflatoxin contamination management in agricultural products.

Results
Effect of resveratrol on the development and aflatoxin production of A. flavus Morphological characteristics and aflatoxin concentrations were quantified to define the effect of resveratrol on the fungal development and toxigenicity of A. flavus. The dry biomass of A. flavus mycelia grown in liquid A&M medium did not show significant difference between the resveratrol treatment (AM-Res) and the control (AM) ( Table 1). A similar result was obtained using solid A&M medium, in which resveratrol treatment did not significantly change the diameter of fungal mycelioid colonies ( Fig. 1 and Table 1). However, obvious variation in morphological characteristics of mycelia in response to resveratrol treatment in shaking liquid culture was observed (Fig. 2). The vegetative mycelial morphologic characteristics of A. flavus were different between the AM and the AM-Res cultures, and the vegetative mycelia of A. flavus appeared some smaller ones in the AM-Res. Resveratrol may affect the germination rate of conidia or growth rate of vegetative mycelia, leading to differences in the length of the mycelia. The results indicated that resveratrol did not affect the biomass accumulation of A. flavus, but it could influence the development of fungal mycelia.
Because no fungal conidium was detected in liquid culture, solid A&M medium was used to investigate the effect of resveratrol on conidia production by A. flavus. The number of conidia was 285.31 × 10 6 and 175.21 × 10 6 per 20 mL medium in the control and resveratrol treatment (3.0 μg/mL), respectively ( Fig. 1 and Table 1). In addition, the number of conidia was 197.83 × 10 6 , 180.83 × 10 6 , 163.06 × 10 6 , and 164.17 × 10 6 per 20 mL medium with resveratrol at concentrations of 0.3, 15.0, 30.0, and 60.0 μg/mL, respectively. Conidia formation of A. flavus was decreased in the resveratrol treatment compared with the control. The aflatoxin concentration was measured after 7 days of growth in liquid culture. The aflatoxin content in the resveratrol-treated culture (3.0 μg/mL) was 203.6 μg/L, which was much less than that in the control (Table 1), indicating that resveratrol inhibited aflatoxin formation. Our previous results also indicated that a resveratrol concentration of 3.0 μg/mL or higher inhibited aflatoxin production [12].

Defining the A. flavus transcriptome
Illumina sequencing of the resveratrol-treated and untreated A. flavus AF2202 generated 52.06 million and 52.72 million raw reads per library, respectively. Among total reads, approximately 51 million reads were purified from each group (Fig. 3), which included 46.00 million (89.42 %, AM-Res) and 45.34 million (89.68 %, AM) reads that were uniquely mapped to the A. flavus reference genome (Additional file 1). The data sets could represent the expressed sequences or transcriptomes of each library. Matching the reads to genes is important to annotate sequences because it could reveal the molecular mechanism behind gene expression [28]. Thus, transcriptomic data obtained in this experiment were used not only to define the differential gene expression in response to resveratrol but also to identify hitherto unknown transcripts including putative isoforms. The genic distribution of reads from mRNA-seq of the resveratrol-treated AF2202 (AM-Res) indicated that most reads (88.0 %) were mapped to protein-coding genes. The others were distributed between introns (0.7 %) and intergenic regions (11.3 %) (Fig. 4). The distribution of reads in the control AF2202 (AM) was similar to that of the resveratrol-treated (AM-Res), but fewer reads in the control were mapped to exons (Fig. 4). In total, 10.2 Gb of valid clean reads were yielded, and 98.82 % (9.13 Gb) of the 9.24 Gb of collected read pairs that passed filtering were mapped uniquely to the genome (Additional file 1), which was higher than the RNA-   seq data published previously for A. flavus [10,26,27] and A. oryzae [13].
In 2005, a 5× sequence coverage of the A. flavus genome was released by J. Craig Venter Research Institute, which has been recently updated at several web sites (http://www.aspergillusflavus.org/genomics/, NCBI). The A. flavus genome is about 37 Mb on 8 chromosomes, and it encodes over 13,487 functional genes (Aspergillus flavus Genome Sequencing Project), which previously were annotated as open reading frames in the genome without unambiguous transcript information. We quantified the overall transcriptional activity of genes in this study by calculating the number of reads per kilobase of exon per million mapped reads (RPKM) [29]. We found that 9,870 (71.14 %) of the 13,875 genes (http://fungi.ensembl.org/Aspergillus_flavus/Info/Annotation/#assembly) in the A. flavus genome database had expression with RPKM ≥ 1 in either of the two samples (Additional files 2 and 3). Besides the small fraction of reads that mapped to rRNA, approximately 11 % of reads could not be mapped to the encoding regions. This might be due to non-coding RNAs or novel exons of known genes. Novel genes were identified by assembling the RNA-seq data from the two samples using Cufflinks 2.1.1 [30]. Among the 14,592 annotated genes expressed in the two samples, 13,875 were previously annotated genes and 717 were novel genes (Additional file 4). According to the improved genome annotation, the transcriptional activity of the genes was estimated.
Functional classification of A. flavus DEGs in response to resveratrol The differently expressed genes between the two libraries provide a clue to the molecular mechanism related to the A. flavus response to resveratrol. By comparing expression levels in resveratrol-treated and untreated samples, we identified 453 genes that were significantly differentially transcribed (log 2 [fold change] = log 2 [AM-Res/AM] > 1, q value < 0.005) between the two samples, including 22 genes (4.86 % of the DEGs) showing a 10fold change (Additional file 5). Of these, 366 genes (80.79 %) were up-regulated and 87 (19.21 %) genes were down-regulated. Genes with altered expression were related to a wide variety of regulatory and metabolic processes.
Functional assignments were defined by Gene Ontology (GO) terms (http://www.genontology.org/), which provided a broad functional classification of genes and gene products for various biological processes, molecular functions, and cellular localizations. GO functional enrichment analysis revealed that these significantly DEGs were mainly involved in fatty acid metabolic process, monocarboxylic acid metabolic process, signaling receptor activity, phosphorelay sensor kinase activity, protein histidine kinase activity, and phosphotransferase activity (nitrogenous group as acceptor) (q value < 0.05) ( Table 2 and Additional file 6). As mentioned above, A. flavus and other Aspergilli preferentially colonize crop seeds with a high oil content because these crops contain high levels of unsaturated fatty acids, linoleic (18:2) and oleic acid (18:1), which are substrates for oxygenases [31,32]. Oxylipins, derived from oxygenases, are a class of oxygenated, unsaturated fatty acids involved in signaling pathways in filamentous fungi, yeasts, oomycetes, plants, and animals [33]. A family of oxylipin-producing oxygenases and their products, which is encoded by ppo (fungi) and lox (plants, animals, and fungi), are involved in regulating sclerotia and conidia formation and secondary metabolism in A. flavus, A. nidulans, and A. parasiticus [5,34]. Oxylipins could potentially interact with G-protein-coupled receptor complexes upstream of the heterotrimeric G-protein complex that regulates aflatoxin/sterigmatocystin production [35,36]. Moreover, the transcription of genes of fatty acid synthases (AFLA_139380 and AFLA_139370), the key genes in the aflatoxin biosynthetic pathway cluster, was significantly down-regulated when A. flavus was treated with resveratrol.
To further investigate the biological functions and interactions of genes, pathway-based analysis was conducted using the Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.genome.ad.jp/kegg) pathway database, which records the network of molecular interactions in the cells and variants specific to particular organisms. However, KEGG metabolic pathway analysis failed to confirm enrichment in any pathway in the RNA-seq data (q value < 0.05) obtained in this study, which is consistent with results of a previous study by Lin [10].

Validation of RNA-seq analysis by quantitative real-time PCR (qRT-PCR)
It is common practice to validate RNA-seq differential analysis with quantitative analysis of a subset of genes by qRT-PCR. Thus, we compared qRT-PCR expression levels of 16 DEGs randomly selected from the resveratroltreated and untreated A. flavus samples (Table 3). Genespecific primer pairs (Table 3) were designed according to The qRT-PCR data mirrored the quantification by RNAseq (Fig. 5), with a significant correlation (R 2 = 0.936), indicating the authenticity of these DEGs. Thus, these comparisons of data from qRT-PCR and RNA-seq of strain AF2202 not only validated the transcriptome analysis, but also reinforced the power of next-generation sequencing in quantifying the entire repertoire of genes.
Gene expression of the secondary metabolism gene clusters in response to resveratrol The A. flavus genome sequence contains hallmark enzymatic genes associated with secondary metabolite synthesis, including 55 putative secondary metabolism gene clusters [24,37,38]. Each secondary metabolism gene cluster contains enzymatic genes and transcription factors for compound synthesis [39]. To evaluate the effect of resveratrol on the regulation of secondary metabolism, we . The backbone genes included a non-ribosomal peptide synthase (NRPS, AFLA_004450), which is an enzyme that produces a polypeptide independent of ribosomes, and an NRPS-like synthase (AFLA_118440), which was named because it does not contain all of the typical domains normally found in an NRPS [2]. In the aflatoxin biosynthetic pathway cluster, 30 genes were identified as cluster #54 [2]. In this study, the transcription of 23 of the 30 genes was depressed to various degrees with resveratrol treatment, and only 7 genes (AFLA_139460, AFLA_139450, AFLA_139390, AFLA_139260, AFLA_139330, AFLA_139220, and AFLA_139210) were not down-regulated (Additional file 7). Interestingly, aflA (AFLA_139380) and aflB (AFLA_139370) were significantly down-regulated (q value < 0.005). aflA and aflB encode the fatty acid synthetase alpha subunit (Fas-2) and beta subunit (Fas-1), respectively, which catalyze the reaction between acetyl-CoA and malonyl-CoA to form the starter unit hexanoate of aflatoxin biosynthesis [40,41]. However, the transcription regulator genes aflR (AFLA_139360) and aflS (AFLA_139340) in the cluster showed only medium or slight down-regulation (Additional file 7), and the transcription of the global secondary metabolite regulators laeA (AFLA_033290) and veA (AFLA_066460) did not show significant depression (Table 4). Moreover, the expression of three genes (AFLA_139500, AFLA_139240, and AFLA_139230) was too low in the control sample to be distinguished by HTSeq (v0.5.3p9).
In cluster #2, only the backbone gene AFLA_004450, encoding a dimethylallyl transferase (DMAT, an indolealkaloids) enzyme, had significantly down-regulated transcription (q value < 0.005), but the transcription of the other 9 genes in this cluster was not significantly different when resveratrol was applied (Additional file 7). Also, in cluster #44, which has an unknown product, only the backbone gene AFLA_118440, encoding a polyketide synthetase (PKS) enzyme, was significantly up-regulated with resveratrol treatment (q value < 0.005) (Additional file 7). A total of 6 genes in cluster #9 (AFLA_017860 and AFLA_017820), #24 (AFLA_069370 and AFLA_069340), and #41 (AFLA_116250 and AFLA_116240) showed significant differential transcription (q value < 0.005) with resveratrol treatment. The products of these three clusters are unknown, but their backbone enzymes were predicted by SMURF. The backbone enzyme of cluster #9 is an NRPS-like enzyme, which is encoded by AFLA_017840. The backbone gene AFLA_069330 of cluster #24 encodes the NRPS Pes1, and the backbone enzyme of cluster #41 is a PKS. Genes in the other secondary metabolism gene clusters did not show significantly altered expression when resveratrol was applied. Only a few genes, such as AFLA_008690 in cluster #6, AFLA_054050 in cluster #16, AFLA_059950 in cluster #18, AFLA_060690 in cluster #19, AFLA_070820 in cluster #25, AFLA_084160 in cluster #30, and AFLA_108440 in cluster #37, were significantly affected (q value < 0.005) (Additional file 7).
Expression analysis of A. flavus growth-and conidial development-related genes in response to resveratrol Using the expression profiles of A. flavus in different treatments, especially the 453 genes that were significantly differentially transcribed (Additional file 5), we observed that some genes involved in mycelial and conidial development were down-regulated when resveratrol was applied (Table 4). Noticeably, the transcription of the socalled fluffy gene family (AFLA_101920, AFLA_137320, AFLA_134030, AFLA_131490, and AFLA_017380), whose null mutants have a fluffy or floccose phenotype [25], was down-regulated to various degrees (Table 4). fluG (AFLA_101920) and flbC (AFLA_137320) showed significantly depressed transcriptional levels (q value < 0.005) when A. flavus was treated with resveratrol. fluG encodes an extracellular development signal protein, which may function as a GSI-related enzyme in the synthesis of a small diffusible factor that acts as an extracellular signal directing asexual sporulation and other aspects of colony growth in A. nidulans [42]. flbC encodes a zinc finger protein (FlbC) that is required for morphogenesis and the elongation of flagellar filaments by facilitating polymerization of the flagellin monomers at the tip of a growing filament in Escherichia coli [43]. Concurrently, some conidia-specific genes (AFLA_098380 and AFLA_014260) were down-regulated slightly ( Table 4). The molecular action of conidial hydrophobins RodA/ RolA (AFLA_098380) and RodB/HypB (AFLA_014260) contributed to the structural integrity of the cell wall [44,45]. In addition, the transcription of genes related to conidial development, such as BrlA (AFLA_082850), AbaA (AFLA_029620), wetA (AFLA_052030), stuA (AFLA_046990), GpaA/FadA (AFLA_018340), and Medusa (AFLA_136410), was decreased to a certain extent (Table 4). Interestingly, transcription of the stuA gene (AFLA_046990), which encodes the APSES transcription factor StuA that affects the orderly differentiation and spatial organization of cell types of the conidiospore, was significantly decreased (the RPKM value decreased from 173.07 to 79.00, q value < 0.005) [46]. We also found that the transcription of the gene AFLA_099000 was significantly increased. In A. oryzae, AFLA_099000 encodes a superoxide dismutase [Cu-Zn] that destroys radicals within the cells and is toxic to biological systems (Table 4). Interestingly, the transcription of genes involved in vacuolar proteins, such as AFLA_027560, AFLA_017730, AFLA_039470, AFLA_060890, AFLA_130130, and Novel 100107, was significantly increased ( Table 4). The vacuolar proteins may take advantage of the exchange of atmosphere and nutrients between the organism and the environment. Nonribosomal peptide synthetase (NRPS) is a key factor responsible for the biosynthesis of bioactive metabolites that potentially contribute to organismal virulence; NRPS is encoded by AFLA_004450, which was significantly down-regulated (Table 4). RPKM: Reads per kilobase per million reads; q value: corrected p value; DEGs: differently expressed genes; AM-Res (Treatment) and AM (Control): A. flavus was cultured in A&M medium with or without resveratrol, respectively; log 2 (fold_change) = log 2 (AM-Res/AM); log 2 (fold_change) > 0: the expression of a gene was up-regulated; log 2 (fold_change) < 0: the expression of a gene was down-regulated

Discussion
Considerable progress has been achieved in reducing aflatoxin contamination in agricultural products [47][48][49][50][51][52] since the discovery that aflatoxin is extremely toxic to humans and animals [2,53]. Many inhibitors [9,54] have been demonstrated to have an inhibitory effect on aflatoxin biosynthesis in A. flavus. Among them, plantderived metabolites are of special interest. Non-host plant metabolites, such as onion, garlic extracts, eugenol [55][56][57], khellin, visnagin [58], caffeine, and piperlongumine [59] inhibit A. flavus growth and/or aflatoxin production. At the same time, host plant-derived metabolites including luteolin [60], eriodictyol [61], and tanning acids [62] also inhibit A. flavus development and/or aflatoxin production. However, most of these compounds can hardly be applied in practice because their biosynthesis pathway and the related biochemical steps are not well understood. Resveratrol, a natural phytoalexin, could protect plants from fungal infection and improve the defense capacity in plants [63]. Some studies indicated that phytoalexins play an important role in resistance to aflatoxin contamination in peanut under drought stress [64], and the isoflavone phytoalexin glyceollin from soybean could inhibit aflatoxin production [65]. Field experiments demonstrated that a correlation existed between increased phytoalexins and decreased aflatoxin B 1 [66]. The results from our previous experiments demonstrated that the resveratrol level was a factor that was relevant to aflatoxin resistance in peanut [12]. Resveratrol could affect A. flavus development and inhibit conidia formation and aflatoxin production. Moreover, the resveratrol biosynthetic pathway is well known [67]. Thus, it is possible to enhance resveratrol biosynthesis and concentration by breeding or genetic engineering, which may provide an applicable means for aflatoxin contamination management. Some transgenic studies indicated that resveratrol could inhibit other pathogenic fungi [68,69]. In addition, resveratrol has been suggested to have a beneficial effect on human health [70,71]. Genetic improvement for an increased resveratrol content in crops could not only ensure food safety, but also enhance the nutritional value of agricultural products.
Aflatoxins are polyketide-derived difurano cumarins, their biosynthesis involves at least 18 enzymatic reactions [72,73]. The first step in the biosynthesis requires the formation of the starter unit hexanoate from acetyl-CoA and malonyl-CoA catalyzed by Fas-1 and Fas-2 [40,41], which is then catalyzed by PksA, a polyketide synthase encoded by aflC, to produce norsolorinic acid anthrone (NAA). The monooxygenase HypC oxidizes the intermediate NAA to form anthraquinone norsolorinic acid, which is the first stable intermediate in aflatoxin biosynthesis [74,75]. In this study, we found that aflA (Fas-2) and aflB (Fas-1) in the aflatoxin biosynthesis pathway cluster were significantly down-regulated when A. flavus was treated with resveratrol. Aflatoxin is a highly oxygenated product, and oxidative stress is a prerequisite for aflatoxin production [76]. Resveratrolinhibited aflatoxin biosynthesis also has antioxidant activity. Superoxide dismutase [Cu-Zn] (AFLA_099000) ( Table 4) had a significantly increased expression level when resveratrol was added to medium; it destroyed radicals that were normally produced within the cells and were toxic to biological systems. Superoxide dismutases inhibit free radical scavenging and protect calcineurin (a Ca 2+ /calmodulin-dependent protein phosphatase) from inactivation [77]. The activity of calcineurin was enhanced at periods corresponding to aflatoxin production, and aflatoxin itself can inhibit calmodulin kinase activity [78]. Calmodulin also activated acetyl-CoA carboxylase, a key enzyme in aflatoxin biosynthesis [79]. The transcription of genes involved in vacuolar proteins significantly increased (Table 4); these proteins may be involved in the exchange of atmosphere and nutrients between the organism and the environment. Moreover, 23 of the 30 genes in the aflatoxin metabolic cluster had depressed transcription levels. These observations suggested that resveratrol directly inhibits the expression of the aflatoxin pathway cluster genes aflA and aflB, which could significantly suppress the formation of the starter unit (hexanoate) of aflatoxin biosynthesis. In addition, resveratrol highly increased the activity of antioxidative enzymes to limit aflatoxin production. There was also a certain amount of aflatoxin produced by A. flavus treated with resveratrol because aflatoxin biosynthesis was only partially inhibited rather than completely suppressed.
Based on our experimental data, resveratrol could influence the mycelia development and conidia formation of A. flavus (Figs. 1 and 2, Table 1), as well as several other biochemical substances previously reported [10,80,81]. Additionally, we have found that some genes involved in mycelial and conidial development were downregulated when A. flavus was treated with resveratrol ( Table 4). The changes in mycelia (the morphological characteristics of mycelia in liquid medium and the myceloid colony on solid medium) was affected by resveratrol, which may be related to the low-level expression of the so-called fluffy gene family. The formation of conidia, the asexual reproductive structure, in A. flavus requires the concerted activity of a number of signaling proteins and transcription factors [25]. Transcription of the conidia-specific genes RodA/RolA and RodB/HypB was down-regulated slightly when resveratrol was applied (Table 4), but the transcription of some genes related to conidial development showed various levels of depression when A. flavus was treated with resveratrol (Table 4). BrlA, AbaA, WetA, and StuA are pivotal regulators of conidiophore development and conidia maturation; they act individually and together to regulate their own expression and that of numerous other sporulation-related genes. The transcriptomic profiling showed that resveratrol significantly depressed the transcriptional levels of stuA, fluG, flbC, and other genes related to conidial development. BrlA, a Cys 2 His 2 -type conidiation transcription factor, mediates the developmental switch from the apical growth pattern of vegetative cells to the budding growth of conidiophores; brlA mutants cannot produce conidia [25]. The transcription factor AbaA controls the temporal and spatial specificity of development in Aspergillus. The expression of abaA leads to the activation of brlA and wetA and to the cessation of growth [82]. The regulatory protein WetA is responsible for activating a set of genes whose products make up the final two conidial wall layers or direct their assembly [82]. The APSES transcription factor StuA (stunted) is required for the orderly differentiation and spatial organization of cell types of the conidiophores. StuA affected the spatial localization of AbaA and BrlA, and the expression of StuA is dependent on a functional BrlA protein [46]. The transcript of GpaA/FadA was involved in conidial development. A previous study demonstrated that ΔgprA and ΔgprD mutants did not exhibit hyper-conidiation when exposed to 0.5 and 1.0 mg fatty acid-soaked discs [83]. The transcriptomic profiling showed that resveratrol significantly depressed the transcriptional levels of stuA, fluG, flbC, and other genes related to conidial development. These results demonstrated that resveratrol greatly reduced conidia formation of A. flavus by affecting gene expression at the transcriptomic level.
Numerous studies indicated that conidial development and secondary metabolism are co-regulated [84][85][86][87]. The heterotrimeric nuclear complex (VelB/VeA/LaeA complex, Velvet complex) was demonstrated to be a global regulator of fungal development and secondary metabolism in Aspergillus [10,84,85,88]. laeA and veA are two critical genes in the formation of the Velvet complex: VeA bridges VelB and LaeA to form the complex [2,88]. In our data, the transcription of laeA and veA were slightly depressed in the A. flavus samples treated with resveratrol (Table 4), which is consistent with the results of previous studies. Microarray data provided evidence for LaeA regulation of 24 clusters, including the aflatoxin biosynthetic cluster, out of 55 secondary metabolism clusters in A. flavus [87]. VeA regulates aflatoxin, aflatrem, and cyclopiazonic acid production [1]. In addition, null mutants of laeA and veA produced fewer conidia and less aflatoxin in host seeds [2,86]. This indicated that the inhibition of laeA and veA transcript levels by resveratrol may lead to a loss of toxin production and conidia formation.

Conclusions
Resveratrol treatment of A. flavus decreased aflatoxin production and conidia formation, and it also caused abnormal mycelia development. First, resveratrol directly inhibited the expression of aflatoxin biosynthetic pathway cluster genes, the aflA and aflB expression at significantly low transcriptional levels could result in an insufficient amount of the starter unit hexanoate for aflatoxin biosynthesis. In addition, resveratrol could significantly enhance the activity of antioxidative enzymes, which destroy radicals produced within the cells that are toxic to biological systems, leading to decreased aflatoxin production. Second, the decreased transcription of fluG and flbC could affect the mycelia (colony) development and asexual sporulation. The altered appearance of mycelia and mycelioid colonies could be a result of the depressed transcription of the so-called fluffy gene family. Third, resveratrol effectively decreased conidia formation of A. flavus, which may result from the suppressed transcription of genes related to conidial development. The low level of transcription of stuA, brlA, wetA, and abaA, the pivotal regulators of conidial development and maturation, as well as other genes related to conidial development, disrupted the cell's orderly differentiation and spatial organization of conidia, blocking conidia formation. Fourth, the inhibition of laeA and veA by resveratrol may lead to a loss of aflatoxin production and conidia formation, which further demonstrated that conidial development and secondary metabolism are co-regulated.

Effects of resveratrol treatment on A. flavus
A toxigenic A. flavus strain (AF2202) isolated from peanut was maintained at Oil Crops Research Institute of Chinese Academy of Agricultural Sciences in glycerol 20 % at −80°C. To prepare the A. flavus inoculum, some conidia of AF2202 were taken from the stored sample and cultured on fresh potato dextrose agar medium at 29 ± 1°C for 7 days. Conidia were then collected and suspended in sterile water containing 0.05 % Tween-80. The concentration of conidia in the suspension was determined using a haemocytometer.
Liquid A&M medium with a final resveratrol (Sigma®) concentration of 3.0 μg/mL [12] was prepared and added to conical 100-mL flasks containing 60 mL medium. Liquid A&M medium without resveratrol was used as a control. Each flask was inoculated with 5 × 10 5 colonyforming units (CFU) A. flavus conidia (freshly prepared). The inoculated flasks were incubated at 30°C with continuous shaking (160 rpm) for 3 or 7 days. After 3 days of incubation, the fungal mycelia were collected by filtration and stored at −80°C for RNA extraction. After 7 days of incubation, fungal mycelia were also collected by filtration, and the biomass was measured after drying overnight at 60°C. The dried mycelia were then ground to fine powder and put back into the corresponding flasks. Then, 25 mL chloroform was added to each flask and incubated with continuous shaking (200 rpm) for 30 min at room temperature. Five milliliters of the chloroform layer was added to a test tube and dried under nitrogen gas. The residue was dissolved in 1.0 mL methanol. Then, aflatoxins B 1 and B 2 were tested using high-performance liquid chromatography as previously described [12].
Solid A&M medium was prepared by adding agar to the above liquid medium with or without resveratrol (1.5 g agar was added to 100 mL liquid medium). The solid medium was used to investigate the effect of resveratrol on conidia formation and mycelia growth of A. flavus. A suspension with 1.6 × 10 5 conidia CFU was inoculated on a 90-mm diameter Petri dish containing 20 mL solid medium. The inoculated Petri dishes were incubated at 30°C in darkness for 7 days. For each dish, the conidia were suspended, the concentration was determined using a haemocytometer, and the colony diameter was measured.
Three replications were performed for each treatment in the above experiments.

RNA isolation and cDNA library construction
Aflatoxin was initially detected at 2 days of A. flavus (AF2202) growth in liquid A&M medium. The aflatoxin level increased at maximum rates between 3 and 4 days and then increased slowly after that point. However, the aflatoxin production rate in the resveratrol-treated culture (3.0 μg/mL) was much lower than that in the control when A. flavus was cultured between 3 and 4 days (Additional file 8). Additionally, the conidia were initially tested at 3 days of AF2202 growth on solid A&M medium. The morphological changes of A. flavus suggested that resveratrol could affect the expression of genes that are related to aflatoxin production and conidia formation. Therefore, we sampled the mycelium to isolate RNA and constructed cDNA libraries when AF2202 was cultured for 3 days.
Total RNA of AF2202 was isolated using the RNeasy® Plant Mini Kit (QIAGEN) according to the manufacturer's protocol. Each pooled RNA sample for cDNA library construction was composed of three individual RNA samples with equal weight extracted from mycelia that were independently cultured. All RNA samples were treated with RNase-free DNase I. The concentration and integrity of the pooled total RNA was checked using a NanoDrop® 1000 spectrophotometer, a Qubit® Fluorometer 2.0, and an Agilent 2100 bioanalyzer to confirm that all samples had an RNA integrity number value above 9.0. Sequencing libraries were constructed using the Illumina TruSeq™ RNA Sample Preparation Kit (Illumina) following the manufacturer's instructions. Four index codes were added to attribute sequences to each sample. The mRNA was isolated from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in Illumina proprietary fragmentation buffer. First-strand cDNA was synthesized using random oligonucleotides and SuperScript II. Second-strand cDNA synthesis was subsequently performed using DNA polymerase I and RNase H. Remaining overhangs were converted into blunt ends via exonuclease/polymerase, and then enzymes were inactivated. After adenylation of 3′ ends of cDNA fragments, Illumina PE adapter oligonucleotides were ligated for hybridization. The cDNA fragments of 200 bp in length were isolated with the AMPure XP system (Beckman Coulter). DNA fragments with ligated adapter molecules on both ends were selectively enriched using Illumina PCR Primer Cocktail in a 15-cycle PCR. Products were purified with the AMPure XP system and quantified using the Agilient high-sensitivity DNA assay with the Agilent Bioanalyzer 2100 system.
Clustering of the index-coded samples was performed with a cBot Cluster Generation System using a TruSeq PE Cluster Kit (v3-cBot-HS, Illumina) according to the manufacturer's instructions. After cluster generation, the prepared libraries were sequenced with the Illumina Hiseq 2000 platform, and then 100-bp paired-end reads were generated.
Mapping reads to the reference genome and quantifying gene expression Raw data (raw reads) of fastq format were first processed using in-house perl scripts. Clean data (clean reads) were obtained by removing reads containing adapter, reads containing poly-N, and low-quality reads from the raw data. The Q20, Q30, GC content, and sequence duplication level of the clean data were calculated (Additional 1 and Fig. 3). All downstream analysis used the clean data with high quality.

Identification and analysis of DEGs
Prior to differential gene expression analysis, the read counts for each sequenced library were adjusted using the edgeR program package (TMM) with one scaling normalization factor. Differential expression analysis of two conditions (with and without resveratrol treatment) was performed using DEGSeq (1.2.2) [91] and edgeR (3.2.4) [92]. The p values were adjusted with the Benjamini and Hochberg method [93], and a corrected p value (q value) of 0.005 and log 2 (fold change) of 1 were set as the threshold for significantly differential expression [94].
GO (http://www.geneontology.org/) is an international standardized gene function classification system that offers a dynamically updated, controlled vocabulary and strictly defined concepts to comprehensively describe properties of genes and their products in any organism. In our experiments, GO enrichment analysis of DEGs was implemented using GOseq (Release 2.12) [95], in which the gene length bias was corrected. GO terms with a corrected p value (q value) less than 0.05 were considered to be significantly enriched by DEGs.
KEGG is the major public pathway-related database [96] used to identify significantly enriched metabolic pathways or signal transduction pathways in DEGs compared to whole-genome background (http://www.genome.jp/kegg/). We used KOBAS software (v 2.0) [97] to test the statistical enrichment of DEGs in KEGG pathways. KEGG terms with a corrected p value (q value) less than 0.05 were considered to be significantly enriched in DEGs.

qRT-PCR analysis
To validate the repeatability and reproducibility of gene expression data obtained by RNA-seq in A. flavus, we randomly selected 16 DEGs for validation by qRT-PCR, as described previously [13,25,98,99]. Gene-specific primer pairs (Table 2) were designed according to the sequences of the 16 genes using the Primer 3-BLAST program available online (bioinfo.ut.ee/primer3-0.4.0/ primer3/, NCBI). To ensure accuracy, each primer was run with three replications on the same plate with a negative control that lacked template cDNA to detect non-specific products. For normalization, the A. flavus 18S rRNA gene was selected as the endogenous reference gene. The relative expression levels of genes were calculated using the 2 -ΔΔCt method [100,101], which represents the C T difference between the reference 18S rRNA and the target gene products. Three replications were performed for each qRT-PCR.