Deep sequencing analysis of transcriptomes in Aspergillus flavus in response to resveratrol
© Wang et al. 2015
Received: 9 January 2015
Accepted: 23 June 2015
Published: 16 September 2015
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.
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.
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.
The widespread fungus Aspergillus flavus is a saprophytic and opportunistic pathogen that is capable of causing diseases in agricultural crops , producing toxic secondary metabolites, and bringing about health problems in animals and humans . Under favorable conditions, A. flavus produces polyketide-derived carcinogenic and mutagenic secondary metabolites, aflatoxins , which are extremely hepatotoxic, immunosuppressive, and antinutritional mycotoxins associated with both acute and chronic toxicity in humans and animals . 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 aflatoxin-contaminated 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–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 . Our previous study found that resveratrol is highly related to resistance to aflatoxin production in peanut seeds and could inhibit aflatoxin production in medium . 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 cost-efficient high-throughput technology for transcriptomic profiling analysis that has been successfully used to interrogate transcriptomes of numerous fungi [10, 13–17]. The application of this technology has greatly accelerated our understanding of the complexity of gene expression, regulation, and networks . Gene expression can be more accurately quantified using RNA-seq approaches than conventional transcriptomic analysis . With higher sensitivity, RNA-seq can also be more efficient in detecting a larger range of dynamically expressed genes than microarrays. Furthermore, RNA-seq 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 . Moreover, the “reference-based” RNA-seq approach was proved to be more sensitive and have higher computational efficiency than a de novo assembly approach . The construction of expressed sequences tags in A. flavus was reported by several researchers , and the whole-genome sequencing of the fungus was completed . 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 , water activity , and temperature  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.
Effect of resveratrol on the development and aflatoxin production of A. flavus
Effect of resveratrol on the development and aflatoxin production of A. flavus
Mycelioid colony diameter
(1 × 106 CFU/20 mL)
0.411 ± 0.027
385.49 ± 17.38
7.1 ± 0.27
285.31 ± 52.50
0.464 ± 0.0067
203.55 ± 31.46 *
7.6 ± 0.35
175.21 ± 40.17 *
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 × 106 and 175.21 × 106 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 × 106, 180.83 × 106, 163.06 × 106, and 164.17 × 106 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 .
Defining the A. flavus transcriptome
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) . 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 . 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 (log2[fold change] = log2[AM-Res/AM] > 1, q value < 0.005) between the two samples, including 22 genes (4.86 % of the DEGs) showing a 10-fold 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.
Gene Ontology (GO) functional enrichment analysis of differentially expressed genes when A. flavus was treated with resveratrol
Fatty acid metabolic process
Monocarboxylic acid metabolic process
signaling receptor activity
Phosphorelay sensor kinase activity
Protein histidine kinase activity
Phosphotransferase activity (nitrogenous group as acceptor)
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 .
Validation of RNA-seq analysis by quantitative real-time PCR (qRT-PCR)
Primer pairs for each target gene for quantitative real-time PCR
Primer sequence (5′-3′)
Annealing temperature (°C)
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 . To evaluate the effect of resveratrol on the regulation of secondary metabolism, we used the information from the web-based software Secondary Metabolite Unknown Regions Finder (SMURF, http://www.jcvi.org/smurf) to identify the putative secondary metabolism gene clusters of A. flavus. Among the 55 clusters in the A. flavus genome, the expression levels of most clusters were not significantly affected by resveratrol (Additional file 7). By comparing the expression levels of secondary metabolic genes in resveratrol-treated and untreated samples, we identified 17 genes (2.45 % of 695 secondary metabolic genes) in 13 secondary metabolism gene clusters that were significantly differentially transcribed (log2[fold change] = log2[AM-Res/AM] > 1, q value < 0.005), with 2 backbone genes (AFLA_004450 and AFLA_118440) (Additional file 7). 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 .
Differentially expressed genes related to the development and conidial formation of A. flavus
Annotated gene function
Extracellular developmental signal biosynthesis protein FluG
Developmental regulator FlbA development.
C2H2 conidiation transcription factor FlbC
Conidial hydrophobin RodA/RolA
Conidial hydrophobin RodB/HypB
G-protein complex alpha subunit GpaA/FadA
APSES transcription factor StuA (stunted)
Transcriptional regulator Medusa(medusa)
C2H2 type conidiation transcription factor BrlA
Transcription factor AbaA
Regulatory protein wetA
Cell pattern formation-associated protein stuA
Regulator of secondary metabolism LaeA
Developmental regulator AflYf / VeA
Nucleoside diphosphatase Gda1/VelB
Vacuolar protein sorting-associated protein 13
Regulator of V-ATPase in vacuolar membrane protein 1
Vacuolar calcium ion transporter
Vacuolar amino acid transporter 3
Putative vacuolar protein sorting-associated protein TDA6
Vacuolar protein sorting-associated protein 13
Superoxide dismutase [Cu-Zn]
Nonribosomal peptide synthetase 12
In cluster #2, only the backbone gene AFLA_004450, encoding a dimethylallyl transferase (DMAT, an indole-alkaloids) 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 so-called fluffy gene family (AFLA_101920, AFLA_137320, AFLA_134030, AFLA_131490, and AFLA_017380), whose null mutants have a fluffy or floccose phenotype , 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 . 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 . 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) .
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).
Considerable progress has been achieved in reducing aflatoxin contamination in agricultural products [47–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, plant-derived metabolites are of special interest. Non-host plant metabolites, such as onion, garlic extracts, eugenol [55–57], khellin, visnagin , caffeine, and piperlongumine  inhibit A. flavus growth and/or aflatoxin production. At the same time, host plant-derived metabolites including luteolin , eriodictyol , and tanning acids  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 . Some studies indicated that phytoalexins play an important role in resistance to aflatoxin contamination in peanut under drought stress , and the isoflavone phytoalexin glyceollin from soybean could inhibit aflatoxin production . Field experiments demonstrated that a correlation existed between increased phytoalexins and decreased aflatoxin B1 . The results from our previous experiments demonstrated that the resveratrol level was a factor that was relevant to aflatoxin resistance in peanut . Resveratrol could affect A. flavus development and inhibit conidia formation and aflatoxin production. Moreover, the resveratrol biosynthetic pathway is well known . 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 . Resveratrol-inhibited 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 Ca2+/calmodulin-dependent protein phosphatase) from inactivation . The activity of calcineurin was enhanced at periods corresponding to aflatoxin production, and aflatoxin itself can inhibit calmodulin kinase activity . Calmodulin also activated acetyl-CoA carboxylase, a key enzyme in aflatoxin biosynthesis . 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 down-regulated 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 . 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 Cys2His2-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 . 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 . 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 . 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 . 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 . 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–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 . VeA regulates aflatoxin, aflatrem, and cyclopiazonic acid production . 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.
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  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 × 105 colony-forming 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 B1 and B2 were tested using high-performance liquid chromatography as previously described .
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 × 105 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.
The A. flavus (NRRL 3357) genome and gene model annotation files were downloaded from the genome website (ftp://ftp.ensemblgenomes.org/pub/release-17/fungi/fasta/Aspergillus_flavus/dna/) directly. The index of the reference genome was built using Bowtie (v2.0.6) , and paired-end clean reads were aligned to the reference genome using Tophat (v2.0.7) with “mismatch 2” as the parameter . HTSeq (v0.5.3p9) (http://www-huber.embl.de/users/anders/HTSeq) was used to count the read numbers mapped to each gene. Then, the RPKM of each gene was calculated based on the length of the gene and read count mapped to this gene .
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)  and edgeR (3.2.4) . The p values were adjusted with the Benjamini and Hochberg method , and a corrected p value (q value) of 0.005 and log2 (fold change) of 1 were set as the threshold for significantly differential expression .
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) , 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  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)  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.
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 CT difference between the reference 18S rRNA and the target gene products. Three replications were performed for each qRT-PCR.
Availability of supporting data
The sequencing data generated in this study have been deposited in NCBI’s Short Read Archive database (SRA, http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/Traces/sra_sub/sub.cgi) and are accessible through SRA series accession number SRP057550 (BioProject ID: PRJNA281697).
This research was funded by the National Basic Research Program of China (2013CB127803), the National Natural Science Foundation of China (No. 31371662 and 31461143022), and the China Agriculture Research System (No. CARS-14).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Duran RM, Cary JW, Calvo AM. Production of cyclopiazonic acid, aflatrem, and aflatoxin by Aspergillus flavus is regulated by veA, a gene necessary for sclerotial formation. Appl Microbiol Biotechnol. 2007;73(5):1158–68.PubMedView ArticleGoogle Scholar
- Amaike S, Keller NP. Aspergillus flavus. Annu Rev Phytopathol. 2011;49:107–33.PubMedView ArticleGoogle Scholar
- Klich MA. Aspergillus flavus: the major producer of aflatoxin. Mol Plant Pathol. 2007;8(6):713–22.PubMedView ArticleGoogle Scholar
- Williams JH, Phillips TD, Jolly PE, Stiles JK, Jolly CM, Aggarwal D. Human aflatoxicosis in developing countries: a review of toxicology, exposure, potential health consequences, and interventions. Am J Clin Nutrit. 2004;80(5):1106–22.PubMedGoogle Scholar
- Burow GB, Nesbitt TC, Dunlap J, Keller NP. Seed lipoxygenase products modulate Aspergillus mycotoxin biosynthesis. Mol Plant Microbe Interact. 1997;10(3):380–7.View ArticleGoogle Scholar
- Goodrichtanrikulu M, Mahoney NE, Rodriguez SB. The plant growth regulator methyl jasmonate inhibits aflatoxin production by Aspergillus flavus. Microbiol Uk. 1995;141:2831–7.View ArticleGoogle Scholar
- Chen ZY, Brown RL, Damann KE, Cleveland TE. Identification of maize kernel endosperm proteins associated with resistance to aflatoxin contamination by Aspergillus flavus. Phytopathology. 2007;97(9):1094–103.PubMedView ArticleGoogle Scholar
- Wang T, Zhang E, Chen X, Li L, Liang X. Identification of seed proteins associated with resistance to pre-harvested aflatoxin contamination in peanut (Arachis hypogaea L). BMC Plant Biol. 2010;10:267–77.PubMed CentralPubMedView ArticleGoogle Scholar
- Holmes RA, Boston RS, Payne GA. Diverse inhibitors of aflatoxin biosynthesis. Appl Microbiol Biotechnol. 2008;78(4):559–72.PubMedView ArticleGoogle Scholar
- Lin J-Q, Zhao X-X, Zhi Q-Q, Zhao M, He Z-M. Transcriptomic profiling of Aspergillus flavus in response to 5-azacytidine. Fungal Genet Biol. 2013;56:78–86.PubMedView ArticleGoogle Scholar
- Chung IM, Park MR, Chun JC, Yun SJ. Resveratrol accumulation and resveratrol synthase gene expression in response to abiotic stresses and hormones in peanut plants. Plant Sci. 2003;164(1):103–9.View ArticleGoogle Scholar
- Wang H, Huang J, Lei Y, Yan L, Wang S, Jiang H, et al. Relationship of resveratrol content and resistance to aflatoxin accumulation caused by Aspergillus flavus in peanut seeds. Acta Agronomica Sinica. 2012;38(10):1875–83.View ArticleGoogle Scholar
- Wang B, Guo G, Wang C, Lin Y, Wang X, Zhao M, et al. Survey of the transcriptome of Aspergillus oryzae via massively parallel mRNA sequencing. Nucleic Acids Res. 2010;38(15):5075–87.PubMed CentralPubMedView ArticleGoogle Scholar
- Nagalakshmi U, Wang Z, Waern K, Shou C, Raha D, Gerstein M, et al. The transcriptional landscape of the yeast genome defined by RNA sequencing. Science. 2008;320(5881):1344–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Wilhelm BT, Marguerat S, Watt S, Schubert F, Wood V, Goodhead I, et al. Dynamic repertoire of a eukaryotic transcriptome surveyed at single-nucleotide resolution. Nature. 2008;453(7199):1239–43.PubMedView ArticleGoogle Scholar
- Linz JE, Wee J, Roze LV. Aspergillus parasiticus SU-1 Genome Sequence, Predicted Chromosome Structure, and Comparative Gene Expression under Aflatoxin-Inducing Conditions: Evidence that Differential Expression Contributes to Species Phenotype. Eukaryot Cell. 2014;13(8):1113–23.PubMed CentralPubMedView ArticleGoogle Scholar
- Chang PK, Scharfenstein LL, Mack B, Yu JJ, Ehrlich KC. Transcriptomic profiles of Aspergillus flavus CA42, a strain that produces small sclerotia, by decanal treatment and after recovery. Fungal Genet Biol. 2014;68:39–47.PubMedView ArticleGoogle Scholar
- McGettigan PA. Transcriptomics in the RNA-seq era. Curr Opin Chem Biol. 2013;17(1):4–11.PubMedView ArticleGoogle Scholar
- Maher CA, Kumar-Sinha C, Cao X, Kalyana-Sundaram S, Han B, Jing X, et al. Transcriptome sequencing to detect gene fusions in cancer. Nature. 2009;458(7234):97–101.PubMed CentralPubMedView ArticleGoogle Scholar
- Martin JA, Wang Z. Next-generation transcriptome assembly. Nat Rev Genet. 2011;12(10):671–82.PubMedView ArticleGoogle Scholar
- Payne GA, Nierman WC, Wortman JR, Pritchard BL, Brown D, Dean RA, et al. Whole genome comparison of Aspergillus flavus and A-oryzae. Med Mycol. 2006;44:S9–S11.View ArticleGoogle Scholar
- Chang P-K, Ehrlich KC. What does genetic diversity of Aspergillus flavus tell us about Aspergillus oryzae? Int J Food Microbiol. 2010;138(3):189–99.PubMedView ArticleGoogle Scholar
- Keller NP, Turner G, Bennett JW. Fungal secondary metabolism - From biochemistry to genomics. Nat Rev Microbiol. 2005;3(12):937–47.PubMedView ArticleGoogle Scholar
- Khaldi N, Seifuddin FT, Turner G, Haft D, Nierman WC, Wolfe KH, et al. SMURF: Genomic mapping of fungal secondary metabolite clusters. Fungal Genet Biol. 2010;47(9):736–41.PubMed CentralPubMedView ArticleGoogle Scholar
- Wilkinson JR, Kale SP, Bhatnagar D, Yu J, Ehrlich KC. Expression Profiling of Non-Aflatoxigenic Aspergillus parasiticus Mutants Obtained by 5-Azacytosine Treatment or Serial Mycelial Transfer. Toxins. 2011;3(8):932–48.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhang F, Guo Z, Zhong H, Wang S, Yang W, Liu Y, et al. RNA-Seq-Based Transcriptome Analysis of Aflatoxigenic Aspergillus flavus in Response to Water Activity. Toxins. 2014;6(11):3187–207.PubMed CentralPubMedView ArticleGoogle Scholar
- Yu J, Fedorova ND, Montalbano BG, Bhatnagar D, Cleveland TE, Bennett JW, et al. Tight control of mycotoxin biosynthesis gene expression in Aspergillus flavus by temperature as revealed by RNA-Seq. Fems Microbiol Lett. 2011;322(2):145–9.PubMedView ArticleGoogle Scholar
- Hegedus Z, Zakrzewska A, Agoston VC, Ordas A, Racz P, Mink M, et al. Deep sequencing of the zebrafish transcriptome response to mycobacterium infection. Mol Immunol. 2009;46(15):2918–30.PubMedView ArticleGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8.PubMedView ArticleGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.PubMed CentralPubMedView ArticleGoogle Scholar
- Calvo AM, Hinze LL, Gardner HW, Keller NP. Sporogenic effect of polyunsaturated fatty acids on development of Aspergillus spp. Appl Environ Microbiol. 1999;65(8):3668–73.PubMed CentralPubMedGoogle Scholar
- Feussner I, Wasternack C. The lipoxygenase pathway. Annu Rev Plant Biol. 2002;53:275–97.PubMedView ArticleGoogle Scholar
- Funk CD. Prostaglandins and leukotrienes: Advances in eicosanoid biology. Science. 2001;294(5548):1871–5.PubMedView ArticleGoogle Scholar
- Brown SH, Zarnowski R, Sharpee WC, Keller NP. Morphological transitions governed by density dependence and lipoxygenase activity in Aspergillus flavus. Appl Environ Microbiol. 2008;74(18):5674–85.View ArticleGoogle Scholar
- Hicks JK, Yu JH, Keller NP, Adams TH. Aspergillus sporulation and mycotoxin production both require inactivation of the FadA G alpha protein-dependent signaling pathway. Embo J. 1997;16(16):4916–23.PubMed CentralPubMedView ArticleGoogle Scholar
- Brodhagen M, Keller NP. Signalling pathways connecting mycotoxin production and sporulation. Mol Plant Pathol. 2006;7(4):285–301.PubMedView ArticleGoogle Scholar
- Cleveland TE, Yu J, Fedorova N, Bhatnagar D, Payne GA, Nierman WC, et al. Potential of Aspergillus flavus genomics for applications in biotechnology. Trends Biotechnol. 2009;27(3):151–7.PubMedView ArticleGoogle Scholar
- Georgianna DR, Fedorova ND, Burroughs JL, Dolezal AL, Bok JW, Horowitz-Brown S, et al. Beyond aflatoxin: four distinct expression patterns and functional roles associated with Aspergillus flavus secondary metabolism gene clusters. Mol Plant Pathol. 2010;11(2):213–26.PubMed CentralPubMedView ArticleGoogle Scholar
- Hoffmeister D, Keller NP. Natural products of filamentous fungi: enzymes, genes, and their regulation. Nat Prod Rep. 2007;24(2):393–416.PubMedView ArticleGoogle Scholar
- Mahanti N, Bhatnagar D, Cary JW, Joubran J, Linz JE. Structure and function of fas-1A, a gene encoding a putative fatty acid synthetase directly involved in aflatoxin biosynthesis in Aspergillus parasiticus. Appl Environ Microbiol. 1996;62(1):191–5.PubMed CentralPubMedGoogle Scholar
- Minto RE, Townsend CA. Enzymology and molecular biology of aflatoxin biosynthesis. Chem Rev. 1997;97(7):2537–55.PubMedView ArticleGoogle Scholar
- Lee BN, Adams TH. The Aspergillus nidulans fluG gene is required for production of an extracellular developmental signal and is related to prokaryotic glutamine synthetase I. Genes Dev. 1994;8(6):641–51.PubMedView ArticleGoogle Scholar
- Kuwajima G, Asaka J, Fujiwara T, Node K, Kondo E. Nucleotide sequence of the hag gene encoding flagellin of Escherichia coli. J Bacteriol. 1986;168(3):1479–83.PubMed CentralPubMedGoogle Scholar
- Stringer MA, Dean RA, Sewall TC, Timberlake WE. Rodletless, a new Aspergillus developmental mutant induced by directed gene inactivation. Genes Dev. 1991;5(7):1161–71.PubMedView ArticleGoogle Scholar
- Fedorova ND, Khaldi N, Joardar VS, Maiti R, Amedeo P, Anderson MJ, et al. Genomic islands in the pathogenic filamentous fungus Aspergillus fumigatus. Plos Genetics. 2008;4(4), e1000046.PubMed CentralPubMedView ArticleGoogle Scholar
- Miller KY, Wu J, Miller BL. StuA is required for cell pattern formation in Aspergillus. Genes Dev. 1992;6(9):1770–82.PubMedView ArticleGoogle Scholar
- Weaver MA, Abbas HK, Falconer LL, Allen TW, Pringle HC, Sciumbato GL. Biological control of aflatoxin is effective and economical in Mississippi field trials. Crop Prot. 2015;69:52–5.View ArticleGoogle Scholar
- Waliyar F, Osiru M, Ntare BR, Kumar KVK, Sudini H, Traore A, et al. Post-harvest management of aflatoxin contamination in groundnut. World Mycotoxin J. 2015;8(2):245–52.View ArticleGoogle Scholar
- Mauro A, Battilani P, Cotty PJ. Atoxigenic Aspergillus flavus endemic to Italy for biocontrol of aflatoxins in maize. Biocontrol. 2015;60(1):125–34.View ArticleGoogle Scholar
- Ehrlich KC, Moore GG, Mellon JE, Bhatnagar D. Challenges facing the biological control strategy for eliminating aflatoxin contamination. World Mycotoxin J. 2015;8(2):225–33.View ArticleGoogle Scholar
- Chen ZY, Rajasekaran K, Brown RL, Sayler RJ, Bhatnagar D. Discovery and confirmation of genes/proteins associated with maize aflatoxin resistance. World Mycotoxin J. 2015;8(2):211–24.View ArticleGoogle Scholar
- Torres AM, Barros GG, Palacios SA, Chulze SN, Battilani P. Review on pre- and post-harvest management of peanuts to minimize aflatoxin contamination. Food Res Int. 2014;62:11–9.View ArticleGoogle Scholar
- Forgacs J, Carll WT. Mycotoxicoses. Adv Vet Sci. 1962;7:273–882.Google Scholar
- Ma JY, Mo HZ, Chen Y, Ding D, Hu LB. Inhibition of Aflatoxin Synthesis in Aspergillus flavus by Three Structurally Modified Lentinans. Int J Mol Sci. 2014;15(3):3860–70.PubMed CentralPubMedView ArticleGoogle Scholar
- Bilgrami KS, Sinha KK, Sinha AK. Inhibition of aflatoxin production and growth of Aspergillus flavus bt eugenol and garlic extracts. Indian J Med Res Section B. 1992;96:171–5.Google Scholar
- Ehrlich K, Ciegler A. Effect of phytate on aflatoxin formation by Aspergillus parasiticus grown on different grains. Mycopathologia. 1985;92(1):3–6.PubMedView ArticleGoogle Scholar
- Mahoney N, Molyneux RJ. Phytochemical inhibition of aflatoxigenicity in Aspergillus flavus by constituents of walnut (Juglans regia). J Agric Food Chem. 2004;52(7):1882–9.PubMedView ArticleGoogle Scholar
- Mabrouk SS, El-Shayeb NMA. Inhibition of aflatoxin production in Aspergillus flavus natural coumarins chromones. World J Microbiol Biotechnol. 1992;8(1):60–2.PubMedView ArticleGoogle Scholar
- Lee SE, Mahoney NE, Campbell BC. Inhibition of aflatoxin B-1 biosynthesis by piperlongumine isolated from Piper longum L. J Microbiol Biotechnol. 2002;12(4):679–82.Google Scholar
- Norton RA. Inhibition of aflatoxin B-1 biosynthesis in Aspergillus flavus by anthocyanidins and related flavonoids. J Agric Food Chem. 1999;47(3):1230–5.PubMedView ArticleGoogle Scholar
- DeLucca II AJ, Palmgren MS, Daigle DJ. Depression of aflatoxin production by flavonoid-type compounds from peanut shells. Phytopathology. 1987;77(11):1560–3.View ArticleGoogle Scholar
- Azaizeh HA, Pettit RE, Sarr BA, Phillips TD. Effect of peanut tannin extracts on growth of Aspergillus parasiticus and aflatoxin production. Mycopathologia. 1990;110(3):125–32.PubMedView ArticleGoogle Scholar
- Hain R, Bieseler B, Kindl H, Schroeder G, Stoecker R. Expression of a stilbene synthase gene in nicotiana-tabacum results in synthesis of the phytoalexin resveratrol. Plant Mol Biol. 1990;15(2):325–36.PubMedView ArticleGoogle Scholar
- Dorner JW, Cole RJ, Sanders TH, Blankenship PD. Interrelationship of kernel water activity, soil temperature, maturity, and phytoalexin production in preharvest aflatoxin contamination of drought-stressed peanuts. Mycopathologia. 1989;105(2):117–28.PubMedView ArticleGoogle Scholar
- Song DK, Karr AL. Soybean phytoalexin, glyceollin, prevents accumulation of aflatoxin B-1 in cultures of Aspergillus flavus. J Chem Ecol. 1993;19(6):1183–94.PubMedView ArticleGoogle Scholar
- Zeringue Jr HJ. Effects of methyl jasmonate on phytoalexin production and aflatoxin control in the developing cotton boll. Biochem Syst Ecol. 2002;30(6):497–503.View ArticleGoogle Scholar
- Halls C, Yu O. Potential for metabolic engineering of resveratrol biosynthesis. Trends Biotechnol. 2008;26(2):77–81.PubMedView ArticleGoogle Scholar
- Fettig S, Hess D. Expression of a chimeric stilbene synthase gene in transgenic wheat lines. Transgenic Res. 1999;8(3):179–89.View ArticleGoogle Scholar
- Thomzik JE, Stenzel K, Stocker R, Schreier PH, Hain R, Stahl DJ. Synthesis of a grapevine phytoalexin in transgenic tomatoes (Lycopersicon esculentum Mill.) conditions resistance against Phytophthora infestans. Physiol Mol Plant Pathol. 1997;51(4):265–78.View ArticleGoogle Scholar
- Athar M, Back JH, Kopelovich L, Bickers DR, Kim AL. Multiple molecular targets of resveratrol: Anti-carcinogenic mechanisms. Arch Biochem Biophys. 2009;486(2):95–102.PubMed CentralPubMedView ArticleGoogle Scholar
- Tili E, Michaille J-J. Resveratrol, microRNAs, inflammation, and cancer. J Nucleic Acids. 2011;2011:102431–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Yabe K, Nakajima H. Enzyme reactions and genes in aflatoxin biosynthesis. Appl Microbiol Biotechnol. 2004;64(6):745–55.PubMedView ArticleGoogle Scholar
- Abrar M, Anjum FM, Butt MS, Pasha I, Randhawa MA, Saeed F, et al. Aflatoxins: biosynthesis, occurrence, toxicity, and remedies. Crit Rev Food Sci Nutr. 2013;53(8):862–74.PubMedView ArticleGoogle Scholar
- Ehrlich KC. Predicted Roles of the Uncharacterized Clustered Genes in Aflatoxin Biosynthesis. Toxins. 2009;1(1):37–58.PubMed CentralPubMedView ArticleGoogle Scholar
- Ehrlich KC, Li P, Scharfenstein L, Chang P-K. HypC, the Anthrone Oxidase Involved in Aflatoxin Biosynthesis. Appl Environ Microbiol. 2010;76(10):3374–7.PubMed CentralPubMedView ArticleGoogle Scholar
- Jayashree T, Subramanyam C. Oxidative stress as a prerequisite for aflatoxin production by Aspergillus parasiticus. Free Radical Biol Med. 2000;29(10):981–5.View ArticleGoogle Scholar
- Wang XT, Culotta VC, Klee CB. Superoxide dismutase protects calcineurin from inactivation. Nature. 1996;383(6599):434–7.PubMedView ArticleGoogle Scholar
- Jayashree T, Rao JP, Subramanyam C. Regulation of aflatoxin production by Ca2+/calmodulin-dependent protein phosphorylation and dephosphorylation. Fems Microbiol Lett. 2000;183(2):215–9.PubMedView ArticleGoogle Scholar
- Rao JP, Subramanyam C. Calmodulin mediated activation of acetyl-CoA carboxylase during aflatoxin production by Aspergillus parasiticus. Lett Appl Microbiol. 2000;30(4):277–81.PubMedView ArticleGoogle Scholar
- Chalfoun SM, Pereira MC, Resende MLV, Angelico CL, da Silva RA. Effect of powdered spice treatments on mycelial growth, sporulation and production of aflatoxins by toxigenic fungi. Ciencia E Agrotecnologia. 2004;28(4):856–62.View ArticleGoogle Scholar
- Greene-McDowelle DM, Ingber B, Wright MS, Zeringue HJ, Bhatnagar D, Cleveland TE. The effects of selected cotton-leaf volatiles on growth, development and aflatoxin production of Aspergillus parasiticus. Toxicon. 1999;37(6):883–93.PubMedView ArticleGoogle Scholar
- Mirabito PM, Adams TH, Timberlake WE. Interactions of three sequentially expressed genes control temporal and spatial specificity in Aspergillus development. Cell. 1989;57(5):859–68.PubMedView ArticleGoogle Scholar
- Affeldt KJ, Brodhagen M, Keller NP. Aspergillus Oxylipin Signaling and Quorum Sensing Pathways Depend on G Protein-Coupled Receptors. Toxins. 2012;4(9):695–717.PubMed CentralPubMedView ArticleGoogle Scholar
- Bok JW, Keller NP. LaeA, a regulator of secondary metabolism in Aspergillus spp. Eukaryot Cell. 2004;3(2):527–35.PubMed CentralPubMedView ArticleGoogle Scholar
- Bok JW, Noordermeer D, Kale SP, Keller NP. Secondary metabolic gene cluster silencing in Aspergillus nidulans. Mol Microbiol. 2006;61(6):1636–45.PubMedView ArticleGoogle Scholar
- Duran RM, Cary JW, Calvo AM. The role of veA in Aspergillus flavus infection of peanut, corn and cotton. Open Mycol J. 2009;3:27–36.View ArticleGoogle Scholar
- Kale SP, Milde L, Trapp MK, Frisvad JC, Keller NP, Bok JW. Requirement of LaeA for secondary metabolism and sclerotial production in Aspergillus flavus. Fungal Genet Biol. 2008;45(10):1422–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Bayram O, Krappmann S, Ni M, Bok JW, Helmstaedt K, Valerius O, et al. VelB/VeA/LaeA complex coordinates light signal with fungal development and secondary metabolism. Science. 2008;320(5882):1504–6.PubMedView ArticleGoogle Scholar
- Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26(1):136–8.PubMedView ArticleGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK. EdgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.PubMed CentralPubMedView ArticleGoogle Scholar
- Benjamini Y, Hochberg Y. Controlling the false discovery rate - a practical and powerful approach to multiple testing. J Royal Stat Soc Series B-Methodological. 1995;57(1):289–300.Google Scholar
- Fang S-M, Hu B-L, Zhou Q-Z, Yu Q-Y, Zhang Z. Comparative analysis of the silk gland transcriptomes between the domestic and wild silkworms. BMC Genomics. 2015;16:60–71.PubMed CentralPubMedView ArticleGoogle Scholar
- Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2):R14–25.PubMed CentralPubMedView ArticleGoogle Scholar
- Saldanha AJ. Java Treeview-extensible visualization of microarray data. Bioinformatics. 2004;20(17):3246–8.PubMedView ArticleGoogle Scholar
- Mao XZ, Cai T, Olyarchuk JG, Wei LP. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics. 2005;21(19):3787–93.PubMedView ArticleGoogle Scholar
- Reverberi M, Punelli M, Scala V, Scarpari M, Uva P, Mentzen WI, et al. Genotypic and Phenotypic Versatility of Aspergillus flavus during Maize Exploitation. Plos One. 2013;8(7), e68735.PubMed CentralPubMedView ArticleGoogle Scholar
- Song Z, Yin Y, Jiang S, Liu J, Chen H, Wang Z. Comparative transcriptome analysis of microsclerotia development in Nomuraea rileyi. BMC Genomics. 2013;14:411–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(−Delta Delta C) method. Methods. 2001;25(4):402–8.PubMedView ArticleGoogle Scholar
- Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29(9), e45.PubMed CentralPubMedView ArticleGoogle Scholar