
==== Front
Commun Biol
Commun Biol
Communications Biology
2399-3642
Nature Publishing Group UK London

6738
10.1038/s42003-024-06738-w
Article
Comparative analysis of the genomes and aflatoxin production patterns of three species within the Aspergillus section Flavi reveals an undescribed chemotype and habitat-specific genetic traits
Schamann Alexandra
Soukup Sebastian T.
Geisen Rolf
Kulling Sabine
http://orcid.org/0000-0003-2844-9873
Schmidt-Heydt Markus Markus.Schmidt-Heydt@mri.bund.de

https://ror.org/045gmmg53 grid.72925.3b 0000 0001 1017 8329 Department of Safety and Quality of Fruit and Vegetables, Max Rubner-Institut, Federal Research Institute of Nutrition and Food, Karlsruhe, Germany
13 9 2024
13 9 2024
2024
7 113418 12 2023
14 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Aflatoxins are the most dangerous mycotoxins for food safety. They are mainly produced by Aspergillus flavus, A. parasiticus, and A. minisclerotigenes. The latter, an understudied species, was the main culprit for outbreaks of fatal aflatoxicosis in Kenya in the past. To determine specific genetic characteristics of these Aspergillus species, their genomes are comparatively analyzed. Differences reflecting the typical habitat are reported, such as an increased number of carbohydrate-active enzymes, including enzymes for lignin degradation, in the genomes of A. minisclerotigenes and A. parasiticus. Further, variations within the aflatoxin gene clusters are described, which are related to different chemotypes of aflatoxin biosynthesis. These include a substitution within the aflL gene of the A. parasiticus isolate, which leads to the translation of a stop codon, thereby switching off the production of the group 1 aflatoxins B1 and G1. In addition, we demonstrate that the inability of the A. minisclerotigenes isolates to produce group G aflatoxins is associated with a 2.2 kb deletion within the aflF and aflU genes. These findings reveal a relatively high genetic homology among the three Aspergillus species investigated. However, they also demonstrate consequential genetic differences that have an important impact on risk-assessment and food safety.

Aspergillus minisclerotigenes, linked to fatal aflatoxicosis in Kenya, shares high genetic similarity with A. flavus and A. parasiticus. However, genome analysis reveals key differences in aflatoxin gene clusters and enzymes, crucial for food safety and risk assessment.

Subject terms

Fungal genomics
Food microbiology
Federal Office for Agriculture and Food (BLE)issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Aspergillus subgenus Circumdati section Flavi comprises many species that are relevant for food and nutrition in contrastive ways. On the one hand, certain species, such as Aspergillus oryzae, are utilized in the production of fermented foods like koji and miso1. On the other hand, many representatives produce mycotoxins and thus pose a threat to food safety. At least 18 different species in the Aspergillus section Flavi synthesize aflatoxins, which have been classified as group 1 carcinogens2,3. Of these, the most important species are A. flavus and A. parasiticus, which have been repeatedly isolated from cereals and nuts4–6. A. parasiticus is described as preferring ground crop hosts such as peanuts, whereas A. flavus infects a wider range of hosts7. It is noted that isolates of A. flavus L strain, which are characterized by the formation of a few large sclerotia (>400 nm in diameter) but numerous conidia, are better adapted to the phyllosphere compared to isolates of A. flavus S strain8,9. In contrast, the latter ones form many small sclerotia (<400 nm) but only a few conidia, indicating that they are better adapted to the soil habitat8,10. A. minisclerotigenes, another highly aflatoxigenic Aspergillus species, was initially isolated from peanuts and has been consistently found in maize samples11–13. However, the preferred habitat of this species has not yet been defined. Knowledge about the preferred habitat and chemotype of fungal species allows for the assessment of health risks associated with a contaminated crop. The presumed chemotype and also the preferred habitat can be deduced to some extent from the interpretation of genomic data. For aflatoxigenic Aspergillus species, for example, different aflatoxin chemotypes exist due to modifications within the aflatoxigenic genes in the form of deletions, insertions, or single nucleotide polymorphisms, despite the structure of the aflatoxin gene cluster, which includes these approximately 30 genes necessary for the biosynthesis of aflatoxins, is mostly conserved14,15. A. parasiticus and A. minisclerotigenes have been described as capable of producing all the four main aflatoxins B1 (AFB1), B2 (AFB2), G1 (AFG1) and G2 (AFG2). In contrast, most aflatoxigenic A. flavus strains only produce AFB1 and AFB23,11,16. This is due to different deletions in A. flavus genomes, including the 5’-end of the coding regions of the aflF and aflU genes, as well as their intergenic region, at the 5’-end of the aflatoxin gene cluster. These two genes are necessary for the production of G-group aflatoxins and code for a cytochrome P450 monooxygenase and an aryl alcohol dehydrogenase17,18. The evolutionary background behind the diversity of various aflatoxin chemotypes and the benefit for the fungus to produce either all or only certain aflatoxins remain unclear. Additionally, the preference for a specific habitat based on genomic data has not been determined until now. Aflatoxins are particularly problematic for food safety in regions with warm and humid climates, like sub-Saharan Africa. In Eastern Kenya, for example, high levels of contamination were detected in maize samples from households. In detail, 75% of the controlled samples contained aflatoxins at levels that exceeded the regulatory limit in Kenya (10 µg/kg)19. The consumption of contaminated food can result in chronic and acute toxic effects20. A significant outbreak of acute aflatoxicosis, resulting in 317 documented cases and 125 deaths, occurred in Eastern Kenya in 2004 due to the consumption of contaminated homegrown maize21. Highly aflatoxigenic Aspergillus sp. strains, which are genetically more closely related to A. minisclerotigenes than to A. flavus and are endemic to Eastern Kenya, were identified as the responsible strains for this outbreak of acute aflatoxicosis22. Later, it was stated that these “lethal aflatoxicosis fungus” strains fell within the genetic diversity of A. minisclerotigenes23.These strains are also characterized by the inability to form G-group aflatoxins due to a deletion of 2.2 kb within the aflatoxin gene cluster, which includes parts of the aflF and aflU genes22. To our knowledge, the whole genome sequence of such a strain has not yet been thoroughly investigated until now.

Some species in the Aspergillus section Flavi produce the mycotoxin cyclopiazonic acid (CPA) in addition to aflatoxins. CPA is an indole-tetramic acid that has been found to cause, among others, neurological disorders and gastrointestinal complaints in animal studies24. In A. flavus, there are three genes (maoA, dmaT, and pks-nrps) that encode the proteins for CPA biosynthesis. Additionally, there is one gene (ctfR1) that appears to regulate the expression of these structural genes. These genes are located in a small gene cluster near the 5’-end of the aflatoxin gene cluster25. A. parasiticus is described as being unable to produce CPA3,26 and the inability to amplify genes of the CPA gene cluster suggests a deletion in this species26. However, the exact extent of this modification has not yet been described. In addition to the aforementioned knowledge gaps, no comprehensive analysis has been conducted at the genomic and analytical levels for important Aspergillus section Flavi species found in Eastern Kenya. However, this analysis is essential in the event of an aflatoxicosis outbreak to quickly predict the severity and extent of the disease and make an appropriate risk assessment in the case of maize infestation with aflatoxin-producing fungi. Thus, in this study, the genomes of strains of A. flavus, A. minisclerotigenes, and A. parasiticus were compared to each other. For this, genomic data of the strains A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410, which were recently sequenced by us27–29, was used and analyzed in detail. Partial gene sequences of calmodulin, nitrate reductase, and ß-tubulin revealed a high homology between strain A. minisclerotigenes MRI390 and strains of the so-called lethal aflatoxicosis fungi (e.g., Aspergillus sp. A1168)28. Indications for the preferred habitat of these different species were investigated at the genomic level. Further, secondary metabolite biosynthetic gene clusters (BGCs) were analyzed with a specific focus on the aflatoxin gene cluster. These clusters were associated with the mycotoxin patterns (chemotypes) produced by the respective strains.

Results

Genome annotation

Characteristics of the genome assemblies of A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410 are published elsewhere27–29. The ab initio eukaryotic gene finding of the genome assemblies using Augustus (v. 3.4.0) predicted the highest number of protein coding genes in the genome of A. parasiticus (11,914). This was followed by A. minisclerotigenes MRI390 (11,551) and MRI400 (11,501), as well as A. flavus (11,489), which were mostly annotated (Table 1).Table 1 Characteristics of the genome assemblies and annotation of A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410

Parameters	A. flavus MRI19	A. minisclerotigenes MRI390	A. minisclerotigenes MRI400	A. parasiticus MRI410	
No. of predicted genes	11,489	11,551	11,501	11,914	
No. of annotated genes (Gene Ontology)	9841	9754	9871	10,109	
Total no. of BUSCO orthologs	1706	1706	1706	1638	
Complete single-copy, complete multicopy, fragmented, and missing orthologs (%)	93.7, 0.9, 0.5, 4.9	93.7, 0.9, 0.5, 4.9	93.4, 0.8, 0.6, 5.2	96.0, 0.9, 0.1, 3.0	

Whole-genome alignment

A whole-genome alignment was performed on the four sequenced strains and in addition, for comparison, with A. flavus AF1330, A. minisclerotigenes CBS 11763531, and A. parasiticus CBS 11761831. A high overall similarity with some larger modifications among the different species was revealed (Fig. 1). In the genome of A. flavus MRI19, two variations were detected in relation to the aforementioned strains. An approximately 360 kb fragment present on chromosome 8 of the other strains was aligned to chromosome 2 in A. flavus MRI19. Further, on chromosome 4, an approximately 210 kb fragment was invertedly incorporated. In A. minisclerotigenes MRI390, approximately 223 kb of chromosome 2 were observed to be inverted compared to the other strains. In addition, smaller modifications were detected. Based on the alignment, an average nucleotide identity of more than 93% was calculated across the whole genomes between all strains (Table 2). The average nucleotide identity between the A. flavus and A. minisclerotigenes isolates was approximately 97.0%, between A. flavus and A. parasiticus approximately 93.9%, and between A. minisclerotigenes and A. parasiticus approximately 94.0%.Fig. 1 Whole-genome sequence-based synteny analysis.

A whole-genome alignment was performed using the CLC Genomics Workbench for the genome assemblies of A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410. Additionally, the genome assemblies of A. flavus AF13 (GenBank ACCN: GCA_014117485.1)30, A. minisclerotigenes CBS 117635 (ACCN: GCA_009176455.1)31, and A. parasiticus CBS 117618 (ACCN: GCA_009176385.1)31 were included for comparison. Modifications mentioned in the text are marked in red.

Table 2 Average nucleotide identity [%] (above the diagonal) and alignment percentage [%] (below the diagonal) between A. flavus, A. minisclerotigenes, and A. parasiticus genomes

Fungal strains	A. flavus MRI19	A. flavus AF13	A. minisclerotigenes MRI390	A. minisclerotigenes MRI400	A. minisclerotigenes CBS 117635	A. parasiticus MRI410	A. parasiticus CBS 117618	
A. flavus MRI19		99.46	97.00	96.98	97.00	93.85	93.88	
A. flavus AF13	96.65		96.97	96.95	96.98	93.83	93.86	
A. minisclerotigenes MRI390	90.80	89.92		99.16	98.88	93.97	93.98	
A. minisclerotigenes MRI400	91.22	90.30	96.53		98.88	93.97	93.98	
A. minisclerotigenes CBS 117635	92.29	91.23	94.95	95.41		93.96	93.97	
A. parasiticus MRI410	87.60	86.73	86.86	87.29	87.99		98.57	
A. parasiticus CBS 117618	88.08	87.09	87.01	87.45	88.45	94.92		

Strain-specific genes

Strain-specific genes were identified among the predicted genes in the genomes of A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410. These genes were defined as having <50% similarity or coverage to the nucleotide sequences of genes from the other sequenced strains. Based on these criteria, more than 10,600 genes were identified as homologous (Fig. 2). The genome of A. parasiticus MRI410 had the highest number of strain-specific genes (533, Fig. 2). An enrichment analysis was conducted to compare the strain-specific genes of an isolate with all predicted genes of the same isolate. The analysis revealed that A. flavus MRI19 is most over-represented in the Gene Ontology (GO) terms “endopeptidase activity” (GO:0004175) and “serine-type endopeptidase activity” (GO:0004252). The genome of A. parasiticus MRI410 is most over-represented in the GO terms “carbon-oxygen lyase activity, acting on phosphates” (GO:00016838) and “pyruvate dehydrogenase activity” (GO:0004738). The most enriched GO terms for the strain-specific genes of A. minisclerotigenes were “anchored component of membrane” (GO:0031225) and “cell adhesion mediator activity” (GO:0098631) for strain MRI390, and “DNA integration” (GO:0015074) and “corticosterone binding” (GO:1903875) for strain MRI400 (Supplementary Data 1).Fig. 2 Venn diagram illustrating the numbers (mean values) of homologous and strain-specific genes identified in the genomes of A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410.

Genes were classified as strain-specific if they had <50% coverage and/or <50% similarity to the nucleotide sequences of all genes in the other strains.

Carbohydrate-active enzymes (CAZymes)

CAZymes, which are enzymes required for the synthesis, modification, and degradation of complex carbohydrates and glycoconjugates, serve as indicators of a fungal species’ ability to metabolize specific carbon sources32,33. CAZymes are also relevant for the degradation of host cell wall components, making them crucial in infection and colonization processes. CAZymes were identified in the genomes of A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410 and were classified into the different modules based on their function, including carbohydrate esterase, auxiliary activity (AA), glycoside hydrolase (GH), glycosyl transferase, polysaccharide lyase, and carbohydrate-binding module32. Overall, a high similarity in the identified families and subfamilies was observed among the sequenced strains (Fig. 3, Supplementary Data 2). The highest deviation was detected for module AA, which includes catalytic enzymes involved in the degradation of plant cell walls. A total of 78 enzymes were detected for A. parasiticus MRI410, 74 and 76 for A. minisclerotigenes MRI390 and MRI400, but only 68 for A. flavus MRI19. Here, the most striking differences were seen in the family AA4 and the subfamily AA3_2, which include lignin-degrading and oxidizing enzymes34. Among these enzymes, the lowest number was found in A. flavus MRI19 (AA4: 0, AA3_2: 24). However, more enzymes were detected in A. minisclerotigenes (MRI390: AA4: 1, AA3_2: 27; MRI400: AA4: 2, AA3_2: 29) and A. parasiticus MRI410 (AA4: 3, AA3_2: 28) (Supplementary Data 2). Further differences include the increased presence of α-rhamnosidases from the family GH78 in the GH module, but a decreased availability of chitinases from the family GH18 and β-1,3-glucanases from the family GH55 in A. flavus and A. parasiticus compared to both strains of A. minisclerotigenes (Supplementary Data 2).Fig. 3 Carbohydrate-active enzymes (CAZymes) were identified in the genomes of A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410 using the dbCAN2 meta server.

These enzymes were classified into different enzyme modules.

BGCs

BGCs were identified in the genomes of A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410 using the fungal version of antiSMASH and were categorized into the following groups based on the key enzyme of the gene cluster: non-ribosomal peptide synthetase (NRPS) clusters, NRPS-like clusters, siderophore clusters, type I polyketide synthase (T1PKS) clusters, type III polyketide synthase (T3PKS) clusters, indole clusters, β-lactone containing protease inhibitor clusters, fungal post-translationally modified peptide (RiPP) clusters, and terpene clusters. Additionally, hybrid clusters have been detected. Differences could be observed among the sequenced Aspergillus species (Fig. 4). The genome of A. parasiticus MRI410 had the highest number of BGCs (55), followed by A. minisclerotigenes MRI390 (54). In the genomes of A. minisclerotigenes strains MRI390 and MRI400, one siderophore cluster was detected on contigs aligned to sequences of the database of the National Center for Biotechnology Information assigned to chromosome 3. The nucleotide sequence was mostly deleted in the other genomes. Within the - compared to the other strains - invertedly incorporated fragment on chromosome 4 of A. flavus MRI19 (see section “Whole-genome alignment”), an unknown BGC with a NRPS as the backbone gene was identified. Further, only in the genome of A. parasiticus MRI410 a RiPP cluster was detected.Fig. 4 Secondary metabolite biosynthetic gene clusters were predicted in the genomes of A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410 using the fungal version of antiSMASH.

Gene clusters were categorized into different groups based on their backbone gene(s) (NRPS, non-ribosomal peptide synthetase; T1PKS, type I polyketide synthase; T3PKS, type III polyketide synthase; RiPP, post-translationally modified peptide).

Aflatoxin gene cluster

The aflatoxin gene cluster, which was identified by antiSMASH as a T1PKS/NRPS/indole hybrid cluster in A. flavus MRI19 and A. minisclerotigenes MRI390 and MRI400, and as a T1PKS cluster in A. parasiticus MRI410, was found in all the four analyzed genomes on a contig that was assigned to chromosome 3. Among the four strains, a nucleotide identity of > 94% for the complete clusters and >96% for only the coding sequences (A. flavus/A. minisclerotigenes > 96.6%, A. flavus/A. parasiticus > 96.6%, A. parasiticus/A. minisclerotigenes > 96.2%, A. minisclerotigenes MRI390/A. minisclerotigenes MRI400 > 98.4%) was calculated using the CLC Genomics Workbench. Further, the gene clusters were visualized using GenVision Pro (Lasergene v. 17, DNASTAR, Inc, Madison, USA) (Fig. 5). All known genes of the aflatoxin gene cluster were present in all the seven analyzed strains, including A. flavus AF1330, A. parasiticus CBS 11761831, and A. minisclerotigenes CBS 11763531 in addition to the aforementioned four strains. The comparison of the aflatoxin gene cluster among the different Aspergillus species revealed that the majority of variations were observed at the 5’-end within the genes aflF and aflU. A deletion of 2176 bp, which includes the 5’-end of both genes and the entire intergenic region, was detected in the genomes of A. minisclerotigenes MRI390 and MRI400 compared to A. minisclerotigenes CBS 117635. Further, in A. flavus MRI19, a fragment of 939 bp including also the 5’-end of both genes and the complete intergenic region was deleted.Fig. 5 Comparison of the aflatoxin gene cluster among Aspergillus species.

Genes of the aflatoxin gene cluster from A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410 were visualized using GenVision Pro. For comparison, the gene sequences of the following strains were used from the database of the NCBI: A. flavus AF13 (GenBank ACCN: GCA_014117485.1)30, A. minisclerotigenes CBS 117635 (ACCN: GCA_009176455.1)31, and A. parasiticus CBS 117618 (ACCN: GCA_009176385.1)31. Modifications mentioned in the text are marked in red.

A further deletion of 23 bp was observed in the transcribed region of aflR in A. minisclerotigenes MRI390 and MRI400, but not in A. minisclerotigenes CBS 117635 or the other fungal species. All additional detected extensive variations were located in intergenic regions, such as a 401 bp deletion in the intergenic region between aflY and nadA in all A. flavus and A. minisclerotigenes strains, in comparison to both A. parasiticus strains. As a further small but highly relevant modification, a substitution (G → T) was detected within the coding sequence of the gene aflL of A. parasiticus MRI410 at the 430th base of this gene. This led to the translation of the stop codon TAG, suggesting that this gene is non-functional. The gene aflL encodes the desaturase that converts versicolorin B to versicolorin A and thus introduces the double bond in the furan ring for the synthesis of group 1 aflatoxins15,35,36. Thus, this step is the branching point for the formation of group 1 and 2 aflatoxins.

CPA gene cluster

The complete CPA gene cluster or its remnants were identified in the genomes of A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410 located near the 5’-end of the aflatoxin gene cluster, and thus also on contigs assigned to chromosome 3. The alignment of this region revealed a high homology between A. flavus MRI19 and A. minisclerotigenes MRI390 and MRI400. These strains contained the complete gene cluster. In the genome of A. parasiticus MRI410, the genes dmaT and maoA were completely deleted as most of its pks/nrps gene (Fig. 6). A deletion of 18,450 bp starting approximately 3200 bp away from the 5’-end of the aflF gene of the aflatoxin gene cluster was detected, spanning over most of the CPA gene cluster. Between two remaining fragments of the pks/nrps gene, a second deletion of 1206 bp and an insertion of 7671 bp were detected in the genome of A. parasiticus MRI410. The inserted sequence contained two genes encoding a putative gluconolactone oxidase and a glucose oxidase precursor. This sequence was also found in the genomes of A. minisclerotigenes and A. flavus. In both A. minisclerotigenes strains, it was, however, located close to the end of the contig representing the 3’-end of chromosome 3. In A. flavus MRI19, it was located approximately 6620 bp away from the 5’-end of chromosome 3. In the genome of A. parasiticus MRI410, no additional copies of the inserted sequence were detected elsewhere.Fig. 6 Schematic presentation of the region between the cyclopiazonic acid (CPA) gene cluster (GC) and the aflatoxin gene cluster of A. parasiticus MRI410.

The lengths of the insertion and the deletions were determined by aligning the region that includes the remnants of the CPA gene cluster and the 5’-end of the aflatoxin gene cluster of A. parasiticus MRI410 with the corresponding sequence of A. minisclerotigenes MRI390. The genes located within this region are indicated by black bars.

Quantitation of aflatoxins

To determine if the different aflatoxin chemotypes correspond to the genetic-structural variations observed in the aflatoxin gene clusters, the fungal strains were incubated for 7 days on malt extract agar (MEA), Czapek yeast autolysate agar (CYA), and yeast extract sucrose agar (YES), and their biosynthesized aflatoxins were measured. Growth rates on laboratory media were similar among the different species, with all reaching a colony diameter of 7.3–7.7 cm on MEA, 7.9–8.8 cm on CYA, and 8.7–8.9 cm on YES. Aflatoxin levels in methanolic extracts were measured using LC-MS with high mass resolution. For the identification of the aflatoxins in the samples, the accurate masses, retention times, isotope ratios, and MS/MS spectra were compared to those of the corresponding standards. Formed aflatoxin levels are shown in Fig. 7 and in the Supplementary Data 3 and 4. The levels varied strongly between extracts of the same stain incubated on different media and between different strains incubated on the same medium. AFB1 was the most abundant aflatoxin produced by A. flavus MRI19 and A. minisclerotigenes MRI390 and MRI400, regardless of the medium on which the three strains were incubated. A. flavus produced the highest amounts of AFB1 on MEA, whereas the A. minisclerotigenes strains produced the highest amounts on YES (Supplementary Data 3 and 4). Except for incubation on MEA, A. minisclerotigenes MRI400 formed clearly higher levels of AFB1 than A. flavus MRI19 and A. minisclerotigenes MRI390. AFB2 and aflatoxin M1 (AFM1) were produced on the three media in low levels of approximately 1–2% compared to AFB1, and even lower levels of aflatoxin M2 (AFM2) but no AFG1 and AFG2 were detected in extracts of the three strains. This fits to their deletion that was observed within the aflF and aflU genes. Further, aflatoxicol (AFL), a hydroxylation product of AFB1, was analyzed. Among the investigated strains, the highest levels of AFL were measured in extracts of A. flavus MRI19, regardless of the medium on which they were incubated. Within the three individual strains (A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400), the highest levels of AFL were detected in extracts of fungi incubated on YES.Fig. 7 Aflatoxin levels formed by Aspergillus flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410 on MEA, CYA, and YES medium.

Aflatoxin levels were measured by LC-MS in extracts of A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410 after 7 days of incubation on MEA (A), CYA (B), and YES (C), respectively. The values are mean values respectively ± standard deviation. Limit of quantitation for all analytes: 0.01 µmol/L. AFB1, aflatoxin B1; AFB2, aflatoxin B2; AFG1, aflatoxin G1; AFG2, aflatoxin G2; AFM1, aflatoxin M1; AFM2, aflatoxin M2; AFL, aflatoxicol.

A. parasiticus MRI410 showed a significantly different chemotype compared to the other strains (Figs. 7 and 8). The aflatoxin produced in the highest quantity by this strain was AFG2, which was measured at the highest levels in extracts of the fungus incubated on YES medium. AFB2 was biosynthesized by A. parasiticus MRI410 in relative amounts of 3.87% on CYA, 16.38% on MEA, and 21.45% on YES in relation to AFG2 levels. Further, AFM2 was measured at relatively low levels in all extracts, regardless of the medium on which A. parasiticus MRI410 was incubated. Low amounts of AFB1 and AFG1 could be detected at concentrations below the limit of quantitation, but not AFM1 or AFL. Additionally, CPA was detected in extracts of A. flavus and A. minisclerotigenes, but not in those of A. parasiticus. However, as CPA was not the focus of this publication, the LC-MS method was not optimized and validated for this mycotoxin and CPA was thus only semi-quantified.Fig. 8 Schematic presentation of a part of the aflatoxin (AF) biosynthesis pathway of A. parasiticus MRI410.

Due to a substitution within the aflL gene, which encodes the desaturase responsible for converting versicolorin B to versicolorin A, versicolorin B is mostly converted into group 2 aflatoxins. The known structures are based on Caceres et al.15. DMST, demethylsterigmatocystin.

Identification of precursors of aflatoxins

The identification of aflatoxin precursors was not the main focus of this publication. However, since full mass scans were generated in high-resolution of the samples, it was also possible to search selected precursors of aflatoxins in extracts of A. parasiticus MRI410 without quantifying these compounds. The identification was based on the accurate masses and isotope patterns, as well as a comparison of the MS/MS spectra with respective data in the literature, if available37,38. No precursors of group 1 aflatoxins, which are typically formed in the biosynthetic pathway after the desaturation of versicolorin B to versicolorin A by AflL, such as sterigmatocystin, O-methylsterigmatocystin, and 11-hydroxy-O-methylsterigmatocystin, were detected. In contrast, versicolorin B and other precursors of group 2 aflatoxins, such as dihydro-O-methylsterigmatocystin and dihydro-11-hydroxy-O-methylsterigmatocystin, were detected (Supplementary Data 5). Additionally, a metabolite with a mass of m/z 329.0665 and the suggested formula C17H12O7 was detected by LC-MS in methanolic extracts of A. parasiticus MRI410, but not in those of the other investigated strains. Based on the evaluation of the respective MS/MS spectrum and the accurate mass, it could be suggested that this analyte may be a metabolite of versicolorin B. The compound differs from versicolorin B (C18H12O7) in accurate mass by 12.0003 Da and in molecular formula most likely by one fewer carbon atom.

Discussion

The discrepancy between the high structural conservation of the aflatoxin gene cluster observed in the species Aspergillus on the one hand and the diversity reflected in the different chemotypes on the other hand, i.e., the evolutionary background of which chemotype generates which benefit for the producing fungus in its specific habitat, can be determined by analyzing the genome. The aim of the study was therefore to comprehensively compare the genomes of several toxicologically relevant strains of Aspergillus section Flavi at both the genetic and analytical levels. These strains have different profiles of aflatoxin biosynthesis. In the genome of the analyzed A. parasiticus strain, a substitution was detected within the aflL gene of the aflatoxin gene cluster. This gene codes for a desaturase, which disables this strain from producing significant levels of group 1 aflatoxins (AFB1, AFG1, AFM1). Further, both analyzed A. minisclerotigenes strains contain a deletion that includes the 5’-end of the coding sequences of aflF and aflU in the aflatoxin gene cluster. This deletion prevents the biosynthesis of G-group aflatoxins. In addition, this study has provided clear references of genetic indicators for the adaption of species to certain habitats. However, it has to be considered that only one respectively two representative strains of each species were included in the study.

A. parasiticus has been described as preferring ground crop hosts and is mostly isolated from soil or peanut samples7. In contrast, L strain isolates of A. flavus appear to prefer the phyllosphere10,39. The preferred habitat of A. minisclerotigenes has not been described yet. Appropriately, an increased number of genes coding for CAZymes classified to the module AA were identified in the genome of A. parasiticus compared to A. flavus. These enzymes are necessary for various functions, including the degradation of lignin and cellulose, the metabolism of metals, and the maintenance of iron homeostasis34. They can thus be seen as an indicator of the rather saprophytic lifestyle of this species in the soil habitat. A further indicator for the soil habitat is the fungal RiPP cluster, which was detected in the genome of A. parasiticus but not in the genomes of A. flavus and A. minisclerotigenes. Fungal RiPP clusters are associated, for example, with competitive defense and mycophagic activities40. Furthermore, various modifications, including two deletions and one insertion, were observed within the CPA gene cluster of A. parasiticus. Consequently, we did not detect any CPA in methanolic extracts of this species. As it can be assumed that multiple events caused these genetic variations, it can be suggested that there was no evolutionary pressure to preserve this cluster. Chalivendra et al.41 reported that A. flavus mutants, which were unable to produce CPA, exhibited lower virulence in causing ear rot disease. On the other hand, non-mutated strains that colonized maize were found to be strong CPA producers. The authors suggest that CPA is produced by fungi as a pathogenicity factor because of its cytotoxicity, which impairs the host’s defense response. Therefore, it can be inferred that the biosynthesis of CPA by A. parasiticus was not essential for this fungus, which is more saprophytic than parasitic. As a result, the maintenance of the CPA gene cluster was obviously discontinued. Alternatively, the inability to produce CPA may have shifted the fungus to the soil habitat. A. minisclerotigenes has already been isolated from soil as well as maize samples12,13,42, indicating its adaptation to both habitats. As for A. parasiticus, an increased number of enzymes necessary for lignin degradation compared to A. flavus, and additionally, magnesium and cobalt transport proteins encoded by species-specific genes of A. minisclerotigenes, might be indicators of adaptation for survival in the soil. In addition, a high CPA biosynthesis was observed for A. minisclerotigenes in this study. As mentioned previously, this could be advantageous for the infection of maize plants. Further, in comparison to the other species, A. minisclerotigenes seemed to be well-equipped for iron uptake. This was evident from the presence of a gene cluster coding for a siderophore, an iron-chelating molecule essential for the fungus’s high-affinity iron uptake43, in the genomes of both sequenced A. minisclerotigenes strains, but not in the genomes of the other species. A high bioavailability of iron, which is essential for the formation of co-factors for numerous enzymes and which is poorly available in most environments, could be a competitive advantage to this species43. A. minisclerotigenes appears to be particularly well adapted to survive in certain regions of Kenya12,13, where strains of this species were most likely responsible for outbreaks of acute aflatoxicosis22. A. minisclerotigenes MRI390 was shown to be genetically very closely related to strains that caused such outbreaks28. Due to the significant impact of the species A. minisclerotigenes on food safety, further research on it is recommended. The third analyzed species, A. flavus, is the most extensively studied species7,44, so the discussion will not focus on it.

In agreement with earlier investigations31, the comparison of the aflatoxin gene cluster generally showed a high degree of homology among the different species. However, relevant structural modifications were identified in all examined Aspergillus strains, and different corresponding chemotypes were observed. While the deletions within aflF and aflU in the genomes of A. flavus and A. minisclerotigenes have already been described17,22,45, the substitution detected within the aflL gene of A. parasiticus MRI410 has not been previously described, to the best of our knowledge. High levels of group 2 aflatoxins (AFB2, AFG2, AFM2) were measured in extracts of this strain, while there were hardly any group 1 aflatoxins (AFB1, AFG1, AFM1). This is highly atypical. Usually, the levels of group 1 aflatoxins are 10–100 times higher than those of group 2 aflatoxins46–48. The aflatoxin precursors detected in methanolic extracts of this strain, i.e., the detection of late precursors of group 2 but not of group 1 aflatoxins, aligns with this aflatoxin profile. The detected levels of AFB1 and AFG1 in methanolic extracts of A. parasiticus, which were found to be below the limit of quantitation, could potentially have been formed by other desaturases or through non-enzymatic processes in order to compensate for the one non-functional desaturation. It is highly likely that all subsequent steps in the pathway for the biosynthesis of group 1 aflatoxins are functional. Possible purposes for the production of aflatoxins for the fungus in general are discussed in the literature, such as their effect against insects or their role in protecting cells against the accumulation of reactive oxygen species49,50. However, the exact roles of the different aflatoxins, such as group 1 and 2 aflatoxins, as well as group B, G, and M aflatoxins, are not yet fully understood. Further research to understand why certain Aspergillus species produce G-group aflatoxins while others are also competitive with the deletion within aflF and aflU, or why A. parasiticus MRI410 is apparently able to survive without producing significant amounts of group 1 aflatoxins, while other Aspergillus strains primarily produce these compounds, would be interesting. The comparative study presented here focusing on three different toxicologically important Aspergillus species contributes to the ongoing research on aflatoxin and shows once again the complexity of the regulation and physiological background of mycotoxin biosynthesis by filamentous fungi.

Materials and methods

Fungal strains, genome sequencing and assembly, and functional annotation

The fungal strains A. flavus MRI19, originally isolated from tiger nuts grown in Spain27, A. minisclerotigenes strains MRI390 and MRI400, isolated from maize samples from Makueni, Kenya28, and A. parasiticus MRI410, isolated from a soil sample of a maize field in Katumani, Kenya29, were used for this study. The genomes of the strains were previously sequenced by us using short-read sequencing on a MiSeq instrument (Illumina, San Diego, USA) and additionally long-read sequencing on a PacBio Sequel instrument (Pacific Biosciences, Menlo Park, USA) for A. flavus MRI19 and A. minisclerotigenes MRI390 and MRI400. For A. parasiticus MRI410, long-read sequencing was performed on a Nanopore MinION Mk1C instrument (Oxford Nanopore Technologies, Oxford, UK) in addition to short-read sequencing. The software SPAdes (v. 3.14.1) was used with default software parameters to generate de novo hybrid assemblies of A. flavus MRI19 and A. minisclerotigenes MRI400, combining data from short and long-read sequencing, resulting in 68 and 52 contigs, respectively27,28,51,52. For the strains A. minisclerotigenes MRI390 and A. parasiticus MRI410, de novo assemblies of the long-reads were generated using the software Flye (v. 2.8.2) with default software parameters, resulting in 28 and 60 contigs, respectively28,29,52,53. Subsequently, the respective MiSeq data was aligned to the contigs using BWA (v. 0.7.17), SAMtools (v. 1.10), and QualiMap (v. 2.2.1), and the alignment was polished using Pilon (v. 1.23) and SAMtools (v. 1.10)52,54–57. Repeat masking was performed on the genome assemblies using RepeatMasker (v. 4.0.9) and Dfam (v. 3.0)58,59. Additionally, ab initio eukaryotic gene finding was conducted on both strands, allowing for the identification of partial genes using Augustus (v. 3.4.0), choosing A. oryzae as the closest species60. Identified genes were searched in the non-redundant protein sequence database performing a BLAST analysis (blastx-fast, e-value threshold 1.0*E-3). The respective taxonomic filter for each species was selected before mapping and annotating blast hits61. For the further functional annotation, InterProScan (v. 5.52-86.0) and EggNOG mapper (v. 2.1.0 with EggNOG v. 5.0.2) were used62,63. All steps were performed using OmicsBox (v. 17.0.2, BioBam Bioinformatics Solutions, Valencia, Spain), and default software parameters were used, unless otherwise specified. BUSCO (Benchmarking Universal Single-Copy Orthologs) analyses were performed using the lineage database ascomycota_odb10 to evaluate the completeness of the genome assembly data sets64.

Whole-genome alignment

A whole-genome alignment was performed between the genome assemblies of A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410 using the Whole Genome Alignment Plugin of the CLC Genomics Workbench (v. 22.0.1, Qiagen, Hilden, Germany) without aligning based on individual chromosomes. The alignment had a minimum initial seed length of 15 nucleotides and a minimum alignment block length of 50 nucleotides. Based on this, the software calculated the average nucleotide identity using a minimum similarity fraction and a minimum length fraction of 0.8. The software further created a tree using the Neighbor-Joining algorithm on the distance matrix calculated from the seed matches between the individual genomes.

Prediction of strain-specific genes and CAZymes

A BLAST database was generated of the nucleotide sequences of the predicted genes from the strains A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410 separately using OmicsBox (v. 17.0.2, BioBam Bioinformatics Solutions, Valencia, Spain). Then, the nucleotide sequences of the genes from each isolate were compared to those of the other three strains, respectively (tblastn, blast expectation value: 1.0*E-3). To predict specific genes of a strain, the results were filtered to include only genes that had either no or <50% similarity or <50% coverage to the nucleotide sequences of the other three strains. An enrichment analysis (Fisher’s exact test) was performed on the resulting genes against all genes of the respective strain (p value: 0.05)65. CAZymes were predicted in the gene sequences of the analyzed genomes using HMMER (e-value < 10-15, coverage > 0.35) against the dbCAN database and DIAMOND (e-value < 10-102) against the CAZy database using the dbCAN meta server (http://bcb.unl.edu/dbCAN2/blast.php)32. Genes identified by both tools were retained and categorized into the following modules based on their functionality: auxiliary activity (AA), carbohydrate-binding module, carbohydrate esterase, glycoside hydrolase (GH), glycosyl transferase, and polysaccharide lyase.

Prediction of BGCs and comparison of the aflatoxin gene cluster

BGCs were predicted using antiSMASH (v. 6.1.0), employing the cluster finder algorithm for BGC border prediction and the relaxed detection strictness66,67. The aflatoxin gene cluster of the strains A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410, as well as that of A. flavus AF1330, A. parasiticus CBS 11761831, and A. minisclerotigenes CBS 11763531 was aligned using the CLC Genomics Workbench (v. 22.0.1, Qiagen, Hilden, Germany) to determine the average nucleotide identity, as described above.

Quantitation of mycotoxins

For the quantitation of mycotoxins produced by Aspergillus sp. strains the following reference standards were used: aflatoxin B1 (AFB1, >99%), aflatoxin B2 (AFB2, >99%), aflatoxin G1 (AFG1, >99%), aflatoxin G2 (AFG2, >99%), and aflatoxin M1 (AFM1, >98%) all dissolved in acetonitrile (Sigma-Aldrich Chemie GmbH, Taufkirchen, Germany), aflatoxin M2 (AFM2) in acetonitrile (>98%), and aflatoxicol (AFL, >99%) (Cfm Oskar Tropitzsch GmbH, Marktredwitz, Germany), and cyclopiazonic acid (CPA) in acetonitrile (>99.9%) (Romer Labs Division Holding GmbH, Tulln, Austria). Spore suspensions containing 1.0 × 104 spores per mL in Tween-80/NaCl-mixture (9 g/L NaCl [Carl Roth, Karlsruhe, Germany], 1 g/L Tween-80 [Serva, Heidelberg, Germany], 1 g/L agar [Agar-Agar Kobe Ι; Carl Roth, Karlsruhe, Germany]) were prepared of A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410. For spore counting, a Thoma cell counting chamber (Paul Marienfeld GmbH & Co. KG, Lauda-Königshofen, Germany) was used. MEA (15 g/L agar [Agar-Agar Kobe Ι], 30 g/L malt extract, 5 g/L tryptone/peptone ex casein [Carl Roth, Karlsruhe, Germany]), CYA (20 g/L agar, 30 g/L saccharose, 5 g/L yeast extract, 3 g/L sodium nitrate, 0.5 g/L magnesium sulfate heptahydrate, 0.01 g/L iron (II) sulfate heptahydrate, 1 g/L dipotassium hydrogen phosphate [Carl Roth, Karlsruhe, Germany], 0.5 g/L potassium chloride [Merck KGaA, Darmstadt, Germany])68, and YES (20 g/L agar, 150 g/L saccharose, 20 g/L yeast extract [Carl Roth, Karlsruhe, Germany])69 plates (Ø = 9 cm, 40 mL medium per Petri dish) were prepared. The plates were inoculated with 50 µL of each spore suspension as a one-point culture. Then, the fungal strains were cultivated for 7 days at 25 °C in the dark. At the end of the incubation period, two agar plugs (Ø = 8 mm) were taken from a colony, one from the center and one from the edge, using a sterile corer (n = 5 per strain and medium). The protocol for the following toxin extraction was based on the DIN EN ISO norm 16050:2011, with modifications as described70. One of the two agar plugs was transferred to a 2 mL tube filled with ceramic beads (Type A, Macherey-Nagel, Düren, Germany) and 1.25 mL of methanol:water (70:30, [v:v]) and 50 mg NaCl were added. The samples were homogenized (4 × 30 s, 6.0 m/s) using the FastPrep-24TM (MP Biomedicals, Eschwege, Germany). Then, the extract was transferred into a new 2 mL tube. Subsequently, the second agar plug of the respective sample was transferred into the emptied 2 mL tube filled with ceramic beads and further 1.25 mL of methanol:water (70:30, [v:v]) and 50 mg NaCl were added. The homogenization procedure was repeated (4 × 30 s, 6.0 m/s) and the extract was transferred once again into a new 2 mL tube. To wash out the emptied 2 mL tube with the ceramic beads, further 1.25 mL of methanol:water (70:30, [v:v]) was pipetted into it and the tube was then processed for 30 s at 6.0 m/s in the FastPrep-24TM. The washing solution was transferred to a third 2 mL tube. The two extracts and the washing solution were centrifuged for 5 min at 16,200 × g. The supernatants from the three corresponding tubes were then combined und filtered using a 0.2 µm PTFE filter (Puradisc-13, WhatmanTM, Merck, Darmstadt, Germany). The concentrations of the target analytes AFB1, AFB2, AFG1, AFG2, AFM1, AFM2, and AFL in the extracted samples were measured on a 1290 Infinity LC system (Agilent Technologies, Waldbronn, Germany) coupled with a Triple TOF 5600 mass spectrometer (AB Sciex, Darmstadt, Germany) using the method described in Schamann et al.71. Briefly, the separation was achieved on a Waters Cortecs UPLC C18 column (2.1 mm × 150 mm, 1.6 µm; Waters, Eschborn, Germany) equipped with a pre-column (Security Guard Ultra UHPLC C18; Phenomenex, Aschaffenburg, Germany). Gradient elution was conducted using aqueous ammonium acetate buffer (10 mmol/L) as eluent A and methanol as eluent B. The ion source of the MS was operated in positive ESI mode. For details, refer to Schamann et al.71. Each sample was measured in two different individual dilutions to ensure that the concentration of each analyte was within the calibrated concentration range (10–640 nmol/L). For data analysis, the software Sciex OS 2.0.1. (AB Sciex, Darmstadt, Germany) was used, and apart from this, it was performed as described in Schamann et al.71. Briefly, analytes were quantified based on standard curves applying linear regression with a weighting of 1/x2. The concentration of the lowest standard (0.01 µmol/L) was set as limit of quantitation of all the target analytes.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Supplementary information

Description of Additional Supplementary Materials

Supplementary Data 1

Supplementary Data 2

Supplementary Data 3

Supplementary Data 4

Supplementary Data 5

Reporting Summary

Supplementary information

The online version contains supplementary material available at 10.1038/s42003-024-06738-w.

Acknowledgements

This research was funded by the Federal Office for Agriculture and Food (BLE), reference AflaZ 2816PROC11. We thank Volker Müller for his excellent technical support as well as Adetoye Wolfgang Adeyemo and Dr. Andreas Dötsch for bioinformatic support. The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of Nature Communications and/or the editor(s). Nature Communications and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Author contributions

“Conceptualization, M.S-H.; methodology, A.S.; software, A.S.; validation, A.S.; formal analysis, A.S., and S.T.S.; investigation, A.S.; resources, M.S-H, R.G. and SK; data curation, A.S., S.T.S., and M.S-H; writing original draft preparation, A.S.; writing review and editing, A.S., M.S-H., S.T.S., R.G., and S.K.; visualization, A.S.; supervision, M.S-H. and S.T.S.; project administration, M.S-H.; funding acquisition, M.S-H. All authors have read and agreed to the published version of the manuscript.”

Peer review

Peer review information

Communications Biology thanks the anonymous reviewers for their contribution to the peer review of this work. Primary Handling Editor: Tobias Goris.

Funding

Open Access funding enabled and organized by Projekt DEAL.

Data availability

The genome assemblies of A. flavus MRI19, A. minisclerotigenes MRI390 and MRI400, and A. parasiticus MRI410 used in this study for comparative genomic analyses have been previously published27–29. Source data underlying Fig. 7 can be found in supplementary data sets 3/4 - currently states availability upon request, as well as the numerical values for the CAZymes screening. Parts of the data presented in this study are available only upon request from the corresponding author due to regulations of the Nagoya Protocol.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Allwood JG Wakeling LT Bean DC Fermentation and the microbial community of Japanese koji and miso: a review J. Food Sci. 2021 86 2194 2207 10.1111/1750-3841.15773 34056716
Allwood, J. G., Wakeling, L. T. & Bean, D. C. Fermentation and the microbial community of Japanese koji and miso: a review. J. Food Sci. 86, 2194–2207 (2021).34056716 10.1111/1750-3841.15773
2. International Agency for Research on Cancer. Aflatoxins. In IARC monographs on the evaluation of carcinogenic risks to humans. Chemical agents and related occupations: a review of human carcinogens. – Vol. 100F. (International Agency for Research on Cancer, Lyon, France, 2012), pp. 225–248.
3. Frisvad JC Taxonomy of Aspergillus section Flavi and their production of aflatoxins, ochratoxins and other mycotoxins Stud. Mycol. 2019 93 1 63 10.1016/j.simyco.2018.06.001 30108412
Frisvad, J. C. et al. Taxonomy of Aspergillus section Flavi and their production of aflatoxins, ochratoxins and other mycotoxins. Stud. Mycol. 93, 1–63 (2019).30108412 10.1016/j.simyco.2018.06.001
4. Wagacha JM Mutegi C Karanja L Kimani J Christie ME Fungal species isolated from peanuts in major Kenyan markets: emphasis on Aspergillus section Flavi Crop. Prot. 2013 52 1 9 10.1016/j.cropro.2013.05.004
Wagacha, J. M., Mutegi, C., Karanja, L., Kimani, J. & Christie, M. E. Fungal species isolated from peanuts in major Kenyan markets: emphasis on Aspergillus section Flavi. Crop. Prot. 52, 1–9 (2013).10.1016/j.cropro.2013.05.004
5. Taniwaki MH Pitt JI Magan N Aspergillus species and mycotoxins: occurrence and importance in major food commodities Curr. Opin. Food Sci. 2018 23 38 43 10.1016/j.cofs.2018.05.008
Taniwaki, M. H., Pitt, J. I. & Magan, N. Aspergillus species and mycotoxins: occurrence and importance in major food commodities. Curr. Opin. Food Sci. 23, 38–43 (2018).10.1016/j.cofs.2018.05.008
6. Del Palacio A Pan D Occurrence and toxigenic potential of Aspergillus section Flavi on wheat and sorghum silages in Uruguay Mycology 2020 11 147 157 10.1080/21501203.2020.1752321 32923022
Del Palacio, A. & Pan, D. Occurrence and toxigenic potential of Aspergillus section Flavi on wheat and sorghum silages in Uruguay. Mycology 11, 147–157 (2020).32923022 10.1080/21501203.2020.1752321
7. Amaike S Keller NP Aspergillus flavus Annu. Rev. Phytopathol. 2011 49 107 133 10.1146/annurev-phyto-072910-095221 21513456
Amaike, S. & Keller, N. P. Aspergillus flavus. Annu. Rev. Phytopathol. 49, 107–133 (2011).21513456 10.1146/annurev-phyto-072910-095221
8. Cotty PJ Virulence and cultural characteristics of two Aspergillus flavus strains pathogenic on cotton Phytopathology 1989 79 808 814 10.1094/Phyto-79-808
Cotty, P. J. Virulence and cultural characteristics of two Aspergillus flavus strains pathogenic on cotton. Phytopathology 79, 808–814 (1989).10.1094/Phyto-79-808
9. Ohkura M Cotty PJ Orbach MJ Comparative genomics of Aspergillus flavus S and L morphotypes yield insights into niche adaptation G3: Genes Genomes Genet. 2018 8 3915 3930 10.1534/g3.118.200553
Ohkura, M., Cotty, P. J. & Orbach, M. J. Comparative genomics of Aspergillus flavus S and L morphotypes yield insights into niche adaptation. G3: Genes Genomes Genet. 8, 3915–3930 (2018).10.1534/g3.118.200553
10. Mehl HL Aspergillus flavus diversity on crops and in the environment can be exploited to reduce aflatoxin exposure and improve health Ann. N. Y. Acad. Sci. 2012 1273 7 17 10.1111/j.1749-6632.2012.06800.x 23230832
Mehl, H. L. et al. Aspergillus flavus diversity on crops and in the environment can be exploited to reduce aflatoxin exposure and improve health. Ann. N. Y. Acad. Sci. 1273, 7–17 (2012).23230832 10.1111/j.1749-6632.2012.06800.x
11. Pildain MB Two novel aflatoxin-producing Aspergillus species from Argentinean peanuts Int. J. Syst. Evol. Microbiol. 2008 58 725 735 10.1099/ijs.0.65123-0 18319485
Pildain, M. B. et al. Two novel aflatoxin-producing Aspergillus species from Argentinean peanuts. Int. J. Syst. Evol. Microbiol. 58, 725–735 (2008).18319485 10.1099/ijs.0.65123-0
12. Okoth S Genetic and toxigenic variability within Aspergillus flavus population isolated from maize in two diverse environments in Kenya Front. Microbiol. 2018 9 57 10.3389/fmicb.2018.00057 29434580
Okoth, S. et al. Genetic and toxigenic variability within Aspergillus flavus population isolated from maize in two diverse environments in Kenya. Front. Microbiol. 9, 57 (2018).29434580 10.3389/fmicb.2018.00057
13. Oloo RD Genetic profiling of Aspergillus isolates with varying aflatoxin production potential from different maize-growing regions of Kenya Toxins 2019 11 467 10.3390/toxins11080467 31404960
Oloo, R. D. et al. Genetic profiling of Aspergillus isolates with varying aflatoxin production potential from different maize-growing regions of Kenya. Toxins 11, 467 (2019).31404960 10.3390/toxins11080467
14. Chang P-K Horn BW Dorner JW Sequence breakpoints in the aflatoxin biosynthesis gene cluster and flanking regions in non-aflatoxigenic Aspergillus flavus isolates Fungal Genet. Biol. 2005 42 914 923 10.1016/j.fgb.2005.07.004 16154781
Chang, P.-K., Horn, B. W. & Dorner, J. W. Sequence breakpoints in the aflatoxin biosynthesis gene cluster and flanking regions in non-aflatoxigenic Aspergillus flavus isolates. Fungal Genet. Biol. 42, 914–923 (2005).16154781 10.1016/j.fgb.2005.07.004
15. Caceres I Aflatoxin biosynthesis and genetic regulation: a review Toxins 2020 12 150 10.3390/toxins12030150 32121226
Caceres, I. et al. Aflatoxin biosynthesis and genetic regulation: a review. Toxins 12, 150 (2020).32121226 10.3390/toxins12030150
16. Varga J Frisvad J Samson R A reappraisal of fungi producing aflatoxins World Mycotoxin J. 2009 2 263 277 10.3920/WMJ2008.1094
Varga, J., Frisvad, J. & Samson, R. A reappraisal of fungi producing aflatoxins. World Mycotoxin J. 2, 263–277 (2009).10.3920/WMJ2008.1094
17. Ehrlich KC Chang P-K Yu J Cotty PJ Aflatoxin biosynthesis cluster gene cypA is required for G aflatoxin formation Appl. Environ. Microbiol. 2004 70 6518 6524 10.1128/AEM.70.11.6518-6524.2004 15528514
Ehrlich, K. C., Chang, P.-K., Yu, J. & Cotty, P. J. Aflatoxin biosynthesis cluster gene cypA is required for G aflatoxin formation. Appl. Environ. Microbiol. 70, 6518–6524 (2004).15528514 10.1128/AEM.70.11.6518-6524.2004
18. Ehrlich KC Scharfenstein LL Montalbano BG Chang P-K Are the genes nadA and norB involved in formation of aflatoxin G1? Int. J. Mol. Sci. 2008 9 1717 1729 10.3390/ijms9091717 19325828
Ehrlich, K. C., Scharfenstein, L. L., Montalbano, B. G. & Chang, P.-K. Are the genes nadA and norB involved in formation of aflatoxin G1? Int. J. Mol. Sci. 9, 1717–1729 (2008).19325828 10.3390/ijms9091717
19. Kagot V Comprehensive analysis of multiple mycotoxins and Aspergillus flavus metabolites in maize from Kenyan households Int. J. Food Microbiol. 2022 363 109502 10.1016/j.ijfoodmicro.2021.109502 34952410
Kagot, V. et al. Comprehensive analysis of multiple mycotoxins and Aspergillus flavus metabolites in maize from Kenyan households. Int. J. Food Microbiol. 363, 109502 (2022).34952410 10.1016/j.ijfoodmicro.2021.109502
20. Benkerroum N Chronic and acute toxicities of aflatoxins: mechanisms of action Int. J. Environ. Res. Public Health 2020 17 423 10.3390/ijerph17020423 31936320
Benkerroum, N. Chronic and acute toxicities of aflatoxins: mechanisms of action. Int. J. Environ. Res. Public Health 17, 423 (2020).31936320 10.3390/ijerph17020423
21. Azziz-Baumgartner E Case–control study of an acute aflatoxicosis outbreak, Kenya, 2004 Environ. Health Perspect. 2005 113 1779 1783 10.1289/ehp.8384 16330363
Azziz-Baumgartner, E. et al. Case–control study of an acute aflatoxicosis outbreak, Kenya, 2004. Environ. Health Perspect. 113, 1779–1783 (2005).16330363 10.1289/ehp.8384
22. Probst C Callicott KA Cotty PJ Deadly strains of Kenyan Aspergillus are distinct from other aflatoxin producers Eur. J. Plant. Pathol. 2012 132 419 429 10.1007/s10658-011-9887-y
Probst, C., Callicott, K. A. & Cotty, P. J. Deadly strains of Kenyan Aspergillus are distinct from other aflatoxin producers. Eur. J. Plant. Pathol. 132, 419–429 (2012).10.1007/s10658-011-9887-y
23. Houbraken J Classification of Aspergillus, Penicillium, Talaromyces and related genera (Eurotiales): an overview of families, genera, subgenera, sections, series and species Stud. Mycol. 2020 95 5 169 10.1016/j.simyco.2020.05.002 32855739
Houbraken, J. et al. Classification of Aspergillus, Penicillium, Talaromyces and related genera (Eurotiales): an overview of families, genera, subgenera, sections, series and species. Stud. Mycol. 95, 5–169 (2020).32855739 10.1016/j.simyco.2020.05.002
24. Ostry V Toman J Grosse Y Malir F Cyclopiazonic acid: 50th anniversary of its discovery World Mycotoxin J. 2018 11 135 148 10.3920/WMJ2017.2243
Ostry, V., Toman, J., Grosse, Y. & Malir, F. Cyclopiazonic acid: 50th anniversary of its discovery. World Mycotoxin J. 11, 135–148 (2018).10.3920/WMJ2017.2243
25. Chang P-K Ehrlich KC Cyclopiazonic acid biosynthesis by Aspergillus flavus Toxin Rev. 2011 30 79 89 10.3109/15569543.2011.576795
Chang, P.-K. & Ehrlich, K. C. Cyclopiazonic acid biosynthesis by Aspergillus flavus. Toxin Rev. 30, 79–89 (2011).10.3109/15569543.2011.576795
26. Chang P-K Horn BW Dorner JW Clustered genes involved in cyclopiazonic acid production are next to the aflatoxin biosynthesis gene cluster in Aspergillus flavus Fungal Genet. Biol. 2009 46 176 182 10.1016/j.fgb.2008.11.002 19038354
Chang, P.-K., Horn, B. W. & Dorner, J. W. Clustered genes involved in cyclopiazonic acid production are next to the aflatoxin biosynthesis gene cluster in Aspergillus flavus. Fungal Genet. Biol. 46, 176–182 (2009).19038354 10.1016/j.fgb.2008.11.002
27. Schamann A Geisen R Schmidt-Heydt M Draft genome sequence of an aflatoxin-producing Aspergillus flavus strain isolated from food Microbiol. Resour. Announc. 2022 11 e00894 21 10.1128/mra.00894-21 35112902
Schamann, A., Geisen, R. & Schmidt-Heydt, M. Draft genome sequence of an aflatoxin-producing Aspergillus flavus strain isolated from food. Microbiol. Resour. Announc. 11, e00894–21 (2022).35112902 10.1128/mra.00894-21
28. Schamann A Geisen R Schmidt-Heydt M Whole-genome sequences of two Kenyan Aspergillus minisclerotigenes strains Microbiol. Resour. Announc. 2023 12 e00219 e00223 37404178
Schamann, A., Geisen, R. & Schmidt-Heydt, M. Whole-genome sequences of two Kenyan Aspergillus minisclerotigenes strains. Microbiol. Resour. Announc. 12, e00219–e00223 (2023).37404178
29. Schamann A Geisen R Schmidt-Heydt M Whole-genome sequence of an Aspergillus parasiticus strain isolated from Kenyan soil Microbiol. Resour. Announc. 2023 12 e00203 e00223 37462395
Schamann, A., Geisen, R. & Schmidt-Heydt, M. Whole-genome sequence of an Aspergillus parasiticus strain isolated from Kenyan soil. Microbiol. Resour. Announc. 12, e00203–e00223 (2023).37462395
30. Fountain JC Two new Aspergillus flavus reference genomes reveal a large insertion potentially contributing to isolate stress tolerance and aflatoxin production G3: Genes Genomes Genet 2020 10 3515 3531 10.1534/g3.120.401405
Fountain, J. C. et al. Two new Aspergillus flavus reference genomes reveal a large insertion potentially contributing to isolate stress tolerance and aflatoxin production. G3: Genes Genomes Genet 10, 3515–3531 (2020).10.1534/g3.120.401405
31. Kjærbølling I A comparative genomics study of 23 Aspergillus species from section Flavi Nat. Commun. 2020 11 1106 10.1038/s41467-019-14051-y 32107379
Kjærbølling, I. et al. A comparative genomics study of 23 Aspergillus species from section Flavi. Nat. Commun. 11, 1106 (2020).32107379 10.1038/s41467-019-14051-y
32. Zhang H dbCAN2: a meta server for automated carbohydrate-active enzyme annotation Nucleic Acids Res. 2018 46 W95 W101 10.1093/nar/gky418 29771380
Zhang, H. et al. dbCAN2: a meta server for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 46, W95–W101 (2018).29771380 10.1093/nar/gky418
33. Cantarel BL The Carbohydrate-Active EnZymes database (CAZy): an expert resource for glycogenomics Nucleic Acids Res. 2009 37 D233 D238 10.1093/nar/gkn663 18838391
Cantarel, B. L. et al. The Carbohydrate-Active EnZymes database (CAZy): an expert resource for glycogenomics. Nucleic Acids Res. 37, D233–D238 (2009).18838391 10.1093/nar/gkn663
34. Levasseur A Drula E Lombard V Coutinho PM Henrissat B Expansion of the enzymatic repertoire of the CAZy database to integrate auxiliary redox enzymes Biotechnol. Biofuels 2013 6 41 10.1186/1754-6834-6-41 23514094
Levasseur, A., Drula, E., Lombard, V., Coutinho, P. M. & Henrissat, B. Expansion of the enzymatic repertoire of the CAZy database to integrate auxiliary redox enzymes. Biotechnol. Biofuels 6, 41 (2013).23514094 10.1186/1754-6834-6-41
35. Kelkar HS Skloss TW Haw JF Keller NP Adams TH Aspergillus nidulans stcL encodes a putative cytochrome P-450 monooxygenase required for bisfuran desaturation during aflatoxin/sterigmatocystin biosynthesis J. Biol. Chem. 1997 272 1589 1594 10.1074/jbc.272.3.1589 8999832
Kelkar, H. S., Skloss, T. W., Haw, J. F., Keller, N. P. & Adams, T. H. Aspergillus nidulans stcL encodes a putative cytochrome P-450 monooxygenase required for bisfuran desaturation during aflatoxin/sterigmatocystin biosynthesis. J. Biol. Chem. 272, 1589–1594 (1997).8999832 10.1074/jbc.272.3.1589
36. Yu J Current understanding on aflatoxin biosynthesis and future perspective in reducing aflatoxin contamination Toxins 2012 4 1024 1057 10.3390/toxins4111024 23202305
Yu, J. Current understanding on aflatoxin biosynthesis and future perspective in reducing aflatoxin contamination. Toxins 4, 1024–1057 (2012).23202305 10.3390/toxins4111024
37. Uka V Secondary metabolite dereplication and phylogenetic analysis identify various emerging mycotoxins and reveal the high intra-species diversity in Aspergillus flavus Front. Microbiol. 2019 10 667 10.3389/fmicb.2019.00667 31024476
Uka, V. et al. Secondary metabolite dereplication and phylogenetic analysis identify various emerging mycotoxins and reveal the high intra-species diversity in Aspergillus flavus. Front. Microbiol. 10, 667 (2019).31024476 10.3389/fmicb.2019.00667
38. Xie H Monitoring metabolite production of aflatoxin biosynthesis by Orbitrap Fusion mass spectrometry and a D-optimal mixture design method Anal. Chem. 2018 90 14331 14338 10.1021/acs.analchem.8b03703 30444348
Xie, H. et al. Monitoring metabolite production of aflatoxin biosynthesis by Orbitrap Fusion mass spectrometry and a D-optimal mixture design method. Anal. Chem. 90, 14331–14338 (2018).30444348 10.1021/acs.analchem.8b03703
39. Sweany RR Damann KE Kaller MD Comparison of soil and corn kernel Aspergillus flavus populations: evidence for niche specialization Phytopathology 2011 101 952 959 10.1094/PHYTO-09-10-0243 21405994
Sweany, R. R., Damann, K. E. & Kaller, M. D. Comparison of soil and corn kernel Aspergillus flavus populations: evidence for niche specialization. Phytopathology 101, 952–959 (2011).21405994 10.1094/PHYTO-09-10-0243
40. Ford RE Foster GD Bailey AM Exploring fungal RiPPs from the perspective of chemical ecology Fungal Biol. Biotechnol. 2022 9 12 10.1186/s40694-022-00144-9 35752794
Ford, R. E., Foster, G. D. & Bailey, A. M. Exploring fungal RiPPs from the perspective of chemical ecology. Fungal Biol. Biotechnol. 9, 12 (2022).35752794 10.1186/s40694-022-00144-9
41. Chalivendra SC DeRobertis C Chang P-K Damann KE Cyclopiazonic acid is a pathogenicity factor for Aspergillus flavus and a promising target for screening germplasm for ear rot resistance Mol. Plant-Microbe Interact. 2017 30 361 373 10.1094/MPMI-02-17-0026-R 28447887
Chalivendra, S. C., DeRobertis, C., Chang, P.-K. & Damann, K. E. Cyclopiazonic acid is a pathogenicity factor for Aspergillus flavus and a promising target for screening germplasm for ear rot resistance. Mol. Plant-Microbe Interact. 30, 361–373 (2017).28447887 10.1094/MPMI-02-17-0026-R
42. Singh P Orbach MJ Cotty PJ Aspergillus texensis: a novel aflatoxin producer with S morphology from the United States Toxins 2018 10 513 10.3390/toxins10120513 30513994
Singh, P., Orbach, M. J. & Cotty, P. J. Aspergillus texensis: a novel aflatoxin producer with S morphology from the United States. Toxins 10, 513 (2018).30513994 10.3390/toxins10120513
43. Misslinger M Hortschansky P Brakhage AA Haas H Fungal iron homeostasis with a focus on Aspergillus fumigatus Biochim. Biophys. Acta – Mol. Cell Res. 2021 1868 118885 10.1016/j.bbamcr.2020.118885 33045305
Misslinger, M., Hortschansky, P., Brakhage, A. A. & Haas, H. Fungal iron homeostasis with a focus on Aspergillus fumigatus. Biochim. Biophys. Acta – Mol. Cell Res. 1868, 118885 (2021).33045305 10.1016/j.bbamcr.2020.118885
44. Amare MG Keller NP Molecular mechanisms of Aspergillus flavus secondary metabolism and development Fungal Genet. Biol. 2014 66 11 18 10.1016/j.fgb.2014.02.008 24613992
Amare, M. G. & Keller, N. P. Molecular mechanisms of Aspergillus flavus secondary metabolism and development. Fungal Genet. Biol. 66, 11–18 (2014).24613992 10.1016/j.fgb.2014.02.008
45. Singh P Callicott KA Orbach MJ Cotty PJ Molecular analysis of S-morphology aflatoxin producers from the United States reveals previously unknown diversity and two new taxa Front. Microbiol. 2020 11 1236 10.3389/fmicb.2020.01236 32625180
Singh, P., Callicott, K. A., Orbach, M. J. & Cotty, P. J. Molecular analysis of S-morphology aflatoxin producers from the United States reveals previously unknown diversity and two new taxa. Front. Microbiol. 11, 1236 (2020).32625180 10.3389/fmicb.2020.01236
46. Kensler TW Roebuck BD Wogan GN Groopman JD Aflatoxin: a 50-year odyssey of mechanistic and translational toxicology Toxicol. Sci. 2011 120 S28 S48 10.1093/toxsci/kfq283 20881231
Kensler, T. W., Roebuck, B. D., Wogan, G. N. & Groopman, J. D. Aflatoxin: a 50-year odyssey of mechanistic and translational toxicology. Toxicol. Sci. 120, S28–S48 (2011).20881231 10.1093/toxsci/kfq283
47. Matumba L Uncommon occurrence ratios of aflatoxin B1, B2, G1, and G2 in maize and groundnuts from Malawi Mycotoxin Res. 2015 31 57 62 10.1007/s12550-014-0209-z 25194830
Matumba, L. et al. Uncommon occurrence ratios of aflatoxin B1, B2, G1, and G2 in maize and groundnuts from Malawi. Mycotoxin Res. 31, 57–62 (2015).25194830 10.1007/s12550-014-0209-z
48. Ting WTE Chang CH Szonyi B Gizachew D Growth and aflatoxin B1, B2, G1, and G2 production by Aspergillus flavus and Aspergillus parasiticus on ground flax seeds (Linum usitatissimum) J. Food Prot. 2020 83 975 983 10.4315/JFP-19-539 32034398
Ting, W. T. E., Chang, C. H., Szonyi, B. & Gizachew, D. Growth and aflatoxin B1, B2, G1, and G2 production by Aspergillus flavus and Aspergillus parasiticus on ground flax seeds (Linum usitatissimum). J. Food Prot. 83, 975–983 (2020).32034398 10.4315/JFP-19-539
49. Drott MT Lazzaro BP Brown DL Carbone I Milgroom MG Balancing selection for aflatoxin in Aspergillus flavus is maintained through interference competition with, and fungivory by insects Proc. Biol. Sci. 2017 284 20172408 29263278
Drott, M. T., Lazzaro, B. P., Brown, D. L., Carbone, I. & Milgroom, M. G. Balancing selection for aflatoxin in Aspergillus flavus is maintained through interference competition with, and fungivory by insects. Proc. Biol. Sci. 284, 20172408 (2017).29263278
50. Kenne GJ Activation of aflatoxin biosynthesis alleviates total ROS in Aspergillus parasiticus Toxins 2018 10 57 10.3390/toxins10020057 29382166
Kenne, G. J. et al. Activation of aflatoxin biosynthesis alleviates total ROS in Aspergillus parasiticus. Toxins 10, 57 (2018).29382166 10.3390/toxins10020057
51. Bankevich A SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing J. Comput. Biol. 2012 19 455 477 10.1089/cmb.2012.0021 22506599
Bankevich, A. et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477 (2012).22506599 10.1089/cmb.2012.0021
52. Gurevich A Saveliev V Vyahhi N Tesler G QUAST: quality assessment tool for genome assemblies Bioinformatics 2013 29 1072 1075 10.1093/bioinformatics/btt086 23422339
Gurevich, A., Saveliev, V., Vyahhi, N. & Tesler, G. QUAST: quality assessment tool for genome assemblies. Bioinformatics 29, 1072–1075 (2013).23422339 10.1093/bioinformatics/btt086
53. Kolmogorov M Yuan J Lin Y Pevzner PA Assembly of long, error-prone reads using repeat graphs Nat. Biotechnol. 2019 37 540 546 10.1038/s41587-019-0072-8 30936562
Kolmogorov, M., Yuan, J., Lin, Y. & Pevzner, P. A. Assembly of long, error-prone reads using repeat graphs. Nat. Biotechnol. 37, 540–546 (2019).30936562 10.1038/s41587-019-0072-8
54. Walker BJ Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement PLoS ONE 2014 9 e112963 10.1371/journal.pone.0112963 25409509
Walker, B. J. et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE 9, e112963 (2014).25409509 10.1371/journal.pone.0112963
55. Li H Durbin R Fast and accurate short read alignment with Burrows-Wheeler transform Bioinformatics 2009 25 1754 1760 10.1093/bioinformatics/btp324 19451168
Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760 (2009).19451168 10.1093/bioinformatics/btp324
56. Li H The sequence alignment/map format and SAMtools Bioinformatics 2009 25 2078 2079 10.1093/bioinformatics/btp352 19505943
Li, H. et al. The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).19505943 10.1093/bioinformatics/btp352
57. Okonechnikov K Conesa A García-Alcalde F Qualimap 2: advanced multi-sample quality control for high-throughput sequencing data Bioinformatics 2016 32 292 294 10.1093/bioinformatics/btv566 26428292
Okonechnikov, K., Conesa, A. & García-Alcalde, F. Qualimap 2: advanced multi-sample quality control for high-throughput sequencing data. Bioinformatics 32, 292–294 (2016).26428292 10.1093/bioinformatics/btv566
58. Smit, A. F., Hubley, R. & Green, P. RepeatMasker Open-4.0. Available online: http://repeatmasker.org/ (accessed on 25 June 2023).
59. Hubley R The Dfam database of repetitive DNA families Nucleic Acids Res. 2016 44 D81 D89 10.1093/nar/gkv1272 26612867
Hubley, R. et al. The Dfam database of repetitive DNA families. Nucleic Acids Res. 44, D81–D89 (2016).26612867 10.1093/nar/gkv1272
60. Hoff KJ Stanke M Predicting genes in single genomes with AUGUSTUS Curr. Protoc. Bioinform. 2019 65 e57 10.1002/cpbi.57
Hoff, K. J. & Stanke, M. Predicting genes in single genomes with AUGUSTUS. Curr. Protoc. Bioinform. 65, e57 (2019).10.1002/cpbi.57
61. Götz S High-throughput functional annotation and data mining with the Blast2GO suite Nucleic Acids Res. 2008 36 3420 3435 10.1093/nar/gkn176 18445632
Götz, S. et al. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 36, 3420–3435 (2008).18445632 10.1093/nar/gkn176
62. Jones P InterProScan 5: genome-scale protein function classification Bioinformatics 2014 30 1236 1240 10.1093/bioinformatics/btu031 24451626
Jones, P. et al. InterProScan 5: genome-scale protein function classification. Bioinformatics 30, 1236–1240 (2014).24451626 10.1093/bioinformatics/btu031
63. Huerta-Cepas J eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses Nucleic Acids Res. 2019 47 D309 D314 10.1093/nar/gky1085 30418610
Huerta-Cepas, J. et al. eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 47, D309–D314 (2019).30418610 10.1093/nar/gky1085
64. Manni M Berkeley MR Seppey M Simão FA Zdobnov EM BUSCO update: novel and streamlined workflows along with broader and deeper phylogenetic coverage for scoring of eukaryotic, prokaryotic, and viral genomes Mol. Biol. Evol. 2021 38 4647 4654 10.1093/molbev/msab199 34320186
Manni, M., Berkeley, M. R., Seppey, M., Simão, F. A. & Zdobnov, E. M. BUSCO update: novel and streamlined workflows along with broader and deeper phylogenetic coverage for scoring of eukaryotic, prokaryotic, and viral genomes. Mol. Biol. Evol. 38, 4647–4654 (2021).34320186 10.1093/molbev/msab199
65. Al-Shahrour F Díaz-Uriarte R Dopazo J FatiGO: a web tool for finding significant associations of Gene Ontology terms with groups of genes Bioinformatics 2004 20 578 580 10.1093/bioinformatics/btg455 14990455
Al-Shahrour, F., Díaz-Uriarte, R. & Dopazo, J. FatiGO: a web tool for finding significant associations of Gene Ontology terms with groups of genes. Bioinformatics 20, 578–580 (2004).14990455 10.1093/bioinformatics/btg455
66. Medema MH antiSMASH: rapid identification, annotation and analysis of secondary metabolite biosynthesis gene clusters in bacterial and fungal genome sequences Nucleic Acids Res. 2011 39 W339 W346 10.1093/nar/gkr466 21672958
Medema, M. H. et al. antiSMASH: rapid identification, annotation and analysis of secondary metabolite biosynthesis gene clusters in bacterial and fungal genome sequences. Nucleic Acids Res. 39, W339–W346 (2011).21672958 10.1093/nar/gkr466
67. Blin K antiSMASH 6.0: improving cluster detection and comparison capabilities Nucleic Acids Res. 2021 49 W29 W35 10.1093/nar/gkab335 33978755
Blin, K. et al. antiSMASH 6.0: improving cluster detection and comparison capabilities. Nucleic Acids Res. 49, W29–W35 (2021).33978755 10.1093/nar/gkab335
68. Pitt, J. I. The Genus Penicillium And Its Teleomorphic States Eupenicillium And Talaromyces. 1st ed. (AP, London, 1979).
69. Hocking, A. D. et al. Selective media and procedures. in Methods for the mycological examination of food, 1st edn. (eds A. D. King, J.I. Pitt, L. R. Beuchat & J. Corry) (Springer, 1986).
70. Deutsches Institut für Normung e. V. Foodstuff - Determination of aflatoxin B 1 , and the total content of aflatoxins B 1 , B 2 , G 1 and G 2 in cereals, nuts and derived products - High performance liquid chromatographic method (ISO 16050:2003) - German version EN ISO 16050:2011. 10.31030/1807322 (Beuth Verlag GmbH, Berlin, 2011).
71. Schamann A Schmidt-Heydt M Geisen R Kulling SE Soukup ST Formation of B- and M-group aflatoxins and precursors by Aspergillus flavus on maize and its implication for food safety Mycotoxin Res. 2022 38 79 92 10.1007/s12550-022-00452-4 35288866
Schamann, A., Schmidt-Heydt, M., Geisen, R., Kulling, S. E. & Soukup, S. T. Formation of B- and M-group aflatoxins and precursors by Aspergillus flavus on maize and its implication for food safety. Mycotoxin Res. 38, 79–92 (2022).35288866 10.1007/s12550-022-00452-4
