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.

Aspergillus flavus La3279, a component strain of the Aflasafe™ biocontrol product, contains a partial aflatoxin biosynthesis gene cluster followed by a genomic region highly variable among A. flavus isolates.

Aspergillus flavus communities in agricultural fields consist of isolates with varying abilities to produce aflatoxins, which are highly toxic and carcinogenic to humans and animals. Biological control using multiple non-aflatoxigenic strains as a formulation to outcompete aflatoxigenic A. flavus has become a mainstream strategy. Aflasafe™ is a biocontrol product composed of four strains, Ka16127, La3279, La3304 and Og0222. It was first developed in Nigeria and is now widely used on maize and groundnut. In this study, phylogenetic analyses based on genome-wide single nucleotide polymorphisms showed that Ka16127 and La3304 were more closely related to each other than both were to La3279, and the three were distantly related to Og0222. Detailed molecular characterization of La3279 indicated that its genome, contradictory to the published report, lacked approximately half of the aflatoxin gene cluster as well as the entire cyclopiazonic acid gene cluster. La3279 was a member of the previously known "pattern E" group, which includes A. flavus and Aspergillus oryzae isolates that have the aforementioned deletion followed by a 3.8-kb "E block" sequence insertion. In comparison to the E block, corresponding regions in typical aflatoxigenic S-morphotype/genotype isolates as well as Ka16127 and La3304 were found to lack 1.1 kb of the 5' portion whereas L-morphotype/genotype isolates contained a complete nonhomologous region characterized by 2.5 copies of A. flavus telomeric repeat sequence at one end. Regions corresponding to the E block were highly variable and were useful for classifying A. flavus isolates into groups that mostly contained both mating types. The presence of both mating-type genes in genetically closely related A. flavus suggests a previously active sexual cycle. It could facilitate the development of a refined biocontrol strategy such as deploying biocontrol strains with the same mating-type that is predominant in a field A. flavus population.

Authentication of Aspergillus parasiticus strains in the genome database of the National Center for Biotechnology Information.

OBJECTIVE: The use of genome sequences from strains authenticated to correct species level is a prerequisite for confidently exploring the evolutionary relationship among related species. Aspergillus strains erroneously curated as Aspergillus oryzae and Aspergillus fumigatus have been noticed in the National Center for Biotechnology Information (NCBI) genome database. Aspergillus parasiticus is one of several aspergilli that produce aflatoxin, the most potent carcinogenic mycotoxin known up to now. To ensure that valid conclusions are drawn by researchers from their genomics-related studies, molecular analyses were carried out to authenticate identities of A. parasiticus strains in the NCBI genome database. RESULTS: Two of the nine supposedly A. parasiticus strains, E1365 and NRRL2999, were found to be misidentified. They turned out to be Aspergillus flavus based on genome-wide single nucleotide polymorphisms (SNPs) and genetic features associated with production of aflatoxin and cyclopiazonic acid. NRRL2999 lacked the additional partial aflatoxin gene cluster known to be present in its equivalent strain, designated as SU-1, and shared a very low total SNPs count specifically with A. flavus NRRL3357 but not with other A. flavus isolates. Therefore, the mislabeled NRRL2999 strain actually is a clonal strain of A. flavus NRRL3357, whose genome was first sequenced in 2005.

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.

Identification of a copper-transporting ATPase involved in biosynthesis of A. flavus conidial pigment.

Conidia are asexual spores and play a crucial role in fungal dissemination. Conidial pigmentation is important for tolerance against UV radiation and contributes to survival of fungi. The molecular basis of conidial pigmentation has been studied in several fungal species. In spite of sharing the initial common step of polyketide formation, other steps for pigment biosynthesis appear to be species-dependent. In this study, we isolated an Aspergillus flavus spontaneous mutant that produced yellow conidia. The underlying genetic defect, a three-nucleotide in-frame deletion in the gene, AFLA_051390, that encodes a copper-transporting ATPase, was identified by a comparative genomics approach. This genetic association was confirmed by disruption of the wild-type gene. When yellow mutants were grown on medium supplemented with copper ions or chloride ions, green conidial color was partially and nearly completely restored, respectively. Further disruption of AFLA_045660, an orthologue of Aspergillus nidulans yA (yellow pigment) that encodes a multicopper oxidase, in wild type and a derived strain producing dark green conidia showed that it yielded mutants that produced gold conidia. The results placed formation of the gold pigment after that of the yellow pigment and before that of the dark green pigment. Using reported inhibitors of DHN-melanin (tricyclazole and phthalide) and DOPA-melanin (tropolone and kojic acid) pathways on a set of conidial color mutants, we investigated the involvement of melanin biosynthesis in A. flavus conidial pigment formation. Results imply that both pathways have no bearing on conidial pigment biosynthesis of A. flavus.

Transcriptional Regulation of Aflatoxin Biosynthesis and Conidiation in Aspergillus flavus by Wickerhamomyces anomalus WRL-076 for Reduction of Aflatoxin Contamination.

Aspergillus flavus is a ubiquitous saprophytic fungus found in soils across the world. The fungus is the major producer of aflatoxin (AF) B1, which is toxic and a potent carcinogen to humans. Aflatoxin B1 (AFB1) is often detected in agricultural crops such as corn, peanut, almond, and pistachio. It is a serious and recurrent problem and causes substantial economic losses. Wickerhamomyces anomalus WRL-076 was identified as an effective biocontrol yeast against A. flavus. In this study, the associated molecular mechanisms of biocontrol were investigated. We found that the expression levels of eight genes, aflR, aflJ, norA, omtA, omtB, pksA, vbs, and ver-1 in the aflatoxin biosynthetic pathway cluster were suppressed. The decreases ranged from several to 10,000 fold in fungal samples co-cultured with W. anomalus. Expression levels of conidiation regulatory genes brlA, abaA, and wetA as well as sclerotial regulatory gene (sclR) were all down regulated. Consistent with the decreased gene expression levels, aflatoxin concentrations in cultural medium were reduced to barely detectable. Furthermore, fungal biomass and conidial number were significantly reduced by 60% and more than 95%, respectively. The results validate the biocontrol efficacy of W. anomalus WRL-076 observed in the field experiments.

Monitoring Metabolite Production of Aflatoxin Biosynthesis by Orbitrap Fusion Mass Spectrometry and a D-Optimal Mixture Design Method.

Aflatoxins, highly toxic and carcinogenic to humans, are synthesized via multiple intermediates by a complex pathway in several Aspergilli, including Aspergillus flavus. Few analytical methods are available for monitoring the changes in metabolite profiles of the aflatoxin biosynthesis pathway under different growth and environmental conditions. In the present study, we developed by a D-optimal mixture design a solvent system, methanol/dichloromethane/ethyl acetate/formic acid (0.36/0.31/0.32/0.01), that was suitable for extracting the pathway metabolites. The matrix effect from dilution of cell extracts was negligible. To facilitate the identification of these metabolites, we constructed a fragmentation ion library. We further employed liquid chromatography coupled with high-resolution mass spectroscopy (UHPLC-HRMS) for simultaneous quantification of the metabolites. The limit of detection (LOD) and limit of quantitation (LOQ) were 0.002-0.016 and 0.008-0.05 µg/kg, respectively. The spiked recovery rates ranged from 81.3 to 100.3% with intraday and interday precision less than 7.6%. Using the method developed to investigate the time-course aflatoxin biosynthesis, we found that precursors, including several possible toxins (with a carcinogenic group similar to aflatoxin B1), occurred together with aflatoxin, and that production increased rapidly at the early growth stage, peaked on day four, and then decreased substantially. The maximum production of aflatoxin B1 and aflatoxin B2 occurred 1 day later. Moreover, the dominant branch pathway was the one for aflatoxin B1 formation. We revealed that the antiaflatoxigenicity mechanism of Leclercia adecarboxylata WT16 was associated with a factor upstream of the aflatoxin biosynthesis pathway. The design strategies can be applied to characterize or detect other secondary metabolites to provide a snapshot of the dynamic changes during their biosynthesis.

Cyclopiazonic Acid Is a Pathogenicity Factor for Aspergillus flavus and a Promising Target for Screening Germplasm for Ear Rot Resistance.

Aspergillus flavus, an opportunistic pathogen, contaminates maize and other key crops with carcinogenic aflatoxins (AFs). Besides AFs, A. flavus makes many more secondary metabolites (SMs) whose toxicity in insects or vertebrates has been studied. However, the role of SMs in the invasion of plant hosts by A. flavus remains to be investigated. Cyclopiazonic acid (CPA), a neurotoxic SM made by A. flavus, is a nanomolar inhibitor of endoplasmic reticulum calcium ATPases (ECAs) and a potent inducer of cell death in plants. We hypothesized that CPA, by virtue of its cytotoxicity, may serve as a key pathogenicity factor that kills plant cells and supports the saprophytic life style of the fungus while compromising the host defense response. This proposal was tested by two complementary approaches. A comparison of CPA levels among A. flavus isolates indicated that CPA may be a determinant of niche adaptation, i.e., isolates that colonize maize make more CPA than those restricted only to the soil. Further, mutants in the CPA biosynthetic pathway are less virulent in causing ear rot than their wild-type parent in field inoculation assays. Additionally, genes encoding ECAs are expressed in developing maize seeds and are induced by A. flavus infection. Building on these results, we developed a seedling assay in which maize roots were exposed to CPA, and cell death was measured as Evans Blue uptake. Among >40 maize inbreds screened for CPA tolerance, inbreds with proven susceptibility to ear rot were also highly CPA sensitive. The publicly available data on resistance to silk colonization or AF contamination for many of the lines was also broadly correlated with their CPA sensitivity. In summary, our studies show that i) CPA serves as a key pathogenicity factor that enables the saprophytic life style of A. flavus and ii) maize inbreds are diverse in their tolerance to CPA. Taking advantage of this natural variation, we are currently pursuing both genome-wide and candidate gene approaches to identify novel components of maize resistance to Aspergillus ear rot.

Development of an Enzyme-Linked Immunosorbent Assay Method Specific for the Detection of G-Group Aflatoxins.

To detect and monitor G-group aflatoxins in agricultural products, we generated class-specific monoclonal antibodies that specifically recognized aflatoxins G1 and G2. Of the final three positive and stable hybridomas obtained, clone 2G6 produced a monoclonal antibody that had equal sensitivity to aflatoxins G1 and G2, and did not cross-react with aflatoxins B1, B2, or M1. Its IC50 values for aflatoxins G1 and G2 were 17.18 ng·mL(-1) and 19.75 ng·mL(-1), respectively. Using this new monoclonal antibody, we developed a competitive indirect enzyme-linked immunosorbent assay (CI-ELISA); the method had a limit of detection of 0.06 ng·mL(-1). To validate this CI-ELISA, we spiked uncontaminated peanut samples with various amounts of aflatoxins G1 and G2 and compared recovery rates with those determined by a standard HPLC method. The recovery rates of the CI-ELISA ranging from 94% to 103% were comparable to those of the HPLC (92% to 102%). We also used both methods to determine the amounts of G-group aflatoxins in five peanut samples contaminated by aflatoxin B1-positive, and their relative standard deviations ranged from 8.4% to 17.7% (under 20%), which demonstrates a good correlation between the two methods. We further used this CI-ELISA to assess the ability of 126 fungal strains isolated from peanuts or field soils to produce G-group aflatoxins. Among these, seven stains producing different amounts of G-group aflatoxins were identified. Our results showed that the monoclonal antibody 2 G6-based CI-ELISA was suitable for the detection of G-group aflatoxins present in peanuts and also those produced by fungi.

Transcriptomic profiles of Aspergillus flavus CA42, a strain that produces small sclerotia, by decanal treatment and after recovery.

Aspergillus flavus is a ubiquitous saprophyte and is capable of producing many secondary metabolites including the carcinogenic aflatoxins. The A. flavus population that produces small sclerotia (S strain) has been implicated as the culprit for persistent aflatoxin contamination in field crops. We investigated how the plant volatile decanal, a C10 fatty aldehyde, affected the growth and development of the S strain A. flavus. Decanal treatment yielded fluffy variants lacking sclerotia and conidia and exhibiting a dosage-dependent radial colony growth. We used RNA-Seq analysis to examine transcriptomic changes caused by decanal and after removal of decanal. Mature sclerotia contained only 80% of the total transcripts detected in all samples in comparison to 94% for the decanal treated culture. Gene ontology (GO) analysis showed that decanal treatment increased expression of genes involved in oxidoreductase activity, cellular carbohydrate metabolism, alcohol metabolism and aflatoxin biosynthesis. The treatment affected cellular components associated with cell wall, and gene expression of glucanases, α-amylases, pectinesterase and peptidase required for its biosynthesis was increased. After decanal was removed, the culture resumed sclerotial production. Moreover, its GO terms significantly overlapped with those of the untreated culture; five of the enriched molecular functions, oxidoreductase activity, monooxygenase activity, electron carrier activity, heme binding, and iron binding were found in the untreated culture. The GO term of cellular component enriched was mainly integral protein constituents of the membrane. The results suggested that decanal halted development at the vegetative state rendering the fungus unable to produce conidia and sclerotia. The induced fluffy phenotype could be related to lower transcript abundance of flbB, flbD, and flbE but not to veA expression. Increased abundance of the laeA transcript in the treated culture correlated with early transcriptional activation of aflatoxin and kojic acid biosynthesis gene clusters. Expression profiles revealed subtle differences in timing of activation of the respective 55 secondary metabolite gene clusters.

Genome-wide analysis of the Zn(II)2Cys6 zinc cluster-encoding gene family in Aspergillus flavus.

Proteins with a Zn(II)2Cys6 domain, Cys-X2-Cys-X6-Cys-X5₋12-Cys-X2-Cys-X6₋9-Cys (hereafter, referred to as the C6 domain), form a subclass of zinc finger proteins found exclusively in fungi and yeast. Genome sequence databases of Saccharomyces cerevisiae and Candida albicans have provided an overview of this family of genes. Annotation of this gene family in most fungal genomes is still far from perfect and refined bioinformatic algorithms are urgently needed. Aspergillus flavus is a saprophytic soil fungus that can produce the carcinogenic aflatoxin. It is the second leading causative agent of invasive aspergillosis. The 37-Mb genome of A. flavus is predicted to encode 12,000 proteins. Two and a half percent of the total proteins are estimated to contain the C6 domain, more than twofold greater than those estimated for yeast, which is about 1 %. The variability in the spacing between cysteines, C3-C4 and C5-C6, in the zinc cluster enables classification of the domains into distinct subgroups, which are also well conserved in Aspergillus nidulans. Sixty-six percent (202/306) of the A. flavus C6 proteins contain a specific transcription factor domain, and 7 % contain a domain of unknown function, DUF3468. Two A. nidulans C6 proteins containing the DUF3468 are involved in asexual conidiation and another two in sexual differentiation. In the anamorphic A. flavus, a homolog of the latter lacks the C6 domain. A. flavus being heterothallic and reproducing mainly through conidiation appears to have lost some components involved in homothallic sexual development. Of the 55 predicted gene clusters thought to be involved in production of secondary metabolites, only about half have a C6-encoding gene in or near the gene clusters. The features revealed by the A. flavus C6 proteins likely are common for other ascomycete fungi.

Characterization of aflatoxigenic and non-aflatoxigenic Aspergillus flavus isolates from pistachio.

Pistachio is a popular snack food. Aflatoxin contamination of pistachio nuts is a serious problem for many producing countries. The development of biological control methods based on ecological parameters is an environmentally friendly approach. Thirty-eight Aspergillus flavus isolates collected from a pistachio orchard in California (CA) were analyzed for production of aflatoxin (AF), cyclopiazonic acid (CPA), vegetative compatibility groups (VCGs), and mating types. All aflatoxigenic isolates produced both AFB1 and CPA. The most toxigenic one was CA28 which produced 164 µg AFB1 per 5 ml PDA fungal culture and small sclerotia (S strain, sclertoium size less than 400 µm). The other aflatoxigenic strains produce AFB1 ranging from 1.2 µg to 80 µg per 5 ml fungal culture. Twenty-one percent of the CA isolates produced AFB1, 84% produced CPA and half formed sclerotia on at least one of three tested media. The 38 CA isolates formed 26 VCGs, 6 of which had two or more isolates and 20 contained single isolates. The S strain isolates belong to 4 different VCGs. Genomic profiling by a retrotransposon DNA probe revealed fingerprint patterns that were highly polymorphic. The predicted VCGs (Pred-VCGs) based on a similarity coefficient >80% matched the VCGs of multiple isolates determined by complementation. All isolates within a VCG had the same mating-type gene of either MAT1-1 or MAT1-2. Uncorrected and VCG-corrected MAT1-1 and MAT1-2 among the isolates were equally distributed.

Loss of msnA, a putative stress regulatory gene, in Aspergillus parasiticus and Aspergillus flavus increased production of conidia, aflatoxins and kojic acid.

Production of the harmful carcinogenic aflatoxins by Aspergillus parasiticus and Aspergillus flavus has been postulated to be a mechanism to relieve oxidative stress. The msnA gene of A. parasiticus and A. flavus is the ortholog of Saccharomyces cerevisiae MSN2 that is associated with multi-stress response. Compared to wild type strains, the msnA deletion (∆msnA) strains of A. parasiticus and A. flavus exhibited retarded colony growth with increased conidiation. The ∆msnA strains also produced slightly higher amounts of aflatoxins and elevated amounts of kojic acid on mixed cereal medium. Microarray assays showed that expression of genes encoding oxidative stress defense enzymes, i.e., superoxide dismutase, catalase, and cytochrome c peroxidase in A. parasiticus ∆msnA, and the catalase A gene in A. flavus ∆msnA, was up-regulated. Both A. parasiticus and A. flavus ∆msnA strains produced higher levels of reactive oxygen species (ROS), and ROS production of A. flavus msnA addback strains was decreased to levels comparable to that of the wild type A. flavus. The msnA gene appears to be required for the maintenance of the normal oxidative state. The impairment of msnA resulted in the aforementioned changes, which might be used to combat the increased oxidative stress in the cells.

Absence of the aflatoxin biosynthesis gene, norA, allows accumulation of deoxyaflatoxin B1 in Aspergillus flavus cultures.

Biosynthesis of the highly toxic and carcinogenic aflatoxins in select Aspergillus species from the common intermediate O-methylsterigmatocystin has been postulated to require only the cytochrome P450 monooxygenase, OrdA (AflQ). We now provide evidence that the aryl alcohol dehydrogenase NorA (AflE) encoded by the aflatoxin biosynthetic gene cluster in Aspergillus flavus affects the accumulation of aflatoxins in the final steps of aflatoxin biosynthesis. Mutants with inactive norA produced reduced quantities of aflatoxin B(1) (AFB(1)), but elevated quantities of a new metabolite, deoxyAFB(1). To explain this result, we suggest that, in the absence of NorA, the AFB(1) reduction product, aflatoxicol, is produced and is readily dehydrated to deoxyAFB(1) in the acidic medium, enabling us to observe this otherwise minor toxin produced in wild-type A. flavus.

Clustered genes involved in cyclopiazonic acid production are next to the aflatoxin biosynthesis gene cluster in Aspergillus flavus.

Cyclopiazonic acid (CPA), an indole-tetramic acid mycotoxin, is produced by many species of Aspergillus and Penicillium. In addition to CPA Aspergillus flavus produces polyketide-derived carcinogenic aflatoxins. Aflatoxin biosynthesis genes form a gene cluster in a subtelomeric region. Isolates of A. flavus lacking aflatoxin production due to the loss of the entire aflatoxin gene cluster and portions of the subtelomeric region are often unable to produce CPA, which suggests a physical link of genes involved in CPA biosynthesis to the aflatoxin gene cluster. Examining the subtelomeric region in A. flavus isolates of different chemotypes revealed a region possibly associated with CPA production. Disruption of three of the four genes present in this region predicted to encode a monoamine oxidase, a dimethylallyl tryptophan synthase, and a hybrid polyketide non-ribosomal peptide synthase abolished CPA production in an aflatoxigenic A. flavus strain. Therefore, some of the CPA biosynthesis genes are organized in a mini-gene cluster that is next to the aflatoxin gene cluster in A. flavus.

Are the genes nadA and norB involved in formation of aflatoxin G(1)?

Aflatoxins, the most toxic and carcinogenic family of fungal secondary metabolites, are frequent contaminants of foods intended for human consumption. Previous studies showed that formation of G-group aflatoxins (AFs) from O-methylsterigmatocystin (OMST) by certain Aspergillus species involves oxidation by the cytochrome P450 monooxygenases, OrdA (AflQ) and CypA (AflU). However, some of the steps in the conversion have not yet been fully defined. Extracts of Aspergillus parasiticus disruption mutants of the OYE-FMN binding domain reductase-encoding gene nadA (aflY) contained a 386 Da AFG(1) precursor. A compound with this mass was predicted as the product of sequential OrdA and CypA oxidation of OMST. Increased amounts of a 362 Da alcohol, the presumptive product of NadA reduction, accumulate in extracts of fungi with disrupted aryl alcohol dehydrogenase-encoding gene norB. These results show that biosynthesis of AFG(1) involves NadA reduction and NorB oxidation.

Genes differentially expressed by Aspergillus flavus strains after loss of aflatoxin production by serial transfers.

Aflatoxins are carcinogenic fungal secondary metabolites produced by Aspergillus flavus and other closely related species. Levels of aflatoxins in agricultural commodities are stringently regulated by many countries because of the health hazard, and thus, aflatoxins are of major concern to both producers and consumers. A cluster of genes responsible for aflatoxin biosynthesis has been identified; however, expression of these genes is a complex and poorly understood phenomenon. To better understand the molecular events that are associated with aflatoxin production, three separate nonaflatoxigenic A. flavus strains were produced through serial transfers of aflatoxigenic parental strains. The three independent aflatoxigenic/nonaflatoxigenic pairs were compared via transcription profiling by microarray analyses. Cross comparisons identified 22 features in common between the aflatoxigenic/nonaflatoxigenic pairs. Physical mapping of the 22 features using the Aspergillus oryzae genome sequence for reference identified 16 unique genes. Aflatoxin biosynthetic and regulatory gene expression levels were not significantly different between the aflatoxigenic/nonaflatoxigenic pairs, which suggests that the inability to produce aflatoxins is not due to decreased expression of known biosynthetic or regulatory genes. Of the 16 in common genes, only one gene homologous to glutathione S-transferase genes showed higher expression in the nonaflatoxigenic progeny relative to the parental strains. This gene, named hcc, was selected for over-expression in an aflatoxigenic A. flavus strain to determine if it was directly responsible for loss of aflatoxin production. Although hcc transformants showed six- to ninefold increase in expression, no discernible changes in colony morphology or aflatoxin production were detected. Possible roles of hcc and other identified genes are discussed in relation to regulation of aflatoxin biosynthesis.

Characterization of AFLAV, a Tf1/Sushi retrotransposon from Aspergillus flavus.

The plasmid, pAF28, a genomic clone from Aspergillus flavus NRRL 6541, has been used as a hybridization probe to fingerprint A. flavus strains isolated in corn and peanut fields. The insert of pAF28 contains a 4.5 kb region which encodes a truncated retrotransposon (AfRTL-1). In search for a full-length and intact copy of retrotransposon, we exploited a novel PCR cloning strategy by amplifying a 3.4 kb region from the genomic DNA of A. flavus NRRL 6541. The fragment was cloned into pCR 4-TOPO. Sequence analysis confirmed that this region encoded putative domains of partial reverse transcriptase, RNase H, and integrase of the predicted retrotransposon. The two flanking long terminal repeats (LTRs) and the sequence between them comprise a putative full-length LTR retrotransposon of 7799 bp in length. This intact retrotransposon sequence is named AFLAV (A. flavus Retrotransposon). The order of the predicted catalytic domains in the polyprotein (Pol) placed AFLAV in the Tf1/sushi subgroup of the Ty3/gypsy retrotransposon family. Primers derived from AFLAV sequence were used to screen this retrotransposon in other strains of A. flavus. More than fifty strains of A. flavus isolated from different geological origins were surveyed and the results show that many strains have extensive deletions in the regions encoding the capsid (Gag) and Pol.

Nonaflatoxigenic Aspergillus flavus TX9-8 competitively prevents aflatoxin accumulation by A. flavus isolates of large and small sclerotial morphotypes.

Aflatoxins are a family of highly toxic and carcinogenic toxins produced by several Aspergillus species. Aflatoxin contamination of agricultural commodities both pre- and postharvest is a serious food safety issue and a significant economic concern. Using nonaflatoxigenic A. flavus isolates to competitively exclude toxigenic A. flavus isolates in agricultural fields has become an adopted approach to reduce aflatoxin contamination. From screening subgroups of nonaflatoxigenic A. flavus, we identified an A. flavus isolate, TX9-8, which competed well with three A. flavus isolates producing low, intermediate, and high levels of aflatoxins, respectively. TX9-8 has a defective polyketide synthase gene (pksA), which is necessary for aflatoxin biosynthesis. Co-inoculating TX9-8 at the same time with large sclerotial (L strain) A. flavus isolates at a ratio of 1:1 or 1:10 (TX9-8:toxigenic) prevented aflatoxin accumulation. The intervention of TX9-8 on small sclerotial (S strain) A. flavus isolates varied and depended on isolate and ratio of co-inoculation. At a ratio of 1:1 TX9-8 prevented aflatoxin accumulation by A. flavus CA28 and reduced aflatoxin accumulation 10-fold by A. flavus CA43. No decrease in aflatoxin accumulation was apparent when TX9-8 was inoculated 24 h after toxigenic L- or S strain A. flavus isolates started growing. The competitive effect likely is due to TX9-8 outgrowing toxigenic A. flavus isolates.

Cladal relatedness among Aspergillus oryzae isolates and Aspergillus flavus S and L morphotype isolates.

Aspergillus flavus is the main etiological agent for aflatoxin contamination of crops. Its close relative, A. oryzae, does not produce aflatoxins and has been widely used to produce fermented foods. We compared the phylogeny of A. oryzae isolates and L- and S-type sclerotial isolates of A. flavus using single nucleotide polymorphisms in the omtA gene in the aflatoxin biosynthesis gene cluster and deletions in and distal to the norB-cypA intergenic region as phylogenetic signals. Aflatoxin-producing ability and sclerotial size also were weighted in the analysis. Like A. flavus, the A. oryzae isolates form a polyphyletic assemblage. A. oryzae isolates in one clade strikingly resemble an A. flavus subgroup of atoxigenic L-type isolates. All toxigenic S-type isolates closely resemble another subgroup of atoxigenic L-type isolates. Because atoxigenic S-type isolates are extremely rare, we hypothesize that loss of aflatoxin production in S-type isolates may occur concomitantly with a change to L-type sclerotia. All toxigenic L-type isolates, unlike A. oryzae, have a 1.0 kb deletion in the norB-cypA region. Although A. oryzae isolates, like S-type, have a 1.5 kb deletion in the norB-cypA region, none were cladally related to S-type A. flavus isolates. Our results show that A. flavus populations are genetically diverse. A. oryzae isolates may descend from certain atoxigenic L-type A. flavus isolates.