Genomic clustering and co-regulation of transcriptional networks in the pathogenic fungus Fusarium graminearum

Trends in biological data integration for the selection of enzymes and transcription factors related to cellulose and hemicellulose degradation in fungi

Fungi are key players in biotechnological applications. Although several studies focusing on fungal diversity and genetics have been performed, many details of fungal biology remain unknown, including how cellulolytic enzymes are modulated within these organisms to allow changes in main plant cell wall compounds, cellulose and hemicellulose, and subsequent biomass conversion. With the advent and consolidation of DNA/RNA sequencing technology, different types of information can be generated at the genomic, structural and functional levels, including the gene expression profiles and regulatory mechanisms of these organisms, during degradation-induced conditions. This increase in data generation made rapid computational development necessary to deal with the large amounts of data generated. In this context, the origination of bioinformatics, a hybrid science integrating biological data with various techniques for information storage, distribution and analysis, was a fundamental step toward the current state-of-the-art in the postgenomic era. The possibility of integrating biological big data has facilitated exciting discoveries, including identifying novel mechanisms and more efficient enzymes, increasing yields, reducing costs and expanding opportunities in the bioprocess field. In this review, we summarize the current status and trends of the integration of different types of biological data through bioinformatics approaches for biological data analysis and enzyme selection.

Identification of Putative Virulence Genes by DNA Methylation Studies in the Cereal Pathogen Fusarium graminearum

DNA methylation mediates organisms’ adaptations to environmental changes in a wide range of species. We investigated if a such a strategy is also adopted by Fusarium graminearum in regulating virulence toward its natural hosts. A virulent strain of this fungus was consecutively sub-cultured for 50 times (once a week) on potato dextrose agar. To assess the effect of subculturing on virulence, wheat seedlings and heads (cv. A416) were inoculated with subcultures (SC) 1, 23, and 50. SC50 was also used to re-infect (three times) wheat heads (SC50×3) to restore virulence. In vitro conidia production, colonies growth and secondary metabolites production were also determined for SC1, SC23, SC50, and SC50×3. Seedling stem base and head assays revealed a virulence decline of all subcultures, whereas virulence was restored in SC50×3. The same trend was observed in conidia production. The DNA isolated from SC50 and SC50×3 was subject to a methylation content-sensitive enzyme and double-digest, restriction-site-associated DNA technique (ddRAD-MCSeEd). DNA methylation analysis indicated 1024 genes, whose methylation levels changed in response to the inoculation on a healthy host after subculturing. Several of these genes are already known to be involved in virulence by functional analysis. These results demonstrate that the physiological shifts following sub-culturing have an impact on genomic DNA methylation levels and suggest that the ddRAD-MCSeEd approach can be an important tool for detecting genes potentially related to fungal virulence.

Recombination suppression and evolutionary strata around mating-type loci in fungi: documenting patterns and understanding evolutionary and mechanistic causes

Genomic regions determining sexual compatibility often display recombination suppression, as occurs in sex chromosomes, plant self-incompatibility loci and fungal mating-type loci. Regions lacking recombination can extend beyond the genes determining sexes or mating types, by several successive steps of recombination suppression. Here we review the evidence for recombination suppression around mating-type loci in fungi, sometimes encompassing vast regions of the mating-type chromosomes. The suppression of recombination at mating-type loci in fungi has long been recognized and maintains the multiallelic combinations required for correct compatibility determination. We review more recent evidence for expansions of recombination suppression beyond mating-type genes in fungi ('evolutionary strata'), which have been little studied and may be more pervasive than commonly thought. We discuss testable hypotheses for the ultimate (evolutionary) and proximate (mechanistic) causes for such expansions of recombination suppression, including (1) antagonistic selection, (2) association of additional functions to mating-type, such as uniparental mitochondria inheritance, (3) accumulation in the margin of nonrecombining regions of various factors, including deleterious mutations or transposable elements resulting from relaxed selection, or neutral rearrangements resulting from genetic drift. The study of recombination suppression in fungi could thus contribute to our understanding of recombination suppression expansion across a broader range of organisms.

"Integrative genomic analysis of the bioprospection of regulators and accessory enzymes associated with cellulose degradation in a filamentous fungus (Trichoderma harzianum)"

BACKGROUND

Unveiling fungal genome structure and function reveals the potential biotechnological use of fungi. Trichoderma harzianum is a powerful CAZyme-producing fungus. We studied the genomic regions in T. harzianum IOC3844 containing CAZyme genes, transcription factors and transporters.

RESULTS

We used bioinformatics tools to mine the T. harzianum genome for potential genomics, transcriptomics, and exoproteomics data and coexpression networks. The DNA was sequenced by PacBio SMRT technology for multiomics data analysis and integration. In total, 1676 genes were annotated in the genomic regions analyzed; 222 were identified as CAZymes in T. harzianum IOC3844. When comparing transcriptome data under cellulose or glucose conditions, 114 genes were differentially expressed in cellulose, with 51 being CAZymes. CLR2, a transcription factor physically and phylogenetically conserved in Trichoderma spp., was differentially expressed under cellulose conditions. The genes induced/repressed under cellulose conditions included those important for plant biomass degradation, including CIP2 of the CE15 family and a copper-dependent LPMO of the AA9 family.

CONCLUSIONS

Our results provide new insights into the relationship between genomic organization and hydrolytic enzyme expression and regulation in T. harzianum IOC3844. Our results can improve plant biomass degradation, which is fundamental for developing more efficient strains and/or enzymatic cocktails to produce hydrolytic enzymes.

The birth, evolution and death of metabolic gene clusters in fungi

Fungi contain a remarkable diversity of both primary and secondary metabolic pathways involved in ecologically specialized or accessory functions. Genes in these pathways are frequently physically linked on fungal chromosomes, forming metabolic gene clusters (MGCs). In this Review, we describe the diversity in the structure and content of fungal MGCs, their population-level and species-level variation, the evolutionary mechanisms that underlie their formation, maintenance and decay, and their ecological and evolutionary impact on fungal populations. We also discuss MGCs from other eukaryotes and the reasons for their preponderance in fungi. Improved knowledge of the evolutionary life cycle of MGCs will advance our understanding of the ecology of specialized metabolism and of the interplay between the lifestyle of an organism and genome architecture.

Comprehensive Description of Fusarium graminearum Pigments and Related Compounds

Several studies have explored in depth the biochemistry and genetics of the pigments present in Fusarium graminearum, but there is a need to discuss their relationship with the mold's observable surface color pattern variation throughout its lifecycle. Furthermore, they require basic cataloguing, including a description of their major features known so far. Colors are a viable alternative to size measurement in growth studies. When grown on yeast extract agar (YEA) at 25 °C, F. graminearum initially exhibits a whitish mycelium, developing into a yellow-orange mold by the sixth day and then turning into wine-red. The colors are likely due to accumulation of the golden yellow polyketide aurofusarin and the red rubrofusarin, but the carotenoid neurosporaxanthin also possibly plays a major role in the yellow or orange coloration. Torulene might contribute to red tones, but it perhaps ends up being converted into neurosporaxanthin. Culmorin is also present, but it does not contribute to the color, though it was initially isolated in pigment studies. Additionally, there is the 5-deoxybostrycoidin-based melanin, but it mostly occurs in the teleomorph's perithecium. There is still a need to chemically quantify the pigments throughout the lifecycle, and analyze their relationships and how much each impacts F. graminearum's surface color.

Whole-Genome Resequencing and Pan-Transcriptome Reconstruction Highlight the Impact of Genomic Structural Variation on Secondary Metabolite Gene Clusters in the Grapevine Esca Pathogen Phaeoacremonium minimum

The Ascomycete fungus Phaeoacremonium minimum is one of the primary causal agents of Esca, a widespread and damaging grapevine trunk disease. Variation in virulence among Pm. minimum isolates has been reported, but the underlying genetic basis of the phenotypic variability remains unknown. The goal of this study was to characterize intraspecific genetic diversity and explore its potential impact on virulence functions associated with secondary metabolism, cellular transport, and cell wall decomposition. We generated a chromosome-scale genome assembly, using single molecule real-time sequencing, and resequenced the genomes and transcriptomes of multiple isolates to identify sequence and structural polymorphisms. Numerous insertion and deletion events were found for a total of about 1 Mbp in each isolate. Structural variation in this extremely gene dense genome frequently caused presence/absence polymorphisms of multiple adjacent genes, mostly belonging to biosynthetic clusters associated with secondary metabolism. Because of the observed intraspecific diversity in gene content due to structural variation we concluded that a transcriptome reference developed from a single isolate is insufficient to represent the virulence factor repertoire of the species. We therefore compiled a pan-transcriptome reference of Pm. minimum comprising a non-redundant set of 15,245 protein-coding sequences. Using naturally infected field samples expressing Esca symptoms, we demonstrated that mapping of meta-transcriptomics data on a multi-species reference that included the Pm. minimum pan-transcriptome allows the profiling of an expanded set of virulence factors, including variable genes associated with secondary metabolism and cellular transport.

Condition-dependent co-regulation of genomic clusters of virulence factors in the grapevine trunk pathogen Neofusicoccum parvum

The ascomycete Neofusicoccum parvum, one of the causal agents of Botryosphaeria dieback, is a destructive wood-infecting fungus and a serious threat to grape production worldwide. The capability to colonize woody tissue, combined with the secretion of phytotoxic compounds, is thought to underlie its pathogenicity and virulence. Here, we describe the repertoire of virulence factors and their transcriptional dynamics as the fungus feeds on different substrates and colonizes the woody stem. We assembled and annotated a highly contiguous genome using single-molecule real-time DNA sequencing. Transcriptome profiling by RNA sequencing determined the genome-wide patterns of expression of virulence factors both in vitro (potato dextrose agar or medium amended with grape wood as substrate) and in planta. Pairwise statistical testing of differential expression, followed by co-expression network analysis, revealed that physically clustered genes coding for putative virulence functions were induced depending on the substrate or stage of plant infection. Co-expressed gene clusters were significantly enriched not only in genes associated with secondary metabolism, but also in those associated with cell wall degradation, suggesting that dynamic co-regulation of transcriptional networks contributes to multiple aspects of N. parvum virulence. In most of the co-expressed clusters, all genes shared at least a common motif in their promoter region, indicative of co-regulation by the same transcription factor. Co-expression analysis also identified chromatin regulators with correlated expression with inducible clusters of virulence factors, suggesting a complex, multi-layered regulation of the virulence repertoire of N. parvum.

A Global Coexpression Network Approach for Connecting Genes to Specialized Metabolic Pathways in Plants

Plants produce diverse specialized metabolites (SMs), but the genes responsible for their production and regulation remain largely unknown, hindering efforts to tap plant pharmacopeia. Given that genes comprising SM pathways exhibit environmentally dependent coregulation, we hypothesized that genes within a SM pathway would form tight associations (modules) with each other in coexpression networks, facilitating their identification. To evaluate this hypothesis, we used 10 global coexpression data sets, each a meta-analysis of hundreds to thousands of experiments, across eight plant species to identify hundreds of coexpressed gene modules per data set. In support of our hypothesis, 15.3 to 52.6% of modules contained two or more known SM biosynthetic genes, and module genes were enriched in SM functions. Moreover, modules recovered many experimentally validated SM pathways, including all six known to form biosynthetic gene clusters (BGCs). In contrast, bioinformatically predicted BGCs (i.e., those lacking an associated metabolite) were no more coexpressed than the null distribution for neighboring genes. These results suggest that most predicted plant BGCs are not genuine SM pathways and argue that BGCs are not a hallmark of plant specialized metabolism. We submit that global gene coexpression is a rich, largely untapped resource for discovering the genetic basis and architecture of plant natural products.

Genome-Wide Association Study Identifies Novel Candidate Genes for Aggressiveness, Deoxynivalenol Production, and Azole Sensitivity in Natural Field Populations of Fusarium graminearum

Genome-wide association studies can identify novel genomic regions and genes that affect quantitative traits. Fusarium head blight is a destructive disease caused by Fusarium graminearum that exhibits several quantitative traits, including aggressiveness, mycotoxin production, and fungicide resistance. Restriction site-associated DNA sequencing was performed for 220 isolates of F. graminearum. A total of 119 isolates were phenotyped for aggressiveness and deoxynivalenol (DON) production under natural field conditions across four environments. The effective concentration of propiconazole that inhibits isolate growth in vitro by 50% was calculated for 220 strains. Approximately 29,000 single nucleotide polymorphism markers were associated to each trait, resulting in 50, 29, and 74 quantitative trait nucleotides (QTNs) that were significantly associated to aggressiveness, DON production, and propiconazole sensitivity, respectively. Approximately 41% of these QTNs caused nonsynonymous substitutions in predicted exons, while the remainder were synonymous substitutions or located in intergenic regions. Three QTNs associated with propiconazole sensitivity were significant after Bonferroni correction. These QTNs were located in genes not previously associated with azole sensitivity. The majority of the detected QTNs were located in genes with predicted regulatory functions, suggesting that nucleotide variation in regulatory genes plays a major role in the corresponding quantitative trait variation.

The Fdb3 transcription factor of the Fusarium Detoxification of Benzoxazolinone gene cluster is required for MBOA but not BOA degradation in Fusarium pseudograminearum

A number of cereals produce the benzoxazolinone class of phytoalexins. Fusarium species pathogenic towards these hosts can typically degrade these compounds via an aminophenol intermediate, and the ability to do so is encoded by a group of genes found in the Fusarium Detoxification of Benzoxazolinone (FDB) cluster. A zinc finger transcription factor encoded by one of the FDB cluster genes (FDB3) has been proposed to regulate the expression of other genes in the cluster and hence is potentially involved in benzoxazolinone degradation. Herein we show that Fdb3 is essential for the ability of Fusarium pseudograminearum to efficiently detoxify the predominant wheat benzoxazolinone, 6-methoxy-benzoxazolin-2-one (MBOA), but not benzoxazoline-2-one (BOA). Furthermore, additional genes thought to be part of the FDB gene cluster, based upon transcriptional response to benzoxazolinones, are regulated by Fdb3. However, deletion mutants for these latter genes remain capable of benzoxazolinone degradation, suggesting that they are not essential for this process.

Degradation of the benzoxazolinone class of phytoalexins is important for virulence of Fusarium pseudograminearum towards wheat

Wheat, maize, rye and certain other agriculturally important species in the Poaceae family produce the benzoxazolinone class of phytoalexins on pest and pathogen attack. Benzoxazolinones can inhibit the growth of pathogens. However, certain fungi can actively detoxify these compounds. Despite this, a clear link between the ability to detoxify benzoxazolinones and pathogen virulence has not been shown. Here, through comparative genome analysis of several Fusarium species, we have identified a conserved genomic region around the FDB2 gene encoding an N-malonyltransferase enzyme known to be involved in benzoxazolinone degradation in the maize pathogen Fusarium verticillioides. Expression analyses demonstrated that a cluster of nine genes was responsive to exogenous benzoxazolinone in the important wheat pathogen Fusarium pseudograminearum. The analysis of independent F. pseudograminearum FDB2 knockouts and complementation of the knockout with FDB2 homologues from F. graminearum and F. verticillioides confirmed that the N-malonyltransferase enzyme encoded by this gene is central to the detoxification of benzoxazolinones, and that Fdb2 contributes quantitatively to virulence towards wheat in head blight inoculation assays. This contrasts with previous observations in F. verticillioides, where no effect of FDB2 mutations on pathogen virulence towards maize was observed. Overall, our results demonstrate that the detoxification of benzoxazolinones is a strategy adopted by wheat-infecting F. pseudograminearum to overcome host-derived chemical defences.

The Fusarium graminearum genome reveals more secondary metabolite gene clusters and hints of horizontal gene transfer

Fungal secondary metabolite biosynthesis genes are of major interest due to the pharmacological properties of their products (like mycotoxins and antibiotics). The genome of the plant pathogenic fungus Fusarium graminearum codes for a large number of candidate enzymes involved in secondary metabolite biosynthesis. However, the chemical nature of most enzymatic products of proteins encoded by putative secondary metabolism biosynthetic genes is largely unknown. Based on our analysis we present 67 gene clusters with significant enrichment of predicted secondary metabolism related enzymatic functions. 20 gene clusters with unknown metabolites exhibit strong gene expression correlation in planta and presumably play a role in virulence. Furthermore, the identification of conserved and over-represented putative transcription factor binding sites serves as additional evidence for cluster co-regulation. Orthologous cluster search provided insight into the evolution of secondary metabolism clusters. Some clusters are characteristic for the Fusarium phylum while others show evidence of horizontal gene transfer as orthologs can be found in representatives of the Botrytis or Cochliobolus lineage. The presented candidate clusters provide valuable targets for experimental examination.