Chromosome assembled and annotated genome sequence of Aspergillus flavus NRRL 3357

zol & fai: large-scale targeted detection and evolutionary investigation of gene clusters

Many universally and conditionally important genes are genomically aggregated within clusters. Here, we introduce fai and zol, which together enable large-scale comparative analysis of different types of gene clusters and mobile-genetic elements (MGEs), such as biosynthetic gene clusters (BGCs) or viruses. Fundamentally, they overcome a current bottleneck to reliably perform comprehensive orthology inference at large scale across broad taxonomic contexts and thousands of genomes. First, fai allows the identification of orthologous instances of a query gene cluster of interest amongst a database of target genomes. Subsequently, zol enables reliable, context-specific inference of ortholog groups for individual protein-encoding genes across gene cluster instances. In addition, zol performs functional annotation and computes a variety of evolutionary statistics for each inferred ortholog group. Importantly, in comparison to tools for visual exploration of homologous relationships between gene clusters, zol can scale to thousands of gene cluster instances and produce detailed reports that are easy to digest. To showcase fai and zol, we apply them for: (i) longitudinal tracking of a virus in metagenomes, (ii) discovering novel population-level genetic insights of two common BGCs in the fungal species Aspergillus flavus, and (iii) uncovering large-scale evolutionary trends of a virulence-associated gene cluster across thousands of genomes from a diverse bacterial genus.

Structure-guided discovery of novel AflG inhibitors for aflatoxin contamination control in aspergillus flavus

Aflatoxins (AFs) are highly carcinogenic metabolites produced by Aspergillus species that can contaminate critical food staples, leading to significant health and economic risks. The cytochrome P450 monooxygenase AflG catalyzes an early step in AF biosynthesis, resulting in the conversion of averantin (AVN) to 5'-hydroxy-averantin. However, the molecular mechanism underlying the AflG-AVN interaction remains unclear. Here, we sought to understand the structural features of AflG in complex with AVN to enable the identification of inhibitors targeting the AflG binding pocket. To achieve this goal, we employed a comprehensive approach combining computational and experimental methods. Structural modeling and microsecond-scale molecular dynamics (MD) simulations yielded new insights into AflG architecture and unveiled unique ligand binding conformations of the AflG-AVN complex. High-throughput virtual screening of more than 1.3 million compounds pinpointed specific subsets with favorable predicted docking scores. The resulting compounds were ranked based on binding free energy calculations and evaluated with MD simulations and in vitro experiments with Aspergillus flavus. Our results revealed two compounds significantly inhibited AF biosynthesis. Comprehensive structural analysis elucidated the binding sites of competitive inhibitors and demonstrated their regulation of AflG dynamics. This structure-guided pipeline successfully enabled the identification of novel AflG inhibitors and provided novel molecular insights that will guide future efforts to develop effective therapeutics that prevent AF contamination.

Pathogenicity is associated with population structure in a fungal pathogen of humans

Aspergillus flavus is a clinically and agriculturally important saprotrophic fungus responsible for severe human infections and extensive crop losses. We analyzed genomic data from 250 (95 clinical and 155 environmental) A. flavus isolates from 9 countries, including 70 newly sequenced clinical isolates, to examine population and pan-genome structure and their relationship to pathogenicity. We identified five A. flavus populations, including a new population, D, corresponding to distinct clades in the genome-wide phylogeny. Strikingly, > 75% of clinical isolates were from population D. Accessory genes, including genes within biosynthetic gene clusters, were significantly more common in some populations but rare in others. Population D was enriched for genes associated with zinc ion binding, lipid metabolism, and certain types of hydrolase activity. In contrast to the major human pathogen Aspergillus fumigatus, A. flavus pathogenicity in humans is strongly associated with population structure, making it a great system for investigating how population-specific genes contribute to pathogenicity.

Aspergillus flavus pangenome (AflaPan) uncovers novel aflatoxin and secondary metabolite associated gene clusters

BACKGROUND

Aspergillus flavus is an important agricultural and food safety threat due to its production of carcinogenic aflatoxins. It has high level of genetic diversity that is adapted to various environments. Recently, we reported two reference genomes of A. flavus isolates, AF13 (MAT1-2 and highly aflatoxigenic isolate) and NRRL3357 (MAT1-1 and moderate aflatoxin producer). Where, an insertion of 310 kb in AF13 included an aflatoxin producing gene bZIP transcription factor, named atfC. Observations of significant genomic variants between these isolates of contrasting phenotypes prompted an investigation into variation among other agricultural isolates of A. flavus with the goal of discovering novel genes potentially associated with aflatoxin production regulation. Present study was designed with three main objectives: (1) collection of large number of A. flavus isolates from diverse sources including maize plants and field soils; (2) whole genome sequencing of collected isolates and development of a pangenome; and (3) pangenome-wide association study (Pan-GWAS) to identify novel secondary metabolite cluster genes.

RESULTS

Pangenome analysis of 346 A. flavus isolates identified a total of 17,855 unique orthologous gene clusters, with mere 41% (7,315) core genes and 59% (10,540) accessory genes indicating accumulation of high genomic diversity during domestication. 5,994 orthologous gene clusters in accessory genome not annotated in either the A. flavus AF13 or NRRL3357 reference genomes. Pan-genome wide association analysis of the genomic variations identified 391 significant associated pan-genes associated with aflatoxin production. Interestingly, most of the significantly associated pan-genes (94%; 369 associations) belonged to accessory genome indicating that genome expansion has resulted in the incorporation of new genes associated with aflatoxin and other secondary metabolites.

CONCLUSION

In summary, this study provides complete pangenome framework for the species of Aspergillus flavus along with associated genes for pathogen survival and aflatoxin production. The large accessory genome indicated large genome diversity in the species A. flavus, however AflaPan is a closed pangenome represents optimum diversity of species A. flavus. Most importantly, the newly identified aflatoxin producing gene clusters will be a new source for seeking aflatoxin mitigation strategies and needs new attention in research.

Genetic regulation of L-tryptophan metabolism in Psilocybe mexicana supports psilocybin biosynthesis

BACKGROUND

Although Basidiomycota produce pharmaceutically and ecologically relevant natural products, knowledge of how they coordinate their primary and secondary metabolism is virtually non-existent. Upon transition from vegetative mycelium to carpophore formation, mushrooms of the genus Psilocybe use L-tryptophan to supply the biosynthesis of the psychedelic tryptamine alkaloid psilocybin with the scaffold, leading to a strongly increased demand for this particular amino acid as this alkaloid may account for up to 2% of the dry mass. Using Psilocybe mexicana as our model and relying on genetic, transcriptomic, and biochemical methods, this study investigated if L-tryptophan biosynthesis and degradation in P. mexicana correlate with natural product formation.

RESULTS

A comparative transcriptomic approach of gene expression in P. mexicana psilocybin non-producing vegetative mycelium versus producing carpophores identified the upregulation of L-tryptophan biosynthesis genes. The shikimate pathway genes trpE1, trpD, and trpB (encoding anthranilate synthase, anthranilate phosphoribosyltransferase, and L-tryptophan synthase, respectively) were upregulated in carpophores. In contrast, genes idoA and iasA, encoding indole-2,3-dioxygenase and indole-3-acetaldehyde synthase, i.e., gateway enzymes for L-tryptophan-consuming pathways, were massively downregulated. Subsequently, IasA was heterologously produced in Escherichia coli and biochemically characterized in vitro. This enzyme represents the first characterized microbial L-tryptophan-preferring acetaldehyde synthase. A comparison of transcriptomic data collected in this study with prior data of Psilocybe cubensis showed species-specific differences in how L-tryptophan metabolism genes are regulated, despite the close taxonomic relationship.

CONCLUSIONS

The upregulated L-tryptophan biosynthesis genes and, oppositely, the concomitant downregulated genes encoding L-tryptophan-consuming enzymes reflect a well-adjusted cellular system to route this amino acid toward psilocybin production. Our study has pilot character beyond the genus Psilocybe and provides, for the first time, insight in the coordination of mushroom primary and secondary metabolism.

Role of RNA Modifications, Especially m6A, in Aflatoxin Biosynthesis of Aspergillus flavus

RNA modifications play key roles in eukaryotes, but the functions in Aspergillus flavus are still unknown. Temperature has been reported previously to be a critical environmental factor that regulates the aflatoxin production of A. flavus, but much remains to be learned about the molecular networks. Here, we demonstrated that 12 kinds of RNA modifications in A. flavus were significantly changed under 29 °C compared to 37 °C incubation; among them, m6A was further verified by a colorimetric method. Then, the transcriptome-wide m6A methylome and m6A-altered genes were comprehensively illuminated through methylated RNA immunoprecipitation sequencing and RNA sequencing, from which 22 differentially methylated and expressed transcripts under 29 °C were screened out. It is especially notable that AFCA_009549, an aflatoxin biosynthetic pathway gene (aflQ), and the m6A methylation of its 332nd adenine in the mRNA significantly affect aflatoxin biosynthesis in A. flavus both on media and crop kernels. The content of sterigmatocystin in both ΔaflQ and aflQA332C strains was significantly higher than that in the WT strain. Together, these findings reveal that RNA modifications are associated with secondary metabolite biosynthesis of A. flavus.

The KdmB-EcoA-RpdA-SntB (KERS) chromatin regulatory complex controls development, secondary metabolism and pathogenicity in Aspergillus flavus

The filamentous fungus Aspergillus flavus is a plant and human pathogen predominantly found in the soil as spores or sclerotia and is capable of producing various secondary metabolites (SM) such as the carcinogenic mycotoxin aflatoxin. Recently, we have discovered a novel nuclear chromatin binding complex (KERS) that contains the JARID1-type histone demethylase KdmB, a putative cohesion acetyl transferase EcoA, a class I type histone deacetylase RpdA and the PHD ring finger reader protein SntB in the model filamentous fungus Aspergillus nidulans. Here, we show the presence of the KERS complex in A. flavus by immunoprecipitation-coupled mass spectrometry and constructed kdmBΔ and rpdAΔ strains to study their roles in fungal development, SM production and histone post-translational modifications (HPTMs). We found that KdmB and RpdA couple the regulation of SM gene clusters with fungal light-responses and HPTMs. KdmB and RpdA have opposing roles in light-induced asexual conidiation, while both factors are positive regulators of sclerotia development through the nsdC and nsdD pathway. KdmB and RpdA are essential for the productions of aflatoxin (similar to findings for SntB) as well as cyclopiazonic acid, ditryptophenaline and leporin B through controlling the respective SM biosynthetic gene clusters. We further show that both KdmB and RpdA regulate H3K4me3 and H3K9me3 levels, while RpdA also acts on H3K14ac levels in nuclear extracts. Therefore, the chromatin modifiers KdmB and RpdA of the KERS complex are key regulators for fungal development and SM metabolism in A. flavus.

Identification of Virulence Factors in Entomopathogenic Aspergillus flavus Isolated from Naturally Infected Rhipicephalus microplus

Aspergillus flavus has been found to be an effective entomopathogenic fungus for various arthropods, including ticks. In particular, natural fungal infections in cattle ticks show promise for biocontrol of the Rhipicephalus (Boophilus) microplus tick, which is a major ectoparasite affecting cattle worldwide. Our study aimed to elucidate the specific entomopathogenic virulence factors encoded in the genome of an A. flavus strain isolated from naturally infected cattle ticks. We performed morphological and biochemical phenotyping alongside complete genome sequencing, which revealed that the isolated fungus was A. flavus related to the L morphotype, capable of producing a range of gene-coded entomopathogenic virulence factors, including ribotoxin, aflatoxin, kojic acid, chitinases, killer toxin, and satratoxin. To evaluate the efficacy of this A. flavus strain against ticks, we conducted experimental bioassays using healthy engorged female ticks. A morbidity rate of 90% was observed, starting at a concentration of 105 conidia/mL. At a concentration of 107 conidia/mL, we observed a 50% mortality rate and a 21.5% inhibition of oviposition. The highest levels of hatch inhibition (30.8%) and estimated reproduction inhibition (34.64%) were achieved at a concentration of 108 conidia/mL. Furthermore, the tick larval progeny that hatched from the infected tick egg masses showed evident symptoms of Aspergillus infection after incubation.

Bifurcate evolution of quinone synthetases in basidiomycetes

BACKGROUND

The terphenylquinones represent an ecologically remarkable class of basidiomycete natural products as they serve as central precursors of pigments and compounds that impact on microbial consortia by modulating bacterial biofilms and motility. This study addressed the phylogenetic origin of the quinone synthetases that assemble the key terphenylquinones polyporic acid and atromentin.

RESULTS

The activity of the Hapalopilus rutilans synthetases HapA1, HapA2 and of Psilocybe cubensis PpaA1 were reconstituted in Aspergilli. Liquid chromatography and mass spectrometry of the culture extracts identified all three enzymes as polyporic acid synthetases. PpaA1 is unique in that it features a C-terminal, yet catalytically inactive dioxygenase domain. Combined with bioinformatics to reconstruct the phylogeny, our results demonstrate that basidiomycete polyporic acid and atromentin synthetases evolved independently, although they share an identical catalytic mechanism and release structurally very closely related products. A targeted amino acid replacement in the substrate binding pocket of the adenylation domains resulted in bifunctional synthetases producing both polyporic acid and atromentin.

CONCLUSIONS

Our results imply that quinone synthetases evolved twice independently in basidiomycetes, depending on the aromatic α-keto acid substrate. Furthermore, key amino acid residues for substrate specificity were identified and changed which led to a relaxed substrate profile. Therefore, our work lays the foundation for future targeted enzyme engineering.

Tractable and Expressive Generative Models of Genetic Variation Data

Population genetic studies often rely on artificial genomes (AGs) simulated by generative models of genetic data. In recent years, unsupervised learning models, based on hidden Markov models, deep generative adversarial networks, restricted Boltzmann machines, and variational autoencoders, have gained popularity due to their ability to generate AGs closely resembling empirical data. These models, however, present a tradeoff between expressivity and tractability. Here, we propose to use hidden Chow-Liu trees (HCLTs) and their representation as probabilistic circuits (PCs) as a solution to this tradeoff. We first learn an HCLT structure that captures the long-range dependencies among SNPs in the training data set. We then convert the HCLT to its equivalent PC as a means of supporting tractable and efficient probabilistic inference. The parameters in these PCs are inferred with an expectation-maximization algorithm using the training data. Compared to other models for generating AGs, HCLT obtains the largest log-likelihood on test genomes across SNPs chosen across the genome and from a contiguous genomic region. Moreover, the AGs generated by HCLT more accurately resemble the source data set in their patterns of allele frequencies, linkage disequilibrium, pairwise haplotype distances, and population structure. This work not only presents a new and robust AG simulator but also manifests the potential of PCs in population genetics.

Whole-genome and comparative genome analysis of Mucor racemosus C isolated from Yongchuan Douchi

Mucor racemosus is the predominant fungal in the zhiqu stage of the fermentation of Yongchuan Douchi (Mucor-type), which plays an important role in the fermentation process of Yongchuan Douchi. However, there is a lack of information on the genetic analysis of M. racemosus. In this study, we isolated and identified M. racemosus C (accession no JAPEHQ000000000) from Yongchuan Douchi and analyzed the physiological indicators, then genomic information of the strain to perform a comprehensive analysis of its fermentation capacity and safety. M. racemosus C had neutral protease activity up to 68.051 U/mL at 30 °C and alkaline protease activity up to 57.367 U/mL at 25 °C. In addition, comparing the genomic data with the COGs database (NCBI), it was predicted that M. racemosus C undergoes extensive amino acid metabolism, making C suitable for the production of fermented foods (e.g., Douchi, Syoyu, and sufu). Finally, we performed virulence genes and resistance genes analysis, hemolysis experiment, aflatoxins assay, antibiotic resistance assay to evaluate the safety of M. racemosus C, and the results showed that M. racemosus C was safe, non-toxin-producing and non-hemolytic.

Hybrid de novo genome assembly and comparative genomics of three different isolates of Gnomoniopsis castaneae

The first genome assemblies of Gnomoniopsis castaneae (syn. G. smithogilvyi), the causal agent of chestnut brown rot of kernels, shoot blight and cankers, are provided here. Specifically, the complete genome of the Italian ex-type MUT401 isolate was compared to the draft genome of a second Italian isolate (GN01) and to the ICMP 14040 isolate from New Zealand. The three genome sequences were obtained through a hybrid assembly using both short Illumina reads and long Nanopore reads, their coding sequences were annotated and compared with each other and with other Diaporthales. The information offered by the genome assembly of the three isolates represents the base of data for further application related to -omics strategies of the fungus and to develop markers for population studies at a local and global scale.

A cell-type-specific multi-protein complex regulates expression of Cyclin B protein in Drosophila male meiotic prophase

During meiosis, germ cell and stage-specific components impose additional layers of regulation on the core cell cycle machinery to set up an extended G2 period termed meiotic prophase. In Drosophila males, meiotic prophase lasts 3.5 days, during which spermatocytes turn up expression of over 3000 genes and grow 25-fold in volume. Previous work showed that the core cell cycle regulator Cyclin B (CycB) is subject to translational repression in immature Drosophila spermatocytes, mediated by the RNA-binding protein Rbp4 and its partner Fest. Here we show that another spermatocyte-specific protein, Lut, is required for translational repression of cycB in an 8-hour window just before spermatocytes are fully mature. In males mutant for rbp4 or lut , spermatocytes enter and exit the meiotic divisions 6-8 hours earlier than in wild-type. In addition, we show that spermatocyte-specific isoforms of Syncrip (Syp) are required for expression of CycB protein and normal entry into the meiotic divisions. Both Lut and Syp interact with Fest in an RNA-independent manner. Thus a complex of spermatocyte-specific regulators choreograph the timing of expression of CycB protein during male meiotic prophase.

SUMMARY STATEMENT

Expression of a conserved cell cycle component, Cyclin B, is regulated by multiple mechanisms in the Drosophila male germline to dictate the correct timing of meiotic division.

Investigating the origin of subtelomeric and centromeric AT-rich elements in Aspergillus flavus

An in silico study of Aspergillus flavus genome stability uncovered significant variations in both coding and non-coding regions. The non-coding insertions uniformly consisted of AT-rich sequences that are evolutionarily maintained, albeit distributed at widely different sites in an array of A. flavus strains. A survey of ≥ 2kb AT-rich elements (AT ≥ 70%; ATEs) in non-centromeric regions uncovered two major categories of ATEs. The first category is composed of homologous insertions at ectopic, non-allelic sites that contain homology to transposable elements (TEs; Classes B, C, D, and E). Strains differed significantly in frequency, position, and TE type, but displayed a common enrichment in subtelomeric regions. The TEs were heavily mutated, with patterns consistent with the ancestral activity of repeat-induced point mutations (RIP). The second category consists of a conserved set of novel subtelomeric ATE repeats (Classes A, G, G, H, I and J) which lack discernible TEs and, unlike TEs, display a constant polarity relative to the telomere. Members of one of these classes are derivatives of a progenitor ATE that is predicted to have undergone extensive homologous recombination during evolution. A third category of ATEs consists of ~100 kb regions at each centromere. Centromeric ATEs and TE clusters within these centromeres display a high level of sequence identity between strains. These studies suggest that transposition and RIP are forces in the evolution of subtelomeric and centromeric structure and function.

Phosphomannose Isomerase Is Involved in Development, Stress Responses, and Pathogenicity of Aspergillus flavus

Aspergillus flavus causes invasive aspergillosis in immunocompromised patients and severe contamination of agriculturally important crops by producing aflatoxins. The fungal cell wall is absent in animals and is structurally different from that of plants, which makes it a potential antifungal drug target due to its essentiality for fungal survival. Mannose is one of the important components in the fungal cell wall, which requires GDP-mannose (GDP-Man) as the primary donor. Three consecutive enzymes, namely, phosphomannose isomerase (PMI), phosphomannose mutase (PMM), and GDP-mannose phosphorylase (GMPP), are required for GDP-Man biosynthesis. Thus, PMI is of prime importance in cell wall biosynthesis and also has an active role in sugar metabolism. Here, we investigated the functional role of PMI in A. flavus by generating a pmiA-deficient strain. The mutant required exogenous mannose to survive and exhibited reduced growth rate, impaired conidiation, early germination, disturbance in stress responses, and defects in colonization of crop seeds. Furthermore, attenuated virulence of the mutant was documented in both Caenorhabditis elegans and Galleria mellonella infection models. Our results suggested that PMI plays an important role in the development, stress responses, and pathogenicity of A. flavus and therefore could serve as a potential target for battling against infection and controlling aflatoxin contamination caused by A. flavus. IMPORTANCE Aspergillus flavus is a common fungal pathogen of humans, animals, and agriculturally important crops. It causes invasive aspergillosis in humans and also produces highly carcinogenic mycotoxins in postharvest crops that threaten food safety worldwide. To alleviate or eliminate the threats posed by A. flavus, it is necessary to identify genes involved in pathogenicity and mycotoxin contamination. However, little progress has been made in this regard. Here, we focused on PMI, which is the first enzyme involved in the biosynthesis pathway of GDP-Man and thus is important for cell wall synthesis and protein glycosylation. Our study revealed that PMI is important for growth of A. flavus. It is also involved in conidiation, germination, morphogenesis, stress responses, and pathogenicity of A. flavus. Thus, PMI is a potent antifungal target to curb the threats posed by A. flavus.

Draft Genome Sequence of an Aflatoxin-Producing Aspergillus flavus Strain Isolated from Food

Aspergillus flavus is the main producer of carcinogenic aflatoxins and thus is one of the most important fungal food contaminants. Here, we report that the genome of A. flavus strain MRI19 was sequenced using MiSeq and PacBio platforms and that a hybrid assembly was generated.