Complete genome of the toxic mold Aspergillus pseudotamarii isolate NRRL 25517 reveals genomic instability of the aflatoxin biosynthesis cluster

Abstract

Fungi can synthesize a broad array of secondary metabolite chemicals. The genes underpinning their biosynthesis are typically arranged in tightly linked clusters in the genome. For example, ∼25 genes responsible for the biosynthesis of carcinogenic aflatoxins by Aspergillus section Flavi species are grouped in a ∼70 Kb cluster. Assembly fragmentation prevents assessment of the role of structural genomic variation in secondary metabolite evolution in this clade. More comprehensive analyses of secondary metabolite evolution will be possible by working with more complete and accurate genomes of taxonomically diverse Aspergillus species. Here, we combined short- and long-read DNA sequencing to generate a highly contiguous genome of the aflatoxigenic fungus, Aspergillus pseudotamarii (isolate NRRL 25517 = CBS 766.97; scaffold N50 = 5.5 Mb). The nuclear genome is 39.4 Mb, encompassing 12,639 putative protein- encoding genes and 74-97 candidate secondary metabolite biosynthesis gene clusters. The circular mitogenome is 29.7 Kb and contains 14 protein-encoding genes that are highly conserved across the genus. This highly contiguous A. pseudotamarii genome assembly enables comparisons of genomic rearrangements between Aspergillus section Flavi series Kitamyces and series Flavi. Although the aflatoxin biosynthesis gene cluster of A. pseudotamarii is conserved with Aspergillus flavus, the cluster has an inverted orientation relative to the telomere and occurs on a different chromosome.