Aspergillus flavus pangenome (AflaPan) uncovers novel aflatoxin and secondary metabolite associated gene clusters.

BACKGROUND: Aspergillus flavus is an important agricultural and food safety threat due to its production of carcinogenic aflatoxins. It has high level of genetic diversity that is adapted to various environments. Recently, we reported two reference genomes of A. flavus isolates, AF13 (MAT1-2 and highly aflatoxigenic isolate) and NRRL3357 (MAT1-1 and moderate aflatoxin producer). Where, an insertion of 310 kb in AF13 included an aflatoxin producing gene bZIP transcription factor, named atfC. Observations of significant genomic variants between these isolates of contrasting phenotypes prompted an investigation into variation among other agricultural isolates of A. flavus with the goal of discovering novel genes potentially associated with aflatoxin production regulation. Present study was designed with three main objectives: (1) collection of large number of A. flavus isolates from diverse sources including maize plants and field soils; (2) whole genome sequencing of collected isolates and development of a pangenome; and (3) pangenome-wide association study (Pan-GWAS) to identify novel secondary metabolite cluster genes. RESULTS: Pangenome analysis of 346 A. flavus isolates identified a total of 17,855 unique orthologous gene clusters, with mere 41% (7,315) core genes and 59% (10,540) accessory genes indicating accumulation of high genomic diversity during domestication. 5,994 orthologous gene clusters in accessory genome not annotated in either the A. flavus AF13 or NRRL3357 reference genomes. Pan-genome wide association analysis of the genomic variations identified 391 significant associated pan-genes associated with aflatoxin production. Interestingly, most of the significantly associated pan-genes (94%; 369 associations) belonged to accessory genome indicating that genome expansion has resulted in the incorporation of new genes associated with aflatoxin and other secondary metabolites. CONCLUSION: In summary, this study provides complete pangenome framework for the species of Aspergillus flavus along with associated genes for pathogen survival and aflatoxin production. The large accessory genome indicated large genome diversity in the species A. flavus, however AflaPan is a closed pangenome represents optimum diversity of species A. flavus. Most importantly, the newly identified aflatoxin producing gene clusters will be a new source for seeking aflatoxin mitigation strategies and needs new attention in research.

Linkage Mapping and Genome-Wide Association Study Identified Two Peanut Late Leaf Spot Resistance Loci, PLLSR-1 and PLLSR-2, Using Nested Association Mapping.

Identification of candidate genes and molecular markers for late leaf spot (LLS) disease resistance in peanut (Arachis hypogaea) has been a focus of molecular breeding for the U.S. industry-funded peanut genome project. Efforts have been hindered by limited mapping resolution due to low levels of genetic recombination and marker density available in traditional biparental mapping populations. To address this, a multi-parental nested association mapping population has been genotyped with the peanut 58K single-nucleotide polymorphism (SNP) array and phenotyped for LLS severity in the field for 3 years. Joint linkage-based quantitative trait locus (QTL) mapping identified nine QTLs for LLS resistance with significant phenotypic variance explained up to 47.7%. A genome-wide association study identified 13 SNPs consistently associated with LLS resistance. Two genomic regions harboring the consistent QTLs and SNPs were identified from 1,336 to 1,520 kb (184 kb) on chromosome B02 and from 1,026.9 to 1,793.2 kb (767 kb) on chromosome B03, designated as peanut LLS resistance loci, PLLSR-1 and PLLSR-2, respectively. PLLSR-1 contains 10 nucleotide-binding site leucine-rich repeat disease resistance genes. A nucleotide-binding site leucine-rich repeat disease resistance gene, Arahy.VKVT6A, was also identified on homoeologous chromosome A02. PLLSR-2 contains five significant SNPs associated with five different genes encoding callose synthase, pollen defective in guidance protein, pentatricopeptide repeat, acyl-activating enzyme, and C2 GRAM domains-containing protein. This study highlights the power of multi-parent populations such as nested association mapping for genetic mapping and marker-trait association studies in peanuts. Validation of these two LLS resistance loci will be needed for marker-assisted breeding.

Macro-micro exploration on dynamic interaction between aflatoxigenic Aspergillus flavus and maize kernels using Vis/NIR hyperspectral imaging and SEM technology.

Aspergillus flavus and its toxic metabolites-aflatoxins infect and contaminate maize kernels, posing a threat to grain safety and human health. Due to the complexity of microbial growth and metabolic processes, dynamic mechanisms among fungal growth, nutrient depletion of maize kernels and aflatoxin production is still unclear. In this study, visible/near infrared (Vis/NIR) hyperspectral imaging (HSI) combined with the scanning electron microscope (SEM) was used to elucidate the critical organismal interaction at kernel (macro-) and microscopic levels. As kernel damage is the main entrance for fungal invasion, maize kernels with gradually aggravated damages from intact to pierced to halved kernels with A. flavus were cultured for 0-120 h. The spectral fingerprints of the A. flavus-maize kernel complex over time were analyzed with principal components analysis (PCA) of hyperspectral images, where the pseudo-color score maps and the loading plots of the first three PCs were used to investigate the dynamic process of fungal infection and to capture the subtle changes in the complex with different hardness of the maize matrix. The dynamic growth process of A. flavus and the interactions of fungus-maize complexes were explained on a microscopic level using SEM. Specifically, fungus morphology, e.g., hyphae, conidia, and conidiophore (stipe) was accurately captured on the microscopic level, and the interaction process between A. flavus and nutrient loss from the maize kernel tissues (i.e., embryo, and endosperm) was described. Furthermore, the growth stage discrimination models based on PLSDA with the results of CCRC = 100 %, CCRV = 97 %, CCRIV = 93 %, and the prediction models of AFB1 based on PLSR with satisfactory performance (R2C = 0.96, R2V = 0.95, R2IV = 0.93 and RPD = 3.58) were both achieved. In conclusion, the results from both macro-level (Vis/NIR-HSI) and micro-level (SEM) assessments revealed the dynamic organismal interactions in A. flavus-maize kernel complex, and the detailed data could be used for modeling, and quantitative prediction of aflatoxin, which would establish a theoretical foundation for the early detection of fungal or toxin contaminated grains to ensure food security.

Precise exogenous insertion and sequence replacements in poplar by simultaneous HDR overexpression and NHEJ suppression using CRISPR-Cas9.

CRISPR-mediated genome editing has become a powerful tool for the genetic modification of biological traits. However, developing an efficient, site-specific, gene knock-in system based on homology-directed DNA repair (HDR) remains a significant challenge in plants, especially in woody species like poplar. Here, we show that simultaneous inhibition of non-homologous end joining (NHEJ) recombination cofactor XRCC4 and overexpression of HDR enhancer factors CtIP and MRE11 can improve HDR efficiency for gene knock-in. Using this approach, the BleoR gene was integrated onto the 3' end of the MKK2 MAP kinase gene to generate a BleoR-MKK2 fusion protein. Based on fully edited nucleotides evaluated by TaqMan real-time PCR, the HDR-mediated knock-in efficiency was up to 48% when using XRCC4 silencing incorporated with a combination of CtIP and MRE11 overexpression compared with no HDR enhancement or NHEJ silencing. Furthermore, this combination of HDR enhancer overexpression and NHEJ repression also increased genome targeting efficiency and gave 7-fold fewer CRISPR-induced insertions and deletions (InDels), resulting in no functional effects on MKK2-based salt stress responses in poplar. Therefore, this approach may be useful not only in poplar and plants or crops but also in mammals for improving CRISPR-mediated gene knock-in efficiency.

Why Do Plant-Pathogenic Fungi Produce Mycotoxins? Potential Roles for Mycotoxins in the Plant Ecosystem.

For many plant-pathogenic or endophytic fungi, production of mycotoxins, which are toxic to humans, may present a fitness gain. However, associations between mycotoxin production and plant pathogenicity or virulence is inconsistent and difficult due to the complexity of these host-pathogen interactions and the influences of environmental and insect factors. Aflatoxin receives a lot of attention due to its potent toxicity and carcinogenicity but the connection between aflatoxin production and pathogenicity is complicated by the pathogenic ability and prevalence of nonaflatoxigenic isolates in crops. Other toxins directly aid fungi in planta, trichothecenes are important virulence factors, and ergot alkaloids limit herbivory and fungal consumption due to insect toxicity. We review a panel discussion at the American Phytopathological Society's Plant Health 2021 conference, which gathered diverse experts representing different research sectors, career stages, ethnicities, and genders to discuss the diverse roles of mycotoxins in the lifestyles of filamentous fungi of the families Clavicipitaceae, Trichocomaceae (Eurotiales), and Nectriaceae (Hypocreales).

Spatio-temporal patterns of Aspergillus flavus infection and aflatoxin B1 biosynthesis on maize kernels probed by SWIR hyperspectral imaging and synchrotron FTIR microspectroscopy.

The dynamics mechanisms regulating the growth and AFB1 production of Aspergillus flavus during its interactions with maize kernels remain unclear. In this study, shortwave infrared hyperspectral imaging (SWIR-HSI) and synchrotron radiation Fourier transform infrared (SR-FTIR) microspectroscopy were combined to investigate chemical and spatial-temporal changes in incremental damaged maize kernels induced by A. flavus infection at macroscopic and microscopic levels. SWIR-HSI was employed to extract spectral information of A. flavus growth and quantitatively detect AFB1 levels. Satisfactory full-spectrum models and simplified multispectral models were obtained respectively by partial least squares regression (PLSR) for three types of samples. Furthermore, SR-FTIR microspectroscopy coupled with two-dimensional correlation spectroscopy (2DCOS) was utilized to reveal the possible sequence of dynamic changes of nutrient loss and trace AFB1 in maize kernels. It exhibited new insights on how to quantify the spatio-temporal patterns of fungal infection and AFB1 accumulation on maize and provided theoretical basis for online sorting.

Global Transcriptome Profiling Identified Transcription Factors, Biological Process, and Associated Pathways for Pre-Harvest Aflatoxin Contamination in Groundnut.

Pre-harvest aflatoxin contamination (PAC) in groundnut is a serious quality concern globally, and drought stress before harvest further exacerbate its intensity, leading to the deterioration of produce quality. Understanding the host-pathogen interaction and identifying the candidate genes responsible for resistance to PAC will provide insights into the defense mechanism of the groundnut. In this context, about 971.63 million reads have been generated from 16 RNA samples under controlled and Aspergillus flavus infected conditions, from one susceptible and seven resistant genotypes. The RNA-seq analysis identified 45,336 genome-wide transcripts under control and infected conditions. This study identified 57 transcription factor (TF) families with major contributions from 6570 genes coding for bHLH (719), MYB-related (479), NAC (437), FAR1 family protein (320), and a few other families. In the host (groundnut), defense-related genes such as senescence-associated proteins, resveratrol synthase, seed linoleate, pathogenesis-related proteins, peroxidases, glutathione-S-transferases, chalcone synthase, ABA-responsive gene, and chitinases were found to be differentially expressed among resistant genotypes as compared to susceptible genotypes. This study also indicated the vital role of ABA-responsive ABR17, which co-regulates the genes of ABA responsive elements during drought stress, while providing resistance against A. flavus infection. It belongs to the PR-10 class and is also present in several plant-pathogen interactions.

Transcriptome Analysis Identified Coordinated Control of Key Pathways Regulating Cellular Physiology and Metabolism upon Aspergillus flavus Infection Resulting in Reduced Aflatoxin Production in Groundnut.

Aflatoxin-affected groundnut or peanut presents a major global health issue to both commercial and subsistence farming. Therefore, understanding the genetic and molecular mechanisms associated with resistance to aflatoxin production during host-pathogen interactions is crucial for breeding groundnut cultivars with minimal level of aflatoxin contamination. Here, we performed gene expression profiling to better understand the mechanisms involved in reduction and prevention of aflatoxin contamination resulting from Aspergillus flavus infection in groundnut seeds. RNA sequencing (RNA-Seq) of 16 samples from different time points during infection (24 h, 48 h, 72 h and the 7th day after inoculation) in U 4-7-5 (resistant) and JL 24 (susceptible) genotypes yielded 840.5 million raw reads with an average of 52.5 million reads per sample. A total of 1779 unique differentially expressed genes (DEGs) were identified. Furthermore, comprehensive analysis revealed several pathways, such as disease resistance, hormone biosynthetic signaling, flavonoid biosynthesis, reactive oxygen species (ROS) detoxifying, cell wall metabolism and catabolizing and seed germination. We also detected several highly upregulated transcription factors, such as ARF, DBB, MYB, NAC and C2H2 in the resistant genotype in comparison to the susceptible genotype after inoculation. Moreover, RNA-Seq analysis suggested the occurrence of coordinated control of key pathways controlling cellular physiology and metabolism upon A. flavus infection, resulting in reduced aflatoxin production.

Draft Genome Sequences of One Aspergillus parasiticus Isolate and Nine Aspergillus flavus Isolates with Varying Stress Tolerance and Aflatoxin Production.

Aspergillus flavus and Aspergillus parasiticus produce carcinogenic aflatoxins during crop infection, with extensive variations in production among isolates, ranging from atoxigenic to highly toxigenic. Here, we report draft genome sequences of one A. parasiticus isolate and nine A. flavus isolates from field environments for use in comparative, functional, and phylogenetic studies.

Two New Aspergillus flavus Reference Genomes Reveal a Large Insertion Potentially Contributing to Isolate Stress Tolerance and Aflatoxin Production.

Efforts in genome sequencing in the Aspergillus genus have led to the development of quality reference genomes for several important species including A. nidulans, A. fumigatus, and A. oryzae However, less progress has been made for A. flavus As part of the effort of the USDA-ARS Annual Aflatoxin Workshop Fungal Genome Project, the isolate NRRL3357 was sequenced and resulted in a scaffold-level genome released in 2005. Our goal has been biologically driven, focusing on two areas: isolate variation in aflatoxin production and drought stress exacerbating aflatoxin production by A. flavus Therefore, we developed two reference pseudomolecule genome assemblies derived from chromosome arms for two isolates: AF13, a MAT1-2, highly stress tolerant, and highly aflatoxigenic isolate; and NRRL3357, a MAT1-1, less stress tolerant, and moderate aflatoxin producer in comparison to AF13. Here, we report these two reference-grade assemblies for these isolates through a combination of PacBio long-read sequencing and optical mapping, and coupled them with comparative, functional, and phylogenetic analyses. This analysis resulted in the identification of 153 and 45 unique genes in AF13 and NRRL3357, respectively. We also confirmed the presence of a unique 310 Kb insertion in AF13 containing 60 genes. Analysis of this insertion revealed the presence of a bZIP transcription factor, named atfC, which may contribute to isolate pathogenicity and stress tolerance. Phylogenomic analyses comparing these and other available assemblies also suggest that the species complex of A. flavus is polyphyletic.

Functional Biology and Molecular Mechanisms of Host-Pathogen Interactions for Aflatoxin Contamination in Groundnut (Arachis hypogaea L.) and Maize (Zea mays L.).

Aflatoxins are secondary metabolites produced by soilborne saprophytic fungus Aspergillus flavus and closely related species that infect several agricultural commodities including groundnut and maize. The consumption of contaminated commodities adversely affects the health of humans and livestock. Aflatoxin contamination also causes significant economic and financial losses to producers. Research efforts and significant progress have been made in the past three decades to understand the genetic behavior, molecular mechanisms, as well as the detailed biology of host-pathogen interactions. A range of omics approaches have facilitated better understanding of the resistance mechanisms and identified pathways involved during host-pathogen interactions. Most of such studies were however undertaken in groundnut and maize. Current efforts are geared toward harnessing knowledge on host-pathogen interactions and crop resistant factors that control aflatoxin contamination. This study provides a summary of the recent progress made in enhancing the understanding of the functional biology and molecular mechanisms associated with host-pathogen interactions during aflatoxin contamination in groundnut and maize.

Nested-association mapping (NAM)-based genetic dissection uncovers candidate genes for seed and pod weights in peanut (Arachis hypogaea).

Multiparental genetic mapping populations such as nested-association mapping (NAM) have great potential for investigating quantitative traits and associated genomic regions leading to rapid discovery of candidate genes and markers. To demonstrate the utility and power of this approach, two NAM populations, NAM_Tifrunner and NAM_Florida-07, were used for dissecting genetic control of 100-pod weight (PW) and 100-seed weight (SW) in peanut. Two high-density SNP-based genetic maps were constructed with 3341 loci and 2668 loci for NAM_Tifrunner and NAM_Florida-07, respectively. The quantitative trait locus (QTL) analysis identified 12 and 8 major effect QTLs for PW and SW, respectively, in NAM_Tifrunner, and 13 and 11 major effect QTLs for PW and SW, respectively, in NAM_Florida-07. Most of the QTLs associated with PW and SW were mapped on the chromosomes A05, A06, B05 and B06. A genomewide association study (GWAS) analysis identified 19 and 28 highly significant SNP-trait associations (STAs) in NAM_Tifrunner and 11 and 17 STAs in NAM_Florida-07 for PW and SW, respectively. These significant STAs were co-localized, suggesting that PW and SW are co-regulated by several candidate genes identified on chromosomes A05, A06, B05, and B06. This study demonstrates the utility of NAM population for genetic dissection of complex traits and performing high-resolution trait mapping in peanut.

Carbohydrate, glutathione, and polyamine metabolism are central to Aspergillus flavus oxidative stress responses over time.

BACKGROUND: The primary and secondary metabolites of fungi are critical for adaptation to environmental stresses, host pathogenicity, competition with other microbes, and reproductive fitness. Drought-derived reactive oxygen species (ROS) have been shown to stimulate aflatoxin production and regulate in Aspergillus flavus, and may function in signaling with host plants. Here, we have performed global, untargeted metabolomics to better understand the role of aflatoxin production in oxidative stress responses, and also explore isolate-specific oxidative stress responses over time. RESULTS: Two field isolates of A. flavus, AF13 and NRRL3357, possessing high and moderate aflatoxin production, respectively, were cultured in medium with and without supplementation with 15 mM H2O2, and mycelia were collected following 4 and 7 days in culture for global metabolomics. Overall, 389 compounds were described in the analysis which encompassed 9 biological super-pathways and 47 sub-pathways. These metabolites were examined for differential accumulation. Significant differences were observed in both isolates in response to oxidative stress and when comparing sampling time points. CONCLUSIONS: The moderately high aflatoxin-producing isolate, NRRL3357, showed extensive stimulation of antioxidant mechanisms and pathways including polyamines metabolism, glutathione metabolism, TCA cycle, and lipid metabolism while the highly aflatoxigenic isolate, AF13, showed a less vigorous response to stress. Carbohydrate pathway levels also imply that carbohydrate repression and starvation may influence metabolite accumulation at the later timepoint. Higher conidial oxidative stress tolerance and antioxidant capacity in AF13 compared to NRRL3357, inferred from their metabolomic profiles and growth curves over time, may be connected to aflatoxin production capability and aflatoxin-related antioxidant accumulation. The coincidence of several of the detected metabolites in H2O2-stressed A. flavus and drought-stressed hosts also suggests their potential role in the interaction between these organisms and their use as markers/targets to enhance host resistance through biomarker selection or genetic engineering.

Mitigating Aflatoxin Contamination in Groundnut through A Combination of Genetic Resistance and Post-Harvest Management Practices.

Aflatoxin is considered a "hidden poison" due to its slow and adverse effect on various biological pathways in humans, particularly among children, in whom it leads to delayed development, stunted growth, liver damage, and liver cancer. Unfortunately, the unpredictable behavior of the fungus as well as climatic conditions pose serious challenges in precise phenotyping, genetic prediction and genetic improvement, leaving the complete onus of preventing aflatoxin contamination in crops on post-harvest management. Equipping popular crop varieties with genetic resistance to aflatoxin is key to effective lowering of infection in farmer's fields. A combination of genetic resistance for in vitro seed colonization (IVSC), pre-harvest aflatoxin contamination (PAC) and aflatoxin production together with pre- and post-harvest management may provide a sustainable solution to aflatoxin contamination. In this context, modern "omics" approaches, including next-generation genomics technologies, can provide improved and decisive information and genetic solutions. Preventing contamination will not only drastically boost the consumption and trade of the crops and products across nations/regions, but more importantly, stave off deleterious health problems among consumers across the globe.

High-density genetic map using whole-genome resequencing for fine mapping and candidate gene discovery for disease resistance in peanut.

Whole-genome resequencing (WGRS) of mapping populations has facilitated development of high-density genetic maps essential for fine mapping and candidate gene discovery for traits of interest in crop species. Leaf spots, including early leaf spot (ELS) and late leaf spot (LLS), and Tomato spotted wilt virus (TSWV) are devastating diseases in peanut causing significant yield loss. We generated WGRS data on a recombinant inbred line population, developed a SNP-based high-density genetic map, and conducted fine mapping, candidate gene discovery and marker validation for ELS, LLS and TSWV. The first sequence-based high-density map was constructed with 8869 SNPs assigned to 20 linkage groups, representing 20 chromosomes, for the 'T' population (Tifrunner × GT-C20) with a map length of 3120 cM and an average distance of 1.45 cM. The quantitative trait locus (QTL) analysis using high-density genetic map and multiple season phenotyping data identified 35 main-effect QTLs with phenotypic variation explained (PVE) from 6.32% to 47.63%. Among major-effect QTLs mapped, there were two QTLs for ELS on B05 with 47.42% PVE and B03 with 47.38% PVE, two QTLs for LLS on A05 with 47.63% and B03 with 34.03% PVE and one QTL for TSWV on B09 with 40.71% PVE. The epistasis and environment interaction analyses identified significant environmental effects on these traits. The identified QTL regions had disease resistance genes including R-genes and transcription factors. KASP markers were developed for major QTLs and validated in the population and are ready for further deployment in genomics-assisted breeding in peanut.

Proteome analysis of Aspergillus flavus isolate-specific responses to oxidative stress in relationship to aflatoxin production capability.

Aspergillus flavus is an opportunistic pathogen of plants such as maize and peanut under conducive conditions such as drought stress resulting in significant aflatoxin production. Drought-associated oxidative stress also exacerbates aflatoxin production by A. flavus. The objectives of this study were to use proteomics to provide insights into the pathogen responses to H2O2-derived oxidative stress, and to identify potential biomarkers and targets for host resistance breeding. Three isolates, AF13, NRRL3357, and K54A with high, moderate, and no aflatoxin production, were cultured in medium supplemented with varying levels of H2O2, and examined using an iTRAQ (Isobaric Tags for Relative and Absolute Quantification) approach. Overall, 1,173 proteins were identified and 220 were differentially expressed (DEPs). Observed DEPs encompassed metabolic pathways including antioxidants, carbohydrates, pathogenicity, and secondary metabolism. Increased lytic enzyme, secondary metabolite, and developmental pathway expression in AF13 was correlated with oxidative stress tolerance, likely assisting in plant infection and microbial competition. Elevated expression of energy and cellular component production in NRRL3357 and K54A implies a focus on oxidative damage remediation. These trends explain isolate-to-isolate variation in oxidative stress tolerance and provide insights into mechanisms relevant to host plant interactions under drought stress allowing for more targeted efforts in host resistance research.

Deciphering drought-induced metabolic responses and regulation in developing maize kernels.

Drought stress conditions decrease maize growth and yield, and aggravate preharvest aflatoxin contamination. While several studies have been performed on mature kernels responding to drought stress, the metabolic profiles of developing kernels are not as well characterized, particularly in germplasm with contrasting resistance to both drought and mycotoxin contamination. Here, following screening for drought tolerance, a drought-sensitive line, B73, and a drought-tolerant line, Lo964, were selected and stressed beginning at 14 days after pollination. Developing kernels were sampled 7 and 14 days after drought induction (DAI) from both stressed and irrigated plants. Comparative biochemical and metabolomic analyses profiled 409 differentially accumulated metabolites. Multivariate statistics and pathway analyses showed that drought stress induced an accumulation of simple sugars and polyunsaturated fatty acids and a decrease in amines, polyamines and dipeptides in B73. Conversely, sphingolipid, sterol, phenylpropanoid and dipeptide metabolites accumulated in Lo964 under drought stress. Drought stress also resulted in the greater accumulation of reactive oxygen species (ROS) and aflatoxin in kernels of B73 in comparison with Lo964 implying a correlation in their production. Overall, field drought treatments disordered a cascade of normal metabolic programming during development of maize kernels and subsequently caused oxidative stress. The glutathione and urea cycles along with the metabolism of carbohydrates and lipids for osmoprotection, membrane maintenance and antioxidant protection were central among the drought stress responses observed in developing kernels. These results also provide novel targets to enhance host drought tolerance and disease resistance through the use of biotechnologies such as transgenics and genome editing.

Oxidative stress and carbon metabolism influence Aspergillus flavus transcriptome composition and secondary metabolite production.

Contamination of crops with aflatoxin is a serious global threat to food safety. Aflatoxin production by Aspergillus flavus is exacerbated by drought stress in the field and by oxidative stress in vitro. We examined transcriptomes of three toxigenic and three atoxigenic isolates of A. flavus in aflatoxin conducive and non-conducive media with varying levels of H2O2 to investigate the relationship of secondary metabolite production, carbon source, and oxidative stress. We found that toxigenic and atoxigenic isolates employ distinct mechanisms to remediate oxidative damage, and that carbon source affected the isolates' expression profiles. Iron metabolism, monooxygenases, and secondary metabolism appeared to participate in isolate oxidative responses. The results suggest that aflatoxin and aflatrem biosynthesis may remediate oxidative stress by consuming excess oxygen and that kojic acid production may limit iron-mediated, non-enzymatic generation of reactive oxygen species. Together, secondary metabolite production may enhance A. flavus stress tolerance, and may be reduced by enhancing host plant tissue antioxidant capacity though genetic improvement by breeding selection.

Responses of Aspergillus flavus to Oxidative Stress Are Related to Fungal Development Regulator, Antioxidant Enzyme, and Secondary Metabolite Biosynthetic Gene Expression.

The infection of maize and peanut with Aspergillus flavus and subsequent contamination with aflatoxin pose a threat to global food safety and human health, and is exacerbated by drought stress. Drought stress-responding compounds such as reactive oxygen species (ROS) are associated with fungal stress responsive signaling and secondary metabolite production, and can stimulate the production of aflatoxin by A. flavus in vitro. These secondary metabolites have been shown to possess diverse functions in soil-borne fungi including antibiosis, competitive inhibition of other microbes, and abiotic stress alleviation. Previously, we observed that isolates of A. flavus showed differences in oxidative stress tolerance which correlated with their aflatoxin production capabilities. In order to better understand these isolate-specific oxidative stress responses, we examined the transcriptional responses of field isolates of A. flavus with varying levels of aflatoxin production (NRRL3357, AF13, and Tox4) to H2O2-induced oxidative stress using an RNA sequencing approach. These isolates were cultured in an aflatoxin-production conducive medium amended with various levels of H2O2. Whole transcriptomes were sequenced using an Illumina HiSeq platform with an average of 40.43 million filtered paired-end reads generated for each sample. The obtained transcriptomes were then used for differential expression, gene ontology, pathway, and co-expression analyses. Isolates which produced higher levels of aflatoxin tended to exhibit fewer differentially expressed genes than isolates with lower levels of production. Genes found to be differentially expressed in response to increasing oxidative stress included antioxidant enzymes, primary metabolism components, antibiosis-related genes, and secondary metabolite biosynthetic components specifically for aflatoxin, aflatrem, and kojic acid. The expression of fungal development-related genes including aminobenzoate degradation genes and conidiation regulators were found to be regulated in response to increasing stress. Aflatoxin biosynthetic genes and antioxidant enzyme genes were also found to be co-expressed and highly correlated with fungal biomass under stress. This suggests that these secondary metabolites may be produced as part of coordinated oxidative stress responses in A. flavus along with antioxidant enzyme gene expression and developmental regulation.