Bisulfite sequencing reveals that Aspergillus flavus holds a hollow in DNA methylation

DNA methylation plays important roles in lifestyle transition of Arthrobotrys oligospora

Trap formation is the key indicator of carnivorous lifestyle transition of nematode-trapping fungi (NTF). Here, the DNA methylation profile was explored during trap induction of Arthrobotrys oligospora, a typical NTF that captures nematodes by developing adhesive networks. Whole-genome bisulfite sequencing identified 871 methylation sites and 1979 differentially methylated regions (DMRs). This first-of-its-kind investigation unveiled the widespread presence of methylation systems in NTF, and suggested potential regulation of ribosomal RNAs through DNA methylation. Functional analysis indicated DNA methylation's involvement in complex gene regulations during trap induction, impacting multiple biological processes like response to stimulus, transporter activity, cell reproduction and molecular function regulator. These findings provide a glimpse into the important roles of DNA methylation in trap induction and offer new insights for understanding the molecular mechanisms driving carnivorous lifestyle transition of NTF.

Methionine biosynthesis enzyme MoMet2 is required for rice blast fungus pathogenicity by promoting virulence gene expression via reducing 5mC modification

The emergence of fungicide resistance severely threatens crop production by limiting the availability and application of established fungicides. Therefore, it is urgent to identify new fungicidal targets for controlling plant diseases. Here, we characterized the function of a conserved homoserine O-acetyltransferase (HOA) from the rice blast fungus Magnaporthe oryzae that could serve as the candidate antifungal target. Deletion of the MoMET2 and MoCYS2 genes encoding HOAs perturbed the biosynthesis of methionine and S-adenyl methionine, a methyl group donor for epigenetic modifications, and severely attenuated the development and virulence of M. oryzae. The ∆Momet2 mutant is significantly increased in 5-methylcytosine (5mC) modification that represses the expression of genes required for pathogenicity, including MoGLIK and MoCDH-CYT. We further showed that host-induced gene silencing (HIGS) targeting MoMET2 and MoCYS2 effectively controls rice blasts. Our studies revealed the importance of HOA in the development and virulence of M. oryzae, which suggests the potential feasibility of HOA as new targets for novel anti-rice blast measurements.

Whole-genome DNA methylome analysis of different developmental stages of the entomopathogenic fungus Beauveria bassiana NCHU-157 by nanopore sequencing

The entomopathogenic fungus (EPF), Beauveria bassiana, is an important and commonly used EPF for microbial control. However, the role of DNA methylation has not been thoroughly studied. Therefore, the whole genomic DNA methylome of one promising EPF isolate, B. bassiana NCHU-157 (Bb-NCHU-157), was investigated by Oxford Nanopore Technologies (ONT). First, the whole genome of Bb-NCHU-157 was sequenced by next-generation sequencing (NGS) and ONT. The genome of Bb-NCHU-157 contains 16 contigs with 34.19 Mb and 50% GC content, which are composed of 10,848 putative protein-coding genes. Two putative DNA methyltransferases (DNMTs) were found, including Dim-2 and C-5 cytosine-specific DNA methylases. Both DNMTs showed higher expression levels in the mycelium stage than in the conidia stage, indicating that development of DNA methylation in Bb-NCHU-157 might occur in the mycelium stage. The global methylation level of the mycelium stage (5 mC = 4.56%, CG = 3.33%, CHG = 0.74%, CHH = 0.49%) was higher than that of the conidial stage (5 mC = 2.99%, CG = 1.99%, CHG = 0.63%, CHH = 0.37%) in both the gene and transposable element (TE) regions. Furthermore, the TE regions showed higher methylation frequencies than the gene regions, especially for CHH site methylation, suggesting regulation of genomic stabilization during mycelium development. In the gene regions, high methylation frequencies were found around the transcription start site (TSS) and transcription end site (TES). Moreover, CG and CHG methylation mainly occur in the promoter and intergenic regions, while CHH methylation occurs in the TE region. Among the methylated regions, 371, 661, and 756 differentially DNA methylated regions (DMRs) were hypermethylated in the mycelium in CG, CHG, and CHH, while only 13 and 7 DMRs were hypomethylated in the mycelium in CHG, and CHH, respectively. Genes located in the DMR shared the GO terms, DNA binding (GO: 0003677), and sequence-specific DNA binding (GO: 0043565) for hypermethylation in the mycelium, suggesting that methylation might regulate gene expression from the initial process. Evaluation of the DNA methylome in Bb-NCHU-157 by ONT provided new insight into this field. These data will be further validated, and epigenetic regulation during the development of B. bassiana will be explored.

The histone code of the fungal genus Aspergillus uncovered by evolutionary and proteomic analyses

Chemical modifications of DNA and histone proteins impact the organization of chromatin within the nucleus. Changes in these modifications, catalysed by different chromatin-modifying enzymes, influence chromatin organization, which in turn is thought to impact the spatial and temporal regulation of gene expression. While combinations of different histone modifications, the histone code, have been studied in several model species, we know very little about histone modifications in the fungal genus Aspergillus, whose members are generally well studied due to their importance as models in cell and molecular biology as well as their medical and biotechnological relevance. Here, we used phylogenetic analyses in 94 Aspergilli as well as other fungi to uncover the occurrence and evolutionary trajectories of enzymes and protein complexes with roles in chromatin modifications or regulation. We found that these enzymes and complexes are highly conserved in Aspergilli, pointing towards a complex repertoire of chromatin modifications. Nevertheless, we also observed few recent gene duplications or losses, highlighting Aspergillus species to further study the roles of specific chromatin modifications. SET7 (KMT6) and other components of PRC2 (Polycomb Repressive Complex 2), which is responsible for methylation on histone H3 at lysine 27 in many eukaryotes including fungi, are absent in Aspergilli as well as in closely related Penicillium species, suggesting that these lost the capacity for this histone modification. We corroborated our computational predictions by performing untargeted MS analysis of histone post-translational modifications in Aspergillus nidulans. This systematic analysis will pave the way for future research into the complexity of the histone code and its functional implications on genome architecture and gene regulation in fungi.

DNA Methyltransferases: From Evolution to Clinical Applications

DNA methylation is an epigenetic mark that living beings have used in different environments. The MTases family catalyzes DNA methylation. This process is conserved from archaea to eukaryotes, from fertilization to every stage of development, and from the early stages of cancer to metastasis. The family of DNMTs has been classified into DNMT1, DNMT2, and DNMT3. Each DNMT has been duplicated or deleted, having consequences on DNMT structure and cellular function, resulting in a conserved evolutionary reaction of DNA methylation. DNMTs are conserved in the five kingdoms of life: bacteria, protists, fungi, plants, and animals. The importance of DNMTs in whether methylate or not has a historical adaptation that in mammals has been discovered in complex regulatory mechanisms to develop another padlock to genomic insurance stability. The regulatory mechanisms that control DNMTs expression are involved in a diversity of cell phenotypes and are associated with pathologies transcription deregulation. This work focused on DNA methyltransferases, their biology, functions, and new inhibitory mechanisms reported. We also discuss different approaches to inhibit DNMTs, the use of non-coding RNAs and nucleoside chemical compounds in recent studies, and their importance in biological, clinical, and industry research.

Post-translational modifications drive secondary metabolite biosynthesis in Aspergillus: a review

Post‐translational modifications (PTMs) are important for protein function and regulate multiple cellular processes and secondary metabolites (SMs) in fungi. Aspergillus species belong to a genus renown for an abundance of bioactive secondary metabolites, many important as toxins, pharmaceuticals and in industrial production. The genes required for secondary metabolites are typically co‐localized in biosynthetic gene clusters (BGCs), which often localize in heterochromatic regions of genome and are ‘turned off’ under laboratory condition. Efforts have been made to ‘turn on’ these BGCs by genetic manipulation of histone modifications, which could convert the heterochromatic structure to euchromatin. Additionally, non‐histone PTMs also play critical roles in the regulation of secondary metabolism. In this review, we collate the known roles of epigenetic and PTMs on Aspergillus SM production. We also summarize the proteomics approaches and bioinformatics tools for PTM identification and prediction and provide future perspectives on the emerging roles of PTM on regulation of SM biosynthesis in Aspergillus and other fungi.

Diversity of Fungal DNA Methyltransferases and Their Association With DNA Methylation Patterns

DNA methyltransferases (DNMTs) are a group of proteins that catalyze DNA methylation by transferring a methyl group to DNA. The genetic variation in DNMTs results in differential DNA methylation patterns associated with various biological processes. In fungal species, DNMTs and their DNA methylation profiles were found to be very diverse and have gained many research interests. We reviewed fungal DNMTs in terms of their biological functions, protein domain structures, and their associated epigenetic regulations compared to those known in plant and animal systems. In addition, we summarized recent reports on potential RNA-directed DNA methylation (RdDM) related to DNMT5 in fungi. We surveyed up to 40 fungal species with published genome-wide DNA methylation profiles (methylomes) and presented the associations between the specific patterns of fungal DNA methylation and their DNMTs based on a phylogenetic tree of protein domain structures. For example, the main DNMTs in Basidiomycota, DNMT1 with RFD domain + DNMT5, contributing to CG methylation preference, were distinct from RID + Dim-2 in Ascomycota, resulting in a non-CG methylation preference. Lastly, we revealed that the dynamic methylation involved in fungal life stage changes was particularly low in mycelium and DNA methylation was preferentially located in transposable elements (TEs). This review comprehensively discussed fungal DNMTs and methylomes and their connection with fungal development and taxonomy to present the diverse usages of DNA methylation in fungal genomes.

DNA Methylation Is Responsive to the Environment and Regulates the Expression of Biosynthetic Gene Clusters, Metabolite Production, and Virulence in Fusarium graminearum

Histone modifications play a significant role in the regulation of biosynthetic gene clusters (BGCs) in the phytopathogen Fusarium graminearum, by contrast, epigenetic regulation by DNA methyltransferases (DNMTs) is less documented. In this study, we characterized two DNMTs (FgDIM-2 and FgRID) in F. graminearum, with homologies to "Deficient in methylation" (DIM-2) and "Repeat-induced point (RIP) deficient" (RID) from Neurospora. The loss of DNMTs resulted in not only a decrease in average methylation density in the nutrient-poor, compared to nutrient-rich conditions, but also differences in the genes expressed between the WT and the DNMT mutant strains, implicating the external environment as an important trigger in altering DNA methylation patterns. Consequently, we observed significant changes in the regulation of multiple BGCs and alterations in the pathogenicity of the fungus.

Repetitive Elements Contribute to the Diversity and Evolution of Centromeres in the Fungal Genus Verticillium

Centromeres are chromosomal regions that are crucial for chromosome segregation during mitosis and meiosis, and failed centromere formation can contribute to chromosomal anomalies. Despite this conserved function, centromeres differ significantly between and even within species. Thus far, systematic studies into the organization and evolution of fungal centromeres remain scarce. In this study, we identified the centromeres in each of the 10 species of the fungal genus Verticillium and characterized their organization and evolution. Chromatin immunoprecipitation of the centromere-specific histone CenH3 (ChIP-seq) and chromatin conformation capture (Hi-C) followed by high-throughput sequencing identified eight conserved, large (∼150-kb), AT-, and repeat-rich regional centromeres that are embedded in heterochromatin in the plant pathogen Verticillium dahliae Using Hi-C, we similarly identified repeat-rich centromeres in the other Verticillium species. Strikingly, a single degenerated long terminal repeat (LTR) retrotransposon is strongly associated with centromeric regions in some but not all Verticillium species. Extensive chromosomal rearrangements occurred during Verticillium evolution, of which some could be linked to centromeres, suggesting that centromeres contributed to chromosomal evolution. The size and organization of centromeres differ considerably between species, and centromere size was found to correlate with the genome-wide repeat content. Overall, our study highlights the contribution of repetitive elements to the diversity and rapid evolution of centromeres within the fungal genus VerticilliumIMPORTANCE The genus Verticillium contains 10 species of plant-associated fungi, some of which are notorious pathogens. Verticillium species evolved by frequent chromosomal rearrangements that contribute to genome plasticity. Centromeres are instrumental for separation of chromosomes during mitosis and meiosis, and failed centromere functionality can lead to chromosomal anomalies. Here, we used a combination of experimental techniques to identify and characterize centromeres in each of the Verticillium species. Intriguingly, we could strongly associate a single repetitive element to the centromeres of some of the Verticillium species. The presence of this element in the centromeres coincides with increased centromere sizes and genome-wide repeat expansions. Collectively, our findings signify a role of repetitive elements in the function, organization, and rapid evolution of centromeres in a set of closely related fungal species.

Epigenetic Regulation of Verticillium dahliae Virulence: Does DNA Methylation Level Play A Role?

Verticillium dahliae is the etiological agent of Verticillium wilt of olive. The virulence of Defoliating V. dahliae isolates usually displays differences and high plasticity. This work studied whether an epigenetic mechanism was involved in this plasticity. An inverse correlation between virulence and DNA methylation of protein-coding genes was found. A set of 831 genes was selected for their highly consistent inverse methylation profile and virulence in the five studied isolates. Of these genes, ATP-synthesis was highly represented, which indicates that the more virulent D isolates are, the more energy requirements they may have. Furthermore, there were numerous genes in the protein biosynthesis process: genes coding for the chromatin structure, which suggests that epigenetic changes may also affect chromatin condensation; many transmembrane transporter genes, which is consistent with denser compounds, traffic through membranes in more virulent isolates; a fucose-specific lectin that may play a role in the attachment to plant cell walls during the host infection process; and pathogenic cutinases that facilitate plant invasion and sporulation genes for rapid spreading alongside plants. Our findings support the notion that differences in the virulence of the Defoliating V. dahliae isolates may be controlled, at least to some extent, by an epigenetic mechanism.

The Pattern and Function of DNA Methylation in Fungal Plant Pathogens

To successfully infect plants and trigger disease, fungal plant pathogens use various strategies that are dependent on characteristics of their biology and genomes. Although pathogenic fungi are different from animals and plants in the genomic heritability, sequence feature, and epigenetic modification, an increasing number of phytopathogenic fungi have been demonstrated to share DNA methyltransferases (MTases) responsible for DNA methylation with animals and plants. Fungal plant pathogens predominantly possess four types of DNA MTase homologs, including DIM-2, DNMT1, DNMT5, and RID. Numerous studies have indicated that DNA methylation in phytopathogenic fungi mainly distributes in transposable elements (TEs), gene promoter regions, and the repetitive DNA sequences. As an important and heritable epigenetic modification, DNA methylation is associated with silencing of gene expression and transposon, and it is responsible for a wide range of biological phenomena in fungi. This review highlights the relevant reports and insights into the important roles of DNA methylation in the modulation of development, pathogenicity, and secondary metabolism of fungal plant pathogens. Recent evidences prove that there are massive links between DNA and histone methylation in fungi, and they commonly regulate fungal development and mycotoxin biosynthesis.

Evaluation of DNA Methylation Changes by CRED-RA Analysis Following Prednisone Treatment of Endophyte, Fusarium oxysporum

Endophytes that represent a sub-set of plant resident microbes are a reservoir of bioactive metabolites. Many of the secondary metabolite biosynthetic gene clusters of endophytes are silent under axenic culture conditions. Epigenetic reprogramming of such cryptic pathways is possible by use of small molecule modulators like prednisone. Methylation changes induced by prednisone, a hypomethylating epigenetic modulator were studied in endophytic Fusarium oxysporum. CRED-RA analysis following exposure to non-cytotoxic dose (300 µM) revealed prednisone as effective in inducing non-methylation and semi-methylation pattern while inhibiting full-methylation of the genome. Effectiveness of prednisone as a DNA methyl transferase inhibitor can be explored in future to study alterations in secondary metabolite gene expression profile in endophytic F. oxysporum.

Revealing mitogenome-wide DNA methylation and RNA editing of three Ascomycotina fungi using SMRT sequencing

Beauveria bassiana, Cordyceps militaris and Ophiocordyceps sinensis (Ascomycotina) are traditional Chinese medicines. Here, mitogenomes of these three Ascomycotina fungi were sequenced and de-novo assembled using single-molecule real-time sequencing. The results showed that their complete mitogenomes were 31,258, 31,854 and 157,584 bp, respectively, with sequencing depth approximately 278,760×, 326,283× and 69,385×. Types of repeat sequences were mainly (AA)n, (AAT)n, (TA)n and (TATT)n. DNA methylation motifs were revealed in DNA modifications of these three fungi. We discovered new models of RNA editing through analysis of transcriptomes from B. bassiana and C. militaris. These data lay a solid foundation for further genetic and biological studies about these three fungi, especially for elucidating the mitogenome evolution and exploring the regulatory mechanism of adapting environment.

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

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⁻¹¹ (A. flavus), 1.1 × 10⁻¹¹ (A. fumigatus) and 4.1 × 10⁻¹¹ (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 end-joining (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⁻¹¹ (A. flavus), 2.2 × 10⁻¹¹ (A. fumigatus) and 4.5 × 10⁻¹¹ (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.

Chromatin-Mediated Regulation of Genome Plasticity in Human Fungal Pathogens

Human fungal pathogens, such as Candida albicans, Aspergillus fumigatus and Cryptococcus neoformans, are a public health problem, causing millions of infections and killing almost half a million people annually. The ability of these pathogens to colonise almost every organ in the human body and cause life-threating infections relies on their capacity to adapt and thrive in diverse hostile host-niche environments. Stress-induced genome instability is a key adaptive strategy used by human fungal pathogens as it increases genetic diversity, thereby allowing selection of genotype(s) better adapted to a new environment. Heterochromatin represses gene expression and deleterious recombination and could play a key role in modulating genome stability in response to environmental changes. However, very little is known about heterochromatin structure and function in human fungal pathogens. In this review, I use our knowledge of heterochromatin structure and function in fungal model systems as a road map to review the role of heterochromatin in regulating genome plasticity in the most common human fungal pathogens: Candida albicans, Aspergillus fumigatus and Cryptococcus neoformans.

Chromatin-dependent regulation of secondary metabolite biosynthesis in fungi: is the picture complete?

Fungal secondary metabolites are small molecules that exhibit diverse biological activities exploited in medicine, industry and agriculture. Their biosynthesis is governed by co-expressed genes that often co-localize in gene clusters. Most of these secondary metabolite gene clusters are inactive under laboratory conditions, which is due to a tight transcriptional regulation. Modifications of chromatin, the complex of DNA and histone proteins influencing DNA accessibility, play an important role in this regulation. However, tinkering with well-characterised chemical and genetic modifications that affect chromatin alters the expression of only few biosynthetic gene clusters, and thus the regulation of the vast majority of biosynthetic pathways remains enigmatic. In the past, attempts to activate silent gene clusters in fungi mainly focused on histone acetylation and methylation, while in other eukaryotes many other post-translational modifications are involved in transcription regulation. Thus, how chromatin regulates the expression of gene clusters remains a largely unexplored research field. In this review, we argue that focusing on only few well-characterised chromatin modifications is significantly hampering our understanding of the chromatin-based regulation of biosynthetic gene clusters. Research on underexplored chromatin modifications and on the interplay between different modifications is timely to fully explore the largely untapped reservoir of fungal secondary metabolites.

Diversity of cytosine methylation across the fungal tree of life

The generation of thousands of fungal genomes is leading to a better understanding of genes and genomic organization within the kingdom. However, the epigenome, which includes DNA and chromatin modifications, remains poorly investigated in fungi. Large comparative studies in animals and plants have deepened our understanding of epigenomic variation, particularly of the modified base 5-methylcytosine (5mC), but taxonomic sampling of disparate groups is needed to develop unifying explanations for 5mC variation. Here we utilize the largest phylogenetic resolution of 5mC methyltransferases (5mC MTases) and genome evolution to better understand levels and patterns of 5mC across fungi. We show that extant 5mC MTase genotypes are descendent from ancestral maintenance and de novo genotypes, whereas the 5mC MTases DIM-2 and RID are more recently derived, and that 5mC levels are correlated with 5mC MTase genotype and transposon content. Our survey also revealed that fungi lack canonical gene body methylation, which distinguishes fungal epigenomes from certain insect and plant species. However, some fungal species possess independently derived clusters of contiguous 5mC encompassing many genes. In some cases, DNA repair pathways and the N6-methyladenine (6mA) DNA modification negatively coevolved with 5mC pathways, which additionally contributed to interspecific epigenomic variation across fungi.

Genome-wide analysis of DNA methylation in subcultured Cordyceps militaris

The entomopathogenic mushroom Cordyceps militaris is a storehouse of various medicinal compounds and pharmacological effects. However, the high frequency of strain degeneration during subculture and preservation severely limits the large-scale production of C. militaris. DNA methylation is an important epigenomic modification involved in gene regulation. In this study, we used bisulfite sequencing for DNA methylation profiling of wild-type and mutant C. militaris. The differentially methylated regions (DMRs) of the two types were analyzed using Gene Ontology (GO) clustering and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. DNA methylation levels of the wild-type and mutant-type C. militaris were 0.48% and 0.56%, respectively. Methylation appeared at CHH dinucleotides in 58.62% and 58.20% of all methylated cytosine sites in the wild and mutant types, respectively. In all, 188 DMRs were identified from the wild and mutant types. Most of the DMRs ranged from 200 to 350 bp in length. KEGG pathways of the expression of DMR-related genes, which are involved in pyruvate metabolism, glycerophospholipid metabolism, DNA replication, and N-glycan biosynthesis. This contributes to the knowledge and understanding of the possible mechanisms of C. militaris strain degeneration.

Epigenetic and Posttranslational Modifications in Regulating the Biology of Aspergillus Species

Epigenetic and posttranslational modifications have been proved to participate in multiple cellular processes and suggested to be an important regulatory mechanism on transcription of genes in eukaryotes. However, our knowledge about epigenetic and posttranslational modifications mainly comes from the studies of yeasts, plants, and animals. Recently, epigenetic and posttranslational modifications have also raised concern for the relevance of regulating fungal biology in Aspergillus. Emerging evidence indicates that these modifications could be a connection between genetic elements and environmental factors, and their combined effects may finally lead to fungal phenotypical changes. This article describes the advances in typical DNA and protein modifications in the genus Aspergillus, focusing on methylation, acetylation, phosphorylation, ubiquitination, sumoylation, and neddylation.

Comparative Methylome Analysis Reveals Perturbation of Host Epigenome in Chestnut Blight Fungus by a Hypovirus

In eukaryotic genomes, DNA methylation is an important type of epigenetic modification that plays crucial roles in many biological processes. To investigate the impact of a hypovirus infection on the methylome of Cryphonectria parasitica, the chestnut blight fungus, whole-genome bisulfite sequencing (WGBS) was employed to generate single-base resolution methylomes of the fungus with/without hypovirus infection. The results showed that hypovirus infection alters methylation in all three contexts (CG, CHG, and CHH), especially in gene promoters. A total of 600 differentially methylated regions (DMRs) were identified, of which 144 could be annotated to functional genes. RNA-seq analysis revealed that DNA methylation in promoter is negatively correlated with gene expression. Among DMRs, four genes were shown to be involved in conidiation, orange pigment production, and virulence. Taken together, our DNA methylomes analysis provide valuable insights into the understanding of the relationship between DNA methylation and hypovirus infection, as well as phenotypic traits in C. parasitica.

Aspergillus flavus Secondary Metabolites: More than Just Aflatoxins

Aspergillus flavus is best known for producing the family of potent carcinogenic secondary metabolites known as aflatoxins. However, this opportunistic plant and animal pathogen also produces numerous other secondary metabolites, many of which have also been shown to be toxic. While about forty of these secondary metabolites have been identified from A. flavus cultures, analysis of the genome has predicted the existence of at least 56 secondary metabolite gene clusters. Many of these gene clusters are not expressed during growth of the fungus on standard laboratory media. This presents researchers with a major challenge of devising novel strategies to manipulate the fungus and its genome so as to activate secondary metabolite gene expression and allow identification of associated cluster metabolites. In this review, we discuss the genetic, biochemical and bioinformatic methods that are being used to identify previously uncharacterized secondary metabolite gene clusters and their associated metabolites. It is important to identify as many of these compounds as possible to determine their bioactivity with respect to fungal development, survival, virulence and especially with respect to any potential synergistic toxic effects with aflatoxin.

Advances in molecular and genomic research to safeguard food and feed supply from aflatoxin contamination

Worldwide recognition that aflatoxin contamination of agricultural commodities by the fungus Aspergillus flavus is a global problem has significantly benefitted from global collaboration for understanding the contaminating fungus, as well as for developing and implementing solutions against the contamination. The effort to address this serious food and feed safety issue has led to a detailed understanding of the taxonomy, ecology, physiology, genomics and evolution of A. flavus, as well as strategies to reduce or control pre-harvest aflatoxin contamination, including (1) biological control, using atoxigenic aspergilli, (2) proteomic and genomic analyses for identifying resistance factors in maize as potential breeding markers to enable development of resistant maize lines, and (3) enhancing host-resistance by bioengineering of susceptible crops, such as cotton, maize, peanut and tree nuts. A post-harvest measure to prevent the occurrence of aflatoxin contamination in storage is also an important component for reducing exposure of populations worldwide to aflatoxins in food and feed supplies. The effect of environmental changes on aflatoxin contamination levels has recently become an important aspect for study to anticipate future contamination levels. The ability of A. flavus to produce dozens of secondary metabolites, in addition to aflatoxins, has created a new avenue of research for understanding the role these metabolites play in the survival and biodiversity of this fungus. The understanding of A. flavus, the aflatoxin contamination problem, and control measures to prevent the contamination has become a unique example for an integrated approach to safeguard global food and feed safety.

Global DNA Methylation in the Chestnut Blight Fungus Cryphonectria parasitica and Genome-Wide Changes in DNA Methylation Accompanied with Sectorization

Mutation in CpBck1, an ortholog of the cell wall integrity mitogen-activated protein kinase kinase kinase (MAPKKK) of Saccharomyces cerevisiae, in the chestnut blight fungus Cryphonectria parasitica resulted in a sporadic sectorization as culture proceeded. The progeny from the sectored area maintained the characteristics of the sector, showing a massive morphogenetic change, including robust mycelial growth without differentiation. Epigenetic changes were investigated as the genetic mechanism underlying this sectorization. Quantification of DNA methylation and whole-genome bisulfite sequencing revealed genome-wide DNA methylation of the wild-type at each nucleotide level and changes in DNA methylation of the sectored progeny. Compared to the wild-type, the sectored progeny exhibited marked genome-wide DNA hypomethylation but increased methylation sites. Expression analysis of two DNA methyltransferases, including two representative types of DNA methyltransferase (DNMTase), demonstrated that both were significantly down-regulated in the sectored progeny. However, functional analysis using mutant phenotypes of corresponding DNMTases demonstrated that a mutant of CpDmt1, an ortholog of RID of Neurospora crassa, resulted in the sectored phenotype but the CpDmt2 mutant did not, suggesting that the genetic basis of fungal sectorization is more complex. The present study revealed that a mutation in a signaling pathway component resulted in sectorization accompanied with changes in genome-wide DNA methylation, which suggests that this signal transduction pathway is important for epigenetic control of sectorization via regulation of genes involved in DNA methylation.

Epigenetic regulation of development and pathogenesis in fungal plant pathogens

Evidently, epigenetics is at forefront in explaining the mechanisms underlying the success of human pathogens and in the identification of pathogen-induced modifications within host plants. However, there is a lack of studies highlighting the role of epigenetics in the modulation of the growth and pathogenicity of fungal plant pathogens. In this review, we attempt to highlight and discuss the role of epigenetics in the regulation of the growth and pathogenicity of fungal phytopathogens using Magnaporthe oryzae, a devastating fungal plant pathogen, as a model system. With the perspective of wide application in the understanding of the development, pathogenesis and control of other fungal pathogens, we attempt to provide a synthesized view of the epigenetic studies conducted on M. oryzae to date. First, we discuss the mechanisms of epigenetic modifications in M. oryzae and their impact on fungal development and pathogenicity. Second, we highlight the unexplored epigenetic mechanisms and areas of research that should be considered in the near future to construct a holistic view of epigenetic functioning in M. oryzae and other fungal plant pathogens. Importantly, the development of a complete understanding of the modulation of epigenetic regulation in fungal pathogens can help in the identification of target points to combat fungal pathogenesis.

Widespread 5-methylcytosine in the genomes of avian Coccidia and other apicomplexan parasites detected by an ELISA-based method

To date, little is known about cytosine methylation in the genomic DNA of apicomplexan parasites, although it has been confirmed that this important epigenetic modification exists in many lower eukaryotes, plants, and animals. In the present study, ELISA-based detection demonstrated that low levels of 5-methylcytosine (5-mC) are present in Eimeria spp., Toxoplasma gondii, Cryptosporidium spp., and Neospora caninum. The proportions of 5-mC in genomic DNA were 0.18 ± 0.02% in E tenella sporulated oocysts, 0.19 ± 0.01% in E. tenella second-generation merozoites, 0.22 ± 0.04% in T. gondii tachyzoites, 0.28 ± 0.03% in N. caninum tachyzoites, and 0.06 ± 0.01, 0.11 ± 0.01, and 0.09 ± 0.01% in C. andersoni, C. baileyi, and C. parvum sporulated oocysts, respectively. In addition, we found that the percentages of 5-mC in E. tenella varied considerably at different life stages, with sporozoites having the highest percentage of 5-mC (0.78 ± 0.10%). Similar stage differences in 5-mC were also found in E. maxima, E. necatrix, and E. acervulina, the levels of 5-mC in their sporozoites being 4.3-, 1.8-, 2.5-, and 2.0-fold higher than that of sporulated oocysts, respectively (p < 0.01). Furthermore, a total DNA methyltransferase-like activity was detected in whole cell extracts prepared from E. tenella sporozoites. In conclusion, genomic DNA methylation is present in these apicomplexan parasites and may play a role in the stage conversion of Eimeria.

Characterization of Cytosine Methylation and the DNA Methyltransferases of Toxoplasma gondii

DNA methylation is a key epigenetic modification which confers phenotypic plasticity and adaptation. Cyst-forming strains of Toxoplasma gondii undergo tachyzoite to bradyzoite conversion after initial acute infection of a host, and the reverse conversion may occur in immune-suppressed hosts. The formation of m⁵C is catalyzed by DNA methyltransferase (DNMT). We identified two functional DNA methyltransferases, TgDNMTa and TgDNMTb, in T. gondii that may mediate DNA methylation. The recombinant proteins showed intrinsic methyltransferase activity; both have higher transcription levels in bradyzoites than that in tachyzoites. We performed genome-wide analysis of DNA methylation in tachyzoites and bradyzoites. The results showed more methylation sites in bradyzoites than that in tachyzoites. The most significantly enriched GO-terms of genes with DNA methylation were associated with basal cellular processes such as energy metabolism and parasite resistance to host immunity. Tachyzoite proliferation in parasitophorous vacuoles (PV) can be inhibited by the DNA methyltransferase inhibitor 5-azacytidine, a chemical analogue of the nucleotide cytosine that can inactivate DNA methyltransferases. These findings provide the first confirmation of DNA methylation in T. gondii.

DNA methyltransferases contribute to the fungal development, stress tolerance and virulence of the entomopathogenic fungus Metarhizium robertsii

DNA methylation is an important epigenetic mark in mammals, plants, and fungi and depends on multiple genetic pathways involving de novo and maintenance DNA methyltransferases (DNMTases). Metarhizium robertsii, a model system for investigating insect-fungus interactions, has been used as an environmentally friendly alternative to chemical insecticides. However, little is known concerning the molecular basis for DNA methylation. Here, we report on the roles of two DNMTases (MrRID and MrDIM-2) by characterizing ΔMrRID, ΔMrDIM-2, and ΔRID/ΔDIM-2 mutants. The results showed that approximately 71, 10, and 8% of ^mC sites remained in the ΔMrRID, ΔMrDIM-2, and ΔRID/ΔDIM-2 strains, respectively, compared with the wild-type (WT) strain. Further analysis showed that MrRID regulates the specificity of DNA methylation and MrDIM-2 is responsible for most DNA methylation, implying an interaction or cooperation between MrRID and MrDIM-2 for DNA methylation. Moreover, the ΔMrDIM-2 and ΔRID/ΔDIM-2 strains showed more defects in radial growth and conidial production compared to the WT. Under ultraviolet (UV) irradiation or heat stress, an obvious reduction in spore viability was observed for all the mutant strains compared to the WT. The spore median lethal times (LT₅₀s) for the ΔMrDIM-2 and ΔRID/ΔDIM-2 strains in the greater wax moth, Galleria mellonella, were decreased by 47.7 and 65.9%, respectively, which showed that MrDIM-2 is required for full fungal virulence. Our data advances the understanding of the function of DNMTase in entomopathogenic fungi, which should contribute to future epigenetic investigations in fungi.

Comparative genomics reveals high biological diversity and specific adaptations in the industrially and medically important fungal genus Aspergillus

BACKGROUND

The fungal genus Aspergillus is of critical importance to humankind. Species include those with industrial applications, important pathogens of humans, animals and crops, a source of potent carcinogenic contaminants of food, and an important genetic model. The genome sequences of eight aspergilli have already been explored to investigate aspects of fungal biology, raising questions about evolution and specialization within this genus.

RESULTS

We have generated genome sequences for ten novel, highly diverse Aspergillus species and compared these in detail to sister and more distant genera. Comparative studies of key aspects of fungal biology, including primary and secondary metabolism, stress response, biomass degradation, and signal transduction, revealed both conservation and diversity among the species. Observed genomic differences were validated with experimental studies. This revealed several highlights, such as the potential for sex in asexual species, organic acid production genes being a key feature of black aspergilli, alternative approaches for degrading plant biomass, and indications for the genetic basis of stress response. A genome-wide phylogenetic analysis demonstrated in detail the relationship of the newly genome sequenced species with other aspergilli.

CONCLUSIONS

Many aspects of biological differences between fungal species cannot be explained by current knowledge obtained from genome sequences. The comparative genomics and experimental study, presented here, allows for the first time a genus-wide view of the biological diversity of the aspergilli and in many, but not all, cases linked genome differences to phenotype. Insights gained could be exploited for biotechnological and medical applications of fungi.

Reduced representation bisulphite sequencing of ten bovine somatic tissues reveals DNA methylation patterns and their impacts on gene expression

Background

As a major epigenetic component, DNA methylation plays important functions in individual development and various diseases. DNA methylation has been well studied in human and model organisms, but only limited data exist in economically important animals like cattle.

Results

Using reduced representation bisulphite sequencing (RRBS), we obtained single-base-resolution maps of bovine DNA methylation from ten somatic tissues. In total, we evaluated 1,868,049 cytosines in CG-enriched regions. While we found slightly low methylation levels (29.87 to 38.06 %) in cattle, the methylation contexts (CGs and non-CGs) of cattle showed similar methylation patterns to other species. Non-CG methylation was detected but methylation levels in somatic tissues were significantly lower than in pluripotent cells. To study the potential function of the methylation, we detected 10,794 differentially methylated cytosines (DMCs) and 836 differentially methylated CG islands (DMIs). Further analyses in the same tissues revealed many DMCs (including non-CGs) and DMIs, which were highly correlated with the expression of genes involved in tissue development.

Conclusions

In summary, our study provides a baseline dataset and essential information for DNA methylation profiles of cattle.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3116-1) contains supplementary material, which is available to authorized users.

The DmtA methyltransferase contributes to Aspergillus flavus conidiation, sclerotial production, aflatoxin biosynthesis and virulence

DNA methylation is essential for epigenetic regulation of gene transcription and development in many animals, plants and fungi. We investigated whether DNA methylation plays a role in the development and secondary metabolism of Aspergillus flavus, identified the DmtA methyltransferase from A. flavus, and produced a dmtA knock-out mutant by replacing the dmtA coding sequence with the pyrG selectable marker. The A. flavus dmtA null mutant lines produced white fluffy mycelium in liquid medium, and displayed a slightly flavescent conidial pigmentation compared with the normal yellow of the wild-type strain when grown on agar. The ΔdmtA lines exhibited decreased conidiation and aflatoxin (AF) biosynthesis, compared with the wild-type line, suggesting that the DmtA knock-out affected the transcriptional level of genes in the AF cluster. In particular, sclerotia development and host colonization were altered in the dmtA null mutants. Green fluorescent protein tagging at the C-terminus of DmtA showed that DmtA localized to the nucleus and cytoplasm. DNA methylation content measurements in the dmtA mutants revealed no widespread DNA methylation in the mutants or wild-type lines. Thus, our findings suggest that DmtA, apart from being a C-5 cytosine methyltransferase in A. flavus, contributes to asexual development, aflatoxin biosynthesis, sclerotial production and virulence.

Systems Biology Approaches for Host-Fungal Interactions: An Expanding Multi-Omics Frontier

Opportunistic fungal infections are an increasing threat for global health, and for immunocompromised patients in particular. These infections are characterized by interaction between fungal pathogen and host cells. The exact mechanisms and the attendant variability in host and fungal pathogen interaction remain to be fully elucidated. The field of systems biology aims to characterize a biological system, and utilize this knowledge to predict the system's response to stimuli such as fungal exposures. A multi-omics approach, for example, combining data from genomics, proteomics, metabolomics, would allow a more comprehensive and pan-optic "two systems" biology of both the host and the fungal pathogen. In this review and literature analysis, we present highly specialized and nascent methods for analysis of multiple -omes of biological systems, in addition to emerging single-molecule visualization techniques that may assist in determining biological relevance of multi-omics data. We provide an overview of computational methods for modeling of gene regulatory networks, including some that have been applied towards the study of an interacting host and pathogen. In sum, comprehensive characterizations of host-fungal pathogen systems are now possible, and utilization of these cutting-edge multi-omics strategies may yield advances in better understanding of both host biology and fungal pathogens at a systems scale.

Epigenetic Segregation of Microbial Genomes from Complex Samples Using Restriction Endonucleases HpaII and McrB

We describe continuing work to develop restriction endonucleases as tools to enrich targeted genomes of interest from diverse populations. Two approaches were developed in parallel to segregate genomic DNA based on cytosine methylation. First, the methyl-sensitive endonuclease HpaII was used to bind non-CG methylated DNA. Second, a truncated fragment of McrB was used to bind CpG methylated DNA. Enrichment levels of microbial genomes can exceed 100-fold with HpaII allowing improved genomic detection and coverage of otherwise trace microbial genomes from sputum. Additionally, we observe interesting enrichment results that correlate with the methylation states not only of bacteria, but of fungi, viruses, a protist and plants. The methods presented here offer promise for testing biological samples for pathogens and global analysis of population methylomes.

COBRA-Seq: Sensitive and Quantitative Methylome Profiling

Combined Bisulfite Restriction Analysis (COBRA) quantifies DNA methylation at a specific locus. It does so via digestion of PCR amplicons produced from bisulfite-treated DNA, using a restriction enzyme that contains a cytosine within its recognition sequence, such as TaqI. Here, we introduce COBRA-seq, a genome wide reduced methylome method that requires minimal DNA input (0.1-1.0 mg) and can either use PCR or linear amplification to amplify the sequencing library. Variants of COBRA-seq can be used to explore CpG-depleted as well as CpG-rich regions in vertebrate DNA. The choice of enzyme influences enrichment for specific genomic features, such as CpG-rich promoters and CpG islands, or enrichment for less CpG dense regions such as enhancers. COBRA-seq coupled with linear amplification has the additional advantage of reduced PCR bias by producing full length fragments at high abundance. Unlike other reduced representative methylome methods, COBRA-seq has great flexibility in the choice of enzyme and can be multiplexed and tuned, to reduce sequencing costs and to interrogate different numbers of sites. Moreover, COBRA-seq is applicable to non-model organisms without the reference genome and compatible with the investigation of non-CpG methylation by using restriction enzymes containing CpA, CpT, and CpC in their recognition site.

Genome-wide analysis of DNA methylation in the sexual stage of the insect pathogenic fungus Cordyceps militaris

DNA methylation is a basic epigenetic mechanism found in eukaryotes, but its patterns and roles vary significantly among diverse taxa. In fungi, DNA methylation has various effects on diverse biological processes. However, its function in the sexual development of fungi remains unclear. Cordyceps militaris, readily performs sexual reproduction and thus provides a remarkably rich model for understanding epigenetic processes in sexual development. Here, we surveyed the methylome of C. militaris at single-base resolution to assess DNA methylation patterns during sexual development using genomic bisulfite sequencing (BS-Seq). The results showed that approximately 0.4 % of cytosines are methylated, similar to the DNA methylation level (0.39 %) during asexual development. Importantly, we found that DNA methylation in the fungi undergoes global reprogramming during fungal development. Moreover, RNA-Seq analysis indicated that the differentially methylated regions (DMRs) have no correlation with the genes that have roles during fungal sexual development in C. militaris. These results provide a comprehensive characterization of DNA methylation in the sexual development of C. militaris, which will contribute to future investigations of epigenetics in fungi.

Transcriptome changes in Fusarium verticillioides caused by mutation in the transporter-like gene FST1

BACKGROUND

Fusarium verticillioides causes an important seed disease on maize and produces the fumonisin group of mycotoxins, which are toxic to humans and livestock. A previous study discovered that a gene (FST1) in the pathogen affects fumonisin production and virulence. Although the predicted amino acid sequence of FST1 is similar to hexose transporters, previous experimental evidence failed to prove function.

RESULTS

Three new phenotypes were identified that are associated with the FST1 mutant of F. verticillioides (Δfst1), namely reduction in macroconidia production, increased sensitivity to hydrogen peroxide, and reduced mycelial hydrophobicity. A transcriptome comparison of the wild type and strain Δfst1 grown on autoclaved maize kernels for six days identified 2677 genes that were differentially expressed. Through gene ontology analysis, 961 genes were assigned to one of 12 molecular function categories. Sets of down-regulated genes in strain Δfst1 were identified that could account for each of the mutant phenotypes.

CONCLUSION

The study provides evidence that disruption of FST1 causes several metabolic and developmental defects in F. verticillioides. FST1 appears to connect the expression of several gene networks, including those involved in secondary metabolism, cell wall structure, conidiogenesis, virulence, and resistance to reactive oxygen species. The results support our hypothesis that FST1 functions within the framework of environmental sensing.

Genome-wide profiling of DNA methylation provides insights into epigenetic regulation of fungal development in a plant pathogenic fungus, Magnaporthe oryzae

DNA methylation is an important epigenetic modification that regulates development of plants and mammals. To investigate the roles of DNA methylation in fungal development, we profiled genome-wide methylation patterns at single-nucleotide resolution during vegetative growth, asexual reproduction, and infection-related morphogenesis in a model plant pathogenic fungus, Magnaporthe oryzae. We found that DNA methylation occurs in and around genes as well as transposable elements and undergoes global reprogramming during fungal development. Such reprogramming of DNA methylation suggests that it may have acquired new roles other than controlling the proliferation of TEs. Genetic analysis of DNA methyltransferase deletion mutants also indicated that proper reprogramming in methylomes is required for asexual reproduction in the fungus. Furthermore, RNA-seq analysis showed that DNA methylation is associated with transcriptional silencing of transposable elements and transcript abundance of genes in context-dependent manner, reinforcing the role of DNA methylation as a genome defense mechanism. This comprehensive approach suggests that DNA methylation in fungi can be a dynamic epigenetic entity contributing to fungal development and genome defense. Furthermore, our DNA methylomes provide a foundation for future studies exploring this key epigenetic modification in fungal development and pathogenesis.

A Bayesian framework to identify methylcytosines from high-throughput bisulfite sequencing data

High-throughput bisulfite sequencing technologies have provided a comprehensive and well-fitted way to investigate DNA methylation at single-base resolution. However, there are substantial bioinformatic challenges to distinguish precisely methylcytosines from unconverted cytosines based on bisulfite sequencing data. The challenges arise, at least in part, from cell heterozygosis caused by multicellular sequencing and the still limited number of statistical methods that are available for methylcytosine calling based on bisulfite sequencing data. Here, we present an algorithm, termed Bycom, a new Bayesian model that can perform methylcytosine calling with high accuracy. Bycom considers cell heterozygosis along with sequencing errors and bisulfite conversion efficiency to improve calling accuracy. Bycom performance was compared with the performance of Lister, the method most widely used to identify methylcytosines from bisulfite sequencing data. The results showed that the performance of Bycom was better than that of Lister for data with high methylation levels. Bycom also showed higher sensitivity and specificity for low methylation level samples (<1%) than Lister. A validation experiment based on reduced representation bisulfite sequencing data suggested that Bycom had a false positive rate of about 4% while maintaining an accuracy of close to 94%. This study demonstrated that Bycom had a low false calling rate at any methylation level and accurate methylcytosine calling at high methylation levels. Bycom will contribute significantly to studies aimed at recalibrating the methylation level of genomic regions based on the presence of methylcytosines.

High-throughput bisulfite sequencing (BS-seq) has advanced tremendously the study of DNA methylation and the determination of methylcytosines at single-base resolution. In BS-seq data analysis, sequencing errors, incomplete bisulfite conversion, and cell heterozygosis affect the accuracy of methylcytosine detection in quite a major way. Simple filtering methods using predefined thresholds have proved to have extremely low efficiency. The commonly used Lister uses binomial distribution to overcome the impacts of non-conversion rate and sequencing errors, but the impact of the cell heterozygosis is not considered. Here, we present Bycom, a novel algorithm based on the Bayesian inference model. To improve the accuracy of methylcytosine calling, Bycom considers sequencing errors, non-conversion rate, and cell heterozygosis integratively to identify methylcytosines from BS-seq data. We evaluated the performance of Bycom using different kinds of BS-seq data. Our results demonstrated that Bycom identified methylcytosines more accurately than Lister, especially in BS-seq data with extremely low genome-wide methylation levels.

Non-exhaustive DNA methylation-mediated transposon silencing in the black truffle genome, a complex fungal genome with massive repeat element content

BACKGROUND

We investigated how an extremely transposon element (TE)-rich organism such as the plant-symbiotic ascomycete truffle Tuber melanosporum exploits DNA methylation to cope with the more than 45,000 repeated elements that populate its genome.

RESULTS

Whole-genome bisulfite sequencing performed on different developmental stages reveals a high fraction of methylated cytosines with a strong preference for CpG sites. The methylation pattern is highly similar among samples and selectively targets TEs rather than genes. A marked trend toward hypomethylation is observed for TEs located within a 1 kb distance from expressed genes, rather than segregated in TE-rich regions of the genome. Approximately 300 hypomethylated or unmethylated TEs are transcriptionally active, with higher expression levels in free-living mycelium compared to fruitbody. Indeed, multiple TE-enriched, copy number variant regions bearing a significant fraction of hypomethylated and expressed TEs are found almost exclusively in free-living mycelium. A reduction of DNA methylation, restricted to non-CpG sites and accompanied by an increase in TE expression, is observed upon treatment of free-living mycelia with 5-azacytidine.

CONCLUSIONS

Evidence derived from analysis of the T. melanosporum methylome indicates that a non-exhaustive, partly reversible, methylation process operates in truffles. This allows for the existence of hypomethylated, transcriptionally active TEs that are associated with copy number variant regions of the genome. Non-exhaustive TE methylation may reflect a role of active TEs in promoting genome plasticity and the ability to adapt to sudden environmental changes.

Small chemical chromatin effectors alter secondary metabolite production in Aspergillus clavatus

The filamentous fungus Aspergillus clavatus is known to produce a variety of secondary metabolites (SM) such as patulin, pseurotin A, and cytochalasin E. In fungi, the production of most SM is strongly influenced by environmental factors and nutrients. Furthermore, it has been shown that the regulation of SM gene clusters is largely based on modulation of a chromatin structure. Communication between fungi and bacteria also triggers chromatin-based induction of silent SM gene clusters. Consequently, chemical chromatin effectors known to inhibit histone deacetylases (HDACs) and DNA-methyltransferases (DNMTs) influence the SM profile of several fungi. In this study, we tested the effect of five different chemicals, which are known to affect chromatin structure, on SM production in A. clavatus using two growth media with a different organic nitrogen source. We found that production of patulin was completely inhibited and cytochalasin E levels strongly reduced, whereas growing A. clavatus in media containing soya-derived peptone led to substantially higher pseurotin A levels. The HDAC inhibitors valproic acid, trichostatin A and butyrate, as well as the DNMT inhibitor 5-azacytidine (AZA) and N-acetyl-d-glucosamine, which was used as a proxy for bacterial fungal co-cultivation, had profound influence on SM accumulation and transcription of the corresponding biosynthetic genes. However, the repressing effect of the soya-based nitrogen source on patulin production could not be bypassed by any of the small chemical chromatin effectors. Interestingly, AZA influenced some SM cluster genes and SM production although no Aspergillus species has yet been shown to carry detectable DNA methylation.

Transcriptomic profiling of Aspergillus flavus in response to 5-azacytidine

Aspergillus flavus is a common saprophyte and opportunistic pathogen producing aflatoxin (AF) and many other secondary metabolites. 5-Azacytidine (5-AC), a derivative of the nucleoside cytidine, is widely used for studies in epigenetics and cancer biology as an inactivator of DNA methyltransferase and is also used for studying secondary metabolism in fungi. Our previous studies showed that 5-AC affects development and inhibits AF production in A. flavus, and that A. flavus lacks DNA methylation. In this study, an RNA-Seq approach was applied to explore the mechanism of 5-AC's effect on A. flavus. We identified 240 significantly differentially expressed (Q-value<0.05) genes after 5-AC treatment, including two backbone genes respectively in secondary metabolite clusters #27 and #35. These two clusters are involved in development or survival of sclerotia. GO functional enrichment analysis showed that these significantly differentially expressed genes were mainly involved in catalytic activity and proteolytic functions. The expressed transcripts of most genes in the AF biosynthetic gene cluster in A. flavus showed no significant changes after treatment with 5-AC and were expressed at low levels, and the transcription regulator genes aflR and aflS in this cluster did not show differential expression relative to the sample without 5-AC treatment. We found that the veA gene, which encodes protein bridges VelB and LaeA, decreased profoundly the expressed transcripts, and brlA, which encodes an early regulator of development, increased its transcripts in A. flavus after 5-AC treatment. Our data support a model whereby 5-AC affects development through increasing the expression of brlA by depressing the expression of veA and AF production through suppressing veA expression and dysregulating carboxypeptidase activity, which then prevents the aflatoxisomes (vesicles) from performing their normal function in AF formation. Furthermore, the suppressed veA expression weakens or even interrupts the connection between VelB and LaeA, leading to dysregulation of the expression pattern of genes involved in development and secondary metabolism in A. flavus. The RNA-seq data presented in this work were also served to improve the annotation of the A. flavus genome. This work provides a comprehensive view of the transcriptome of A. flavus responsive to 5-AC and supports the conclusion that fungal development and secondary metabolism are co-regulated.

Aflatoxins, fumonisins, and trichothecenes: a convergence of knowledge

Plant pathogenic fungi Aspergillus flavus, Fusarium verticillioides, and Fusarium graminearum infect seeds of the most important food and feed crops, including maize, wheat, and barley. More importantly, these fungi produce aflatoxins, fumonisins, and trichothecenes, respectively, which threaten health and food security worldwide. In this review, we examine the molecular mechanisms and environmental factors that regulate mycotoxin biosynthesis in each fungus, and discuss the similarities and differences in the collective body of knowledge. Whole-genome sequences are available for these fungi, providing reference databases for genomic, transcriptomic, and proteomic analyses. It is well recognized that genes responsible for mycotoxin biosynthesis are organized in clusters. However, recent research has documented the intricate transcriptional and epigenetic regulation that affects these gene clusters. Significantly, molecular networks that respond to environmental factors, namely nitrogen, carbon, and pH, are connected to components regulating mycotoxin production. Furthermore, the developmental status of seeds and specific tissue types exert conditional influences during fungal colonization. A comparison of the three distinct mycotoxin groups provides insight into new areas for research collaborations that will lead to innovative strategies to control mycotoxin contamination of grain.