Genome sequencing of evolved aspergilli populations reveals robust genomes, transversions in A. flavus, and sexual aberrancy in non-homologous end-joining mutants

Abstract

Background: Aspergillus spp. comprises a very diverse group of lower eukaryotes with a high relevance for industrial applications and clinical implications. These multinucleate species are often cultured for many generations in the laboratory, which can unknowingly propagate hidden genetic mutations. To assess the likelihood of such events, we studied the genome stability of aspergilli by using a combination of mutation accumulation (MA) lines and whole genome sequencing. Results: We sequenced the whole genomes of 30 asexual and 10 sexual MA lines of three Aspergillus species (A. flavus, A. fumigatus and A. nidulans) and estimated that each MA line accumulated mutations for over 4000 mitoses during asexual cycles. We estimated mutation rates of 4.2 × 10−11 (A. flavus), 1.1 × 10−11 (A. fumigatus) and 4.1 × 10−11 (A. nidulans) per site per mitosis, suggesting that the genomes are very robust. Unexpectedly, we found a very high rate of GC→TA transversions only in A. flavus. In parallel, 30 asexual lines of the non-homologous endjoining (NHEJ) mutants of the three species were also allowed to accumulate mutations for the same number of mitoses. Sequencing of these NHEJ MA lines gave an estimated mutation rate of 5.1 × 10−11 (A. flavus), 2.2 × 10−11 (A. fumigatus) and 4.5 × 10−11 (A. nidulans) per base per mitosis, which is slightly higher than in the wild-type strains and some ~ 5–6 times lower than in the yeasts. Additionally, in A. nidulans, we found a NHEJ-dependent interference of the sexual cycle that is independent of the accumulation of mutations. Conclusions: We present for the first time direct counts of the mutation rate of filamentous fungal species and find that Aspergillus genomes are very robust. Deletion of the NHEJ machinery results in a slight increase in the mutation rate, but at a rate we suggest is still safe to use for biotechnology purposes. Unexpectedly, we found GC→TA transversions predominated only in the species A. flavus, which could be generated by the hepatocarcinogen secondary metabolite aflatoxin. Lastly, a strong effect of the NHEJ mutation in self-crossing was observed and an increase in the mutations of the asexual lines was quantified