Next Article in Journal
Indoxyl Sulfate, a Tubular Toxin, Contributes to the Development of Chronic Kidney Disease
Next Article in Special Issue
Assessment of Citrinin in Spices and Infant Cereals Using Immunoaffinity Column Clean-Up with HPLC-Fluorescence Detection
Previous Article in Journal
Transcriptomic and Proteomic Analysis Reveals Mechanisms of Patulin-Induced Cell Toxicity in Human Embryonic Kidney Cells
Previous Article in Special Issue
Transcriptomic Insights into the Antifungal Effects of Magnolol on the Growth and Mycotoxin Production of Alternaria alternata
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dimethylformamide Inhibits Fungal Growth and Aflatoxin B1 Biosynthesis in Aspergillus flavus by Down-Regulating Glucose Metabolism and Amino Acid Biosynthesis

1
Institute of Food Science and Technology, Chinese Academy of Agricultural Sciences/Key Laboratory of Agro-products Quality and Safety Control in Storage and Transport Process, Ministry of Agriculture, Yuanmingyuan West Road, Haidian District, Beijing 100193, China
2
College of Food Science and Engineering, Qingdao Agricultural University, Qingdao 266109, China
*
Author to whom correspondence should be addressed.
Toxins 2020, 12(11), 683; https://doi.org/10.3390/toxins12110683
Submission received: 5 October 2020 / Revised: 22 October 2020 / Accepted: 27 October 2020 / Published: 29 October 2020
(This article belongs to the Special Issue Detection and Prevention Technologies for Toxins)

Abstract

:
Aflatoxins (AFs) are secondary metabolites produced by plant fungal pathogens infecting crops with strong carcinogenic and mutagenic properties. Dimethylformamide (DMF) is an excellent solvent widely used in biology, medicine and other fields. However, the effect and mechanism of DMF as a common organic solvent against fungal growth and AFs production are not clear. Here, we discovered that DMF had obvious inhibitory effect against A. flavus, as well as displayed complete strong capacity to combat AFs production. Hereafter, the inhibition mechanism of DMF act on AFs production was revealed by the transcriptional expression analysis of genes referred to AFs biosynthesis. With 1% DMF treatment, two positive regulatory genes of AFs biosynthetic pathway aflS and aflR were down-regulated, leading to the suppression of the structural genes in AFs cluster like aflW, aflP. These changes may be due to the suppression of VeA and the subsequent up-regulation of FluG. Exposure to DMF caused the damage of cell wall and the dysfunction of mitochondria. In particular, it is worth noting that most amino acid biosynthesis and glucose metabolism pathway were down-regulated by 1% DMF using Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Taken together, these RNA-Seq data strongly suggest that DMF inhibits fungal growth and aflatoxin B1 (AFB1) production by A. flavus via the synergistic interference of glucose metabolism, amino acid biosynthesis and oxidative phosphorylation.
Key Contribution: The inhibitory mechanism of action of DMF on fungal growth and AFs biosynthesis at the transcriptomic level is elucidated. In terms of agricultural applications; this research may provide a basis for synergistic antifungal and antitoxigenic effect between DMF and other fungicides.

1. Introduction

Aspergillus flavus as plant-invasive fungal pathogens cause enormous losses in the yield and quality of field crops worldwide [1]. Under the suitable environmental conditions, A. flavus is prone to produce a series of strong carcinogenic and mutagenic secondary metabolites aflatoxins (AFs) during the process of infecting food and feed [2]. AFs are the predominant and most carcinogenic naturally occurring compounds which inevitably result in health complications, including hepatocellular carcinoma, acute intoxication, immune system disorder and growth retardation in children [3,4]. In 1993, AFs were classified as a Class I carcinogens by the International Agency for Research on Cancer (IARC) [5,6]. Among AFs, aflatoxin B1 (AFB1) is the most toxic and carcinogenic compound known. AFs are commonly relevant to the cereals, nuts and a scope of their agricultural products, especially peanuts, maize and rice [7,8]. In the United States, the Food and Drug Administration (FDA) has set the limiting value at 20 μg/kg for total AFs (B1, B2, G1, G2) for all foods, and 100 μg/kg for peanut and corn feed products [9,10]. In the European Union, the European Commission (EC) set the upper limit at 2 μg/kg for AFB1 and 4 μg/kg for total AFs [11]. The appearance of antifungal resistance of chemical fungicides and the safety requirements of practical application in crops globally have incurred the discovery of novel antifungal agents and other antifungal substances [12].
Many strategies have been used to reduce AFs contamination. Safe and efficient natural substances for preventing and controlling A. flavus growth and AFs production are necessary. Essential and plant hormone possessing potent anti-microbial, antioxidant activities, were applied to agricultural industry [13,14]. However, some fungicides have low solubility in general solvents. As a universal solvent, dimethylformamide (DMF) is effective to dissolve several high-efficiency antifungal agents. Therefore, to determine the inhibitory effect and mechanism of DMF on A. flavus is necessary for basic research. Although DMF have toxicity, previous research showed only inhibited the cell viability in cells exposed to 160 mM DMF (78.7%, p < 0.01) [15].
Transcriptional sequencing (RNA-Seq) has been widely applied to study lots of eukaryotic transcriptomes [16,17,18]. It also has been used to decipher the inhibitory mechanism of eugenol [19], gallic acid [20], and cinnamaldehyde [21] on A. flavus growth and AFs formation. The objective of this study is to investigate the effect of DMF on A. flavus growth and AFs production, and to determine transcriptomic changes in A. flavus treated with DMF. In particular, the inhibitory mechanism of action of DMF on AFs biosynthesis at the transcriptomic level is elucidated. In terms of agricultural applications, this research may provide a basis for synergistic antifungal effect between DMF and new fungicides.

2. Results

2.1. Inhibitory Potential of DMF Acts on Growth and Toxicity by A. flavus

The antifungal effect of DMF on A. flavus is shown in Figure 1. After DMF treatment, the mycelia growth of A. flavus was significantly suppressed in a dose-dependent manner. Under 4% DMF concentration, the maximum growth inhibition was observed (Figure 1A). However, the complete inhibition was not obtained with all the tested concentrations. With 0.25–1% of DMF treatment, the colony growth was retarded compared to the control group (Figure 1B).
Similarly, as shown in Figure 2, DMF significantly inhibited the AFB1 production in a dose-dependent manner in the yeast extract sucrose (YES) broth at the tested levels. Moreover, the production of AFB1 was completely suppressed by 2% and 4% DMF. Obviously, the difference compared with the growth data, DMF has a significant inhibitory effect on AFB1. Taken together, these results indicated that DMF significantly suppressed A. flavus growth and AFB1 production in a dose-dependent manner.

2.2. Changes on Gene Expression Profile of A. flavus Treated with DMF

The transcriptomes of all treatment and control groups were sequenced to obtain a comprehensive overview of the transcriptional response of A. flavus to DMF. Using RNA sequencing, averagely 46.35 million, 48.23 million raw reads were gained from control and 1% of DMF treatment samples, respectively. After percolating, about 44.39 million and 46.39 million clean reads were obtained from control and treatment transcriptome. According to fragments per kilobase per million mapped fragments (FPKM) values with FDR ≤ 0.05 and log2Ratio ≥1 or ≤−1, differentially expression genes (DEGs) between the control and treatment samples were identified. Compared with control, a total of 2353 genes were significantly differentially expressed in the DMF group. Among them, 1204 (51.17%) genes were up-regulated and 1149 (48.83%) genes were down-regulated. These results suggested that DMF effectively influenced the expression of large number of genes.

2.3. Functions and Involved Pathways of Significant DEGs

GO functional enrichment analysis revealed that these significantly differentially expressed genes were mainly involved in molecular function (MF), cellular component (CC), and biological process (BP). Figure 3 shows the top 30 terms of the most obvious enrichment in three gene categories. For the up-regulated genes (Figure 3A), flavin adenine dinucleotide binding, inorganic molecular entity transmembrane, metal ion transmembrane transporter activity were the staple terms in molecular function. Lipid catabolic process and mental iron transport were the most affected terms belonging to the biological process, followed by signaling, signal transduction and intracellular signal transduction. The main terms in cellular component were Golgi-associated vesicle membrane, coated vesicle membrane, vesicle membrane, cytoplasmic vesicle membrane, vesicle coat. The result indicated that the treatment with 1% DMF mainly up-regulated the genes involved in the vesicle and related cellular component. However, most of down-regulated genes were enriched in biological process (BP) (Figure 3B). For the down-regulated genes, organic acid metabolic process, oxoacid metabolic process, carboxylic acid metabolic process, small molecule biosynthetic process and cellular amino acid metabolic process were the abundant terms in biological process. Oxidoreductase activity, protein dimerization activity, NAD binding and electron transfer activity were the predominant terms in molecular function. The most terms belonging to cellular component were mitochondrion, mitochondrial part, organelle envelope and envelope.
By Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, the top 20 enriched pathways of significant DEGs in A. flavus treated with 1% DMF treatment were shown in Figure 4. For up-regulated DEGs with 1% DMF (Figure 4A), the most abundant genes (26 DEGs) were enriched protein processing in endoplasmic reticulum (afv04141), and 20 DEGs, 16 DEGs, 13 DEGs, 11 DEGs were enriched in RNA transport (afv03013), MAPK signaling pathway (afv04011), phenylalanine metabolism (afv00360), beta-alanine metabolism (afv00410), respectively. For down-regulated DEGs (Figure 4B), the most down-regulated genes (72 DEGs) were enriched in biosynthesis of amino acids (afv01230), and 38 DEGs, 31 DEGs, 14 DEGs were enriched in carbon metabolism, 2-Oxocarboxylic acid metabolism (afv01210) and valine, leucine and isoleucine biosynthesis (afv00290). In general, the most important amino acid (such as valine, leucine, isoleucine, arginine and lysine) biosynthesis pathways were suppressed, while some amino acid (alanine, phenylalanine) metabolism and degradation pathways were improved. In addition, the genes involved in fungal basal metabolism (carbon metabolism and nitrogen metabolism) and mitochondrial were also down-regulated.

2.4. Genes Related to Pigment Biosynthesis and Fungal Development

Transcriptional activity of genes involved in A. flavus pigment and development was presented at Table 1. The genes related to pigment biosynthesis AFLA_016120 (O-methyltransferase family protein), AFLA_016130 (a hypothetical protein), AFLA_016140 (conidial pigment biosynthesis scytalone dehydratase Arp1) were all significantly up-regulated. Obviously, consistent with the growth phenotype with exposure to 1% DMF treatment, most gene involved in fungal development were not significantly influenced. A FluG family protein (AFLA_039530), a conidiation-specific family protein (AFLA_044790) and a conidiation protein Con−6 (AFLA_044800) was conspicuously stimulated. APSES (ASM-1, Phd1, StuA, EFG1, and Sok2) transcription factor stuA, developmental regulator flbA, C2H2 conidiation transcription facto flbC and brlA were up-regulated slightly. However, sexual development transcription factor steA, sexual development transcription factor nsdD, G-protein complex alpha subunit gpaA/fadA were down-regulated. LaeA which regulated secondary metabolism was slightly down-regulated. Interestingly, the developmental regulator VeA show slight down-regulated level, as the velvet regulator VosA which plays a pivotal role in spore survival and metabolism in Aspergillus [22] was tightly up-regulated.

2.5. Genes Related to the Biosynthesis of Aflatrem, Aflatoxins, and Cyclopiazonic Acid

The transcription activities of the genes involved in the biosynthesis of aflatrem (#15), aflatoxins (#54), and cyclopiazonic acid (#55) were shown in Table 2. In pathway #15, most of genes were expressed at very low levels, and the expression of AFLA_045470 and AFLA_045530 was undetected. It is worth mentioning that AFLA_045570 encoding a putative MFS multidrug transporter and AFLA_045460 encoding a putative acetyl xylan esterase were significantly down-regulated. However, most of other genes were slightly up-regulated, including the genes encoding the hybrid PKS/NRPS enzyme, integral membrane protein and cytochrome P450. In AFs biosynthesis pathway (54#), all genes were down-regulated by with 1% DMF except aflNa encoding a hypothetical protein and aflA encoding fatty acid synthase alpha subunit were slightly up-regulated. The key regulator genes aflR and aflS were both slightly suppressed with log2 D1/CK values −0.263 and −0.419, respectively. Interestingly, all genes referred to lipid redox were significantly down-regulated, including aflYa, aflY, aflX, aflW, aflP, aflO, aflM, aflJ and aflH. In Pathway 55, AFLA_139470 encoding a FAD dependent oxidoreductase, AFLA_139480 encoding a tryptophan dimethylallyl transferase and AFLA_139480 encoding a hybrid PKS/NRPS enzyme were significantly up-regulated with 1% DMF treatment, whereas AFLA_139460 encoding an MFS multidrug transporter was obviously down-regulated.

2.6. Genes Involved in Cell Wall

The transcription activities of the genes involved in cell wall were shown in Table 3. AFLA_038420 encoding an endo-chitosanase B and AFLA_024770 encoding a putative symbiotic chitinase were significantly up-regulated with log2 D1/CK values 4.6754 and 2.2492, respectively. Rather than, the genes related to chitin hydrolysis and chitin synthesis were both significantly down-regulated. It is worthy that the down-regulated transcriptional levels of gene involved in chitinase were much higher than chitin synthase. All genes related to glucan synthesis were up-regulated, AFLA_023460 encoding an alpha−1,3-glucan synthase Ags1 and AFLA_134100 encoding an alpha−1,3-glucan synthase Ags2 were up-regulated with log2 D1/CK values 1.8014 and 0.7643, respectively. All genes related to glucan hydrolysis were significantly down-regulated, including AFLA_095680 encoding a putative alpha−1,3-glucanase, AFLA_029950 encoding a putative endo−1, 3(4)-beta-glucanase, AFLA_045290 encoding a putative extracellular endoglucanase/cellulase, AFLA_102640 encoding a putative exo-beta−1, 3-glucanase, AFLA_053390 encoding a GPI-anchored cell wall beta−1,3-endoglucanase EglC, AFLA_068300 encoding a 1, 3-beta-glucanosyltransferase Bgt1.

2.7. Genes Involved in Glucose Metabolism Pathway

The transcription activities of the genes involved in glucose metabolism pathway was shown in Table S1. All genes involved in glucose metabolism were down-regulated at different degrees. In glycolysis pathway, AFLA_101470 encoding a putative glyceraldehyde−3-phosphate dehydrogenase, AFLA_085400 encoding a 2,3-bisphosphoglycerate-independent phosphoglycerate mutase, AFLA_069370 encoding a putative phosphoglycerate kinase PgkA and AFLA_119290 encoding a putative phosphofructokinase were significantly down-regulated with log2 D1/CK values −5.1528, −1.1066, −1.0872 and −1.0641, respectively. In the tricarboxylic acid (TCA) cycle and glyoxylic acid cycle, AFLA_052400 encoding a isocitrate lyase AcuD, AFLA_086400 encoding a putative isocitrate dehydrogenase Idp1 and AFLA_069370 encoding a putative phosphoglycerate kinase PgkA were obviously down-regulated with log2 D1/CK values −1.6033, −1.1747 and −1.0872, respectively.

2.8. Genes Involved in Oxidative Phosphorylation and Amino Acid Biosynthesis/Metabolism

AFs biosynthesis may be controlled by the energy state of specific subcellular compartments, and the production of these secondary metabolites may be affected by the synthesis of ATP in mitochondria [23]. As shown in Figure 5, there are five complexes involved in the oxidative phosphorylation, including complex I, II, III, IV, and V. In each complex, several genes encoding NADH dehydrogenase, succinate dehydrogenase, cytochrome oxidase and ATPase, were down-regulated at different degrees by 1% DMF. In addition, the biosynthesis pathways and metabolism pathways of almost all 20 essential amino acid were affected by 1% DMF, such as phenylalanine, tryptophan, tyrosine, valine, leucine, isoleucine, alanine, glycine, serine, threonine, and cysteine. It is noteworthy that most of amino acids metabolism related genes were up-regulated, while most of amino acids biosynthesis genes were down-regulated in A. flavus treated with 1% DMF. It indicates that the amino acids cannot be synthesized, then some important intermediate products cannot be accumulated with putting a lot of pressure on the fungal cells.

2.9. Genes Involved in MAPK Pathway, Oxylipins, GPCRs and Oxidative Stress Response

As shown in Table S2, the expression changes of most genes in oxidative stress response (OSR), MAPK pathway, oxylipins and GPCRs were slightly changed after 1% DMF treatment. The spore-specific catalase CatA was significantly up-regulated, while catalase Cat was suppressed. For MAP kinase, the sakA1 and sakA2 were both up-regulated by DMF with Log2 (FPKM) values 0.678 and 1.837, respectively. The fatty acid oxygenase ppoA, ppoB and ppoC were all significantly up-regulated. Interestingly, gfdB encoding a glycerol 3-phosphate dehydrogenase, was clearly down-regulated. The G protein-coupled receptor gene gprG encoding a PQ loop repeat protein was significantly down-regulated after DMF treatment, while other receptor genes gprH, gprK, gprM were significantly up-regulated.

3. Discussion

The biosynthesis of toxic and carcinogenic AFs involves multiple biochemical reactions, which require the activity of more than various 27 enzymes [24]. These enzymes are coded by the genes grouped in a cluster of aflatoxin pathway, as their expression regulated by cluster-specific transcription activator aflR and transcription enhancer aflS [25,26]. After DMF treatment, all genes in the cluster were down-regulated at different degrees except for aflA encoding fatty acid synthase alpha subunit. However, DMF could not completely inhibit any genes in the cluster. Two crucial regulator genes aflS and aflR were slightly repressed in A. flavus with DMF treatment, accompanying with significant reduction in the expression of AFs structural genes such as aflD, aflH, aflI, aflJ, aflL, aflM, aflO, aflP, aflQ, aflW, aflX and aflY, leading to a consecutive loss in the ability to synthesize AFs intermediates [24]. These findings suggest that the expression changes of structural genes are more significant compared with the key regulators aflR and aflS with inhibitor treatment.
Similar findings were obtained in A.flavus treated with eugenol [17,19], piperine [27], cinnamaldehyde [21,28], ethanol [29], 5-Azacytidine (5-AC) [30] and gallic acid [20]. In A. flavus treated with 5-AC, the expressions of aflR and aflS were basically unchanged. However, at least three structural genes including aflQ, aflI and aflLa were completely inhibited, and five structural genes were suppressed with high or middle levels by 5-AC, especially aflG and aflX [30]. After treatment with eugenol, the expression of aflR did not change obviously and the expression of aflS was slightly up-regulated. The expression of most structural genes was down-regulated and the most strongly down-regulated gene was aflMa, followed by aflI, aflJ, aflCa, aflH, aflNa, aflE. However, the expression of these genes was not completely inhibited [31]. Exposed to cinnamaldehyde, aflR and aflS showed slightly up-regulation. Excepting the up-regulated expression of aflF, all the structural genes in the cluster were down-regulated. The most strongly down-regulation gene was aflD, followed by the key structural genes aflG, aflH, aflP, aflM, aflI, aflL and aflE [21]. When A. flavus treated with antioxidant gallic acid, aflR and aflS were slightly up-regulated while structure gene showed down-regulation [20]. After treatment with 3.5% ethanol, aflR, aflS were all down-regulated significantly and aflS/aflR showed the up-regulation. At least three structural genes including aflK, aflLa, aflL, aflG and aflM were completely inhibited following the up-regulation of aflS/aflR [28]. These similar results showed that after treatment with different anti-aflatoxigenic compounds, the AFs cluster transcription factors aflR and aflS were not significantly changed although many structural genes were significantly repressed, suggesting the stable expression of the two key regulator genes in A. flavus.
Sugars metabolism produce the basic substance unit acetyl-CoA, is prerequisite for all known polyketides compounds especially AFs [32,33]. Acetyl-CoA is mainly generated by glycolysis pathway in cytoplasm and fatty acid β-oxidation pathway in peroxisomes [34,35]. The reduction of AFs biosynthesis was associated with the decrease of tricarboxylic acid (TCA) cycle intermediates, the suppression of fatty acid biosynthesis, and the increase of pentose phosphate pathway substance, as reported [36]. The numerous DEGs were enriched in carboxylic acid metabolic process. The utilization of a given carbon source in A. flavus involves the sugar cluster, of which aflYe, aflYd, aflYc, aflYb were up-regulated. Additionally, a zinc finger transcription repressor creA was slightly down-regulated, creA suppresses the expression of aclR and the latter is a positive regulatory factor for the genes aldA [37,38]. It is showed that the DMF treatment contributes to the significant increase of pentose phosphate pathway of genes expression, such as AFLA_079220 (glucose dehydrogenase), AFLA_026950 (3-ketoacyl-CoA thiolase peroxisomal A precursor) from our research. The similar result was reported that gallic acid inhibited the AFs formation via up-regulation of pentose phosphate pathway [20]. Cinnamaldehyde, cultured in AFs inhibitory medium, the pentose phosphate pathway was accelerated, leading to NADPH accumulation and AFs reduction [21,28]. On the other hand, a large number of genes involved in TCA cycle-related and glycolysis pathway showed significant down-regulation. Interestingly, the glyoxylic acid-related genes especially AFLA_052400 encoding an isocitrate lyase AcuD, AFLA_049390 (malate synthase AcuE) and the TCA cycle-related genes especially AFLA_086400 (socitrate dehydrogenase Idp1) were significantly inhibited, leading to the accumulation of isocitrate and the subsequent depletion of acetyl-CoA. Taken together, the inhibition of AFs may be due to the depletion of acetyl-CoA and the lack of NADPH.
One important factor that has been found to affect AFs biosynthesis is amino acid catabolism and biosynthesis [39]. Wilkinson et al. [39] reported that glutamine and tyrosine favor AFs production in A. flavus, while tryptophan seem more complicated. Adye et al. and Naik, M et al. [40,41] reported that tyrosine, tryptophan, phenylalanine and methionine were easily absorbed into AFs biosynthesis pathway of A. flavus. Similarly, Wilkinson et al. [39] found the supplementation of amino acid in YES media could positively modulate the AFs biosynthesis in A. flavus and A. parasitiucs. In addition to these aromatic amino acids which easily influenced the AFs biosynthesis, Payne et al. [42] unearthed that proline and asparagine can increase more AFs production than tryptophan or methionine in A. flavus. Roze et al. [43] found that veA negatively regulated catabolism of branched chain amino acids and the synthesis of ethanol in A. parasiticus. The analogous results existing in our study were that veA was slightly suppressed, and the up-regulated and down-regulated DEGs were abundantly enriched in metabolism and biosynthesis pathways, respectively. However, 2-PE at a low level rendered the decreased activities in the metabolism of branched-chain amino acid, of which may be necessary to activate the AFs pathway by providing building blocks and energy regeneration [44]. These results implied that amino acids played complex roles in AFs production.
Oxidative phosphorylation in mitochondria can convert the energy released by organisms into ATP during the decomposition process. Chung et al. [45] reported that DMF influenced electron-transfer reactivity of cytochrome b5. As a powerful inhibitor of electron transport, DMF was found to have a marked inhibiting effect of the phosphorylation reaction at concentrations lower than 6.0% (v/v) [46]. Consistent with the result, our RNA-seq data showed that the expressional levels of several genes in oxidative phosphorylation were consistently down-regulated, including complex I (NADH complex), complex II (dehydrogenase complex), complex III (cytochrome complex) and complex IV (cytochrome oxidase). All genes related to electron transport component were repressed at different degrees, especially AFLA_129610 encoding a putative subunit G of NADH-ubiquinone oxidoreductase was obviously down-regulated. Inhibition of oxidative phosphorylation by the application of exogenous compounds such as resveratrol has been shown to compromise fungal oxidative stress tolerance by altering mitochondrial respiration and oxidative stress as a prerequisite for AFs production by Aspergillus parasiticus [47,48]. These results suggested that oxidative phosphorylation dysfunction might be associated with the reduction of AFs biosynthesis.
The cell wall acts as a protective barrier for fungi against environmental factors and is essential for the survival of the fungus during development and reproduction [49]. The cell wall plays an important molecular target for various antifungal compounds [50]. The important components of fungal cell wall, alpha−1,3-glucan and beta−1,3-glucan, play a crucial role to maintain the normal morphology of fungal cell wall. The alpha−1,3-glucan synthase encoded by ags1, ags2, ags3 is essential for the formation of alpha−1,3-glucan [51]. fksP encoding the beta−1,3-glucan synthase participates in the synthesis of beta−1,3-glucan [52]. In our RNA-seq data, ags1, ags2, ags3 and fksP were significantly up-regulated to resist DMF destruction. As the important constituent ingredient of the cell wall, chitin biosynthesis is directly controlled by chitin synthase [53]. In our current study, all chitin synthase was significantly depressed with DMF treatment. Chitinase facilitates the separation of its cells during fungal growth and reproduction [54]. Wang et al. [55] found the glucanase gene crh11 was down-regulated, could cause the obstruction of fungal reproduction. Similarly, our RNA-seq data displayed that all glucanase genes related to cell wall were down-regulated at different degrees. These results suggest that DMF attacks the cell wall of A. flavus, with destroying the main components of the fungal cell wall, chitin and structural polysaccharides. Consequently, the cell will generate a stress response to maintain the basic structure of the cell wall by over-expressing alpha/beta-1,3-glucan synthase genes to resist external stimuli.
Figure 6 showed an elementary diagram illustrating the antifungal effect of dimethylformamide act on A. flavus NRRL3357. To sum up in Figure 6, dimethylformamide inhibits the AFs biosynthesis and fungal growth of A. flavus via (1) attacking the cell wall by regulating the expression of cell wall integrity (CWI) related genes, and then cause the disorder of related protein of the cell membrane, and the spread to oxylipins genes by signal transduction; (2) increasing the depletion of acetyl-CoA and suppressing the NADPH accumulation by glucose metabolism; (3) disturbing the function of oxidative phosphorylation, then reducing ATP deemed as key elements for fungal cell to perform various reactions; (4) weakening amino acid synthesis which are indispensable for AFs accumulation. In brief, the dysfunction of cell wall integrity, glucose metabolism, amino acid biosynthesis, oxidative phosphorylation tight-knit interact with AFs production.

4. Conclusions

This study provides new insights of mechanism to interpret the inhibition of transcriptional regulation with DMF against AFs synthesis via a large number of comparative RNA sequencing. Based on existing research, we conclude that (1) DMF attacks the cell wall of A. flavus, with destroying the fungal cell wall integrity and the cell will subsequently generate a stress response to maintain the basic structure of the cell wall by over-expressing glucan synthase genes to resist external stimuli; (2) in the presence of DMF, the most intuitive performance of the decrease of AFs production is following increased expression of their specific regulators aflS/aflR and down-regulation of AFs cluster genes; (3) the down-regulation of the global regulator VeA and up-regulation of FluG is associated with the increase of conidiophore development; (4) glucose metabolism pathways are greatly interfered including TCA cycle, pentose phosphate pathway, glyoxylic acid pathway; (5) oxidative phosphorylation is in disorder, with reducing ATP required by fungal organisms to perform various reactions; and (6) the biosynthesis and metabolism of most amino acid is affected, this indicates that most amino acids cannot be synthesized and some important intermediate products cannot be accumulated, which puts a lot of pressure on the fungal cells. In general, these results strongly suggest that DMF disturb a variety of cellular reactions in A. flavus, thereby interfering fungal growth and metabolic function.

5. Materials and Methods

5.1. Fungal Strain, Chemicals and Treatment

DMF (DMF, 100% purity) was purchased from Beijing Chemical Works (Beijing, China). Chromatographic grade methanol was purchased from Thermo Fisher Scientific (Waltham, MA, USA). The AFB1 standard was purchased from Sigma-Aldrich (Sigma-Aldrich Chemicals, St. Louis, MO, USA).
The A. flavus strain NRRL3357 used in this study [29] was maintained in the dark condition on potato dextrose agar (PDA) medium at 4 °C as reserving. A conidia inoculum was prepared by washing PDA surface culture and adjusted to 107 conidia/mL with 0.1% Tween-80 solution. The AFB1 standard was dissolved in 70% methanol.

5.2. Determination of Fungal Growth and AFB1 Production

Different treatment of DMF was added to the sterilized PDA medium at final concentrations of 0.25%, 0.5%, 1%, 2%, and 4% respectively. Then, 10 µL of 107 conidia/mL suspension was inoculated on the central of PDA medium and incubated at 28 °C for 7 days. Determination of A. flavus growth indexes were by measuring colony diameters.
Similarly, the different treatment concentrations of DMF were added to yeast extract sucrose (YES, Hopebio, Qingdao, China) broth to obtain the concentrations of 0.25%, 0.5%, 1.0%, 2.0% and 4.0% DMF. The control cultures (CK) were treated similarly but without DMF. Then, 100 µL of 107 conidia/mL suspension was inoculated to 100 mL YES broth. Fungal mycelia were collected and weighed as the method described by Yamazaki et al. [56] after incubation at 28 °C and 180 rpm in the dark for 5 days. Extraction and quantification of AFB1 by high-performance liquid chromatography (HPLC) was conducted according to the previous reference [57]. AFB1 was extracted with acetonitrile: water (84:16) solution from YES broth and purified by a ToxinFast immunoaffinity column (Huaan Magnech Biotech, Beijing, China). An Agilent 1220 Infinity II system coupled with a fluorescence detector (Santa Clara, CA, USA), an Agilent TC-C18 column (250 mm × 4.6 mm, 5 μm particle size) and a post-column derivation system (Huaan Magnech, Beijing, China) was used to quantify the AFB1 concentrations. Each treatment was performed in triplicate.

5.3. Construction of cDNA Libraries and RNA Sequencing

Total RNA extraction, cDNA libraries construction was prepared according to the methods described by Lin et al. [30]. Illumina® HiSeq 4000TM system (San Diego, CA, USA) was used to sequence the cDNA libraries. The original RNA-seq data have been uploaded in the NCBI Sequence Read Archive (SRA) with accession code SUB8228352.

5.4. Analysis of Sequence Data

After removing the false reads and deleting the end-sequence with low quality, the reads shorter than 50 bp was discarded. The remaining reads were mapped to the A. flavus genome (http://www.ncbi.nlm.nih.gov/nuccore/AAIH00000000). Then, the high-quality reads were assembled into unigenes by using the method described by Grabherr et al. [58]. The transcriptional levels of genes in A. flavus were represented using Fragments per kilobase per million mapped fragments (FPKM) values [59]. The mean FPKM of triplicate samples was analyzed by using DEseq software [60] for the differential expression of genes. The significant differential expression genes were identified as log2Ratio ≥ 1 and q < 0.05 between these compared samples [20]. For the DEGs, gene ontology (GO) functional analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment were conducted using FungiFun and KAAs, respectively [61,62,63].

5.5. RT-qPCR Analysis of AFs Biosynthesis Genes

All genes in the AFs biosynthesis cluster were chosen for RT-qPCR validation of the RNA-Seq results according to the methods described by Ren et al. [29]. All the data generated from real-time PCR were analyzed using SPSS software version 16.0 with one-way ANOVA method. The significance level of 0.05 has been indicated with lowercases. The gene was defined as significantly up- or down-regulated only if the relative expression level was more than two-fold and showed significant at 0.05 level compared to the control group [64].

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6651/12/11/683/s1, Table S1. Transcriptional activity of genes involved in A. flavus carbon metabolism; Table S2. Transcriptional activity of genes involved in A. flavus MAPK pathway, Oxylipins, GPCRs and OSR

Author Contributions

Conceptualization, P.C. and F.X.; data curation, L.P.; formal analysis, L.P., P.C., J.J. and Q.Y.; methodology, J.J. and F.X.; writing—original draft, L.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Key R&D Program of China (No. 2016YFD0400105), National Peanut Industrial Technology System (CARS−13) and the National Natural Science Foundation of China (31972179). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rajasekaran, K.; Sayler, R.J.; Sickler, C.M.; Majumdar, R.; Jaynes, J.; Cary, J.W. Control of Aspergillus flavus growth and aflatoxin production in transgenic maize kernels expressing a tachyplesin-derived synthetic peptide, AGM182. Plant Sci. 2018, 270, 150–156. [Google Scholar] [CrossRef] [PubMed]
  2. Klich, M.A. Environmental and developmental factors influencing aflatoxin production by Aspergillus flavus and Aspergillus parasiticus. Mycoscience 2007, 48, 71–80. [Google Scholar] [CrossRef]
  3. Groopman, J.D.; Kensler, T.W.; Wild, C.P. Protective interventions to prevent aflatoxin-induced carcinogenesis in developing countries. Annu. Rev. Public Health 2008, 29, 187–203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Reddy, K.; Salleh, B.; Saad, B.; Abbas, H.; Abel, C.; Shier, W. An overview of mycotoxin contamination in foods and its implications for human health. Toxin Rev. 2010, 29, 3–26. [Google Scholar] [CrossRef]
  5. International Agency for Research on Cancer [IARC]. Some naturally occurring substances: Food items and constituents, heterocyclic aromatic amines and mycotoxins. IARC Monogr. Eval. Carcinog. Risks Hum. 1993, 56, 245–395. [Google Scholar]
  6. International Agency for Research on Cancer [IARC]. Some traditional herbal medicines, some mycotoxins, naphthalene and styrene. IARC Monogr. Eval. Carcinog. Risks Hum. 2002, 82, 1–556. [Google Scholar]
  7. Luttfullah, G.; Hussain, A. Studies on contamination level of aflatoxins in some dried fruits and nuts of pakistan. Food Control. 2011, 22, 426–429. [Google Scholar] [CrossRef]
  8. Atayde, D.D.; Reis, T.A.; Ignácio, J.G.; Zorzete, P.; Reis, G.M.; Benedit, C. Mycobiota and aflatoxins in a peanut variety grown in different regions in the state of so paulo, brazil. Crop Prot. 2012, 33, 7–12. [Google Scholar] [CrossRef]
  9. The U.S. Food and Drug Administration [FDA]. CPG Sec. 570.375 Aflatoxin in Peanuts and Peanut Products; 2005. Available online: https://www.fda.gov/media/72073/download (accessed on 29 October 2020).
  10. The U.S. Food and Drug Administration [FDA]. CPG Sec 555.400 Foods-Adulteration with Aflatoxin; 2005. Available online: https://www.fda.gov/regulatory-information/search-fda-guidance-documents/cpg-sec-555400-foods-adulteration-aflatoxin (accessed on 29 October 2020).
  11. Commission Regulation (EU) No 165/2010 of 26 February 2010 Amending Regulation (EC) No 1881/2006 Setting Maximum Levels for Certain Contaminants in Foodstuffs as Regards Aflatoxins (Text with EEA Relevance) OJ L 50, 27.2. 2010, pp. 8–12 (BG, ES, CS, DA, DE, ET, EL, EN, FR, IT, LV, LT, HU, MT, NL, PL, PT, RO, SK, SL, FI, SV) Special Edition in Croatian: Chapter 13 Volume 038 pp. 244–248. Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/?uri=CELEX%3A32010R0165&qid=1603901359517 (accessed on 29 October 2020).
  12. Bhatnagar-Mathur, P.; Sunkara, S.; Bhatnagar-Panwar, M.; Waliyar, F.; Sharma, K.K. Biotechnological advances for combating Aspergillus flavus and aflatoxin contamination in crops. Plant Sci. 2015, 234, 119–132. [Google Scholar] [CrossRef] [Green Version]
  13. Ceker, S.; Agar, G.; Alpsoy, L.; Nardemir, G.; Kizil, H.E. Antagonistic effects of Satureja hortensis essential oil against AFB1 on human lymphocytes in vitro. Cytol. Genet. 2014, 48, 327–332. [Google Scholar] [CrossRef]
  14. Abass, M.H. In vitro antifungal activity of different plant hormones on the growth and toxicity of nigrospora spp. on date palm (Phoenix Dactylifera L.). Open Plant Sci. J. 2017, 10, 10–20. [Google Scholar] [CrossRef]
  15. Xu, L.; Zhao, Q.W.; Luo, J.; Ma, W.L.; Jin, Y.; Li, C.H.; Hou, Y.F.; Feng, M.Y.; Wang, Y.; Chen, J.; et al. Integration of proteomics, lipidomics, and metabolomics reveals novel metabolic mechanisms underlying N, N-dimethylformamide induced hepatotoxicity. Ecotoxicol. Environ. Saf. 2020, 205, 111166. [Google Scholar] [CrossRef] [PubMed]
  16. Watt, S.; Schubert, F.; Wood, V.; Goodhead, I.; Penkett, C.J.; Rogers, J.; Bhler, J.; Wilhelm, B.T.; Marguerat, S. Dynamic repertoire of a eukaryotic transcriptome surveyed at single-nucleotide resolution. Nature 2008, 453, 1239–1243. [Google Scholar]
  17. Wang, Z.; Gerstein, M.; Snyder, M. RNA-Seq: A revolutionary tool for transcriptomics. Nat. Rev. Genet. 2009, 10, 57–63. [Google Scholar] [CrossRef]
  18. Borries, A.; Vogel, J.; Sharma, C.M. Differential RNA Sequencing(dRNA-Seq): Deep-Sequencing-Based Analysis of Primary Transcriptomes. Tag-Based Next Generation Sequencing; Wiley-VCH Verlag Gmbh&Co. KGaA: Weinheim, Germany, 2012. [Google Scholar]
  19. Lv, C.; Wang, P.; Ma, L.; Zheng, M.; Liu, Y.; Xing, F. Large-scale comparative analysis of eugenol-induced/repressed genes expression in Aspergillus flavus using RNA-seq. Front. Microbiol. 2018, 9, 1116. [Google Scholar] [CrossRef]
  20. Zhao, X.; Zhi, Q.Q.; Li, J.Y.; Keller, N.P.; He, Z.M. The antioxidant gallic acid inhibits aflatoxin formation in Aspergillus flavus by modulating transcription factors FarB and CreA. Toxins 2018, 10, 270. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Wang, P.; Ma, L.; Jin, J.; Zheng, M.; Pan, L.; Zhao, Y. The antiaflatoxigenic mechanism of cinnamaldehyde in Aspergillus flavus. Sci. Rep. 2019, 9, 10499. [Google Scholar] [CrossRef] [PubMed]
  22. Sarikaya-Bayram, Ö.; Palmer, J.M.; Keller, N.; Braus, G.H.; Bayram, Ö. One Juliet and four Romeos: VeA and its methyltransferases. Front. Microbiol. 2015, 6, 1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Buchanan, R.L.; Jones, S.B.; Gerasimowicz, W.V.; Zaica, L.L.; Stahl, H.G.; Ocker, L.A. Regulation of aflatoxin biosynthesis: Assessment of the role of cellular energy status as a regulator of the induction of aflatoxin production. Appl. Environ. Microbiol. 1987, 53, 1224–1231. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Ehrlich, K.C. Predicted roles of the uncharacterized clustered genes in aflatoxin biosynthesis. Toxins 2009, 1, 37–58. [Google Scholar] [CrossRef] [Green Version]
  25. Chang, P. The Aspergillus parasiticus protein AFLJ interacts with the aflatoxin pathway-specific regulator AFLR. Mol. Genet. Genom. 2003, 268, 711–719. [Google Scholar] [CrossRef] [PubMed]
  26. Price, M.S.; Yu, J.; Nierman, W.C.; Kim, H.; Pritchard, B.; Jacobus, C.A.; Bhatnagar, D.; Cleveland, T.E.; Payne, G.A. The aflatoxin pathway regulator AFLR induces gene transcription inside and outside of the aflatoxin biosynthetic cluster. FEMS Microbiol. Lett. 2006, 255, 275–279. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Caceres, I.; El Khoury, R.; Bailly, S.; Oswald, I.P.; Puel, O.; Bailly, J.D. Piperine inhibits aflatoxin B1 production in Aspergillus flavus by modulating fungal oxidative stress response. Fungal Genet. Biol. 2017, 107, 77–85. [Google Scholar] [CrossRef]
  28. Yin, H.B.; Chen, C.H.; Kollanoor-Johny, A.; Darre, M.J.; Venkitanarayanan, K. Controlling Aspergillus flavus and Aspergillus parasiticus growth and aflatoxin production in poultry feed using carvacrol and trans-cinnamaldehyde. Poult. Sci. 2015, 94, 2183–2190. [Google Scholar] [CrossRef]
  29. Ren, Y.; Jin, J.; Zheng, M.; Yang, Q.; Xing, F. Ethanol inhibits aflatoxin B1 biosynthesis in Aspergillus flavus by up-regulating oxidative stress-related genes. Front. Microbiol. 2020, 10, 2946. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. 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. [Google Scholar] [CrossRef]
  31. Georgianna, D.R.; Payne, G.A. Genetic regulation of aflatoxin biosynthesis: From gene to genome. Fungal Genet. Biol. 2009, 46, 113–125. [Google Scholar] [CrossRef]
  32. Davis, N.D.; Diener, U.L. Growth and aflatoxin production by Aspergillus parasiticus from various carbon sources. Appl. Microbiol. 1968, 16, 158–159. [Google Scholar] [CrossRef] [Green Version]
  33. Maggio-Hall, L.A.; Wilson, R.A.; Keller, N.P. Fundamental contribution of beta-oxidation to polyketiden mycotoxin production in planta. Mol. Plant Microbe Interact. 2005, 18, 783–793. [Google Scholar] [CrossRef] [Green Version]
  34. Maggio-Hall, L.A.; Keller, N.P. Mitochondrial β-oxidation in Aspergillus nidulans. Mol. Microbiol. 2004, 54, 1173–1185. [Google Scholar] [CrossRef]
  35. Reverberi, M.; Punelli, M.; Smith, C.A.; Zjalic, S.; Scarpari, M.; Scala, V.; Cardinali, G.; Aspite, N.; Pinzari, F.; Payne, G.A.; et al. How peroxisomes affect aflatoxin biosynthesis in Aspergillus flavus. PLoS ONE 2012, 7, e48097. [Google Scholar] [CrossRef] [PubMed]
  36. Yan, S.; Liang, Y.; Zhang, J.; Liu, C. Aspergillus flavus grown in peptone as the carbon source exhibits spore density-and peptone concentration-dependent aflatoxin biosynthesis. BMC Microbiol. 2012, 12, 106. [Google Scholar] [CrossRef] [Green Version]
  37. Ruijter, G.J.G.; Visser, J. Carbon repression in Aspergilli. FEMS Microbiol. Lett. 1997, 151, 103–114. [Google Scholar] [CrossRef] [PubMed]
  38. Shroff, R.A.; Lockington, R.A.; Kelly, J.M. Analysis of mutations in the creA gene involved in carbon catabolite repression in Aspergillus nidulans. Can. J. Microbiol. 1996, 42, 950–959. [Google Scholar] [CrossRef]
  39. Wilkinson, J.R.; Yu, J.; Bland, J.M.; Nierman, W.C.; Bhatnagar, D.; Cleveland, T.E. Amino acid supplementation reveals differential regulation of aflatoxin biosynthesis in Aspergillus flavus NRRL 3357 and Aspergillus parasiticus SRRC 143. Appl. Microbiol. Biotechnol. 2007, 74, 1308–1319. [Google Scholar] [CrossRef]
  40. Adye, J.C.; Mateles, R.I. Incorporation of labeled compound into aflatoxin. Biochim. Et Biophys. Acta 1964, 86, 418–420. [Google Scholar] [CrossRef]
  41. Naik, M.; Modi, V.V.; Patel, N.C. Studies on aflatoxin synthesis in Aspergillus flavus. Indian J. Exp. Biol. 1970, 8, 345–346. [Google Scholar] [PubMed]
  42. Payne, G.A.; Hagler, W.M. Effect of specific amino acids on growth and aflatoxin production by Aspergillus parasiticus and Aspergillus flavus in defined media. Appl. Environ. Microbiol. 1983, 46, 805–812. [Google Scholar] [CrossRef] [Green Version]
  43. Roze, L.V.; Chanda, A.; Laivenieks, M.; Beaudry, R.M.; Artymovich, K.A.; Koptina, A.V.; Awad, D.W.; Valeeva, D.; Jones, A.D.; Linz, J.E. Volatile profiling reveals intracellular metabolic changes in Aspergillus parasiticus: veA regulates branched chain amino acid and ethanol metabolism. BMC Biochem. 2010, 11, 33. [Google Scholar] [CrossRef] [Green Version]
  44. Chang, P.K.; Hu, S.S.; Sarreal, S.B.; Li, R.W. Suppression of aflatoxin biosynthesis in Aspergillus flavus by 2-phenylethanol is associated with stimulated growth and decreased degradation of branched-chain amino acids. Toxins 2015, 7, 3887–3902. [Google Scholar] [CrossRef] [Green Version]
  45. Chung, H.; Lu, J.; Zhang, W.; Suzukizuki, I.; Li, G. Enhanced electron-transfer reactivity of cytochrome b5 by Dimethylsulfoxide and N, N′-Dimethylformamide. Anal. Sci. 2002, 18, 1031–1033. [Google Scholar]
  46. Pe~Na-Valdivia, C.B.; Rodriguez-Flores, L.; Gómez-Puyou, M.T.; Lotina-Hennsen, B. Inhibition of photophosphorylation and electron transport by N, N-dimethylformamide. Biophys. Chem. 1991, 41, 169–174. [Google Scholar] [CrossRef]
  47. Madrigal-Perez, L.A.; Ramos-Gomez, M. Resveratrol inhibition of cellular respiration: New paradigm for an old mechanism. Int. J. Mol. Sci. 2016, 17, 368. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Jayashree, T.; Subramanyam, C. Oxidative stress as a prerequisite for aflatoxin production by Aspergillus parasiticus. Free. Radic. Biol. Med. 2000, 29, 981–985. [Google Scholar] [CrossRef]
  49. Lee, M.J.; Sheppard, D.C. Recent advances in the understanding of the Aspergillus fumigatus cell wall. J. Microbiol. 2016, 54, 232–242. [Google Scholar] [CrossRef]
  50. Dichtl, K.; Helmschrott, C.; Dirr, F.; Wagener, J. Deciphering cell wall integrity signalling in Aspergillus fumigatus: Identification and functional characterization of cell wall stress sensors and relevant Rho GTPases. Mol. Microbiol. 2012, 83, 506–519. [Google Scholar] [CrossRef]
  51. Henry, C.; Latgé, J.P.; Beauvais, A. α-1, 3-glucans are dispensable in Aspergillus fumigatus. Eukaryot. Cell 2012, 11, 26–29. [Google Scholar] [CrossRef] [Green Version]
  52. Meetei, P.A.; Rathore, R.S.; Prabhu, N.P.; Vindal, V. In silico screening for identification of novel β-1, 3-glucansynthase inhibitors using pharmacophore and 3D-QSAR methodologies. Springerplus 2016, 5, 965. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Henar, M.V.; Durán, A.; Roncero, C. Chitin synthases in yeast and fungi. EXS 1999, 87, 55–69. [Google Scholar]
  54. Yang, M.; Lu, L.; Li, S.; Zhang, J.; Li, Z.; Wu, S.; Guo, Q.; Liu, H.; Wu, C. Transcriptomic insights into benzenamine effects on the development, aflatoxin biosynthesis, and virulence of Aspergillus flavus. Toxins 2019, 11, 70. [Google Scholar] [CrossRef] [Green Version]
  55. Wang, Y.; Feng, K.; Yang, H.; Zhang, Z.; Yuan, Y.; Yue, T. Effect of cinnamaldehyde and citral combination on transcriptional profile, growth, oxidative damage and patulin biosynthesis of Penicillium expansum. Front. Microbiol. 2018, 9, 597. [Google Scholar] [CrossRef] [PubMed]
  56. Yamazaki, H.; Yamazaki, D.; Takaya, N.; Takagi, M.; Ohta, A.; Horiuchi, H. A chitinase gene, chiB, involved in the autolytic process of Aspergillus nidulans. Curr. Genet. 2007, 51, 89. [Google Scholar] [CrossRef]
  57. Liang, D.; Xing, F.; Selvaraj, J.N.; Liu, X.; Wang, L.; Hua, H.; Zhou, L.; Zhao, Y.; Wang, Y.; Liu, Y. Inhibitory effect of cinnamaldehyde, citral and eugenol on aflatoxin biosynthetic gene expression and aflatoxin B1 biosynthesis in Aspergillus flavus. J. Food Sci. 2015, 80, M2917–M2924. [Google Scholar] [CrossRef]
  58. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Mortazavi, A.; Williams, B.A.; McCue, K.; Schaeffer, L.; Wold, B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 2008, 5, 621–628. [Google Scholar] [CrossRef] [PubMed]
  60. Anders, S.; Huber, W. Differential expression analysis for sequence count data. Genome Biol. 2010, 11, R106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Kanehisa, M.; Araki, M.; Goto, S.; Hattori, M.; Hirakawa, M.; Itoh, M.; Katayama, T.; Kawashima, S.; Okuda, S.; Tokimatsu, T.; et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, 36, D480–D484. [Google Scholar] [CrossRef] [PubMed]
  62. Priebe, S.; Linde, J.; Albrecht, D.; Guthke, R.; Brakhage, A.A. FungiFun: A web-based application for functional categorization of fungal genes and proteins. Fungal Genet. Biol. 2011, 48, 353–358. [Google Scholar] [CrossRef]
  63. Li, X.; Jiang, Y.; Ma, L.; Ma, X.; Liu, Y.; Shan, J.; Ma, K.; Xing, F. Comprehensive transcriptome and proteome analyses reveal the modulation of aflatoxin production by Aspergillus flavus on different crop substrates. Front. Microbiol. 2020, 11, 1497. [Google Scholar] [CrossRef]
  64. Lu, L.X.; Zhou, F.; Zhou, Y.; Fan, X.L.; Ye, S.F.; Wang, L.; Chen, H.; Lin, Y.J. Expression profile analysis of the polygalacturonase-inhibiting protein genes in rice and their responses to phytohormones and fungal infection. Plant Cell Rep. 2012, 31, 1173–1187. [Google Scholar] [CrossRef]
Figure 1. Inhibitory effect of dimethylformamide on fungal growth of A. flavus NRRL3357. (A) After 6 days of inoculation with A. flavus conidia suspension (107), the morphology of A. flavus colony on PDA medium under different concentrations (0% to 4%) of dimethylformamide. (B) The colony diameter of A. flavus treated with dimethylformamide (0 to 4%). CK: control group. Compared with CK, * p < 0.05, ** p < 0.01.
Figure 1. Inhibitory effect of dimethylformamide on fungal growth of A. flavus NRRL3357. (A) After 6 days of inoculation with A. flavus conidia suspension (107), the morphology of A. flavus colony on PDA medium under different concentrations (0% to 4%) of dimethylformamide. (B) The colony diameter of A. flavus treated with dimethylformamide (0 to 4%). CK: control group. Compared with CK, * p < 0.05, ** p < 0.01.
Toxins 12 00683 g001
Figure 2. Inhibitory effect of dimethylformamide on AFB1 production and fungal growth of A. flavus NRRL3357. (A) The AFB1 production of A. flavus and the inhibition rate of AFB1 in YES broth at 120 h post-treatment. (B) The mycelium weight of A. flavus in YES broth at 120 h post-treatment. CK: control group. Compared with CK, * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 2. Inhibitory effect of dimethylformamide on AFB1 production and fungal growth of A. flavus NRRL3357. (A) The AFB1 production of A. flavus and the inhibition rate of AFB1 in YES broth at 120 h post-treatment. (B) The mycelium weight of A. flavus in YES broth at 120 h post-treatment. CK: control group. Compared with CK, * p < 0.05, ** p < 0.01, *** p < 0.001.
Toxins 12 00683 g002
Figure 3. Go functional enrichment of up-regulated (A) and down-regulated (B) differentially expression genes (DEGs) with 1% dimethylformamide. The ordinate means the -log10 of the control and 1% dimethylformamide treatment. The size of the plot represents the number of DEGs in one GO term; the color of the plot close to red represents more significant enrichment.
Figure 3. Go functional enrichment of up-regulated (A) and down-regulated (B) differentially expression genes (DEGs) with 1% dimethylformamide. The ordinate means the -log10 of the control and 1% dimethylformamide treatment. The size of the plot represents the number of DEGs in one GO term; the color of the plot close to red represents more significant enrichment.
Toxins 12 00683 g003
Figure 4. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment of up-regulated (A) and down-regulated (B) DEGs with 1% dimethylformamide. The ordinate represents the KEGG classification. The size of the plot represents the number of DEGs; the color of the plot close to red represents more significant enrichment in one KEGG term.
Figure 4. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment of up-regulated (A) and down-regulated (B) DEGs with 1% dimethylformamide. The ordinate represents the KEGG classification. The size of the plot represents the number of DEGs; the color of the plot close to red represents more significant enrichment in one KEGG term.
Toxins 12 00683 g004
Figure 5. The diagram of oxidative phosphorylation disorder including all enzymes involved in the oxidative phosphorylation, including complexes I, II, III, IV, V. The green box represents the down-regulation expression of the gene.
Figure 5. The diagram of oxidative phosphorylation disorder including all enzymes involved in the oxidative phosphorylation, including complexes I, II, III, IV, V. The green box represents the down-regulation expression of the gene.
Toxins 12 00683 g005
Figure 6. An elementary diagram illustrating the antifungal effect of dimethylformamide act on A. flavus NRRL3357. Up- or down-regulation expression of cell substance with dimethylformamide exposure is expressed using red and blue arrow, respectively.
Figure 6. An elementary diagram illustrating the antifungal effect of dimethylformamide act on A. flavus NRRL3357. Up- or down-regulation expression of cell substance with dimethylformamide exposure is expressed using red and blue arrow, respectively.
Toxins 12 00683 g006
Table 1. Transcriptional level of genes involved in A. flavus Pigment, development.
Table 1. Transcriptional level of genes involved in A. flavus Pigment, development.
Gene IDCK * (FPKM)D1 * (FPKM)Log2 D1/CKAnnotated Gene Function
AFLA_0161204.0644.503.45O-methyltransferase family protein
AFLA_0161304.7045.643.28hypothetical protein
AFLA_0161402.2516.602.88Arp1 conidial pigment biosynthesis scytalone dehydratase
AFLA_0061800.841.080.34Arb2/brown2 conidial pigment biosynthesis oxidase
AFLA_0093402.114.241.01Mod-A developmental regulator, putative
AFLA_0142600.621.991.67RodB/HypB conidial hydrophobin
AFLA_0983800.110.130.28RodA/RolA conidial hydrophobin
AFLA_0395304.1419.882.26FluG
AFLA_044790154.53520.231.75conidiation-specific family protein
AFLA_04480018.5864.491.80conidiation protein Con-6, putative
AFLA_046990166.39188.530.18StuA APSES transcription factor
AFLA_018340147.87136.85−0.11GpaA/FadA G-protein complex alpha subunit
AFLA_08149036.6423.97−0.61Gda1/VelB nucleoside diphosphatase
AFLA_0210902.064.231.04sporulation associated protein
AFLA_02489025.6832.990.36Fsr1/Pro1 cell differentiation and development protein
AFLA_0296202.754.350.66AbaA transcription factor
AFLA_02690017.7032.200.86VosA developmental regulator
AFLA_066460135.21112.55−0.26VeA developmental regulator
AFLA_03329040.1627.03−0.57LaeA regulator of secondary metabolism
AFLA_13403025.7716.89−0.61developmental regulator FlbA
AFLA_136410158.96139.49−0.19transcriptional regulator Medusa
AFLA_137320122.3473.75−0.73C2H2 conidiation transcription factor FlbC
AFLA_05203011.1417.110.62WetA developmental regulatory protein
AFLA_071090291.36406.280.48EsdC GTP-binding protein
AFLA_07971054.7657.870.08HymA conidiophore development protein
AFLA_0801705.127.750.60FlbD MYB family conidiophore development protein, putative
AFLA_0828501.682.200.39BrlA C2H2 type conidiation transcription factor
AFLA_08311034.5550.120.54conidiation-specific protein (Con-10), putative
AFLA_1019208.4114.18−0.75FluG extracellular developmental signal biosynthesis protein
AFLA_13149062.2762.350.00conserved hypothetical protein
* CK = Control; D1 = 1% dimethylformamide.
Table 2. Transcriptional level of genes involved in the biosynthesis of Aflatrem (#15), Aflatoxins (#54), and Cyclopiazonic Acid (#55).
Table 2. Transcriptional level of genes involved in the biosynthesis of Aflatrem (#15), Aflatoxins (#54), and Cyclopiazonic Acid (#55).
Cluster IDGene IDCK * (FPKM)D1 * (FPKM)Log2 D1/CKAnnotated Gene Function
15AFLA_0454602.990.90−1.73MFS multidrug transporter, putative
15AFLA_04547000NAFAD dependent oxidoreductase, putative
15AFLA_0454800.060.070.20dimethylallyl tryptophan synthase, putative
15AFLA_04549000.03Uphybrid PKS/NRPS enzyme, putative
15AFLA_0455000.040.070.85cytochrome P450, putative
15AFLA_04551000.02Upintegral membrane protein
15AFLA_04552000.07Upintegral membrane protein
15AFLA_04553000NAhypothetical protein
15AFLA_0455400.040.060.76cytochrome P450, putative
15AFLA_0455503.646.430.82hypothetical protein
15AFLA_0455604.576.630.54carboxylic acid transport protein
15AFLA_0455702.230.62−1.85acetyl xylan esterase, putative
54AFLA_1391006.517.620.23aflYe/orf/Ser -Thr protein phosphatase family protein
54AFLA_1391105.187.600.55aflYd/sugR/sugar regulator
54AFLA_1391203.947.220.87aflYc/glcA/glucosidase
54AFLA_1391304.605.980.38aflYb/hxtA/putative hexose transporter
54AFLA_1391402.681.05−1.34aflYa/nadA/NADH oxidase
54AFLA_1391509.792.77−1.82aflY/hypA/hypP/hypothetical protein
54AFLA_13916010.324.95−1.06aflX/ordB/monooxygenase/oxidase
54AFLA_13917015.595.08−1.62aflW/moxY/monooxygenase
54AFLA_13918011.525.79−0.99aflV/cypX/cytochrome P450 monooxygenase
54AFLA_13919010.625.74−0.89aflK/vbs/VERB synthase
54AFLA_1392003.911.43−1.45aflQ/ordA/ord-1/oxidoreductase/cytochrome P450 monooxigenase
54AFLA_13921016.725.02−1.73aflP/omtA/omt-1/O-methyltransferase A
54AFLA_13922027.508.70−1.66aflO/omtB/dmtA/O-methyltransferase B
54AFLA_1392301.470.41−1.81aflI/avfA/cytochrome P450 monooxygenase
54AFLA_1392406.383.60−0.82aflLa/hypB/hypothetical protein
54AFLA_13925010.044.37−1.19aflL/verB/desaturase/P450 monooxygenase
54AFLA_1392606.543.59−0.86aflG/avnA/ord-1/cytochrome P450 monooxygenase
54AFLA_139270176.54181.240.03aflNa/hypD/hypothetical protein
54AFLA_1392804.944.62−0.09aflN/verA/monooxygenase
54AFLA_13929013.897.10−0.96aflMa/hypE/hypothetical protein
54AFLA_13930055.6317.29−1.68aflM/ver-1/dehydrogenase/ketoreductase
54AFLA_13931015.868.02−0.98aflE/norA/aad/adh-2/NOR reductase/dehydrogenase
54AFLA_13932033.0213.07−1.33aflJ/estA/esterase
54AFLA_13933028.7713.33−1.10aflH/adhA/short chain alcohol dehydrogenase
54AFLA_139340108.5390.40−0.26aflS/pathway regulator
54AFLA_13936076.4957.18−0.41aflR/apa-2/afl-2/transcription activator
54AFLA_1393708.477.73−0.13aflB/fas-1/fatty acid synthase beta subunit
54AFLA_1393806.897.990.21aflA/fas-2/hexA/fatty acid synthase alpha subunit
54AFLA_13939039.8519.86−1.00aflD/nor-1/reductase
54AFLA_13940013.558.20−0.72aflCa/hypC/hypothetical protein
54AFLA_13941010.416.66−0.64aflC/pksA/pksL1/polyketide synthase
54AFLA_139420130.55123.18−0.08aflT/aflT/transmembrane protein
54AFLA_13943021.5616.39−0.39aflU/cypA/P450 monooxygenase
54AFLA_13944022.4916.67−0.43aflF/norB/dehydrogenase
55AFLA_139460485.20216.59−1.16MFS multidrug transporter, putative
55AFLA_13947064.20249.281.95FAD dependent oxidoreductase, putative
55AFLA_139480155.11556.661.84dimethylallyl tryptophan synthase, putative
55AFLA_1394900.956.292.72hybrid PKS/NRPS enzyme, putative
* CK = Control; D1 = 1% dimethylformamide; NA = Not applicable, means the FKPM value of the gene in CK group and D1 group were both zero; UP means the FKPM value of the gene in CK group was zero and the transcriptional level of the gene in D1 group was up-regulated compared with CK group.
Table 3. Transcriptional activity of genes involved in A. flavus cell wall.
Table 3. Transcriptional activity of genes involved in A. flavus cell wall.
Gene IDCK * (FPKM)D1 * (FPKM)Log2 D1/CKAnnotated Gene Function
AFLA_0384200.02 0.59 4.68 endo-chitosanase B
AFLA_0247700.89 4.22 2.25 symbiotic chitinase, putative
AFLA_0234605.01 17.45 1.80 alpha-1,3-glucan synthase Ags1
AFLA_1341000.05 0.09 0.76 alpha-1,3-glucan synthase Ags2
AFLA_052800293.07 314.70 0.10 1,3-beta-glucan synthase catalytic subunit FksP
AFLA_0410600.05 0.01 −2.88 cell wall associated protein, putative
AFLA_1046800.05 0.01 −2.27 class V chitinase ChiB1
AFLA_0132800.75 0.22 −1.76 class V chitinase, putative
AFLA_03138088.88 34.17 −1.38 class V chitinase, putative
AFLA_0544700.31 0.25 −0.30 class V chitinase Chi100
AFLA_11476024.13 10.06 −1.26 chitin synthase B
AFLA_0860700.02 0.01 −0.94 chitin synthase, putative
AFLA_06753043.44 29.37 −0.56 chitin biosynthesis protein (Chs7), putative
AFLA_1372001.12 0.81 −0.47 chitin synthase, putative
AFLA_01369078.45 58.64 −0.42 chitin synthase C
AFLA_09130062.44 50.43 −0.31 chitin biosynthesis protein (Chs5), putative
AFLA_0527805.03 3.14 −0.68 cell wall glucanase (Scw4), putative
AFLA_0966801.53 0.41 −1.89 glucan endo-1,3-alpha-glucosidase agn1 precursor, putative
AFLA_0956801.08 0.34 −1.67 alpha-1,3-glucanase, putative
AFLA_0299506.84 2.55 −1.43 endo-1,3(4)-beta-glucanase, putative
AFLA_045290570.20 213.72 −1.42 extracellular endoglucanase/cellulase, putative
AFLA_1026401.53 0.60 −1.36 exo-beta-1,3-glucanase, putative
AFLA_0533901894.71 823.30 −1.20 GPI-anchored cell wall beta-1,3-endoglucanase EglC
AFLA_0683004467.43 2012.10 −1.15 1,3-beta-glucanosyltransferase Bgt1
AFLA_129440704.92 365.98 −0.95 1,3-beta-glucanosyltransferase, putative
AFLA_0349204.00 2.11 −0.92 glucan endo-1,3-alpha-glucosidase agn1 precursor, putative
AFLA_0584804287.31 2388.65 −0.84 1,3-beta-glucanosyltransferase Gel1
AFLA_08787034.41 20.02 −0.78 Endoglucanase, putative
AFLA_1119703.86 2.26 −0.77 Endoglucanase, putative
AFLA_1264104.24 2.61 −0.70 endoglucanase-1 precursor, putative
AFLA_0527805.03 3.14 −0.68 cell wall glucanase (Scw4), putative
* CK = Control; D1 = 1% dimethylformamide.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pan, L.; Chang, P.; Jin, J.; Yang, Q.; Xing, F. Dimethylformamide Inhibits Fungal Growth and Aflatoxin B1 Biosynthesis in Aspergillus flavus by Down-Regulating Glucose Metabolism and Amino Acid Biosynthesis. Toxins 2020, 12, 683. https://doi.org/10.3390/toxins12110683

AMA Style

Pan L, Chang P, Jin J, Yang Q, Xing F. Dimethylformamide Inhibits Fungal Growth and Aflatoxin B1 Biosynthesis in Aspergillus flavus by Down-Regulating Glucose Metabolism and Amino Acid Biosynthesis. Toxins. 2020; 12(11):683. https://doi.org/10.3390/toxins12110683

Chicago/Turabian Style

Pan, Lin, Peng Chang, Jing Jin, Qingli Yang, and Fuguo Xing. 2020. "Dimethylformamide Inhibits Fungal Growth and Aflatoxin B1 Biosynthesis in Aspergillus flavus by Down-Regulating Glucose Metabolism and Amino Acid Biosynthesis" Toxins 12, no. 11: 683. https://doi.org/10.3390/toxins12110683

APA Style

Pan, L., Chang, P., Jin, J., Yang, Q., & Xing, F. (2020). Dimethylformamide Inhibits Fungal Growth and Aflatoxin B1 Biosynthesis in Aspergillus flavus by Down-Regulating Glucose Metabolism and Amino Acid Biosynthesis. Toxins, 12(11), 683. https://doi.org/10.3390/toxins12110683

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

Article Metrics

Back to TopTop