Global population genomics redefines domestication and clinical diversity in the Aspergillus flavus–oryzae complex

1 Global population genomics redefines domestication and clinical diversity in the Aspergillus flavus–oryzae complex Walter P. Pfliegler1, Bálint Németh1,2 , Veronika Bodnár1,2 , Tünde Pusztahelyi3, Ignazio Carbone4, István Pócsi1,5 1 Department of Molecular Biotechnology and Microbiology, Institute of Biotechnology, Faculty of Science and Technology, University of Debrecen, H-4032, Egyetem tér 1, Debrecen, Hungary 2 Doctoral School of Nutrition and Food Sciences, University of Debrecen, H-4032, Egyetem tér 1, Debrecen, Hungary 3 Central Laboratory of Agricultural and Food Products, Faculty of Agricultural and Food Sciences and Environmental Management, University of Debrecen; H-4032 Böszörményi str. 138. Debrecen, Hungary 4 Center for Integrated Fungal Research, Department of Entomology and Plant Pathology, North Carolina State University, Raleigh, NC 27695, USA 5 FoodChainSafetyLaboratoryDirectorate,NationalFoodChainSafetyOffice,H-1095Mesterstr.81,Budapest,Hungary Corresponding author: Walter P. Pfliegler (pfliegler[email protected]) Copyright: This is an open access article distributed under the terms of the CC0 Public Domain Dedication. Research Article Abstract Aspergillus flavus is a globally important human pathogen and agricultural contaminant, while its domesticated relative A. oryzae is widely used in food fermentation and biotechnology. Despite their importance, the evolutionary relationship, population structure and domestication history of these fungi remain unresolved. Here, we present the first global population genomic analysis of 639 A. flavus and A. oryzae isolates from clinical, environmental and food-fermentation sources across multiple continents. Our analyses reveal a complex evolutionary landscape comprising wellseparated clades interspersed with highly admixed mosaic groups and potential evidence for multiple independent domestication events giving rise to A. oryzae. Clinical A. flavus isolates are distributed across several clades and mosaic groups, some overlapping with fermentation strains, highlighting an apparent role of domestication and admixture in shaping pathogen diversity. These results challenge current species boundaries and provide a framework for understanding evolutionary history, taxonomy and pangenomic architecture in these fungi, with broad implications for pathogenicity, food safety, biocontrol and metagenomic surveillance. Key words: Aspergillosis, biocontrol, microbe domestication, mycobiome, plant pathogen, traditional fermentation Introduction Fungi play diverse, significant and sometimes contrasting roles in agriculture, the food-chain and human health and disease (Hyde et al. 2018; Han et al. 2024; Hill and Round 2024). Despite the immense diversity of the Fungal Kingdom, it is quite rare for a single genus to be both widely used in food production and regarded as a major pathogen of plants and humans. The species-rich Aspergillus genus (Houbraken et al. 2020) is probably the most widely known and Academic editor: Janneke Aylward Received: 17 September 2025 Accepted: 4 December 2025 Published: 23 December 2025 Citation: Pfliegler WP, Németh B, Bodnár V, Pusztahelyi T, Carbone I, Pócsi I (2025) Global population genomics redefines domestication and clinical diversity in the Aspergillus flavus–oryzae complex. IMA Fungus 16: e172343. https://doi.org/10.3897/ imafungus.16.172343 IMA Fungus 16: e172343 (2025) DOI: 10.3897/imafungus.16.172343 2 IMA Fungus 16: e172343 (2025), DOI: 10.3897/imafungus.16.172343 Walter P. Pfliegler et al.: Aspergillus flavus and oryzae global phylogenomics studied example of such a fungal taxon: its so-called section Flavi (Han et al. 2024) includes, amongst others, two widely known, closely-related and often hard-to-distinguish taxa: A. flavus, described by Link in the beginning of the 19th century (Link 1809) and A. oryzae, described at the end of the 19th century (Korschelt 1878; Cohn 1884), originally as Eurotium oryzae. The species A. flavus is an important crop and post-harvest pathogen with wide-ranging effects on the global food chain, not just due to yield loss, but as the producer of the aflatoxins B1, B2 and G (AFs), indol-tetramic acid and cyclopiazonic acid (CPA) (Horn 2003). Aflatoxin is significant due to its high toxicity and potent hepatocarcinogenic effects (Peles et al. 2019; Rushing and Selim 2019). A. flavus is also the second most important Aspergillus species causing human infections, following A. fumigatus (Dabas et al. 2022; Kanaujia et al. 2023). However, its non-aflatoxinogenic strains have been found to be effective biocontrol agents against aflatoxin-producing strains (Khan et al. 2021; Molo et al. 2022). The most well-studied such biocontrol agents are Afla-Guard (strain NRRL 21882) and AF36 Prevail (AF36, strain NRRL 18543), used in North America, that competitively interfere with native aflatoxigenic strains (Dorner 2004), while new strains and mixes of strains are also constantly being tested in various geographical settings including Asian (e.g. Rasheed et al. (2024); Xu et al. (2025)), African (e.g. Augusto et al. (2024)) and European countries (Alvarez et al. 2022). On the other hand, A. oryzae is an important microbe in traditional and modern food technology, used in the solid-state fermentation steps of producing soy sauces, fermented soybean pastes, rice wines, baiju and various rice vinegars (Bal et al. 2016; Yun et al. 2020; Chin et al. 2024; Peng et al. 2025), taking advantage of its high hydrolytic enzyme activities. Its most common colloquial name is the koji fungus, originating from the yellow kōji of Japanese sake brewing. Its use has recently expanded into biocontrol of both insects and toxigenic fungi (e.g. Alshannaq et al. (2018); You et al. (2024)), postbiotics and animal feed supplements (Seidler et al. 2024) and commercial enzyme production, recently by invoking the toolset of synthetic biology as well (Sun et al. 2024). These two closely-related species and the problems of their delimitation have generated interest since the advent of molecular genetic techniques. Based on early molecular genetic findings, the koji fungus was reclassified as a variety of A. flavus (Kurtzman et al. 1986). The single-species concept for the two fungi was not accepted in a recent large-scale revision of section Flavi by Frisvad et al. (2019). However, the authors note that A. oryzae and A. sojae appear to be domesticated forms of the aflatoxigenic species A. flavus and A. parasiticus, respectively. In line with this, the origin of A. oryzae from within the species A. flavus through domestication has generally been supported by the literature of the past decades, yet the species status of both is upheld (Chang et al. 2006; Gibbons et al. 2012; Watarai et al. 2019; Han et al. 2024; Jeong and Seo 2025). The regulatory confusion that conspecificity would generate is often cited as a pragmatic reason to maintain their separate species status (Geiser et al. 1998; Houbraken et al. 2014). A few analyses have challenged the view of an extremely close relationship between the two “species”. For example, Kjærbølling et al. identified A. minisclerotigenes and A. aflatoxiformans as the closest relatives of the koji fungus (Kjærbølling et al. 2020). 3 IMA Fungus 16: e172343 (2025), DOI: 10.3897/imafungus.16.172343 Walter P. Pfliegler et al.: Aspergillus flavus and oryzae global phylogenomics The genome of A. oryzae was one of the first in the genus to have a full assembly (Machida et al. 2005), while a chromosome-level assembly of A. flavus strain NRRL 3357 has recently been published (Skerker et al. 2021) as well. Recent advances in genome sequencing technology have enabled the sequencing and analysis of hundreds of strains and isolates, resulting in significant insights into the population structure and genome evolution of A. flavus and A. oryzae (Moore et al. 2017; Gell et al. 2020). In the case of A. flavus, isolates were first distinguished by sclerotial morphotypes into the S (small sclerotia; should not be confused with another species characterised by small sclerotia, namely A. minisclerotigenes) and L (large sclerotia, with high variation on toxigenicity) groups and some small sclerotial isolates from Thailand were described as separate variant of the species, as A. flavus var. parvisclerotigenus (Saito and Tsuruta 1993). Later, based on phylogenetics of a limited number of sequenced genes, groups I (Sand L-strains that produce AFB1) and II (S-strains that produce AFB1 and/or G-aflatoxins) were distinguished (Geiser et al. 1998). Importantly, A. oryzae strains fell into group I. Group I was subsequently divided into three lineages: IA, IB and IC (Geiser et al. 2000). Group IB was associated with A. oryzae strains, while the sclerotial morphotypes did not cluster with the newly-established genetic groups, as both group IA and group II included strains with the S morphotype. Later works using phylogenetic methods revealed the importance of recombination events amongst the lineages and chemotypes (aflatoxigenic and non-aflatoxigenic types) of group I (Moore et al. 2009, 2013). Moore et al. (2017) examined a worldwide collection of 1,304 A. flavus, A. oryzae, A. nomius, A. parasiticus and A. sojae isolates with multi-locus sequence typing and identified seven distinct genetic clusters. A. flavus was grouped into three clusters: L-type (group I), SB and SBG (group II). The latter two represent distinct chemotypes of S-type isolates, producing only aflatoxin B (SB) or both aflatoxins B and G (SBG). The SBG cluster was proposed to be of hybrid origin, containing haplotypes from both A. flavus and A. parasiticus. The first study to compare three S-type and three L-type genomes of the species also found no clear phylogenetic separation between the two morphotypes (Ohkura et al. 2018). In a study of 94 A. flavus whole-genome sequences, primarily from the continental United States, Drott et al. (2020) identified three genetically distinct populations of L-type isolates (A, B and C), based on phylogenomic networks, along with a fourth, then-unnamed group comprising two S-type isolates. Population A was characterised by relatively frequent sexual processes, based on linkage disequilibrium decay analysis and higher aflatoxin production. Population B had variable aflatoxigenicity and included non-aflatoxigenic isolates and was determined to be sympatric to population A. Population C was mainly non-aflatoxigenic (Drott et al. 2020). The genomes were later used in a study describing the secondary metabolic gene clusters in the pangenome of the species (Drott et al. 2021) and in a phylogeny for twelve newly-sequenced isolates from Pakistan (Ajmal et al. 2022). The latter work recovered the same three populations. Gangurde et al. (2024) inferred a phylogeny of 346 A. flavus isolates, including 225 newly-sequenced ones from the US and recovered the same populations, with a focus on the pangenome of the species. Most recently, 70 predominantly European clinical isolates were sequenced and phylogenomic analysis revealed a fourth population greatly enriched in clinical isolates, named Population D (Hatmaker et al. 2025). The group of small sclerotial isolates was also recovered 4 IMA Fungus 16: e172343 (2025), DOI: 10.3897/imafungus.16.172343 Walter P. Pfliegler et al.: Aspergillus flavus and oryzae global phylogenomics and designated as a population of its own, named S-type. A recent study by Molo et al. (2022) utilised a ddRADseq-based dataset of 815 samples, which were subjected to phylogenetic, admixture and cluster analyses. The authors reported that population A from Drott et al. (2020) largely corresponds to lineage IC, population B to lineage IB, while population C showed affinities to IC. Based on these results, two distinct evolutionary lineages, IB and IC, were delineated (lineages IA and II strains were not included in this analysis). Lineage IB was genetically more homogeneous, whereas lineage IC exhibited a highly reticulate pattern in the neighbour-net analysis and greater variation in the SNPbased principal component analysis. Notably, using 353 unique genome-wide haplotypes, Molo et al. (2022) also demonstrated that SNP-based analyses were largely unaffected by changing the reference genome from A. flavus NRRL3357 to the A. oryzae RIB40 assembly, providing further evidence of the extremely close relatedness between the two taxa. In the case of A. oryzae, eight main clades have been described and designated by letters A – H based on a large-scale analysis of 82 genomes (Watarai et al. 2019) with evidence of past recombination events amongst clades. The clades were also recovered in a subsequent analysis of 91 genomes (Chacón-Vargas et al. 2021). It is notable that the A. oryzae literature typically refers to “clades”, whereas A. flavus is more often described in terms of lineages or groups (IA, IB, IC and II) or populations (A, B, C, D and S-type) in literature, as outlined above. Most other members of section Flavi still lack comprehensive population-level genome sequencing studies (Kjærbølling et al. 2020). Despite the steadily increasing number of available assemblies and whole genome sequencing projects, only recently has a phylogenomic analysis been published that included a considerable number of genomes from both species. Han et al. (2024) highlighted the genomic similarities between the two Aspergillus species and recovered A. oryzae as a clade nested within the genetic diversity of A. flavus. However, the phylogeny, based on amino acid sequences of all single-copy orthologous genes in 93 assemblies, did not include the A. flavus populations of Drott et al. (2020) or the previously described lineages I and II. Despite the recent advancements in the phylogenomics of these Aspergillus species, clades of A. oryzae and populations and lineages of A. flavus have not been directly compared. Additionally, the use of the A. flavus and/or A. oryzae reference genomes is not consistent in the literature. There is, thus, a need for more robust taxonomy of these important fungi as proposed earlier (Arias et al. 2021). Our aim in this study was to collect available short-read sequencing data on these two species, summarise information on the strains and isolates they represent and subject them to phylogenomic and population genomic analysis. Our focus in this is on the phylogenetic relationships of clades and potential mosaic groups, and their geographic distribution. Materials and methods Overview of taxonomy, literature and databases To place phylogenomic and population genomic findings into both a taxonomic and applied context, we reviewed: (1) literature on the taxonomic descriptions and type strains along with described infraspecific taxa of the two species; 5 IMA Fungus 16: e172343 (2025), DOI: 10.3897/imafungus.16.172343 Walter P. Pfliegler et al.: Aspergillus flavus and oryzae global phylogenomics (2) their population structures, clades, lineages and assemblies; and (3) as well as published BioProjects with sequences. We also queried major fungal taxonomic databases (Index and Species Fungorum, Mycobank) and databases frequently used in molecular genetic species identification (NCBI Taxonomy Database, FungiDB, SILVA or UNITE) to check the status of A. flavus and A. oryzae and infraspecific taxa associated with them. Novel isolates and isolation of genomic DNA Genomic DNA was isolated from three Hungarian isolates from the collection of the Central Laboratory of Agricultural and Food Products (Faculty of Agricultural and Food Sciences and Environmental Management, University of Debrecen) and our laboratory’s stock of the reference strain NRRL 3357 as a control strain for genomic analyses. Library preparation was performed using tagmentation with the Nextera DNA Flex Library Prep kit (Illumina, San Diego, CA, USA) according to the manufacturer’s protocol. Sequencing was performed using 150 bp paired-end reads on an Illumina NextSeq 500 system, with approximately 140× coverage of the nuclear genome. Raw reads were deposited in NCBI SRA under BioProject PRJNA1307175. Accession numbers are listed in Suppl. material 2. Published genomes A comprehensive list of published genomes sequenced using short-read technology was compiled for both species and additional genomes were identified by searching NCBI BioProjects with the species’ names (accessed 16 January 2025). Metadata associated with the sequencing files were collected from ENA and SRA and information on the isolates and strains was compiled from publications, including an ontology of habitats. Habitat categories included Agriculture (with subcategories of plant types or field soil), Environmental (e.g. soil samples where agricultural origin was not specifically stated), Human (clinical isolates) and Food industry (e.g. koji, meju, nuruk or qu). Illumina sequencing runs were downloaded from SRA. The Illumina FASTQ sequencing files were trimmed and filtered using fastp for further analysis (Chen et al. 2018). In the list of genomes, geographic areas were grouped into larger categories following the United Nations geoscheme (accessed 01 June 2025), while isolation sources were grouped into the following major categories: Environmental, Agriculture, Human, Food fermentation and Artificial Environment. Phylogenomic analysis Illumina reads were mapped to the ASM901741v1 (GenBank accession: GCA_009017415.1) A. flavus strain NRRL 3357 (Skerker et al. 2021) reference genome using BWA 0.7.17. (Li and Durbin 2010). As the reference contained no mitochondrial assembly, the complete mitochondrial sequence of the isolate CA14 (GenBank accession: CP061812.1) was appended to the file. It has been demonstrated that the two commonly used assemblies in section Flavi are interchangeable (Molo et al. 2022), thus, A. oryzae genomes were also mapped to this former A. flavus reference. Sorted BAM files from reference-based alignments were obtained using Samtools 1.7 (Li and Durbin 2010) and Picard-tools 6 IMA Fungus 16: e172343 (2025), DOI: 10.3897/imafungus.16.172343 Walter P. Pfliegler et al.: Aspergillus flavus and oryzae global phylogenomics 2.23.8 was used to mark duplicated reads (Van der Auwera et al. 2013). BAM files were subjected to local re-alignment around indels and joint variant calling using GATK 4.1.9.0. (Van der Auwera et al. 2013; Poplin et al. 2018) with ploidy set to 1. In the first step, genomic VCF files were obtained with GATK Haplotype Caller and joint genotyping of the gVCF files was applied. Subsequent filtering of resulting VCF files was based on either SNPs or INDELS. SNPs were filtered according to the following parameters: QD < 5.0; QUAL < 30.0; SOR > 3.0; FS > 60.0; MQ < 40.0; MQRankSum < -12.5; ReadPosRankSum < -8.0. INDELS were filtered according to the parameters QD < 5.0; QUAL < 30.0; FS > 60.0; ReadPosRankSum < -20.0. INDELS were then left-aligned. For the final VCF files, INDELS and SNPs were merged, filtered and non-variant sites were removed. Calling was performed for every chromosome and the mitochondrial genome separately. Two outgroup genomes were also included in the cohort variant calling workflow: A. minisclerotigenes E1406 and A. parasiticus AA-M17-157. Initial tests on the dataset (SNP calling as described here and neighbour-net graphs as described below) showed that a few A. flavus genomes from literature (Arias et al. 2020) clustered with A. minisclerotigenes and these were excluded from the final dataset (samples E1288, E1293, E1316, E1376). The initial analyses also showed that mapping the A. minisclerotigenes and A. parasiticus reads to the NRRL 3357 reference results in a high proportion of callable sites despite their divergence (gap and ambiguity: 9.36% and 14.22% of total 3,469,698 called variant sites, respectively). Thus, they are appropriate as outgroups in the SNP matrix-based analyses. The SNP VCF files from cohort calling were used to produce genotype matrices using vcf2phylip (Ortiz 2019). Matrices were used as input in the software SplitsTree 4.15.1. (Huson and Bryant 2006) to create a phylogenomic neighbour-net graph using uncorrected P-distance. This was also converted to a neighbour-joining tree and visualised in iTOL, along with additional data (e.g. population genomics, see below) (Letunic and Bork 2019). Linkage disequilibrium pruning was not applied and, therefore, all SNPs were considered for building the network. The neighbour-net graph was generated both with and without the outgroup species to improve the visibility of the resulting figures. Identity-by-state (IBS) analysis IBS analysis was carried out using SNP data with the R package “SNPRelate” 1.24.0 (Zheng et al. 2012) using the function “snpgdsIBSNum” with default settings, calculating the minor allele frequency and missing rate for each SNP over all the samples. The analysis was based on converting the input VCF file with the snpgdsVCF2GDS function and creating n-by-n matrices. Due to the method used in the package (.gds files are converted from .vcf files), only biallelic SNPs were kept. Since linkage disequilibrium pruning was not applied for this analysis, all 1,764,492 biallelic SNPs were considered. SNPRelate expects diploid calls, thus the genotypes in the vcf file were corrected as if they were diploid homozygous calls using awk. This way, two informative datasets were obtained: the number of identical genotypes between each pair of genomes (analogous to IBS2 of diploids) and the number of non-identical genotypes (analogous to IBS0 of diploids). Instead of applying clustering, the order of genomes was changed to fit that of the neighbour-joining tree. Results were visualised with the R package “heatmaply” (Galili et al. 2018) with default settings. Zoomable .html files 7 IMA Fungus 16: e172343 (2025), DOI: 10.3897/imafungus.16.172343 Walter P. Pfliegler et al.: Aspergillus flavus and oryzae global phylogenomics for heatmaps are uploaded to FigShare (doi: 10.6084/m9.figshare.29896676). After clades and the genomes belonging to them had been established, the table-format data generated by the IBS analysis was used again, to calculate SNP identity amongst all possible pairs of genomes in each pair of groups. These were visualized as violin plots for all possible group comparisons (StatisticsKingdom 2017). Outgroup species were excluded from this analysis. Population genomic analysis was performed using fastStructure 1.0 (Raj et al. 2014) run on SNPs filtered for biallelic variants with minor allele frequency greater than 5%. Linkage disequilibrium pruning was conducted with plink 1.9 (Chang et al. 2015) (plink -indep-pairwise 50 1 0.5). Populations were specified from 2 to 10. The chooseK command was used to identify the best number of populations. The population genomic analysis after pruning relied on 850,234 variants (i.e. ~ 48% of all biallelic SNPs). Population genomic data was visualied in iTOL as a dataset added to a dendrogram. Outgroup species were excluded from this analysis. A subset of the genomes was chosen for Maximum Likelihood analysis using IQTree 2.3.6. (Minh et al. 2020) and ModelFinder (Kalyaanamoorthy et al. 2017). From each clade and mosaic group, five representative genomes were chosen randomly and only these were kept in the genotype matrix. The matrix was used as an input for finding the best model and to obtain a dendrogram with 1000 ultrafast bootstrap approximations, which was visualised in iTOL. A single A. minisclerotigenes outgroup (E1406) was used to root the dendrogram. The above mentioned phylogenomic and population genomic analyses were run on the individual chromosomes and on the mitochondrial genome, as well as on a combined dataset as well, representing the whole genome. However, the A. oryzae and A. flavus reference genomes are not completely syntenic, therefore, relying on the chromosomal boundaries of the NRRL 3357 genome could introduce uncertainties. Thus, we limited detailed analyses to the combined dataset. Results Neotypes, reference genomes and the status of A. flavus and A. oryzae In this work, we first reviewed the original descriptions of A. flavus and A. oryzae. We also tracked how their taxonomic status changed over time, which neotype strains were assigned to them and which of these has an available genome sequence. A. flavus Link was described in the beginning of the 19th century (Link 1809). The first mention of A. oryzae is as Eurotium Oryzae Ahlburg [sic] in a publication on Japanese sake-brewing dating back to 1878 (Korschelt 1878). In this publication, Korschelt scientifically described the koji fungus in the name and in honour of Ahlburg (Korschelt 1878). Years later, it was later reclassified as A. Oryzae [sic] by Cohn (1884), its full botanical name becoming A. oryzae (Ahlb.) Cohn. Several variations and forms of A. oryzae have been described in the 20th century and many of them have been subsequently synonymised (Raper and Fennell 1973). Mycobank lists 18 such names (Mycobank, accessed 16.01.2025). It is worth noting that the publication of A. oryzae var. basidiferus Costantin & Lucet in the early 20th century (Costantin and Lucet 1905) resulted in the previously described original koji fungus being recognised as A. oryzae var. oryzae (Ahlb.) Cohn, in accordance with the rules of botanical nomenclature. 8 IMA Fungus 16: e172343 (2025), DOI: 10.3897/imafungus.16.172343 Walter P. Pfliegler et al.: Aspergillus flavus and oryzae global phylogenomics Fungal nomenclature relies on type strains or herbarium type material deposited in collections and these types must be taken into account when discussing the taxonomic status of the two fungi. Due to A. flavus and A. oryzae being first described in the 19th century, original type strains for the species A. flavus and A. oryzae have not been designated and deposited, further complicating their taxonomic status as only neotypes are available. The neotype of A. flavus is available in multiple collections (NRRL 1957, ATCC 16883, CBS 569.65, CBS 100927, IMI 124930). Interestingly, this neotype is far from a typical representative of the species: it was collected from the cellophane diaphragm at an unknown location in the South Pacific (ARS Culture Collection, accessed 16.01.2025). Its genome was recently sequenced, but it is not considered the commonly used reference. The widely adopted reference genome is NRRL 3357 (Skerker et al. 2021), which includes complete chromosomes, but lacks an assembled mitochondrial genome. A neotype for the taxon E. oryzae Ahlb., collected in Japan, was designated as IMI 16266 (also deposited as CBS 102.07 and CBS 110.47). Furthermore, another type strain was designated for A. oryzae var. oryzae (Ahlb.) Cohn, which was collected in 1969 in Japan and currently deposited in multiple collections (ATCC 42149, CBS 466.91, IFO 6215, RIB40, WB 5032). The assembled genome of this strain is widely used as the reference genome for A. oryzae (Machida et al. 2005). Major taxonomic databases now regard E. oryzae Ahlb. and A. oryzae Cohn as a variety of A. flavus Link. Its full botanical name is thus listed as A. flavus var. oryzae (Ahlb.) Kurtzman, M.J. Smiley, Robnett & Wicklow. Whenever this infra-specific status is accepted, the non-oryzae A. flavus taxon has to be specified as A. flavus var. flavus. In contrast, the taxon A. oryzae var. oryzae (Ahlb.) Cohn 1884 is accepted as a legitimate and the current name for the koji fungus as it has never been synonymised in literature (Mycobank and Index Fungorum, both accessed 16.01.2025). The NCBI Taxonomic Database and, hence, the NCBI RefSeq database often used for mapping-based metagenomic analysis (e.g. Yan et al. (2024)) (accessed 15.03.2025) accepts the two species as separate, while UNITE (accessed 15.03.2025) lists only A. flavus barcode sequences. The FungiDB database (Basenko et al. 2024) and the SILVA ribosomal RNA reference database (Glöckner et al. 2017) (both accessed 15.03.2025) lists both species separately. Combined phylogenetic and population genomics of 639 A. oryzae and A. flavus genomes re-defines clades, mosaic groups and the origin of the koji fungus We identified, tabulated and downloaded all available A. oryzae and A. flavus short-read genome sequences from previously published papers and additional BioProjects containing publicly accessible sequencing files (NCBI SRA, accessed 16 January 2025). Three additional isolates from Hungary were also sequenced in this study and added to the list, which contains 639 genomes from 35 studies and Bioprojects. A small number of genomes above this number uploaded into SRA as S-type A. flavus clustered with A. minisclerotigenes (E1288, E1316, E1293, E1376, E1406) (Arias et al. 2020) in a preliminary phylogeny. These were discarded from further analyses, except for E1406 which was used 9 IMA Fungus 16: e172343 (2025), DOI: 10.3897/imafungus.16.172343 Walter P. Pfliegler et al.: Aspergillus flavus and oryzae global phylogenomics as an outgroup, along with the A. parasiticus strain AA-M17-157 (Hatmaker et al. 2024, 2025). A strain in SRA described as A. parasiticus NRRL 2999 proved to be identical to A. flavus NRRL 33757, as shown previously (Fountain et al. 2020; Houbraken et al. 2021) and was included with an identifier ‘NRR33572999’. A limited number of strains were sequenced by more than one study due to their importance in biocontrol or as a reference genome (Suppl. material 2). Thus, the total amount of unique strains and isolates (excluding outgroups) in the collection presented in this study from the past 20 years amounts to 632, collected from 17 countries and regions. The 639 genomes, including the replicates, were included in the variant cohort calling pipeline and in subsequent phylogenomic and population genomic analyses. Based on an SNP-matrix containing 1,667,702 parsimony-informative and 1,801,996 singleton sites (as determined by running IQTree on the whole genotype matrix) from all called variant sites of the genome, we first generated a distance matrix of the genomes and, subsequently, a neighbour-net graph based on this, following earlier studies (Drott et al. 2020; Molo et al. 2022). Using the NRRL 3357 reference genome for all genomes and the approach of reference-based variant calling resulted in a very high proportion of comparable variants. Only four samples had more than 10% gaps and unknown variants across the dataset and the mean percentage of such undetermined sites was 5.45%. Our dataset recovered groups within both A. flavus and A. oryzae: the “populations” described for A. flavus and the clades described for A. oryzae. Inclusion of both fungi in a single phylogenomic network enabled us to compare the diversity within groups of the two species. The number of described clades for the koji fungus proved to be too high to enable meaningful comparisons with the more diverse “populations” of A. flavus. The described A. oryzae clades (Watarai et al. 2019; Chacón-Vargas et al. 2021) were thus grouped into two superclades, from here on referred to as clade Oryzae A-B (former clades A and B) and clade Oryzae C-H (former clades C to H) (Watarai et al. 2019; Chacón-Vargas et al. 2021) (Fig. 1). For the sake of consistency, the “populations” of A. flavus are also referred to as clades in our analysis in cases when population genomic analysis does not point to mosaic genomic origin. The S-type “population” of A. flavus, from here on referred to Flavus S clade, had the fewest genomes and formed a distinct group connected by long parallel edges to the rest of the network. This included the genome AF70 and, thus, corresponds to the lineage IA of A. flavus. The “population B”, from here on referred to as the Flavus B clade, was also connected by long parallel edges to the rest of the groups and was characteristically compact, suggesting lower diversity. The Flavus B encompassed a subgroup which consisted of the Afla-Guard biocontrol strain and many highly similar, almost indistinguishable field isolates (Fig. 1). This group thus corresponds to lineage IB, based on the genomes clustering into it that are classified into lineage IB in literature. The “population A”, named here Flavus A clade, was characterised by higher diversity, unclear within-clade connections and formed a broad sub-network. This included the AF36 biocontrol strain and several extremely similar field isolates connected to the Flavus A subnetwork by long parallel edges (thus, clear separation and unclear within-group relations due to very high similarity). Based on the list of genomes clustering to this Flavus A clade, it largely corresponds to lineage IC in literature. 16 IMA Fungus 16: e172343 (2025), DOI: 10.3897/imafungus.16.172343 Walter P. Pfliegler et al.: Aspergillus flavus and oryzae global phylogenomics Many isolates (382) originated in northern America. This region was dominated by Flavus A and Flavus B genomes, along with large number of isolates in Flavus C and D clades. Eastern Asian samples predominantly originated from food and beverage fermentations and belonged to the two Oryzae clades, as well as the Flavus A lineage. Western European samples predominantly represented the Flavus D and Flavus A clades, the two mosaic groups (mostly in the Figure 5. Analysis of non-identical and identical SNP alleles in haploid calls in a heatmap format amongst the genomes. Samples are listed according to the NJ dendrogram (shown on the left and top edges, along with population genomic data with k = 10). Heatmap scales are on the right side. Four genomes with a high proportion of missing data are also included, resulting in uniform lines on the heatmaps. a number of non-identical variants b number of identical variants. Note the high divergence of the B clade resulting in a higher number of non-identical SNPs compared to all other clades and the high uniformity in A. orzyae clades. Figure 6. a Isolation sources of the 639 genomes and their placement into the various clades and mosaic groups b Geographic origin of the 639 genomes and their placement into the various clades and mosaic groups. Total number of genomes in each source type or geographic region is indicated above the stacked bar charts. Clades and mosaic groups are colour-coded. Note that the Human source contains 11 isogenic Flavus A genomes collected from a single patient (Suppl. material 2). 17 IMA Fungus 16: e172343 (2025), DOI: 10.3897/imafungus.16.172343 Walter P. Pfliegler et al.: Aspergillus flavus and oryzae global phylogenomics form of human samples), while a smaller number of Flavus B and Oryzae C-H samples were also found amongst them (in the latter case, these were human samples). Eastern Africa and southern Asia were represented by substantially smaller sample sizes and were dominated by Flavus C and Flavus A and by Flavus A samples, respectively. Discussion In this study, we provide a phylogenomic and population genomic overview using all available short-read genome sequence data for A. flavus and A. oryzae. We applied an approach relying on reference-based variant-calling that enables comparing their genomes even in cases where high-quality assemblies are not available. The 639 genomes included in the current dataset significantly exceed the sample sizes used in previous genome-wide SNP phylogenies. For comparison, Drott et al. (2020) examined 94 genomes, Gangurde et al. (2024) included 346 and Hatmaker et al. (2025) examined 300. However, the current dataset was smaller than the 815 genomes analysed by Molo et al. (2022) using ddRADSeq, which samples only a fraction of the genome. That study was also focused on a more restricted geographic area. Similar to Molo et al. (2022), our analyses showed that the A. flavus and A. oryzae reference genomes are interchangeable for population and phylogenomic studies using a reference-guided variant discovery pipeline. This fact has important implications in metagenomics as well. Most mycobiome studies perform read mapping against a collection of fungal reference genomes, often using combined reference datasets available from NCBI. This contains several A. flavus and A. oryzae references and, if no similarity-based clustering method is applied prior to read mapping (e.g. as applied in Yan et al. (2024)), reads from metagenomic libraries may map to both. This would indicate the simultaneous presence of both lineages even if only one is present. Similarly, several databases used for ITS metabarcoding were found to contain both A. flavus and A. oryzae. This may also result in unreliable mapping of amplicons and the simultaneous identification of both species. This may cause important uncertainties, especially in cases where the presence of both is not unexpected. For example, Peng et al. (2025) identified both species in some traditional solid-state fermentations, based on high-throughput amplicon sequencing. Our comprehensive database on available sequenced strains and isolates of both species along with their collection metadata (Suppl. material 1) aimed to provide information for further studies, facilitating large-scale phylogenies and the correct placement of newly-sequenced genomes. However, apparent biases could be observed in the list of genomes, with samples from the USA over-represented in comparison with other geographic areas where A. flavus is a common agricultural problem. Another apparent bias in the current dataset was the presence of large clusters of the most used biocontrol agents, Afla-Guard and AF36. The presence of these highly similar genomes can bias population genomic algorithms, as shown here (Figs 1, 3). The African continent was clearly under-represented in the dataset and several described variants of the two species lack sequenced type or neotype strains. Samples from natural ecosystems were almost completely missing from the global dataset, with only two samples worldwide. Sampling efforts and sequencing novel 18 IMA Fungus 16: e172343 (2025), DOI: 10.3897/imafungus.16.172343 Walter P. Pfliegler et al.: Aspergillus flavus and oryzae global phylogenomics isolates is paramount to understand the global diversity of these species. This was perhaps best illustrated by the fact that, in this study, while only three Central European samples were sequenced, they represented three distinct major clades: Flavus B, Flavus A and Flavus D. More sequences would certainly reveal more about the biogeography of the clades. It is especially important to note that the entire lineage II (that likely corresponds to the cluster named SBG, considered to be an interspecific mosaic) is missing from our present analyses as none of its members has been sequenced so far. There are also no available sequences for A. flavus var. parvisclerotigenus and, therefore, its placement and relatedness to A. flavus and A. oryzae remains unclear. The clades delineated here based on NJ-based network, ML phylogeny and population genomic analysis corresponded to the previously described ones (in the case of A. flavus, originally referred to as populations) (Drott et al. 2020; Hatmaker et al. 2025) and the fact that A. oryzae was analysed together with the A. flavus genomes in a comprehensive phylogeny enabled an unprecedented overview of their clades’ relationships. Reconstruction of the genome’s evolutionary history recapitulated the existence of A. flavus evolutionary lineages 1A (= Flavus S), IB (= Flavus B) and IC (= Flavus A) and identified A. flavus populations that were intermixed with A. oryzae (Flavus C and D). This further supports the important role of genetic exchange and recombination in driving further clade diversification. For example, the two newly-delineated mosaic groups showed genomic ancestry intermediate between well-established clades of both species, such as Flavus D and Orzyzae C-H and Flavus C and Oryzae C-H. Ongoing genetic exchange and recombination may also confound phylogenetic inferences, which assume a strictly bifurcating evolutionary history. We showed that the mosaic genomes may be misplaced in Maximum Likelihood-based phylogenies (Fig. 4) and only the integration of network and population-genomic analyses can reliably reveal their admixed nature. Relatively high levels of introgression were also detected in the more well-defined clades, for example in Oryzae A-B and in the highly diverse Flavus A clade that was characterised by a highly heterogeneous within-clade structure (Figs 1–4). In the case of the latter, FastStructure identified two to three presumed populations with diverse forms of admixture amongst these, without clear subclades (apart from the biocontrol AF36 and its very closely-related isolates). Lineage IC (= Flavus A) appears to be the evolutionary cradle of A. flavus, harbouring more genetic and functional diversity (e.g. aflatoxin chemotypes) than any other species in section Flavi (Carbone et al. 2007; Moore et al. 2009). While introgression from Flavus B and Flavus D into Flavus A was evident (Figs 2, 3), Flavus B itself did not display such patterns. These patterns are reminiscent of the recent findings of asymmetric introgression into Flavus A observed in a substantial number of North American isolates (Molo et al. 2022) and of the low levels of admixture in Flavus B reported by Hatmaker et al. (Hatmaker et al. 2025). Furthermore, a recently published and subsequently withdrawn preprint described phylogenomic and population genomic analyses of more than 550 unpublished A. flavus genomes from China. These genomes remain inaccessible (Xie et al. 2023). Our understanding of the genetic relatedness amongst agricultural, clinical and domesticated strains of A. flavus remains limited. The Flavus D clade, described recently as an A. flavus population of infectious human isolates, is related to Flavus C (Hatmaker et al. 2025) and was found to be much more closely 19 IMA Fungus 16: e172343 (2025), DOI: 10.3897/imafungus.16.172343 Walter P. Pfliegler et al.: Aspergillus flavus and oryzae global phylogenomics related to A. oryzae clades than to Flavus C in our extended dataset that also included genomes of the koji fungi (Figs 1, 2). This raises the question of whether domestication and the colonisation or infection of the human host are related processes in section Flavi, similar to what has been observed in Saccharomyces cerevisiae baker’s yeasts. In this latter, widely used species, both non-admixed and highly admixed clades adapted to fermentation niches proved to be associated with isolates colonising and infecting humans, while wild clades in natural niches not (Peter et al. 2018; Loegler et al. 2025; Rácz et al. 2025). In the current Aspergillus dataset, human isolates represented a substantial proportion (~ 17%) of all genomes and these isolates predominantly grouped with the Flavus D and Oryzae C-H clades, along with the two mosaic groups related to these clades (and to Flavus C in the case of one mosaic group). The remaining human isolates fell almost exclusively within the Flavus A clade, including a clonal cluster of 11 samples from a single patient from the Netherlands (Buil et al. 2021). As expected, the Oryzae A-B and Oryzae C-H clades were highly enriched in isolates from solid state fermentations for traditional food products and alcoholic beverages and are considered domesticated (Watarai et al. 2019; Chacón-Vargas et al. 2021). Several hypotheses about the relationship between domestication and human infections may be proposed, but their evaluation will require further study. For example, domestication and frequent use could merely have promoted increased opportunities for human infections or, alternatively, domestication may have contributed pre-adaptive traits for the infectious groups of these fungi. This needs to be further elucidated by studies with broader sampling, thereby lowering the potential effect of sampling bias. Amongst the currently available clinical isolate genomes, Europe is vastly over-represented, while agricultural isolates are extremely low in numbers from the continent. This may obscure our understanding on the potential origin of clinical isolates. It is also possible that a single earlier clade bifurcated into what are now predominantly clinical and food fermentation clades with contrasting evolutionary histories. Another important observation about the clinical isolates is that they are not only common in the Flavus D clade (Hatmaker et al. 2025), but also in the Flavus D/Oryzae C-H and the Flavus C/Oryzae C-H mosaic groups, where clinical samples are found almost exclusively, even in higher proportions than in Flavus D (Fig. 2). These findings indicate that genetic admixture may be a significant evolutionary force in the origin of the approximately 30 clinical genomes in these groups. Alternatively, the apparent lack of non-clinical strains in these mosaic groups may simply reflect sampling bias and increased sampling may later clarify whether admixture contributed to pathogenicity in humans. The most striking observations made when comparing all available shortread sequenced genomes are related to the status of the koji fungus and its treatment in regard to the species delimitation of A. flavus. We found that the neotype strain of A. flavus, NRRL 1957, was placed in the Oryzae C-H clade closely related to the neotype strain of A. oryzae var. oryzae, namely RIB40. Similarly, Hatmaker et al. (2025) placed this genome in population D, which, in the present study, was divided into distinct clades and mosaic groups due to broader sampling and the inclusion of koji fungus genomes. This unexpected placement has major taxonomic implications, as it is impossible to split the members of the two Aspergillus species into any two monophyletic groups that: (1) would together encompass all samples and (2) contain only one of the neotypes in each. 20 IMA Fungus 16: e172343 (2025), DOI: 10.3897/imafungus.16.172343 Walter P. Pfliegler et al.: Aspergillus flavus and oryzae global phylogenomics In other words, A. oryzae var. oryzae should not just be treated as a junior synonym of A. flavus because the koji fungus originated from within the diversity of A. flavus, but also because the very clade that contains its neotype already holds the neotype that represents A. flavus itself. A similar conclusion was reached in a preprint, based on phylogenomic analysis of 313 A. flavus and A. oryzae genomes during the review process of this current work (Drott et al. 2025). Furthermore, the koji fungus was split into two well-separated clades that did not form a single monophyletic group, suggesting at least two independent domestication events and very limited subsequent admixture from the Oryzae C-H to Oryzae A-B clades (Figs 2, 3). The Oryzae A-B consisted of two groups of apparently clonal genomes, both with very low genetic diversity (Figs 1, 5). The origin of this clade remains unclear as the number of variants it shared with other clades and mosaic groups (Fig. 5, Suppl. material 1) and its population genomic relationships (Fig. 3) placed it well-separated from others. Additionally, the larger clonal group, namely Oryzae A of Watarai et al. (2019) showed no sign of admixture. Importantly, based on two fully-assembled genomes and their structural differences, it has already been shown that there are differences in alpha-amylase gene duplications, in gene presence/absence and in high-impact mutations in biosynthetic genes in the two A. oryzae clades (Chacón-Vargas et al. 2021). These suggested independent domestication events. Gell et al. (2020) constructed genetic maps, based on three A. flavus mapping populations to gain insights into the genomic organisation of A. flavus and how it differs from A. oryzae RIB40. They found an inverted reciprocal translocation between chromosomes 2 and 6 that was present in all six parental strains, which included strains from both lineages IB (= Flavus B) and lineage IC (= Flavus A), when these genomes were compared to A. oryzae RIB40. They hypothesised that the inverted reciprocal translocation event most likely occurred between chromosomes 2 and 6 in the RIB40 strain, potentially during the process of domestication. It is possible that domestication has led to additional translocation events, occurring independently and multiple times across different species, evolutionary lineages and populations; these are mostly undetectable using short-read sequencing and reference-guided variant discovery methods. Future work should include extensive de novo long-read sequencing and assembly of mosaic and non-mosaic genomes to determine whether structural variants are confounding evolutionary relationships within these groups. These revised phylogenetic relationships have profound implications for several applied fields, including food fermentations and safety, animal nutrition and human clinical diagnostics, as well as metagenomics and metabarcoding. It also complicates efforts to align taxonomy with applied mycological concepts and practices in food technology and safety. A clearer, global overview of the genomic diversity and population structure of these fungi can help in supporting regulatory decisions and diagnostics, especially as agricultural A. flavus and A. oryzae strains are increasingly used as biocontrol agents against toxigenic strains (Lewis et al. 2019; Molo et al. 2019; Khan et al. 2021) and insect pests (You et al. 2024). Similarly, the koji fungus is widely sought after for at-home fermentations, gaining attention as a source for health-promoting foods and post-biotics (Choi et al. 2024; Seidler et al. 2024) and feedstuff supplements (Uwineza et al. 2024). These uses may facilitate the spread of several clades or mosaic groups into novel ecological niches and geographic regions. 21 IMA Fungus 16: e172343 (2025), DOI: 10.3897/imafungus.16.172343 Walter P. Pfliegler et al.: Aspergillus flavus and oryzae global phylogenomics Conclusions Here, for the first time, we analysed all available A. flavus and A. oryzae shortread genomes together. This joint analysis not only confirmed the origin of the domesticated form (A. oryzae) from A. flavus, but also revealed an unexpectedly complex population structure, with mostly well-separated clades interspersed with highly admixed mosaic groups. The distinction between these fungi as two separate species is complicated by the paraphyly of A. flavus when A. oryzae is excluded and by evidence of independent domestication events that gave rise to the koji fungus. Clarifying their taxonomic status will require thorough and, potentially, controversial decisions, but doing so would ultimately improve safety assessments in food production, provide information for biocontrol applications and increase the accuracy of automated metagenomic analyses. Our phylogenomic and population genomic results substantially altered the context in which the origin and diversity of A. flavus clinical isolates are understood, relative to recent work (Hatmaker et al. 2025). Human-infecting isolates are not merely enriched in a single distinct clade of A. flavus (Flavus D), but instead occur across a broader set of clades and mosaic groups that also include traditional food fermentation isolates of A. oryzae. Future studies should, therefore, consider domestication history when analysing genomic markers of pathogenicity. As domestication and admixture can move accessory genes between lineages, the pangenome and a comprehensive analysis of structural variations in the genome is the logical next frontier. Although we did not analyse the pangenome of the 639 samples here, subsequent work should clarify which accessory genes in A. flavus mosaic groups may derive from A. oryzae and, thus, from domestication. These efforts could even adopt a super-pangenome framework, recently demonstrated in plant genomics where domesticated lineages and related wild species are subjected to an integrative comparative analysis (Khan et al. 2020). Such a framework will require expanded geographic sampling (notably, Africa), sequencing of previously described, but unsampled groups (e.g. lineage II) and inclusion of closely-related wild and domesticated species, such as A. minisclerotigenes, A. parasiticus and A. sojae. This has the potential to resolve the evolutionary origins of accessory genes important in virulence, secondary metabolite production or carbohydrate metabolism in domesticated and non-domesticated lineages and would clarify the roles of domestication and admixture played in shaping adaptation to specific niches, like the human host or food fermentation. Acknowledgements We thank anonymous reviewers and the editor Dr. Janneke Aylward for valuable comments and suggestions. Additional information Conflict of interest The authors have declared that no competing interests exist. Ethical statement No ethical statement was reported. 22 IMA Fungus 16: e172343 (2025), DOI: 10.3897/imafungus.16.172343 Walter P. Pfliegler et al.: Aspergillus flavus and oryzae global phylogenomics Use of AI No use of AI was reported. Adherence to national and international regulations All the fungal strains used in this study have been legally obtained, respecting the Convention on Biological Diversity (Rio Convention). Funding This work was supported by Project no. 2018-1.2.1-NKP-2018-00002 which has been implemented with the support provided from the National Research, Development and Innovation Fund of Hungary, financed under the 2018-1.2.1-NKP funding scheme, by Project no. TKP2021-EGA-20 (Biotechnology), which has been implemented with the support provided from the National Research, Development and Innovation Fund of Hungary, financed under the TKP2021-EGA funding scheme, by the National Research, Development and Innovation Office of Hungary with the grant K142801. The study was also funded by the European Union (project 101181208 - MYMATCH). Views and opinions expressed are, however, those of the author(s) only and do not necessarily reflect those of the European Union or CINEA. Neither the European Union nor the granting authority can be held responsible for them. Author contributions Conceptualisation, W.P.P., I.P.; methodology, W.P.P., B.N., V.B.; formal analysis, W.P.P., B.N., V.B.; resources, T.P., I.P.; data curation, W.P.P., V.B.; writing—original draft preparation, W.P.P.; writing—review and editing, W.P.P., I.P, T.P., N.B., I.C.; visualisation, N.B., W.P.P.; supervision, W.P.P., I.P.; project administration, I.P.; funding acquisition, I.P. Author ORCIDs Walter P. Pfliegler https://orcid.org/0000-0001-6723-4416 Bálint Németh https://orcid.org/0009-0004-4654-6073 Veronika Bodnár https://orcid.org/0009-0008-3241-8889 Tünde Pusztahelyi https://orcid.org/0000-0002-5495-6273 Ignazio Carbone https://orcid.org/0000-0003-2721-6656 István Pócsi https://orcid.org/0000-0003-2692-6453 Data availability The datasets generated during the current study are available in the FigShare repository (https://figshare.com/projects/The_clades_mosaic_groups_domestication_history_and_the_global_diversity_of_the_Aspergillus_flavus_oryzae_complex/260192) and in NCBI SRA (PRJNA1307175). References Ajmal M, Alshannaq AF, Moon H et al. (2022) Characterization of 260 isolates of Aspergillus section Flavi obtained from sesame seeds in Punjab, Pakistan. Toxins (Basel) 14: 117. https://doi.org/10.3390/TOXINS14020117 Alshannaq AF, Gibbons JG, Lee MK et al. (2018) Controlling aflatoxin contamination and propagation of Aspergillus flavus by a soy-fermenting Aspergillus oryzae strain. Scientific Reports 8: 1–14. https://doi.org/10.1038/s41598-018-35246-1 23 IMA Fungus 16: e172343 (2025), DOI: 10.3897/imafungus.16.172343 Walter P. Pfliegler et al.: Aspergillus flavus and oryzae global phylogenomics Alvarez F, Arena M, Auteri D et al. (2022) Peer review of the pesticide risk assessment of the active substance Aspergillus flavus strain MUCL54911. EFSA journal. European Food Safety Authority 20: e07202. https://doi.org/10.2903/J.EFSA.2022.7202 Arias RS, Orner VA, Martinez-Castillo J et al. (2021) Aspergillus section Flavi, need for a robust taxonomy. Microbiology Resource Announcements 10: 99–101. https://doi. org/10.1128/mra.00784-21 Arias RS, Mohammed A, Orner VA et al. (2020) Sixteen draft genome sequences representing the genetic diversity of Aspergillus flavus and Aspergillus parasiticus colonizing peanut seeds in Ethiopia. Microbiology Resource Announcements 9: e00591-20. https://doi.org/10.1128/MRA.00591-20 Augusto J, Atehnkeng J, Ortega-Beltran A et al. (2024) Keeping toxigenic Aspergillus section Flavi and aflatoxin contamination at bay by deploying atoxigenic-based biocontrol products during production of groundnut and maize in Mozambique. Frontiers in Microbiology 15: 1501924. https://doi.org/10.3389/FMICB.2024.1501924 Van der Auwera GA, Carneiro MO, Hartl C et al. (2013) From fastQ data to high-confidence variant calls: The genome analysis toolkit best practices pipeline. Current Protocols in Bioinformatics 43: 11.10.1-11.10.33. https://doi.org/10.1002/0471250953. bi1110s43 Bal J, Yun SH, Chun J et al. (2016) Taxonomic characterization, evaluation of toxigenicity, and saccharification capability of Aspergillus section Flavi isolates from Korean traditional wheat-based fermentation starter nuruk. Mycobiology 44: 155–161. https://doi.org/10.5941/MYCO.2016.44.3.155 Basenko EY, Shanmugasundram A, Böhme U et al. (2024) What is new in FungiDB: a webbased bioinformatics platform for omics-scale data analysis for fungal and oomycete species. Genetics 227: iyae035. https://doi.org/10.1093/GENETICS/IYAE035 Buil JB, Houbraken J, Reijers MH et al. (2021) Genetic and phenotypic characterization of in-host developed azole-resistant Aspergillus flavus isolates. Journal of Fungi (Basel) 7: 164. https://doi.org/10.3390/JOF7030164 Carbone I, Ramirez-Prado JH, Jakobek JL et al. (2007) Gene duplication, modularity and adaptation in the evolution of the aflatoxin gene cluster. BMC Evolutionary Biology 7: 111. https://doi.org/10.1186/1471-2148-7-111 Chacón-Vargas K, McCarthy CO, Choi D et al. (2021) Comparison of two Aspergillus oryzae genomes from different clades reveals independent evolution of alpha-amylase duplication, variation in secondary metabolism genes, and differences in primary metabolism. Frontiers in Microbiology 12: 691296. https://doi.org/10.3389/ FMICB.2021.691296 Chang CC, Chow CC, Tellier LCAM et al. (2015) Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience 4: 7. https://doi.org/10.1186/ S13742-015-0047-8 Chang PK, Ehrlich KC, Hua SST (2006) Cladal relatedness among Aspergillus oryzae isolates and Aspergillus flavus S and L morphotype isolates. International Journal of Food Microbiology 108: 172–177. https://doi.org/10.1016/J.IJFOODMICRO.2005.11.008 Chen S, Zhou Y, Chen Y et al. (2018) fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34: i884–i890. https://doi.org/10.1093/bioinformatics/bty560 Chin YW, Hong SP, Lim SD et al. (2024) Investigation of microbial community of Korean soy sauce (ganjang) using shotgun metagenomic sequencing and its relationship with sensory characteristics. Microorganisms 12: 2559. https://doi.org/10.3390/microorganisms12122559 24 IMA Fungus 16: e172343 (2025), DOI: 10.3897/imafungus.16.172343 Walter P. Pfliegler et al.: Aspergillus flavus and oryzae global phylogenomics Choi D, Alshannaq AF, Bok Y et al. (2024) Broad-spectrum antimicrobial activities of a food fermentate of Aspergillus oryzae. Microbiology Spectrum 12: e0185424. https:// doi.org/10.1128/spectrum.01854-24 Cohn F (1884) Über Schimmelpilze aIs Gährungserreger. Jahresbericht der Schlesischen Gesellschaft für Vaterländische Kultur 61: 226–229. Costantin M, Lucet M (1905) Quelques Aspergillus pathogènes. Annales des Sciences Naturelles. Botanique 9: 119–171. Dabas Y, Xess I, Pandey M et al. (2022) Epidemiology and antifungal susceptibility patterns of invasive fungal infections (IFIs) in India: A prospective observational study. Journal of Fungi (Basel) 8: 33. https://doi.org/10.3390/JOF8010033 Dorner JW (2004) Biological control of aflatoxin contamination of crops. Journal of Toxicology: Toxin Reviews 23: 425–450. https://doi.org/10.1081/TXR-200027877 Drott MT, Rush TA, Satterlee TR et al. (2021) Microevolution in the pansecondary metabolome of Aspergillus flavus and its potential macroevolutionary implications for filamentous fungi. Proceedings of the National Academy of Sciences of the United States of America 118: e2021683118. https://doi.org/10.1073/pnas.202168311 Drott MT, Wang Y-W, Hatmaker EA et al. (2025) Population-specific transcriptomic shifts underlie secondary metabolic diversification in Aspergillus flavus and the domestication of Aspergillus oryzae. bioRxiv 2025.10.02.680074. https://doi. org/10.1101/2025.10.02.680074 Drott MT, Satterlee TR, Skerker JM et al. (2020) The frequency of sex: Population genomics reveals differences in recombination and population structure of the aflatoxin-producing fungus Aspergillus flavus. mBio 11: 1–13. https://doi.org/10.1128/ MBIO.00963-20 Fountain JC, Clevenger JP, Nadon B et al. (2020) Draft genome sequences of one Aspergillus parasiticus isolate and nine Aspergillus flavus isolates with varying stress tolerance and aflatoxin production. Microbiology Resource Announcements 9: e0047820. https://doi.org/10.1128/MRA.00478-20 Frisvad JC, Hubka V, Ezekiel CN et al. (2019) Taxonomy of Aspergillus section Flavi and their production of aflatoxins, ochratoxins and other mycotoxins. Studies in Mycology 93: 1–63. https://doi.org/10.1016/j.simyco.2018.06.001 Galili T, O’Callaghan A, Sidi J et al. (2018) heatmaply: an R package for creating interactive cluster heatmaps for online publishing. Bioinformatics 34: 1600–1602. https:// doi.org/10.1093/BIOINFORMATICS/BTX657 Gangurde SS, Korani W, Bajaj P et al. (2024) Aspergillus flavus pangenome (AflaPan) uncovers novel aflatoxin and secondary metabolite associated gene clusters. BMC Plant Biology 24: 1–17. https://doi.org/10.1186/s12870-024-04950-8 Geiser DM, Pitt JI, Taylor JW (1998) Cryptic speciation and recombination in the aflatoxin-producing fungus Aspergillus flavus. Proceedings of the National Academy of Sciences of the United States of America 95: 388–393. https://doi.org/10.1073/PNAS.95.1.388 Geiser DM, Dorner JW, Horn BW et al. (2000) The phylogenetics of mycotoxin and sclerotium production in Aspergillus flavus and Aspergillus oryzae. Fungal Genetics and Biology 31: 169–179. https://doi.org/10.1006/FGBI.2000.1215 Gell RM, Horn BW, Carbone I (2020) Genetic map and heritability of Aspergillus flavus. Fungal Genetics and Biology 144: 103478. https://doi.org/10.1016/J.FGB.2020.103478 Gibbons JG, Salichos L, Slot JC et al. (2012) The evolutionary imprint of domestication on genome variation and function of the filamentous fungus Aspergillus oryzae. Current Biology 22: 1403–1409. https://doi.org/10.1016/J.CUB.2012.05.033 25 IMA Fungus 16: e172343 (2025), DOI: 10.3897/imafungus.16.172343 Walter P. Pfliegler et al.: Aspergillus flavus and oryzae global phylogenomics Glöckner FO, Yilmaz P, Quast C et al. (2017) 25 years of serving the community with ribosomal RNA gene reference databases and tools. Journal of Biotechnology 261: 169–176. https://doi.org/10.1016/J.JBIOTEC.2017.06.1198 Han DM, Baek JH, Choi DG et al. (2024) Comparative pangenome analysis of Aspergillus flavus and Aspergillus oryzae reveals their phylogenetic, genomic, and metabolic homogeneity. Food Microbiology 119: 104435. https://doi.org/10.1016/J. FM.2023.104435 Hatmaker EA, Barber AE, Drott MT et al. (2025) Population structure in a fungal human pathogen is potentially linked to pathogenicity. Nature Communications 16: 1–12. https://doi.org/10.1038/s41467-025-62777-9 Hatmaker EA, Barber AE, Drott MT et al. (2024) Pathogenicity is associated with population structure in a fungal pathogen of humans. bioRxiv 2024.07.05.602241 https:// doi.org/10.1101/2024.07.05.602241 Hill JH, Round JL (2024) Intestinal fungal-host interactions in promoting and maintaining health. Cell Host & Microbe 32: 1668–1680. https://doi.org/10.1016/J. CHOM.2024.09.010 Horn BW (2003) Ecology and population biology of aflatoxigenic fungi in soil. Journal of Toxicology: Toxin Reviews 22: 351–379. https://doi.org/10.1081/TXR-120024098 Houbraken J, Kocsubé S, Visagie CM et al. (2020) Classification of Aspergillus, Penicillium, Talaromyces and related genera (Eurotiales): An overview of families, genera, subgenera, sections, series and species. Studies in Mycology 95: 5–169. https://doi. org/10.1016/J.SIMYCO.2020.05.002 Houbraken J, de Vries RP, Samson RA (2014) Modern taxonomy of biotechnologically important Aspergillus and Penicillium species. Advances in Applied Microbiology 86: 199–249. https://doi.org/10.1016/B978-0-12-800262-9.00004-4 Houbraken J, Visagie CM, Frisvad JC (2021) Recommendations to prevent taxonomic misidentification of genome-sequenced fungal strains. Microbiology Resource Announcements 10: e0107420. https://doi.org/10.1128/MRA.01074-20 Huson DH, Bryant D (2006) Application of phylogenetic networks in evolutionary studies. Molecular Biology and Evolution 23: 254–267. https://doi.org/10.1093/MOLBEV/MSJ030 Hyde KD, Al-Hatmi AMS, Andersen B et al. (2018) The world’s ten most feared fungi. Fungal Diversity 93: 161–194. https://doi.org/10.1007/s13225-018-0413-9 Jeong E, Seo JA (2025) Comparative genome analyses of Aspergillus oryzae and Aspergillus flavus originated from a Korean fermentation starter, nuruk. Food Microbiology 131. https://doi.org/10.1016/j.fm.2025.104807 Kalyaanamoorthy S, Minh BQ, Wong TKF et al. (2017) ModelFinder: Fast model selection for accurate phylogenetic estimates. Nature Methods 14: 587. https://doi. org/10.1038/NMETH.4285 Kanaujia R, Singh S, Rudramurthy SM (2023) Aspergillosis: an update on clinical spectrum, diagnostic schemes, and management. Current Fungal Infection Reports 17: 144–155. https://doi.org/10.1007/S12281-023-00461-5 Khan AW, Garg V, Roorkiwal M et al. (2020) Super-pangenome by integrating the wild side of a species for accelerated crop improvement. Trends in Plant Science 25: 148–158. https://doi.org/10.1016/j.tplants.2019.10.012 Khan R, Ghazali FM, Mahyudin NA et al. (2021) Biocontrol of aflatoxins using non-aflatoxigenic Aspergillus flavus: A literature review. Journal of Fungi (Basel) 7: 381. https://doi.org/10.3390/JOF7050381